Binding Affinity

The dissociation constant (Kd) is the gold-standard measure of affinity, representing the ligand concentration at which half the target's binding sites are occupied. A lower Kd indicates higher affinity.

• Formula: Kd = [P][L] / [PL]
• Units: Molar (M), typically ranging from mM (weak) to pM (ultra-tight)
• Relationship: Kd is the reciprocal of the association constant (Ka)
• Measurement: Determined via Surface Plasmon Resonance (SPR), Isothermal Titration Calorimetry (ITC), or radioligand binding assays
• Example: Biotin-streptavidin complex exhibits a Kd of ~10⁻¹⁴ M, one of the strongest known non-covalent interactions

The Gibbs free energy change (ΔG) upon binding determines spontaneity and is directly related to Kd. A more negative ΔG signifies a more favorable, higher-affinity interaction.

• Equation: ΔG = -RT ln(Ka) = RT ln(Kd)
• Components: ΔG = ΔH - TΔS, where ΔH is enthalpy change and ΔS is entropy change
• Energy Scale: A ΔG of -10 kcal/mol corresponds to a Kd of ~50 nM at 298 K
• Enthalpy-Driven: Dominated by favorable hydrogen bonds, van der Waals contacts, and electrostatic interactions
• Entropy-Driven: Driven by the hydrophobic effect and release of ordered water molecules from the binding interface

Enthalpy-entropy compensation is the observed phenomenon where favorable enthalpic gains from forming intermolecular bonds are often offset by unfavorable entropic losses from conformational restriction, and vice versa.

• Mechanism: Stronger binding interactions (favorable ΔH) rigidify the complex (unfavorable -TΔS)
• Solvent Reorganization: Release of structured water from hydrophobic surfaces is entropically favorable but enthalpically neutral or slightly unfavorable
• Drug Design Implication: Optimizing for enthalpy-driven binders often yields higher selectivity, as entropic contributions from desolvation are less target-specific
• ITC Deconvolution: Isothermal Titration Calorimetry uniquely measures ΔH and ΔS independently, revealing the thermodynamic signature of binding

Binding affinity is a ratio of kinetic rate constants. The association rate (kon) and dissociation rate (koff) define how quickly a complex forms and decays, with Kd = koff / kon.

• kon: Diffusion-limited maximum ~10⁸ to 10⁹ M⁻¹s⁻¹; influenced by electrostatic steering and conformational selection
• koff: Ranges from milliseconds (weak binders) to days (tight binders); the primary determinant of residence time
• Residence Time (τ): τ = 1/koff; often correlates better with in vivo efficacy than Kd alone
• SPR Biosensors: Surface Plasmon Resonance provides real-time kon and koff measurements without labels
• Clinical Relevance: Drugs with long residence times (slow koff) can maintain target engagement even after plasma clearance

Desolvation is the energetic cost of stripping water molecules from the binding interfaces of both protein and ligand before complex formation. This process profoundly influences the net binding free energy.

• Hydrophobic Effect: Burial of non-polar surface area releases ordered water into bulk solvent, providing a favorable entropic driving force
• Polar Desolvation Penalty: Removing water from charged or polar groups is energetically costly and must be compensated by new hydrogen bonds in the complex
• Water Networks: Structured water molecules bridging protein and ligand can contribute favorably to binding enthalpy if optimally positioned
• Computational Prediction: Implicit solvation models (GBSA, PBSA) approximate these effects, while explicit solvent molecular dynamics provide detailed water maps

IC50 (half-maximal inhibitory concentration) and EC50 (half-maximal effective concentration) are functional measures of a compound's potency, which depend on both binding affinity and downstream cellular context.

• IC50: Concentration required to inhibit a biological process by 50%; commonly used in enzyme and cell-based assays
• EC50: Concentration producing 50% of the maximal response; used for agonists and activators
• Cheng-Prusoff Equation: Relates IC50 to Ki (inhibition constant) for competitive inhibitors: Ki = IC50 / (1 + [S]/Km)
• Assay Dependency: IC50 values vary with substrate concentration, incubation time, and cell type, unlike the system-independent Kd
• Caution: IC50 is not a direct measure of binding affinity; it conflates target engagement with functional pathway amplification

A mathematical function that approximates the binding free energy (ΔG) of a protein-ligand complex. It evaluates non-covalent interactions—such as van der Waals forces, hydrogen bonds, and electrostatic complementarity—to rank docking poses.

• Force-field-based: Calculates physical interaction energies (e.g., DOCK, AutoDock).
• Empirical: Sums weighted terms fitted to experimental affinity data (e.g., ChemScore).
• Knowledge-based: Derives statistical potentials from observed atom-pair frequencies in structural databases (e.g., PMF).
• Consensus scoring combines multiple functions to improve hit-rate enrichment.

A rigorous alchemical free energy calculation method based on statistical mechanics. FEP computes the relative binding affinity between two related ligands by simulating a non-physical thermodynamic path where one molecule is gradually 'mutated' into another.

• Achieves chemical accuracy (~1 kcal/mol) when properly converged.
• Requires extensive molecular dynamics sampling, making it computationally expensive.
• Modern GPU-accelerated implementations (e.g., Schrödinger FEP+) have made it viable for lead optimization campaigns.

Molecular Mechanics/Generalized Born Surface Area is an end-point free energy method that estimates binding affinity from only the initial and final states of a simulation—the free ligand, free protein, and the complex.

• Computationally efficient compared to FEP, requiring only a few hundred snapshots from an MD trajectory.
• Decomposes binding energy into per-residue contributions, enabling identification of hot-spot residues.
• Less accurate for absolute binding free energies but effective for ranking congeneric series.

A binary or count-based vector encoding the specific intermolecular contacts between a protein's residues and a bound ligand. Each bit represents a defined interaction type at a specific residue position.

• Interaction types: Hydrogen bond donor/acceptor, hydrophobic contact, π-π stacking, salt bridge, halogen bond.
• Used as features for machine learning models predicting binding affinity or classifying active vs. inactive compounds.
• Enables interaction pattern similarity searching to find ligands with analogous binding modes regardless of chemical scaffold.

The reciprocal of the dissociation rate constant (koff), representing the average duration a ligand remains bound to its target. While binding affinity (Kd) reflects the equilibrium between association and dissociation, residence time captures the kinetic stability of the complex.

• τ = 1/koff: A drug with a residence time of hours may maintain target occupancy long after plasma concentration drops.
• Correlates with prolonged pharmacodynamic effects and improved in vivo efficacy.
• Optimized through structure-kinetic relationships (SKRs) rather than purely thermodynamic considerations.

A retrospective performance metric quantifying how effectively a virtual screening method identifies active compounds early in a ranked list. It measures the fold-improvement over random selection.

• EFx% = (Actives found in top x%) / (Expected actives by random chance).
• An EF1% of 50 means the method found 50 times more actives in the top 1% of the ranked database than random screening would.
• Evaluated using decoy sets (presumed inactives with similar physical properties to known actives) to prevent artificial enrichment from trivial property discrimination.