1
|
Ali AE, Li LL, Courtney MJ, Pentikäinen OT, Postila PA. Atomistic simulations reveal impacts of missense mutations on the structure and function of SynGAP1. Brief Bioinform 2024; 25:bbae458. [PMID: 39311700 PMCID: PMC11418247 DOI: 10.1093/bib/bbae458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2024] [Revised: 07/20/2024] [Accepted: 09/04/2024] [Indexed: 09/26/2024] Open
Abstract
De novo mutations in the synaptic GTPase activating protein (SynGAP) are associated with neurological disorders like intellectual disability, epilepsy, and autism. SynGAP is also implicated in Alzheimer's disease and cancer. Although pathogenic variants are highly penetrant in neurodevelopmental conditions, a substantial number of them are caused by missense mutations that are difficult to diagnose. Hence, in silico mutagenesis was performed for probing the missense effects within the N-terminal region of SynGAP structure. Through extensive molecular dynamics simulations, encompassing three 150-ns replicates for 211 variants, the impact of missense mutations on the protein fold was assessed. The effect of the mutations on the folding stability was also quantitatively assessed using free energy calculations. The mutations were categorized as potentially pathogenic or benign based on their structural impacts. Finally, the study introduces wild-type-SynGAP in complex with RasGTPase at the inner membrane, while considering the potential effects of mutations on these key interactions. This study provides structural perspective to the clinical assessment of SynGAP missense variants and lays the foundation for future structure-based drug discovery.
Collapse
Affiliation(s)
- Aliaa E Ali
- MedChem.fi, Institute of Biomedicine, Integrative Physiology and Pharmacology, University of Turku, FI-20014 Turku, Finland
- InFLAMES Research Flagship, University of Turku, 20014 Turku, Finland
| | - Li-Li Li
- Neuronal Signalling Laboratory and Turku Screening Unit, Turku Bioscience Centre, University of Turku and Åbo Akademi University, Turku, Finland
| | - Michael J Courtney
- Neuronal Signalling Laboratory and Turku Screening Unit, Turku Bioscience Centre, University of Turku and Åbo Akademi University, Turku, Finland
| | - Olli T Pentikäinen
- MedChem.fi, Institute of Biomedicine, Integrative Physiology and Pharmacology, University of Turku, FI-20014 Turku, Finland
- InFLAMES Research Flagship, University of Turku, 20014 Turku, Finland
| | - Pekka A Postila
- MedChem.fi, Institute of Biomedicine, Integrative Physiology and Pharmacology, University of Turku, FI-20014 Turku, Finland
- InFLAMES Research Flagship, University of Turku, 20014 Turku, Finland
| |
Collapse
|
2
|
Roterman I, Slupina M, Stapor K, Konieczny L, Gądek K, Nowakowski P. Chameleon Sequences-Structural Effects in Proteins Characterized by Hydrophobicity Disorder. ACS OMEGA 2024; 9:38506-38522. [PMID: 39310170 PMCID: PMC11411663 DOI: 10.1021/acsomega.4c03658] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/16/2024] [Revised: 08/20/2024] [Accepted: 08/23/2024] [Indexed: 09/25/2024]
Abstract
Repeated protein folding processes both in vivo and in vitro leading to the same structure for a specific amino acid sequence prove that the amino acid sequence determines protein structuring. This is also evidenced by the variability of structuring, dependent on the introduced mutations. An important phenomenon in this regard is the presence of a differentiated secondary structure for chain fragments of identical sequence representing distinct forms of the secondary-order structure. Proteins termed chameleon proteins contain polypeptide chain fragments of identical sequence (length 6-12 aa) showing structural differentiation: helix versus β-structure. In the present paper, it was shown that these fragments represent components matching the structural status dictated by the physicochemical properties of the entire structural unit. This structural matching is related to achieving the goal of the biological function of the structural unit. The corresponding secondary structure represents a means to achieving this goal, not an end in itself. A selected set of proteins from the ChSeq database have been analyzed using a fuzzy oil drop model (FOD-M) identifying the uniqueness of the hydrophobicity distribution taken as a medium for recording the specificity of a given protein and a given chameleon section in particular. It was shown that in the vast majority, the status of chameleon sections turns out to be comparable regardless of the represented secondary structure.
Collapse
Affiliation(s)
- Irena Roterman
- Department
of Bioinformatics and Telemedicine, Jagiellonian
University—Medical College, Medyczna 7, 30-688 Krakow, Poland
| | - Mateusz Slupina
- ALSTOM
ZWUS Sp. z o.o., Modelarska
12, 40-142 Katowice, Poland
| | - Katarzyna Stapor
- Faculty
of Automatic, Electronics and Computer Science, Department of Applied
Informatics, Silesian University of Technology, Akademicka 16, 44-100 Gliwice, Poland
| | - Leszek Konieczny
- Chair
of Medical Biochemistry, Jagiellonian University—Medical
College, Kopernika 7, 31-034 Krakow, Poland
| | - Krzysztof Gądek
- AGH
Cyfronet, SANO SCIENCE, Nawojki 11, 30-950 Kraków, Poland
| | | |
Collapse
|
3
|
Moses K, Van Tassel PR. Polyelectrolyte Influence on Beta-Hairpin Peptide Stability: A Simulation Study. J Phys Chem B 2023; 127:359-370. [PMID: 36574611 DOI: 10.1021/acs.jpcb.2c06641] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Assemblies of proteins and charged macromolecules (polyelectrolytes) find important applications as pharmaceutical formulations, biocatalysts, and cell-contacting substrates. A key question is how the polymer component influences the structure and function of the protein. The present paper addresses the influence of charged polymers on the thermal stability of two model beta-hairpin-forming peptides through an all-atom, replica exchange molecular dynamics simulation. The (negatively charged) peptides consist of the terminal 16 amino acids of the B1 domain of Protein G (GB1) and a variant with three of the GB1 residues substituted with tryptophan (Tryptophan Zipper 4, or TZ4). A (cationic) lysine polymer is seen to thermally stabilize TZ4 and destabilize GB1, while a (also cationic) chitosan polymer slightly stabilizes GB1 but has essentially no effect on TZ4. Free energy profiles reveal folded and unfolded conformations to be separated by kinetic barriers generally acting in the direction of the thermodynamically favored state. Through application of an Ising-like statistical mechanical model, a mechanism is proposed based on competition between (indirect) entropic stabilization of folded versus unfolded states and (direct) competition for hydrogen-bonding and hydrophobic interactions. These findings have important implications to the design of polyelectrolyte-based materials for biomedical and biotechnological applications.
Collapse
Affiliation(s)
- Kevin Moses
- Dept. of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520-8286, United States
| | - Paul R Van Tassel
- Dept. of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520-8286, United States
| |
Collapse
|
4
|
Bandyopadhyay S, Mondal J. A deep autoencoder framework for discovery of metastable ensembles in biomacromolecules. J Chem Phys 2021; 155:114106. [PMID: 34551528 DOI: 10.1063/5.0059965] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
Biomacromolecules manifest dynamic conformational fluctuation and involve mutual interconversion among metastable states. A robust mapping of their conformational landscape often requires the low-dimensional projection of the conformational ensemble along optimized collective variables (CVs). However, the traditional choice for the CV is often limited by user-intuition and prior knowledge about the system, and this lacks a rigorous assessment of their optimality over other candidate CVs. To address this issue, we propose an approach in which we first choose the possible combinations of inter-residue Cα-distances within a given macromolecule as a set of input CVs. Subsequently, we derive a non-linear combination of latent space embedded CVs via auto-encoding the unbiased molecular dynamics simulation trajectories within the framework of the feed-forward neural network. We demonstrate the ability of the derived latent space variables in elucidating the conformational landscape in four hierarchically complex systems. The latent space CVs identify key metastable states of a bead-in-a-spring polymer. The combination of the adopted dimensional reduction technique with a Markov state model, built on the derived latent space, reveals multiple spatially and kinetically well-resolved metastable conformations for GB1 β-hairpin. A quantitative comparison based on the variational approach-based scoring of the auto-encoder-derived latent space CVs with the ones obtained via independent component analysis (principal component analysis or time-structured independent component analysis) confirms the optimality of the former. As a practical application, the auto-encoder-derived CVs were found to predict the reinforced folding of a Trp-cage mini-protein in aqueous osmolyte solution. Finally, the protocol was able to decipher the conformational heterogeneities involved in a complex metalloenzyme, namely, cytochrome P450.
Collapse
Affiliation(s)
- Satyabrata Bandyopadhyay
- Tata Institute of Fundamental Research, Center for Interdisciplinary Sciences, Hyderabad 500046, India
| | - Jagannath Mondal
- Tata Institute of Fundamental Research, Center for Interdisciplinary Sciences, Hyderabad 500046, India
| |
Collapse
|
5
|
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
|
6
|
Zerze GH, Stillinger FH, Debenedetti PG. Effect of heterochiral inversions on the structure of a β-hairpin peptide. Proteins 2019; 87:569-578. [PMID: 30811673 DOI: 10.1002/prot.25680] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2018] [Accepted: 02/24/2019] [Indexed: 01/25/2023]
Abstract
We study computationally a family of β-hairpin peptides with systematically introduced chiral inversions, in explicit water, and we investigate the extent to which the backbone structure is able to fold in the presence of heterochiral perturbations. In contrast to the recently investigated case of a helical peptide, we do not find a monotonic change in secondary structure content as a function of the number of L- to D-inversions. The effects of L- to D-inversions are instead found to be highly position-specific. Additionally, in contrast to the helical peptide, some inversions increase the stability of the folded peptide: in such cases, we compute an increase in β-sheet content in the aqueous solution equilibrium ensemble. However, the tertiary structures of the stable (folded) configurations for peptides for which inversions cause an increase in β-sheet content show differences from one another, as well as from the native fold of the nonchirally perturbed β-hairpin. Our results suggest that although some chiral perturbations can increase folding stability, chirally perturbed proteins may still underperform functionally, given the relationship between structure and function.
Collapse
Affiliation(s)
- Gül H Zerze
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey
| | | | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey
| |
Collapse
|
7
|
Ahalawat N, Mondal J. Assessment and optimization of collective variables for protein conformational landscape: GB1 β-hairpin as a case study. J Chem Phys 2018; 149:094101. [PMID: 30195312 DOI: 10.1063/1.5041073] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Collective variables (CVs), when chosen judiciously, can play an important role in recognizing rate-limiting processes and rare events in any biomolecular systems. However, high dimensionality and inherent complexities associated with such biochemical systems render the identification of an optimal CV a challenging task, which in turn precludes the elucidation of an underlying conformational landscape in sufficient details. In this context, a relevant model system is presented by a 16-residue β-hairpin of GB1 protein. Despite being the target of numerous theoretical and computational studies for understanding the protein folding, the set of CVs optimally characterizing the conformational landscape of the β-hairpin of GB1 protein has remained elusive, resulting in a lack of consensus on its folding mechanism. Here we address this by proposing a pair of optimal CVs which can resolve the underlying free energy landscape of the GB1 hairpin quite efficiently. Expressed as a linear combination of a number of traditional CVs, the optimal CV for this system is derived by employing the recently introduced time-structured independent component analysis approach on a large number of independent unbiased simulations. By projecting the replica-exchange simulated trajectories along these pair of optimized CVs, the resulting free energy landscape of this system is able to resolve four distinct well-separated metastable states encompassing the extensive ensembles of folded, unfolded, and molten globule states. Importantly, the optimized CVs were found to be capable of automatically recovering a novel partial helical state of this protein, without needing to explicitly invoke helicity as a constituent CV. Furthermore, a quantitative sensitivity analysis of each constituent in the optimized CV provided key insights on the relative contributions of the constituent CVs in the overall free energy landscapes. Finally, the kinetic pathways connecting these metastable states, constructed using a Markov state model, provide an optimum description of the underlying folding mechanism of the peptide. Taken together, this work offers a quantitatively robust approach toward comprehensive mapping of the underlying folding landscape of a quintessential model system along its optimized CV.
Collapse
Affiliation(s)
- Navjeet Ahalawat
- Tata Institute of Fundamental Research, Center for Interdisciplinary Sciences, Hyderabad 500107, India
| | - Jagannath Mondal
- Tata Institute of Fundamental Research, Center for Interdisciplinary Sciences, Hyderabad 500107, India
| |
Collapse
|
8
|
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
|
9
|
Gupta M, Khatua P, Chakravarty C, Bandyopadhyay S. Hydration Behavior along the Folding Pathways of Trpzip4, Trpzip5 and Trpzip6. J Phys Chem B 2018; 122:1560-1572. [DOI: 10.1021/acs.jpcb.7b10135] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
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, Kharagpur 721302, India
| | - Charusita Chakravarty
- Department
of Chemistry, Indian Institute of Technology-Delhi, New Delhi 110016, India
| | - Sanjoy Bandyopadhyay
- Molecular
Modeling Laboratory, Department of Chemistry, Indian Institute of Technology-Kharagpur, Kharagpur 721302, India
| |
Collapse
|
10
|
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
|
11
|
Xiao J, Li Y, Huang Q. Application of Monte Carlo simulation in addressing key issues of complex coacervation formed by polyelectrolytes and oppositely charged colloids. Adv Colloid Interface Sci 2017; 239:31-45. [PMID: 27265512 DOI: 10.1016/j.cis.2016.05.010] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2016] [Revised: 05/18/2016] [Accepted: 05/21/2016] [Indexed: 10/21/2022]
Abstract
This paper reviews the recent advance of Monte Carlo (MC) simulation in addressing key issues of complex coacervation between polyelectrolytes and oppositely charged colloids. Readers were first supplied with a brief overview of current knowledge and experimental strategies in the study of complex coacervation. In the next section, the general MC simulation procedures as well as representative strategies applied in complex coacervation were summarized. The unique contributions of MC simulation in either capturing delicate features, easing the experimental trials or proving the concept were then elucidated through the following aspects: i) identify phase boundary and decouple interaction contributions; ii) clarify composition distribution and internal structure; iii) predict the influences of physicochemical conditions on complex coacervation; iv) delineate the mechanisms for "binding on the wrong side of the isoelectric point". Finally, current challenges as well as prospects of MC simulation in complex coacervation are also discussed. The ultimate goal of this review is to provide readers with basic guideline for synergistic design of experiments in combination with MC simulation, and deliver convincing interpretation and reliable prediction for the structure and behavior in polyelectrolyte-macroion complex coacervation.
Collapse
|
12
|
Yu TQ, Lu J, Abrams CF, Vanden-Eijnden E. Multiscale implementation of infinite-swap replica exchange molecular dynamics. Proc Natl Acad Sci U S A 2016; 113:11744-11749. [PMID: 27698148 PMCID: PMC5081654 DOI: 10.1073/pnas.1605089113] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Replica exchange molecular dynamics (REMD) is a popular method to accelerate conformational sampling of complex molecular systems. The idea is to run several replicas of the system in parallel at different temperatures that are swapped periodically. These swaps are typically attempted every few MD steps and accepted or rejected according to a Metropolis-Hastings criterion. This guarantees that the joint distribution of the composite system of replicas is the normalized sum of the symmetrized product of the canonical distributions of these replicas at the different temperatures. Here we propose a different implementation of REMD in which (i) the swaps obey a continuous-time Markov jump process implemented via Gillespie's stochastic simulation algorithm (SSA), which also samples exactly the aforementioned joint distribution and has the advantage of being rejection free, and (ii) this REMD-SSA is combined with the heterogeneous multiscale method to accelerate the rate of the swaps and reach the so-called infinite-swap limit that is known to optimize sampling efficiency. The method is easy to implement and can be trivially parallelized. Here we illustrate its accuracy and efficiency on the examples of alanine dipeptide in vacuum and C-terminal β-hairpin of protein G in explicit solvent. In this latter example, our results indicate that the landscape of the protein is a triple funnel with two folded structures and one misfolded structure that are stabilized by H-bonds.
Collapse
Affiliation(s)
- Tang-Qing Yu
- Courant Institute of Mathematical Sciences, New York University, New York, NY 10012
| | - Jianfeng Lu
- Department of Mathematics, Duke University, Durham, NC 27708; Department of Physics, Duke University, Durham, NC 27708; Department of Chemistry, Duke University, Durham, NC 27708
| | - Cameron F Abrams
- Department of Chemical and Biological Engineering, Drexel University, Philadelphia, PA 19104
| | - Eric Vanden-Eijnden
- Courant Institute of Mathematical Sciences, New York University, New York, NY 10012;
| |
Collapse
|
13
|
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
|
14
|
Paul S, Taraphder S. Determination of the Reaction Coordinate for a Key Conformational Fluctuation in Human Carbonic Anhydrase II. J Phys Chem B 2015; 119:11403-15. [DOI: 10.1021/acs.jpcb.5b03655] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Sanjib Paul
- Department of Chemistry, Indian Institute of Technology, Kharagpur 721302, India
| | - Srabani Taraphder
- Department of Chemistry, Indian Institute of Technology, Kharagpur 721302, India
| |
Collapse
|
15
|
Zerze GH, Uz B, Mittal J. Folding thermodynamics ofβ-hairpins studied by replica-exchange molecular dynamics simulations. Proteins 2015; 83:1307-15. [DOI: 10.1002/prot.24827] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2015] [Revised: 04/24/2015] [Accepted: 04/29/2015] [Indexed: 01/21/2023]
Affiliation(s)
- Gül H. Zerze
- Department of Chemical and Biomolecular Engineering; Lehigh University; Bethlehem Pennsylvania 18015
| | - Bilge Uz
- Department of Chemical and Biomolecular Engineering; Lehigh University; Bethlehem Pennsylvania 18015
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering; Lehigh University; Bethlehem Pennsylvania 18015
| |
Collapse
|
16
|
Razavi AM, Voelz VA. Kinetic Network Models of Tryptophan Mutations in β-Hairpins Reveal the Importance of Non-Native Interactions. J Chem Theory Comput 2015; 11:2801-12. [DOI: 10.1021/acs.jctc.5b00088] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Asghar M. Razavi
- Department
of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Vincent A. Voelz
- Department
of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
17
|
Peter EK, Pivkin IV, Shea JE. A kMC-MD method with generalized move-sets for the simulation of folding of α-helical and β-stranded peptides. J Chem Phys 2015; 142:144903. [DOI: 10.1063/1.4915919] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- Emanuel K. Peter
- Faculty of Informatics, Institute of Computational Science, University of Lugano, Lugano, Switzerland
| | - Igor V. Pivkin
- Faculty of Informatics, Institute of Computational Science, University of Lugano, Lugano, Switzerland
| | - Joan-Emma Shea
- Department of Chemistry and Biochemistry, and Department of Physics, University of California, Santa Barbara, California 93106, USA
| |
Collapse
|
18
|
Ardevol A, Tribello GA, Ceriotti M, Parrinello M. Probing the Unfolded Configurations of a β-Hairpin Using Sketch-Map. J Chem Theory Comput 2015; 11:1086-93. [DOI: 10.1021/ct500950z] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Affiliation(s)
- Albert Ardevol
- Computational
Science, Department of Chemistry and Applied Biosciences, ETH Zurich, USI-Campus, Via Giuseppe Buffi 13, C-6900 Lugano, Switzerland
| | - Gareth A. Tribello
- Atomistic
Simulation Centre, School of Mathematics and Physics, Queen’s University Belfast, Belfast BT7 1NN, United Kingdom
| | - Michele Ceriotti
- Laboratory
of Computational Science and Modelling, EPFL, CH-1015 Lausanne, Switzerland
| | - Michele Parrinello
- Computational
Science, Department of Chemistry and Applied Biosciences, ETH Zurich, USI-Campus, Via Giuseppe Buffi 13, C-6900 Lugano, Switzerland
| |
Collapse
|
19
|
Voelz VA, Elman B, Razavi AM, Zhou G. Surprisal Metrics for Quantifying Perturbed Conformational Dynamics in Markov State Models. J Chem Theory Comput 2014; 10:5716-28. [DOI: 10.1021/ct500827g] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Vincent A. Voelz
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Brandon Elman
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Asghar M. Razavi
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Guangfeng Zhou
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
20
|
Narayanan C, Dias CL. Exploring the free energy landscape of a model β-hairpin peptide and its isoform. Proteins 2014; 82:2394-402. [PMID: 24825659 DOI: 10.1002/prot.24601] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2013] [Revised: 03/21/2014] [Accepted: 04/29/2014] [Indexed: 12/16/2022]
Abstract
Secondary structural transitions from α-helix to β-sheet conformations are observed in several misfolding diseases including Alzheimer's and Parkinson's. Determining factors contributing favorably to the formation of each of these secondary structures is therefore essential to better understand these disease states. β-hairpin peptides form basic components of anti-parallel β-sheets and are suitable model systems for characterizing the fundamental forces stabilizing β-sheets in fibrillar structures. In this study, we explore the free energy landscape of the model β-hairpin peptide GB1 and its E2 isoform that preferentially adopts α-helical conformations at ambient conditions. Umbrella sampling simulations using all-atom models and explicit solvent are performed over a large range of end-to-end distances. Our results show the strong preference of GB1 and the E2 isoform for β-hairpin and α-helical conformations, respectively, consistent with previous studies. We show that the unfolded states of GB1 are largely populated by misfolded β-hairpin structures which differ from each other in the position of the β-turn. We discuss the energetic factors contributing favorably to the formation of α-helix and β-hairpin conformations in these peptides and highlight the energetic role of hydrogen bonds and non-bonded interactions.
Collapse
Affiliation(s)
- Chitra Narayanan
- Department of Physics, New Jersey Institute of Technology, University Heights, Newark, New Jersey, 07102-1982
| | | |
Collapse
|
21
|
Lee J. Exact partition function zeros of the Wako-Saitô-Muñoz-Eaton β hairpin model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:022710. [PMID: 24032867 DOI: 10.1103/physreve.88.022710] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/01/2013] [Indexed: 06/02/2023]
Abstract
I compute exact partition function zeros of β hairpins, using both analytic and numerical methods, extending previous work [J. Lee, Phys. Rev. Lett. 110, 248101 (2013)] where only a restricted class of hairpins was considered. The zeros of β hairpins with an odd number of peptide bonds are computed and the difference of the distribution of zeros from those for an even number of peptide bonds is explained in terms of additional entropy of liberating the extra bond at the turn region. Upon the introduction of a hydrophobic core in the central region of the hairpin, the zeros are distributed uniformly on two concentric circles corresponding to the hydrophobic collapse and the transition to the fully folded conformation. One of the circles dissolves as the core moves toward the turn or the tip region, which is explained in terms of the similarity of the intermediate state with the folded or unfolded states. The exact partition function zeros for a hairpin with a more complex structure of native contacts, the 16 C-terminal residues of streptococcal protein G B1, are numerically computed and their loci are closely approximated by concentric circles.
Collapse
Affiliation(s)
- Julian Lee
- Department of Bioinformatics and Life Science, Soongsil University, Seoul 156-743, Korea
| |
Collapse
|
22
|
Tsai MY, Yuan JM, Teranishi Y, Lin SH. Thermodynamics of protein folding using a modified Wako-Saitô-Muñoz-Eaton model. J Biol Phys 2012; 38:543-71. [PMID: 24615219 PMCID: PMC3473134 DOI: 10.1007/s10867-012-9271-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2011] [Accepted: 05/07/2012] [Indexed: 10/28/2022] Open
Abstract
Herein, we propose a modified version of the Wako-Saitô-Muñoz-Eaton (WSME) model. The proposed model introduces an empirical temperature parameter for the hypothetical structural units (i.e., foldons) in proteins to include site-dependent thermodynamic behavior. The thermodynamics for both our proposed model and the original WSME model were investigated. For a system with beta-hairpin topology, a mathematical treatment (contact-pair treatment) to facilitate the calculation of its partition function was developed. The results show that the proposed model provides better insight into the site-dependent thermodynamic behavior of the system, compared with the original WSME model. From this site-dependent point of view, the relationship between probe-dependent experimental results and model's thermodynamic predictions can be explained. The model allows for suggesting a general principle to identify foldon behavior. We also find that the backbone hydrogen bonds may play a role of structural constraints in modulating the cooperative system. Thus, our study may contribute to the understanding of the fundamental principles for the thermodynamics of protein folding.
Collapse
Affiliation(s)
- Min-Yeh Tsai
- National Chiao Tung University, 1001 Ta Hsuen Road, Hsinchu, Taiwan, Republic of China,
| | | | | | | |
Collapse
|
23
|
|
24
|
Barducci A, Bonomi M, Derreumaux P. Assessing the Quality of the OPEP Coarse-Grained Force Field. J Chem Theory Comput 2011; 7:1928-34. [DOI: 10.1021/ct100646f] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Alessandro Barducci
- Computational Science, Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, via Buffi 13, CH-6900 Lugano, Switzerland
| | - Massimiliano Bonomi
- Computational Science, Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, via Buffi 13, CH-6900 Lugano, Switzerland
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, Institut de Biologie Physico-Chimique and Université Paris Diderot, Paris 7, Institut Universitaire de France, 13 rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
25
|
Hayre NR, Singh RRP, Cox DL. Evaluating force field accuracy with long-time simulations of a β-hairpin tryptophan zipper peptide. J Chem Phys 2011; 134:035103. [DOI: 10.1063/1.3532931] [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
|
26
|
Zou X, Liu Y, Chen Z, Cárdenas-Jirón GI, Schulten K. Flow-induced beta-hairpin folding of the glycoprotein Ibalpha beta-switch. Biophys J 2010; 99:1182-91. [PMID: 20713002 DOI: 10.1016/j.bpj.2010.05.035] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2010] [Revised: 05/25/2010] [Accepted: 05/27/2010] [Indexed: 11/19/2022] Open
Abstract
Flow-induced shear has been identified as a regulatory driving force in blood clotting. Shear induces beta-hairpin folding of the glycoprotein Ibalpha beta-switch which increases affinity for binding to the von Willebrand factor, a key step in blood clot formation and wound healing. Through 2.1-micros molecular dynamics simulations, we investigate the kinetics of flow-induced beta-hairpin folding. Simulations sampling different flow velocities reveal that under flow, beta-hairpin folding is initiated by hydrophobic collapse, followed by interstrand hydrogen-bond formation and turn formation. Adaptive biasing force simulations are employed to determine the free energy required for extending the unfolded beta-switch from a loop to an elongated state. Lattice and freely jointed chain models illustrate how the folding rate depends on the entropic and enthalpic energy, the latter controlled by flow. The results reveal that the free energy landscape of the beta-switch has two stable conformations imprinted on it, namely, loop and hairpin--with flow inducing a transition between the two.
Collapse
Affiliation(s)
- Xueqing Zou
- School of Physics, Peking University, Beijing, China; Beckman Institute, University of Illinois, Urbana, Illinois, USA
| | | | | | | | | |
Collapse
|