Abstract:Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.

DeltaGzip: Computing Biopolymer-Ligand Binding Affinity via Kolmogorov Complexity and Lossless Compression

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities

Analysis of Compression Techniques for DNA Sequence Data

A Rigorous Framework for Calculating Protein-Protein Binding Affinities in Membranes

Gap-Δenergy, a New Metric of the Bond Energy State, Assisting to Predict Molecular Toxicity

Recent theoretical and computational advances for modeling protein–ligand binding affinities

DeltaDock: A Unified Framework for Accurate, Efficient, and Physically Reliable Molecular Docking

Computing Ligand Binding Free Energy in a Large Flexible Pocket of a Large Protein

Estimating protein–ligand interactions with geometric deep learning and mixture density models

Improved Modeling of Peptide-Protein Binding Through Global Docking and Accelerated Molecular Dynamics Simulations

Ligand Gaussian accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides

Massively parallel free energy calculations for affinity maturation of designed miniproteins

DynamicBind: predicting ligand-specific protein-ligand complex structure with a deep equivariant generative model

Free-energy calculations for semi-flexible macromolecules: Applications to DNA knotting and looping

Achieving Rapid Free Energy Estimates from λ-Dynamics with Bias Updated Gibbs Sampling

Predictions Of Binding Of A Diverse Set Of Ligands To Gelatinase-A By A Combination Of Molecular Dynamics And Continuum Solvent Models

Prediction of atomic structure from sequence for double helical DNA oligomers

A critical analysis of computational protein design with sparse residue interaction graphs

Rapid calculation of side chain packing and free energy with applications to protein molecular dynamics

Insights into the Molecular Mechanisms of Protein-Ligand Interactions by Molecular Docking and Molecular Dynamics Simulation: A Case of Oligopeptide Binding Protein