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.

Free-energy landscapes of transmembrane homodimers by bias-exchange adaptively biased molecular dynamics

Free energy, rates, and mechanism of transmembrane dimerization in lipid bilayers from dynamically unbiased molecular dynamics simulations

Mapping transiently formed and sparsely populated conformations on a complex energy landscape

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

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

Exploration of Multi-State Conformational Dynamics and Underlying Global Functional Landscape of Maltose Binding Protein

Protein Conformational Transitions from All-Atom Adaptively Biased Path Optimization

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

Balancing Force Field Protein–Lipid Interactions To Capture Transmembrane Helix–Helix Association

Binding Energy Distribution Analysis Method: Hamiltonian Replica Exchange with Torsional Flattening for Binding Mode Prediction and Binding Free Energy Estimation

Efficiently determining membrane-bound conformations of peripheral membrane proteins using replica exchange with hybrid tempering

Massively parallel free energy calculations for affinity maturation of designed miniproteins

Thermodynamic free energy methods to investigate shape transitions in bilayer membranes

Critical residues and motifs for homodimerization of the first transmembrane domain of the plasma membrane glycoprotein CD36

Exploring Free Energy Landscapes for Protein Partitioning into Membrane Domains in All-Atom and Coarse-Grained Simulations

All-Atom Molecular Dynamics Elucidating Molecular Mechanisms of Single-Transmembrane Model Peptide Dimerization in a Lipid Bilayer

Assessing the performance of the MM/PBSA and MM/GBSA methods. 10. Impacts of enhanced sampling and variable dielectric model on protein-protein Interactions

Application of the Accelerated Molecular Dynamics Simulations to the Folding of A Small Protein

Computing committors via Mahalanobis diffusion maps with enhanced sampling data

Elucidating Protein Structures in the Gas Phase: Traversing Configuration Space with Biasing Methods