Mechanistic picture for conformational transition of a membrane transporter at atomic resolution

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. 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.

[1]  A. M. George,et al.  Mechanism of ABC transporters: A molecular dynamics simulation of a well characterized nucleotide-binding subunit , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[2]  C. Roland,et al.  Calculating relative transition rates with driven nonequilibrium simulations , 2011 .

[3]  Jinhui Dong,et al.  Structural Basis of Energy Transduction in the Transport Cycle of MsbA , 2005, Science.

[4]  N. Mousseau,et al.  Molecular dynamics simulations of the bacterial ABC transporter SAV1866 in the closed form. , 2012, The journal of physical chemistry. B.

[5]  R. Dawson,et al.  Structure of a bacterial multidrug ABC transporter , 2006, Nature.

[6]  Yue Weng,et al.  Structure of P-Glycoprotein Reveals a Molecular Basis for Poly-Specific Drug Binding , 2009, Science.

[7]  C. Klug,et al.  Effects of the L511P and D512G mutations on the Escherichia coli ABC transporter MsbA. , 2011, Biochemistry.

[8]  Jue Chen,et al.  Structural insights into ABC transporter mechanism. , 2008, Current opinion in structural biology.

[9]  C Jarzynski,et al.  How does a system respond when driven away from thermal equilibrium? , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[10]  Y. Sugita,et al.  Multidimensional replica-exchange method for free-energy calculations , 2000, cond-mat/0009120.

[11]  E. Tajkhorshid,et al.  Driven Metadynamics: Reconstructing Equilibrium Free Energies from Driven Adaptive-Bias Simulations , 2013, The journal of physical chemistry letters.

[12]  Marc-Etienne Rousseau,et al.  Cysteines 431 and 1074 Are Responsible for Inhibitory Disulfide Cross-linking between the Two Nucleotide-binding Sites in Human P-glycoprotein* , 2001, The Journal of Biological Chemistry.

[13]  H. Mchaourab,et al.  Conformational cycle of the ABC transporter MsbA in liposomes: detailed analysis using double electron-electron resonance spectroscopy. , 2009, Journal of molecular biology.

[14]  Mahmoud Moradi,et al.  Reaction path ensemble of the B–Z-DNA transition: a comprehensive atomistic study , 2012, Nucleic acids research.

[15]  D. Clarke,et al.  The Packing of the Transmembrane Segments of Human Multidrug Resistance P-glycoprotein Is Revealed by Disulfide Cross-linking Analysis* , 2000, The Journal of Biological Chemistry.

[16]  C. Roland,et al.  Conformations and free energy landscapes of polyproline peptides , 2009, Proceedings of the National Academy of Sciences.

[17]  C. Jarzynski Nonequilibrium Equality for Free Energy Differences , 1996, cond-mat/9610209.

[18]  Jue Chen,et al.  Crystal structure of the multidrug transporter P-glycoprotein from Caenorhabditis elegans , 2012 .

[19]  Berthold K. P. Horn,et al.  Closed-form solution of absolute orientation using unit quaternions , 1987 .

[20]  G. Hummer,et al.  Free energy reconstruction from nonequilibrium single-molecule pulling experiments , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[21]  N. Dokholyan,et al.  ATP hydrolysis at one of the two sites in ABC transporters initiates transport related conformational transitions. , 2011, Biochimica et biophysica acta.

[22]  K. Dill,et al.  Using quaternions to calculate RMSD , 2004, J. Comput. Chem..

[23]  Hassane S. Mchaourab,et al.  On the Origin of Large Flexibility of P-glycoprotein in the Inward-facing State* , 2013, The Journal of Biological Chemistry.

[24]  E. Tajkhorshid,et al.  Conformational coupling of the nucleotide-binding and the transmembrane domains in ABC transporters. , 2009, Biophysical journal.

[25]  K. Schulten,et al.  Molecular dynamics study of unbinding of the avidin-biotin complex. , 1997, Biophysical journal.

[26]  G. Crooks Path-ensemble averages in systems driven far from equilibrium , 1999, cond-mat/9908420.

[27]  C. Domene,et al.  Dynamics of a bacterial multidrug ABC transporter in the inward- and outward-facing conformations , 2012, Proceedings of the National Academy of Sciences.

[28]  C. Soares,et al.  Conformational changes induced by ATP‐hydrolysis in an ABC transporter: A molecular dynamics study of the Sav1866 exporter , 2011, Proteins.

[29]  Giacomo Fiorin,et al.  Using collective variables to drive molecular dynamics simulations , 2013 .

[30]  Geoffrey Chang,et al.  Flexibility in the ABC transporter MsbA: Alternating access with a twist , 2007, Proceedings of the National Academy of Sciences.

[31]  Jue Chen,et al.  Structure, Function, and Evolution of Bacterial ATP-Binding Cassette Systems , 2008, Microbiology and Molecular Biology Reviews.

[32]  Mark S. P. Sansom,et al.  Asymmetric Switching in a Homodimeric ABC Transporter: A Simulation Study , 2010, PLoS Comput. Biol..

[33]  Jürgen Schlitter,et al.  Targeted Molecular Dynamics Simulation of Conformational Change-Application to the T ↔ R Transition in Insulin , 1993 .

[34]  S. Wilkens,et al.  P-glycoprotein Retains Drug-stimulated ATPase Activity upon Covalent Linkage of the Two Nucleotide Binding Domains at Their C-terminal Ends* , 2011, The Journal of Biological Chemistry.

[35]  E. Pardon,et al.  Structures of P-glycoprotein reveal its conformational flexibility and an epitope on the nucleotide-binding domain , 2013, Proceedings of the National Academy of Sciences.

[36]  C. Higgins,et al.  A bacterial antibiotic-resistance gene that complements the human multidrug-resistance P-glycoprotein gene , 1998, Nature.

[37]  Emad Tajkhorshid,et al.  Dimer opening of the nucleotide binding domains of ABC transporters after ATP hydrolysis. , 2008, Biophysical journal.

[38]  K. Locher,et al.  Structure of the multidrug ABC transporter Sav1866 from Staphylococcus aureus in complex with AMP‐PNP , 2007, FEBS letters.

[39]  C. Roland,et al.  A classical molecular dynamics investigation of the free energy and structure of short polyproline conformers. , 2010, The Journal of chemical physics.

[40]  Feifei Gu,et al.  Thermal stability of apolipoprotein A-I in high-density lipoproteins by molecular dynamics. , 2009, Biophysical journal.

[41]  C. Klug,et al.  Characterization of the E506Q and H537A dysfunctional mutants in the E. coli ABC transporter MsbA. , 2011, Biochemistry.

[42]  Wonpil Im,et al.  Transmembrane helix assembly by window exchange umbrella sampling. , 2012, Physical review letters.

[43]  K. Schulten,et al.  Energetics of glycerol conduction through aquaglyceroporin GlpF , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[44]  A. M. George,et al.  Molecular-dynamics simulations of the ATP/apo state of a multidrug ATP-binding cassette transporter provide a structural and mechanistic basis for the asymmetric occluded state. , 2011, Biophysical journal.

[45]  Andrew J Ballard,et al.  Replica exchange with nonequilibrium switches , 2009, Proceedings of the National Academy of Sciences.

[46]  M. Prévost,et al.  Dynamics and structural changes induced by ATP binding in SAV1866, a bacterial ABC exporter. , 2010, The journal of physical chemistry. B.

[47]  R. Dawson,et al.  Structure and mechanism of ABC transporter proteins. , 2007, Current opinion in structural biology.