The native-to-loop (N-L) unfolding transition of Trp-cage protein was studied via optimized forward flux sampling (FFS) methods with trajectories evolved using molecular dynamics. The rate constant calculated from our simulations is in good agreement with the experimental value for the native-to-unfolded transition of this protein; furthermore, the trajectories sampled a phase region consistent with that reported in previous studies for the N-L transition using transition path sampling and transition interface sampling. A new variant of FFS is proposed and implemented that allows a better control of a constant flux of partial paths. A reaction coordinate model was obtained, at no extra cost, from the transition path ensemble generated by FFS, through iterative use of the FFS-least-square estimation method [E. E. Borrero and F. A. Escobedo, J. Chem. Phys. 127, 164101 (2007)] and an adaptive staging optimization algorithm [E. E. Borrero and F. A. Escobedo, J. Chem. Phys. 129, 024115 (2008)]. Finally, we further elucidate the unfolding mechanism by correlating the unfolding progress with changes in the root mean square deviation from the α carbons of the native state, the root mean square deviation from an ideal α-helix, and other structural properties of the protein.
Download full-text PDF |
Source |
|---|---|
| http://dx.doi.org/10.1063/1.3474803 | DOI Listing |
Publication Analysis
Top Keywords
Similar Publications
Estimation of vibrational spectra of Trp-cage protein from nonequilibrium metadynamics simulations.
Biophys J
October 2024
Department of Chemical & Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina. Electronic address:
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.
View Article and Find Full Text PDFA New Transfer Free Energy Based Implicit Solvation Model for the Description of Disordered and Folded Proteins.
J Phys Chem B
August 2022
Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States.
Most biological events occur on time scales that are difficult to access using conventional all-atom molecular dynamics simulations in explicit solvent. Implicit solvent techniques offer a promising solution to this problem, alleviating the computational cost associated with the simulation of large systems and accelerating the sampling compared to explicit solvent models. The substitution of water molecules by a mean field, however, introduces simplifications that may penalize accuracy and impede the prediction of certain physical properties.
View Article and Find Full Text PDFTopological bio-scaling analysis as a universal measure of protein folding.
R Soc Open Sci
July 2022
Infochemistry Scientific Center (ISC), ITMO University, 9 Lomonosova St., St Petersburg 191002, Russia.
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.
View Article and Find Full Text PDFAccurate Machine Learning Prediction of Protein Circular Dichroism Spectra with Embedded Density Descriptors.
JACS Au
December 2021
Hefei National Laboratory for Physical Sciences at the Microscale, Collaborative Innovation Center of Chemistry for Energy Materials, School of Chemistry and Materials Science, University of Science and Technology of China, Hefei, Anhui 230026, P. R. China.
A data-driven approach to simulate circular dichroism (CD) spectra is appealing for fast protein secondary structure determination, yet the challenge of predicting electric and magnetic transition dipole moments poses a substantial barrier for the goal. To address this problem, we designed a new machine learning (ML) protocol in which ordinary pure geometry-based descriptors are replaced with alternative embedded density descriptors and electric and magnetic transition dipole moments are successfully predicted with an accuracy comparable to first-principle calculation. The ML model is able to not only simulate protein CD spectra nearly 4 orders of magnitude faster than conventional first-principle simulation but also obtain CD spectra in good agreement with experiments.
View Article and Find Full Text PDFResidue Folding Degree-Relationship to Secondary Structure Categories and Use as Collective Variable.
Int J Mol Sci
December 2021
Center for Computational Sciences, University of Tsukuba, Tsukuba 305-8577, Ibaraki, Japan.
Recently, we have shown that the residue folding degree, a network-based measure of folded content in proteins, is able to capture backbone conformational transitions related to the formation of secondary structures in molecular dynamics (MD) simulations. In this work, we focus primarily on developing a collective variable (CV) for MD based on this residue-bound parameter to be able to trace the evolution of secondary structure in segments of the protein. We show that this CV can do just that and that the related energy profiles (potentials of mean force, PMF) and transition barriers are comparable to those found by others for particular events in the folding process of the model mini protein Trp-cage.
View Article and Find Full Text PDFWant AI Summaries of new PubMed Abstracts delivered to your In-box?
Enter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!