During their transport cycle, ATP-binding cassette (ABC) transporters undergo large-scale conformational changes between inward- and outward-facing states. Using an approach based on designing system-specific reaction coordinates and using nonequilibrium work relations, we have performed extensive all-atom molecular dynamics simulations in the presence of explicit membrane/solvent to sample a large number of mechanistically distinct pathways for the conformational transition of MsbA, a bacterial ABC exporter whose structure has been solved in multiple functional states. The computational approach developed here is based on (i) extensive exploration of system-specific biasing protocols (e.g., using collective variables designed based on available low-resolution crystal structures) and (ii) using nonequilibrium work relations for comparing the relevance of the transition pathways. The most relevant transition pathway identified using this approach involves several distinct stages reflecting the complex nature of the structural changes associated with the function of the protein. The opening of the cytoplasmic gate during the outward- to inward-facing transition of apo MsbA is found to be disfavored when the periplasmic gate is open and facilitated by a twisting motion of the nucleotide-binding domains that involves a dramatic change in their relative orientation. These results highlight the cooperativity between the transmembrane and the nucleotide-binding domains in the conformational transition of ABC exporters. The approach introduced here provides a framework to study large-scale conformational changes of other membrane transporters whose computational investigation at an atomic resolution may not be currently feasible using conventional methods. Significance Membrane transporters rely on large-scale conformational changes for their function. Despite considerable effort, however, details of such structural transitions are largely unknown, leaving many fundamental questions regarding the transport mechanism in this important family of membrane proteins unanswered. Here, we present a nonequilibrium approach to characterize the conformational transition of MsbA, a member of the ATP-binding cassette exporter family, which is involved in transport of diverse substrates across the membrane. The design of the study is based on complex, system-specific reaction coordinates and protocols resulting in an unprecedented level of detail on the nature of conformational coupling and the mechanism of transport. The presented approach opens opportunities for investigating large-scale conformational changes of other membrane transporters.

Membrane transporters actively translocate their substrate by undergoing large-scale structural transitions between inward- (IF) and outward-facing (OF) states (‘alternating-access’ mechanism). Despite extensive structural studies, atomic-level mechanistic details of such structural transitions, and as importantly, their coupling to chemical events supplying the energy, remain amongst the most elusive aspects of the function of these proteins. Here we present a quantitative, atomic-level description of the functional thermodynamic cycle for the glycerol-3-phosphate:phosphate antiporter GlpT by using a novel approach in reconstructing the free energy landscape governing the IF↔OF transition along a cyclic transition pathway involving both apo and substrate-bound states. Our results provide a fully atomic description of the complete transport process, offering a structural model for the alternating-access mechanism and substantiating the close coupling between global structural transitions and local chemical events.

The structure of the proline amino acid allows folded polyproline peptides to exist as both left- (PPII) and right-handed (PPI) helices. We have characterized the free energy landscapes of hexamer, nanomer, and tridecamer polyproline peptides in gas phase and implicit water as well as explicit hexane and 1-propanol for the nanomer. To enhance the sampling provided by regular molecular dynamics, we used the recently developed adaptively biased molecular dynamics method, which describes Landau free energy maps in terms of relevant collective variables. These maps, as a function of the collective variables of handedness, radius of gyration, and three others based on the peptide torsion angle ω, were used to determine the relative stability of the different structures, along with an estimate of the transition pathways connecting the different minima. Results show the existence of several metastable isomers and therefore provide a complementary view to experimental conclusions based on photo-induced electron transfer experiments with regard to the existence of stable heterogeneous subpopulations in PPII polyproline.

Characterizing large-scale structural transitions in biomolecular systems poses major technical challenges to both experimental and computational approaches. On the computational side, efficient sampling of the configuration space along the transition pathway remains the most daunting challenge. Recognizing this issue, we introduce a knowledge-based computational approach toward describing large-scale conformational transitions using (i) nonequilibrium, driven simulations combined with work measurements and (ii) free energy calculations using empirically optimized biasing protocols. The first part is based on designing mechanistically relevant, system-specific reaction coordinates whose usefulness and applicability in inducing the transition of interest are examined using knowledge-based, qualitative assessments along with nonequilirbrium work measurements which provide an empirical framework for optimizing the biasing protocol. The second part employs the optimized biasing protocol resulting from the first part to initiate free energy calculations and characterize the transition quantitatively. Using a biasing protocol fine-tuned to a particular transition not only improves the accuracy of the resulting free energies but also speeds up the convergence. The efficiency of the sampling will be assessed by employing dimensionality reduction techniques to help detect possible flaws and provide potential improvements in the design of the biasing protocol. Structural transition of a membrane transporter will be used as an example to illustrate the workings of the proposed approach.

The protein–ligand binding affinity quantifies the binding strength between a protein and its ligand. Computer modeling and simulations can be used to estimate the binding affinity or binding free energy using data- or physics-driven methods or a combination thereof. Here we discuss a purely physics-based sampling approach based on biased molecular dynamics simulations. Our proposed method generalizes and simplifies previously suggested stratification strategies that use umbrella sampling or other enhanced sampling simulations with additional collective-variable-based restraints. The approach presented here uses a flexible scheme that can be easily tailored for any system of interest. We estimate the binding affinity of human fibroblast growth factor 1 to heparin hexasaccharide based on the available crystal structure of the complex as the initial model and four different variations of the proposed method to compare against the experimentally determined binding affinity obtained from isothermal titration calorimetry experiments.