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.

Calculations of the Binding Affinities of Protein-Protein Complexes with the Fast Multipole Method

Calculation of Relative Binding Free Energy for Mutations in Protein Complexes: The Alchemical Path

A Discontinuous Potential Model for Protein-Protein Interactions.

Calculation of Protein–ligand Binding Affinities Based on a Fragment Quantum Mechanical Method

Binding Free Energy Estimation for Protein-Ligand Complex Based on MM-PBSA with Various Partial Charge Models.

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Exploring the Binding Mechanism Between Human Profilin (PFN1) and Polyproline-10 Through Binding Mode Screening

Accurate and Efficient Calculation of Protein-Protein Binding Free Energy-Interaction Entropy with Residue Type-Specific Dielectric Constants.

Fragment-based Quantum Mechanical Calculation of Protein-Protein Binding Affinities

Calculation of Absolute Protein-Ligand Binding Affinity Using Path and Endpoint Approaches

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

Computational prediction of protein–protein binding affinities

A Computational Analysis of the Binding Affinities of FKBP12 Inhibitors Using the MM‐PB/SA Method

Robust prediction of relative binding energies for protein-protein complex mutations using free energy perturbation calculations

Efficient Calculation of Protein-Ligand Binding Free Energy Using GFN Methods: the Power of the Cluster Model.

Binding free-energy calculation of an ion-peptide complex by constrained dynamics.

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

Calculating Structure Factors of Protein Solutions by Atomistic Modeling of Protein-Protein Interactions

Binding free energy of protein/ligand complexes calculated using dissociation Parallel Cascade Selection Molecular Dynamics and Markov state model

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