Module 3: Fundamentals of Ligand-Receptor Interactions
Analyze the energetic drivers behind ligand binding. Explore Gibbs free energy, calculate potential energy curves, and inspect the structural basis of the hydrophobic effect.
Learning outcomes
- Decompose binding free energy into enthalpic and entropic contributions.
- Identify the major non-covalent interactions and their characteristic geometries.
- Explain the hydrophobic effect and why desolvation can dominate an affinity change.
- Use ligand efficiency and LLE to compare hits of different sizes and lipophilicities.
- Name the four ways the pairwise-contact picture of binding breaks down.
1. Molecular Recognition Models
How do a drug molecule and its target protein recognize each other in a crowded biological environment? Two classical theories explain this process:
Lock-and-Key Model
Formulated by Emil Fischer in 1894. Suggests the receptor and ligand possess complementary, rigid geometries. It explains high specificity but fails to account for structural plasticity.
Induced-Fit Model
Proposed by Daniel Koshland in 1958. Proposes that ligand binding triggers conformational rearrangements in the receptor pocket to optimize contacts, matching thermodynamic realities.
2. Thermodynamics of Binding
The affinity of a ligand for its receptor is defined by the change in Gibbs Free Energy (ΔG) upon complex formation. A more negative ΔG corresponds to tighter binding (higher affinity):
Gibbs Free Energy change is directly related to the binding dissociation constant (K_d) by the fundamental thermodynamic equilibrium equation:
Where R is the gas constant and T is the absolute temperature. Because of this logarithmic relationship, small, linear changes in Gibbs Free Energy translate into exponential improvements in binding affinity:
| ΔG (at 298 K) | Affinity Constant (K_d) | Qualitative Description |
|---|---|---|
| -4.1 kcal/mol | 1.0 mM (Millimolar) | Very weak binding, typical of small fragments. |
| -8.2 kcal/mol | 1.0 µM (Micromolar) | Moderate binding, typical of high-throughput screen hits. |
| -12.3 kcal/mol | 1.0 nM (Nanomolar) | Tightly bound drug candidate. |
| -16.4 kcal/mol | 1.0 pM (Picomolar) | Exceptional affinity binding, rare and highly optimized. |
Enthalpy (ΔH)
Refers to the release of heat resulting from the formation of specific, directional non-covalent contacts (hydrogen bonds, ionic pairs, van der Waals, and halogen bonds) between the ligand and target.
Entropy (-TΔS)
Represents the change in disorder. Complexation restricts ligand and side-chain rotation, costing conformational entropy. However, this penalty is offset by the hydrophobic effect: water displacement from hydrophobic surfaces into bulk.
Interactive Playground: Distance-Dependent H-Bond Potential
Adjust the distance slider to bring the Hydrogen Bond donor atom closer to the acceptor. Watch the energy shift along the potential curve.
-2.17kcal/mol
This 12-6 radial potential is a teaching approximation for short-range repulsion and an attractive well. Real hydrogen-bond strength also depends on donor-acceptor angle, protonation, solvent, and the surrounding electrostatic field.
3. From Affinity to Design Efficiency
Potency is essential, but it does not show what a molecule had to become to achieve that potency. Medicinal chemists use efficiency metrics and thermodynamic measurements to compare compounds while tracking size, lipophilicity, and binding mechanism.
| Measure | What it adds | Important limit |
|---|---|---|
| Ligand efficiency (LE) | Relates binding free energy to non-hydrogen atom count, helping compare differently sized hits. | Size normalization has known biases; compare related series and use the same affinity endpoint. |
| Lipophilic ligand efficiency (LLE or LipE) | pActivity - logD asks whether potency rises faster than lipophilicity. | State the assay endpoint, pH, and lipophilicity measure. The octanol-water reference is useful, not a universal measure of specificity. |
| Isothermal titration calorimetry (ITC) | A titration can estimate affinity, stoichiometry, and binding enthalpy; entropy is inferred from the free-energy relationship. | Buffer ionization, proton transfer, concentration accuracy, and coupled conformational changes can affect the observed heat. |
Enthalpy-entropy compensation
Two ligands can have similar affinity with different enthalpic and entropic profiles. Solvent reorganization, protonation, and conformational change couple the terms, so an apparently favorable enthalpy does not automatically identify a better lead.
Water is part of the mechanism
Displacing an unfavorable water can help binding, while disrupting a stable bridging network can hurt it. Inspect water networks and receptor state instead of treating every buried water release as an automatic gain.
4. Types of Non-Covalent Interactions
Each interaction below is a physical phenomenon first and a line of arithmetic second. Module 4 shows how a force field turns them into computable terms — hydrogen bonds and salt bridges fall out of the Coulomb term, dispersion and steric clash out of the Lennard-Jones term — and every docking score (Module 6) and free-energy estimate (Module 10) you meet later is built from exactly these contributions.
Hydrogen Bonds
Formed between a hydrogen atom covalently bound to an electronegative atom (Donor: O-H, N-H) and another electronegative atom with lone pairs (Acceptor: O, N). Strong, highly directional, with typical optimal donor-acceptor distances of 2.7–3.2 Å and bond angles close to 180°.
Halogen Bonds (σ-Hole Interactions)
An interaction between an electronegative atom and the electropositive region on the tip of a halogen atom (Cl, Br, or I) bound to carbon. This positive region, known as the σ-hole, renders halogen bonds highly directional.
Electrostatic Interactions
Salt bridges formed between oppositely charged functional groups (e.g. protonated amine on ligand and carboxylate side-chain of Aspartate or Glutamate on receptor). These interactions are long-range (V ∝ 1/r).
Cation-π Interactions
A special ion-dipole interaction between a cation (e.g. a protonated amine, Lys-NH₃⁺, or Arg guanidinium) and the electron-rich π face of an aromatic ring. Because the effect depends on ring electron density, electron-donating substituents (e.g. -NH₂) strengthen it while electron-withdrawing groups (e.g. -CN) weaken it. This explains why electron-rich Trp and Tyr engage in cation-π contacts far more often than Phe.
Van der Waals / London Dispersion & π-π Stacking
Weak (~0.5–1 kcal/mol each), short-range forces from transient, induced dipoles between all atoms in close contact. Individually negligible, but summed over a well-packed binding pocket they contribute substantially to affinity — this is the energetic basis of shape complementarity. π-π stacking between aromatic rings (parallel-displaced or T-shaped, ~3.5–4.0 Å) is a directional special case.
Interactive Playground: Hydrophobic Effect & Desolvation
Toggle the "Bind Ligand" button to push a lipophilic compound into a hydrophobic pocket. The slate circles illustrate interfacial waters that can be released into bulk solution when nonpolar surfaces are buried.
Thermodynamic Analysis
Both the hydrophobic pocket and the lipophilic ligand are surrounded by highly organized, 'caged' water structures. This localized structure is translationally and rotationally restricted, resulting in low system entropy.
Interactive Playground: The Supramolecular Binding Sandbox
In physical drug discovery, binding affinity depends on matching multiple functional groups simultaneously. Drag the central ligand core (or use the coordinates sliders) to align the three chemical arms with their target pocket residues. Watch the Gibbs Free Energy update live. Note: Aligning all three perfectly is restricted by scaffold geometry; you must resolve the conformational strain trade-off to minimize ΔG!
Thermodynamic Calculations
5. Where This Picture Breaks Down
Everything above treats binding as a sum of pairwise contacts between a ligand and a rigid pocket. That picture is useful enough to design drugs with, and wrong in four specific ways that resurface throughout the course.
Interactions are not additive
Adding a hydrogen bond worth 2 kcal/mol to a ligand rarely buys 2 kcal/mol of affinity. The new contact may cost desolvation, restrict a rotatable bond, or strain the pose. Non-additivity is precisely why medicinal chemistry still requires synthesis rather than arithmetic — and why scoring functions (Module 6) fail in the way they do.
The pocket is not rigid
Induced fit above is a cartoon of a much larger effect: side chains rotate, loops close, and some pockets do not exist until a ligand arrives. Anything computed on one fixed structure inherits that structure's assumptions (Modules 6 and 10).
Water is a participant, not a background
Every contact you form must first break a contact with water, and a few ordered waters in a pocket can be worth more than a whole substituent. Desolvation is the single most commonly underestimated term in this module.
Affinity is not the objective
ΔG tells you how tightly a ligand binds its target — not whether it is selective, absorbed, metabolically stable, or safe. Ligand efficiency and LLE exist to keep potency honest, and Modules 12 and 15 supply the constraints that ultimately decide whether a tight binder becomes a drug.
Every one of these is a physical shortcoming of the pairwise model, and each is the reason a later technique exists. Module 4 turns these same interactions into computable energy terms; keep the four caveats in mind, because they explain most of what goes wrong afterwards.