1
|
Fischer SA, Roeters SJ, Meuzelaar H, Woutersen S, Weidner T, Pfaendtner J. Estimation of vibrational spectra of Trp-cage protein from nonequilibrium metadynamics simulations. Biophys J 2024:S0006-3495(24)00562-9. [PMID: 39180216 DOI: 10.1016/j.bpj.2024.08.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2024] [Revised: 06/28/2024] [Accepted: 08/21/2024] [Indexed: 08/26/2024] Open
Abstract
The development of methods that allow a structural interpretation of linear and nonlinear vibrational spectra is of great importance, both for spectroscopy and for optimizing force field quality. The experimentally measured signals are ensemble averages over all accessible configurations, which complicates spectral calculations. To account for this, we present a recipe for calculating vibrational amide-I spectra of proteins based on metadynamics molecular dynamics simulations. For each frame, a one-exciton Hamiltonian is set up for the backbone amide groups, in which the couplings are estimated with the transition-charge coupling model for nonnearest neighbors, and with a parametrized map of ab initio calculations that give the coupling as a function of the dihedral angles for nearest neighbors. The local-mode frequency variations due to environmental factors such as hydrogen bonds are modeled by exploiting the linear relationship between the amide C-O bond length and the amide-I frequency. The spectra are subsequently calculated while taking into account the equilibrium statistical weights of the frames that are determined using a previously published reweighting procedure. By implementing all these steps in an efficient Fortran code, the spectra can be averaged over very large amounts of structures, thereby extensively covering the phase space of proteins. Using this recipe, the spectral responses of 2.5 million frames of a metadynamics simulation of the miniprotein Trp-cage are averaged to reproduce the experimental temperature-dependent IR spectra very well. The spectral calculations provide new insight into the origin of the various spectral signatures (which are typically challenging to disentangle in the congested amide-I region), and allow for a direct structural interpretation of the experimental spectra and for validation of the molecular dynamics simulations of ensembles.
Collapse
Affiliation(s)
- Sean A Fischer
- Chemistry Division, Naval Research Laboratory, Washington, District of Columbia
| | - Steven J Roeters
- Department of Chemistry, Aarhus University, Aarhus C, Denmark; Amsterdam UMC, Department of Anatomy and Neurosciences, Vrije Universiteit, Amsterdam, the Netherlands
| | - Heleen Meuzelaar
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, the Netherlands
| | - Sander Woutersen
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, the Netherlands
| | - Tobias Weidner
- Department of Chemistry, Aarhus University, Aarhus C, Denmark; Department of Chemical & Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina.
| | - Jim Pfaendtner
- Department of Chemical & Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina.
| |
Collapse
|
2
|
Feng RR, Wang M, Zhang W, Gai F. Unnatural Amino Acids for Biological Spectroscopy and Microscopy. Chem Rev 2024; 124:6501-6542. [PMID: 38722769 DOI: 10.1021/acs.chemrev.3c00944] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/23/2024]
Abstract
Due to advances in methods for site-specific incorporation of unnatural amino acids (UAAs) into proteins, a large number of UAAs with tailored chemical and/or physical properties have been developed and used in a wide array of biological applications. In particular, UAAs with specific spectroscopic characteristics can be used as external reporters to produce additional signals, hence increasing the information content obtainable in protein spectroscopic and/or imaging measurements. In this Review, we summarize the progress in the past two decades in the development of such UAAs and their applications in biological spectroscopy and microscopy, with a focus on UAAs that can be used as site-specific vibrational, fluorescence, electron paramagnetic resonance (EPR), or nuclear magnetic resonance (NMR) probes. Wherever applicable, we also discuss future directions.
Collapse
Affiliation(s)
- Ran-Ran Feng
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Manxi Wang
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Wenkai Zhang
- Department of Physics and Applied Optics Beijing Area Major Laboratory, Beijing Normal University, Beijing 100875, China
| | - Feng Gai
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
3
|
Maruyama Y, Mitsutake A. Effect of Main and Side Chains on the Folding Mechanism of the Trp-Cage Miniprotein. ACS OMEGA 2023; 8:43827-43835. [PMID: 38027385 PMCID: PMC10666239 DOI: 10.1021/acsomega.3c05809] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/07/2023] [Revised: 09/19/2023] [Accepted: 10/27/2023] [Indexed: 12/01/2023]
Abstract
Proteins that do not fold into their functional native state have been linked to diseases. In this study, the influence of the main and side chains of individual amino acids on the folding of the tryptophan cage (Trp-cage), a designed 20-residue miniprotein, was analyzed. For this purpose, we calculated the solvation free energy (SFE) contributions of individual atoms by using the 3D-reference interaction site model with the atomic decomposition method. The mechanism by which the Trp-cage is stabilized during the folding process was examined by calculating the total energy, which is the sum of the conformational energy and SFE. The folding process of the Trp-cage resulted in a stable native state, with a total energy that was 62.4 kcal/mol lower than that of the unfolded state. The solvation entropy, which is considered to be responsible for the hydrophobic effect, contributed 31.3 kcal/mol to structural stabilization. In other words, the contribution of the solvation entropy accounted for approximately half of the total contribution to Trp-cage folding. The hydrophobic core centered on Trp6 contributed 15.6 kcal/mol to the total energy, whereas the solvation entropy contribution was 6.3 kcal/mol. The salt bridge formed by the hydrophilic side chains of Asp9 and Arg16 contributed 10.9 and 5.0 kcal/mol, respectively. This indicates that not only the hydrophobic core but also the salt bridge of the hydrophilic side chains gain solvation entropy and contribute to stabilizing the native structure of the Trp-cage.
Collapse
Affiliation(s)
- Yutaka Maruyama
- Data
Science Center for Creative Design and Manufacturing, The Institute of Statistical Mathematics, 10-3 Midori-cho, Tachikawa, Tokyo 190-8562, Japan
- Department
of Physics, School of Science and Technology, Meiji University, 1-1-1
Higashi-Mita, Tama-ku, Kawasaki-shi, Kanagawa 214-8571, Japan
| | - Ayori Mitsutake
- Department
of Physics, School of Science and Technology, Meiji University, 1-1-1
Higashi-Mita, Tama-ku, Kawasaki-shi, Kanagawa 214-8571, Japan
| |
Collapse
|
4
|
Bajpai S, Petkov BK, Tong M, Abreu CRA, Nair NN, Tuckerman ME. An interoperable implementation of collective-variable based enhanced sampling methods in extended phase space within the OpenMM package. J Comput Chem 2023; 44:2166-2183. [PMID: 37464902 DOI: 10.1002/jcc.27182] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2023] [Revised: 05/30/2023] [Accepted: 06/06/2023] [Indexed: 07/20/2023]
Abstract
Collective variable (CV)-based enhanced sampling techniques are widely used today for accelerating barrier-crossing events in molecular simulations. A class of these methods, which includes temperature accelerated molecular dynamics (TAMD)/driven-adiabatic free energy dynamics (d-AFED), unified free energy dynamics (UFED), and temperature accelerated sliced sampling (TASS), uses an extended variable formalism to achieve quick exploration of conformational space. These techniques are powerful, as they enhance the sampling of a large number of CVs simultaneously compared to other techniques. Extended variables are kept at a much higher temperature than the physical temperature by ensuring adiabatic separation between the extended and physical subsystems and employing rigorous thermostatting. In this work, we present a computational platform to perform extended phase space enhanced sampling simulations using the open-source molecular dynamics engine OpenMM. The implementation allows users to have interoperability of sampling techniques, as well as employ state-of-the-art thermostats and multiple time-stepping. This work also presents protocols for determining the critical parameters and procedures for reconstructing high-dimensional free energy surfaces. As a demonstration, we present simulation results on the high dimensional conformational landscapes of the alanine tripeptide in vacuo, tetra-N-methylglycine (tetra-sarcosine) peptoid in implicit solvent, and the Trp-cage mini protein in explicit water.
Collapse
Affiliation(s)
- Shitanshu Bajpai
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur, India
| | - Brian K Petkov
- Department of Chemistry, New York University (NYU), New York, New York, USA
| | - Muchen Tong
- Department of Chemistry, New York University (NYU), New York, New York, USA
| | - Charlles R A Abreu
- Chemical Engineering Department, Escola de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
| | - Nisanth N Nair
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur, India
| | - Mark E Tuckerman
- Department of Chemistry, New York University (NYU), New York, New York, USA
- Courant Institute of Mathematical Sciences, New York University (NYU), New York, New York, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China
- Simons Center for Computational Physical Chemistry, New York University, New York, New York, USA
| |
Collapse
|
5
|
Sun Q, He X, Fu Y. The "Beacon" Structural Model of Protein Folding: Application for Trp-Cage in Water. Molecules 2023; 28:5164. [PMID: 37446826 DOI: 10.3390/molecules28135164] [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: 05/30/2023] [Revised: 06/30/2023] [Accepted: 06/30/2023] [Indexed: 07/15/2023] Open
Abstract
Protein folding is a process in which a polypeptide must undergo folding process to obtain its three-dimensional structure. Thermodynamically, it is a process of enthalpy to overcome the loss of conformational entropy in folding. Folding is primarily related to hydrophobic interactions and intramolecular hydrogen bondings. During folding, hydrophobic interactions are regarded to be the driving forces, especially in the initial structural collapse of a protein. Additionally, folding is guided by the strong interactions within proteins, such as intramolecular hydrogen bondings related to the α-helices and β-sheets of proteins. Therefore, a protein is divided into the folding key (FK) regions related to intramolecular hydrogen bondings and the non-folding key (non-FK) regions. Various conformations are expected for FK and non-FK regions. Different from non-FK regions, it is necessary for FK regions to form the specific conformations in folding, which are regarded as the necessary folding pathways (or "beacons"). Additionally, sequential folding is expected for the FK regions, and the intermediate state is found during folding. They are reflected on the local basins in the free energy landscape (FEL) of folding. To demonstrate the structural model, molecular dynamics (MD) simulations are conducted on the folding pathway of the TRP-cage in water.
Collapse
Affiliation(s)
- Qiang Sun
- Key Laboratory of Orogenic Belts and Crustal Evolution, Ministry of Education, The School of Earth and Space Sciences, Peking University, Beijing 100871, China
| | - Xian He
- Key Laboratory of Orogenic Belts and Crustal Evolution, Ministry of Education, The School of Earth and Space Sciences, Peking University, Beijing 100871, China
| | - Yanfang Fu
- Key Laboratory of Orogenic Belts and Crustal Evolution, Ministry of Education, The School of Earth and Space Sciences, Peking University, Beijing 100871, China
| |
Collapse
|
6
|
Gerig JT. Examination of Solvent Interactions with Trp-Cage in 1,1,1,3,3,3-Hexafluoro-2-propanol-water at 298 K through MD Simulations and Intermolecular Nuclear Overhauser Effects. J Phys Chem B 2023; 127:5062-5071. [PMID: 37249321 PMCID: PMC10258800 DOI: 10.1021/acs.jpcb.3c01029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 05/08/2023] [Indexed: 05/31/2023]
Abstract
MD simulations of the peptide Trp-cage dissolved in 28% hexafluoroisopropanol (HFIP)-water have been carried out at 298 K with the goal of exploring peptide hydrogen-solvent fluorine nuclear spin cross-relaxation. The work was motivated by the observation that most experimental fluoroalcohol-peptide cross-relaxation terms at 298 K are small, both positive and negative, and not always well predicted from simulations. The cross-relaxation terms for hydrogens of the caged tryptophan residue of Trp-cage are substantially negative, a result consistent with simulations. It was concluded that hexafluoroisopropanol interactions near this part of the peptide are particularly long-lived. While both HFIP and water are present in all regions of the simulation box, the composition of the solvent mixture is not homogeneous throughout the system. HFIP generally accumulates near the peptide surface, while water molecules are preferentially found in regions that are more than 1.5 nm from the surface of the peptide. However, some water remains in higher-than-expected amounts in the solvent layer surrounding 6Trp, 9Asp, Ser13, and Ser14 residues in the helical region of Trp-cage. As observed in simulations of this system at 278 K, HFIP molecules aggregate into clusters that continually form and re-form. Translational diffusion of both HFIP and water appears to be slowed near the surface of the peptide with reduction in diffusion near the 6Trp residue 2- to 3-fold larger than calculated for solvent interactions with other regions of Trp-cage.
Collapse
Affiliation(s)
- J. T. Gerig
- Department of Chemistry &
Biochemistry, University of California,
Santa Barbara, Santa
Barbara, California 93106, United States
| |
Collapse
|
7
|
Shityakov S, Skorb EV, Nosonovsky M. Folding-unfolding asymmetry and a RetroFold computational algorithm. ROYAL SOCIETY OPEN SCIENCE 2023; 10:221594. [PMID: 37153361 PMCID: PMC10154942 DOI: 10.1098/rsos.221594] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/23/2022] [Accepted: 03/30/2023] [Indexed: 05/09/2023]
Abstract
We treat protein folding as molecular self-assembly, while unfolding is viewed as disassembly. Fracture is typically a much faster process than self-assembly. Self-assembly is often an exponentially decaying process, since energy relaxes due to dissipation, while fracture is a constant-rate process as the driving force is opposed by damping. Protein folding takes two orders of magnitude longer than unfolding. We suggest a mathematical transformation of variables, which makes it possible to view self-assembly as time-reversed disassembly, thus folding can be studied as reversed unfolding. We investigate the molecular dynamics modelling of folding and unfolding of the short Trp-cage protein. Folding time constitutes about 800 ns, while unfolding (denaturation) takes only about 5.0 ns and, therefore, fewer computational resources are needed for its simulation. This RetroFold approach can be used for the design of a novel computation algorithm, which, while approximate, is less time-consuming than traditional folding algorithms.
Collapse
Affiliation(s)
- Sergey Shityakov
- Infochemistry Scientific Center (ISC), ITMO University, 9 Lomonosova Street, St. Petersburg 191002, Russia
| | - Ekaterina V. Skorb
- Infochemistry Scientific Center (ISC), ITMO University, 9 Lomonosova Street, St. Petersburg 191002, Russia
| | - Michael Nosonovsky
- Infochemistry Scientific Center (ISC), ITMO University, 9 Lomonosova Street, St. Petersburg 191002, Russia
- College of Engineering and Applied Science, University of Wisconsin-Milwaukee, Milwaukee, WI 53211, USA
| |
Collapse
|
8
|
Franke L, Peter C. Visualizing the Residue Interaction Landscape of Proteins by Temporal Network Embedding. J Chem Theory Comput 2023; 19:2985-2995. [PMID: 37122117 DOI: 10.1021/acs.jctc.2c01228] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
Characterizing the structural dynamics of proteins with heterogeneous conformational landscapes is crucial to understanding complex biomolecular processes. To this end, dimensionality reduction algorithms are used to produce low-dimensional embeddings of the high-dimensional conformational phase space. However, identifying a compact and informative set of input features for the embedding remains an ongoing challenge. Here, we propose to harness the power of Residue Interaction Networks (RINs) and their centrality measures, established tools to provide a graph theoretical view on molecular structure. Specifically, we combine the closeness centrality, which captures global features of the protein conformation at residue-wise resolution, with EncoderMap, a hybrid neural-network autoencoder/multidimensional-scaling like dimensionality reduction algorithm. We find that the resulting low-dimensional embedding is a meaningful visualization of the residue interaction landscape that resolves structural details of the protein behavior while retaining global interpretability. This feature-based graph embedding of temporal protein graphs makes it possible to apply the general descriptive power of RIN formalisms to the analysis of protein simulations of complex processes such as protein folding and multidomain interactions requiring no protein-specific input. We demonstrate this on simulations of the fast folding protein Trp-Cage and the multidomain signaling protein FAT10. Due to its generality and modularity, the presented approach can easily be transferred to other protein systems.
Collapse
Affiliation(s)
- Leon Franke
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, Konstanz 78457, Germany
- Konstanz Research School Chemical Biology, University of Konstanz, Universitätsstraße 10, Konstanz 78457, Germany
| | - Christine Peter
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, Konstanz 78457, Germany
| |
Collapse
|
9
|
Chan AM, Nijhawan AK, Hsu DJ, Leshchev D, Rimmerman D, Kosheleva I, Kohlstedt KL, Chen LX. The Role of Transient Intermediate Structures in the Unfolding of the Trp-Cage Fast-Folding Protein: Generating Ensembles from Time-Resolved X-ray Solution Scattering with Genetic Algorithms. J Phys Chem Lett 2023; 14:1133-1139. [PMID: 36705525 PMCID: PMC10167713 DOI: 10.1021/acs.jpclett.2c03680] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
The Trp-cage miniprotein is one of the smallest systems to exhibit a stable secondary structure and fast-folding dynamics, serving as an apt model system to study transient intermediates with both experimental and computational analyses. Previous spectroscopic characterizations that have been done on Trp-cage have inferred a single stable intermediate on a pathway from folded to unfolded basins. We aim to bridge the understanding of Trp-cage structural folding dynamics on microsecond-time scales, by utilizing time-resolved X-ray solution scattering to probe the temperature-induced unfolding pathway. Our results indicate the formation of a conformationally extended intermediate on the time scale of 1 μs, which undergoes complete unfolding within 5 μs. We further investigated the atomistic structural details of the unfolding pathway using a genetic algorithm to generate ensemble model fits to the scattering profiles. This analysis paves the way for direct benchmarking of theoretical models of protein folding ensembles produced with molecular dynamics simulations.
Collapse
Affiliation(s)
- Arnold M Chan
- Department of Chemistry, Northwestern University, Evanston, Illinois60208, United States
| | - Adam K Nijhawan
- Department of Chemistry, Northwestern University, Evanston, Illinois60208, United States
| | - Darren J Hsu
- Department of Chemistry, Northwestern University, Evanston, Illinois60208, United States
| | - Denis Leshchev
- Department of Chemistry, Northwestern University, Evanston, Illinois60208, United States
| | - Dolev Rimmerman
- Department of Chemistry, Northwestern University, Evanston, Illinois60208, United States
| | - Irina Kosheleva
- Center for Advanced Radiation Sources, The University of Chicago, Chicago, Illinois60637, United States
| | - Kevin L Kohlstedt
- Department of Chemistry, Northwestern University, Evanston, Illinois60208, United States
| | - Lin X Chen
- Department of Chemistry, Northwestern University, Evanston, Illinois60208, United States
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois60439, United States
| |
Collapse
|
10
|
Shityakov S, Skorb EV, Nosonovsky M. Topological bio-scaling analysis as a universal measure of protein folding. ROYAL SOCIETY OPEN SCIENCE 2022; 9:220160. [PMID: 35845855 PMCID: PMC9277272 DOI: 10.1098/rsos.220160] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/21/2022] [Accepted: 06/14/2022] [Indexed: 05/24/2023]
Abstract
Scaling relationships for polymeric molecules establish power law dependencies between the number of molecular segments and linear dimensions, such as the radius of gyration. They also establish spatial topological properties of the chains, such as their dimensionality. In the spatial domain, power exponents α = 1 (linear stretched molecule), α = 0.5 (the ideal chain) and α = 0.333 (compact globule) are significant. During folding, the molecule undergoes the transition from the one-dimensional linear to the three-dimensional globular state within a very short time. However, intermediate states with fractional dimensions can be stabilized by modifying the solubility (e.g. by changing the solution temperature). Topological properties, such as dimension, correlate with the interaction energy, and thus by tuning the solubility one can control molecular interaction. We investigate these correlations using the example of a well-studied short model of Trp-cage protein. The radius of gyration is used to estimate the fractal dimension of the chain at different stages of folding. It is expected that the same principle is applicable to much larger molecules and that topological (dimensional) characteristics can provide insights into molecular folding and interactions.
Collapse
Affiliation(s)
- Sergey Shityakov
- Infochemistry Scientific Center (ISC), ITMO University, 9 Lomonosova St., St Petersburg 191002, Russia
| | - Ekaterina V. Skorb
- Infochemistry Scientific Center (ISC), ITMO University, 9 Lomonosova St., St Petersburg 191002, Russia
| | - Michael Nosonovsky
- Infochemistry Scientific Center (ISC), ITMO University, 9 Lomonosova St., St Petersburg 191002, Russia
| |
Collapse
|
11
|
Ghorbani M, Prasad S, Klauda JB, Brooks BR. GraphVAMPNet, using graph neural networks and variational approach to Markov processes for dynamical modeling of biomolecules. J Chem Phys 2022; 156:184103. [PMID: 35568532 PMCID: PMC9094994 DOI: 10.1063/5.0085607] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2022] [Accepted: 04/22/2022] [Indexed: 11/14/2022] Open
Abstract
Finding a low dimensional representation of data from long-timescale trajectories of biomolecular processes, such as protein folding or ligand-receptor binding, is of fundamental importance, and kinetic models, such as Markov modeling, have proven useful in describing the kinetics of these systems. Recently, an unsupervised machine learning technique called VAMPNet was introduced to learn the low dimensional representation and the linear dynamical model in an end-to-end manner. VAMPNet is based on the variational approach for Markov processes and relies on neural networks to learn the coarse-grained dynamics. In this paper, we combine VAMPNet and graph neural networks to generate an end-to-end framework to efficiently learn high-level dynamics and metastable states from the long-timescale molecular dynamics trajectories. This method bears the advantages of graph representation learning and uses graph message passing operations to generate an embedding for each datapoint, which is used in the VAMPNet to generate a coarse-grained dynamical model. This type of molecular representation results in a higher resolution and a more interpretable Markov model than the standard VAMPNet, enabling a more detailed kinetic study of the biomolecular processes. Our GraphVAMPNet approach is also enhanced with an attention mechanism to find the important residues for classification into different metastable states.
Collapse
Affiliation(s)
| | - Samarjeet Prasad
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, Maryland 20824, USA
| | - Jeffery B. Klauda
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| | - Bernard R. Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, Maryland 20824, USA
| |
Collapse
|
12
|
Bò L, Milanetti E, Chen CG, Ruocco G, Amadei A, D’Abramo M. Computational Modeling of the Thermodynamics of the Mesophilic and Thermophilic Mutants of Trp-Cage Miniprotein. ACS OMEGA 2022; 7:13448-13454. [PMID: 35559192 PMCID: PMC9088802 DOI: 10.1021/acsomega.1c06206] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Accepted: 03/10/2022] [Indexed: 06/15/2023]
Abstract
We characterize the folding-unfolding thermodynamics of two mutants of the miniprotein Trp-cage by combining extended molecular dynamics simulations and an advanced statistical-mechanical-based approach. From a set of molecular dynamics simulations in an explicit solvent performed along a reference isobar, we evaluated the structural and thermodynamic behaviors of a mesophilic and a thermophilic mutant of the Trp-cage and their temperature dependence. In the case of the thermophilic mutant, computational data confirm that our theoretical-computational approach is able to reproduce the available experimental estimate with rather good accuracy. On the other hand, the mesophilic mutant does not show a clear two-state (folded and unfolded) behavior, preventing us from reconstructing its thermodynamics; thus, an analysis of its structural behavior along a reference isobar is presented. Our results show that an extended sampling of these kinds of systems coupled to an advanced statistical-mechanical-based treatment of the data can provide an accurate description of the folding-unfolding thermodynamics along a reference isobar, rationalizing the discrepancies between the simulated and experimental systems.
Collapse
Affiliation(s)
- Leonardo Bò
- Department
of Physics, Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
| | - Edoardo Milanetti
- Department
of Physics, Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
- Center
for Life Nano & Neuroscience, Italian
Institute of Technology, Viale Regina Elena 291, 00161 Rome, Italy
| | - Cheng Giuseppe Chen
- Department
of Chemistry, Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
| | - Giancarlo Ruocco
- Department
of Physics, Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
- Center
for Life Nano & Neuroscience, Italian
Institute of Technology, Viale Regina Elena 291, 00161 Rome, Italy
| | - Andrea Amadei
- Department
of Chemical Sciences and Technology, Universitá
degli Studi di Roma Tor Vergata, Via della ricerca scientifica 00133 Rome, Italy
| | - Marco D’Abramo
- Department
of Chemistry, Sapienza University, Piazzale Aldo Moro 5, 00185 Rome, Italy
| |
Collapse
|
13
|
van Adrichem KE, Jansen TLC. AIM: A Mapping Program for Infrared Spectroscopy of Proteins. J Chem Theory Comput 2022; 18:3089-3098. [PMID: 35387451 PMCID: PMC9097285 DOI: 10.1021/acs.jctc.2c00113] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
Here, we present
a new analysis program, AIM, that allows extracting
the vibrational amide-I Hamiltonian using molecular dynamics trajectories
for protein infrared spectroscopy modeling. The constructed Hamiltonians
can be used as input for spectral calculations allowing the calculation
of infrared absorption spectra, vibrational circular dichroism, and
two-dimensional infrared spectra. These spectroscopies allow the study
of the structure and dynamics of proteins. We will explain the essence
of how AIM works and give examples of the information and spectra
that can be obtained with the program using the Trypsin Inhibitor
as an example. AIM is freely available from GitHub, and the package
contains a demonstration allowing easy introduction to the use of
the program.
Collapse
Affiliation(s)
- Kim E van Adrichem
- Zernike Institute for Advanced Materials, University of Groningen, 9747 AG Groningen, The Netherlands
| | - Thomas L C Jansen
- Zernike Institute for Advanced Materials, University of Groningen, 9747 AG Groningen, The Netherlands
| |
Collapse
|
14
|
Kapakayala AB, Nair NN. Boosting the conformational sampling by combining replica exchange with solute tempering and well-sliced metadynamics. J Comput Chem 2021; 42:2233-2240. [PMID: 34585768 DOI: 10.1002/jcc.26752] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2021] [Revised: 08/30/2021] [Accepted: 09/12/2021] [Indexed: 01/22/2023]
Abstract
Methods that combine collective variable (CV) based enhanced sampling and global tempering approaches are used in speeding-up the conformational sampling and free energy calculation of large and soft systems with a plethora of energy minima. In this paper, a new method of this kind is proposed in which the well-sliced metadynamics approach (WSMTD) is united with replica exchange with solute tempering (REST2) method. WSMTD employs a divide-and-conquer strategy wherein high-dimensional slices of a free energy surface are independently sampled and combined. The method enables one to accomplish a controlled exploration of the CV-space with a restraining bias as in umbrella sampling, and enhance-sampling of one or more orthogonal CVs using a metadynamics like bias. The new hybrid method proposed here enables boosting the sampling of more slow degrees of freedom in WSMTD simulations, without the need to specify associated CVs, through a replica exchange scheme within the framework of REST2. The high-dimensional slices of the probability distributions of CVs computed from the united WSMTD and REST2 simulations are subsequently combined using the weighted histogram analysis method to obtain the free energy surface. We show that the new method proposed here is accurate, improves the conformational sampling, and achieves quick convergence in free energy estimates. We demonstrate this by computing the conformational free energy landscapes of solvated alanine tripeptide and Trp-cage mini protein in explicit water.
Collapse
Affiliation(s)
- Anji Babu Kapakayala
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, India.,School of Pharmacy and Biomedical Sciences, Curtin University, Perth, Australia
| | - Nisanth N Nair
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, India
| |
Collapse
|
15
|
Ghorbani M, Prasad S, Klauda JB, Brooks BR. Variational embedding of protein folding simulations using Gaussian mixture variational autoencoders. J Chem Phys 2021; 155:194108. [PMID: 34800961 DOI: 10.1063/5.0069708] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Conformational sampling of biomolecules using molecular dynamics simulations often produces a large amount of high dimensional data that makes it difficult to interpret using conventional analysis techniques. Dimensionality reduction methods are thus required to extract useful and relevant information. Here, we devise a machine learning method, Gaussian mixture variational autoencoder (GMVAE), that can simultaneously perform dimensionality reduction and clustering of biomolecular conformations in an unsupervised way. We show that GMVAE can learn a reduced representation of the free energy landscape of protein folding with highly separated clusters that correspond to the metastable states during folding. Since GMVAE uses a mixture of Gaussians as its prior, it can directly acknowledge the multi-basin nature of the protein folding free energy landscape. To make the model end-to-end differentiable, we use a Gumbel-softmax distribution. We test the model on three long-timescale protein folding trajectories and show that GMVAE embedding resembles the folding funnel with folded states down the funnel and unfolded states outside the funnel path. Additionally, we show that the latent space of GMVAE can be used for kinetic analysis and Markov state models built on this embedding produce folding and unfolding timescales that are in close agreement with other rigorous dynamical embeddings such as time independent component analysis.
Collapse
Affiliation(s)
- Mahdi Ghorbani
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, Maryland 20824, USA
| | - Samarjeet Prasad
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, Maryland 20824, USA
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, Maryland 20824, USA
| |
Collapse
|
16
|
Nijhawan AK, Chan AM, Hsu DJ, Chen LX, Kohlstedt KL. Resolving Dynamics in the Ensemble: Finding Paths through Intermediate States and Disordered Protein Structures. J Phys Chem B 2021; 125:12401-12412. [PMID: 34748336 PMCID: PMC9096987 DOI: 10.1021/acs.jpcb.1c05820] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Proteins have been found to inhabit a diverse set of three-dimensional structures. The dynamics that govern protein interconversion between structures happen over a wide range of time scales─picoseconds to seconds. Our understanding of protein functions and dynamics is largely reliant upon our ability to elucidate physically populated structures. From an experimental structural characterization perspective, we are often limited to measuring the ensemble-averaged structure both in the steady-state and time-resolved regimes. Generating kinetic models and understanding protein structure-function relationships require atomistic knowledge of the populated states in the ensemble. In this Perspective, we present ensemble refinement methodologies that integrate time-resolved experimental signals with molecular dynamics models. We first discuss integration of experimental structural restraints to molecular models in disordered protein systems that adhere to the principle of maximum entropy for creating a complete set of ensemble structures. We then propose strategies to find kinetic pathways between the refined structures, using time-resolved inputs to guide molecular dynamics trajectories and the use of inference to generate tailored stimuli to prepare a desired ensemble of protein states.
Collapse
Affiliation(s)
- Adam K Nijhawan
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Arnold M Chan
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Darren J Hsu
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Lin X Chen
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Kevin L Kohlstedt
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| |
Collapse
|
17
|
Peter EK, Manstein DJ, Shea JE, Schug A. CORE-MD II: A fast, adaptive, and accurate enhanced sampling method. J Chem Phys 2021; 155:104114. [PMID: 34525829 DOI: 10.1063/5.0063664] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we present a fast and adaptive correlation guided enhanced sampling method (CORE-MD II). The CORE-MD II technique relies, in part, on partitioning of the entire pathway into short trajectories that we refer to as instances. The sampling within each instance is accelerated by adaptive path-dependent metadynamics simulations. The second part of this approach involves kinetic Monte Carlo (kMC) sampling between the different states that have been accessed during each instance. Through the combination of the partition of the total simulation into short non-equilibrium simulations and the kMC sampling, the CORE-MD II method is capable of sampling protein folding without any a priori definitions of reaction pathways and additional parameters. In the validation simulations, we applied the CORE-MD II on the dialanine peptide and the folding of two peptides: TrpCage and TrpZip2. In a comparison with long time equilibrium Molecular Dynamics (MD), 1 µs replica exchange MD (REMD), and CORE-MD I simulations, we find that the level of convergence of the CORE-MD II method is improved by a factor of 8.8, while the CORE-MD II method reaches acceleration factors of ∼120. In the CORE-MD II simulation of TrpZip2, we observe the formation of the native state in contrast to the REMD and the CORE-MD I simulations. The method is broadly applicable for MD simulations and is not restricted to simulations of protein folding or even biomolecules but also applicable to simulations of protein aggregation, protein signaling, or even materials science simulations.
Collapse
Affiliation(s)
- Emanuel K Peter
- Institute for Biophysical Chemistry, Fritz-Hartmann-Centre for Medical Research, Hannover Medical School, Carl-Neuberg-Str. 1, Hannover 30625, Germany
| | - Dietmar J Manstein
- Institute for Biophysical Chemistry, Fritz-Hartmann-Centre for Medical Research, Hannover Medical School, Carl-Neuberg-Str. 1, Hannover 30625, Germany
| | - Joan-Emma Shea
- Department of Chemistry and Biochemistry, Department of Physics, University of California, Santa Barbara, California 93106, USA
| | - Alexander Schug
- John von Neumann Institute for Computing and Jülich Supercomputing Centre, Institute for Advanced Simulation, Forschungszentrum Jülich, 52425 Jülich, Germany
| |
Collapse
|
18
|
Gerig JT. Examination of Interactions of Hexafluoro-2-propanol with Trp-Cage in Hexafluoro-2-propanol-Water by MD Simulations and Intermolecular Nuclear Overhauser Effects. J Phys Chem B 2020; 124:9793-9802. [PMID: 33095591 DOI: 10.1021/acs.jpcb.0c06476] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
All-atom molecular dynamic simulations of the peptide Trp-cage in 30% hexafluoro-2-propanol- water (V/V) at 278 K have been carried out with the goal of exploring peptide hydrogen-solvent fluorine nuclear spin cross relaxation. Force field parameters for HFIP reported by Fioroni et al. along with the fluorine parameters of the TFE5 model reported by this lab were used. Water was represented by the TIP5P-Ew model. Peptide modeling used the AMBER99SB-ILDN force field. Translational diffusion coefficients of solution components at 278 K were predicted to within 35% of experimental values using these parameter sets. The simulations indicate that the solvent mixture is not homogeneous, with HFIP molecules clustered into aggregates as large as 53 fluoroalcohol molecules. The solvent environment of surface atoms of Trp-cage fluctuates between being HFIP-rich and more water-rich about every 10 ns. In accord with previous studies by other groups, the average concentration of HFIP near the surface of the peptide is significantly enhanced over the concentration of HFIP in the bulk solvent. In the simulations, ∼7% of the initial contacts between HFIP molecules and Trp-cage develop into peptide-fluoroalcohol interactions that persist for times as long as 8 ns. Most of the available experimental nuclear spin cross-relaxation rates (ΣHF) for hydrogens of the Trp-cage in 30% HFIP-water are reproduced from the MD trajectories to within uncertainties of the experimental data and the simulations. However, a few calculated ΣHF values for hydrogens of the Trp-cage do not agree with experiment. These tend to be situations where long-lived peptide-HFIP interactions are predicted. The disagreements between observed and calculated ΣHF in these instances signal defects in the modeling parameters and procedures that are presently unrecognized.
Collapse
Affiliation(s)
- J T Gerig
- Department of Chemistry & Biochemistry University of California, Santa Barbara Santa Barbara, California 93106, United States
| |
Collapse
|
19
|
Peter EK, Shea JE, Schug A. CORE-MD, a path correlated molecular dynamics simulation method. J Chem Phys 2020; 153:084114. [DOI: 10.1063/5.0015398] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Emanuel K. Peter
- John von Neumann Institute for Computing and Julich Supercomputing Centre, Institute for Advanced Simulation, Forschungszentrum Jülich, Jülich, Germany
| | - Joan-Emma Shea
- Department of Chemistry and Biochemistry, Department of Physics, University of California, Santa Barbara, Santa Barbara, California 93106, USA
| | - Alexander Schug
- John von Neumann Institute for Computing and Julich Supercomputing Centre, Institute for Advanced Simulation, Forschungszentrum Jülich, Jülich, Germany
- Faculty of Biology, University of Duisburg-Essen, Duisburg, Germany
| |
Collapse
|
20
|
Preußke N, Moormann W, Bamberg K, Lipfert M, Herges R, Sönnichsen FD. Visible-light-driven photocontrol of the Trp-cage protein fold by a diazocine cross-linker. Org Biomol Chem 2020; 18:2650-2660. [PMID: 32207764 DOI: 10.1039/c9ob02442e] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Diazocines are characterized by extraordinary photochemical properties rendering them of particular interest for switching the conformation of biomolecules with visible light. Current developments afford synthetic access to unprecedented diazocine derivatives promising particular opportunities in photocontrol of proteins and biological systems. In this work, the well-established approach of photocontrolling the secondary structure of α-helices was exploited using a diazocine to reversibly fold and unfold the tertiary structure of a small protein. The protein of choice was the globulary folded Trp-cage, a widely used model system for the elucidation of protein folding pathways. A specifically designed, short and rigid dicarboxy-functionalized diazocine-based cross-linker was attached to two solvent-exposed side chains at the α-helix of the miniprotein through the use of a primary amine-selective active ester. This cross-linking strategy is orthogonal to the common cysteine-based chemistry. The cross-linked Trp-cage was successfully photoisomerized and exhibited a strong correlation between protein fold and diazocine isomeric state. As determined by NMR spectroscopy, the cis-isomer stabilized the fold, while the trans-isomer led to complete protein unfolding. The successful switching of the protein fold in principle demonstrates the ability to control protein function, as the activity depends on their structural integrity.
Collapse
Affiliation(s)
- Nils Preußke
- Otto-Diels-Institute for Organic Chemistry, Christian-Albrechts-University of Kiel, Otto-Hahn-Platz 4, 24118 Kiel, Germany.
| | | | | | | | | | | |
Collapse
|
21
|
Chalyavi F, Schmitz AJ, Tucker MJ. Unperturbed Detection of the Dynamic Structure in the Hydrophobic Core of Trp-Cage via Two-Dimensional Infrared Spectroscopy. J Phys Chem Lett 2020; 11:832-837. [PMID: 31931573 PMCID: PMC7026909 DOI: 10.1021/acs.jpclett.9b03706] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
The tyrosine ring mode is an intrinsic non-perturbing site-specific infrared reporter for conformational dynamics within protein systems. This transition is influenced by direct and indirect interactions associated with the electron-donating ability and the hydrophobicity of the surrounding molecules. Utilizing an intrinsic tyrosine moiety, two-dimensional infrared spectra of Trp-cage, often called the "hydrogen atom" of protein folding, were measured in the folded and denatured states to uncover the dynamics of the hydrophobic core. The vibrational lifetimes and the correlation decays of the tyrosine ring mode showed significant changes upon both temperature and chemical denaturation of the Trp-cage miniprotein, indicating important structural features of the hydrophobic core and its dynamics. The observed Trp6-Tyr3 interactions are in good agreement with the prior studies of the folded state, but they reach beyond the static structure. These stacking interactions and orientations fluctuate on the picosecond time scale as measured through the spectral dephasing within a dehydrated environment.
Collapse
Affiliation(s)
- Farzaneh Chalyavi
- Department of Chemistry , University of Nevada, Reno , Reno , Nevada 89557 , United States
| | - Andrew J Schmitz
- Department of Chemistry , University of Nevada, Reno , Reno , Nevada 89557 , United States
| | - Matthew J Tucker
- Department of Chemistry , University of Nevada, Reno , Reno , Nevada 89557 , United States
| |
Collapse
|
22
|
Feng CJ, Dhayalan B, Tokmakoff A. Refinement of Peptide Conformational Ensembles by 2D IR Spectroscopy: Application to Ala‒Ala‒Ala. Biophys J 2019; 114:2820-2832. [PMID: 29925019 DOI: 10.1016/j.bpj.2018.05.003] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2018] [Revised: 04/27/2018] [Accepted: 05/02/2018] [Indexed: 10/28/2022] Open
Abstract
Characterizing ensembles of intrinsically disordered proteins is experimentally challenging because of the ill-conditioned nature of ensemble determination with limited data and the intrinsic fast dynamics of the conformational ensemble. Amide I two-dimensional infrared (2D IR) spectroscopy has picosecond time resolution to freeze structural ensembles as needed for probing disordered-protein ensembles and conformational dynamics. Also, developments in amide I computational spectroscopy now allow a quantitative and direct prediction of amide I spectra based on conformational distributions drawn from molecular dynamics simulations, providing a route to ensemble refinement against experimental spectra. We performed a Bayesian ensemble refinement method on Ala-Ala-Ala against isotope-edited Fourier-transform infrared spectroscopy and 2D IR spectroscopy and tested potential factors affecting the quality of ensemble refinements. We found that isotope-edited 2D IR spectroscopy provides a stringent constraint on Ala-Ala-Ala conformations and returns consistent conformational ensembles with the dominant ppII conformer across varying prior distributions from many molecular dynamics force fields and water models. The dominant factor influencing ensemble refinements is the systematic frequency uncertainty from spectroscopic maps. However, the uncertainty of conformer populations can be significantly reduced by incorporating 2D IR spectra in addition to traditional Fourier-transform infrared spectra. Bayesian ensemble refinement against isotope-edited 2D IR spectroscopy thus provides a route to probe equilibrium-complex protein ensembles and potentially nonequilibrium conformational dynamics.
Collapse
Affiliation(s)
- Chi-Jui Feng
- Department of Chemistry, James Franck Institute and Institute for Biophysical Dynamics, University of Chicago, Chicago, Illinois
| | - Balamurugan Dhayalan
- Department of Chemistry, James Franck Institute and Institute for Biophysical Dynamics, University of Chicago, Chicago, Illinois
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute and Institute for Biophysical Dynamics, University of Chicago, Chicago, Illinois.
| |
Collapse
|
23
|
Graham KA, Byrne A, Mason M, Andersen NH. Optimizing the fold stability of the circularly permuted Trp-cage motif. Biopolymers 2019; 110:e23327. [PMID: 31479150 DOI: 10.1002/bip.23327] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2019] [Revised: 08/11/2019] [Accepted: 08/13/2019] [Indexed: 11/09/2022]
Abstract
Through optimization of the linker region and key stabilizing mutations, it has been possible to improve the stability of the circularly permuted (cp) Trp-cage miniprotein. However, even the most stable Trp-cage circular permutants are still less stable than the analogous standard topology (std) Trp-cages. Extending mutational studies of Trp-cage fold stability to cp-species, including analogs lacking chain terminal charges, has uncovered and quantitated some additional stabilizing and destabilizing interactions. Upon protonation, the circular permutants are destabilized to a much greater extent than the standard topology series. End effects, particularly Coulombic interactions, appear to be more important for the cp-series while the Y10/P4 interaction in the cp-series is not as significant a stabilizing feature as the corresponding Y3/P19 in the standard topology series.
Collapse
Affiliation(s)
- Katherine A Graham
- Department of Chemistry, University of Washington, Seattle, Washington, United States of America
| | - Aimee Byrne
- Department of Chemistry, University of Washington, Seattle, Washington, United States of America
| | - Micheal Mason
- Department of Chemistry, University of Washington, Seattle, Washington, United States of America
| | - Niels H Andersen
- Department of Chemistry, University of Washington, Seattle, Washington, United States of America
| |
Collapse
|
24
|
Sidky H, Chen W, Ferguson AL. High-Resolution Markov State Models for the Dynamics of Trp-Cage Miniprotein Constructed Over Slow Folding Modes Identified by State-Free Reversible VAMPnets. J Phys Chem B 2019; 123:7999-8009. [DOI: 10.1021/acs.jpcb.9b05578] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Hythem Sidky
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Wei Chen
- Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States
| | - Andrew L. Ferguson
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
25
|
Shin K, Tran DP, Takemura K, Kitao A, Terayama K, Tsuda K. Enhancing Biomolecular Sampling with Reinforcement Learning: A Tree Search Molecular Dynamics Simulation Method. ACS OMEGA 2019; 4:13853-13862. [PMID: 31497702 PMCID: PMC6714528 DOI: 10.1021/acsomega.9b01480] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/21/2019] [Accepted: 08/01/2019] [Indexed: 05/13/2023]
Abstract
This paper proposes a novel molecular simulation method, called tree search molecular dynamics (TS-MD), to accelerate the sampling of conformational transition pathways, which require considerable computation. In TS-MD, a tree search algorithm, called upper confidence bounds for trees, which is a type of reinforcement learning algorithm, is applied to sample the transition pathway. By learning from the results of the previous simulations, TS-MD efficiently searches conformational space and avoids being trapped in local stable structures. TS-MD exhibits better performance than parallel cascade selection molecular dynamics, which is one of the state-of-the-art methods, for the folding of miniproteins, Chignolin and Trp-cage, in explicit water.
Collapse
Affiliation(s)
- Kento Shin
- Graduate School
of Frontier Sciences, The University of
Tokyo, 5-1-5 Kashiwa-no-ha, Kashiwa, Chiba 277-8561, Japan
| | - Duy Phuoc Tran
- Graduate School
of Frontier Sciences, The University of
Tokyo, 5-1-5 Kashiwa-no-ha, Kashiwa, Chiba 277-8561, Japan
| | - Kazuhiro Takemura
- School
of Life Sciences and Technology, Tokyo Institute
of Technology, 2-12-1, Ookayama, Meguro-ku, Tokyo 152-8550, Japan
| | - Akio Kitao
- School
of Life Sciences and Technology, Tokyo Institute
of Technology, 2-12-1, Ookayama, Meguro-ku, Tokyo 152-8550, Japan
| | - Kei Terayama
- RIKEN Center for
Advanced Intelligence Project, 1-4-1 Nihonbashi, Chuo-ku, Tokyo 103-0027, Japan
- Medical Sciences
Innovation Hub Program, RIKEN Cluster for Science, Technology and
Innovation Hub, Kanagawa 230-0045, Japan
- Department
of Biomedical Data Intelligence, Graduate School of Medicine, Kyoto University, 54 Shogoin-Kawaharacho, Sakyo-ku, Kyoto 606-8507, Japan
- E-mail: (Kei Terayama)
| | - Koji Tsuda
- Graduate School
of Frontier Sciences, The University of
Tokyo, 5-1-5 Kashiwa-no-ha, Kashiwa, Chiba 277-8561, Japan
- RIKEN Center for
Advanced Intelligence Project, 1-4-1 Nihonbashi, Chuo-ku, Tokyo 103-0027, Japan
- Research
and Services Division of Materials Data and Integrated System, National Institute for Materials Science, Ibaraki 305-0047, Japan
- E-mail: . Phone: +81(4)-7136-3983. Fax: +81(4)-7136-3975 (Koji Tsuda)
| |
Collapse
|
26
|
Peter EK, Černý J. A Hybrid Hamiltonian for the Accelerated Sampling along Experimental Restraints. Int J Mol Sci 2019; 20:E370. [PMID: 30654563 PMCID: PMC6359555 DOI: 10.3390/ijms20020370] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2018] [Revised: 01/07/2019] [Accepted: 01/10/2019] [Indexed: 12/25/2022] Open
Abstract
In this article, we present an enhanced sampling method based on a hybrid Hamiltonian which combines experimental distance restraints with a bias dependent from multiple path-dependent variables. This simulation method determines the bias-coordinates on the fly and does not require a priori knowledge about reaction coordinates. The hybrid Hamiltonian accelerates the sampling of proteins, and, combined with experimental distance information, the technique considers the restraints adaptively and in dependency of the system's intrinsic dynamics. We validate the methodology on the dipole relaxation of two water models and the conformational landscape of dialanine. Using experimental NMR-restraint data, we explore the folding landscape of the TrpCage mini-protein and in a second example apply distance restraints from chemical crosslinking/mass spectrometry experiments for the sampling of the conformation space of the Killer Cell Lectin-like Receptor Subfamily B Member 1A (NKR-P1A). The new methodology has the potential to adaptively introduce experimental restraints without affecting the conformational space of the system along an ergodic trajectory. Since only a limited number of input- and no-order parameters are required for the setup of the simulation, the method is broadly applicable and has the potential to be combined with coarse-graining methods.
Collapse
Affiliation(s)
- Emanuel K Peter
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Průmyslová 595, 252 50 Vestec, Prague West, Czech Republic.
| | - Jiří Černý
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Průmyslová 595, 252 50 Vestec, Prague West, Czech Republic.
| |
Collapse
|
27
|
Peter EK, Cerny J. Enriched Conformational Sampling of DNA and Proteins with a Hybrid Hamiltonian Derived from the Protein Data Bank. Int J Mol Sci 2018; 19:E3405. [PMID: 30380800 PMCID: PMC6274895 DOI: 10.3390/ijms19113405] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2018] [Revised: 10/22/2018] [Accepted: 10/27/2018] [Indexed: 12/26/2022] Open
Abstract
In this article, we present a method for the enhanced molecular dynamics simulation of protein and DNA systems called potential of mean force (PMF)-enriched sampling. The method uses partitions derived from the potentials of mean force, which we determined from DNA and protein structures in the Protein Data Bank (PDB). We define a partition function from a set of PDB-derived PMFs, which efficiently compensates for the error introduced by the assumption of a homogeneous partition function from the PDB datasets. The bias based on the PDB-derived partitions is added in the form of a hybrid Hamiltonian using a renormalization method, which adds the PMF-enriched gradient to the system depending on a linear weighting factor and the underlying force field. We validated the method using simulations of dialanine, the folding of TrpCage, and the conformational sampling of the Dickerson⁻Drew DNA dodecamer. Our results show the potential for the PMF-enriched simulation technique to enrich the conformational space of biomolecules along their order parameters, while we also observe a considerable speed increase in the sampling by factors ranging from 13.1 to 82. The novel method can effectively be combined with enhanced sampling or coarse-graining methods to enrich conformational sampling with a partition derived from the PDB.
Collapse
Affiliation(s)
- Emanuel K Peter
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Průmyslová 595, 252 50 Vestec, Czech Republic.
| | - Jiri Cerny
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Průmyslová 595, 252 50 Vestec, Czech Republic.
| |
Collapse
|
28
|
Chen W, Ferguson AL. Molecular enhanced sampling with autoencoders: On-the-fly collective variable discovery and accelerated free energy landscape exploration. J Comput Chem 2018; 39:2079-2102. [PMID: 30368832 DOI: 10.1002/jcc.25520] [Citation(s) in RCA: 113] [Impact Index Per Article: 18.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2017] [Accepted: 06/14/2018] [Indexed: 01/08/2023]
Abstract
Macromolecular and biomolecular folding landscapes typically contain high free energy barriers that impede efficient sampling of configurational space by standard molecular dynamics simulation. Biased sampling can artificially drive the simulation along prespecified collective variables (CVs), but success depends critically on the availability of good CVs associated with the important collective dynamical motions. Nonlinear machine learning techniques can identify such CVs but typically do not furnish an explicit relationship with the atomic coordinates necessary to perform biased sampling. In this work, we employ auto-associative artificial neural networks ("autoencoders") to learn nonlinear CVs that are explicit and differentiable functions of the atomic coordinates. Our approach offers substantial speedups in exploration of configurational space, and is distinguished from existing approaches by its capacity to simultaneously discover and directly accelerate along data-driven CVs. We demonstrate the approach in simulations of alanine dipeptide and Trp-cage, and have developed an open-source and freely available implementation within OpenMM. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Wei Chen
- Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois, 61801
| | - Andrew L Ferguson
- Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois, 61801.,Department of Materials Science and Engineering, University of Illinois at Urbana-Champaign, 1304 W Green Street, Urbana, Illinois, 61801.,Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois, 61801
| |
Collapse
|
29
|
Borgohain G, Paul S. Atomistic level understanding of the stabilization of protein Trp cage in denaturing and mixed osmolyte solutions. COMPUT THEOR CHEM 2018. [DOI: 10.1016/j.comptc.2018.03.013] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
30
|
Salamatova E, Cunha AV, Bloem R, Roeters SJ, Woutersen S, Jansen TLC, Pshenichnikov MS. Hydrophobic Collapse in N-Methylacetamide-Water Mixtures. J Phys Chem A 2018; 122:2468-2478. [PMID: 29425450 PMCID: PMC6028151 DOI: 10.1021/acs.jpca.8b00276] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2018] [Revised: 02/02/2018] [Indexed: 11/28/2022]
Abstract
Aqueous N-methylacetamide solutions were investigated by polarization-resolved pump-probe and 2D infrared spectroscopy (2D IR), using the amide I mode as a reporter. The 2D IR results are compared with molecular dynamics simulations and spectral calculations to gain insight into the molecular structures in the mixture. N-Methylacetamide and water molecules tend to form clusters with "frozen" amide I dynamics. This is driven by a hydrophobic collapse as the methyl groups of the N-methylacetamide molecules cluster in the presence of water. Since the studied system can be considered as a simplified model for the backbone of proteins, the present study forms a convenient basis for understanding the structural and vibrational dynamics in proteins. It is particularly interesting to find out that a hydrophobic collapse as the one driving protein folding is observed in such a simple system.
Collapse
Affiliation(s)
- Evgeniia Salamatova
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Ana V. Cunha
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Robbert Bloem
- Van
’t Hoff Institute for Molecular Sciences, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Steven J. Roeters
- Van
’t Hoff Institute for Molecular Sciences, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Sander Woutersen
- Van
’t Hoff Institute for Molecular Sciences, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Thomas L. C. Jansen
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Maxim S. Pshenichnikov
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
31
|
Zavadlav J, Marrink SJ, Praprotnik M. Multiscale Simulation of Protein Hydration Using the SWINGER Dynamical Clustering Algorithm. J Chem Theory Comput 2018; 14:1754-1761. [PMID: 29439560 DOI: 10.1021/acs.jctc.7b01129] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
To perform computationally efficient concurrent multiscale simulations of biological macromolecules in solution, where the all-atom (AT) models are coupled to supramolecular coarse-grained (SCG) solvent models, previous studies resorted to modified AT water models, such as the bundled-simple point charge (SPC) models, that use semiharmonic springs to restrict the relative movement of water molecules within a cluster. Those models can have a significant impact on the simulated biomolecules and can lead, for example, to a partial unfolding of a protein. In this work, we employ the recently developed alternative approach with a dynamical clustering algorithm, SWINGER, which enables a direct coupling of original unmodified AT and SCG water models. We perform an adaptive resolution molecular dynamics simulation of a Trp-Cage miniprotein in multiscale water, where the standard SPC water model is interfaced with the widely used MARTINI SCG model, and demonstrate that, compared to the corresponding full-blown AT simulations, the structural and dynamic properties of the solvated protein and surrounding solvent are well reproduced by our approach.
Collapse
Affiliation(s)
- Julija Zavadlav
- Computational Science & Engineering Laboratory , ETH Zurich , Clausiusstrasse 33 , CH-8092 Zurich , Switzerland
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials , University of Groningen , Nijenborgh 7 , 9747 AG Groningen , The Netherlands
| | - Matej Praprotnik
- Laboratory for Molecular Modeling , National Institute of Chemistry , Hajdrihova 19 , SI-1001 Ljubljana , Slovenia.,Department of Physics, Faculty of Mathematics and Physics , University of Ljubljana , Jadranska 19 , SI-1000 Ljubljana , Slovenia
| |
Collapse
|
32
|
Peter EK, Shea JE. An adaptive bias - hybrid MD/kMC algorithm for protein folding and aggregation. Phys Chem Chem Phys 2018. [PMID: 28650060 DOI: 10.1039/c7cp03035e] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
In this paper, we present a novel hybrid Molecular Dynamics/kinetic Monte Carlo (MD/kMC) algorithm and apply it to protein folding and aggregation in explicit solvent. The new algorithm uses a dynamical definition of biases throughout the MD component of the simulation, normalized in relation to the unbiased forces. The algorithm guarantees sampling of the underlying ensemble in dependency of one average linear coupling factor 〈α〉τ. We test the validity of the kinetics in simulations of dialanine and compare dihedral transition kinetics with long-time MD-simulations. We find that for low 〈α〉τ values, kinetics are in good quantitative agreement. In folding simulations of TrpCage and TrpZip4 in explicit solvent, we also find good quantitative agreement with experimental results and prior MD/kMC simulations. Finally, we apply our algorithm to study growth of the Alzheimer Amyloid Aβ 16-22 fibril by monomer addition. We observe two possible binding modes, one at the extremity of the fibril (elongation) and one on the surface of the fibril (lateral growth), on timescales ranging from ns to 8 μs.
Collapse
Affiliation(s)
- Emanuel K Peter
- Department of Pharmacy and Chemistry, Institute of Physical and Theoretical Chemistry, University of Regensburg, Germany
| | | |
Collapse
|
33
|
Chen G, Miao M, Jiang B, Jin J, Campanella OH, Feng B. Effects of high hydrostatic pressure on Rhizopus chinensis lipase: II. Intermediate states during unfolding. INNOV FOOD SCI EMERG 2018. [DOI: 10.1016/j.ifset.2017.08.018] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
|
34
|
Goyal S, Chattopadhyay A, Kasavajhala K, Priyakumar UD. Role of Urea–Aromatic Stacking Interactions in Stabilizing the Aromatic Residues of the Protein in Urea-Induced Denatured State. J Am Chem Soc 2017; 139:14931-14946. [DOI: 10.1021/jacs.7b05463] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Siddharth Goyal
- Center for Computational
Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| | - Aditya Chattopadhyay
- Center for Computational
Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| | - Koushik Kasavajhala
- Center for Computational
Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| | - U. Deva Priyakumar
- Center for Computational
Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| |
Collapse
|
35
|
Andryushchenko VA, Chekmarev SF. Temperature evolution of Trp-cage folding pathways: An analysis by dividing the probability flux field into stream tubes. J Biol Phys 2017; 43:565-583. [PMID: 28983809 DOI: 10.1007/s10867-017-9470-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2016] [Accepted: 09/01/2017] [Indexed: 11/25/2022] Open
Abstract
Owing to its small size and very fast folding rate, the Trp-cage miniprotein has become a benchmark system to study protein folding. Two folding pathways were found to be characteristic of this protein: pathway I, in which the hydrophobic collapse precedes the formation of α-helix, and pathway II, in which the events occur in the reverse order. At the same time, the relative contribution of these pathways at different temperatures as well as the nature of transition from one pathway to the other remain unclear. To gain insight into this issue, we employ a recently proposed hydrodynamic description of protein folding, in which the process of folding is considered as a motion of a "folding fluid" (Chekmarev et al., Phys. Rev. Lett. 100(1), 018107 2008). Using molecular dynamics simulations, we determine the field of probability fluxes of transitions in a space of collective variables and divide it into stream tubes. Each tube contains a definite fraction of the total folding flow and can be associated with a certain pathway. Specifically, three temperatures were considered, T = 285K, T = 315K, and T = 325K. We have found that as the temperature increases, the contribution of pathway I, which is approximately 90% of the total folding flow at T = 285K, decreases to approximately 10% at T = 325K, i.e., pathway II becomes dominant. At T = 315K, both pathways contribute approximately equally. All these temperatures are found below the calculated melting point, which suggests that the Trp-cage folding mechanism is determined by kinetic factors rather than thermodynamics.
Collapse
Affiliation(s)
- Vladimir A Andryushchenko
- Institute of Thermophysics, SB RAS, 630090, Novosibirsk, Russia
- Department of Physics, Novosibirsk State University, 630090, Novosibirsk, Russia
| | - Sergei F Chekmarev
- Institute of Thermophysics, SB RAS, 630090, Novosibirsk, Russia.
- Department of Physics, Novosibirsk State University, 630090, Novosibirsk, Russia.
| |
Collapse
|
36
|
Koole C, Reynolds CA, Mobarec JC, Hick C, Sexton PM, Sakmar TP. Genetically encoded photocross-linkers determine the biological binding site of exendin-4 peptide in the N-terminal domain of the intact human glucagon-like peptide-1 receptor (GLP-1R). J Biol Chem 2017; 292:7131-7144. [PMID: 28283573 PMCID: PMC5409479 DOI: 10.1074/jbc.m117.779496] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2017] [Revised: 03/09/2017] [Indexed: 12/25/2022] Open
Abstract
The glucagon-like peptide-1 receptor (GLP-1R) is a key therapeutic target in the management of type II diabetes mellitus, with actions including regulation of insulin biosynthesis and secretion, promotion of satiety, and preservation of β-cell mass. Like most class B G protein-coupled receptors (GPCRs), there is limited knowledge linking biological activity of the GLP-1R with the molecular structure of an intact, full-length, and functional receptor·ligand complex. In this study, we have utilized genetic code expansion to site-specifically incorporate the photoactive amino acid p-azido-l-phenylalanine (azF) into N-terminal residues of a full-length functional human GLP-1R in mammalian cells. UV-mediated photolysis of azF was then carried out to induce targeted photocross-linking to determine the proximity of the azido group in the mutant receptor with the peptide exendin-4. Cross-linking data were compared directly with the crystal structure of the isolated N-terminal extracellular domain of the GLP-1R in complex with exendin(9-39), revealing both similarities as well as distinct differences in the mode of interaction. Generation of a molecular model to accommodate the photocross-linking constraints highlights the potential influence of environmental conditions on the conformation of the receptor·peptide complex, including folding dynamics of the peptide and formation of dimeric and higher order oligomeric receptor multimers. These data demonstrate that crystal structures of isolated receptor regions may not give a complete reflection of peptide/receptor interactions and should be combined with additional experimental constraints to reveal peptide/receptor interactions occurring in the dynamic, native, and full-length receptor state.
Collapse
Affiliation(s)
- Cassandra Koole
- From the Laboratory of Chemical Biology and Signal Transduction, The Rockefeller University, New York, New York 10065
| | - Christopher A Reynolds
- the School of Biological Sciences, University of Essex, Wivenhoe Park, Colchester CO4 3SQ, United Kingdom
| | - Juan C Mobarec
- the School of Biological Sciences, University of Essex, Wivenhoe Park, Colchester CO4 3SQ, United Kingdom
| | - Caroline Hick
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria 3052, Australia, and
| | - Patrick M Sexton
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Victoria 3052, Australia, and
| | - Thomas P Sakmar
- From the Laboratory of Chemical Biology and Signal Transduction, The Rockefeller University, New York, New York 10065,
| |
Collapse
|
37
|
Meshkin H, Zhu F. Thermodynamics of Protein Folding Studied by Umbrella Sampling along a Reaction Coordinate of Native Contacts. J Chem Theory Comput 2017; 13:2086-2097. [PMID: 28355066 DOI: 10.1021/acs.jctc.6b01171] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Spontaneous transitions between the native and non-native protein conformations are normally rare events that hardly take place in typical unbiased molecular dynamics simulations. It was recently demonstrated that such transitions can be well described by a reaction coordinate, Q, that represents the collective fraction of the native contacts between the protein atoms. Here we attempt to use this reaction coordinate to enhance the conformational sampling. We perform umbrella sampling simulations with biasing potentials on Q for two model proteins, Trp-Cage and BBA, using the CHARMM force field. Hamiltonian replica exchange is implemented in these simulations to further facilitate the sampling. The simulations appear to have reached satisfactory convergence, resulting in unbiased free energies as a function of Q. In addition to the native structure, multiple folded conformations are identified in the reconstructed equilibrium ensemble. Some conformations without any native contacts nonetheless have rather compact geometries and are stabilized by hydrogen bonds not present in the native structure. Whereas the enhanced sampling along Q reasonably reproduces the equilibrium conformational space, we also find that the folding of an α-helix in Trp-Cage is a slow degree of freedom orthogonal to Q and therefore cannot be accelerated by biasing the reaction coordinate. Overall, we conclude that whereas Q is an excellent parameter to analyze the simulations, it is not necessarily a perfect reaction coordinate for enhanced sampling, and better incorporation of other slow degrees of freedom may further improve this reaction coordinate.
Collapse
Affiliation(s)
- Hamed Meshkin
- Department of Physics, Indiana University Purdue University Indianapolis , 402 North Blackford Street, Indianapolis, Indiana 46202, United States
| | - Fangqiang Zhu
- Department of Physics, Indiana University Purdue University Indianapolis , 402 North Blackford Street, Indianapolis, Indiana 46202, United States
| |
Collapse
|
38
|
Ota M, Ikeguchi M, Kidera A. Itinerary profiling to analyze a large number of protein-folding trajectories. Biophys Physicobiol 2016; 13:295-304. [PMID: 28409081 PMCID: PMC5221515 DOI: 10.2142/biophysico.13.0_295] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2016] [Accepted: 09/05/2016] [Indexed: 12/01/2022] Open
Abstract
Understanding how proteins fold through a vast number of unfolded states is a major subject in the study of protein folding. Herein, we present itinerary profiling as a simple method to analyze molecular dynamics trajectories, and apply this method to Trp-cage. In itinerary profiling, structural clusters included in a trajectory are represented by a bit sequence, and a number of trajectories, as well as the structural clusters, can be compared and classified. As a consequence, the structural clusters that characterize the foldability of trajectories were able to be identified. The connections between the clusters were then illustrated as a network and the structural features of the clusters were examined. We found that in the true folding funnel, Trp-cage formed a left-handed main-chain topology and the Trp6 side-chain was located at the front of the main-chain ring, even in the initial unfolded states. In contrast, in the false folding funnel of the pseudo-native states, in which the Trp6 side-chain is upside down in the protein core, Trp-cage had a right-handed main-chain topology and the Trp side-chain was at the back. The initial topological partition, as determined by the main-chain handedness and the location of the Trp residue, predetermines Trp-cage foldability and the destination of the trajectory to the native state or the pseudo-native states.
Collapse
Affiliation(s)
- Motonori Ota
- Graduate School of Information Science, Nagoya University, Nagoya, Aichi 464-8601, Japan
| | - Mitsunori Ikeguchi
- Graduate School of Medical Life Science, Yokohama City University, Yokohama, Kanagawa 230-0045, Japan
| | - Akinori Kidera
- Graduate School of Medical Life Science, Yokohama City University, Yokohama, Kanagawa 230-0045, Japan
| |
Collapse
|
39
|
Meuzelaar H, Panman MR, van Dijk CN, Woutersen S. Folding of a Zinc-Finger ββα-Motif Investigated Using Two-Dimensional and Time-Resolved Vibrational Spectroscopy. J Phys Chem B 2016; 120:11151-11158. [PMID: 27723346 DOI: 10.1021/acs.jpcb.6b08883] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Small proteins provide good model systems for studying the fundamental forces that control protein folding. Here, we investigate the folding dynamics of the 28-residue zinc-finger mutant FSD-1, which is designed to form a metal-independent folded ββα-motif, and which provides a testing ground for proteins containing a mixed α/β fold. Although the folding of FSD-1 has been actively studied, the folding mechanism remains largely unclear. In particular, it is unclear in what stage of folding the α-helix is formed. To address this issue we investigate the folding mechanism of FSD-1 using a combination of temperature-dependent UV circular dichroism (UV-CD), Fourier transform infrared (FTIR) spectroscopy, two-dimensional infrared (2D-IR) spectroscopy, and temperature-jump (T-jump) transient-IR spectroscopy. Our UV-CD and FTIR data show different thermal melting transitions, indicating multistate folding behavior. Temperature-dependent 2D-IR spectra indicate that the α-helix is the most stable structural element of FSD-1. To investigate the folding/unfolding re-equilibration dynamics of FSD-1, the conformational changes induced by a nanosecond T-jump are probed with transient-IR and transient dispersed-pump-probe (DPP) IR spectroscopy. We observe biexponential T-jump relaxation kinetics (with time constants of 80 ± 13 ns and 1300 ± 100 ns at 322 K), confirming that the folding involves an intermediate state. The IR and dispersed-pump-probe IR spectra associated with the two kinetic components suggest that the folding of FSD-1 involves early formation of the α-helix, followed by the formation of the β-hairpin and hydrophobic contacts.
Collapse
Affiliation(s)
- Heleen Meuzelaar
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Matthijs R Panman
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Chris N van Dijk
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Sander Woutersen
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098 XH Amsterdam, The Netherlands
| |
Collapse
|
40
|
Bondarenko AS, Jansen TLC. Application of two-dimensional infrared spectroscopy to benchmark models for the amide I band of proteins. J Chem Phys 2016; 142:212437. [PMID: 26049457 DOI: 10.1063/1.4919716] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
In this paper, we present a novel benchmarking method for validating the modelling of vibrational spectra for the amide I region of proteins. We use the linear absorption spectra and two-dimensional infrared spectra of four experimentally well-studied proteins as a reference and test nine combinations of molecular dynamics force fields, vibrational frequency mappings, and coupling models. We find that two-dimensional infrared spectra provide a much stronger test of the models than linear absorption does. The best modelling approach in the present study still leaves significant room for future improvement. The presented benchmarking scheme, thus, provides a way of validating future protocols for modelling the amide I band in proteins.
Collapse
Affiliation(s)
- Anna S Bondarenko
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Thomas L C Jansen
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
41
|
Zhang H, Lv Y, Tan T, van der Spoel D. Atomistic Simulation of Protein Encapsulation in Metal–Organic Frameworks. J Phys Chem B 2016; 120:477-84. [DOI: 10.1021/acs.jpcb.5b10437] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Haiyang Zhang
- Beijing
Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yongqin Lv
- Beijing
Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
| | - Tianwei Tan
- Beijing
Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
| | - David van der Spoel
- Uppsala
Center for Computational Chemistry, Science for Life Laboratory, Department
of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box
596, SE-75124 Uppsala, Sweden
| |
Collapse
|
42
|
Liang C. Catching protein structural dynamics by two-dimensional infrared spectroscopy. Biophys J 2016; 108:1577-1579. [PMID: 25863046 DOI: 10.1016/j.bpj.2015.03.007] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2015] [Revised: 02/26/2015] [Accepted: 03/04/2015] [Indexed: 02/01/2023] Open
Affiliation(s)
- Chungwen Liang
- Department of Chemistry and Biochemistry, University of California, San Diego, San Diego, California.
| |
Collapse
|
43
|
Peter EK, Shea JE, Pivkin IV. Coarse kMC-based replica exchange algorithms for the accelerated simulation of protein folding in explicit solvent. Phys Chem Chem Phys 2016; 18:13052-65. [DOI: 10.1039/c5cp06867c] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
In this paper, we present a coarse replica exchange molecular dynamics (REMD) approach, based on kinetic Monte Carlo (kMC).
Collapse
Affiliation(s)
- Emanuel K. Peter
- Institute of Computational Science
- Faculty of Informatics
- University of Lugano
- Switzerland
| | - Joan-Emma Shea
- Department of Chemistry and Biochemistry
- Department of Physics
- University of California
- Santa Barbara
- USA
| | - Igor V. Pivkin
- Institute of Computational Science
- Faculty of Informatics
- University of Lugano
- Switzerland
- Swiss Institute of Bioinformatics
| |
Collapse
|
44
|
Harada R, Nakamura T, Shigeta Y. Sparsity-weighted outlier FLOODing (OFLOOD) method: Efficient rare event sampling method using sparsity of distribution. J Comput Chem 2015; 37:724-38. [DOI: 10.1002/jcc.24255] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2015] [Revised: 10/06/2015] [Accepted: 10/29/2015] [Indexed: 01/18/2023]
Affiliation(s)
- Ryuhei Harada
- Department of Physics Graduate School of Pure and Applied Sciences; University of Tsukuba; 1-1-1 Tennodai Tsukuba Ibaraki 305-8571 Japan
- Division of Life Sciences; Center for Computational Sciences, University of Tsukuba; 1-1-1 Tennodai Tsukuba Ibaraki 305-8571 Japan
- Computational Engineering Application Unit, RIKEN Advanced Center for Computing and Communication; 2-1, Hirosawa Wako Saitama 351-0198 Japan
| | | | - Yasuteru Shigeta
- Department of Physics Graduate School of Pure and Applied Sciences; University of Tsukuba; 1-1-1 Tennodai Tsukuba Ibaraki 305-8571 Japan
- Division of Life Sciences; Center for Computational Sciences, University of Tsukuba; 1-1-1 Tennodai Tsukuba Ibaraki 305-8571 Japan
- JST, CREST; 4-1-8 Honcho Kawaguchi Saitama 332-0012 Japan
| |
Collapse
|
45
|
A hydrodynamic view of the first-passage folding of Trp-cage miniprotein. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2015; 45:229-43. [PMID: 26559408 DOI: 10.1007/s00249-015-1089-7] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/24/2015] [Revised: 09/27/2015] [Accepted: 10/09/2015] [Indexed: 12/11/2022]
Abstract
We study folding of Trp-cage miniprotein in the conditions when the native state of the protein is stable and unfolding events are improbable, which corresponds to physiological conditions. Using molecular dynamics simulations with an implicit solvent model, an ensemble of folding trajectories from unfolded (practically extended) states of the protein to the native state was generated. To get insight into the folding kinetics, the free energy surface and kinetic network projected on this surface were constructed. This, "conventional" analysis of the folding reaction was followed by a recently proposed hydrodynamic description of protein folding (Chekmarev et al. in Phys Rev Lett 100(1):018107, 2008), in which the process of the first-passage folding is viewed as a stationary flow of a folding "fluid" from the unfolded to native state. This approach is conceptually different from the previously used approaches and thus allows an alternative view of the folding dynamics and kinetics of Trp-cage, the conclusions about which are very diverse. In agreement with most previous studies, we observed two characteristic folding pathways: in one pathway (I), the collapse of the hydrophobic core precedes the formation of the [Formula: see text]-helix, and in the other pathway (II), these events occur in the reverse order. We found that although pathway II is complicated by a repeated partial protein unfolding, it contributes to the total folding flow as little as ≈10%, so that the folding kinetics remain essentially single-exponential.
Collapse
|
46
|
Zhou CY, Jiang F, Wu YD. Folding Thermodynamics and Mechanism of Five Trp-Cage Variants from Replica-Exchange MD Simulations with RSFF2 Force Field. J Chem Theory Comput 2015; 11:5473-80. [PMID: 26574335 DOI: 10.1021/acs.jctc.5b00581] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
To test whether our recently developed residue-specific force field RSFF2 can reproduce the mutational effect on the thermal stability of Trp-cage mini-protein and decipher its detailed folding mechanism, we carried out long-time replica-exchange molecular dynamics (REMD) simulations on five Trp-cage variants, including TC5b and TC10b. Initiated from their unfolded structures, the simulations not only well-reproduce their experimental structures but also their melting temperatures and folding enthalpies reasonably well. For each Trp-cage variant, the overall folding free energy landscape is apparently two-state, but some intermediate states can be observed when projected on more detailed coordinates. We also found different variants have the same major folding pathway, including the well formed PII-helix in the unfolded state, the formation of W6-P12/P18/P19 contacts and the α-helix before the transition state, the following formation of most native contacts, and the final native loop formation. The folding mechanism derived here is consistent with many previous simulations and experiments.
Collapse
Affiliation(s)
- Chen-Yang Zhou
- Laboratory of Computational Chemistry and Drug Design, Laboratory of Chemical Genomics, Peking University Shenzhen Graduate School , Shenzhen 518055, China.,College of Chemistry and Molecular Engineering, Peking University , Beijing 100871, China
| | - Fan Jiang
- Laboratory of Computational Chemistry and Drug Design, Laboratory of Chemical Genomics, Peking University Shenzhen Graduate School , Shenzhen 518055, China
| | - Yun-Dong Wu
- Laboratory of Computational Chemistry and Drug Design, Laboratory of Chemical Genomics, Peking University Shenzhen Graduate School , Shenzhen 518055, China.,College of Chemistry and Molecular Engineering, Peking University , Beijing 100871, China
| |
Collapse
|
47
|
Baker JL, Furbish J, Lindberg GE. Influence of the ionic liquid [C4mpy][Tf2N] on the structure of the miniprotein Trp-cage. J Mol Graph Model 2015; 62:202-212. [PMID: 26479192 DOI: 10.1016/j.jmgm.2015.10.003] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2015] [Revised: 09/14/2015] [Accepted: 10/06/2015] [Indexed: 10/22/2022]
Abstract
We examine the effect of the ionic liquid [C4mpy][Tf2N] on the structure of the miniprotein Trp-cage and contrast these results with the behavior of Trp-cage in water. We find the ionic liquid has a dramatic effect on Trp-cage, though many similarities with aqueous Trp-cage are observed. We assess Trp-cage folding by monitoring root mean square deviation from the crystallographic structure, radius of gyration, proline cis/trans isomerization state, protein secondary structure, amino acid contact formation and distance, and native and non-native contact formation. Starting from an unfolded configuration, Trp-cage folds in water at 298 K in less than 500 ns of simulation, but has very little mobility in the ionic liquid at the same temperature, which can be ascribed to the higher ionic liquid viscosity. At 365 K, the mobility of the ionic liquid is increased and initial stages of Trp-cage folding are observed, however Trp-cage does not reach the native folded state in 2 μs of simulation in the ionic liquid. Therefore, in addition to conventional molecular dynamics, we also employ scaled molecular dynamics to expedite sampling, and we demonstrate that Trp-cage in the ionic liquid does closely approach the aqueous folded state. Interestingly, while the reduced mobility of the ionic liquid is found to restrict Trp-cage motion, the ionic liquid does facilitate proline cis/trans isomerization events that are not seen in our aqueous simulations.
Collapse
Affiliation(s)
- Joseph L Baker
- Department of Chemistry, The College of New Jersey, 2000 Pennington Road, Ewing, NJ 08628, United States.
| | - Jeffrey Furbish
- Department of Chemistry, The College of New Jersey, 2000 Pennington Road, Ewing, NJ 08628, United States
| | - Gerrick E Lindberg
- Department of Chemistry and Biochemistry, Northern Arizona University, 700 S. Osbourne Drive, Flagstaff, AZ 86011, United States.
| |
Collapse
|
48
|
Zang T, Yu L, Zhang C, Ma J. Parallel continuous simulated tempering and its applications in large-scale molecular simulations. J Chem Phys 2015; 141:044113. [PMID: 25084887 DOI: 10.1063/1.4890038] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
In this paper, we introduce a parallel continuous simulated tempering (PCST) method for enhanced sampling in studying large complex systems. It mainly inherits the continuous simulated tempering (CST) method in our previous studies [C. Zhang and J. Ma, J. Chem. Phys. 130, 194112 (2009); C. Zhang and J. Ma, J. Chem. Phys. 132, 244101 (2010)], while adopts the spirit of parallel tempering (PT), or replica exchange method, by employing multiple copies with different temperature distributions. Differing from conventional PT methods, despite the large stride of total temperature range, the PCST method requires very few copies of simulations, typically 2-3 copies, yet it is still capable of maintaining a high rate of exchange between neighboring copies. Furthermore, in PCST method, the size of the system does not dramatically affect the number of copy needed because the exchange rate is independent of total potential energy, thus providing an enormous advantage over conventional PT methods in studying very large systems. The sampling efficiency of PCST was tested in two-dimensional Ising model, Lennard-Jones liquid and all-atom folding simulation of a small globular protein trp-cage in explicit solvent. The results demonstrate that the PCST method significantly improves sampling efficiency compared with other methods and it is particularly effective in simulating systems with long relaxation time or correlation time. We expect the PCST method to be a good alternative to parallel tempering methods in simulating large systems such as phase transition and dynamics of macromolecules in explicit solvent.
Collapse
Affiliation(s)
- Tianwu Zang
- Applied Physics Program and Department of Bioengineering, Rice University, Houston, Texas 77005, USA
| | - Linglin Yu
- Applied Physics Program and Department of Bioengineering, Rice University, Houston, Texas 77005, USA
| | - Chong Zhang
- Applied Physics Program and Department of Bioengineering, Rice University, Houston, Texas 77005, USA
| | - Jianpeng Ma
- Applied Physics Program and Department of Bioengineering, Rice University, Houston, Texas 77005, USA
| |
Collapse
|
49
|
Kamo F, Ishizuka R, Matubayasi N. Correlation analysis for heat denaturation of Trp-cage miniprotein with explicit solvent. Protein Sci 2015; 25:56-66. [PMID: 26189564 DOI: 10.1002/pro.2754] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2015] [Accepted: 07/15/2015] [Indexed: 11/10/2022]
Abstract
Energetics was analyzed for Trp-cage miniprotein in water to elucidate the solvation effect in heat denaturation. The solvation free energy was computed for a set of protein structures at room and high temperatures with all-atom molecular dynamics simulation combined with the solution theory in the energy representation, and its correlations were investigated against the intramolecular (structural) energy of the protein and the average interaction energy of the protein with the solvent water. It was observed both at room and high temperatures that the solvation free energy is anticorrelated to the structural energy and varies in parallel to the electrostatic component of the protein-water interaction energy without correlations to the van der Waals and excluded-volume components. When the set of folded structures sampled at room temperature was compared with the set of unfolded ones at high temperature, it was found that the preference order of the two sets is in correspondence to the van der Waals and excluded-volume components in the sum of the protein intramolecular and protein-water intermolecular interactions and is not distinguished by the electrostatic component.
Collapse
Affiliation(s)
- Fumitaka Kamo
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka, 560-8531, Japan
| | - Ryosuke Ishizuka
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka, 560-8531, Japan.,Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Katsura, Kyoto, 615-8520, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka, 560-8531, Japan.,Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Katsura, Kyoto, 615-8520, Japan
| |
Collapse
|
50
|
Miao Y, Feixas F, Eun C, McCammon JA. Accelerated molecular dynamics simulations of protein folding. J Comput Chem 2015; 36:1536-49. [PMID: 26096263 DOI: 10.1002/jcc.23964] [Citation(s) in RCA: 110] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2015] [Revised: 05/11/2015] [Accepted: 05/19/2015] [Indexed: 02/02/2023]
Abstract
Folding of four fast-folding proteins, including chignolin, Trp-cage, villin headpiece and WW domain, was simulated via accelerated molecular dynamics (aMD). In comparison with hundred-of-microsecond timescale conventional molecular dynamics (cMD) simulations performed on the Anton supercomputer, aMD captured complete folding of the four proteins in significantly shorter simulation time. The folded protein conformations were found within 0.2-2.1 Å of the native NMR or X-ray crystal structures. Free energy profiles calculated through improved reweighting of the aMD simulations using cumulant expansion to the second-order are in good agreement with those obtained from cMD simulations. This allows us to identify distinct conformational states (e.g., unfolded and intermediate) other than the native structure and the protein folding energy barriers. Detailed analysis of protein secondary structures and local key residue interactions provided important insights into the protein folding pathways. Furthermore, the selections of force fields and aMD simulation parameters are discussed in detail. Our work shows usefulness and accuracy of aMD in studying protein folding, providing basic references in using aMD in future protein-folding studies.
Collapse
Affiliation(s)
- Yinglong Miao
- Howard Hughes Medical Institute, University of California at San Diego, La Jolla, California
| | - Ferran Feixas
- Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California.,Department of Pharmacology, University of California at San Diego, La Jolla, California
| | - Changsun Eun
- Howard Hughes Medical Institute, University of California at San Diego, La Jolla, California
| | - J Andrew McCammon
- Howard Hughes Medical Institute, University of California at San Diego, La Jolla, California.,Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California.,Department of Pharmacology, University of California at San Diego, La Jolla, California
| |
Collapse
|