1
|
Yuan S, Han X, Zhang J, Xie Z, Fan C, Xiao Y, Gao YQ, Yang YI. Generating High-Precision Force Fields for Molecular Dynamics Simulations to Study Chemical Reaction Mechanisms Using Molecular Configuration Transformer. J Phys Chem A 2024; 128:4378-4390. [PMID: 38759697 DOI: 10.1021/acs.jpca.4c01267] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/19/2024]
Abstract
Theoretical studies on chemical reaction mechanisms have been crucial in organic chemistry. Traditionally, calculating the manually constructed molecular conformations of transition states for chemical reactions using quantum chemical calculations is the most commonly used method. However, this way is heavily dependent on individual experience and chemical intuition. In our previous study, we proposed a research paradigm that used enhanced sampling in molecular dynamics simulations to study chemical reactions. This approach can directly simulate the entire process of a chemical reaction. However, the computational speed limited the use of high-precision potential energy functions for simulations. To address this issue, we presented a scheme for training high-precision force fields for molecular modeling using a previously developed graph-neural-network-based molecular model, molecular configuration transformer. This potential energy function allowed for highly accurate simulations at a low computational cost, leading to more precise calculations of the mechanism of chemical reactions. We applied this approach to study a Claisen rearrangement reaction and a carbonyl insertion reaction catalyzed by manganese.
Collapse
Affiliation(s)
- Sihao Yuan
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
| | - Xu Han
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jun Zhang
- Changping Laboratory, Beijing 102200, China
| | - Zhaoxin Xie
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
| | - Cheng Fan
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
| | - Yunlong Xiao
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Yi Qin Gao
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
- Changping Laboratory, Beijing 102200, China
| | - Yi Isaac Yang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
| |
Collapse
|
2
|
Kong X, Huang Z, Zhang C, Hua Y, Chen Y, Li X. Phenolic compounds in walnut pellicle improve walnut (Juglans regia L.) protein solubility under pH-shifting condition. Food Res Int 2023; 163:112156. [PMID: 36596107 DOI: 10.1016/j.foodres.2022.112156] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Revised: 10/30/2022] [Accepted: 11/15/2022] [Indexed: 11/21/2022]
Abstract
This study focused on the interaction of walnut protein with phenolic extracts of walnut pellicle (PEWP) under alkaline condition, leading to enhancement of protein solubility under neutral condition. First, the change of PEWP under alkaline condition was determined by RP-HPLC and mass spectrometry, and the results showed that most ellagitannins in PEWP could be retained under alkaline condition within 3 h. Interaction between PEWP and walnut protein under pH-shifting condition resulted in the remarkable increase of protein solubility (above 90%) at neutral pH. The results from SDS-PAGE and SEC showed that the improved solubility lied in the formation of large and soluble protein aggregates due to the covalent interaction among walnut protein and polyphenols. A significant change in tertiary structure of protein-phenolic complex was witnessed by fluorescence spectrum and near-UV circular dichroism. Meanwhile, walnut protein-polyphenol interaction led to a slight increase in β-turn while a slight decrease in β-sheet. Combined with amino acid composition, it could be illustrated that the covalent bonding for walnut protein with polyphenol mainly occurred at Lysine residues.
Collapse
Affiliation(s)
- Xiangzhen Kong
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, PR China; School of Food Science and Technology, Jiangnan University, Wuxi, PR China.
| | - Zilin Huang
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, PR China; School of Food Science and Technology, Jiangnan University, Wuxi, PR China
| | - Caimeng Zhang
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, PR China; School of Food Science and Technology, Jiangnan University, Wuxi, PR China
| | - Yufei Hua
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, PR China; School of Food Science and Technology, Jiangnan University, Wuxi, PR China
| | - Yeming Chen
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, PR China; School of Food Science and Technology, Jiangnan University, Wuxi, PR China
| | - Xingfei Li
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, PR China; School of Food Science and Technology, Jiangnan University, Wuxi, PR China
| |
Collapse
|
3
|
Su C, He Z, Li H. Covalent interactions between rabbit myofibrillar proteins and quercetin: A promising approach to enhance protein antioxidant capacity and thermal stability. Lebensm Wiss Technol 2022. [DOI: 10.1016/j.lwt.2022.114132] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
4
|
Zhang K, Wen Q, Li T, Zhang Y, Huang J, Huang Q, Gao L. Effect of covalent conjugation with chlorogenic acid and luteolin on allergenicity and functional properties of wheat gliadin. J Cereal Sci 2022. [DOI: 10.1016/j.jcs.2022.103484] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
5
|
Yang YI, Shao Q, Zhang J, Yang L, Gao YQ. Enhanced sampling in molecular dynamics. J Chem Phys 2019; 151:070902. [PMID: 31438687 DOI: 10.1063/1.5109531] [Citation(s) in RCA: 172] [Impact Index Per Article: 34.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
Although molecular dynamics simulations have become a useful tool in essentially all fields of chemistry, condensed matter physics, materials science, and biology, there is still a large gap between the time scale which can be reached in molecular dynamics simulations and that observed in experiments. To address the problem, many enhanced sampling methods were introduced, which effectively extend the time scale being approached in simulations. In this perspective, we review a variety of enhanced sampling methods. We first discuss collective-variables-based methods including metadynamics and variationally enhanced sampling. Then, collective variable free methods such as parallel tempering and integrated tempering methods are presented. At last, we conclude with a brief introduction of some newly developed combinatory methods. We summarize in this perspective not only the theoretical background and numerical implementation of these methods but also the new challenges and prospects in the field of the enhanced sampling.
Collapse
Affiliation(s)
- Yi Isaac Yang
- Institute of Systems Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
| | - Qiang Shao
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Jun Zhang
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, Berlin 14195, Germany
| | - Lijiang Yang
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Yi Qin Gao
- Institute of Systems Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
| |
Collapse
|
6
|
Yang C, Wang B, Wang J, Xia S, Wu Y. Effect of pyrogallic acid (1,2,3-benzenetriol) polyphenol-protein covalent conjugation reaction degree on structure and antioxidant properties of pumpkin (Cucurbita sp.) seed protein isolate. Lebensm Wiss Technol 2019. [DOI: 10.1016/j.lwt.2019.04.034] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
|
7
|
Shao Q, Zhu W. Effective Conformational Sampling in Explicit Solvent with Gaussian Biased Accelerated Molecular Dynamics. J Chem Theory Comput 2017; 13:4240-4252. [DOI: 10.1021/acs.jctc.7b00242] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Qiang Shao
- Drug
Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi
Road, Shanghai, 201203, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Weiliang Zhu
- Drug
Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi
Road, Shanghai, 201203, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| |
Collapse
|
8
|
Simulation of the T-jump triggered unfolding and thermal unfolding vibrational spectroscopy related to polypeptides conformation fluctuation. Sci China Chem 2017. [DOI: 10.1007/s11426-016-9055-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
9
|
A test of AMBER force fields in predicting the secondary structure of α-helical and β-hairpin peptides. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.04.074] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
|
10
|
Shao Q, Shi J, Zhu W. Determining Protein Folding Pathway and Associated Energetics through Partitioned Integrated-Tempering-Sampling Simulation. J Chem Theory Comput 2017; 13:1229-1243. [DOI: 10.1021/acs.jctc.6b00967] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Qiang Shao
- Drug
Discovery and Design Center, CAS Key Laboratory of Receptor Research,
Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Jiye Shi
- UCB Biopharma
SPRL, Chemin du Foriest, 1420 Braine-l’Alleud, Belgium
| | - Weiliang Zhu
- Drug
Discovery and Design Center, CAS Key Laboratory of Receptor Research,
Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| |
Collapse
|
11
|
Gupta M, Khatua P, Chakravarty C, Bandyopadhyay S. The sensitivity of folding free energy landscapes of trpzips to mutations in the hydrophobic core. Phys Chem Chem Phys 2017; 19:22813-22825. [DOI: 10.1039/c7cp03825a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The sensitivity of the stability of folded states and free energy landscapes to the differences in the hydrophobic content of the core residues has been studied for the set of 16-residue trpzips, namely, Trpzip4, Trpzip5 and Trpzip6.
Collapse
Affiliation(s)
- Madhulika Gupta
- Department of Chemistry
- Indian Institute of Technology-Delhi
- New Delhi 110016
- India
| | - Prabir Khatua
- Molecular Modeling Laboratory
- Department of Chemistry
- Indian Institute of Technology
- Kharagpur 721302
- India
| | | | - Sanjoy Bandyopadhyay
- Molecular Modeling Laboratory
- Department of Chemistry
- Indian Institute of Technology
- Kharagpur 721302
- India
| |
Collapse
|
12
|
Xie L, Chen ZN. Enhanced molecular dynamics simulation of the transformation between α-helix and β-hairpin structures for peptide. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1161249] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Liangxu Xie
- Department of Chemistry, The University of Hong Kong, Hong Kong SAR, China
| | - Zhe-Ning Chen
- Department of Chemistry, The University of Hong Kong, Hong Kong SAR, China
| |
Collapse
|
13
|
Makwana KM, Mahalakshmi R. Nature of aryl-tyrosine interactions contribute to β-hairpin scaffold stability: NMR evidence for alternate ring geometry. Phys Chem Chem Phys 2016; 17:4220-30. [PMID: 25569770 DOI: 10.1039/c4cp04991h] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The specific contribution of the acidic-aromatic β-sheet favouring amino acid tyrosine to the stability of short octapeptide β-hairpin structures is presented here. Solution NMR analysis in near-apolar environments suggests the energetically favourable mode of interaction to be T-shaped face-to-edge (FtE) and that a Trp-Tyr interacting pair is the most stabilizing. Alternate aryl geometries also exist in solution, which readily equilibrate between a preferred π···π conformation to an aromatic-amide conformation, without any change in the backbone structure. While the phenolic ring is readily accommodated at the "edge" of FtE aryl interactions, it exhibits an overall lowered contribution to scaffold stability in the "face" orientation. Such differential tyrosine interactions are key to its dual nature in proteins.
Collapse
Affiliation(s)
- Kamlesh Madhusudan Makwana
- Molecular Biophysics Laboratory, Department of Biological Sciences, Indian Institute of Science Education and Research, Bhopal-462023, India.
| | | |
Collapse
|
14
|
Correct folding of an α-helix and a β-hairpin using a polarized 2D torsional potential. Sci Rep 2015; 5:10359. [PMID: 26039188 PMCID: PMC5380191 DOI: 10.1038/srep10359] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2014] [Accepted: 04/09/2015] [Indexed: 01/12/2023] Open
Abstract
A new modification to the AMBER force field that incorporates the coupled two-dimensional main chain torsion energy has been evaluated for the balanced representation of secondary structures. In this modified AMBER force field (AMBER03(2D)), the main chain torsion energy is represented by 2-dimensional Fourier expansions with parameters fitted to the potential energy surface generated by high-level quantum mechanical calculations of small peptides in solution. Molecular dynamics simulations are performed to study the folding of two model peptides adopting either α-helix or β-hairpin structures. Both peptides are successfully folded into their native structures using an AMBER03(2D) force field with the implementation of a polarization scheme (AMBER03(2D)p). For comparison, simulations using a standard AMBER03 force field with and without polarization, as well as AMBER03(2D) without polarization, fail to fold both peptides successfully. The correction to secondary structure propensity in the AMBER03 force field and the polarization effect are critical to folding Trpzip2; without these factors, a helical structure is obtained. This study strongly suggests that this new force field is capable of providing a more balanced preference for helical and extended conformations. The electrostatic polarization effect is shown to be indispensable to the growth of secondary structures.
Collapse
|
15
|
Yang L, Liu CW, Shao Q, Zhang J, Gao YQ. From thermodynamics to kinetics: enhanced sampling of rare events. Acc Chem Res 2015; 48:947-55. [PMID: 25781363 DOI: 10.1021/ar500267n] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
Despite great advances in molecular dynamics simulations, there remain large gaps between the simulations and experimental observations in terms of the time and length scales that can be approached. Developing fast and accurate algorithms and methods is of ultimate importance to bridge these gaps. In this Account, we briefly summarize recent efforts in such directions. In particular, we focus on integrated tempering sampling. The efficiency of this sampling method has been demonstrated by applications to a range of chemical and biological problems: protein folding, molecular cluster structure searches, and chemical reactions. The combination of integrated tempering sampling and a trajectory sampling method allows the calculation of rate constants and reaction pathways without predefined collective coordinates.
Collapse
Affiliation(s)
| | | | - Qiang Shao
- Drug
Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | | | | |
Collapse
|
16
|
Shao Q. Important roles of hydrophobic interactions in folding and charge interactions in misfolding of α-helix bundle protein. RSC Adv 2015. [DOI: 10.1039/c4ra14265a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Abstract
An enhanced-sampling molecular dynamics simulation is presented to quantitatively demonstrate the important roles of hydrophobic and charge interactions in the folding and misfolding of α-helix bundle protein, respectively.
Collapse
Affiliation(s)
- Qiang Shao
- Drug Discovery and Design Center
- Shanghai Institute of Materia Medica
- Chinese Academy of Sciences
- Shanghai
- China
| |
Collapse
|
17
|
Wiebe H, Weinberg N. Theoretical volume profiles as a tool for probing transition states: folding kinetics. J Chem Phys 2014; 140:124105. [PMID: 24697422 DOI: 10.1063/1.4868549] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
The mechanism by which conformational changes, particularly folding and unfolding, occur in proteins and other biopolymers has been widely discussed in the literature. Molecular dynamics (MD) simulations of protein folding present a formidable challenge since these conformational changes occur on a time scale much longer than what can be afforded at the current level of computational technology. Transition state (TS) theory offers a more economic description of kinetic properties of a reaction system by relating them to the properties of the TS, or for flexible systems, the TS ensemble (TSE). The application of TS theory to protein folding is limited by ambiguity in the definition of the TSE for this process. We propose to identify the TSE for conformational changes in flexible systems by comparison of its experimentally determined volumetric property, known as the volume of activation, to the structure-specific volume profile of the process calculated using MD. We illustrate this approach by its successful application to unfolding of a model chain system.
Collapse
Affiliation(s)
- H Wiebe
- Department of Chemistry, Simon Fraser University, Burnaby, British Columbia V5A 1S6, Canada
| | - N Weinberg
- Department of Chemistry, Simon Fraser University, Burnaby, British Columbia V5A 1S6, Canada
| |
Collapse
|
18
|
Markiewicz BN, Yang L, Culik RM, Gao YQ, Gai F. How quickly can a β-hairpin fold from its transition state? J Phys Chem B 2014; 118:3317-25. [PMID: 24611730 PMCID: PMC3969101 DOI: 10.1021/jp500774q] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
![]()
Understanding the structural nature
of the free energy bottleneck(s)
encountered in protein folding is essential to elucidating the underlying
dynamics and mechanism. For this reason, several techniques, including
Φ-value analysis, have previously been developed to infer the
structural characteristics of such high free-energy or transition
states. Herein we propose that one (or few) appropriately placed backbone
and/or side chain cross-linkers, such as disulfides, could be used
to populate a thermodynamically accessible conformational state that
mimics the folding transition state. Specifically, we test this hypothesis
on a model β-hairpin, Trpzip4, as its folding mechanism has
been extensively studied and is well understood. Our results show
that cross-linking the two β-strands near the turn region increases
the folding rate by an order of magnitude, to about (500 ns)−1, whereas cross-linking the termini results in a hyperstable β-hairpin
that has essentially the same folding rate as the uncross-linked peptide.
Taken together, these findings suggest that cross-linking is not only
a useful strategy to manipulate folding free energy barriers, as shown
in other studies, but also, in some cases, it can be used to stabilize
a folding transition state analogue and allow for direct assessment
of the folding process on the downhill side of the free energy barrier.
The calculated free energy landscape of the cross-linked Trpzip4 also
supports this picture. An empirical analysis further suggests, when
folding of β-hairpins does not involve a significant free energy
barrier, the folding time (τ) follows a power law dependence
on the number of hydrogen bonds to be formed (nH), namely, τ = τ0nHα, with
τ0 = 20 ns and α = 2.3.
Collapse
Affiliation(s)
- Beatrice N Markiewicz
- Department of Chemistry and ‡Department of Biochemistry and Biophysics, University of Pennsylvania , Philadelphia, Pennsylvania 19104, United States
| | | | | | | | | |
Collapse
|
19
|
Shao Q, Wang J, Shi J, Zhu W. The universality of β-hairpin misfolding indicated by molecular dynamics simulations. J Chem Phys 2013; 139:165103. [DOI: 10.1063/1.4826461] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
20
|
Tsai M, Yuan J, Yamaki M, Lin C, Lin SH. Molecular Dynamics Insight into the Diverse Thermodynamic Behavior of a Beta‐Hairpin Peptide. J CHIN CHEM SOC-TAIP 2013. [DOI: 10.1002/jccs.201300173] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Min‐Yeh Tsai
- Department of Applied Chemistry, National Chiao Tung University, Hsinchu, Taiwan 30010, ROC
- Department of Chemistry, National Taiwan University, Taipei, Taiwan 10617, ROC
| | - Jian‐Min Yuan
- Department of Physics, Drexel University, Philadelphia, Pennsylvania 19104, USA
| | - Masahiro Yamaki
- Department of Applied Chemistry, National Chiao Tung University, Hsinchu, Taiwan 30010, ROC
| | - Chih‐Kai Lin
- Department of Applied Chemistry, National Chiao Tung University, Hsinchu, Taiwan 30010, ROC
| | - Sheng Hsien Lin
- Department of Applied Chemistry, National Chiao Tung University, Hsinchu, Taiwan 30010, ROC
| |
Collapse
|
21
|
Lin C, Culik RM, Gai F. Using VIPT-jump to distinguish between different folding mechanisms: application to BBL and a Trpzip. J Am Chem Soc 2013; 135:7668-73. [PMID: 23642153 PMCID: PMC3706100 DOI: 10.1021/ja401473m] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Protein folding involves a large number of sequential molecular steps or conformational substates. Thus, experimental characterization of the underlying folding energy landscape for any given protein is difficult. Herein, we present a new method that can be used to determine the major characteristics of the folding energy landscape in question, e.g., to distinguish between activated and barrierless downhill folding scenarios. This method is based on the idea that the conformational relaxation kinetics of different folding mechanisms at a given final condition will show different dependences on the initial condition. We show, using both simulation and experiment, that it is possible to differentiate between disparate kinetic folding models by comparing temperature jump (T-jump) relaxation traces obtained with a fixed final temperature and varied initial temperatures, which effectively varies the initial potential (VIP) of the system of interest. We apply this method (hereafter refer to as VIPT-jump) to two model systems, tryptophan zipper (Trpzip)-2c and BBL, and our results show that BBL exhibits characteristics of barrierless downhill folding, whereas Trpzip-2c folding encounters a free energy barrier. In addition, using the T-jump data of BBL we are able to provide, via Langevin dynamics simulations, a realistic estimate of its conformational diffusion coefficient.
Collapse
Affiliation(s)
- Chun–Wei Lin
- Department of Chemistry, University of Pennsylvania, Philadelphia, PA 19104
| | - Robert M. Culik
- Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, PA 19104
| | - Feng Gai
- Department of Chemistry, University of Pennsylvania, Philadelphia, PA 19104
| |
Collapse
|
22
|
Shao Q, Shi J, Zhu W. Molecular dynamics simulation indicating cold denaturation of β-hairpins. J Chem Phys 2013; 138:085102. [DOI: 10.1063/1.4792299] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
|
23
|
Lin Z, van Gunsteren WF. Combination of Enveloping Distribution Sampling (EDS) of a Soft-Core Reference-State Hamiltonian with One-Step Perturbation to Predict the Effect of Side Chain Substitution on the Relative Stability of Right- and Left-Helical Folds of β-Peptides. J Chem Theory Comput 2012; 9:126-34. [PMID: 26589016 DOI: 10.1021/ct300929q] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Folding free enthalpies of many not too different polypeptides can be efficiently and accurately predicted with the one-step perturbation (OSP) method using only one or a few molecular dynamics (MD) simulations. In this article, we introduce a combination of enveloping distribution sampling (EDS) and the OSP method (EDS-OSP) and apply it to predict the free enthalpy differences between a right-handed 2.710/12-helix and a left-handed 314-helix for 16 β-peptides with slightly different side-chain substitution patterns. An EDS simulation of a designed soft-core reference-state peptide was carried out in which both helices were sampled. Then, the soft-core atoms were perturbed into physical atoms. Thus, free enthalpy differences between the two helices for the 16 β-peptides can be predicted from only one simulation. The results predicted by EDS-OSP and a previous OSP study are very similar, i.e., the deviations between the results of the 16 peptides are mostly within the order of kBT, and the average absolute deviation is 1.2 kJ mol(-1). Together with the EDS parameter update simulation, about 128 ns of MD simulations needed to be carried out using the EDS-OSP method, while 700 ns of MD simulations were required in the previous OSP study where two separate reference-state simulations and an additional long time MD simulation of one of the 16 β-peptides were carried out. Thus, the computational effort was significantly reduced, i.e., by more than a factor of 5, using the EDS-OSP method. Hence, we consider this method an efficient tool to predict conformational free enthalpy differences from MD simulations.
Collapse
Affiliation(s)
- Zhixiong Lin
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093 Zürich, Switzerland
| | - Wilfred F van Gunsteren
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093 Zürich, Switzerland
| |
Collapse
|
24
|
Shao Q, Zhu W, Gao YQ. Robustness in Protein Folding Revealed by Thermodynamics Calculations. J Phys Chem B 2012; 116:13848-56. [DOI: 10.1021/jp307684h] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Qiang Shao
- Institute of Theoretical and
Computational Chemistry, College of Chemistry and Molecular Engineering,
Beijing National Laboratory of Molecular Sciences, Peking University, Beijing 100871, China
- Drug Discovery and Design Center,
Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203,
China
| | - Weiliang Zhu
- Drug Discovery and Design Center,
Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203,
China
| | - Yi Qin Gao
- Institute of Theoretical and
Computational Chemistry, College of Chemistry and Molecular Engineering,
Beijing National Laboratory of Molecular Sciences, Peking University, Beijing 100871, China
| |
Collapse
|
25
|
Shao Q, Shi J, Zhu W. Enhanced sampling molecular dynamics simulation captures experimentally suggested intermediate and unfolded states in the folding pathway of Trp-cage miniprotein. J Chem Phys 2012; 137:125103. [DOI: 10.1063/1.4754656] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Wu X, Yang G, Zu Y, Zhou L. Molecular dynamics studies of β-hairpin folding with the presence of the sodium ion. Comput Biol Chem 2012; 38:1-9. [PMID: 22487489 DOI: 10.1016/j.compbiolchem.2012.02.003] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2011] [Revised: 01/20/2012] [Accepted: 02/14/2012] [Indexed: 11/19/2022]
Abstract
Metal ions are ubiquitous in protein systems and play a significant role during their folding processes. Nineteen independent structures were determined for the Na(+)/β-hairpin interacting systems, and their folding pathways are different. (i) For Na(S47), the turn is rapidly shaped with the help of Na(+) and acts as the folding nucleus for the rest regions. Two intermediate states are observed and the resulted structure is the most folded. (ii) For Na(B41), Na(B52), Na(B54), Na(S55) and Na(B56), the inclusive Na(+) ions are anchored by β-strands. The local structures around the Na(+) ions and the turn regions fold simultaneously and serve as two independent folding nuclei. (iii) The other systems have no folding nuclei and correspond to low-folded structures. Long-range electrostatic interactions contribute a lot to the folding, especially from the four negatively charged residues (Glu42, Asp46, Asp47 and Glu56). The initial positions of the Na(+) ions are largely responsible for the different folding behaviors. The interactions with sidechain- rather than backbone-O atoms generally lead to more compact structures. Another factor affecting the folding is whether the O atoms are associated with native H-bonds, and those involved show decreased affinities to metal ions. The addition of water solvent does not induce obvious folding and conformational transitions to the Na(+)/β-hairpin interacting systems.
Collapse
Affiliation(s)
- Xiaomin Wu
- Key Laboratory of Forest Plant Ecology, Ministry of Education, Northeast Forestry University, Harbin 150040, PR China
| | | | | | | |
Collapse
|
27
|
Kaur H, Sasidhar YU. For the Sequence YKGQ, the Turn and Extended Conformational Forms Are Separated by Small Barriers and the Turn Propensity Persists Even at High Temperatures: Implications for Protein Folding. J Phys Chem B 2012; 116:3850-60. [DOI: 10.1021/jp210227s] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Affiliation(s)
- Harpreet Kaur
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400 076, India
| | - Yellamraju U. Sasidhar
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400 076, India
| |
Collapse
|
28
|
Shao Q, Yang L, Gao YQ. Structure change of β-hairpin induced by turn optimization: An enhanced sampling molecular dynamics simulation study. J Chem Phys 2011; 135:235104. [DOI: 10.1063/1.3668288] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
29
|
Xu Y, Du D, Oyola R. Infrared study of the stability and folding kinetics of a series of β-hairpin peptides with a common NPDG turn. J Phys Chem B 2011; 115:15332-8. [PMID: 22136248 DOI: 10.1021/jp2046867] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The thermal stability and folding kinetics of a series of 15-residue β-hairpins with a common Type I [3:5] NPDG turn were studied using Fourier transform infrared spectroscopy (FTIR) and laser-induced temperature jump (T-jump) with infrared detection, respectively. Mutations at positions 3, 5, or 13 in the peptide sequence SEXYXNPDGTWTXTE, where X represents the position of mutation, were performed to study the roles of hydrophobic interactions in determining the thermodynamic and kinetic properties of β-hairpin folding. The thermal stability studies show a broad thermal folding/unfolding transition for all the peptides. T-jump studies indicate that these β-hairpin peptides fold in less than 2 μs. In addition, both folding and unfolding rate constants decrease with increasing strength of hydrophobic interactions. Kinetically, the hydrophobic interactions have more significant influence on the unfolding rate than the folding rate. Φ-value analysis indicates that the hydrophobic interactions between the side chains are mainly formed at the latter part of the transition-state region during the folding process. In summary, the results suggest that the formation of the native structure of these β-hairpins depends on the correct topology of the hydrophobic cluster. Besides the formation of the turn region as a key process for folding as suggested by previous studies, a hydrophobic collapse process may also play a crucial role during β-hairpin folding.
Collapse
Affiliation(s)
- Yao Xu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | | | | |
Collapse
|
30
|
Shao Q, Qin Gao Y. The relative helix and hydrogen bond stability in the B domain of protein A as revealed by integrated tempering sampling molecular dynamics simulation. J Chem Phys 2011; 135:135102. [DOI: 10.1063/1.3630127] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
|
31
|
Shao Q, Gao YQ. Temperature Dependence of Hydrogen-Bond Stability in β-Hairpin Structures. J Chem Theory Comput 2010. [DOI: 10.1021/ct100436r] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Qiang Shao
- College of Chemistry and Molecular Engineering, National Laboratory of Molecular Sciences, Peking University, Beijing, China
| | - Yi Qin Gao
- College of Chemistry and Molecular Engineering, National Laboratory of Molecular Sciences, Peking University, Beijing, China
| |
Collapse
|