Exploring the Conformational Transitions of Biomolecular Systems Using a Simple Two-State Anisotropic Network Model
Identifieur interne : 000255 ( Pmc/Corpus ); précédent : 000254; suivant : 000256Exploring the Conformational Transitions of Biomolecular Systems Using a Simple Two-State Anisotropic Network Model
Auteurs : Avisek Das ; Mert Gur ; Mary Hongying Cheng ; Sunhwan Jo ; Ivet Bahar ; Benoît RouxSource :
- PLoS Computational Biology [ 1553-734X ] ; 2014.
Abstract
Biomolecular conformational transitions are essential to biological functions. Most experimental methods report on the long-lived functional states of biomolecules, but information about the transition pathways between these stable states is generally scarce. Such transitions involve short-lived conformational states that are difficult to detect experimentally. For this reason, computational methods are needed to produce plausible hypothetical transition pathways that can then be probed experimentally. Here we propose a simple and computationally efficient method, called
Url:
DOI: 10.1371/journal.pcbi.1003521
PubMed: 24699246
PubMed Central: 3974643
Links to Exploration step
PMC:3974643Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Exploring the Conformational Transitions of Biomolecular Systems Using a Simple Two-State Anisotropic Network Model</title>
<author><name sortKey="Das, Avisek" sort="Das, Avisek" uniqKey="Das A" first="Avisek" last="Das">Avisek Das</name>
<affiliation><nlm:aff id="aff1"><addr-line>Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Gur, Mert" sort="Gur, Mert" uniqKey="Gur M" first="Mert" last="Gur">Mert Gur</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Cheng, Mary Hongying" sort="Cheng, Mary Hongying" uniqKey="Cheng M" first="Mary Hongying" last="Cheng">Mary Hongying Cheng</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Jo, Sunhwan" sort="Jo, Sunhwan" uniqKey="Jo S" first="Sunhwan" last="Jo">Sunhwan Jo</name>
<affiliation><nlm:aff id="aff1"><addr-line>Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Bahar, Ivet" sort="Bahar, Ivet" uniqKey="Bahar I" first="Ivet" last="Bahar">Ivet Bahar</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Roux, Benoit" sort="Roux, Benoit" uniqKey="Roux B" first="Benoît" last="Roux">Benoît Roux</name>
<affiliation><nlm:aff id="aff1"><addr-line>Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">24699246</idno>
<idno type="pmc">3974643</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3974643</idno>
<idno type="RBID">PMC:3974643</idno>
<idno type="doi">10.1371/journal.pcbi.1003521</idno>
<date when="2014">2014</date>
<idno type="wicri:Area/Pmc/Corpus">000255</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000255</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Exploring the Conformational Transitions of Biomolecular Systems Using a Simple Two-State Anisotropic Network Model</title>
<author><name sortKey="Das, Avisek" sort="Das, Avisek" uniqKey="Das A" first="Avisek" last="Das">Avisek Das</name>
<affiliation><nlm:aff id="aff1"><addr-line>Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Gur, Mert" sort="Gur, Mert" uniqKey="Gur M" first="Mert" last="Gur">Mert Gur</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Cheng, Mary Hongying" sort="Cheng, Mary Hongying" uniqKey="Cheng M" first="Mary Hongying" last="Cheng">Mary Hongying Cheng</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Jo, Sunhwan" sort="Jo, Sunhwan" uniqKey="Jo S" first="Sunhwan" last="Jo">Sunhwan Jo</name>
<affiliation><nlm:aff id="aff1"><addr-line>Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Bahar, Ivet" sort="Bahar, Ivet" uniqKey="Bahar I" first="Ivet" last="Bahar">Ivet Bahar</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Roux, Benoit" sort="Roux, Benoit" uniqKey="Roux B" first="Benoît" last="Roux">Benoît Roux</name>
<affiliation><nlm:aff id="aff1"><addr-line>Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">PLoS Computational Biology</title>
<idno type="ISSN">1553-734X</idno>
<idno type="eISSN">1553-7358</idno>
<imprint><date when="2014">2014</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>Biomolecular conformational transitions are essential to biological functions. Most experimental methods report on the long-lived functional states of biomolecules, but information about the transition pathways between these stable states is generally scarce. Such transitions involve short-lived conformational states that are difficult to detect experimentally. For this reason, computational methods are needed to produce plausible hypothetical transition pathways that can then be probed experimentally. Here we propose a simple and computationally efficient method, called <italic>ANMPathway</italic>
, for constructing a physically reasonable pathway between two endpoints of a conformational transition. We adopt a coarse-grained representation of the protein and construct a two-state potential by combining two elastic network models (ENMs) representative of the experimental structures resolved for the endpoints. The two-state potential has a cusp hypersurface in the configuration space where the energies from both the ENMs are equal. We first search for the minimum energy structure on the cusp hypersurface and then treat it as the transition state. The continuous pathway is subsequently constructed by following the steepest descent energy minimization trajectories starting from the transition state on each side of the cusp hypersurface. Application to several systems of broad biological interest such as adenylate kinase, ATP-driven calcium pump SERCA, leucine transporter and glutamate transporter shows that <italic>ANMPathway</italic>
yields results in good agreement with those from other similar methods and with data obtained from all-atom molecular dynamics simulations, in support of the utility of this simple and efficient approach. Notably the method provides experimentally testable predictions, including the formation of non-native contacts during the transition which we were able to detect in two of the systems we studied. An open-access web server has been created to deliver <italic>ANMPathway</italic>
results.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Henzler Wildman, Ka" uniqKey="Henzler Wildman K">KA Henzler-Wildman</name>
</author>
<author><name sortKey="Thai, V" uniqKey="Thai V">V Thai</name>
</author>
<author><name sortKey="Lei, M" uniqKey="Lei M">M Lei</name>
</author>
<author><name sortKey="Ott, M" uniqKey="Ott M">M Ott</name>
</author>
<author><name sortKey="Wolf Watz, M" uniqKey="Wolf Watz M">M Wolf-Watz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lei, M" uniqKey="Lei M">M Lei</name>
</author>
<author><name sortKey="Velos, J" uniqKey="Velos J">J Velos</name>
</author>
<author><name sortKey="Gardino, A" uniqKey="Gardino A">A Gardino</name>
</author>
<author><name sortKey="Kivenson, A" uniqKey="Kivenson A">A Kivenson</name>
</author>
<author><name sortKey="Karplus, M" uniqKey="Karplus M">M Karplus</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lacroix, Jj" uniqKey="Lacroix J">JJ Lacroix</name>
</author>
<author><name sortKey="Pless, Sa" uniqKey="Pless S">SA Pless</name>
</author>
<author><name sortKey="Maragliano, L" uniqKey="Maragliano L">L Maragliano</name>
</author>
<author><name sortKey="Campos, Fv" uniqKey="Campos F">FV Campos</name>
</author>
<author><name sortKey="Galpin, Jd" uniqKey="Galpin J">JD Galpin</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Shaw, De" uniqKey="Shaw D">DE Shaw</name>
</author>
<author><name sortKey="Deneroff, Mm" uniqKey="Deneroff M">MM Deneroff</name>
</author>
<author><name sortKey="Dror, Ro" uniqKey="Dror R">RO Dror</name>
</author>
<author><name sortKey="Kuskin, Js" uniqKey="Kuskin J">JS Kuskin</name>
</author>
<author><name sortKey="Larson, Rh" uniqKey="Larson R">RH Larson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chandler, D" uniqKey="Chandler D">D Chandler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dellago, C" uniqKey="Dellago C">C Dellago</name>
</author>
<author><name sortKey="Bolhuis, Pg" uniqKey="Bolhuis P">PG Bolhuis</name>
</author>
<author><name sortKey="Csajka, Fs" uniqKey="Csajka F">FS Csajka</name>
</author>
<author><name sortKey="Chandler, D" uniqKey="Chandler D">D Chandler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bolhuis, Pg" uniqKey="Bolhuis P">PG Bolhuis</name>
</author>
<author><name sortKey="Chandler, D" uniqKey="Chandler D">D Chandler</name>
</author>
<author><name sortKey="Dellago, C" uniqKey="Dellago C">C Dellago</name>
</author>
<author><name sortKey="Geissler, Pl" uniqKey="Geissler P">PL Geissler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Du, Wn" uniqKey="Du W">WN Du</name>
</author>
<author><name sortKey="Marino, Ka" uniqKey="Marino K">KA Marino</name>
</author>
<author><name sortKey="Bolhuis, Pg" uniqKey="Bolhuis P">PG Bolhuis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Elber, R" uniqKey="Elber R">R Elber</name>
</author>
<author><name sortKey="Karplus, M" uniqKey="Karplus M">M Karplus</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Elber, R" uniqKey="Elber R">R Elber</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Maragliano, L" uniqKey="Maragliano L">L Maragliano</name>
</author>
<author><name sortKey="Fischer, A" uniqKey="Fischer A">A Fischer</name>
</author>
<author><name sortKey="Vanden Eijnden, E" uniqKey="Vanden Eijnden E">E Vanden-Eijnden</name>
</author>
<author><name sortKey="Ciccotti, G" uniqKey="Ciccotti G">G Ciccotti</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pan, Ac" uniqKey="Pan A">AC Pan</name>
</author>
<author><name sortKey="Sezer, D" uniqKey="Sezer D">D Sezer</name>
</author>
<author><name sortKey="Roux, B" uniqKey="Roux B">B Roux</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gan, W" uniqKey="Gan W">W Gan</name>
</author>
<author><name sortKey="Yang, S" uniqKey="Yang S">S Yang</name>
</author>
<author><name sortKey="Roux, B" uniqKey="Roux B">B Roux</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Vashisth, H" uniqKey="Vashisth H">H Vashisth</name>
</author>
<author><name sortKey="Maragliano, L" uniqKey="Maragliano L">L Maragliano</name>
</author>
<author><name sortKey="Abrams, Cf" uniqKey="Abrams C">CF Abrams</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Vashisth, H" uniqKey="Vashisth H">H Vashisth</name>
</author>
<author><name sortKey="Abrams, Cf" uniqKey="Abrams C">CF Abrams</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Matsunaga, Y" uniqKey="Matsunaga Y">Y Matsunaga</name>
</author>
<author><name sortKey="Fujisaki, H" uniqKey="Fujisaki H">H Fujisaki</name>
</author>
<author><name sortKey="Terada, T" uniqKey="Terada T">T Terada</name>
</author>
<author><name sortKey="Furuta, T" uniqKey="Furuta T">T Furuta</name>
</author>
<author><name sortKey="Moritsugu, K" uniqKey="Moritsugu K">K Moritsugu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stober, St" uniqKey="Stober S">ST Stober</name>
</author>
<author><name sortKey="Abrams, Cf" uniqKey="Abrams C">CF Abrams</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jo, S" uniqKey="Jo S">S Jo</name>
</author>
<author><name sortKey="Rui, H" uniqKey="Rui H">H Rui</name>
</author>
<author><name sortKey="Lim, Jb" uniqKey="Lim J">JB Lim</name>
</author>
<author><name sortKey="Klauda, Jb" uniqKey="Klauda J">JB Klauda</name>
</author>
<author><name sortKey="Im, W" uniqKey="Im W">W Im</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ovchinnikov, V" uniqKey="Ovchinnikov V">V Ovchinnikov</name>
</author>
<author><name sortKey="Karplus, M" uniqKey="Karplus M">M Karplus</name>
</author>
<author><name sortKey="Vanden Eijnden, E" uniqKey="Vanden Eijnden E">E Vanden-Eijnden</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kirmizialtin, S" uniqKey="Kirmizialtin S">S Kirmizialtin</name>
</author>
<author><name sortKey="Nguyen, V" uniqKey="Nguyen V">V Nguyen</name>
</author>
<author><name sortKey="Johnson, Ka" uniqKey="Johnson K">KA Johnson</name>
</author>
<author><name sortKey="Elber, R" uniqKey="Elber R">R Elber</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huber, G" uniqKey="Huber G">G Huber</name>
</author>
<author><name sortKey="Kim, S" uniqKey="Kim S">S Kim</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhang, Bw" uniqKey="Zhang B">BW Zhang</name>
</author>
<author><name sortKey="Jasnow, D" uniqKey="Jasnow D">D Jasnow</name>
</author>
<author><name sortKey="Zuckerman, Dm" uniqKey="Zuckerman D">DM Zuckerman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bhatt, D" uniqKey="Bhatt D">D Bhatt</name>
</author>
<author><name sortKey="Zhang, Bw" uniqKey="Zhang B">BW Zhang</name>
</author>
<author><name sortKey="Zuckerman, Dm" uniqKey="Zuckerman D">DM Zuckerman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bhatt, D" uniqKey="Bhatt D">D Bhatt</name>
</author>
<author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Woolf, Tb" uniqKey="Woolf T">TB Woolf</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zuckerman, Dm" uniqKey="Zuckerman D">DM Zuckerman</name>
</author>
<author><name sortKey="Woolf, Tb" uniqKey="Woolf T">TB Woolf</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Perilla, Jr" uniqKey="Perilla J">JR Perilla</name>
</author>
<author><name sortKey="Beckstein, O" uniqKey="Beckstein O">O Beckstein</name>
</author>
<author><name sortKey="Denning, Ej" uniqKey="Denning E">EJ Denning</name>
</author>
<author><name sortKey="Woolf, Tb" uniqKey="Woolf T">TB Woolf</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Grubmuller, H" uniqKey="Grubmuller H">H Grubmüller</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Laio, A" uniqKey="Laio A">A Laio</name>
</author>
<author><name sortKey="Parrinello, M" uniqKey="Parrinello M">M Parrinello</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hamelberg, D" uniqKey="Hamelberg D">D Hamelberg</name>
</author>
<author><name sortKey="Mongan, J" uniqKey="Mongan J">J Mongan</name>
</author>
<author><name sortKey="Mccammon, Ja" uniqKey="Mccammon J">JA McCammon</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Seeliger, D" uniqKey="Seeliger D">D Seeliger</name>
</author>
<author><name sortKey="Haas, J" uniqKey="Haas J">J Haas</name>
</author>
<author><name sortKey="De Groot, Bl" uniqKey="De Groot B">BL de Groot</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
<author><name sortKey="Atilgan, Ar" uniqKey="Atilgan A">AR Atilgan</name>
</author>
<author><name sortKey="Erman, B" uniqKey="Erman B">B Erman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Haliloglu, T" uniqKey="Haliloglu T">T Haliloglu</name>
</author>
<author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
<author><name sortKey="Erman, B" uniqKey="Erman B">B Erman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tirion, Mm" uniqKey="Tirion M">MM Tirion</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Atilgan, A" uniqKey="Atilgan A">A Atilgan</name>
</author>
<author><name sortKey="Durell, S" uniqKey="Durell S">S Durell</name>
</author>
<author><name sortKey="Jernigan, R" uniqKey="Jernigan R">R Jernigan</name>
</author>
<author><name sortKey="Demirel, M" uniqKey="Demirel M">M Demirel</name>
</author>
<author><name sortKey="Keskin, O" uniqKey="Keskin O">O Keskin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Thomas, A" uniqKey="Thomas A">A Thomas</name>
</author>
<author><name sortKey="Hinsen, K" uniqKey="Hinsen K">K Hinsen</name>
</author>
<author><name sortKey="Field, Mj" uniqKey="Field M">MJ Field</name>
</author>
<author><name sortKey="Perahia, D" uniqKey="Perahia D">D Perahia</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tama, F" uniqKey="Tama F">F Tama</name>
</author>
<author><name sortKey="Sanejouand, Yh" uniqKey="Sanejouand Y">YH Sanejouand</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Reuter, N" uniqKey="Reuter N">N Reuter</name>
</author>
<author><name sortKey="Hinsen, K" uniqKey="Hinsen K">K Hinsen</name>
</author>
<author><name sortKey="Lacapre, Jj" uniqKey="Lacapre J">JJ Lacapre</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
<author><name sortKey="Rader, A" uniqKey="Rader A">A Rader</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
<author><name sortKey="Lezon, Tr" uniqKey="Lezon T">TR Lezon</name>
</author>
<author><name sortKey="Yang, Lw" uniqKey="Yang L">LW Yang</name>
</author>
<author><name sortKey="Eyal, E" uniqKey="Eyal E">E Eyal</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
<author><name sortKey="Lezon, Tr" uniqKey="Lezon T">TR Lezon</name>
</author>
<author><name sortKey="Bakan, A" uniqKey="Bakan A">A Bakan</name>
</author>
<author><name sortKey="Shrivastava, Ih" uniqKey="Shrivastava I">IH Shrivastava</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gur, M" uniqKey="Gur M">M Gur</name>
</author>
<author><name sortKey="Zomot, E" uniqKey="Zomot E">E Zomot</name>
</author>
<author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kim, Mk" uniqKey="Kim M">MK Kim</name>
</author>
<author><name sortKey="Chirikjian, Gs" uniqKey="Chirikjian G">GS Chirikjian</name>
</author>
<author><name sortKey="Jernigan, Rl" uniqKey="Jernigan R">RL Jernigan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kim, Mk" uniqKey="Kim M">MK Kim</name>
</author>
<author><name sortKey="Jernigan, Rl" uniqKey="Jernigan R">RL Jernigan</name>
</author>
<author><name sortKey="Chirikjian, Gs" uniqKey="Chirikjian G">GS Chirikjian</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Schuyler, Ad" uniqKey="Schuyler A">AD Schuyler</name>
</author>
<author><name sortKey="Chirikjian, Gs" uniqKey="Chirikjian G">GS Chirikjian</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kim, Mk" uniqKey="Kim M">MK Kim</name>
</author>
<author><name sortKey="Jernigan, Rl" uniqKey="Jernigan R">RL Jernigan</name>
</author>
<author><name sortKey="Chirikjian, Gs" uniqKey="Chirikjian G">GS Chirikjian</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Schuyler, Ad" uniqKey="Schuyler A">AD Schuyler</name>
</author>
<author><name sortKey="Jernigan, Rl" uniqKey="Jernigan R">RL Jernigan</name>
</author>
<author><name sortKey="Qasba, Pk" uniqKey="Qasba P">PK Qasba</name>
</author>
<author><name sortKey="Ramakrishnan, B" uniqKey="Ramakrishnan B">B Ramakrishnan</name>
</author>
<author><name sortKey="Chirikjian, Gs" uniqKey="Chirikjian G">GS Chirikjian</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Miyashita, O" uniqKey="Miyashita O">O Miyashita</name>
</author>
<author><name sortKey="Onuchic, Jn" uniqKey="Onuchic J">JN Onuchic</name>
</author>
<author><name sortKey="Wolynes, Pg" uniqKey="Wolynes P">PG Wolynes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Maragakis, P" uniqKey="Maragakis P">P Maragakis</name>
</author>
<author><name sortKey="Karplus, M" uniqKey="Karplus M">M Karplus</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yang, Z" uniqKey="Yang Z">Z Yang</name>
</author>
<author><name sortKey="Majek, P" uniqKey="Majek P">P Majek</name>
</author>
<author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Best, Rb" uniqKey="Best R">RB Best</name>
</author>
<author><name sortKey="Chen, Yg" uniqKey="Chen Y">YG Chen</name>
</author>
<author><name sortKey="Hummer, G" uniqKey="Hummer G">G Hummer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zheng, W" uniqKey="Zheng W">W Zheng</name>
</author>
<author><name sortKey="Brooks, Br" uniqKey="Brooks B">BR Brooks</name>
</author>
<author><name sortKey="Hummer, G" uniqKey="Hummer G">G Hummer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhu, F" uniqKey="Zhu F">F Zhu</name>
</author>
<author><name sortKey="Hummer, G" uniqKey="Hummer G">G Hummer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yang, S" uniqKey="Yang S">S Yang</name>
</author>
<author><name sortKey="Roux, B" uniqKey="Roux B">B Roux</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="G, N" uniqKey="G N">N Gō</name>
</author>
<author><name sortKey="Abe, H" uniqKey="Abe H">H Abe</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Abe, H" uniqKey="Abe H">H Abe</name>
</author>
<author><name sortKey="G, N" uniqKey="G N">N Gō</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chu, Jw" uniqKey="Chu J">JW Chu</name>
</author>
<author><name sortKey="Voth, Ga" uniqKey="Voth G">GA Voth</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Franklin, J" uniqKey="Franklin J">J Franklin</name>
</author>
<author><name sortKey="Koehl, P" uniqKey="Koehl P">P Koehl</name>
</author>
<author><name sortKey="Doniach, S" uniqKey="Doniach S">S Doniach</name>
</author>
<author><name sortKey="Delarue, M" uniqKey="Delarue M">M Delarue</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Muller, C" uniqKey="Muller C">C Müller</name>
</author>
<author><name sortKey="Schlauderer, G" uniqKey="Schlauderer G">G Schlauderer</name>
</author>
<author><name sortKey="Reinstein, J" uniqKey="Reinstein J">J Reinstein</name>
</author>
<author><name sortKey="Schulz, G" uniqKey="Schulz G">G Schulz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Muller, Cw" uniqKey="Muller C">CW Müller</name>
</author>
<author><name sortKey="Schulz, Ge" uniqKey="Schulz G">GE Schulz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gur, M" uniqKey="Gur M">M Gur</name>
</author>
<author><name sortKey="Madura, J" uniqKey="Madura J">J Madura</name>
</author>
<author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kantarci Carsibasi, N" uniqKey="Kantarci Carsibasi N">N Kantarci-Carsibasi</name>
</author>
<author><name sortKey="Haliloglu, T" uniqKey="Haliloglu T">T Haliloglu</name>
</author>
<author><name sortKey="Doruker, P" uniqKey="Doruker P">P Doruker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Arora, K" uniqKey="Arora K">K Arora</name>
</author>
<author><name sortKey="Brooks, Cl" uniqKey="Brooks C">CL Brooks</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Daily, Md" uniqKey="Daily M">MD Daily</name>
</author>
<author><name sortKey="Jr, Gnp" uniqKey="Jr G">GNP Jr</name>
</author>
<author><name sortKey="Cui, Q" uniqKey="Cui Q">Q Cui</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Daily, Md" uniqKey="Daily M">MD Daily</name>
</author>
<author><name sortKey="Phillips, Gn" uniqKey="Phillips G">GN Phillips</name>
</author>
<author><name sortKey="Cui, Q" uniqKey="Cui Q">Q Cui</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Whitford, Pc" uniqKey="Whitford P">PC Whitford</name>
</author>
<author><name sortKey="Miyashita, O" uniqKey="Miyashita O">O Miyashita</name>
</author>
<author><name sortKey="Levy, Y" uniqKey="Levy Y">Y Levy</name>
</author>
<author><name sortKey="Onuchic, Jn" uniqKey="Onuchic J">JN Onuchic</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Whitford, Pc" uniqKey="Whitford P">PC Whitford</name>
</author>
<author><name sortKey="Gosavi, S" uniqKey="Gosavi S">S Gosavi</name>
</author>
<author><name sortKey="Onuchic, Jn" uniqKey="Onuchic J">JN Onuchic</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Brokaw, Jb" uniqKey="Brokaw J">JB Brokaw</name>
</author>
<author><name sortKey="Chu, Jw" uniqKey="Chu J">JW Chu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Beckstein, O" uniqKey="Beckstein O">O Beckstein</name>
</author>
<author><name sortKey="Denning, Ej" uniqKey="Denning E">EJ Denning</name>
</author>
<author><name sortKey="Perilla, Jr" uniqKey="Perilla J">JR Perilla</name>
</author>
<author><name sortKey="Woolf, Tb" uniqKey="Woolf T">TB Woolf</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yamashita, A" uniqKey="Yamashita A">A Yamashita</name>
</author>
<author><name sortKey="Singh, Sk" uniqKey="Singh S">SK Singh</name>
</author>
<author><name sortKey="Kawate, T" uniqKey="Kawate T">T Kawate</name>
</author>
<author><name sortKey="Jin, Y" uniqKey="Jin Y">Y Jin</name>
</author>
<author><name sortKey="Gouaux, E" uniqKey="Gouaux E">E Gouaux</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Krishnamurthy, H" uniqKey="Krishnamurthy H">H Krishnamurthy</name>
</author>
<author><name sortKey="Gouaux, E" uniqKey="Gouaux E">E Gouaux</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cheng, Mh" uniqKey="Cheng M">MH Cheng</name>
</author>
<author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yernool, D" uniqKey="Yernool D">D Yernool</name>
</author>
<author><name sortKey="Boudker, O" uniqKey="Boudker O">O Boudker</name>
</author>
<author><name sortKey="Jin, Y" uniqKey="Jin Y">Y Jin</name>
</author>
<author><name sortKey="Gouaux, E" uniqKey="Gouaux E">E Gouaux</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Reyes, N" uniqKey="Reyes N">N Reyes</name>
</author>
<author><name sortKey="Ginter, C" uniqKey="Ginter C">C Ginter</name>
</author>
<author><name sortKey="Boudker, O" uniqKey="Boudker O">O Boudker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Verdon, G" uniqKey="Verdon G">G Verdon</name>
</author>
<author><name sortKey="Boudker, O" uniqKey="Boudker O">O Boudker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jiang, J" uniqKey="Jiang J">J Jiang</name>
</author>
<author><name sortKey="Shrivastava, Ih" uniqKey="Shrivastava I">IH Shrivastava</name>
</author>
<author><name sortKey="Watts, Sd" uniqKey="Watts S">SD Watts</name>
</author>
<author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
<author><name sortKey="Amara, Sg" uniqKey="Amara S">SG Amara</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zomot, E" uniqKey="Zomot E">E Zomot</name>
</author>
<author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shrivastava, Ih" uniqKey="Shrivastava I">IH Shrivastava</name>
</author>
<author><name sortKey="Jiang, J" uniqKey="Jiang J">J Jiang</name>
</author>
<author><name sortKey="Amara, Sg" uniqKey="Amara S">SG Amara</name>
</author>
<author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lezon, Tr" uniqKey="Lezon T">TR Lezon</name>
</author>
<author><name sortKey="Bahar, I" uniqKey="Bahar I">I Bahar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Toyoshima, C" uniqKey="Toyoshima C">C Toyoshima</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Toyoshima, C" uniqKey="Toyoshima C">C Toyoshima</name>
</author>
<author><name sortKey="Nakasako, M" uniqKey="Nakasako M">M Nakasako</name>
</author>
<author><name sortKey="Nomura, H" uniqKey="Nomura H">H Nomura</name>
</author>
<author><name sortKey="Ogawa, H" uniqKey="Ogawa H">H Ogawa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sorensen, Tlm" uniqKey="Sorensen T">TLM Sorensen</name>
</author>
<author><name sortKey="Moller, Jv" uniqKey="Moller J">JV Moller</name>
</author>
<author><name sortKey="Nissen, P" uniqKey="Nissen P">P Nissen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Toyoshima, C" uniqKey="Toyoshima C">C Toyoshima</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Espinoza Fonseca, Lm" uniqKey="Espinoza Fonseca L">LM Espinoza-Fonseca</name>
</author>
<author><name sortKey="Thomas, Dd" uniqKey="Thomas D">DD Thomas</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Winters, Dl" uniqKey="Winters D">DL Winters</name>
</author>
<author><name sortKey="Autry, Jm" uniqKey="Autry J">JM Autry</name>
</author>
<author><name sortKey="Svensson, B" uniqKey="Svensson B">B Svensson</name>
</author>
<author><name sortKey="Thomas, Dd" uniqKey="Thomas D">DD Thomas</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tekpinar, M" uniqKey="Tekpinar M">M Tekpinar</name>
</author>
<author><name sortKey="Zheng, W" uniqKey="Zheng W">W Zheng</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Krger, Dm" uniqKey="Krger D">DM Krger</name>
</author>
<author><name sortKey="Ahmed, A" uniqKey="Ahmed A">A Ahmed</name>
</author>
<author><name sortKey="Gohlke, H" uniqKey="Gohlke H">H Gohlke</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sfriso, P" uniqKey="Sfriso P">P Sfriso</name>
</author>
<author><name sortKey="Emperador, A" uniqKey="Emperador A">A Emperador</name>
</author>
<author><name sortKey="Orellana, L" uniqKey="Orellana L">L Orellana</name>
</author>
<author><name sortKey="Hospital, A" uniqKey="Hospital A">A Hospital</name>
</author>
<author><name sortKey="Gelpi, Jl" uniqKey="Gelpi J">JL Gelpí</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="E, W" uniqKey="E W">W E</name>
</author>
<author><name sortKey="Ren, W" uniqKey="Ren W">W Ren</name>
</author>
<author><name sortKey="Vanden Eijnden, E" uniqKey="Vanden Eijnden E">E Vanden-Eijnden</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Adelman, Jl" uniqKey="Adelman J">JL Adelman</name>
</author>
<author><name sortKey="Dale, Al" uniqKey="Dale A">AL Dale</name>
</author>
<author><name sortKey="Zwier, Mc" uniqKey="Zwier M">MC Zwier</name>
</author>
<author><name sortKey="Bhatt, D" uniqKey="Bhatt D">D Bhatt</name>
</author>
<author><name sortKey="Chong, Lt" uniqKey="Chong L">LT Chong</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS Comput Biol</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS Comput. Biol</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">ploscomp</journal-id>
<journal-title-group><journal-title>PLoS Computational Biology</journal-title>
</journal-title-group>
<issn pub-type="ppub">1553-734X</issn>
<issn pub-type="epub">1553-7358</issn>
<publisher><publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">24699246</article-id>
<article-id pub-id-type="pmc">3974643</article-id>
<article-id pub-id-type="publisher-id">PCOMPBIOL-D-13-01573</article-id>
<article-id pub-id-type="doi">10.1371/journal.pcbi.1003521</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
<subj-group subj-group-type="Discipline-v2"><subject>Biology and Life Sciences</subject>
<subj-group><subject>Biophysics</subject>
</subj-group>
<subj-group><subject>Computational Biology</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group><article-title>Exploring the Conformational Transitions of Biomolecular Systems Using a Simple Two-State Anisotropic Network Model</article-title>
<alt-title alt-title-type="running-head">Conformational Transitions of Biomolecular Systems</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author" equal-contrib="yes"><name><surname>Das</surname>
<given-names>Avisek</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author" equal-contrib="yes"><name><surname>Gur</surname>
<given-names>Mert</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author" equal-contrib="yes"><name><surname>Cheng</surname>
<given-names>Mary Hongying</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Jo</surname>
<given-names>Sunhwan</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Bahar</surname>
<given-names>Ivet</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Roux</surname>
<given-names>Benoît</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
<xref ref-type="corresp" rid="cor1"><sup>*</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label>
<addr-line>Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America</addr-line>
</aff>
<aff id="aff2"><label>2</label>
<addr-line>Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America</addr-line>
</aff>
<contrib-group><contrib contrib-type="editor"><name><surname>de Groot</surname>
<given-names>Bert L.</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>Max Planck Institute for Biophysical Chemistry, Germany</addr-line>
</aff>
<author-notes><corresp id="cor1">* E-mail: <email>roux@uchicago.edu</email>
</corresp>
<fn fn-type="conflict"><p>The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con"><p>Conceived and designed the experiments: AD IB BR. Performed the experiments: AD MG MHC. Analyzed the data: AD MG MHC SJ IB BR. Wrote the paper: AD MG MHC SJ IB BR.</p>
</fn>
</author-notes>
<pub-date pub-type="collection"><month>4</month>
<year>2014</year>
</pub-date>
<pub-date pub-type="epub"><day>3</day>
<month>4</month>
<year>2014</year>
</pub-date>
<volume>10</volume>
<issue>4</issue>
<elocation-id>e1003521</elocation-id>
<history><date date-type="received"><day>3</day>
<month>9</month>
<year>2013</year>
</date>
<date date-type="accepted"><day>1</day>
<month>2</month>
<year>2014</year>
</date>
</history>
<permissions><copyright-statement>© 2014 Das et al</copyright-statement>
<copyright-year>2014</copyright-year>
<copyright-holder>Das et al</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are properly credited.</license-p>
</license>
</permissions>
<abstract><p>Biomolecular conformational transitions are essential to biological functions. Most experimental methods report on the long-lived functional states of biomolecules, but information about the transition pathways between these stable states is generally scarce. Such transitions involve short-lived conformational states that are difficult to detect experimentally. For this reason, computational methods are needed to produce plausible hypothetical transition pathways that can then be probed experimentally. Here we propose a simple and computationally efficient method, called <italic>ANMPathway</italic>
, for constructing a physically reasonable pathway between two endpoints of a conformational transition. We adopt a coarse-grained representation of the protein and construct a two-state potential by combining two elastic network models (ENMs) representative of the experimental structures resolved for the endpoints. The two-state potential has a cusp hypersurface in the configuration space where the energies from both the ENMs are equal. We first search for the minimum energy structure on the cusp hypersurface and then treat it as the transition state. The continuous pathway is subsequently constructed by following the steepest descent energy minimization trajectories starting from the transition state on each side of the cusp hypersurface. Application to several systems of broad biological interest such as adenylate kinase, ATP-driven calcium pump SERCA, leucine transporter and glutamate transporter shows that <italic>ANMPathway</italic>
yields results in good agreement with those from other similar methods and with data obtained from all-atom molecular dynamics simulations, in support of the utility of this simple and efficient approach. Notably the method provides experimentally testable predictions, including the formation of non-native contacts during the transition which we were able to detect in two of the systems we studied. An open-access web server has been created to deliver <italic>ANMPathway</italic>
results.</p>
</abstract>
<abstract abstract-type="summary"><title>Author Summary</title>
<p>Many biomolecules are like tiny molecular machines that need to change their shapes and visit many states to perform their biological functions. For a complete molecular understanding of a biological process, one needs to have information on the relevant stable states of the system in question, as well as the pathways by which the system travels from one state to another. We report here an efficient computational method that uses the knowledge of experimental structures of a pair of stable states in order to construct an energetically favoravle pathway between them. We adopt a simple representation of the molecular system by replacing the atoms with beads connected by springs and constructing an energy function with two minima around the end-states. We searched for the structure with highest energy that the system is most likely to visit during the transition and created two paths starting from this structure and proceeding toward the end-states. The combined result of these two paths is the minimum energy pathway between the two stable states. We apply this method to study important structural changes in one enzyme and three large proteins that transport small molecules and ions across the cell membrane.</p>
</abstract>
<funding-group><funding-statement>This work was carried out in the context of the Membrane Protein Structural Dynamics Consortium which is funded by grant U54-GM087519 from the National Institute of Health (NIH). Additional support from NIH through grants (to I.B.) R01-GM099738, R01GM086238 and P41GM103712, and computing award from the NSF TeraGrid resources (TG-MCB130006) is gratefully acknowledged. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement>
</funding-group>
<counts><page-count count="17"></page-count>
</counts>
</article-meta>
</front>
<body><sec id="s1"><title>Introduction</title>
<p>Complex macromolecular systems such as enzymes, channels, transporters and pumps need to change their shapes and visit many conformational states in order to perform their functions. Experimental data from functional, biochemical, spectroscopic and structural techniques often inform us on the long-lived stable functional states of macromolecular systems. Accordingly, the average structures of thousands of important biomolecules have been determined using X-ray crystallography or NMR. For many well-studied systems, hundreds of structures have been resolved in the presence of different ligands, or under different conditions or functional states. In contrast, for most systems, little or no experimental data are often available on the intermediate structures along the conformational transition pathway associated with a function. To understand the molecular mechanism of a specific biological process, one needs to go beyond the static information and determine how macromolecules change their conformations as a function of time. In practice, however, obtaining direct structural data about a transition pathway is exceedingly difficult, because intermediate conformations are transient and usually short-lived compared to the timescale of the whole process.</p>
<p>Computational methods can help generate physically plausible pathways for conformational transitions, which can then serve as “hypotheses” to be tested and refined experimentally <xref rid="pcbi.1003521-HenzlerWildman1" ref-type="bibr">[1]</xref>
–<xref rid="pcbi.1003521-Lacroix1" ref-type="bibr">[3]</xref>
. The relevance of any <italic>in silico</italic>
pathway lies in its ability to predict the occurrence of intermediates, which can sometimes be detected using X-ray crystallography <xref rid="pcbi.1003521-HenzlerWildman1" ref-type="bibr">[1]</xref>
or indirectly inferred by perturbing the system via site-directed mutagenesis <xref rid="pcbi.1003521-Lacroix1" ref-type="bibr">[3]</xref>
. Computation and experimental validation thus offers a powerful combination to study the mechanisms of complex biomolecular events.</p>
<p>All-atom molecular dynamics (MD) simulation, arguably, provides the most realistic representation of biomolecular dynamics <xref rid="pcbi.1003521-Allen1" ref-type="bibr">[4]</xref>
, <xref rid="pcbi.1003521-Frenkel1" ref-type="bibr">[5]</xref>
. If one could simulate the system of interest for sufficiently long time-scales then the trajectory could provide the information required to understand a conformational transition (albeit based on a virtual model). However, brute-force MD is often impractical since most large-scale conformational changes take place over timescales ranging from milliseconds to seconds, which are far beyond the reach of the most powerful supercomputers. Special-purpose hardware and software like Anton <xref rid="pcbi.1003521-Shaw1" ref-type="bibr">[6]</xref>
are pushing the limits of current molecular simulations; however, they still fall short of accessing the relevant time-scales for the cooperative structural changes of large biological systems. Statistical mechanical methods have also been developed specifically to simulate rare dynamical events <xref rid="pcbi.1003521-Chandler1" ref-type="bibr">[7]</xref>
–<xref rid="pcbi.1003521-Du1" ref-type="bibr">[10]</xref>
, though their application to study large-scale conformational transitions in biological macromolecules remains challenging.</p>
<p>One alternative strategy has been to formulate the problem of the conformational reaction pathway as a “chain-of-state”, i.e., a sequence of configurations representing the progress of the system between two known end-states in the multi-dimensional conformation space <xref rid="pcbi.1003521-Elber1" ref-type="bibr">[11]</xref>
–<xref rid="pcbi.1003521-Pan1" ref-type="bibr">[15]</xref>
. For example, the so-called “string method” based on all-atom MD simulations has been employed successfully to study functionally important conformational transitions in a variety of biological systems, including Src kinases <xref rid="pcbi.1003521-Gan1" ref-type="bibr">[16]</xref>
, insulin receptor kinase <xref rid="pcbi.1003521-Vashisth1" ref-type="bibr">[17]</xref>
, <xref rid="pcbi.1003521-Vashisth2" ref-type="bibr">[18]</xref>
, adenylate kinase <xref rid="pcbi.1003521-Matsunaga1" ref-type="bibr">[19]</xref>
, amyloidogenic isomerization of 2-microglobulin <xref rid="pcbi.1003521-Stober1" ref-type="bibr">[20]</xref>
, cholesterol flip in membranes <xref rid="pcbi.1003521-Jo1" ref-type="bibr">[21]</xref>
, myosin VI <xref rid="pcbi.1003521-Ovchinnikov1" ref-type="bibr">[22]</xref>
, DNA polymerase <xref rid="pcbi.1003521-Kirmizialtin1" ref-type="bibr">[23]</xref>
, and voltage-gated K<sup>+</sup>
channels <xref rid="pcbi.1003521-Lacroix1" ref-type="bibr">[3]</xref>
. Other notable methods that seek to shed some light on the important intermediate structures between experimentally known intermediates include the weighted ensemble method <xref rid="pcbi.1003521-Huber1" ref-type="bibr">[24]</xref>
–<xref rid="pcbi.1003521-Bhatt2" ref-type="bibr">[27]</xref>
and dynamic importance sampling <xref rid="pcbi.1003521-Woolf1" ref-type="bibr">[28]</xref>
–<xref rid="pcbi.1003521-Perilla1" ref-type="bibr">[30]</xref>
. Several enhanced sampling methods such as conformational flooding <xref rid="pcbi.1003521-Grubmller1" ref-type="bibr">[31]</xref>
, metadynamics <xref rid="pcbi.1003521-Laio1" ref-type="bibr">[32]</xref>
and accelerated molecular dynamics <xref rid="pcbi.1003521-Hamelberg1" ref-type="bibr">[33]</xref>
have also been used for similar purpose even though these methods are not designed for searching for transition pathways between two known endpoints. A very different approach based on shapes of biomolecules rather than detailed energetics of the system has been used to search for transition pathways in the tCONCOORD method of Seeliger <italic>et al.</italic>
<xref rid="pcbi.1003521-Seeliger1" ref-type="bibr">[34]</xref>
. However, despite these promising advances, the investigation of large-scale transitions of multimeric systems at atomic details remains prohibitively expensive. It is, therefore, desirable to dispose of simpler models and computationally efficient methods to generate a putative pathway with qualitatively reasonable features that can be tested by experiments.</p>
<p>A straightforward way to reduce the computational cost is to simplify the atomistic system by constructing a coarse-grained (CG) model and adopting a simple potential function that uses knowledge of resolved structures. We adopt here a broadly used/tested structure-based CG model, the elastic network model (ENM) <xref rid="pcbi.1003521-Bahar1" ref-type="bibr">[35]</xref>
–<xref rid="pcbi.1003521-Tirion1" ref-type="bibr">[37]</xref>
, a powerful example of which is the anisotropic network model (ANM) <xref rid="pcbi.1003521-Atilgan1" ref-type="bibr">[38]</xref>
.</p>
<p>In the ANM, the protein is represented by a set of CG sites (nodes) placed at the positions of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e001.jpg"></inline-graphic>
</inline-formula>
atoms of all the residues and the energy function is a pairwise additive harmonic potential where each site interacts with all the sites within a cut-off distance. ENMs are often used in conjunction with normal mode analysis (NMA) where one diagonalizes the Hessian matrix of the potential constructed around an experimental structure and studies the deformation of the system along the low frequency normal modes. The simplified potential function in the ANM presents the advantage of yielding an analytical expression for the Hessian, directly expressed in terms of the known structure coordinates <xref rid="pcbi.1003521-Atilgan1" ref-type="bibr">[38]</xref>
, which is readily decomposed to obtain the ANM (normal) modes. It has been found for a wide variety of large biomolecular systems that collective motions relevant to function occur along the low energy normal modes of motions predicted by ENMs <xref rid="pcbi.1003521-Thomas1" ref-type="bibr">[39]</xref>
–<xref rid="pcbi.1003521-Gur1" ref-type="bibr">[45]</xref>
, suggesting that native contact topology accounted for by the network model is a major determinant of accessible modes of function.</p>
<p>Even though ENMs coupled with NMA have been successful in providing insights into important conformational transitions, they explore, by definition, the neighborhood of a given energy minimum and as such they are not adequate for constructing a transition pathway between two endpoints (minima on conformational energy landscape). However, ENM-based approaches have been very influential in the development of a series of methods that aim at providing plausible intermediate structures along a transition. One of the early studies along these lines is that of Jernigan, Chirikjian and coworkers <xref rid="pcbi.1003521-Kim1" ref-type="bibr">[46]</xref>
, <xref rid="pcbi.1003521-Kim2" ref-type="bibr">[47]</xref>
who have used an interpolation technique with distance constraints to avoid steric clashes. They also showed that normal mode calculations could be accelerated by dividing the system into rigid clusters connected by elastic springs <xref rid="pcbi.1003521-Schuyler1" ref-type="bibr">[48]</xref>
, and employed cluster-NMA for constructing pathways by successively creating new structures from an end-state <xref rid="pcbi.1003521-Kim3" ref-type="bibr">[49]</xref>
, <xref rid="pcbi.1003521-Schuyler2" ref-type="bibr">[50]</xref>
. Miyashita <italic>et al.</italic>
<xref rid="pcbi.1003521-Miyashita1" ref-type="bibr">[51]</xref>
started from one stable state, performed successive normal mode calculations and for each new set of normal modes used a small subset based on the overlap with the other end structure to create an intermediate structure. In the plastic network model (PNM), Maragakis and Karplus <xref rid="pcbi.1003521-Maragakis1" ref-type="bibr">[52]</xref>
constructed a two-state elastic network potential by mixing two ENMs, one for each end-state, and then the pathway was constructed in two steps: identification of a saddle point and two steepest descent minimizations. Yang <italic>et al.</italic>
<xref rid="pcbi.1003521-Yang1" ref-type="bibr">[53]</xref>
used the same two-state potential to start from both end structures and used well-defined criteria for recruiting small subsets of normal modes to create a series of intermediate conformers via an <italic>a</italic>
daptive ANM (<italic>a</italic>
ANM) methodology until the two intermediates merged within a predefined root-mean-square-deviation (RMSD). Hummer and co-workers <xref rid="pcbi.1003521-Best1" ref-type="bibr">[54]</xref>
–<xref rid="pcbi.1003521-Zhu1" ref-type="bibr">[56]</xref>
also constructed a two-state potential by mixing two ENM surfaces using an exponential mixing rule and constructed the pathway on this surface using saddle point search. Yang and Roux <xref rid="pcbi.1003521-Yang2" ref-type="bibr">[57]</xref>
have used a two-state G<inline-formula><inline-graphic xlink:href="pcbi.1003521.e002.jpg"></inline-graphic>
</inline-formula>
model <xref rid="pcbi.1003521-G1" ref-type="bibr">[58]</xref>
, <xref rid="pcbi.1003521-Abe1" ref-type="bibr">[59]</xref>
and extensive CG simulations in conjunction with clustering methods to investigate pathways in conformational transition of Src-kinase. Chu and Voth <xref rid="pcbi.1003521-Chu1" ref-type="bibr">[60]</xref>
used a more complicated two-state potential by representing each pairwise interaction as a double well and used a saddle point search algorithm to construct the pathway. Their double-well network model has more frustration than a two-state elastic network model and captures complexity of a transition that are not present in models with smoother potential energy functions. Franklin <italic>et al.</italic>
<xref rid="pcbi.1003521-Franklin1" ref-type="bibr">[61]</xref>
used two ENM surfaces in an entirely different way to construct a pathway method. In their MinActionPath method, they developed an algorithm based on the minimization of the Onsager-Machlup action to construct the path with minimum resistance between two stable states. Even though the problem of construction of pathway between two stable states described by simple CG models has attracted a lot of attention, there is still enough room and need for development of new methods that can address the scalability problem in particular, and help efficiently calculate pathways for large systems.</p>
<p>We propose a simple and efficient method, called <italic>ANMPathway</italic>
, and apply the method to understand conformational transitions of several important globular and membrane proteins. We adopt a simple ENM representation for each of the end-states, which accounts for the topology of inter-residue contacts in the structure. We construct a very simple two-state potential by mixing these two ENMs. Our potential has a cusp hypersurface where the energies from both the ENMs are same. We search for a minimum energy structure on the cusp hypersurface and treat it as the transition state. We then start from the transition state and perform two separate steepest descent minimizations to connect the end-states. Conformers collected from two steepest descent paths along with the transition state provide a pathway. Even though the existence of a cusp hypersurface in our potential is somewhat unphysical, we demonstrate, by way of applications to several systems (adenylate kinase (AK), ATP driven calcium pump SERCA, leucine transporter (LeuT) and glutamate transporter (<inline-formula><inline-graphic xlink:href="pcbi.1003521.e003.jpg"></inline-graphic>
</inline-formula>
)), that <italic>ANMPathway</italic>
gives physically meaningful pathways and helps generate experimentally testable hypotheses.</p>
</sec>
<sec sec-type="methods" id="s2"><title>Methods</title>
<p>The goal of the <italic>ANMPathway</italic>
method is to construct a transition pathway between two end-states of a conformational transition. As in the string method, the pathway is represented by a chain of equidistant states (conformers/images) <xref rid="pcbi.1003521-Elber1" ref-type="bibr">[11]</xref>
–<xref rid="pcbi.1003521-Maragliano1" ref-type="bibr">[14]</xref>
. The macromolecular structure is described by a CG model where interaction sites are placed at the positions of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e004.jpg"></inline-graphic>
</inline-formula>
atoms, which serve as the set of collective variables for the string. Conceptually, the string is the minimum free energy pathway on the potential of mean force (PMF) of the system with respect to those collective variables <xref rid="pcbi.1003521-Maragliano1" ref-type="bibr">[14]</xref>
. Assuming a Euclidian metric in the cartesian space of the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e005.jpg"></inline-graphic>
</inline-formula>
atoms, the equidistant conditions implies that neighboring images along the string are separated by a fixed RMSD. However, in practice, the ENM energy function employed in the CG model is a knowledge-based construct that is only applicable near the experimental structure used to construct the model. For describing a conformational transition between two stable states, we need an approximation to the PMF that is applicable for large distortions from the experimental structures.</p>
<p>In the presence of structural data on the end-states, it is reasonable to construct an effective energy function with two minima centered around the endpoints of the transition. One obvious route to such two-state potentials involves creating two separate energy surfaces that are defined around each of the end-states and then combine these surfaces by an empirical rule. We have adopted this strategy and used two ANMs <xref rid="pcbi.1003521-Atilgan1" ref-type="bibr">[38]</xref>
and a very simple mixing rule to construct an energy function with two minima.</p>
<p>For a protein with <italic>N</italic>
residues, the configuration of the system is denoted by a 3<italic>N</italic>
dimensional vector <inline-formula><inline-graphic xlink:href="pcbi.1003521.e006.jpg"></inline-graphic>
</inline-formula>
where, <inline-formula><inline-graphic xlink:href="pcbi.1003521.e007.jpg"></inline-graphic>
</inline-formula>
is a three-dimensional vector giving the position of the <italic>i</italic>
th site (<inline-formula><inline-graphic xlink:href="pcbi.1003521.e008.jpg"></inline-graphic>
</inline-formula>
atom of the <italic>i</italic>
th residue). ANM is an elastic network model defined around an experimental structure (e.g. crystal or NMR structure) <inline-formula><inline-graphic xlink:href="pcbi.1003521.e009.jpg"></inline-graphic>
</inline-formula>
with the following energy function<disp-formula id="pcbi.1003521.e010"><graphic xlink:href="pcbi.1003521.e010.jpg" position="anchor" orientation="portrait"></graphic>
<label>(1)</label>
</disp-formula>
Here <inline-formula><inline-graphic xlink:href="pcbi.1003521.e011.jpg"></inline-graphic>
</inline-formula>
is the distance between nodes <italic>i</italic>
and <italic>j</italic>
, <italic>k</italic>
is the uniform force constant, <italic>C<sub>ij</sub>
</italic>
is an element of the contact matrix defined by<disp-formula id="pcbi.1003521.e012"><graphic xlink:href="pcbi.1003521.e012.jpg" position="anchor" orientation="portrait"></graphic>
<label>(2)</label>
</disp-formula>
<disp-formula id="pcbi.1003521.e013"><graphic xlink:href="pcbi.1003521.e013.jpg" position="anchor" orientation="portrait"></graphic>
<label>(3)</label>
</disp-formula>
<inline-formula><inline-graphic xlink:href="pcbi.1003521.e014.jpg"></inline-graphic>
</inline-formula>
is the cut-off distance and <inline-formula><inline-graphic xlink:href="pcbi.1003521.e015.jpg"></inline-graphic>
</inline-formula>
is the energy of the system at the reference state. The advantage of including the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e016.jpg"></inline-graphic>
</inline-formula>
term is that it allows us to create energy difference between the end-states when more than one ENM are included in the model. In order to construct the potential function we first define two ANM energy functions, <inline-formula><inline-graphic xlink:href="pcbi.1003521.e017.jpg"></inline-graphic>
</inline-formula>
and <inline-formula><inline-graphic xlink:href="pcbi.1003521.e018.jpg"></inline-graphic>
</inline-formula>
, centered around the end structures <inline-formula><inline-graphic xlink:href="pcbi.1003521.e019.jpg"></inline-graphic>
</inline-formula>
and <inline-formula><inline-graphic xlink:href="pcbi.1003521.e020.jpg"></inline-graphic>
</inline-formula>
and combine them by the following mixing rule,<disp-formula id="pcbi.1003521.e021"><graphic xlink:href="pcbi.1003521.e021.jpg" position="anchor" orientation="portrait"></graphic>
<label>(4)</label>
</disp-formula>
The energy difference between the end-states of this two-state ENM is <inline-formula><inline-graphic xlink:href="pcbi.1003521.e022.jpg"></inline-graphic>
</inline-formula>
.</p>
<p>The two-state potential based on <xref ref-type="disp-formula" rid="pcbi.1003521.e021">Eq. (4)</xref>
has a cusp hypersurface in the 3<italic>N</italic>
-dimensional configuration space. Even though potential energy functions developed for real systems are differentiable everywhere, we will show that the simple two-state potential is a reasonable first approximation and is capable of capturing important qualitative features of the conformational transition in question. Both the ANM energy functions are 3<italic>N</italic>
-dimensional harmonic surfaces and the hypersurface where they intersect (i.e. where the energies from both ANM surfaces are same) is another harmonic surface of dimension <inline-formula><inline-graphic xlink:href="pcbi.1003521.e023.jpg"></inline-graphic>
</inline-formula>
. We define the transition state as the minimum energy structure on the cusp hypersurface. Given a sequence of conformers that linearly interpolates the Cartesian distance between two conformers that reside on the opposite sides of the cusp hypersurface, it is possible to identify a conformer that has equal energies, within a tolerance, from both the ANM surfaces. This conformer, by construction, resides on the cusp hypersurface. This simple observation allows us to devise an algorithm to search for the energy minimum on the cusp hypersurface i.e. the transition state. Once we have identified the transition state, we can start from there and slide down the harmonic surfaces until we reach the endpoints, by performing two separate steepest descent minimizations. In the end, we collect all the conformers in proper order to construct the transition pathway. The pathway obtained by <italic>ANMPathway</italic>
can be regarded as the minimum energy path between the end structures since it is the combination of two steepest descent paths on two surfaces joined at the transition state which is the minimum energy conformer on the cusp hypersurface. A detailed description of the algorithm is given below.</p>
<list list-type="order"><list-item><p>Two end structures are represented by the positions of their <inline-formula><inline-graphic xlink:href="pcbi.1003521.e024.jpg"></inline-graphic>
</inline-formula>
atoms. These structures are aligned and <inline-formula><inline-graphic xlink:href="pcbi.1003521.e025.jpg"></inline-graphic>
</inline-formula>
new intermediate conformers/images are generated by linearly interpolating between the end structures. The value of <italic>M</italic>
is chosen by the user and is dependent on the value of the tolerance parameter <inline-formula><inline-graphic xlink:href="pcbi.1003521.e026.jpg"></inline-graphic>
</inline-formula>
, which is the smallest energy difference between two conformers that are considered to be different.</p>
</list-item>
<list-item><p>For each image the energy is determined using the two-state potential defined in <xref ref-type="disp-formula" rid="pcbi.1003521.e021">Eq. (4)</xref>
. We identify the conformer <inline-formula><inline-graphic xlink:href="pcbi.1003521.e027.jpg"></inline-graphic>
</inline-formula>
for which energies from both the surfaces (i.e. <inline-formula><inline-graphic xlink:href="pcbi.1003521.e028.jpg"></inline-graphic>
</inline-formula>
and <inline-formula><inline-graphic xlink:href="pcbi.1003521.e029.jpg"></inline-graphic>
</inline-formula>
) are equal within the tolerance parameter <inline-formula><inline-graphic xlink:href="pcbi.1003521.e030.jpg"></inline-graphic>
</inline-formula>
.</p>
</list-item>
<list-item><p>Starting from <inline-formula><inline-graphic xlink:href="pcbi.1003521.e031.jpg"></inline-graphic>
</inline-formula>
the transition state is searched by the following iterative procedure:</p>
<list list-type="alpha-lower"><list-item><p>With appropriate choices of step-sizes <inline-formula><inline-graphic xlink:href="pcbi.1003521.e032.jpg"></inline-graphic>
</inline-formula>
and <inline-formula><inline-graphic xlink:href="pcbi.1003521.e033.jpg"></inline-graphic>
</inline-formula>
and knowledge of transition state for the present iteration <inline-formula><inline-graphic xlink:href="pcbi.1003521.e034.jpg"></inline-graphic>
</inline-formula>
, one step of steepest descent minimization is carried out on each surface using the force of the respective surface and two new sets of coordinates, <inline-formula><inline-graphic xlink:href="pcbi.1003521.e035.jpg"></inline-graphic>
</inline-formula>
and <inline-formula><inline-graphic xlink:href="pcbi.1003521.e036.jpg"></inline-graphic>
</inline-formula>
are generated, where<disp-formula id="pcbi.1003521.e037"><graphic xlink:href="pcbi.1003521.e037.jpg" position="anchor" orientation="portrait"></graphic>
<label>(5)</label>
</disp-formula>
<disp-formula id="pcbi.1003521.e038"><graphic xlink:href="pcbi.1003521.e038.jpg" position="anchor" orientation="portrait"></graphic>
<label>(6)</label>
</disp-formula>
and<disp-formula id="pcbi.1003521.e039"><graphic xlink:href="pcbi.1003521.e039.jpg" position="anchor" orientation="portrait"></graphic>
<label>(7)</label>
</disp-formula>
</p>
</list-item>
<list-item><p>A linear interpolation is performed between <inline-formula><inline-graphic xlink:href="pcbi.1003521.e040.jpg"></inline-graphic>
</inline-formula>
and <inline-formula><inline-graphic xlink:href="pcbi.1003521.e041.jpg"></inline-graphic>
</inline-formula>
to find out the conformer that resides on the cusp hypersurface. This is the new approximation for transition state i.e. <inline-formula><inline-graphic xlink:href="pcbi.1003521.e042.jpg"></inline-graphic>
</inline-formula>
</p>
</list-item>
<list-item><p>We iterate steps (3a) and (3b) until the energy difference between two transition state conformers, obtained in two successive iterations, is less than a tolerance <inline-formula><inline-graphic xlink:href="pcbi.1003521.e043.jpg"></inline-graphic>
</inline-formula>
.</p>
</list-item>
</list>
</list-item>
<list-item><p>Two separate steepest descent minimizations are performed, one on each surface, starting from the final transition state conformation <inline-formula><inline-graphic xlink:href="pcbi.1003521.e044.jpg"></inline-graphic>
</inline-formula>
and conformers separated by a user-defined RMSD are collected.</p>
</list-item>
<list-item><p>Conformers are indexed in the following sequence to construct a pathway: end structure <italic>A</italic>
, conformers collected on surface <italic>A</italic>
with increasing RMSD from the end structure <italic>A</italic>
, transition state conformer <inline-formula><inline-graphic xlink:href="pcbi.1003521.e045.jpg"></inline-graphic>
</inline-formula>
, conformers collected on surface <italic>B</italic>
with decreasing RMSD from the end structure <italic>B</italic>
, end structure <italic>B</italic>
.</p>
</list-item>
</list>
<p>Several parameters listed above need to be specified before performing the calculation. The two-state potential function is characterized by the force constants and cut-off distances of ANMs. The ANM force constant does not affect the qualitative results (or the shape of conformational change driven by the normal modes), but uniformly scales the absolute size of motions. Our choice of force constants is inconsequential, since the absolute size of the motion is adjusted by the step sizes <inline-formula><inline-graphic xlink:href="pcbi.1003521.e046.jpg"></inline-graphic>
</inline-formula>
and <inline-formula><inline-graphic xlink:href="pcbi.1003521.e047.jpg"></inline-graphic>
</inline-formula>
. The cut-off distance is usually selected in the range <inline-formula><inline-graphic xlink:href="pcbi.1003521.e048.jpg"></inline-graphic>
</inline-formula>
Å and the overall qualitative features of the pathways were found to be quite insensitive to the choice of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e049.jpg"></inline-graphic>
</inline-formula>
within this range for all the systems we studied. If desired the force constants can be estimated by fitting the crystallographic B-factors although the B-factors themselves may be biased by the crystallization conditions and crystal contacts. The energy offsets can be tuned if there are experimental information on the relative energies of the end-states. The value of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e050.jpg"></inline-graphic>
</inline-formula>
was chosen to be in the range between 10<sup>−4</sup>
and 10<sup>−5</sup>
which could be achieved by setting <italic>M</italic>
= 100 (step 1 of the algorithm). The most important parameters for an efficient implementation of the algorithm turned out to be the step-sizes involved in the transition state search on the cusp hypersurface (<italic>s<sub>A</sub>
</italic>
and <italic>s<sub>B</sub>
</italic>
in step 3a). If step-sizes are too large then the resultant movement of the transition state structure on the cusp hypersurface is large and the minimization algorithm does not work. On the other hand, if the chosen values are too small then the convergence becomes slow. For optimal values of step-sizes, short trial runs were performed for several choices, starting from large values and systematically decreasing them at each trial run until the energy of the transition state conformer decreased monotonically for the entire duration of the trial run. The starting values of <italic>s<sub>A</sub>
</italic>
and <italic>s<sub>B</sub>
</italic>
were chosen between 0.8 and 0.4 with ANM force constants set at 0.1 kcal/(mol Å<sup>2</sup>
). These values need to be adjusted if the force constants are changed by maintaining the inverse proportionality between the step-size and the force constant. Our experience shows that a few very short trial runs are sufficient for finding the optimal values of <italic>s<sub>A</sub>
</italic>
and <italic>s<sub>B</sub>
</italic>
and the overall procedure is extremely efficient. The convergence criterion <inline-formula><inline-graphic xlink:href="pcbi.1003521.e051.jpg"></inline-graphic>
</inline-formula>
was selected between 10<sup>−4</sup>
and 10<sup>−5</sup>
. The number of iterations needed for convergence ranged from 200 to 1000 for the systems studied in this paper.</p>
<p>The pathway is constituted of equally spaced structures obtained by the above mentioned algorithm between the two end structures. We have calculated several quantities to analyse the pathway and understand the conformational transition in terms of collective coordinates. The change in energy of the system along the pathway illustrates the shape of the harmonic surfaces used to describe the system. For example, if the structural change involves movements along the low frequency modes, the energy changes are smaller for a given deformation, compared to those involved in movements along high frequency (more local) movements. The cusp region along the pathway, which is easily identified as the place where the system hops from one surface to another, does not necessarily fall in the middle of the two end structures if one endpoint is more compact than the other. In order to understand the importance of the normal modes (ANM modes) of one end structure in describing the conformational transition, we calculated the cumulative correlation cosine, defined below, of few selected structures along the pathway<disp-formula id="pcbi.1003521.e052"><graphic xlink:href="pcbi.1003521.e052.jpg" position="anchor" orientation="portrait"></graphic>
<label>(8)</label>
</disp-formula>
where <inline-formula><inline-graphic xlink:href="pcbi.1003521.e053.jpg"></inline-graphic>
</inline-formula>
is the displacement of a selected conformer/image from one of the reference end structures, <inline-formula><inline-graphic xlink:href="pcbi.1003521.e054.jpg"></inline-graphic>
</inline-formula>
is the <italic>q</italic>
th normal mode of that structure and <italic>m</italic>
is the total number of modes (starting from the lowest frequency modes) used for evaluating the cumulative correlation cosine.</p>
<p>It is difficult to validate the results of our method with direct experimental observations. Many of the intermediate structures in a pathway are short-lived and may not be amenable to experimental detection. However, it is reasonable to expect that some predictions can be indirectly verified. In order to make closer connections to experiments, we have looked at the possible formation of close non-native contacts along the pathway. The hope is that some of these predictions can be tested in cross-linking experiments. We looked for pairs of residues that are far apart (>10 Å) in both the native states but come close (<5–7 Å) somewhere along the pathway. We were able to find such pairs in two of the four systems we studied.</p>
</sec>
<sec id="s3"><title>Results</title>
<sec id="s3a"><title>Adenylate Kinase</title>
<p>Adenylate kinase (AK) is an enzyme that catalyzes the transfer of a single phosphoryl group from ATP to AMP via the reversible reaction <inline-formula><inline-graphic xlink:href="pcbi.1003521.e055.jpg"></inline-graphic>
</inline-formula>
. The structure of AK consists of three domains: the AMP-binding domain (NMP), the ATP-binding domain (LID) and the CORE domain (<xref ref-type="fig" rid="pcbi-1003521-g001">Fig. (1A)</xref>
). The phosphoryl transfer reaction involves a large-scale conformational transition in AK. In the open (O) state, the NMP and LID are farther apart; and in the closed state, they are tightly packed (right and left structures in <xref ref-type="fig" rid="pcbi-1003521-g001">Fig. (1A)</xref>
). We have applied <italic>ANMPathway</italic>
on the open (O) to closed (C) transition in AK. The end-states were obtained from the crystal structures (PDB IDs: 4AKE <xref rid="pcbi.1003521-Mller1" ref-type="bibr">[62]</xref>
and 1AKE <xref rid="pcbi.1003521-Mller2" ref-type="bibr">[63]</xref>
for the O and C states, respectively). The pathway has 100 images with an RMSD of 0.1 Å between two consecutive images and the transition state corresponds to image 89 (<xref ref-type="fig" rid="pcbi-1003521-g001">Fig. (1B)</xref>
) which is situated almost at the end of the open to close transition. The transition between the functional substates of AK comprises large scale hinge-like motions of NMP and LID with respect to a rigid CORE. At the initial stage only the LID moves like a rigid body and the rest of the protein is almost unchanged. This motion corresponds to the slow rise in energy (<xref ref-type="fig" rid="pcbi-1003521-g001">Fig. (1B)</xref>
). Then the NMP starts to move and the energy rises as local structural rearrangements take place. Finally the CORE domain undergoes some changes and the transition is complete (movie S1 in Supplementary Material (SM)). The overall result is a two step transition mechanism: LID closing followed by NMP closing (or in the reverse direction: NMP opening followed by LID opening).</p>
<fig id="pcbi-1003521-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1003521.g001</object-id>
<label>Figure 1</label>
<caption><title>Conformational transition between the open and the closed states of adenylate kinase.</title>
<p><bold>A.</bold>
Structures of the open (<italic>left</italic>
, PDB ID: 4AKE) and the closed (<italic>right</italic>
, PDB ID: 1AKE) states. The LID and the NMP domains are shown in <italic>red</italic>
and <italic>orange</italic>
respectively. The CORE domain and the rest of the protein are shown in <italic>blue</italic>
. The central structure is the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e056.jpg"></inline-graphic>
</inline-formula>
trace of the transition state produced by <italic>ANMPathway</italic>
. <bold>B.</bold>
The energy of the system along the transition. Total number of images in the pathway is 100, RMSD between two consecutive images is ∼0.1 Å. The transition state corresponds to image 89. <bold>C.</bold>
Cumulative squared cosines between ANM modes and the change in structure between the initial state and a few selected conformers (images) along the transition pathway. The modes were calculated for the starting structure (open state). Here and in the counterparts generated for other test cases, the lowest frequency (slowest modes) end of the graph is enlarged in the inset. The force constants and cut-offs for both the end-states were set to 0.1 kcal/(mol Angstrom) and 15 Å in all applications in the present study, and no energy offsets were used for either of the end structures.</p>
</caption>
<graphic xlink:href="pcbi.1003521.g001"></graphic>
</fig>
<p>Because of the functional importance of domain opening/closing, it is natural to expect that the transition can be described by a small number of normal modes, as shown in <xref ref-type="fig" rid="pcbi-1003521-g001">Fig. (1C)</xref>
. For image 40, only two modes are sufficient to represent 90% of the displacement from the starting (O) structure. As we move away from the reference state, more modes are needed but the number of modes increases slowly. For attaining the other end structure (C, black curve) with a correlation cosine of more than 0.90 no more than 10–15 modes are needed, which is only ∼2% of the total number of available modes. Therefore the normal mode picture is extremely useful for studying this transition.</p>
<p>In recent years, several computational studies have examined the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e057.jpg"></inline-graphic>
</inline-formula>
allosteric transition in AK. These studies revealed a multiplicity of pathways, as well as their dependence on the initial conformers. Among them, two types of transitions appear to be consistently observed in independent studies: (i) LID closing followed by NMP closing along the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e058.jpg"></inline-graphic>
</inline-formula>
transition <xref rid="pcbi.1003521-Maragakis1" ref-type="bibr">[52]</xref>
, <xref rid="pcbi.1003521-Gur2" ref-type="bibr">[64]</xref>
–<xref rid="pcbi.1003521-Brokaw1" ref-type="bibr">[71]</xref>
and (ii) LID opening followed by NMP opening along the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e059.jpg"></inline-graphic>
</inline-formula>
transition <xref rid="pcbi.1003521-Gur2" ref-type="bibr">[64]</xref>
, <xref rid="pcbi.1003521-Whitford2" ref-type="bibr">[70]</xref>
–<xref rid="pcbi.1003521-Beckstein1" ref-type="bibr">[72]</xref>
. <xref ref-type="fig" rid="pcbi-1003521-g002">Fig. (2)</xref>
illustrates these pathways in the conformational space defined by LID-CORE and NMP-CORE angles (see <xref ref-type="fig" rid="pcbi-1003521-g002">Fig. (2)</xref>
caption for the definition of the angles), with the intermediates obtained by a recently introduced hybrid methodology, <italic>coMD</italic>
<xref rid="pcbi.1003521-Gur2" ref-type="bibr">[64]</xref>
, and by <italic>ANMPathway</italic>
, as labeled. As can be seen, mechanisms (i) and (ii) are predicted by <italic>coMD</italic>
provided that the starting points are the O and C states, respectively. The pathway (iii), on the other hand, is obtained by conducting two parallel runs, starting from both ends, and generating intermediates until the two paths merge. The initial steps conform to paths (i) and (ii) in this case. Similar transition mechanisms were recently reported by Kidera and coworkers <xref rid="pcbi.1003521-Matsunaga1" ref-type="bibr">[19]</xref>
in further support to the LID-movement-first behavior in both directions. This behavior is also in agreement with the free energy surface obtained by Woolf and coworkers <xref rid="pcbi.1003521-Beckstein1" ref-type="bibr">[72]</xref>
as a function of LID-CORE and NMP-CORE angles, although other mechanisms, such as NMP-first closure for unligated AK <xref rid="pcbi.1003521-Daily1" ref-type="bibr">[67]</xref>
or monotonous LID/NMP closing <xref rid="pcbi.1003521-Beckstein1" ref-type="bibr">[72]</xref>
have also been reported. <italic>ANMPathway</italic>
yielded a transition state very close to the closed state, i.e. by definition, most of the change in structure (89 steps out of 100) proceeded in the energy well near the open state, hence its consistency with modes accessible to end-state O, or pathway (i).</p>
<fig id="pcbi-1003521-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1003521.g002</object-id>
<label>Figure 2</label>
<caption><title>Transition pathways displayed as a function of LID-CORE and NMP-CORE angles.</title>
<p>Closed (<italic>bottom</italic>
, <italic>square</italic>
) and Open (<italic>top</italic>
, <italic>diamond</italic>
) substates of AK for the crystal structures (<italic>black</italic>
) and their MD equilibrated (<italic>magenta</italic>
) conformers. <italic>ANMPathway</italic>
and three <italic>coMD</italic>
<xref rid="pcbi.1003521-Gur2" ref-type="bibr">[64]</xref>
pathways are depicted as <italic>black</italic>
, <italic>blue</italic>
, <italic>magenta</italic>
and <italic>green</italic>
dots, respectively, with the arrows indicating the direction of reconfiguration in each case. The relative positions of the domains are defined similar to Beckstein <italic>et al.</italic>
<xref rid="pcbi.1003521-Beckstein1" ref-type="bibr">[72]</xref>
and Gur <italic>et al.</italic>
<xref rid="pcbi.1003521-Gur2" ref-type="bibr">[64]</xref>
. The NMP-CORE angle is the angle between the centers of mass of three segments I90-G100, L115-V125 and L35-A55, based on <inline-formula><inline-graphic xlink:href="pcbi.1003521.e060.jpg"></inline-graphic>
</inline-formula>
atoms and the LID-CORE angle is the angle between centers of mass of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e061.jpg"></inline-graphic>
</inline-formula>
atoms of L115-V125, I179-E185 and V125-L153.</p>
</caption>
<graphic xlink:href="pcbi.1003521.g002"></graphic>
</fig>
</sec>
<sec id="s3b"><title>Leucine Transporter (LeuT)</title>
<p>Neurotransmitter sodium symporters (NSS) are integral membrane proteins responsible for secondary transport of glycine, <italic>γ</italic>
-amino butyric acid and biogenic amines across the plasma membrane. LeuT is a bacterial orthologue of eukaryotic NSS. The protein consists of twelve transmembrane (TM) helices (TM1-12), extracellular (EC) (EL2, EL3, EL4a, El4b) and intracellular (IC) loops (IL1, IL5) and two <italic>β</italic>
-sheets. The crystal structures of LeuT in the outward-facing (OF) occluded (PDB ID: 2A65 <xref rid="pcbi.1003521-Yamashita1" ref-type="bibr">[73]</xref>
) and inward-facing (IF) open (PDB ID: 3TT3 <xref rid="pcbi.1003521-Krishnamurthy1" ref-type="bibr">[74]</xref>
) states have been resolved (left and right structures of <xref ref-type="fig" rid="pcbi-1003521-g003">Fig. (3A)</xref>
, see figure caption for color codes).</p>
<fig id="pcbi-1003521-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1003521.g003</object-id>
<label>Figure 3</label>
<caption><title>Conformational transition between the outward-facing occluded/closed (OFc) and the inward-facing open (IFo) states of leucine transporter, LeuT.</title>
<p><bold>A.</bold>
Structures of the OFc state (<italic>left</italic>
, PDB ID: 2A65) and the IFo state (<italic>right</italic>
, PDB ID: 3TT3). The scaffold domain, which does not undergo significant conformational changes, is shown in <italic>gray</italic>
and the rest of the protein is shown in <italic>blue</italic>
(EL4), <italic>red</italic>
(TM1a), <italic>orange</italic>
(TM1b), <italic>purple</italic>
(TM2 and TM7), <italic>yellow</italic>
(TM5), <italic>green</italic>
(TM6a) and <italic>lime</italic>
(TM6b). The central structure is the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e062.jpg"></inline-graphic>
</inline-formula>
trace of the transition state produced by the <italic>ANMPathway</italic>
method. <bold>B.</bold>
The energy of the system along the transition. Total number of images in the pathway is 95, RMSD between two consecutive images is ∼0.05 Å. The transition state corresponds to image 40. <bold>C.</bold>
Cumulative squared cosines between the structural change to reach a few selected conformers/images along the transition pathway and the ANM modes accessible to the starting (outward-facing occluded) structure.</p>
</caption>
<graphic xlink:href="pcbi.1003521.g003"></graphic>
</fig>
<p>We note that several residues were not resolved in these structures. We have built the ANMs based on the residues commonly resolved in the two structures and used the scaffold region (TM3, TM4, TM8 and TM9) <xref rid="pcbi.1003521-Krishnamurthy1" ref-type="bibr">[74]</xref>
for structural alignment of all conformers along the pathway.</p>
<p>The results are presented in <xref ref-type="fig" rid="pcbi-1003521-g003">Fig. (3B–C)</xref>
for the transition from the OF occluded to IF open state. The pathway is composed of 95 images with an RMSD of 0.05 Å between two consecutive images and the transition state is located at image 40 (closer to the OF occluded state). In order to analyze the transition we have looked at the rearrangements of EL4a, EL4b, TM1a, TM1b, TM2, TM5, TM6a, TM6b and TM7 which play important roles in substrate-binding and gating at the EC and IC regions <xref rid="pcbi.1003521-Krishnamurthy1" ref-type="bibr">[74]</xref>
. The motions are much more subtle compared to AK. At the initial stage there seems to be an almost rigid-body rotation in the aforementioned domains (movies S2 and S3 in SM). Then a downward movement of EL4a and EL4b closes the opening at the EC side. Subsequently, concerted motions of TM5 and TM1a occur. At the later phase of the transition, the main event is the movement of TM1a, which, along with other domains, creates an opening into the IC region. Barring the initial stage, various motions involve intra-domain movements and do not follow a strict rigid-body character. TM2 and TM7 move together for the entire duration of the conformational transition.</p>
<p>The complexity of the motions and lack of rigid-body character are reflected on the normal mode projections of various conformers along the pathway (<xref ref-type="fig" rid="pcbi-1003521-g003">Fig. (3C)</xref>
). At the initial stage when the motions are rigid-body like, a few normal modes are sufficient to describe the structural change. But as the transition progresses many more normal modes are needed to represent the displacement from the reference (OF occluded) state. For the later stage of the transition and for the end-state almost 500 modes out of possible ∼1550 modes are required to attain a cumulative correlation of 0.8. This is in sharp contrast to the transition in AK where far fewer percentages of modes were sufficient in describing the structural changes.</p>
<p>To gain more insights and validate the analysis, we compared the pathway from <italic>ANMPathway</italic>
to a 235 ns long all-atom MD trajectory using as initial structure an intermediate close to the OF occluded state <xref rid="pcbi.1003521-Cheng1" ref-type="bibr">[75]</xref>
. The details of the simulation protocols are described in the <xref ref-type="supplementary-material" rid="pcbi.1003521.s005">text S1</xref>
of the Supplemental Material (SM). The time evolution of the structure during this transition was probed by monitoring a relevant order parameter, namely the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e063.jpg"></inline-graphic>
</inline-formula>
distance between the binding-site residues N21 and S256 shown in <xref ref-type="fig" rid="pcbi-1003521-g004">Fig. (4A)</xref>
. Notably, a spontaneous transition to IF open state was observed in this conventional (unbiased) full atomistic MD simulation. The MD trajectory thus provides an important data-set for benchmarking the <italic>ANMPathway</italic>
method. <xref ref-type="fig" rid="pcbi-1003521-g004">Fig. (4B)</xref>
compares the projection from the <italic>ANMPathway</italic>
(white curve) and the MD trajectory on the space of two order parameters, the N21-S256 <inline-formula><inline-graphic xlink:href="pcbi.1003521.e064.jpg"></inline-graphic>
</inline-formula>
distance and the RMSD from the OF occluded state. The MD data are represented as a crude free energy calculated by taking the logarithm of the two-dimensional histogram of the above mentioned two order parameters shown in <xref ref-type="fig" rid="pcbi-1003521-g004">Fig. (4B)</xref>
. The pathway predicted by <italic>ANMPathway</italic>
is constructed on a smooth potential energy function and has no thermal fluctuations. As such, it is representative of an average pathway of the real system and it should go through the low energy regions of the free energy landscape obtained from the MD simulation of LeuT embedded in fully solvated membrane lipids at finite temperature. This is exactly the behavior we observe in <xref ref-type="fig" rid="pcbi-1003521-g004">Fig. (4B)</xref>
for the most parts where all-atom MD data are available. This agreement is satisfactory, given the minimal computational cost required by <italic>ANMPathway</italic>
compared to that (several orders of magnitude larger) required for the full scale all-atom MD simulation.</p>
<fig id="pcbi-1003521-g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1003521.g004</object-id>
<label>Figure 4</label>
<caption><title>Comparison with all-atom simulation results.</title>
<p><bold>A.</bold>
Time trace of the distance between <inline-formula><inline-graphic xlink:href="pcbi.1003521.e065.jpg"></inline-graphic>
</inline-formula>
atoms of residues 21 and 256 from a 235 ns long conventional MD simulation of the fully solvated system. The simulation was started from a structure obtained from a targeted MD simulation originated from the OFc state. The system undergoes a spontaneous transition to the IFo state. <bold>B.</bold>
Comparison of the <italic>ANMPathway</italic>
method and all-atom MD in the space of two order parameters. The all-atom MD results are shown as a pseudo free energy landscape <inline-formula><inline-graphic xlink:href="pcbi.1003521.e066.jpg"></inline-graphic>
</inline-formula>
, where <italic>P</italic>
is the 2D distribution. The color-scale goes from <italic>blue</italic>
(low energy) to <italic>red</italic>
(high energy). The pathway predicted by <italic>ANMPathway</italic>
(<italic>white</italic>
line) goes mostly through the low energy regions of the free energy landscape.</p>
</caption>
<graphic xlink:href="pcbi.1003521.g004"></graphic>
</fig>
<p>There is another crystal structure of LeuT which models the outward-facing open state (PDB ID: 3TT1 <xref rid="pcbi.1003521-Krishnamurthy1" ref-type="bibr">[74]</xref>
). The sequence of functional states in the reaction cycle is: OF open (PDB ID: 3TT1)<inline-formula><inline-graphic xlink:href="pcbi.1003521.e067.jpg"></inline-graphic>
</inline-formula>
OF occluded (PDB ID: 2A65)<inline-formula><inline-graphic xlink:href="pcbi.1003521.e068.jpg"></inline-graphic>
</inline-formula>
IF open (PDB ID: 3TT3). There are important differences in the helical orientations of several TM helices between the OF open and OF occluded structures (<xref ref-type="fig" rid="pcbi-1003521-g005">Fig. (5A)</xref>
) even though their overall architectures are quite similar. It is natural to ask whether the <italic>ANMPathway</italic>
could predict the existence of OF closed state along a transition pathway calculated between OF open and IF open states. We indeed found a conformer which is very close to the OF closed state (RMSD from 2A65: ∼1.0 Å) along the pathway between the crystal structures of OF open and IF open states. The detection of the occluded intermediate is studied by monitoring order parameters that describe the instantaneous conformations of the TM helices responsible for gating and binding of ions. The order parameters are the center of mass (COM) distances between the pairs of helices TM1a-TM10, TM6b-TM10 and TM1a-TM6b. The results are shown in <xref ref-type="fig" rid="pcbi-1003521-g005">Fig. (5B)</xref>
. In all three distance profiles the intermediate is detected as indicated by the yellow arrows. Given the simplicity of the potential energy function, it is quite remarkable that the method is capable of detecting a functionally relevant and experimentally observed intermediate state. This fact highlights the usefulness of the method as well as significance of global modes of motion in facilitating the conformational transition of transporters.</p>
<fig id="pcbi-1003521-g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1003521.g005</object-id>
<label>Figure 5</label>
<caption><title>Occurrence of outward-facing occluded/closed (OFc) state along the computed transition pathway between outward-facing open (OFo) and inward-facing open (IFo) states of LeuT.</title>
<p><bold>A.</bold>
Relative orientations of TM1, TM6 and TM10 of LeuT in the OFo (<italic>orange</italic>
), OFc (<italic>yellow</italic>
) and IFo (<italic>cyan</italic>
) crystal structures. The OFo crystal structure is shown in transparent cartoon. Inward tilting of the TM1b and TM6a segments contributes to the closure of the extracellular vestibule, indicating the decreased distances of TM1b-TM10 and TM6a-TM10. Outward tilting of TM1a dominates the opening of the intracellular vestibule, resulting in the increased distance between TM1a andTM6b. <bold>B.</bold>
Variation of the center of mass (COM) distances of TM1a-TM10 (<italic>red</italic>
+), TM6b-TM10 (<italic>blue</italic>
+), and TM1a-TM6b (<italic>gray</italic>
+) as the LeuT undergoes transition from the OFo to IFo states. For COM distance calculations, TM1a (R11 to A22), TM1b (L25 to A35), TM6a (G242 to L255), TM6b (F259 to Y268), and TM10 (K398 to V412) <inline-formula><inline-graphic xlink:href="pcbi.1003521.e069.jpg"></inline-graphic>
</inline-formula>
atoms are taken from <italic>ANMPathway</italic>
calculations of LeuT from IFo to OFo. <italic>Yellow arrows</italic>
point to the values found in the OFc crystal structure. Clearly, occluded intermediates are identified by <italic>ANMPathway</italic>
calculations (highlighted by the horizontal <italic>blue arrow</italic>
).</p>
</caption>
<graphic xlink:href="pcbi.1003521.g005"></graphic>
</fig>
</sec>
<sec id="s3c"><title>Glutamate Transporter (Glt<sub>Ph</sub>
)</title>
<p>Excitatory amino acid transporters (EAATS) constitute a class of integral membrane proteins that are responsible for secondary active transport of amino acids like glutamate and aspartate across the plasma membrane. The aspartate transporter <inline-formula><inline-graphic xlink:href="pcbi.1003521.e070.jpg"></inline-graphic>
</inline-formula>
, an archaeal orthologue of eukaryotic EAAT, is broadly used as a structural prototype, being functionally resolved in multiple states. The protein is a homotrimer. Each protomer consists of eight TM helices (TM1-8) and two helix-turn-helix motifs (HP1 and HP2) at the substrate-binding core <xref rid="pcbi.1003521-Yernool1" ref-type="bibr">[76]</xref>
. According to the alternating access mechanism that enables the transport of substrate, the trimer alternates between OF and IF states and <italic>vice versa</italic>
, via structural changes in all three protomers, to expose the substrate-binding site to the EC and IC regions, respectively. Crystal structures of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e071.jpg"></inline-graphic>
</inline-formula>
in OF (PDB ID: 1XFH <xref rid="pcbi.1003521-Yernool1" ref-type="bibr">[76]</xref>
), IF (PDB ID: 3KBC <xref rid="pcbi.1003521-Reyes1" ref-type="bibr">[77]</xref>
) states, as well a mixed intermediate state (iOF) with two protomers in IF conformation and one in an intermediate between OF and IF conformations (PDB ID: 3V8G <xref rid="pcbi.1003521-Verdon1" ref-type="bibr">[78]</xref>
) have been determined (<xref ref-type="fig" rid="pcbi-1003521-g006">Fig. (6A)</xref>
). The iOF state has been suggested to be relevant to uncoupled anion permeation during the transport process <xref rid="pcbi.1003521-Verdon1" ref-type="bibr">[78]</xref>
. This asymmetric structure closely approximates the intermediate predicted earlier by a combined experimental and computational study <xref rid="pcbi.1003521-Jiang1" ref-type="bibr">[79]</xref>
.</p>
<fig id="pcbi-1003521-g006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1003521.g006</object-id>
<label>Figure 6</label>
<caption><title>Conformational transition between states with two protomers facing inward and all protomers facing inward of the glutamate transporter, Glt<sub>Ph</sub>
.</title>
<p><bold>A.</bold>
<italic>Left</italic>
: Structure of the intermediate state (iOF) state where two protomers are in IF and the third (<italic>green</italic>
) protomer is in the OF conformations (PDB ID: 3V8G); <italic>right</italic>
: IF state where all protomers are in IF conformations (PDB ID: 3KBC). The central structure is the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e072.jpg"></inline-graphic>
</inline-formula>
trace of the transition state produced by <italic>ANMPathway</italic>
. <bold>B.</bold>
The energy of the system along the transition. Total number of images in the pathway is 79, RMSD between two consecutive images is ∼0.1 Å. The transition state corresponds to image 47. <bold>C.</bold>
Cumulative squared cosines between the structural change to reach a few selected conformers/images along the transition pathway and the ANM modes accessible to the starting (iOF) structure.</p>
</caption>
<graphic xlink:href="pcbi.1003521.g006"></graphic>
</fig>
<p>We examined the transition between the iOF and IF states, using <italic>ANMPathway</italic>
. The crystal structures reveal that protomers can be divided into two domains: a trimerization domain (TM1, TM2, TM4 and TM5) that closely maintains its internal conformation during the transition, and a transport domain (TM3, TM6, HP1, TM7, HP2 and TM8) that practically undergoes a downward rigid-body movement (perpendicular to the membrane) relative to the trimerization domain. Our pathway is made of 79 images with an RMSD of 0.1 Å between two consecutive images and the transition state is located at image 47 (<xref ref-type="fig" rid="pcbi-1003521-g006">Fig. (6B)</xref>
). In accord with experimental data, the trimerization domain does not exhibit significant change in its internal structure during the entire transition while the transport domain of the reconfiguring protomer (colored green in <xref ref-type="fig" rid="pcbi-1003521-g006">Fig. (6A)</xref>
) first moves like a rigid-body and then undergoes intra-domain rearrangements. The short segment of (the broken helix) TM8 moves first, then the lower longer segment and at the end the entire helix moves. Similarly, the upward part of TM7 moves initially followed by the rest of the helix. The movements of HP1 and HP2 are particularly important since they are involved in substrate gating <xref rid="pcbi.1003521-Reyes1" ref-type="bibr">[77]</xref>
. Along the inward transition of the protomer, HP2 moves first followed by HP1 toward the end, consistent with the higher (initial) mobility of HP2 observed in all-atom MD simulations <xref rid="pcbi.1003521-Zomot1" ref-type="bibr">[80]</xref>
, <xref rid="pcbi.1003521-Shrivastava1" ref-type="bibr">[81]</xref>
. See the movies S4 and S5 in SM. The long flexible loop in the extracellular part also undergoes significant movements during the transition. A simpler view of these local rearrangements can be obtained if one constructs two blocks as suggested by Reyes <italic>et al.</italic>
<xref rid="pcbi.1003521-Reyes1" ref-type="bibr">[77]</xref>
and predicted by ANM analysis to constitute two distinctive substructures subject to anticorrelated motions <xref rid="pcbi.1003521-Lezon1" ref-type="bibr">[82]</xref>
. Block 1 is composed of HP1 and the lower part of TM7; and block 2 consists of HP2 and the upper part of TM8 up to the point where the helix is broken. Block 2 moves first followed by block 1 (movie S6 in SM).</p>
<p><xref ref-type="fig" rid="pcbi-1003521-g006">Fig. (6C)</xref>
displays the number of normal modes needed for describing the transition to several images/conformers starting from the iOF state. At the initial stage when motions are more rigid-body like, very few modes are sufficient to describe the structural changes. However, similar to LeuT, this number increases as more localized events that do not conform to <italic>en bloc</italic>
movements of low frequency modes become important. These involve flexible regions within the transport domain. However, the total number of modes for a reasonable description of the overall structural change remains significantly smaller than those available (e.g. <100 modes accomplish a cumulative squared cosine of 0.8 with the reference state), supporting the utility of low frequency modes for efficiently mapping the transition pathway.</p>
<p>The complex structural changes result in the formation of eight non-native contacts along the transition. The formation of non-native contact is defined as two residues that are more than 10 Å apart in the end-states but come closer by less than 7 Å during the transition. The non-native contact forming pairs are Val58-Ala358, Val62-Ala353, Leu152-Leu347, Gln220-Met385, Val274-Ala391, Thr275-Gly357, Gly280-Pro356 and Val355-Ile361. All the pairs belong to the protomer that is undergoing the transition to the IF conformation (chain C, colored green in <xref ref-type="fig" rid="pcbi-1003521-g006">Fig. (6A)</xref>
). It is worth noting that none of these pairs corresponds to the cysteine cross-link present in the crystal structure of the IF state <xref rid="pcbi.1003521-Reyes1" ref-type="bibr">[77]</xref>
. <xref ref-type="fig" rid="pcbi-1003521-g007">Fig. (7)</xref>
shows the distance profiles for three of these pairs along the transition. Most of these contacts involve pairs where one residue is in the trimerization domain and the other in the transport domain (movies S7 to S9 in the SM). These observations provide a route to test the predictions of the <italic>ANMPathway</italic>
method against experiments.</p>
<fig id="pcbi-1003521-g007" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1003521.g007</object-id>
<label>Figure 7</label>
<caption><title>Occurrence of transient non-native contacts along the conformational transition of Glt<sub>Ph</sub>
between iOF and IF states.</title>
<p><bold>A.</bold>
Distances of three contact-forming pairs plotted against the index of the conformers along the pathway. All contacts are formed between residues belonging to the same protomer. Two residues are said to form a non-native contact if they are more than 10 Å apart in the reference states but somewhere along the transition, the pairwise distance between them becomes less than 7 Å. <bold>B.</bold>
Conformer at the point of closest contact between Val62 and Ala353 (lowest point of the <italic>red</italic>
curve in <bold>A</bold>
, corresponding to the image 46). The contact forming residues, Val62 and Ala353 are shown as <italic>blue</italic>
and <italic>red</italic>
spheres respectively. The color code is: trimerization domain, <italic>cyan</italic>
; TM3, <italic>deep sky blue</italic>
; TM6, <italic>lime</italic>
; HP1, <italic>yellow</italic>
; TM7, <italic>orange</italic>
; HP2, <italic>orange</italic>
; TM8, <italic>magenta</italic>
and the rest, <italic>gray</italic>
.</p>
</caption>
<graphic xlink:href="pcbi.1003521.g007"></graphic>
</fig>
<p>We have compared the <italic>ANMPathway</italic>
results with the data collected from all-atom simulations. It is very challenging to simulate a spontaneous transition by straightforward conventional MD simulations. In order to perform a qualitative comparison, we have adopted the following protocol. First, a targeted MD (TMD) pathway is generated between the end-states (iOF and IF) with targeting forces acting on the backbone atoms only (see SM <xref ref-type="supplementary-material" rid="pcbi.1003521.s005">text S1</xref>
for details). Then, we launched a series of conventional MD runs from various intermediates visited during the TMD. These runs might be expected to follow the local free energy gradient and proceed along the result from <italic>ANMPathway</italic>
provided that the latter offers a reasonable approximation to the actual transition pathway. In order to understand the transition in terms of a simple order parameters we have used, as order parameter, the <italic>z</italic>
component of the distance vector between the COM of the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e073.jpg"></inline-graphic>
</inline-formula>
atoms of the transport domain of a conformer and that of the crystal structure of the iOF state. <xref ref-type="fig" rid="pcbi-1003521-g008">Fig. (8)</xref>
shows the projection of the predicted pathway and snapshots from unbiased MD runs initiated from various intermediate structures (shown in different colors) on the space spanned by the above-mentioned order parameter and by the RMSD from the end structure (iOF state). Except for one of the trajectories (shown in red points, near the starting point), the MD runs yielded snapshots in accord with the transition pathway predicted by <italic>ANMPathway</italic>
. The iOF state of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e074.jpg"></inline-graphic>
</inline-formula>
represents an intermediate between the OF and IF states, it is conceivable that the MD runs starting from this intermediate (shown in red) tend to go back to the more stable OF state, instead of drifting toward the IF state. This is primarily due to the proximity of the initial structure to the deep free energy basin of the more stable (OF) end-state. Overall, these data validate the ability of <italic>ANMPathway</italic>
to provide a meaningful description of the structural changes involved in the global transition of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e075.jpg"></inline-graphic>
</inline-formula>
protomers as they reconfigure from OF to IF states.</p>
<fig id="pcbi-1003521-g008" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1003521.g008</object-id>
<label>Figure 8</label>
<caption><title>Comparison of <italic>ANMPathway</italic>
with all-atom simulation results for the Glt<sub>Ph</sub>
transition.</title>
<p>The <italic>black</italic>
line is the projection of the pathway on a space spanned by two order parameters, RMSD from the endpoint (iOF state) and the <italic>z</italic>
component of the displacement of the center of mass (COM) of the transport domain (based on <inline-formula><inline-graphic xlink:href="pcbi.1003521.e076.jpg"></inline-graphic>
</inline-formula>
atoms) with respect to the initial crystal structure. Negative values along the ordinate point to the IC region. The dots of different colors are the projection of conformers sampled in unbiased all-atom MD simulation runs initiated from various points along a TMD path between the two end-states. For the middle section the data clouds from the unbiased MD simulations cluster around the pathway predicted by <italic>ANMPathway</italic>
.</p>
</caption>
<graphic xlink:href="pcbi.1003521.g008"></graphic>
</fig>
</sec>
<sec id="s3d"><title>Sarcoplasmic Reticulum Ca<sup>2+</sup>
-ATPase (SERCA)</title>
<p>Calcium transporting pump of sarco/endoplasmic reticulum (SERCA) is an integral membrane protein that pumps <inline-formula><inline-graphic xlink:href="pcbi.1003521.e077.jpg"></inline-graphic>
</inline-formula>
ions from calcium-poor cytoplasm of the muscle cell to the calcium-rich lumen of the sarcoplasmic reticulum at the expense of ATP hydrolysis. This process gets rid of the excess <inline-formula><inline-graphic xlink:href="pcbi.1003521.e078.jpg"></inline-graphic>
</inline-formula>
ions in the cytoplasm caused by their release from the lumen during muscle contraction and reverts the muscle to relaxed state. The protein is composed of a single polypeptide chain of 994 amino acids that form three cytoplasmic domains (nucleotide-binding domain N, phosphorylation domain P, actuator domain A) and ten TM helices (M1–10) (<xref ref-type="fig" rid="pcbi-1003521-g009">Fig. (9A)</xref>
). Extensive structural studies have revealed atomic models of various functionally relevant states in the pumping cycle <xref rid="pcbi.1003521-Toyoshima1" ref-type="bibr">[83]</xref>
. We have used the <italic>ANMPathway</italic>
method to explore the transition between the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e079.jpg"></inline-graphic>
</inline-formula>
(PDB ID: 1SU4 <xref rid="pcbi.1003521-Toyoshima2" ref-type="bibr">[84]</xref>
) and E1.ATP (PDB ID: 1T5S <xref rid="pcbi.1003521-Sorensen1" ref-type="bibr">[85]</xref>
) states. In the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e080.jpg"></inline-graphic>
</inline-formula>
state the calcium ions can dissociate from the transmembrane binding sites but, in the E1P state (i.e. phosphorylated <inline-formula><inline-graphic xlink:href="pcbi.1003521.e081.jpg"></inline-graphic>
</inline-formula>
state) they are occluded and can not go back to the cytoplasmic side. The architecture of the E1.ATP state is almost identical to that of E1P or the transition state analog E1∼P.ADP state (<inline-formula><inline-graphic xlink:href="pcbi.1003521.e082.jpg"></inline-graphic>
</inline-formula>
RMSD<0.5 Å). Therefore, at the level of a <inline-formula><inline-graphic xlink:href="pcbi.1003521.e083.jpg"></inline-graphic>
</inline-formula>
CG model, a transition pathway between the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e084.jpg"></inline-graphic>
</inline-formula>
and E1.ATP states can provide important insights into the large scale motions responsible for occlusion of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e085.jpg"></inline-graphic>
</inline-formula>
ions in the transmembrane binding sites.</p>
<fig id="pcbi-1003521-g009" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1003521.g009</object-id>
<label>Figure 9</label>
<caption><title>Conformational transition between E1.2CA<sup>2+</sup>
and E1.ATP states of the Sarcoplasmic Reticulum Ca<sup>2+</sup>
-ATPase (SERCA).</title>
<p><bold>A.</bold>
Structures of the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e086.jpg"></inline-graphic>
</inline-formula>
(<italic>left</italic>
, PDB ID: 1SU4) and the E1.ATP (∼E1P) (<italic>right</italic>
, PDB ID: 1T5S) states. The P, A and N domains are shown in <italic>blue</italic>
, <italic>red</italic>
and <italic>green</italic>
respectively. Rest of the protein including the transmembrane domain is shown in <italic>gray</italic>
. The central structure is the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e087.jpg"></inline-graphic>
</inline-formula>
trace of the transition state produced by <italic>ANMPathway</italic>
. <bold>B.</bold>
The energy of the system along the transition. Total number of images in the pathway is 106, RMSD between two consecutive images is ∼0.2 Å. The transition state corresponds to image 74. <bold>C.</bold>
Cumulative squared cosines between the structural change to reach a few selected conformers/images along the transition pathway and the ANM modes accessible to the starting (<inline-formula><inline-graphic xlink:href="pcbi.1003521.e088.jpg"></inline-graphic>
</inline-formula>
state) structure.</p>
</caption>
<graphic xlink:href="pcbi.1003521.g009"></graphic>
</fig>
<p>The pathway consists of 106 images/conformers with an RMSD of 0.2 Å between two consecutive images and the transition state is located at image 74 which is closer to the E1.ATP state (<xref ref-type="fig" rid="pcbi-1003521-g009">Fig. (9B)</xref>
). Careful examination of the pathway obtained by <italic>ANMPathway</italic>
reveals that at the initial stage only the N domain moves while the rest of the protein remains fixed. This is reflected in the initial slow rise in the energy as shown in <xref ref-type="fig" rid="pcbi-1003521-g009">Fig. (9B)</xref>
. As the N domain comes closer to the P and A domains, first the P domain undergoes intra-domain changes and then the A domain rotates. This rotation causes upward motions of M1 and M2 and in this process M1 helix adopts a bent conformation. The energy rises quickly as the cytoplasmic domains close and local structural rearrangements take place. The bending of the M1 helix is responsible for shielding the TM calcium-binding sites from the cytoplasmic side and helps to form the occluded state <xref rid="pcbi.1003521-Toyoshima1" ref-type="bibr">[83]</xref>
, <xref rid="pcbi.1003521-Toyoshima3" ref-type="bibr">[86]</xref>
(movie S10 in SM). These observations can be regarded as a crude qualitative picture of the conformational couplings among various parts of SERCA that give rise to the occluded state. The initial motion of the N domain which costs little energy suggests that the protein can exist in alternative conformations where the N domain may be much closer to other cytoplasmic domains compared to the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e089.jpg"></inline-graphic>
</inline-formula>
crystal structure. This observation is validated by other simulation and experimental studies <xref rid="pcbi.1003521-EspinozaFonseca1" ref-type="bibr">[87]</xref>
, <xref rid="pcbi.1003521-Winters1" ref-type="bibr">[88]</xref>
. <xref ref-type="fig" rid="pcbi-1003521-g009">Fig. (9C)</xref>
shows that only a handful of modes are sufficient for describing the transition. This is due to the initial rigid-body motion of the N domain, in good agreement with the slow modes predicted by the ANM based on the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e090.jpg"></inline-graphic>
</inline-formula>
structure. The entire pathway can be projected onto ∼50 modes to take into account more than 80% of the displacement. The normal mode picture is therefore very useful for exploring the transition of this system.</p>
<p>The distance profiles of contact-forming residues along the transition are shown in <xref ref-type="fig" rid="pcbi-1003521-g010">Fig. (10)</xref>
. The formation of non-native contact is defined as two residues that are more than 10 Å apart in the reference states but come close by less 5 Å during the transition. All contact-forming pairs involve one residue in the M1 helix. The other residue is either on the M4 helix or on the cytoplasmic loop region (see movies S11 to S13 in the SM). These contacts form as the result of bending of the M1 helix which is thought to be responsible for the formation of the occluded state. Therefore the non-native contacts can be probed experimentally to validate the predictions of method as well as to establish structural changes that have functional consequence.</p>
<fig id="pcbi-1003521-g010" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1003521.g010</object-id>
<label>Figure 10</label>
<caption><title>Occurrence of transient non-native contacts along the conformational transition of SERCA.</title>
<p><bold>A.</bold>
Distances between contact-forming residue pairs plotted against the index of the conformers along the pathway. Two residues are said to form non-native contact if they are more than 10 Å apart in the reference states but somewhere along the transition, the pairwise distance between them becomes less than 5 Å. <bold>B.</bold>
Intermediate conformer at the point of closest contact between Leu49 (cyan sphere) and Lys246 (<italic>purple</italic>
sphere) (minimum of the red curve in <bold>A</bold>
, corresponding to the image index 79). The color code is same as in <xref ref-type="fig" rid="pcbi-1003521-g009">Fig. (9)</xref>
.</p>
</caption>
<graphic xlink:href="pcbi.1003521.g010"></graphic>
</fig>
</sec>
</sec>
<sec id="s4"><title>Discussion</title>
<p>We presented a new computational method, <italic>ANMPathway</italic>
, for constructing the most probable (lowest energy) transition pathway between two stable endpoints of a conformational transition. Conceptually, <italic>ANMPathway</italic>
represents a direct application of the string method to a two-state CG system approximated by ANM energy surfaces. For this reason, the method is simple and efficient; it can produce a completely optimized pathway for a 1000 residue protein in about one hour on a single CPU of a standard desktop computer. We have implemented the method as an open access web server with user-friendly features at <ext-link ext-link-type="uri" xlink:href="http://anmpathway.lcrc.anl.gov/anmpathway.cgi">http://anmpathway.lcrc.anl.gov/anmpathway.cgi</ext-link>
. Although there exist other servers for exploring the transition paths between pairs of known endpoints <xref rid="pcbi.1003521-Yang1" ref-type="bibr">[53]</xref>
, <xref rid="pcbi.1003521-Zheng1" ref-type="bibr">[55]</xref>
, <xref rid="pcbi.1003521-Franklin1" ref-type="bibr">[61]</xref>
, <xref rid="pcbi.1003521-Tekpinar1" ref-type="bibr">[89]</xref>
–<xref rid="pcbi.1003521-Sfriso1" ref-type="bibr">[91]</xref>
, <italic>ANMPathway</italic>
is perhaps one of the simplest approaches that efficiently provides a unique solution for the most probable pathway, with a minimal number of parameters and no biased simulations using the ANMs for the two endpoints. The resulting pathway can be interpreted as the minimum energy pathway on the two-state elastic network surface with a cusp hypersurface. The presence of cusp hypersurface does not seem to influence the qualitative nature of the pathway as evident from the comparison of our results on adenylate kinase with those from other methods. Franklin <italic>et al.</italic>
<xref rid="pcbi.1003521-Franklin1" ref-type="bibr">[61]</xref>
in their MinActionPath method effectively used a similar two-state potential with cusp hypersurface and their results for the AK system are in good agreement with ours. We note that the presence of cusp hypersurface will have noticeable effects on the quantitative details of the pathway, especially near the transition state region. The energetic cost of breaking non-bonded contacts increases as one moves away from the native state which is reflected in the rapid change in energy near the transition state. A smooth surface, like the one adopted in other similar studies <xref rid="pcbi.1003521-Maragakis1" ref-type="bibr">[52]</xref>
, <xref rid="pcbi.1003521-Zheng1" ref-type="bibr">[55]</xref>
, <xref rid="pcbi.1003521-Tekpinar1" ref-type="bibr">[89]</xref>
, will have a better representation of the underlying physics compared to our present two-state potential with a cusp hypersurface. This is the sacrifice we have to make in order to exploit the simplicity and algorithmic efficiency imparted by the two-state potential. However, the influence of this drawback on the pathways produced by <italic>ANMPathway</italic>
was verified to be minimal and localized in the transition state region. This is clearly demonstrated by the refinement of the pathways by systematic smoothing of the two-state potential as described below. The method of Zheng <italic>et al.</italic>
<xref rid="pcbi.1003521-Zheng1" ref-type="bibr">[55]</xref>
had been implemented as a web server called AD-ENM which allowed us to directly compare our results with that of a similar method based on a more realistic energy surface. The pathways produced by the AD-ENM server had different numbers of images than the corresponding pathways generated by <italic>ANMPathway</italic>
. For the purpose of comparison, we projected the pathways using the same sets of order parameters already used to analyze <italic>ANMPathway</italic>
results in the previous section. For the AK case, it seems that AD-ENM produced a qualitatively different pathway than <italic>ANMPathway</italic>
(SM <xref ref-type="supplementary-material" rid="pcbi.1003521.s001">Fig. S1</xref>
). The AD-ENM pathway is closer to the scenario (iii) in <xref ref-type="fig" rid="pcbi-1003521-g002">Fig. (2)</xref>
, whereas <italic>ANMPathway</italic>
result resembles scenario (i) in the same figure. Both methods produced qualitatively similar pathways for LeuT and <inline-formula><inline-graphic xlink:href="pcbi.1003521.e091.jpg"></inline-graphic>
</inline-formula>
as illustrated by the projections presented in SM <xref ref-type="supplementary-material" rid="pcbi.1003521.s002">Figs. S2</xref>
and <xref ref-type="supplementary-material" rid="pcbi.1003521.s003">S3</xref>
respectively. The AD-ENM server makes small modifications to the end-states which can explain the discrepancy found near one of the endpoints in SM <xref ref-type="supplementary-material" rid="pcbi.1003521.s002">Fig. S2</xref>
. The <italic>ANMPathway</italic>
method is also closely related to the PNM method of Maragakis and Karplus <xref rid="pcbi.1003521-Maragakis1" ref-type="bibr">[52]</xref>
; our two-state potential is a limiting form of their smooth two-state ENM, when the mixing parameter goes to zero. We use a novel and efficient algorithm for searching for the minimum energy structure on the cusp hypersurface, which is equivalent to the saddle point search in the PNM method. If a sufficiently small value of mixing parameter is used, the PNM result will be very similar to ours and this is indeed observed for the AK system.</p>
<p>The effect of the cusp hypersurface on the potential energy surface was further examined by refining the pathway on a smoothed surface based on the following potential function <xref rid="pcbi.1003521-Best1" ref-type="bibr">[54]</xref>
, <xref rid="pcbi.1003521-Yang2" ref-type="bibr">[57]</xref>
,<disp-formula id="pcbi.1003521.e092"><graphic xlink:href="pcbi.1003521.e092.jpg" position="anchor" orientation="portrait"></graphic>
<label>(9)</label>
</disp-formula>
where <inline-formula><inline-graphic xlink:href="pcbi.1003521.e093.jpg"></inline-graphic>
</inline-formula>
and <inline-formula><inline-graphic xlink:href="pcbi.1003521.e094.jpg"></inline-graphic>
</inline-formula>
are ANM potentials defined in <xref ref-type="disp-formula" rid="pcbi.1003521.e010">Eq. (1)</xref>
. The parameter <inline-formula><inline-graphic xlink:href="pcbi.1003521.e095.jpg"></inline-graphic>
</inline-formula>
determines the amount of mixing and the height of the barrier. In the limit <inline-formula><inline-graphic xlink:href="pcbi.1003521.e096.jpg"></inline-graphic>
</inline-formula>
, the potential defined above becomes the two-state potential with a cusp hypersurface defined previously in <xref ref-type="disp-formula" rid="pcbi.1003521.e021">Eq. (4)</xref>
i. e. <inline-formula><inline-graphic xlink:href="pcbi.1003521.e097.jpg"></inline-graphic>
</inline-formula>
. Therefore, if we refine the pathway after smoothing the potential using <xref ref-type="disp-formula" rid="pcbi.1003521.e092">Eq. (9)</xref>
with a high value of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e098.jpg"></inline-graphic>
</inline-formula>
, the resulting pathway should not be very different. We tested this hypothesis by performing zero temperature string method <xref rid="pcbi.1003521-E1" ref-type="bibr">[92]</xref>
calculation on the smoothed energy surface starting from the final path produced by the <italic>ANMPathway</italic>
method. The results for AK are shown in <xref ref-type="fig" rid="pcbi-1003521-g011">Fig. (11)</xref>
. It is noteworthy that, for <inline-formula><inline-graphic xlink:href="pcbi.1003521.e099.jpg"></inline-graphic>
</inline-formula>
, the RMSDs between corresponding images along the two pathways are extremely small (compared to the the resolution of the structures) and the position of the transition state is almost unchanged (the projections of the pathways on the space spanned by the two order parameters described in <xref ref-type="fig" rid="pcbi-1003521-g002">Fig. (2)</xref>
are shown in the SM <xref ref-type="supplementary-material" rid="pcbi.1003521.s004">Fig. S4</xref>
). As the value of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e100.jpg"></inline-graphic>
</inline-formula>
is decreased the RMSDs increase and the position of the transition state moves slowly to the right. These observations demonstrate that the algorithm used in <italic>ANMPathway</italic>
has the ability to find a proper transition state on the two-state surface with the cusp hypersurface and also illustrate that presence of cusp hypersurface does not have a significant impact on the over all features of the two-state potential.</p>
<fig id="pcbi-1003521-g011" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1003521.g011</object-id>
<label>Figure 11</label>
<caption><title>Effect of smoothing of the two-state potential on the adenylate kinase transition pathway.</title>
<p><bold>A.</bold>
Energy of the system along the pathway for several pathways obtained by refining the original pathway by zero temperature string method calculation on the smoothed potential defined in <xref ref-type="disp-formula" rid="pcbi.1003521.e092">Eq. (9)</xref>
. <bold>B.</bold>
RMSD of the refined pathways from the pathway produced by <italic>ANMPathway</italic>
.</p>
</caption>
<graphic xlink:href="pcbi.1003521.g011"></graphic>
</fig>
<p>A desirable feature of any pathway method is the ability to detect functionally relevant intermediate structures. The two-state potential energy function used in the <italic>ANMPathway</italic>
may be too simple for such purpose. If an intermediate has a topology which is very different from both the end-states then it is unlikely to be detected by this method. However, it is reasonable to expect that an intermediate that maintains the overall fold and much of the secondary structure shared by the two endpoints would show up along the computed transition pathway. This was verified in the transition pathway of LeuT from OF open to IF open state, via the intermediate OF occluded state (see <xref ref-type="fig" rid="pcbi-1003521-g005">Fig. (5)</xref>
) consistent with experimental data. We note that Adelman <italic>et al.</italic>
<xref rid="pcbi.1003521-Adelman1" ref-type="bibr">[93]</xref>
also observed the passage over an occluded intermediate in their weighted-ensemble simulations of the transition of Mhp1 (another transporter that shares LeuT fold) between its OF and IF states. In their case, a significantly more detailed two-state G<inline-formula><inline-graphic xlink:href="pcbi.1003521.e101.jpg"></inline-graphic>
</inline-formula>
potential (with several parameters) is used for each endpoint, together with a computationally expensive Monte Carlo sampling algorithm (compared to <italic>ANMPathway</italic>
). The present approach offers the multiple advantages of being simpler and more efficient, while detecting at the same level of resolution an experimentally validated intermediate.</p>
<p>The effective two-state CG potential used by <italic>ANMPathway</italic>
is meant to approximate the true PMF of the system with respect to the <inline-formula><inline-graphic xlink:href="pcbi.1003521.e102.jpg"></inline-graphic>
</inline-formula>
coordinates chosen as collective variables <xref rid="pcbi.1003521-Maragliano1" ref-type="bibr">[14]</xref>
. The basic assumption of such structure-based CG models is that the information about the native contact topology is sufficient for describing large scale transitions. Yet, the effective CG potential can produce pathways with various degrees of complexity depending on the system. For AK, the main structural difference between the end-states is rigid movements of large domains and the transition also involves largely rigid-body motions of various domains. But for <inline-formula><inline-graphic xlink:href="pcbi.1003521.e103.jpg"></inline-graphic>
</inline-formula>
, even though the inspection of end-states reveals a rigid-body movement of the transport domain with respect to the trimerization domain, the structural changes along the pathway are much more complex. This fact highlights that simple ANM-based surfaces are capable of capturing significant complexity of biomolecular conformational change. For systems involving large-scale rigid-body domain motions, a few ANM modes predicted for the initial structure enable displacements far away from the starting state toward the end-state, as observed in transitions in AK and SERCA. However, if the intermediate structures involve more subtle changes, including localized motions, many more modes may be necessary (e.g. the case of LeuT). Notably, ANM-based surfaces are capable of capturing the complexity of biomolecular conformational change by inclusion of a higher number of modes. These observations validate the conventional picture of correlation between initial direction of structural change and displacements along low frequency normal modes <xref rid="pcbi.1003521-Bahar2" ref-type="bibr">[42]</xref>
–<xref rid="pcbi.1003521-Bahar4" ref-type="bibr">[44]</xref>
and will have important ramifications for pathway methods based on subsets of normal modes.</p>
<p>From a computational point of view, <italic>ANMPathway</italic>
can be used as a first approximation to simulate the real transition in all-atom models. Our simple CG model does not have any residue-specific information or side chains and for membrane protein systems we introduce a further simplification by excluding the membrane. Therefore, results from <italic>ANMPathway</italic>
should be interpreted accordingly and should only be used for a quick assessment of a likely path, or for studying relatively large-scale movements such as those of entire domains or subunits. The study demonstrates that <italic>ANMPathway</italic>
can yield very good starting point for building all-atom pathways and refining them with more advanced methods.</p>
<p>Experimental validation of any conformational transition pathway produced via computations is obviously difficult. Owing to their stability, the end-states are often amenable to direct observation by scattering or spectroscopic methods. In contrast, the short-lived intermediate structures occurring transiently along the pathway are much more difficult to detect directly. The most one can hope for is that the predicted pathway could be validated indirectly. In this regard, a particularly interesting possibility is to engineer cross-links between pairs of residues that come in close contact somewhere along the pathway but are otherwise far away from one another in the two end-states. In practice, we were able to find such pairs of residues in <inline-formula><inline-graphic xlink:href="pcbi.1003521.e104.jpg"></inline-graphic>
</inline-formula>
and in SERCA which are farther than 10 Å in both end-states, but transiently come within 5–6 Å from one another along the transition pathway. This is encouraging and perhaps experimentally probing pairs of residues satisfying this criterion could become a routine endeavor to validate and test computational pathways.</p>
</sec>
<sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title>
<supplementary-material content-type="local-data" id="pcbi.1003521.s001"><label>Figure S1</label>
<caption><p>Comparison of AD-ENM <xref rid="pcbi.1003521-Zheng1" ref-type="bibr">[55]</xref>
and <italic>ANMPathway</italic>
paths of AK.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pcbi.1003521.s001.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s002"><label>Figure S2</label>
<caption><p>Comparison of AD-ENM <xref rid="pcbi.1003521-Zheng1" ref-type="bibr">[55]</xref>
and <italic>ANMPathway</italic>
paths of LeuT.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pcbi.1003521.s002.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s003"><label>Figure S3</label>
<caption><p>Comparison of AD-ENM <xref rid="pcbi.1003521-Zheng1" ref-type="bibr">[55]</xref>
and <italic>ANMPathway</italic>
paths of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e105.jpg"></inline-graphic>
</inline-formula>
.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pcbi.1003521.s003.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s004"><label>Figure S4</label>
<caption><p>Projection of refined pathways on the smooth two-state potential defined in <xref ref-type="disp-formula" rid="pcbi.1003521.e092">Eq. (9)</xref>
on the space of order parameters used in <xref ref-type="fig" rid="pcbi-1003521-g002">Fig. (2)</xref>
.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pcbi.1003521.s004.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s005"><label>Text S1</label>
<caption><p>Protocols for all-atom molecular dynamics simulations of LeuT and <inline-formula><inline-graphic xlink:href="pcbi.1003521.e106.jpg"></inline-graphic>
</inline-formula>
.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pcbi.1003521.s005.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s006"><label>Video S1</label>
<caption><p>Transition pathway between open and closed states of AK.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s006.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s007"><label>Video S2</label>
<caption><p>Transition pathway between the outward-facing occluded state and inward-facing open state of LeuT.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s007.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s008"><label>Video S3</label>
<caption><p>Structural changes in functionally important domains of LeuT. For clarity only the EL4a (<italic>blue</italic>
), EL4b (<italic>blue</italic>
), TM1a (<italic>red</italic>
), TM1b (<italic>orange</italic>
), TM2 (<italic>purple</italic>
), TM5 (<italic>yellow</italic>
), TM6a (<italic>green</italic>
), TM6b (<italic>lime</italic>
) and TM7 (<italic>purple</italic>
) domains are shown.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s008.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s009"><label>Video S4</label>
<caption><p>Transition pathway between iOF and IF states of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e107.jpg"></inline-graphic>
</inline-formula>
. All three protomers are shown, the <italic>green</italic>
protomer undergoes transition.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s009.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s010"><label>Video S5</label>
<caption><p>Only the protomer undergoing transition in <inline-formula><inline-graphic xlink:href="pcbi.1003521.e108.jpg"></inline-graphic>
</inline-formula>
is shown with the following coloring scheme: trimerization domain, <italic>cyan</italic>
; TM3, <italic>deep sky blue</italic>
; TM6, <italic>lime</italic>
; HP1, <italic>yellow</italic>
; TM7, <italic>orange</italic>
; HP2, <italic>orange</italic>
; TM8, <italic>magenta</italic>
and the rest, <italic>gray</italic>
.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s010.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s011"><label>Video S6</label>
<caption><p>A simpler view of the transition in <inline-formula><inline-graphic xlink:href="pcbi.1003521.e109.jpg"></inline-graphic>
</inline-formula>
in terms of two blocks. The color code is: block 1, <italic>blue</italic>
; block 2, <italic>green</italic>
; rest of TM7, <italic>orange</italic>
; rest of TM8, <italic>magenta</italic>
.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s011.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s012"><label>Video S7</label>
<caption><p>Contact formation between Val62 (<italic>blue</italic>
sphere) and Ala353 (<italic>red</italic>
sphere) of chain C during the transition between iOF and IF states of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e110.jpg"></inline-graphic>
</inline-formula>
.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s012.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s013"><label>Video S8</label>
<caption><p>Contact formation between Gln220 (<italic>blue</italic>
sphere) and Met385 (<italic>red</italic>
sphere) of chain C during the transition between iOF and IF states of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e111.jpg"></inline-graphic>
</inline-formula>
.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s013.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s014"><label>Video S9</label>
<caption><p>Contact formation between Val274 (<italic>blue</italic>
sphere) and Ala391 (<italic>red</italic>
sphere) of chain C during the transition between iOF and IF states of <inline-formula><inline-graphic xlink:href="pcbi.1003521.e112.jpg"></inline-graphic>
</inline-formula>
.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s014.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s015"><label>Video S10</label>
<caption><p>Transtion pathway between <inline-formula><inline-graphic xlink:href="pcbi.1003521.e113.jpg"></inline-graphic>
</inline-formula>
and E1.ATP states of SERCA.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s015.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s016"><label>Video S11</label>
<caption><p>Contact formation between Leu49 (<italic>cyan</italic>
sphere) and Lys246 (<italic>purple</italic>
sphere) during the transition between <inline-formula><inline-graphic xlink:href="pcbi.1003521.e114.jpg"></inline-graphic>
</inline-formula>
and E1.ATP states of SERCA.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s016.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s017"><label>Video S12</label>
<caption><p>Contact formation between Trp50 (<italic>cyan</italic>
sphere) and Thr316 (<italic>purple</italic>
sphere) during the transition between <inline-formula><inline-graphic xlink:href="pcbi.1003521.e115.jpg"></inline-graphic>
</inline-formula>
and E1.ATP states of SERCA.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s017.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1003521.s018"><label>Video S13</label>
<caption><p>Contact formation between Ile54 (<italic>cyan</italic>
sphere) and Pro312 (<italic>purple</italic>
sphere) during the transition between <inline-formula><inline-graphic xlink:href="pcbi.1003521.e116.jpg"></inline-graphic>
</inline-formula>
and E1.ATP states of SERCA.</p>
<p>(MP4)</p>
</caption>
<media xlink:href="pcbi.1003521.s018.mp4"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><ref-list><title>References</title>
<ref id="pcbi.1003521-HenzlerWildman1"><label>1</label>
<mixed-citation publication-type="journal"><name><surname>Henzler-Wildman</surname>
<given-names>KA</given-names>
</name>
, <name><surname>Thai</surname>
<given-names>V</given-names>
</name>
, <name><surname>Lei</surname>
<given-names>M</given-names>
</name>
, <name><surname>Ott</surname>
<given-names>M</given-names>
</name>
, <name><surname>Wolf-Watz</surname>
<given-names>M</given-names>
</name>
, <etal>et al</etal>
(<year>2007</year>
) <article-title>Intrinsic motions along an enzymatic reaction trajectory</article-title>
. <source>Nature</source>
<volume>450</volume>
: <fpage>838</fpage>
–<lpage>844</lpage>
.<pub-id pub-id-type="pmid">18026086</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Lei1"><label>2</label>
<mixed-citation publication-type="journal"><name><surname>Lei</surname>
<given-names>M</given-names>
</name>
, <name><surname>Velos</surname>
<given-names>J</given-names>
</name>
, <name><surname>Gardino</surname>
<given-names>A</given-names>
</name>
, <name><surname>Kivenson</surname>
<given-names>A</given-names>
</name>
, <name><surname>Karplus</surname>
<given-names>M</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Segmented transition pathway of the signaling protein nitrogen regulatory protein C</article-title>
. <source>J Mol Biol</source>
<volume>392</volume>
: <fpage>823</fpage>
–<lpage>836</lpage>
.<pub-id pub-id-type="pmid">19576227</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Lacroix1"><label>3</label>
<mixed-citation publication-type="journal"><name><surname>Lacroix</surname>
<given-names>JJ</given-names>
</name>
, <name><surname>Pless</surname>
<given-names>SA</given-names>
</name>
, <name><surname>Maragliano</surname>
<given-names>L</given-names>
</name>
, <name><surname>Campos</surname>
<given-names>FV</given-names>
</name>
, <name><surname>Galpin</surname>
<given-names>JD</given-names>
</name>
, <etal>et al</etal>
(<year>2012</year>
) <article-title>Intermediate state trapping of a voltage sensor</article-title>
. <source>J Gen Physiol</source>
<volume>140</volume>
: <fpage>635</fpage>
–<lpage>652</lpage>
.<pub-id pub-id-type="pmid">23183699</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Allen1"><label>4</label>
<mixed-citation publication-type="book">Allen MP, Tildesley DP (1987) Computer Simulation of Liquids. Oxford University Press: Oxford.</mixed-citation>
</ref>
<ref id="pcbi.1003521-Frenkel1"><label>5</label>
<mixed-citation publication-type="book">Frenkel M, Smit B (2002) Understanding Molecular Simulation: From Algorithms to Applications. Academic Press.</mixed-citation>
</ref>
<ref id="pcbi.1003521-Shaw1"><label>6</label>
<mixed-citation publication-type="journal"><name><surname>Shaw</surname>
<given-names>DE</given-names>
</name>
, <name><surname>Deneroff</surname>
<given-names>MM</given-names>
</name>
, <name><surname>Dror</surname>
<given-names>RO</given-names>
</name>
, <name><surname>Kuskin</surname>
<given-names>JS</given-names>
</name>
, <name><surname>Larson</surname>
<given-names>RH</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>Anton, a special-purpose machine for molecular dynamics simulation</article-title>
. <source>Commun ACM</source>
<volume>51</volume>
: <fpage>91</fpage>
–<lpage>97</lpage>
.</mixed-citation>
</ref>
<ref id="pcbi.1003521-Chandler1"><label>7</label>
<mixed-citation publication-type="journal"><name><surname>Chandler</surname>
<given-names>D</given-names>
</name>
(<year>1978</year>
) <article-title>Statistical mechanics of isomerization dynamics in liquids and the transition state approximation</article-title>
. <source>J Chem Phys</source>
<volume>68</volume>
: <fpage>2959</fpage>
–<lpage>2970</lpage>
.</mixed-citation>
</ref>
<ref id="pcbi.1003521-Dellago1"><label>8</label>
<mixed-citation publication-type="journal"><name><surname>Dellago</surname>
<given-names>C</given-names>
</name>
, <name><surname>Bolhuis</surname>
<given-names>PG</given-names>
</name>
, <name><surname>Csajka</surname>
<given-names>FS</given-names>
</name>
, <name><surname>Chandler</surname>
<given-names>D</given-names>
</name>
(<year>1998</year>
) <article-title>Transition path sampling and the calculation of rate constants</article-title>
. <source>J Chem Phys</source>
<volume>108</volume>
: <fpage>1964</fpage>
–<lpage>1977</lpage>
.</mixed-citation>
</ref>
<ref id="pcbi.1003521-Bolhuis1"><label>9</label>
<mixed-citation publication-type="journal"><name><surname>Bolhuis</surname>
<given-names>PG</given-names>
</name>
, <name><surname>Chandler</surname>
<given-names>D</given-names>
</name>
, <name><surname>Dellago</surname>
<given-names>C</given-names>
</name>
, <name><surname>Geissler</surname>
<given-names>PL</given-names>
</name>
(<year>2002</year>
) <article-title>Transition path sampling: throwing ropes over rough mountain passes, in the dark</article-title>
. <source>Annu Rev Phys Chem</source>
<volume>53</volume>
: <fpage>291</fpage>
–<lpage>318</lpage>
.<pub-id pub-id-type="pmid">11972010</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Du1"><label>10</label>
<mixed-citation publication-type="journal"><name><surname>Du</surname>
<given-names>WN</given-names>
</name>
, <name><surname>Marino</surname>
<given-names>KA</given-names>
</name>
, <name><surname>Bolhuis</surname>
<given-names>PG</given-names>
</name>
(<year>2011</year>
) <article-title>Multiple state transition interface sampling of alanine dipeptide in explicit solvent</article-title>
. <source>J Chem Phys</source>
<volume>135</volume>
: <fpage>145102</fpage>
.<pub-id pub-id-type="pmid">22010733</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Elber1"><label>11</label>
<mixed-citation publication-type="journal"><name><surname>Elber</surname>
<given-names>R</given-names>
</name>
, <name><surname>Karplus</surname>
<given-names>M</given-names>
</name>
(<year>1987</year>
) <article-title>A method for determining reaction paths in large molecules - application to myoglobin</article-title>
. <source>Chem Phys Lett</source>
<volume>139</volume>
: <fpage>375</fpage>
–<lpage>80</lpage>
.</mixed-citation>
</ref>
<ref id="pcbi.1003521-Jonsson1"><label>12</label>
<mixed-citation publication-type="other">Jonsson H, Mills G, Jacobsen KW (1998) In: Berne BJ, Ciccotti G, Coker DF, editors, Classical and Quantum Dynamics in Condensed Phase Simulations, World Scientific. pp. 385–404.</mixed-citation>
</ref>
<ref id="pcbi.1003521-Elber2"><label>13</label>
<mixed-citation publication-type="journal"><name><surname>Elber</surname>
<given-names>R</given-names>
</name>
(<year>2005</year>
) <article-title>Long-timescale simulation methods</article-title>
. <source>Curr Opin Struct Biol</source>
<volume>15</volume>
: <fpage>151</fpage>
–<lpage>156</lpage>
.<pub-id pub-id-type="pmid">15837172</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Maragliano1"><label>14</label>
<mixed-citation publication-type="journal"><name><surname>Maragliano</surname>
<given-names>L</given-names>
</name>
, <name><surname>Fischer</surname>
<given-names>A</given-names>
</name>
, <name><surname>Vanden-Eijnden</surname>
<given-names>E</given-names>
</name>
, <name><surname>Ciccotti</surname>
<given-names>G</given-names>
</name>
(<year>2006</year>
) <article-title>String method in collective variables: Minimum free energy paths and isocommittor surfaces</article-title>
. <source>J Chem Phys</source>
<volume>125</volume>
: <fpage>024106</fpage>
.</mixed-citation>
</ref>
<ref id="pcbi.1003521-Pan1"><label>15</label>
<mixed-citation publication-type="journal"><name><surname>Pan</surname>
<given-names>AC</given-names>
</name>
, <name><surname>Sezer</surname>
<given-names>D</given-names>
</name>
, <name><surname>Roux</surname>
<given-names>B</given-names>
</name>
(<year>2008</year>
) <article-title>Finding transition pathways using the string method with swarms of trajectories</article-title>
. <source>J Phys Chem B</source>
<volume>112</volume>
: <fpage>3432</fpage>
–<lpage>3440</lpage>
.<pub-id pub-id-type="pmid">18290641</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Gan1"><label>16</label>
<mixed-citation publication-type="journal"><name><surname>Gan</surname>
<given-names>W</given-names>
</name>
, <name><surname>Yang</surname>
<given-names>S</given-names>
</name>
, <name><surname>Roux</surname>
<given-names>B</given-names>
</name>
(<year>2009</year>
) <article-title>Atomistic View of the Conformational Activation of Src Kinase Using the String Method with Swarms-of-Trajectories</article-title>
. <source>Biophys J</source>
<volume>97</volume>
: <fpage>L8</fpage>
–<lpage>L10</lpage>
.<pub-id pub-id-type="pmid">19686639</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Vashisth1"><label>17</label>
<mixed-citation publication-type="journal"><name><surname>Vashisth</surname>
<given-names>H</given-names>
</name>
, <name><surname>Maragliano</surname>
<given-names>L</given-names>
</name>
, <name><surname>Abrams</surname>
<given-names>CF</given-names>
</name>
(<year>2012</year>
) <article-title>“DFG-flip” in the insulin receptor kinase is facilitated by a helical intermediate state of the activation loop</article-title>
. <source>Biophys J</source>
<volume>102</volume>
: <fpage>1979</fpage>
–<lpage>1987</lpage>
.<pub-id pub-id-type="pmid">22768955</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Vashisth2"><label>18</label>
<mixed-citation publication-type="journal"><name><surname>Vashisth</surname>
<given-names>H</given-names>
</name>
, <name><surname>Abrams</surname>
<given-names>CF</given-names>
</name>
(<year>2013</year>
) <article-title>All-atom structural models of insulin binding to the insulin receptor in the presence of a tandem hormone-binding element</article-title>
. <source>Proteins</source>
<volume>81</volume>
: <fpage>1017</fpage>
–<lpage>1030</lpage>
.<pub-id pub-id-type="pmid">23348915</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Matsunaga1"><label>19</label>
<mixed-citation publication-type="journal"><name><surname>Matsunaga</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Fujisaki</surname>
<given-names>H</given-names>
</name>
, <name><surname>Terada</surname>
<given-names>T</given-names>
</name>
, <name><surname>Furuta</surname>
<given-names>T</given-names>
</name>
, <name><surname>Moritsugu</surname>
<given-names>K</given-names>
</name>
, <etal>et al</etal>
(<year>2012</year>
) <article-title>Minimum free energy path of ligand-induced transition in adenylate kinase</article-title>
. <source>PLoS Comput Biol</source>
<volume>8</volume>
: <fpage>e1002555</fpage>
.<pub-id pub-id-type="pmid">22685395</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Stober1"><label>20</label>
<mixed-citation publication-type="journal"><name><surname>Stober</surname>
<given-names>ST</given-names>
</name>
, <name><surname>Abrams</surname>
<given-names>CF</given-names>
</name>
(<year>2012</year>
) <article-title>Energetics and mechanism of the normal-to-amyloidogenic isomerization of 2-microglobulin: On-the-fly string method calculations</article-title>
. <source>J Phys Chem B</source>
<volume>116</volume>
: <fpage>9371</fpage>
–<lpage>9375</lpage>
.<pub-id pub-id-type="pmid">22793795</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Jo1"><label>21</label>
<mixed-citation publication-type="journal"><name><surname>Jo</surname>
<given-names>S</given-names>
</name>
, <name><surname>Rui</surname>
<given-names>H</given-names>
</name>
, <name><surname>Lim</surname>
<given-names>JB</given-names>
</name>
, <name><surname>Klauda</surname>
<given-names>JB</given-names>
</name>
, <name><surname>Im</surname>
<given-names>W</given-names>
</name>
(<year>2010</year>
) <article-title>Cholesterol flip-flop: Insights from free energy simulation studies</article-title>
. <source>J Phys Chem B</source>
<volume>114</volume>
: <fpage>13342</fpage>
–<lpage>13348</lpage>
.<pub-id pub-id-type="pmid">20923227</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Ovchinnikov1"><label>22</label>
<mixed-citation publication-type="journal"><name><surname>Ovchinnikov</surname>
<given-names>V</given-names>
</name>
, <name><surname>Karplus</surname>
<given-names>M</given-names>
</name>
, <name><surname>Vanden-Eijnden</surname>
<given-names>E</given-names>
</name>
(<year>2011</year>
) <article-title>Free energy of conformational transition paths in biomolecules: The string method and its application to myosin VI</article-title>
. <source>J Chem Phys</source>
<volume>134</volume>
: <fpage>085103</fpage>
.<pub-id pub-id-type="pmid">21361558</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Kirmizialtin1"><label>23</label>
<mixed-citation publication-type="journal"><name><surname>Kirmizialtin</surname>
<given-names>S</given-names>
</name>
, <name><surname>Nguyen</surname>
<given-names>V</given-names>
</name>
, <name><surname>Johnson</surname>
<given-names>KA</given-names>
</name>
, <name><surname>Elber</surname>
<given-names>R</given-names>
</name>
(<year>2012</year>
) <article-title>How conformational dynamics of DNA polymerase select correct substrates: experiments and simulations</article-title>
. <source>Structure</source>
<volume>20</volume>
: <fpage>618</fpage>
–<lpage>627</lpage>
.<pub-id pub-id-type="pmid">22483109</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Huber1"><label>24</label>
<mixed-citation publication-type="journal"><name><surname>Huber</surname>
<given-names>G</given-names>
</name>
, <name><surname>Kim</surname>
<given-names>S</given-names>
</name>
(<year>1996</year>
) <article-title>Weighted-ensemble Brownian dynamics simulations for protein association reactions</article-title>
. <source>Biophys J</source>
<volume>70</volume>
: <fpage>97</fpage>
–<lpage>110</lpage>
.<pub-id pub-id-type="pmid">8770190</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Zhang1"><label>25</label>
<mixed-citation publication-type="journal"><name><surname>Zhang</surname>
<given-names>BW</given-names>
</name>
, <name><surname>Jasnow</surname>
<given-names>D</given-names>
</name>
, <name><surname>Zuckerman</surname>
<given-names>DM</given-names>
</name>
(<year>2010</year>
) <article-title>The “weighted ensemble” path sampling method is statistically exact for a broad class of stochastic processes and binning procedures</article-title>
. <source>J Chem Phys</source>
<volume>132</volume>
: <fpage>054107</fpage>
.<pub-id pub-id-type="pmid">20136305</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Bhatt1"><label>26</label>
<mixed-citation publication-type="journal"><name><surname>Bhatt</surname>
<given-names>D</given-names>
</name>
, <name><surname>Zhang</surname>
<given-names>BW</given-names>
</name>
, <name><surname>Zuckerman</surname>
<given-names>DM</given-names>
</name>
(<year>2010</year>
) <article-title>Steady-state simulations using weighted ensemble path sampling</article-title>
. <source>J Chem Phys</source>
<volume>133</volume>
: <fpage>014110</fpage>
.<pub-id pub-id-type="pmid">20614962</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Bhatt2"><label>27</label>
<mixed-citation publication-type="journal"><name><surname>Bhatt</surname>
<given-names>D</given-names>
</name>
, <name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
(<year>2012</year>
) <article-title>An adaptive weighted ensemble procedure for efficient computation of free energies and first passage rates</article-title>
. <source>J Chem Phys</source>
<volume>137</volume>
: <fpage>104101</fpage>
.<pub-id pub-id-type="pmid">22979844</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Woolf1"><label>28</label>
<mixed-citation publication-type="journal"><name><surname>Woolf</surname>
<given-names>TB</given-names>
</name>
(<year>1998</year>
) <article-title>Path corrected functionals of stochastic trajectories: towards relative free energy and reaction coordinate calculations</article-title>
. <source>Chem Phys Lett</source>
<volume>289</volume>
: <fpage>433</fpage>
–<lpage>441</lpage>
.</mixed-citation>
</ref>
<ref id="pcbi.1003521-Zuckerman1"><label>29</label>
<mixed-citation publication-type="journal"><name><surname>Zuckerman</surname>
<given-names>DM</given-names>
</name>
, <name><surname>Woolf</surname>
<given-names>TB</given-names>
</name>
(<year>1999</year>
) <article-title>Dynamic reaction paths and rates through importance-sampled stochastic dynamics</article-title>
. <source>J Chem Phys</source>
<volume>111</volume>
: <fpage>9475</fpage>
–<lpage>9484</lpage>
.</mixed-citation>
</ref>
<ref id="pcbi.1003521-Perilla1"><label>30</label>
<mixed-citation publication-type="journal"><name><surname>Perilla</surname>
<given-names>JR</given-names>
</name>
, <name><surname>Beckstein</surname>
<given-names>O</given-names>
</name>
, <name><surname>Denning</surname>
<given-names>EJ</given-names>
</name>
, <name><surname>Woolf</surname>
<given-names>TB</given-names>
</name>
(<year>2011</year>
) <article-title>Computing ensembles of transitions from stable states: Dynamic importance sampling</article-title>
. <source>J Comput Chem</source>
<volume>32</volume>
: <fpage>196</fpage>
–<lpage>209</lpage>
.<pub-id pub-id-type="pmid">21132840</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Grubmller1"><label>31</label>
<mixed-citation publication-type="journal"><name><surname>Grubmüller</surname>
<given-names>H</given-names>
</name>
(<year>1995</year>
) <article-title>Predicting slow structural transitions in macromolecular systems: Conformational flooding</article-title>
. <source>Phys Rev E</source>
<volume>52</volume>
: <fpage>2893</fpage>
–<lpage>2906</lpage>
.</mixed-citation>
</ref>
<ref id="pcbi.1003521-Laio1"><label>32</label>
<mixed-citation publication-type="journal"><name><surname>Laio</surname>
<given-names>A</given-names>
</name>
, <name><surname>Parrinello</surname>
<given-names>M</given-names>
</name>
(<year>2002</year>
) <article-title>Escaping free-energy minima</article-title>
. <source>Proc Natl Acad Sci U S A</source>
<volume>99</volume>
: <fpage>12562</fpage>
–<lpage>12566</lpage>
.<pub-id pub-id-type="pmid">12271136</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Hamelberg1"><label>33</label>
<mixed-citation publication-type="journal"><name><surname>Hamelberg</surname>
<given-names>D</given-names>
</name>
, <name><surname>Mongan</surname>
<given-names>J</given-names>
</name>
, <name><surname>McCammon</surname>
<given-names>JA</given-names>
</name>
(<year>2004</year>
) <article-title>Accelerated molecular dynamics: A promising and efficient simulation method for biomolecules</article-title>
. <source>J Chem Phys</source>
<volume>120</volume>
: <fpage>11919</fpage>
–<lpage>29</lpage>
.<pub-id pub-id-type="pmid">15268227</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Seeliger1"><label>34</label>
<mixed-citation publication-type="journal"><name><surname>Seeliger</surname>
<given-names>D</given-names>
</name>
, <name><surname>Haas</surname>
<given-names>J</given-names>
</name>
, <name><surname>de Groot</surname>
<given-names>BL</given-names>
</name>
(<year>2007</year>
) <article-title>Geometry-based sampling of conformational transitions in proteins</article-title>
. <source>Structure</source>
<volume>15</volume>
: <fpage>1482</fpage>
–<lpage>1492</lpage>
.<pub-id pub-id-type="pmid">17997973</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Bahar1"><label>35</label>
<mixed-citation publication-type="journal"><name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
, <name><surname>Atilgan</surname>
<given-names>AR</given-names>
</name>
, <name><surname>Erman</surname>
<given-names>B</given-names>
</name>
(<year>1997</year>
) <article-title>Direct evaluation of thermal fluctuations in proteins using a single-parameter harmonic potential</article-title>
. <source>Fold Des</source>
<volume>2</volume>
: <fpage>173</fpage>
–<lpage>181</lpage>
.<pub-id pub-id-type="pmid">9218955</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Haliloglu1"><label>36</label>
<mixed-citation publication-type="journal"><name><surname>Haliloglu</surname>
<given-names>T</given-names>
</name>
, <name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
, <name><surname>Erman</surname>
<given-names>B</given-names>
</name>
(<year>1997</year>
) <article-title>Gaussian dynamics of folded proteins</article-title>
. <source>Phys Rev Lett</source>
<volume>79</volume>
: <fpage>3090</fpage>
–<lpage>3093</lpage>
.</mixed-citation>
</ref>
<ref id="pcbi.1003521-Tirion1"><label>37</label>
<mixed-citation publication-type="journal"><name><surname>Tirion</surname>
<given-names>MM</given-names>
</name>
(<year>1996</year>
) <article-title>Large amplitude elastic motions in proteins from a single-parameter, atomic analysis</article-title>
. <source>Phys Rev Lett</source>
<volume>77</volume>
: <fpage>1905</fpage>
–<lpage>1908</lpage>
.<pub-id pub-id-type="pmid">10063201</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Atilgan1"><label>38</label>
<mixed-citation publication-type="journal"><name><surname>Atilgan</surname>
<given-names>A</given-names>
</name>
, <name><surname>Durell</surname>
<given-names>S</given-names>
</name>
, <name><surname>Jernigan</surname>
<given-names>R</given-names>
</name>
, <name><surname>Demirel</surname>
<given-names>M</given-names>
</name>
, <name><surname>Keskin</surname>
<given-names>O</given-names>
</name>
, <etal>et al</etal>
(<year>2001</year>
) <article-title>Anisotropy of fluctuation dynamics of proteins with an elastic network model</article-title>
. <source>Biophys J</source>
<volume>80</volume>
: <fpage>505</fpage>
–<lpage>515</lpage>
.<pub-id pub-id-type="pmid">11159421</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Thomas1"><label>39</label>
<mixed-citation publication-type="journal"><name><surname>Thomas</surname>
<given-names>A</given-names>
</name>
, <name><surname>Hinsen</surname>
<given-names>K</given-names>
</name>
, <name><surname>Field</surname>
<given-names>MJ</given-names>
</name>
, <name><surname>Perahia</surname>
<given-names>D</given-names>
</name>
(<year>1999</year>
) <article-title>Tertiary and quaternary conformational changes in aspartate transcarbamylase: a normal mode study</article-title>
. <source>Proteins</source>
<volume>34</volume>
: <fpage>96</fpage>
–<lpage>112</lpage>
.<pub-id pub-id-type="pmid">10336386</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Tama1"><label>40</label>
<mixed-citation publication-type="journal"><name><surname>Tama</surname>
<given-names>F</given-names>
</name>
, <name><surname>Sanejouand</surname>
<given-names>YH</given-names>
</name>
(<year>2001</year>
) <article-title>Conformational change of proteins arising from normal mode calculations</article-title>
. <source>Protein Eng</source>
<volume>14</volume>
: <fpage>1</fpage>
–<lpage>6</lpage>
.<pub-id pub-id-type="pmid">11287673</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Reuter1"><label>41</label>
<mixed-citation publication-type="journal"><name><surname>Reuter</surname>
<given-names>N</given-names>
</name>
, <name><surname>Hinsen</surname>
<given-names>K</given-names>
</name>
, <name><surname>Lacapre</surname>
<given-names>JJ</given-names>
</name>
(<year>2003</year>
) <article-title>Transconformations of the SERCA1 Ca-ATPase: A normal mode study</article-title>
. <source>Biophys J</source>
<volume>85</volume>
: <fpage>2186</fpage>
–<lpage>2197</lpage>
.<pub-id pub-id-type="pmid">14507684</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Bahar2"><label>42</label>
<mixed-citation publication-type="journal"><name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
, <name><surname>Rader</surname>
<given-names>A</given-names>
</name>
(<year>2005</year>
) <article-title>Coarse-grained normal mode analysis in structural biology</article-title>
. <source>Curr Opin Struct Biol</source>
<volume>15</volume>
: <fpage>586</fpage>
–<lpage>592</lpage>
.<pub-id pub-id-type="pmid">16143512</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Bahar3"><label>43</label>
<mixed-citation publication-type="journal"><name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
, <name><surname>Lezon</surname>
<given-names>TR</given-names>
</name>
, <name><surname>Yang</surname>
<given-names>LW</given-names>
</name>
, <name><surname>Eyal</surname>
<given-names>E</given-names>
</name>
(<year>2010</year>
) <article-title>Global dynamics of proteins: Bridging between structure and function</article-title>
. <source>Annu Rev Biophys</source>
<volume>39</volume>
: <fpage>23</fpage>
–<lpage>42</lpage>
.<pub-id pub-id-type="pmid">20192781</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Bahar4"><label>44</label>
<mixed-citation publication-type="journal"><name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
, <name><surname>Lezon</surname>
<given-names>TR</given-names>
</name>
, <name><surname>Bakan</surname>
<given-names>A</given-names>
</name>
, <name><surname>Shrivastava</surname>
<given-names>IH</given-names>
</name>
(<year>2010</year>
) <article-title>Normal mode analysis of biomolecular structures: Functional mechanisms of membrane proteins</article-title>
. <source>Chem Rev</source>
<volume>110</volume>
: <fpage>1463</fpage>
–<lpage>1497</lpage>
.<pub-id pub-id-type="pmid">19785456</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Gur1"><label>45</label>
<mixed-citation publication-type="journal"><name><surname>Gur</surname>
<given-names>M</given-names>
</name>
, <name><surname>Zomot</surname>
<given-names>E</given-names>
</name>
, <name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
(<year>2013</year>
) <article-title>Global motions exhibited by proteins in micro- to milliseconds simulations concur with anisotropic network model predictions</article-title>
. <source>J Chem Phys</source>
<volume>139</volume>
: <fpage>121912</fpage>
.<pub-id pub-id-type="pmid">24089724</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Kim1"><label>46</label>
<mixed-citation publication-type="journal"><name><surname>Kim</surname>
<given-names>MK</given-names>
</name>
, <name><surname>Chirikjian</surname>
<given-names>GS</given-names>
</name>
, <name><surname>Jernigan</surname>
<given-names>RL</given-names>
</name>
(<year>2002</year>
) <article-title>Elastic models of conformational transitions in macromolecules</article-title>
. <source>J Mol Graph Model</source>
<volume>21</volume>
: <fpage>151</fpage>
–<lpage>160</lpage>
.<pub-id pub-id-type="pmid">12398345</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Kim2"><label>47</label>
<mixed-citation publication-type="journal"><name><surname>Kim</surname>
<given-names>MK</given-names>
</name>
, <name><surname>Jernigan</surname>
<given-names>RL</given-names>
</name>
, <name><surname>Chirikjian</surname>
<given-names>GS</given-names>
</name>
(<year>2002</year>
) <article-title>Efficient generation of feasible pathways for protein conformational transitions</article-title>
. <source>Biophys J</source>
<volume>83</volume>
: <fpage>1620</fpage>
–<lpage>1630</lpage>
.<pub-id pub-id-type="pmid">12202386</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Schuyler1"><label>48</label>
<mixed-citation publication-type="journal"><name><surname>Schuyler</surname>
<given-names>AD</given-names>
</name>
, <name><surname>Chirikjian</surname>
<given-names>GS</given-names>
</name>
(<year>2005</year>
) <article-title>Efficient determination of low-frequency normal modes of large protein structures by cluster-NMA</article-title>
. <source>J Mol Graph Model</source>
<volume>24</volume>
: <fpage>46</fpage>
–<lpage>58</lpage>
.<pub-id pub-id-type="pmid">15990344</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Kim3"><label>49</label>
<mixed-citation publication-type="journal"><name><surname>Kim</surname>
<given-names>MK</given-names>
</name>
, <name><surname>Jernigan</surname>
<given-names>RL</given-names>
</name>
, <name><surname>Chirikjian</surname>
<given-names>GS</given-names>
</name>
(<year>2005</year>
) <article-title>Rigid-cluster models of conformational transitions in macromolecular machines and assemblies</article-title>
. <source>Biophys J</source>
<volume>89</volume>
: <fpage>43</fpage>
–<lpage>55</lpage>
.<pub-id pub-id-type="pmid">15833998</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Schuyler2"><label>50</label>
<mixed-citation publication-type="journal"><name><surname>Schuyler</surname>
<given-names>AD</given-names>
</name>
, <name><surname>Jernigan</surname>
<given-names>RL</given-names>
</name>
, <name><surname>Qasba</surname>
<given-names>PK</given-names>
</name>
, <name><surname>Ramakrishnan</surname>
<given-names>B</given-names>
</name>
, <name><surname>Chirikjian</surname>
<given-names>GS</given-names>
</name>
(<year>2009</year>
) <article-title>Iterative cluster- NMA: A tool for generating conformational transitions in proteins</article-title>
. <source>Proteins</source>
<volume>74</volume>
: <fpage>760</fpage>
–<lpage>776</lpage>
.<pub-id pub-id-type="pmid">18712827</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Miyashita1"><label>51</label>
<mixed-citation publication-type="journal"><name><surname>Miyashita</surname>
<given-names>O</given-names>
</name>
, <name><surname>Onuchic</surname>
<given-names>JN</given-names>
</name>
, <name><surname>Wolynes</surname>
<given-names>PG</given-names>
</name>
(<year>2003</year>
) <article-title>Nonlinear elasticity, proteinquakes, and the energy landscapes of functional transitions in proteins</article-title>
. <source>Proc Natl Acad Sci U S A</source>
<volume>100</volume>
: <fpage>12570</fpage>
–<lpage>12575</lpage>
.<pub-id pub-id-type="pmid">14566052</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Maragakis1"><label>52</label>
<mixed-citation publication-type="journal"><name><surname>Maragakis</surname>
<given-names>P</given-names>
</name>
, <name><surname>Karplus</surname>
<given-names>M</given-names>
</name>
(<year>2005</year>
) <article-title>Large amplitude conformational change in proteins explored with a plastic network model: Adenylate kinase</article-title>
. <source>J Mol Biol</source>
<volume>352</volume>
: <fpage>807</fpage>
–<lpage>822</lpage>
.<pub-id pub-id-type="pmid">16139299</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Yang1"><label>53</label>
<mixed-citation publication-type="journal"><name><surname>Yang</surname>
<given-names>Z</given-names>
</name>
, <name><surname>Majek</surname>
<given-names>P</given-names>
</name>
, <name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
(<year>2009</year>
) <article-title>Allosteric transitions of supramolecular systems explored by network models: Application to chaperonin GroEL</article-title>
. <source>PLoS Comput Biol</source>
<volume>5</volume>
: <fpage>e1000360</fpage>
.<pub-id pub-id-type="pmid">19381265</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Best1"><label>54</label>
<mixed-citation publication-type="journal"><name><surname>Best</surname>
<given-names>RB</given-names>
</name>
, <name><surname>Chen</surname>
<given-names>YG</given-names>
</name>
, <name><surname>Hummer</surname>
<given-names>G</given-names>
</name>
(<year>2005</year>
) <article-title>Slow protein conformational dynamics from multiple experimental structures: the helix/sheet transition of arc repressor</article-title>
. <source>Structure</source>
<volume>13</volume>
: <fpage>1755</fpage>
–<lpage>1763</lpage>
.<pub-id pub-id-type="pmid">16338404</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Zheng1"><label>55</label>
<mixed-citation publication-type="journal"><name><surname>Zheng</surname>
<given-names>W</given-names>
</name>
, <name><surname>Brooks</surname>
<given-names>BR</given-names>
</name>
, <name><surname>Hummer</surname>
<given-names>G</given-names>
</name>
(<year>2007</year>
) <article-title>Protein conformational transitions explored by mixed elastic network models</article-title>
. <source>Proteins</source>
<volume>69</volume>
: <fpage>43</fpage>
–<lpage>57</lpage>
.<pub-id pub-id-type="pmid">17596847</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Zhu1"><label>56</label>
<mixed-citation publication-type="journal"><name><surname>Zhu</surname>
<given-names>F</given-names>
</name>
, <name><surname>Hummer</surname>
<given-names>G</given-names>
</name>
(<year>2009</year>
) <article-title>Gating transition of pentameric ligand-gated ion channels</article-title>
. <source>Biophys J</source>
<volume>97</volume>
: <fpage>2456</fpage>
–<lpage>2463</lpage>
.<pub-id pub-id-type="pmid">19883588</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Yang2"><label>57</label>
<mixed-citation publication-type="journal"><name><surname>Yang</surname>
<given-names>S</given-names>
</name>
, <name><surname>Roux</surname>
<given-names>B</given-names>
</name>
(<year>2008</year>
) <article-title>Src kinase conformational activation: Thermodynamics, pathways, and mechanisms</article-title>
. <source>PLoS Comput Biol</source>
<volume>4</volume>
: <fpage>e1000047</fpage>
.<pub-id pub-id-type="pmid">18369437</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-G1"><label>58</label>
<mixed-citation publication-type="journal"><name><surname>Gō</surname>
<given-names>N</given-names>
</name>
, <name><surname>Abe</surname>
<given-names>H</given-names>
</name>
(<year>1981</year>
) <article-title>Noninteracting local-structure model of folding and unfolding transition in globular proteins. i. formulation</article-title>
. <source>Biopolymers</source>
<volume>20</volume>
: <fpage>991</fpage>
–<lpage>1011</lpage>
.<pub-id pub-id-type="pmid">7225531</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Abe1"><label>59</label>
<mixed-citation publication-type="journal"><name><surname>Abe</surname>
<given-names>H</given-names>
</name>
, <name><surname>Gō</surname>
<given-names>N</given-names>
</name>
(<year>1981</year>
) <article-title>Noninteracting local-structure model of folding and unfolding transition in globular proteins. ii. applications to two-diemnsional lattice proteins</article-title>
. <source>Biopolymers</source>
<volume>20</volume>
: <fpage>1013</fpage>
–<lpage>1031</lpage>
.<pub-id pub-id-type="pmid">7225529</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Chu1"><label>60</label>
<mixed-citation publication-type="journal"><name><surname>Chu</surname>
<given-names>JW</given-names>
</name>
, <name><surname>Voth</surname>
<given-names>GA</given-names>
</name>
(<year>2007</year>
) <article-title>Coarse-grained free energy functions for studying protein conformational changes: A double-well network model</article-title>
. <source>Biophys J</source>
<volume>93</volume>
: <fpage>3860</fpage>
–<lpage>3871</lpage>
.<pub-id pub-id-type="pmid">17704151</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Franklin1"><label>61</label>
<mixed-citation publication-type="journal"><name><surname>Franklin</surname>
<given-names>J</given-names>
</name>
, <name><surname>Koehl</surname>
<given-names>P</given-names>
</name>
, <name><surname>Doniach</surname>
<given-names>S</given-names>
</name>
, <name><surname>Delarue</surname>
<given-names>M</given-names>
</name>
(<year>2007</year>
) <article-title>MinActionPath: maximum likelihood trajectory for large-scale structural transitions in a coarse-grained locally harmonic energy landscape</article-title>
. <source>Nucleic Acids Res</source>
<volume>35</volume>
: <fpage>W477</fpage>
–<lpage>W482</lpage>
.<pub-id pub-id-type="pmid">17545201</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Mller1"><label>62</label>
<mixed-citation publication-type="journal"><name><surname>Müller</surname>
<given-names>C</given-names>
</name>
, <name><surname>Schlauderer</surname>
<given-names>G</given-names>
</name>
, <name><surname>Reinstein</surname>
<given-names>J</given-names>
</name>
, <name><surname>Schulz</surname>
<given-names>G</given-names>
</name>
(<year>1996</year>
) <article-title>Adenylate kinase motions during catalysis: an energetic counterweight balancing substrate binding</article-title>
. <source>Structure</source>
<volume>4</volume>
: <fpage>147</fpage>
–<lpage>156</lpage>
.<pub-id pub-id-type="pmid">8805521</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Mller2"><label>63</label>
<mixed-citation publication-type="journal"><name><surname>Müller</surname>
<given-names>CW</given-names>
</name>
, <name><surname>Schulz</surname>
<given-names>GE</given-names>
</name>
(<year>1992</year>
) <article-title>Structure of the complex between adenylate kinase from escherichia coli and the inhibitor ap5a refined at 1.9 Å resolution: A model for a catalytic transition state</article-title>
. <source>J Mol Biol</source>
<volume>224</volume>
: <fpage>159</fpage>
–<lpage>177</lpage>
.<pub-id pub-id-type="pmid">1548697</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Gur2"><label>64</label>
<mixed-citation publication-type="journal"><name><surname>Gur</surname>
<given-names>M</given-names>
</name>
, <name><surname>Madura</surname>
<given-names>J</given-names>
</name>
, <name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
(<year>2013</year>
) <article-title>Global transitions of proteins explored by a multiscale hybrid methodology: Application to adenylate kinase</article-title>
. <source>Biophys J</source>
<volume>105</volume>
: <fpage>1643</fpage>
–<lpage>1652</lpage>
.<pub-id pub-id-type="pmid">24094405</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-KantarciCarsibasi1"><label>65</label>
<mixed-citation publication-type="journal"><name><surname>Kantarci-Carsibasi</surname>
<given-names>N</given-names>
</name>
, <name><surname>Haliloglu</surname>
<given-names>T</given-names>
</name>
, <name><surname>Doruker</surname>
<given-names>P</given-names>
</name>
(<year>2008</year>
) <article-title>Conformational transition pathways explored by Monte Carlo simulation integrated with collective modes</article-title>
. <source>Biophys J</source>
<volume>95</volume>
: <fpage>5862</fpage>
–<lpage>5873</lpage>
.<pub-id pub-id-type="pmid">18676657</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Arora1"><label>66</label>
<mixed-citation publication-type="journal"><name><surname>Arora</surname>
<given-names>K</given-names>
</name>
, <name><surname>Brooks</surname>
<given-names>CL</given-names>
</name>
(<year>2007</year>
) <article-title>Large-scale allosteric conformational transitions of adenylate kinase appear to involve a population-shift mechanism</article-title>
. <source>Proc Natl Acad Sci U S A</source>
<volume>104</volume>
: <fpage>18496</fpage>
–<lpage>18501</lpage>
.<pub-id pub-id-type="pmid">18000050</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Daily1"><label>67</label>
<mixed-citation publication-type="journal"><name><surname>Daily</surname>
<given-names>MD</given-names>
</name>
, <name><surname>Jr</surname>
<given-names>GNP</given-names>
</name>
, <name><surname>Cui</surname>
<given-names>Q</given-names>
</name>
(<year>2010</year>
) <article-title>Many local motions cooperate to produce the adenylate kinase conformational transition</article-title>
. <source>J Mol Biol</source>
<volume>400</volume>
: <fpage>618</fpage>
–<lpage>631</lpage>
.<pub-id pub-id-type="pmid">20471396</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Daily2"><label>68</label>
<mixed-citation publication-type="journal"><name><surname>Daily</surname>
<given-names>MD</given-names>
</name>
, <name><surname>Phillips</surname>
<given-names>GN</given-names>
<suffix>Jr</suffix>
</name>
, <name><surname>Cui</surname>
<given-names>Q</given-names>
</name>
(<year>2011</year>
) <article-title>Interconversion of functional motions between mesophilic and thermophilic adenylate kinases</article-title>
. <source>PLoS Comput Biol</source>
<volume>7</volume>
: <fpage>e1002103</fpage>
.<pub-id pub-id-type="pmid">21779157</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Whitford1"><label>69</label>
<mixed-citation publication-type="journal"><name><surname>Whitford</surname>
<given-names>PC</given-names>
</name>
, <name><surname>Miyashita</surname>
<given-names>O</given-names>
</name>
, <name><surname>Levy</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Onuchic</surname>
<given-names>JN</given-names>
</name>
(<year>2007</year>
) <article-title>Conformational transitions of adenylate kinase: Switching by cracking</article-title>
. <source>J Mol Biol</source>
<volume>366</volume>
: <fpage>1661</fpage>
–<lpage>1671</lpage>
.<pub-id pub-id-type="pmid">17217965</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Whitford2"><label>70</label>
<mixed-citation publication-type="journal"><name><surname>Whitford</surname>
<given-names>PC</given-names>
</name>
, <name><surname>Gosavi</surname>
<given-names>S</given-names>
</name>
, <name><surname>Onuchic</surname>
<given-names>JN</given-names>
</name>
(<year>2008</year>
) <article-title>Conformational transitions in adenylate kinase: Allosteric communication reduces misligation</article-title>
. <source>J Biol Chem</source>
<volume>283</volume>
: <fpage>2042</fpage>
–<lpage>2048</lpage>
.<pub-id pub-id-type="pmid">17998210</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Brokaw1"><label>71</label>
<mixed-citation publication-type="journal"><name><surname>Brokaw</surname>
<given-names>JB</given-names>
</name>
, <name><surname>Chu</surname>
<given-names>JW</given-names>
</name>
(<year>2010</year>
) <article-title>On the roles of substrate binding and hinge unfolding in conformational changes of adenylate kinase</article-title>
. <source>Biophys J</source>
<volume>99</volume>
: <fpage>3420</fpage>
–<lpage>3429</lpage>
.<pub-id pub-id-type="pmid">21081091</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Beckstein1"><label>72</label>
<mixed-citation publication-type="journal"><name><surname>Beckstein</surname>
<given-names>O</given-names>
</name>
, <name><surname>Denning</surname>
<given-names>EJ</given-names>
</name>
, <name><surname>Perilla</surname>
<given-names>JR</given-names>
</name>
, <name><surname>Woolf</surname>
<given-names>TB</given-names>
</name>
(<year>2009</year>
) <article-title>Zipping and unzipping of adenylate kinase: Atomistic insights into the ensemble of open ↔ closed transitions</article-title>
. <source>J Mol Biol</source>
<volume>394</volume>
: <fpage>160</fpage>
–<lpage>176</lpage>
.<pub-id pub-id-type="pmid">19751742</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Yamashita1"><label>73</label>
<mixed-citation publication-type="journal"><name><surname>Yamashita</surname>
<given-names>A</given-names>
</name>
, <name><surname>Singh</surname>
<given-names>SK</given-names>
</name>
, <name><surname>Kawate</surname>
<given-names>T</given-names>
</name>
, <name><surname>Jin</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Gouaux</surname>
<given-names>E</given-names>
</name>
(<year>2005</year>
) <article-title>Crystal structure of a bacterial homologue of Na<sup>+</sup>
/Cl<sup>−</sup>
-dependent neurotransmitter transporters</article-title>
. <source>Nature</source>
<volume>437</volume>
: <fpage>215</fpage>
–<lpage>223</lpage>
.<pub-id pub-id-type="pmid">16041361</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Krishnamurthy1"><label>74</label>
<mixed-citation publication-type="journal"><name><surname>Krishnamurthy</surname>
<given-names>H</given-names>
</name>
, <name><surname>Gouaux</surname>
<given-names>E</given-names>
</name>
(<year>2012</year>
) <article-title>X-ray structures of LeuT in substrate-free outward-open and apo inward-open states</article-title>
. <source>Nature</source>
<volume>481</volume>
: <fpage>469</fpage>
–<lpage>474</lpage>
.<pub-id pub-id-type="pmid">22230955</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Cheng1"><label>75</label>
<mixed-citation publication-type="journal"><name><surname>Cheng</surname>
<given-names>MH</given-names>
</name>
, <name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
<article-title>Coupled global and local changes direct substrate translocation by neurotransmitter-sodium symporter ortholog</article-title>
. <source>Biophys J</source>
<volume>105</volume>
: <fpage>630</fpage>
–<lpage>639</lpage>
.<pub-id pub-id-type="pmid">23931311</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Yernool1"><label>76</label>
<mixed-citation publication-type="journal"><name><surname>Yernool</surname>
<given-names>D</given-names>
</name>
, <name><surname>Boudker</surname>
<given-names>O</given-names>
</name>
, <name><surname>Jin</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Gouaux</surname>
<given-names>E</given-names>
</name>
(<year>2004</year>
) <article-title>Structure of a glutamate transporter homologue from pyrococcus horikoshii</article-title>
. <source>Nature</source>
<volume>431</volume>
: <fpage>811</fpage>
–<lpage>818</lpage>
.<pub-id pub-id-type="pmid">15483603</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Reyes1"><label>77</label>
<mixed-citation publication-type="journal"><name><surname>Reyes</surname>
<given-names>N</given-names>
</name>
, <name><surname>Ginter</surname>
<given-names>C</given-names>
</name>
, <name><surname>Boudker</surname>
<given-names>O</given-names>
</name>
(<year>2009</year>
) <article-title>Transport mechanism of a bacterial homologue of glutamate transporters</article-title>
. <source>Nature</source>
<volume>462</volume>
: <fpage>880</fpage>
–<lpage>885</lpage>
.<pub-id pub-id-type="pmid">19924125</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Verdon1"><label>78</label>
<mixed-citation publication-type="journal"><name><surname>Verdon</surname>
<given-names>G</given-names>
</name>
, <name><surname>Boudker</surname>
<given-names>O</given-names>
</name>
(<year>2012</year>
) <article-title>Crystal structure of an asymmetric trimer of a bacterial glutamate transporter homolog</article-title>
. <source>Nat Struct Mol Biol</source>
<volume>19</volume>
: <fpage>355</fpage>
–<lpage>357</lpage>
.<pub-id pub-id-type="pmid">22343718</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Jiang1"><label>79</label>
<mixed-citation publication-type="journal"><name><surname>Jiang</surname>
<given-names>J</given-names>
</name>
, <name><surname>Shrivastava</surname>
<given-names>IH</given-names>
</name>
, <name><surname>Watts</surname>
<given-names>SD</given-names>
</name>
, <name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
, <name><surname>Amara</surname>
<given-names>SG</given-names>
</name>
(<year>2011</year>
) <article-title>Large collective motions regulate the functional properties of glutamate transporter trimers</article-title>
. <source>Proc Natl Acad Sci U S A</source>
<volume>108</volume>
: <fpage>15141</fpage>
–<lpage>15146</lpage>
.<pub-id pub-id-type="pmid">21876140</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Zomot1"><label>80</label>
<mixed-citation publication-type="journal"><name><surname>Zomot</surname>
<given-names>E</given-names>
</name>
, <name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
(<year>2013</year>
) <article-title>Intracellular gating in an inward-facing state of aspartate transporter Glt<sub>Ph</sub>
is regulated by the movements of the helical hairpin HP2</article-title>
. <source>J Biol Chem</source>
<volume>288</volume>
: <fpage>8231</fpage>
–<lpage>8237</lpage>
.<pub-id pub-id-type="pmid">23386619</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Shrivastava1"><label>81</label>
<mixed-citation publication-type="journal"><name><surname>Shrivastava</surname>
<given-names>IH</given-names>
</name>
, <name><surname>Jiang</surname>
<given-names>J</given-names>
</name>
, <name><surname>Amara</surname>
<given-names>SG</given-names>
</name>
, <name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
(<year>2008</year>
) <article-title>Time-resolved mechanism of extracellular gate opening and substrate binding in a glutamate transporter</article-title>
. <source>J Biol Chem</source>
<volume>283</volume>
: <fpage>28680</fpage>
–<lpage>28690</lpage>
.<pub-id pub-id-type="pmid">18678877</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Lezon1"><label>82</label>
<mixed-citation publication-type="journal"><name><surname>Lezon</surname>
<given-names>TR</given-names>
</name>
, <name><surname>Bahar</surname>
<given-names>I</given-names>
</name>
(<year>2012</year>
) <article-title>Constraints imposed by the membrane selectively guide the alternating access dynamics of the glutamate transporter Glt<sub>Ph</sub>
</article-title>
. <source>Biophys J</source>
<volume>102</volume>
: <fpage>1331</fpage>
–<lpage>1340</lpage>
.<pub-id pub-id-type="pmid">22455916</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Toyoshima1"><label>83</label>
<mixed-citation publication-type="journal"><name><surname>Toyoshima</surname>
<given-names>C</given-names>
</name>
(<year>2008</year>
) <article-title>Structural aspects of ion pumping by Ca<sup>2+</sup>
-ATPase of sarcoplasmic reticulum</article-title>
. <source>Arch Biochem Biophys</source>
<volume>476</volume>
: <fpage>3</fpage>
–<lpage>11</lpage>
.<pub-id pub-id-type="pmid">18455499</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Toyoshima2"><label>84</label>
<mixed-citation publication-type="journal"><name><surname>Toyoshima</surname>
<given-names>C</given-names>
</name>
, <name><surname>Nakasako</surname>
<given-names>M</given-names>
</name>
, <name><surname>Nomura</surname>
<given-names>H</given-names>
</name>
, <name><surname>Ogawa</surname>
<given-names>H</given-names>
</name>
(<year>2000</year>
) <article-title>Crystal structure of the calcium pump of sarcoplasmic reticulum at 2.6 Å resolution</article-title>
. <source>Nature</source>
<volume>405</volume>
: <fpage>647</fpage>
–<lpage>655</lpage>
.<pub-id pub-id-type="pmid">10864315</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Sorensen1"><label>85</label>
<mixed-citation publication-type="journal"><name><surname>Sorensen</surname>
<given-names>TLM</given-names>
</name>
, <name><surname>Moller</surname>
<given-names>JV</given-names>
</name>
, <name><surname>Nissen</surname>
<given-names>P</given-names>
</name>
(<year>2004</year>
) <article-title>Phosphoryl transfer and calcium ion occlusion in the calcium pump</article-title>
. <source>Science</source>
<volume>304</volume>
: <fpage>1672</fpage>
–<lpage>1675</lpage>
.<pub-id pub-id-type="pmid">15192230</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Toyoshima3"><label>86</label>
<mixed-citation publication-type="journal"><name><surname>Toyoshima</surname>
<given-names>C</given-names>
</name>
(<year>2009</year>
) <article-title>How Ca<sup>2+</sup>
-ATPase pumps ions across the sarcoplasmic reticulum membrane</article-title>
. <source>Biochim Biophys Acta</source>
<volume>1793</volume>
: <fpage>941</fpage>
–<lpage>946</lpage>
.<pub-id pub-id-type="pmid">19010358</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-EspinozaFonseca1"><label>87</label>
<mixed-citation publication-type="journal"><name><surname>Espinoza-Fonseca</surname>
<given-names>LM</given-names>
</name>
, <name><surname>Thomas</surname>
<given-names>DD</given-names>
</name>
(<year>2011</year>
) <article-title>Atomic-level characterization of the activation mechanism of SERCA by calcium</article-title>
. <source>PLoS One</source>
<volume>6</volume>
: <fpage>e26936</fpage>
.<pub-id pub-id-type="pmid">22046418</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Winters1"><label>88</label>
<mixed-citation publication-type="journal"><name><surname>Winters</surname>
<given-names>DL</given-names>
</name>
, <name><surname>Autry</surname>
<given-names>JM</given-names>
</name>
, <name><surname>Svensson</surname>
<given-names>B</given-names>
</name>
, <name><surname>Thomas</surname>
<given-names>DD</given-names>
</name>
(<year>2008</year>
) <article-title>Interdomain fluorescence resonance energy transfer in SERCA probed by cyan-fluorescent protein fused to the actuator domain</article-title>
. <source>Biochemistry</source>
<volume>47</volume>
: <fpage>4246</fpage>
–<lpage>4256</lpage>
.<pub-id pub-id-type="pmid">18338856</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Tekpinar1"><label>89</label>
<mixed-citation publication-type="journal"><name><surname>Tekpinar</surname>
<given-names>M</given-names>
</name>
, <name><surname>Zheng</surname>
<given-names>W</given-names>
</name>
(<year>2010</year>
) <article-title>Predicting order of conformational changes during protein conformational transitions using an interpolated elastic network model</article-title>
. <source>Proteins</source>
<volume>78</volume>
: <fpage>2469</fpage>
–<lpage>2481</lpage>
.<pub-id pub-id-type="pmid">20602461</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Krger1"><label>90</label>
<mixed-citation publication-type="journal"><name><surname>Krger</surname>
<given-names>DM</given-names>
</name>
, <name><surname>Ahmed</surname>
<given-names>A</given-names>
</name>
, <name><surname>Gohlke</surname>
<given-names>H</given-names>
</name>
(<year>2012</year>
) <article-title>NMSim Web Server: integrated approach for normal modebased geometric simulations of biologically relevant conformational transitions in proteins</article-title>
. <source>Nucleic Acids Res</source>
<volume>40</volume>
: <fpage>W310</fpage>
–<lpage>W316</lpage>
.<pub-id pub-id-type="pmid">22669906</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Sfriso1"><label>91</label>
<mixed-citation publication-type="journal"><name><surname>Sfriso</surname>
<given-names>P</given-names>
</name>
, <name><surname>Emperador</surname>
<given-names>A</given-names>
</name>
, <name><surname>Orellana</surname>
<given-names>L</given-names>
</name>
, <name><surname>Hospital</surname>
<given-names>A</given-names>
</name>
, <name><surname>Gelpí</surname>
<given-names>JL</given-names>
</name>
, <etal>et al</etal>
(<year>2012</year>
) <article-title>Finding conformational transition pathways from discrete molecular dynamics simulations</article-title>
. <source>J Chem Theory Comput</source>
<volume>8</volume>
: <fpage>4707</fpage>
–<lpage>4718</lpage>
.<pub-id pub-id-type="pmid">26605625</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-E1"><label>92</label>
<mixed-citation publication-type="journal"><name><surname>E</surname>
<given-names>W</given-names>
</name>
, <name><surname>Ren</surname>
<given-names>W</given-names>
</name>
, <name><surname>Vanden-Eijnden</surname>
<given-names>E</given-names>
</name>
(<year>2007</year>
) <article-title>Simplified and improved string method for computing the minimum energy paths in barrier-crossing events</article-title>
. <source>J Chem Phys</source>
<volume>126</volume>
: <fpage>164103</fpage>
.<pub-id pub-id-type="pmid">17477585</pub-id>
</mixed-citation>
</ref>
<ref id="pcbi.1003521-Adelman1"><label>93</label>
<mixed-citation publication-type="journal"><name><surname>Adelman</surname>
<given-names>JL</given-names>
</name>
, <name><surname>Dale</surname>
<given-names>AL</given-names>
</name>
, <name><surname>Zwier</surname>
<given-names>MC</given-names>
</name>
, <name><surname>Bhatt</surname>
<given-names>D</given-names>
</name>
, <name><surname>Chong</surname>
<given-names>LT</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Simulations of the alternating access mechanism of the sodium symporter mhp1</article-title>
. <source>Biophys J</source>
<volume>101</volume>
: <fpage>2399</fpage>
–<lpage>2407</lpage>
.<pub-id pub-id-type="pmid">22098738</pub-id>
</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Amérique/explor/PittsburghV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000255 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000255 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Wicri/Amérique |area= PittsburghV1 |flux= Pmc |étape= Corpus |type= RBID |clé= PMC:3974643 |texte= Exploring the Conformational Transitions of Biomolecular Systems Using a Simple Two-State Anisotropic Network Model }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:24699246" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a PittsburghV1
![]() | This area was generated with Dilib version V0.6.38. | ![]() |