1
|
Zschau RL, Zacharias M. Mechanism of β-hairpin formation in AzoChignolin and Chignolin. J Comput Chem 2023; 44:988-1001. [PMID: 36575994 DOI: 10.1002/jcc.27059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2022] [Revised: 11/23/2022] [Accepted: 11/30/2022] [Indexed: 12/29/2022]
Abstract
AzoChignolin is a photoswitchable variant of the mini-protein Chignolin with an azobenzene (AMPP) replacing the central loop. AzoChignolin is unfolded with AMPP in the trans-isomer. Transition to the cis-isomer causes β-hairpin folding similar to Chignolin. The AzoChignolin system is excellently suited for comprehensive analysis of folding nucleation kinetics. Utilizing multiple long-time MD simulations of AzoChignolin and Chignolin in MeOH and water, we estimated Markov models to examine folding kinetics of both peptides. We show that while AzoChignolin mimics Chignolin's structure well, the folding kinetics are quite different. Not only folding times but also intermediate states differ, particularly Chignolin is able to fold in MeOH into an α-helical intermediate which is impossible to form in AzoChignolin. The Markov models demonstrate that AzoChignolin's kinetics are generally faster, specifically when comparing the two main microfolding processes of hydrophobic collapse and turn formation. Photoswitchable loops are used frequently to understand the kinetics of elementary protein folding nucleation. However, our results indicate that intermediates and folding kinetics may differ between natural loops and photoswitchable variants.
Collapse
Affiliation(s)
- Richard L Zschau
- Physics Department and Center of Protein Assemblies, Technical University of Munich, Garching, Germany
| | - Martin Zacharias
- Physics Department and Center of Protein Assemblies, Technical University of Munich, Garching, Germany
| |
Collapse
|
2
|
Schmid F. Understanding and Modeling Polymers: The Challenge of Multiple Scales. ACS POLYMERS AU 2022. [DOI: 10.1021/acspolymersau.2c00049] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Friederike Schmid
- Institut für Physik, Johannes Gutenberg-Universität Mainz, Staudingerweg 9, 55128Mainz, Germany
| |
Collapse
|
3
|
Beyerle ER, Guenza MG. Identifying the leading dynamics of ubiquitin: A comparison between the tICA and the LE4PD slow fluctuations in amino acids' position. J Chem Phys 2021; 155:244108. [PMID: 34972386 DOI: 10.1063/5.0059688] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Molecular Dynamics (MD) simulations of proteins implicitly contain the information connecting the atomistic molecular structure and proteins' biologically relevant motion, where large-scale fluctuations are deemed to guide folding and function. In the complex multiscale processes described by MD trajectories, it is difficult to identify, separate, and study those large-scale fluctuations. This problem can be formulated as the need to identify a small number of collective variables that guide the slow kinetic processes. The most promising method among the ones used to study the slow leading processes in proteins' dynamics is the time-structure based on time-lagged independent component analysis (tICA), which identifies the dominant components in a noisy signal. Recently, we developed an anisotropic Langevin approach for the dynamics of proteins, called the anisotropic Langevin Equation for Protein Dynamics or LE4PD-XYZ. This approach partitions the protein's MD dynamics into mostly uncorrelated, wavelength-dependent, diffusive modes. It associates with each mode a free-energy map, where one measures the spatial extension and the time evolution of the mode-dependent, slow dynamical fluctuations. Here, we compare the tICA modes' predictions with the collective LE4PD-XYZ modes. We observe that the two methods consistently identify the nature and extension of the slowest fluctuation processes. The tICA separates the leading processes in a smaller number of slow modes than the LE4PD does. The LE4PD provides time-dependent information at short times and a formal connection to the physics of the kinetic processes that are missing in the pure statistical analysis of tICA.
Collapse
Affiliation(s)
- E R Beyerle
- Institute for Fundamental Science and Department of Chemistry and Biochemistry, University of Oregon, Eugene, Oregon 97403, USA
| | - M G Guenza
- Institute for Fundamental Science and Department of Chemistry and Biochemistry, University of Oregon, Eugene, Oregon 97403, USA
| |
Collapse
|
4
|
Mardt A, Noé F. Progress in deep Markov state modeling: Coarse graining and experimental data restraints. J Chem Phys 2021; 155:214106. [PMID: 34879670 DOI: 10.1063/5.0064668] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Recent advances in deep learning frameworks have established valuable tools for analyzing the long-timescale behavior of complex systems, such as proteins. In particular, the inclusion of physical constraints, e.g., time-reversibility, was a crucial step to make the methods applicable to biophysical systems. Furthermore, we advance the method by incorporating experimental observables into the model estimation showing that biases in simulation data can be compensated for. We further develop a new neural network layer in order to build a hierarchical model allowing for different levels of details to be studied. Finally, we propose an attention mechanism, which highlights important residues for the classification into different states. We demonstrate the new methodology on an ultralong molecular dynamics simulation of the Villin headpiece miniprotein.
Collapse
Affiliation(s)
- Andreas Mardt
- Department of Mathematics and Computer Science, Freie Universität Berlin, Berlin, Germany
| | - Frank Noé
- Department of Mathematics and Computer Science, Freie Universität Berlin, Berlin, Germany
| |
Collapse
|
5
|
Kraml J, Hofer F, Quoika PK, Kamenik AS, Liedl KR. X-Entropy: A Parallelized Kernel Density Estimator with Automated Bandwidth Selection to Calculate Entropy. J Chem Inf Model 2021; 61:1533-1538. [PMID: 33719418 PMCID: PMC8154256 DOI: 10.1021/acs.jcim.0c01375] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
X-Entropy is a Python package used to calculate the entropy of a given distribution, in this case, based on the distribution of dihedral angles. The dihedral entropy facilitates an alignment-independent measure of local protein flexibility. The key feature of our approach is a Gaussian kernel density estimation (KDE) using a plug-in bandwidth selection, which is fully implemented in a C++ backend and parallelized with OpenMP. We further provide a Python frontend, with predefined wrapper functions for classical coordinate-based dihedral entropy calculations, using a 1D approximation. This makes the package very straightforward to include in any Python-based analysis workflow. Furthermore, the frontend allows full access to the C++ backend, so that the KDE can be used on any binnable one-dimensional input data. In this application note, we discuss implementation and usage details and illustrate potential applications. In particular, we benchmark the performance of our module in calculating the entropy of samples drawn from a Gaussian distribution and the analytical solution thereof. Further, we analyze the computational performance of this module compared to well-established python libraries that perform KDE analyses. X-Entropy is available free of charge on GitHub (https://github.com/liedllab/X-Entropy).
Collapse
Affiliation(s)
- Johannes Kraml
- Institute
for General, Inorganic and Theoretical Chemistry, Center for Molecular
Biosciences Innsbruck (CMBI), University
of Innsbruck, A-6020 Innsbruck, Austria
| | - Florian Hofer
- Institute
for General, Inorganic and Theoretical Chemistry, Center for Molecular
Biosciences Innsbruck (CMBI), University
of Innsbruck, A-6020 Innsbruck, Austria
| | - Patrick K. Quoika
- Institute
for General, Inorganic and Theoretical Chemistry, Center for Molecular
Biosciences Innsbruck (CMBI), University
of Innsbruck, A-6020 Innsbruck, Austria
| | - Anna S. Kamenik
- Institute
for General, Inorganic and Theoretical Chemistry, Center for Molecular
Biosciences Innsbruck (CMBI), University
of Innsbruck, A-6020 Innsbruck, Austria
| | - Klaus R. Liedl
- Institute
for General, Inorganic and Theoretical Chemistry, Center for Molecular
Biosciences Innsbruck (CMBI), University
of Innsbruck, A-6020 Innsbruck, Austria,
| |
Collapse
|
6
|
Hofer F, Kraml J, Kahler U, Kamenik AS, Liedl KR. Catalytic Site p Ka Values of Aspartic, Cysteine, and Serine Proteases: Constant pH MD Simulations. J Chem Inf Model 2020; 60:3030-3042. [PMID: 32348143 PMCID: PMC7312390 DOI: 10.1021/acs.jcim.0c00190] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
![]()
Enzymatic function and activity of
proteases is closely controlled
by the pH value. The protonation states of titratable residues in
the active site react to changes in the pH value, according to their
pKa, and thereby determine the functionality
of the enzyme. Knowledge of the titration behavior of these residues
is crucial for the development of drugs targeting the active site
residues. However, experimental pKa data
are scarce, since the systems’ size and complexity make determination
of these pKa values inherently difficult.
In this study, we use single pH constant pH MD simulations as a fast
and robust tool to estimate the active site pKa values of a set of aspartic, cysteine, and serine proteases.
We capture characteristic pKa shifts of
the active site residues, which dictate the experimentally determined
activity profiles of the respective protease family. We find clear
differences of active site pKa values
within the respective families, which closely match the experimentally
determined pH preferences of the respective proteases. These shifts
are caused by a distinct network of electrostatic interactions characteristic
for each protease family. While we find convincing agreement with
experimental data for serine and aspartic proteases, we observe clear
deficiencies in the description of the titration behavior of cysteines
within the constant pH MD framework and highlight opportunities for
improvement. Consequently, with this work, we provide a concise set
of active site pKa values of aspartic
and serine proteases, which could serve as reference for future theoretical
as well as experimental studies.
Collapse
Affiliation(s)
- Florian Hofer
- Institute for General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Johannes Kraml
- Institute for General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Ursula Kahler
- Institute for General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Anna S Kamenik
- Institute for General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Klaus R Liedl
- Institute for General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| |
Collapse
|
7
|
Dreßler C, Kabbe G, Brehm M, Sebastiani D. Exploring non-equilibrium molecular dynamics of mobile protons in the solid acid CsH2PO4 at the micrometer and microsecond scale. J Chem Phys 2020; 152:164110. [DOI: 10.1063/5.0002167] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Christian Dreßler
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Gabriel Kabbe
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Martin Brehm
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Daniel Sebastiani
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| |
Collapse
|
8
|
Dreßler C, Kabbe G, Brehm M, Sebastiani D. Dynamical matrix propagator scheme for large-scale proton dynamics simulations. J Chem Phys 2020; 152:114114. [PMID: 32199428 DOI: 10.1063/1.5140635] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
We derive a matrix formalism for the simulation of long range proton dynamics for extended systems and timescales. On the basis of an ab initio molecular dynamics simulation, we construct a Markov chain, which allows us to store the entire proton dynamics in an M × M transition matrix (where M is the number of oxygen atoms). In this article, we start from common topology features of the hydrogen bond network of good proton conductors and utilize them as constituent constraints of our dynamic model. We present a thorough mathematical derivation of our approach and verify its uniqueness and correct asymptotic behavior. We propagate the proton distribution by means of transition matrices, which contain kinetic data from both ultra-short (sub-ps) and intermediate (ps) timescales. This concept allows us to keep the most relevant features from the microscopic level while effectively reaching larger time and length scales. We demonstrate the applicability of the transition matrices for the description of proton conduction trends in proton exchange membrane materials.
Collapse
Affiliation(s)
- Christian Dreßler
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Gabriel Kabbe
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Martin Brehm
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Daniel Sebastiani
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| |
Collapse
|
9
|
Hofer F, Dietrich V, Kamenik AS, Tollinger M, Liedl KR. pH-Dependent Protonation of the Phl p 6 Pollen Allergen Studied by NMR and cpH-aMD. J Chem Theory Comput 2019; 15:5716-5726. [PMID: 31476118 PMCID: PMC6994067 DOI: 10.1021/acs.jctc.9b00540] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We use state-of-the-art NMR experiments to measure apparent pKa values in the native protein environment and employ a cutting-edge combination of enhanced sampling and constant pH molecular dynamics (MD) simulations to rationalize strong pKa shifts. The major timothy grass pollen allergen Phl p 6 serves as an ideal model system for both methods due to its high number of titratable residues despite its comparably small size. We present a proton transition analysis as intuitive tool to depict the captured protonation state ensemble in atomistic detail. Combining microscopic structural details from MD simulations and macroscopic ensemble averages from NMR shifts leads to a comprehensive view on pH dependencies of protonation states and tautomers. Overall, we find striking agreement between simulation-based pKa predictions and experiment. However, our analyses suggest subtle differences in the underlying molecular origin of the observed pKa shifts. From accelerated constant pH MD simulations, we identify immediate proximity of opposite charges, followed by vicinity of equal charges as major driving forces for pKa shifts. NMR experiments on the other hand, suggest only a weak relation of pKa shifts and close contacts to charged residues, while the strongest influence derives from the dipolar character of α helices. The presented study hence pinpoints opportunities for improvements concerning the theoretical description of protonation state and tautomer probabilities. However, the coherence in the resulting apparent pKa values from simulations and experiment affirms cpH-aMD as a reliable tool to study allergen dynamics at varying pH levels.
Collapse
Affiliation(s)
- Florian Hofer
- †Institute
for General, Inorganic and Theoretical Chemistry and ‡Institute for Organic Chemistry,
Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Valentin Dietrich
- †Institute
for General, Inorganic and Theoretical Chemistry and ‡Institute for Organic Chemistry,
Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Anna S. Kamenik
- †Institute
for General, Inorganic and Theoretical Chemistry and ‡Institute for Organic Chemistry,
Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Martin Tollinger
- †Institute
for General, Inorganic and Theoretical Chemistry and ‡Institute for Organic Chemistry,
Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Klaus R. Liedl
- †Institute
for General, Inorganic and Theoretical Chemistry and ‡Institute for Organic Chemistry,
Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria,E-mail:
| |
Collapse
|
10
|
Remington JM, McCullagh M, Kohler B. Molecular Dynamics Simulations of 2-Aminopurine-Labeled Dinucleoside Monophosphates Reveal Multiscale Stacking Kinetics. J Phys Chem B 2019; 123:2291-2304. [PMID: 30767498 DOI: 10.1021/acs.jpcb.8b12172] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics (MD) simulations of 2-aminopurine (2Ap)-labeled DNA dinucleoside monophosphates (DNMPs) were performed to investigate the hypothesis that base stacking dynamics occur on timescales sufficiently rapid to influence the emission signals measured in time-resolved fluorescence experiments. Analysis of multiple microsecond-length trajectories shows that the DNMPs sample all four coplanar stacking motifs. In addition, three metastable unstacked conformations are detected. A hidden Markov-state model (HMSM) was applied to the simulations to estimate transition rates between the stacked and unstacked states. Transitions between different stacked states generally occur at higher rates when the number of nucleobase faces requiring desolvation is minimized. Time constants for structural relaxation range between 1.6 and 25 ns, suggesting that emission from photoexcited 2Ap, which has an excited-state lifetime of 10 ns, is sensitive to base stacking kinetics. A master equation model for the excited-state population of 2Ap predicts multiexponential emission decays that reproduce the sub-10 ns emission decay lifetimes and amplitudes seen in experiments. Combining MD simulations with HMSM analysis is a powerful way to understand the dynamics that influence 2Ap excited-state relaxation and represents an important step toward using observed emission signals to validate MD simulations.
Collapse
Affiliation(s)
- Jacob M Remington
- Department of Chemistry and Biochemistry , Montana State University , Bozeman , Montana 59717 , United States
| | - Martin McCullagh
- Department of Chemistry , Colorado State University , Fort Collins , Colorado 80523 , United States
| | - Bern Kohler
- Department of Chemistry and Biochemistry , The Ohio State University , 100 West 18th Avenue , Columbus , Ohio 43210 , United States
| |
Collapse
|
11
|
Abstract
Dynamic neutron scattering directly probes motions in biological systems on femtosecond to microsecond timescales. When combined with molecular dynamics simulation and normal mode analysis, detailed descriptions of the forms and frequencies of motions can be derived. We examine vibrations in proteins, the temperature dependence of protein motions, and concepts describing the rich variety of motions detectable using neutrons in biological systems at physiological temperatures. New techniques for deriving information on collective motions using coherent scattering are also reviewed.
Collapse
Affiliation(s)
- Jeremy C Smith
- UT/ORNL Center for Molecular Biophysics, Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6309, USA; .,Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Pan Tan
- School of Physics and Astronomy and Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Loukas Petridis
- UT/ORNL Center for Molecular Biophysics, Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6309, USA; .,Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Liang Hong
- School of Physics and Astronomy and Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, China
| |
Collapse
|
12
|
Affiliation(s)
- Brooke E. Husic
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Vijay S. Pande
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
13
|
McKiernan KA, Husic BE, Pande VS. Modeling the mechanism of CLN025 beta-hairpin formation. J Chem Phys 2017; 147:104107. [PMID: 28915754 PMCID: PMC5597441 DOI: 10.1063/1.4993207] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2017] [Accepted: 08/24/2017] [Indexed: 01/26/2023] Open
Abstract
Beta-hairpins are substructures found in proteins that can lend insight into more complex systems. Furthermore, the folding of beta-hairpins is a valuable test case for benchmarking experimental and theoretical methods. Here, we simulate the folding of CLN025, a miniprotein with a beta-hairpin structure, at its experimental melting temperature using a range of state-of-the-art protein force fields. We construct Markov state models in order to examine the thermodynamics, kinetics, mechanism, and rate-determining step of folding. Mechanistically, we find the folding process is rate-limited by the formation of the turn region hydrogen bonds, which occurs following the downhill hydrophobic collapse of the extended denatured protein. These results are presented in the context of established and contradictory theories of the beta-hairpin folding process. Furthermore, our analysis suggests that the AMBER-FB15 force field, at this temperature, best describes the characteristics of the full experimental CLN025 conformational ensemble, while the AMBER ff99SB-ILDN and CHARMM22* force fields display a tendency to overstabilize the native state.
Collapse
Affiliation(s)
- Keri A McKiernan
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Brooke E Husic
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Vijay S Pande
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
14
|
Wu H, Nüske F, Paul F, Klus S, Koltai P, Noé F. Variational Koopman models: Slow collective variables and molecular kinetics from short off-equilibrium simulations. J Chem Phys 2017; 146:154104. [DOI: 10.1063/1.4979344] [Citation(s) in RCA: 69] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Affiliation(s)
- Hao Wu
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Feliks Nüske
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Fabian Paul
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Stefan Klus
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Péter Koltai
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Frank Noé
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| |
Collapse
|
15
|
Olsson S, Noé F. Mechanistic Models of Chemical Exchange Induced Relaxation in Protein NMR. J Am Chem Soc 2016; 139:200-210. [DOI: 10.1021/jacs.6b09460] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Simon Olsson
- Computational Molecular
Biology,
FB Mathematik und Informatik, Freie Universität Berlin, Berlin 14195, Germany
| | - Frank Noé
- Computational Molecular
Biology,
FB Mathematik und Informatik, Freie Universität Berlin, Berlin 14195, Germany
| |
Collapse
|
16
|
Pérez-Hernández G, Noé F. Hierarchical Time-Lagged Independent Component Analysis: Computing Slow Modes and Reaction Coordinates for Large Molecular Systems. J Chem Theory Comput 2016; 12:6118-6129. [PMID: 27792332 DOI: 10.1021/acs.jctc.6b00738] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Analysis of molecular dynamics, for example using Markov models, often requires the identification of order parameters that are good indicators of the rare events, i.e. good reaction coordinates. Recently, it has been shown that the time-lagged independent component analysis (TICA) finds the linear combinations of input coordinates that optimally represent the slow kinetic modes and may serve in order to define reaction coordinates between the metastable states of the molecular system. A limitation of the method is that both computing time and memory requirements scale with the square of the number of input features. For large protein systems, this exacerbates the use of extensive feature sets such as the distances between all pairs of residues or even heavy atoms. Here we derive a hierarchical TICA (hTICA) method that approximates the full TICA solution by a hierarchical, divide-and-conquer calculation. By using hTICA on distances between heavy atoms we identify previously unknown relaxation processes in the bovine pancreatic trypsin inhibitor.
Collapse
Affiliation(s)
- Guillermo Pérez-Hernández
- Department of Mathematics and Computer Science, Freie Universitat Berlin , Arnimallee 6, Berlin, Germany 14195
| | - Frank Noé
- Department of Mathematics and Computer Science, Freie Universitat Berlin , Arnimallee 6, Berlin, Germany 14195
| |
Collapse
|
17
|
Lemke O, Keller BG. Density-based cluster algorithms for the identification of core sets. J Chem Phys 2016; 145:164104. [DOI: 10.1063/1.4965440] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Affiliation(s)
- Oliver Lemke
- Department of Biology, Chemistry, Pharmacy, Freie Universität Berlin, Takustraße 3, D-14195 Berlin, Germany
| | - Bettina G. Keller
- Department of Biology, Chemistry, Pharmacy, Freie Universität Berlin, Takustraße 3, D-14195 Berlin, Germany
| |
Collapse
|
18
|
Loring RF. Lattice model of spatial correlations in catalysis. J Chem Phys 2016; 145:134508. [DOI: 10.1063/1.4964282] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Roger F. Loring
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
19
|
Witek J, Keller BG, Blatter M, Meissner A, Wagner T, Riniker S. Kinetic Models of Cyclosporin A in Polar and Apolar Environments Reveal Multiple Congruent Conformational States. J Chem Inf Model 2016; 56:1547-62. [DOI: 10.1021/acs.jcim.6b00251] [Citation(s) in RCA: 73] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jagna Witek
- Laboratory
of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Bettina G. Keller
- Department
of Biology, Chemistry, Pharmacy, Freie Universität Berlin, Takustrasse 3, 14195 Berlin, Germany
| | - Markus Blatter
- Novartis
Institutes for BioMedical Research, Novartis Pharma AG, Novartis Campus, 4002 Basel, Switzerland
| | - Axel Meissner
- Novartis
Institutes for BioMedical Research, Novartis Pharma AG, Novartis Campus, 4002 Basel, Switzerland
| | - Trixie Wagner
- Novartis
Institutes for BioMedical Research, Novartis Pharma AG, Novartis Campus, 4002 Basel, Switzerland
| | - Sereina Riniker
- Laboratory
of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
20
|
Vural D, Hu X, Lindner B, Jain N, Miao Y, Cheng X, Liu Z, Hong L, Smith JC. Quasielastic neutron scattering in biology: Theory and applications. Biochim Biophys Acta Gen Subj 2016; 1861:3638-3650. [PMID: 27316321 DOI: 10.1016/j.bbagen.2016.06.015] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2016] [Revised: 06/08/2016] [Accepted: 06/09/2016] [Indexed: 02/03/2023]
Abstract
Neutrons scatter quasielastically from stochastic, diffusive processes, such as overdamped vibrations, localized diffusion and transitions between energy minima. In biological systems, such as proteins and membranes, these relaxation processes are of considerable physical interest. We review here recent methodological advances and applications of quasielastic neutron scattering (QENS) in biology, concentrating on the role of molecular dynamics simulation in generating data with which neutron profiles can be unambiguously interpreted. We examine the use of massively-parallel computers in calculating scattering functions, and the application of Markov state modeling. The decomposition of MD-derived neutron dynamic susceptibilities is described, and the use of this in combination with NMR spectroscopy. We discuss dynamics at very long times, including approximations to the infinite time mean-square displacement and nonequilibrium aspects of single-protein dynamics. Finally, we examine how neutron scattering and MD can be combined to provide information on lipid nanodomains. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.
Collapse
Affiliation(s)
- Derya Vural
- Center for Molecular Biophysics, Oak Ridge National Laboratory, TN 37831, USA; Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA
| | - Xiaohu Hu
- Center for Molecular Biophysics, Oak Ridge National Laboratory, TN 37831, USA; Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA
| | - Benjamin Lindner
- Institute of Natural Sciences & Department of Physics and Astronomy, Shanghai Jiao Tong University, 200240, China
| | - Nitin Jain
- Center for Molecular Biophysics, Oak Ridge National Laboratory, TN 37831, USA; Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA
| | - Yinglong Miao
- Center for Molecular Biophysics, Oak Ridge National Laboratory, TN 37831, USA; Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA
| | - Xiaolin Cheng
- Center for Molecular Biophysics, Oak Ridge National Laboratory, TN 37831, USA; Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA
| | - Zhuo Liu
- Institute of Natural Sciences & Department of Physics and Astronomy, Shanghai Jiao Tong University, 200240, China
| | - Liang Hong
- Institute of Natural Sciences & Department of Physics and Astronomy, Shanghai Jiao Tong University, 200240, China
| | - Jeremy C Smith
- Center for Molecular Biophysics, Oak Ridge National Laboratory, TN 37831, USA; Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA.
| |
Collapse
|
21
|
Reppert M, Tokmakoff A. Computational Amide I 2D IR Spectroscopy as a Probe of Protein Structure and Dynamics. Annu Rev Phys Chem 2016; 67:359-86. [DOI: 10.1146/annurev-physchem-040215-112055] [Citation(s) in RCA: 75] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Mike Reppert
- Department of Chemistry, James Franck Institute, Institute for Biophysical Dynamics, University of Chicago, Chicago, Illinois 60637;
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute, Institute for Biophysical Dynamics, University of Chicago, Chicago, Illinois 60637;
| |
Collapse
|
22
|
Rudzinski JF, Kremer K, Bereau T. Communication: Consistent interpretation of molecular simulation kinetics using Markov state models biased with external information. J Chem Phys 2016; 144:051102. [DOI: 10.1063/1.4941455] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, 55128 Mainz, Germany
| | - Tristan Bereau
- Max Planck Institute for Polymer Research, 55128 Mainz, Germany
| |
Collapse
|
23
|
Boninsegna L, Gobbo G, Noé F, Clementi C. Investigating Molecular Kinetics by Variationally Optimized Diffusion Maps. J Chem Theory Comput 2015; 11:5947-60. [DOI: 10.1021/acs.jctc.5b00749] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Lorenzo Boninsegna
- Center
for Theoretical Biological Physics and Department of Chemistry, Rice University, 6100 Main Street, Houston, Texas 77005, United States
| | - Gianpaolo Gobbo
- Maxwell
Institute for Mathematical Sciences and School of Mathematics, The University of Edinburgh, Peter Guthrie Tait Road, Edinburgh EH9 3FD, United Kingdom
| | - Frank Noé
- Department
of Mathematics, Computer Science and Bioinformatics, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Cecilia Clementi
- Center
for Theoretical Biological Physics and Department of Chemistry, Rice University, 6100 Main Street, Houston, Texas 77005, United States
| |
Collapse
|
24
|
Trendelkamp-Schroer B, Wu H, Paul F, Noé F. Estimation and uncertainty of reversible Markov models. J Chem Phys 2015; 143:174101. [DOI: 10.1063/1.4934536] [Citation(s) in RCA: 79] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
| | - Hao Wu
- Institut für Mathematik und Informatik, FU Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Fabian Paul
- Institut für Mathematik und Informatik, FU Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Frank Noé
- Institut für Mathematik und Informatik, FU Berlin, Arnimallee 6, 14195 Berlin, Germany
| |
Collapse
|
25
|
Scherer MK, Trendelkamp-Schroer B, Paul F, Pérez-Hernández G, Hoffmann M, Plattner N, Wehmeyer C, Prinz JH, Noé F. PyEMMA 2: A Software Package for Estimation, Validation, and Analysis of Markov Models. J Chem Theory Comput 2015; 11:5525-42. [PMID: 26574340 DOI: 10.1021/acs.jctc.5b00743] [Citation(s) in RCA: 700] [Impact Index Per Article: 77.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Markov (state) models (MSMs) and related models of molecular kinetics have recently received a surge of interest as they can systematically reconcile simulation data from either a few long or many short simulations and allow us to analyze the essential metastable structures, thermodynamics, and kinetics of the molecular system under investigation. However, the estimation, validation, and analysis of such models is far from trivial and involves sophisticated and often numerically sensitive methods. In this work we present the open-source Python package PyEMMA ( http://pyemma.org ) that provides accurate and efficient algorithms for kinetic model construction. PyEMMA can read all common molecular dynamics data formats, helps in the selection of input features, provides easy access to dimension reduction algorithms such as principal component analysis (PCA) and time-lagged independent component analysis (TICA) and clustering algorithms such as k-means, and contains estimators for MSMs, hidden Markov models, and several other models. Systematic model validation and error calculation methods are provided. PyEMMA offers a wealth of analysis functions such that the user can conveniently compute molecular observables of interest. We have derived a systematic and accurate way to coarse-grain MSMs to few states and to illustrate the structures of the metastable states of the system. Plotting functions to produce a manuscript-ready presentation of the results are available. In this work, we demonstrate the features of the software and show new methodological concepts and results produced by PyEMMA.
Collapse
Affiliation(s)
- Martin K Scherer
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | | | - Fabian Paul
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Guillermo Pérez-Hernández
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Moritz Hoffmann
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Nuria Plattner
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Christoph Wehmeyer
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Jan-Hendrik Prinz
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Frank Noé
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| |
Collapse
|
26
|
Wu Y, Stauffer SR, Stanfield RL, Tapia PH, Ursu O, Fisher GW, Szent-Gyorgyi C, Evangelisti A, Waller A, Strouse JJ, Carter MB, Bologa C, Gouveia K, Poslusney M, Waggoner AS, Lindsley CW, Jarvik JW, Sklar LA. Discovery of Small-Molecule Nonfluorescent Inhibitors of Fluorogen-Fluorogen Activating Protein Binding Pair. ACTA ACUST UNITED AC 2015; 21:74-87. [PMID: 26442911 DOI: 10.1177/1087057115609145] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2015] [Accepted: 09/09/2015] [Indexed: 11/17/2022]
Abstract
A new class of biosensors, fluorogen activating proteins (FAPs), has been successfully used to track receptor trafficking in live cells. Unlike the traditional fluorescent proteins (FPs), FAPs do not fluoresce unless bound to their specific small-molecule fluorogens, and thus FAP-based assays are highly sensitive. Application of the FAP-based assay for protein trafficking in high-throughput flow cytometry resulted in the discovery of a new class of compounds that interferes with the binding between fluorogens and FAP, thus blocking the fluorescence signal. These compounds are high-affinity, nonfluorescent analogs of fluorogens with little or no toxicity to the tested cells and no apparent interference with the normal function of FAP-tagged receptors. The most potent compound among these, N,4-dimethyl-N-(2-oxo-2-(4-(pyridin-2-yl)piperazin-1-yl)ethyl)benzenesulfonamide (ML342), has been investigated in detail. X-ray crystallographic analysis revealed that ML342 competes with the fluorogen, sulfonated thiazole orange coupled to diethylene glycol diamine (TO1-2p), for the same binding site on a FAP, AM2.2. Kinetic analysis shows that the FAP-fluorogen interaction is more complex than a homogeneous one-site binding process, with multiple conformational states of the fluorogen and/or the FAP, and possible dimerization of the FAP moiety involved in the process.
Collapse
Affiliation(s)
- Yang Wu
- Department of Pathology, University of New Mexico, Albuquerque, NM, USA Center for Molecular Discovery, University of New Mexico, Albuquerque, NM, USA
| | - Shaun R Stauffer
- Vanderbilt Specialized Chemistry Center, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Robyn L Stanfield
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, CA, USA
| | - Phillip H Tapia
- Center for Molecular Discovery, University of New Mexico, Albuquerque, NM, USA
| | - Oleg Ursu
- Center for Molecular Discovery, University of New Mexico, Albuquerque, NM, USA
| | - Gregory W Fisher
- Molecular Biosensor and Imaging Center, Carnegie Mellon University, Pittsburgh, PA, USA
| | | | - Annette Evangelisti
- Center for Molecular Discovery, University of New Mexico, Albuquerque, NM, USA
| | - Anna Waller
- Center for Molecular Discovery, University of New Mexico, Albuquerque, NM, USA
| | - J Jacob Strouse
- Center for Molecular Discovery, University of New Mexico, Albuquerque, NM, USA
| | - Mark B Carter
- Center for Molecular Discovery, University of New Mexico, Albuquerque, NM, USA
| | - Cristian Bologa
- Center for Molecular Discovery, University of New Mexico, Albuquerque, NM, USA
| | - Kristine Gouveia
- Center for Molecular Discovery, University of New Mexico, Albuquerque, NM, USA
| | - Mike Poslusney
- Vanderbilt Specialized Chemistry Center, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Alan S Waggoner
- Molecular Biosensor and Imaging Center, Carnegie Mellon University, Pittsburgh, PA, USA Department of Biological Science, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Craig W Lindsley
- Vanderbilt Specialized Chemistry Center, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Jonathan W Jarvik
- Molecular Biosensor and Imaging Center, Carnegie Mellon University, Pittsburgh, PA, USA Department of Biological Science, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Larry A Sklar
- Department of Pathology, University of New Mexico, Albuquerque, NM, USA Center for Molecular Discovery, University of New Mexico, Albuquerque, NM, USA
| |
Collapse
|
27
|
Vitalini F, Noé F, Keller BG. A Basis Set for Peptides for the Variational Approach to Conformational Kinetics. J Chem Theory Comput 2015; 11:3992-4004. [DOI: 10.1021/acs.jctc.5b00498] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- F. Vitalini
- Department
of Biology, Chemistry, Pharmacy, Freie Universität Berlin, Takustraße
3, D-14195 Berlin, Germany
| | - F. Noé
- Department
of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, D-14195 Berlin, Germany
| | - B. G. Keller
- Department
of Biology, Chemistry, Pharmacy, Freie Universität Berlin, Takustraße
3, D-14195 Berlin, Germany
| |
Collapse
|
28
|
Rakers C, Bermudez M, Keller BG, Mortier J, Wolber G. Computational close up on protein-protein interactions: how to unravel the invisible using molecular dynamics simulations? WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2015. [DOI: 10.1002/wcms.1222] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Christin Rakers
- Institute of Pharmacy; Freie Universität Berlin; Berlin Germany
| | - Marcel Bermudez
- Institute of Pharmacy; Freie Universität Berlin; Berlin Germany
| | - Bettina G. Keller
- Institute for Chemistry and Biochemistry; Freie Universität Berlin; Berlin Germany
| | - Jérémie Mortier
- Institute of Pharmacy; Freie Universität Berlin; Berlin Germany
| | - Gerhard Wolber
- Institute of Pharmacy; Freie Universität Berlin; Berlin Germany
| |
Collapse
|
29
|
Razavi AM, Voelz VA. Kinetic Network Models of Tryptophan Mutations in β-Hairpins Reveal the Importance of Non-Native Interactions. J Chem Theory Comput 2015; 11:2801-12. [DOI: 10.1021/acs.jctc.5b00088] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Asghar M. Razavi
- Department
of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Vincent A. Voelz
- Department
of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
30
|
Wang H, Schütte C. Building Markov State Models for Periodically Driven Non-Equilibrium Systems. J Chem Theory Comput 2015; 11:1819-31. [DOI: 10.1021/ct500997y] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Affiliation(s)
- Han Wang
- CAEP Software Center for High Performance Numerical Simulation, Beijing, China
- Zuse Institute Berlin (ZIB), Berlin, Germany
| | - Christof Schütte
- Zuse Institute Berlin (ZIB), Berlin, Germany
- Institute
for Mathematics, Freie Universität Berlin, Berlin, Germany
| |
Collapse
|
31
|
Vitalini F, Mey ASJS, Noé F, Keller BG. Dynamic properties of force fields. J Chem Phys 2015; 142:084101. [DOI: 10.1063/1.4909549] [Citation(s) in RCA: 60] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Affiliation(s)
- F. Vitalini
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, D-14195 Berlin, Germany
- Department of Biology, Chemistry, Pharmacy, Freie Universität Berlin, Takustraße 3, D-14195 Berlin, Germany
| | - A. S. J. S. Mey
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, D-14195 Berlin, Germany
| | - F. Noé
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, D-14195 Berlin, Germany
| | - B. G. Keller
- Department of Biology, Chemistry, Pharmacy, Freie Universität Berlin, Takustraße 3, D-14195 Berlin, Germany
| |
Collapse
|
32
|
Noé F, Wu H, Prinz JH, Plattner N. Projected and hidden Markov models for calculating kinetics and metastable states of complex molecules. J Chem Phys 2014; 139:184114. [PMID: 24320261 DOI: 10.1063/1.4828816] [Citation(s) in RCA: 115] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Markov state models (MSMs) have been successful in computing metastable states, slow relaxation timescales and associated structural changes, and stationary or kinetic experimental observables of complex molecules from large amounts of molecular dynamics simulation data. However, MSMs approximate the true dynamics by assuming a Markov chain on a clusters discretization of the state space. This approximation is difficult to make for high-dimensional biomolecular systems, and the quality and reproducibility of MSMs has, therefore, been limited. Here, we discard the assumption that dynamics are Markovian on the discrete clusters. Instead, we only assume that the full phase-space molecular dynamics is Markovian, and a projection of this full dynamics is observed on the discrete states, leading to the concept of Projected Markov Models (PMMs). Robust estimation methods for PMMs are not yet available, but we derive a practically feasible approximation via Hidden Markov Models (HMMs). It is shown how various molecular observables of interest that are often computed from MSMs can be computed from HMMs/PMMs. The new framework is applicable to both, simulation and single-molecule experimental data. We demonstrate its versatility by applications to educative model systems, a 1 ms Anton MD simulation of the bovine pancreatic trypsin inhibitor protein, and an optical tweezer force probe trajectory of an RNA hairpin.
Collapse
Affiliation(s)
- Frank Noé
- Department of Mathematics and Computer Science, FU Berlin, Arnimallee 6, 14159 Berlin, Germany
| | | | | | | |
Collapse
|
33
|
Chodera JD, Noé F. Markov state models of biomolecular conformational dynamics. Curr Opin Struct Biol 2014; 25:135-44. [PMID: 24836551 DOI: 10.1016/j.sbi.2014.04.002] [Citation(s) in RCA: 498] [Impact Index Per Article: 49.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2014] [Revised: 04/08/2014] [Accepted: 04/12/2014] [Indexed: 10/25/2022]
Abstract
It has recently become practical to construct Markov state models (MSMs) that reproduce the long-time statistical conformational dynamics of biomolecules using data from molecular dynamics simulations. MSMs can predict both stationary and kinetic quantities on long timescales (e.g. milliseconds) using a set of atomistic molecular dynamics simulations that are individually much shorter, thus addressing the well-known sampling problem in molecular dynamics simulation. In addition to providing predictive quantitative models, MSMs greatly facilitate both the extraction of insight into biomolecular mechanism (such as folding and functional dynamics) and quantitative comparison with single-molecule and ensemble kinetics experiments. A variety of methodological advances and software packages now bring the construction of these models closer to routine practice. Here, we review recent progress in this field, considering theoretical and methodological advances, new software tools, and recent applications of these approaches in several domains of biochemistry and biophysics, commenting on remaining challenges.
Collapse
Affiliation(s)
- John D Chodera
- Computational Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065, USA.
| | - Frank Noé
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany.
| |
Collapse
|
34
|
Nüske F, Keller BG, Pérez-Hernández G, Mey ASJS, Noé F. Variational Approach to Molecular Kinetics. J Chem Theory Comput 2014; 10:1739-52. [DOI: 10.1021/ct4009156] [Citation(s) in RCA: 210] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Feliks Nüske
- Department for Mathematics
and Computer Science, Freie Universität Berlin, 14195 Berlin, Germany
| | - Bettina G. Keller
- Department for Mathematics
and Computer Science, Freie Universität Berlin, 14195 Berlin, Germany
| | | | - Antonia S. J. S. Mey
- Department for Mathematics
and Computer Science, Freie Universität Berlin, 14195 Berlin, Germany
| | - Frank Noé
- Department for Mathematics
and Computer Science, Freie Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
35
|
Noé F, Prinz JH. Analysis of Markov Models. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2014; 797:75-90. [DOI: 10.1007/978-94-007-7606-7_6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/16/2023]
|
36
|
Yi Z, Lindner B, Prinz JH, Noé F, Smith JC. Dynamic neutron scattering from conformational dynamics. II. Application using molecular dynamics simulation and Markov modeling. J Chem Phys 2013; 139:175102. [DOI: 10.1063/1.4824071] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
|
37
|
Lindner B, Yi Z, Prinz JH, Smith JC, Noé F. Dynamic neutron scattering from conformational dynamics. I. Theory and Markov models. J Chem Phys 2013; 139:175101. [DOI: 10.1063/1.4824070] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
|
38
|
Pérez-Hernández G, Paul F, Giorgino T, De Fabritiis G, Noé F. Identification of slow molecular order parameters for Markov model construction. J Chem Phys 2013; 139:015102. [DOI: 10.1063/1.4811489] [Citation(s) in RCA: 605] [Impact Index Per Article: 55.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
|
39
|
Trendelkamp-Schroer B, Noé F. Efficient Bayesian estimation of Markov model transition matrices with given stationary distribution. J Chem Phys 2013; 138:164113. [DOI: 10.1063/1.4801325] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
|
40
|
Prigozhin MB, Gruebele M. Microsecond folding experiments and simulations: a match is made. Phys Chem Chem Phys 2013; 15:3372-88. [PMID: 23361200 PMCID: PMC3632410 DOI: 10.1039/c3cp43992e] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
For the past two decades, protein folding experiments have been speeding up from the second or millisecond time scale to the microsecond time scale, and full-atom simulations have been extended from the nanosecond to the microsecond and even millisecond time scale. Where the two meet, it is now possible to compare results directly, allowing force fields to be validated and refined, and allowing experimental data to be interpreted in atomistic detail. In this perspective we compare recent experiments and simulations on the microsecond time scale, pointing out the progress that has been made in determining native structures from physics-based simulations, refining experiments and simulations to provide more quantitative underlying mechanisms, and tackling the problems of multiple reaction coordinates, downhill folding, and complex underlying structure of unfolded or misfolded states.
Collapse
Affiliation(s)
- M. B. Prigozhin
- Department of Chemistry, Center for Biophsyics and Computational Biology, 600 South Mathews Ave. Box 5–6, Urbana IL 61801, USA
| | - M. Gruebele
- Department of Chemistry, Center for Biophsyics and Computational Biology, 600 South Mathews Ave. Box 5–6, Urbana IL 61801, USA
- Department of Physics, Center for Biophsyics and Computational Biology, 600 South Mathews Ave. Box 5–6, Urbana IL 61801, USA
| |
Collapse
|
41
|
Steger K, Bollmann S, Noé F, Doose S. Systematic evaluation of fluorescence correlation spectroscopy data analysis on the nanosecond time scale. Phys Chem Chem Phys 2013; 15:10435-45. [DOI: 10.1039/c3cp50644d] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
|
42
|
Held M, Imhof P, Keller BG, Noé F. Modulation of a ligand's energy landscape and kinetics by the chemical environment. J Phys Chem B 2012; 116:13597-607. [PMID: 23025812 DOI: 10.1021/jp3006684] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Understanding how the chemical environment modulates the predominant conformations and kinetics of flexible molecules is a core interest of biochemistry and a prerequisite for the rational design of synthetic catalysts. This study combines molecular dynamics simulation and Markov state models (MSMs) to a systematic computational strategy for investigating the effect of the chemical environment of a molecule on its conformations and kinetics. MSMs allow quantities to be computed that are otherwise difficult to access, such as the metastable sets, their free energies, and the relaxation time scales related to the rare transitions between metastable states. Additionally, MSMs are useful to identify observables that may act as sensors for the conformational or binding state of the molecule, thus guiding the design of experiments. In the present study, the conformation dynamics of UDP-GlcNAc are studied in vacuum, water, water + Mg(2+), and in the protein UDP-GlcNAc 2-epimerase. It is found that addition of Mg(2+) significantly affects the conformational stability, thermodynamics, and kinetics of UDP-GlcNAc. In particular, the slowest structural process, puckering of the GlcNAc sugar, depends on the overall conformation of UDP-GlcNAc and may thus act as a sensor of whether Mg(2+) is bound or not. Interestingly, transferring the molecule from vacuum to water makes the protein-binding conformations UDP-GlcNAc first accessible, while adding Mg(2+) further stabilizes them by specifically associating to binding-competent conformations. While Mg(2+) is not cocrystallized in the UDP-GlcNAc 2-epimerase complex, the selectively stabilized Mg(2+)/UDP-GlcNAc complex may be a template for the bound state, and Mg(2+) may accompany the binding-competent ligand conformation to the binding pocket. This serves as a possible explanation of the enhanced epimerization rate in the presence of Mg(2+). This role of Mg(2+) has previously not been described and opens the question whether "binding co-factors" may be a concept of general relevance for protein-ligand binding.
Collapse
Affiliation(s)
- Martin Held
- Institute of Mathematics, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | | | | | | |
Collapse
|
43
|
Shalashilin DV, Beddard GS, Paci E, Glowacki DR. Peptide kinetics from picoseconds to microseconds using boxed molecular dynamics: Power law rate coefficients in cyclisation reactions. J Chem Phys 2012; 137:165102. [DOI: 10.1063/1.4759088] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Caliandro R, Rossetti G, Carloni P. Local Fluctuations and Conformational Transitions in Proteins. J Chem Theory Comput 2012; 8:4775-85. [DOI: 10.1021/ct300610y] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Rocco Caliandro
- CNR—Institute of Crystallography,
via Amendola 122/o, I-70126, Bari, Italy
| | - Giulia Rossetti
- Institute for Research in Biomedicine
and Barcelona Supercomputing Center, Joint Research Program on Computational
Biology, Baldiri I Reixac 10, 08028, Barcelona, Spain
- Jülich Supercomputing Centre,
Institute for Advanced Simulation, Forschungszentrum Jülich,
D-52425 Jülich, Germany
- Computational Biophysics, German
Research School for Simulation Sciences 1, D-52425 Jülich,
Germany, and Institute for Advanced Simulation, Forschungszentrum
Jülich, D-52425 Jülich, Germany
| | - Paolo Carloni
- Computational Biophysics, German
Research School for Simulation Sciences 1, D-52425 Jülich,
Germany, and Institute for Advanced Simulation, Forschungszentrum
Jülich, D-52425 Jülich, Germany
- Statistical and Biological Physics
Sector, International School for Advanced Studies (SISSA), Trieste,
Italy
| |
Collapse
|
45
|
Senne M, Trendelkamp-Schroer B, Mey ASJS, Schütte C, Noé F. EMMA: A Software Package for Markov Model Building and Analysis. J Chem Theory Comput 2012; 8:2223-38. [PMID: 26588955 DOI: 10.1021/ct300274u] [Citation(s) in RCA: 124] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
The study of folding and conformational changes of macromolecules by molecular dynamics simulations often requires the generation of large amounts of simulation data that are difficult to analyze. Markov (state) models (MSMs) address this challenge by providing a systematic way to decompose the state space of the molecular system into substates and to estimate a transition matrix containing the transition probabilities between these substates. This transition matrix can be analyzed to reveal the metastable, i.e., long-living, states of the system, its slowest relaxation time scales, and transition pathways and rates, e.g., from unfolded to folded, or from dissociated to bound states. Markov models can also be used to calculate spectroscopic data and thus serve as a way to reconcile experimental and simulation data. To reduce the technical burden of constructing, validating, and analyzing such MSMs, we provide the software framework EMMA that is freely available at https://simtk.org/home/emma .
Collapse
Affiliation(s)
- Martin Senne
- Department for Mathematics and Computer Science, FU Berlin
| | | | | | | | - Frank Noé
- Department for Mathematics and Computer Science, FU Berlin
| |
Collapse
|