1
|
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
|
2
|
Man VH, He X, Wang J. Stable Cavitation Interferes with Aβ 16-22 Oligomerization. J Chem Inf Model 2022; 62:3885-3895. [PMID: 35920625 DOI: 10.1021/acs.jcim.2c00764] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Ultrasound and microbubbles are used for many medical applications nowadays. Scanning ultrasound can remove amyloid-β (Aβ) aggregates in the mouse brain and restores memory in an Alzheimer's disease mouse model. In vitro studies showed that amyloid fibrils are fragmented due to the ultrasound-induced bubble inertial cavitation, and ultrasonic pulses accelerate the depolymerization of Aβ fibrils into monomers at 1 μM of concentration. Under applied ultrasound, microbubbles can be in a stable oscillating state or unstable inertial cavitation state. The latter occurs when ultrasound causes a dramatic change of bubble sizes above a certain acoustic pressure. We have developed and implemented a nonequilibrium molecular dynamics simulation algorithm to the AMBER package, to facilitate the investigation of the molecular mechanism of Aβ oligomerization under stable cavitation. Our results indicated that stable cavitation not only inhibited oligomeric formation, but also prevented the formation of β-rich oligomers. The network analysis of state transitions revealed that stable cavitation altered the oligomerization pathways of Aβ16-22 peptides. Our simulation tool may be applied to optimize the experimental conditions to achieve the best therapeutical effect.
Collapse
Affiliation(s)
- Viet Hoang Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
3
|
Kamenik AS, Linker SM, Riniker S. Enhanced sampling without borders: on global biasing functions and how to reweight them. Phys Chem Chem Phys 2022; 24:1225-1236. [PMID: 34935813 PMCID: PMC8768491 DOI: 10.1039/d1cp04809k] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Accepted: 12/14/2021] [Indexed: 12/17/2022]
Abstract
Molecular dynamics (MD) simulations are a powerful tool to follow the time evolution of biomolecular motions in atomistic resolution. However, the high computational demand of these simulations limits the timescales of motions that can be observed. To resolve this issue, so called enhanced sampling techniques are developed, which extend conventional MD algorithms to speed up the simulation process. Here, we focus on techniques that apply global biasing functions. We provide a broad overview of established enhanced sampling methods and promising new advances. As the ultimate goal is to retrieve unbiased information from biased ensembles, we also discuss benefits and limitations of common reweighting schemes. In addition to concisely summarizing critical assumptions and implications, we highlight the general application opportunities as well as uncertainties of global enhanced sampling.
Collapse
Affiliation(s)
- Anna S Kamenik
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
| | - Stephanie M Linker
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
| | - Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland.
| |
Collapse
|
4
|
Wang D, Wang Y, Chang J, Zhang L, Wang H, E W. Efficient sampling of high-dimensional free energy landscapes using adaptive reinforced dynamics. NATURE COMPUTATIONAL SCIENCE 2022; 2:20-29. [PMID: 38177702 DOI: 10.1038/s43588-021-00173-1] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/04/2021] [Accepted: 11/15/2021] [Indexed: 01/06/2024]
Abstract
Enhanced sampling methods such as metadynamics and umbrella sampling have become essential tools for exploring the configuration space of molecules and materials. At the same time, they have long faced a number of issues such as the inefficiency when dealing with a large number of collective variables (CVs) or systems with high free energy barriers. Here we show that, with clustering and adaptive tuning techniques, the reinforced dynamics (RiD) scheme can be used to efficiently explore the configuration space and free energy landscapes with a large number of CVs or systems with high free energy barriers. We illustrate this by studying various representative and challenging examples. First we demonstrate the efficiency of adaptive RiD compared with other methods and construct the nine-dimensional (9D) free energy landscape of a peptoid trimer, which has energy barriers of more than 8 kcal mol-1. We then study the folding of the protein chignolin using 18 CVs. In this case, both the folding and unfolding rates are observed to be 4.30 μs-1. Finally, we propose a protein structure refinement protocol based on RiD. This protocol allows us to efficiently employ more than 100 CVs for exploring the landscape of protein structures and it gives rise to an overall improvement of 14.6 units over the initial global distance test-high accuracy (GDT-HA) score.
Collapse
Affiliation(s)
- Dongdong Wang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA
- DP Technology, Beijing, People's Republic of China
| | - Yanze Wang
- DP Technology, Beijing, People's Republic of China
- College of Chemistry and Molecular Engineering, Peking University, Beijing, People's Republic of China
| | - Junhan Chang
- DP Technology, Beijing, People's Republic of China
- College of Chemistry and Molecular Engineering, Peking University, Beijing, People's Republic of China
| | - Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA.
- DP Technology, Beijing, People's Republic of China.
| | - Han Wang
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing, People's Republic of China.
| | - Weinan E
- School of Mathematical Sciences, Peking University, Beijing, People's Republic of China
- Department of Mathematics and Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA
- Beijing Institute of Big Data Research, Beijing, People's Republic of China
| |
Collapse
|
5
|
Nguyen PH, Derreumaux P. Computer Simulations Aimed at Exploring Protein Aggregation and Dissociation. Methods Mol Biol 2022; 2340:175-196. [PMID: 35167075 DOI: 10.1007/978-1-0716-1546-1_9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.
Collapse
Affiliation(s)
- Phuong H Nguyen
- Laboratoire de Biochimie Théorique, UPR 9080, CNRS, Université de Paris, Paris, France
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, UPR 9080, CNRS, Université de Paris, Paris, France.
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris, France.
| |
Collapse
|
6
|
Wada H, Okamoto Y. Two-dimensional simulated tempering for the isobaric–isothermal ensemble with fast on-the-fly weight determination. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1904156] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Hiromune Wada
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Japan
| | - Yuko Okamoto
- Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Japan
- Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya, Japan
- Information Technology Center, Nagoya University, Nagoya, Japan
| |
Collapse
|
7
|
Man VH, He X, Gao J, Wang J. Effects of All-Atom Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of PHF6 Peptide of Tau Protein. J Chem Theory Comput 2021; 17:6458-6471. [PMID: 34491058 DOI: 10.1021/acs.jctc.1c00028] [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/05/2023]
Abstract
Molecular dynamics (MD) simulations play a vital role in revealing the mechanism of amyloid aggregation that is crucial to the therapeutic agent development for Alzheimer's Disease. However, the accuracy of MD simulation results strongly depends on the force field employed. In our previous benchmark for 17 all-atom force fields on modeling of amyloid aggregation using the Aβ16-22 dimer, we showed that AMBER14SB and CHARMM36m are suitable force fields for amyloid aggregation simulation, while GROMOS54a7 and OPLSAA are not good for the task. In this work, we continue assessing the applicability of atomistic force fields on amyloid aggregation using the VQIVYK (PHF6) peptide which is essential for tau-protein aggregation. Although, both Aβ16-22 and PHF6 peptides formed fibrils in vitro, the PHF6 fibrils are parallel β-sheets, while the Aβ16-22 fibrils are antiparallel β-sheets. We performed an all-atom large-scale MD simulation in explicit water on the PHF6 dimer and octa-peptides systems using five mainstream force fields, including AMBER99SB-disp, AMBER14SB, CHARMM36m, GROMOS54a7, and OPLSAA. The accumulated simulation time is 0.2 ms. Our result showed that the β-sheet structures of PHF6 peptides sampled by AMBER99SB-disp, AMBER14SB, GROMOS54a7, and OPLSAA are in favor of the antiparallel β-sheets, while the dominant type of β-sheet structures is parallel β-sheet by using CHARMM36m. Among the five force fields, CHARMM36m provides the strongest CH-π interaction that was observed in an NMR study. The comparison between our results and experimental observation indicates that CHARMM36m achieved the best performance on modeling the aggregation of PHF6 peptides. In summary, CHARMM36m is currently the most suitable force field for studying the aggregation of both amyloid-β and Tau through MD simulations.
Collapse
Affiliation(s)
- Viet Hoang Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Jie Gao
- Department of Neuroscience, The Ohio State University Wexner Medical Center, Columbus, Ohio 43210, United States
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
8
|
Atazanavir, Alone or in Combination with Ritonavir, Inhibits SARS-CoV-2 Replication and Proinflammatory Cytokine Production. Antimicrob Agents Chemother 2020; 64:AAC.00825-20. [PMID: 32759267 PMCID: PMC7508582 DOI: 10.1128/aac.00825-20] [Citation(s) in RCA: 80] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2020] [Accepted: 08/02/2020] [Indexed: 12/15/2022] Open
Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is already responsible for far more deaths than previous pathogenic coronaviruses (CoVs) from 2002 and 2012. The identification of clinically approved drugs to be repurposed to combat 2019 CoV disease (COVID-19) would allow the rapid implementation of potentially life-saving procedures. The major protease (Mpro) of SARS-CoV-2 is considered a promising target, based on previous results from related CoVs with lopinavir (LPV), an HIV protease inhibitor. However, limited evidence exists for other clinically approved antiretroviral protease inhibitors. Extensive use of atazanavir (ATV) as antiretroviral and previous evidence suggesting its bioavailability within the respiratory tract prompted us to study this molecule against SARS-CoV-2. Our results show that ATV docks in the active site of SARS-CoV-2 Mpro with greater strength than LPV, blocking Mpro activity. We confirmed that ATV inhibits SARS-CoV-2 replication, alone or in combination with ritonavir (RTV) in Vero cells and a human pulmonary epithelial cell line. ATV/RTV also impaired virus-induced enhancement of interleukin 6 (IL-6) and tumor necrosis factor alpha (TNF-α) levels. Together, our data strongly suggest that ATV and ATV/RTV should be considered among the candidate repurposed drugs undergoing clinical trials in the fight against COVID-19.
Collapse
|
9
|
Nguyen PH, Sterpone F, Derreumaux P. Aggregation of disease-related peptides. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:435-460. [PMID: 32145950 DOI: 10.1016/bs.pmbts.2019.12.002] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Protein misfolding and aggregation of amyloid proteins is the fundamental cause of more than 20 diseases. Molecular mechanisms of the self-assembly and the formation of the toxic aggregates are still elusive. Computer simulations have been intensively used to study the aggregation of amyloid peptides of various amino acid lengths related to neurodegenerative diseases. We review atomistic and coarse-grained simulations of short amyloid peptides aimed at determining their transient oligomeric structures and the early and late aggregation steps.
Collapse
Affiliation(s)
- Phuong H Nguyen
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | - Fabio Sterpone
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | - Philippe Derreumaux
- Laboratory of Theoretical Chemistry, Ton Duc Thang University, Ho Chi Minh City, Vietnam; Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam.
| |
Collapse
|
10
|
Calero-Rubio C, Saluja A, Sahin E, Roberts CJ. Predicting High-Concentration Interactions of Monoclonal Antibody Solutions: Comparison of Theoretical Approaches for Strongly Attractive Versus Repulsive Conditions. J Phys Chem B 2019; 123:5709-5720. [PMID: 31241333 DOI: 10.1021/acs.jpcb.9b03779] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Nonspecific protein-protein interactions of a monoclonal antibody were quantified experimentally using light scattering from low to high protein concentrations (c2) and compared with prior work for a different antibody that yielded qualitatively different behavior. The c2 dependence of the excess Rayleigh ratio (Rex) provided the osmotic second virial coefficient (B22) at low c2 and the static structure factor (Sq=0) at high c2, as a function of solution pH, total ionic strength (TIS), and sucrose concentration. Net repulsive interactions were observed at pH 5, with weaker repulsions at higher TIS. Conversely, attractive electrostatic interactions were observed at pH 6.5, with weaker attractions at higher TIS. Refined coarse-grained models were used to fit model parameters using experimental B22 versus TIS data. The parameters were used to predict high-c2 Rex values via Monte Carlo simulations and separately with Mayer-sampling calculations of higher-order virial coefficients. For both methods, predictions for repulsive to mildly attractive conditions were quantitatively accurate. However, only qualitatively accurate predictions were practical for strongly attractive conditions. An alternative, higher resolution model was used to show semiquantitatively and quantitatively accurate predictions of strong electrostatic attractions at low c2 and low ionic strength.
Collapse
Affiliation(s)
- Cesar Calero-Rubio
- Department of Chemical and Biomolecular Engineering , University of Delaware , Newark , Delaware 19716 , United States
| | - Atul Saluja
- Drug Product Science and Technology , Bristol-Myers Squibb , New Brunswick , New Jersey 08901 , United States
| | - Erinc Sahin
- Drug Product Science and Technology , Bristol-Myers Squibb , New Brunswick , New Jersey 08901 , United States
| | - Christopher J Roberts
- Department of Chemical and Biomolecular Engineering , University of Delaware , Newark , Delaware 19716 , United States
| |
Collapse
|
11
|
Pasquali S, Frezza E, Barroso da Silva FL. Coarse-grained dynamic RNA titration simulations. Interface Focus 2019; 9:20180066. [PMID: 31065339 DOI: 10.1098/rsfs.2018.0066] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/13/2019] [Indexed: 12/11/2022] Open
Abstract
Electrostatic interactions play a pivotal role in many biomolecular processes. The molecular organization and function in biological systems are largely determined by these interactions. Owing to the highly negative charge of RNA, the effect is expected to be more pronounced in this system. Moreover, RNA base pairing is dependent on the charge of the base, giving rise to alternative secondary and tertiary structures. The equilibrium between uncharged and charged bases is regulated by the solution pH, which is therefore a key environmental condition influencing the molecule's structure and behaviour. By means of constant-pH Monte Carlo simulations based on a fast proton titration scheme, coupled with the coarse-grained model HiRE-RNA, molecular dynamic simulations of RNA molecules at constant pH enable us to explore the RNA conformational plasticity at different pH values as well as to compute electrostatic properties as local pK a values for each nucleotide.
Collapse
Affiliation(s)
- S Pasquali
- Laboratoire de Cristallographie et RMN Biologiques, CNRS UMR 8015, Université Paris Descartes, Paris 75006, France
| | - E Frezza
- Laboratoire de Cristallographie et RMN Biologiques, CNRS UMR 8015, Université Paris Descartes, Paris 75006, France
| | - F L Barroso da Silva
- Departamento de Física e Química, Faculdade de Ciência s Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Av. do café, s/no, Ribeirão Preto, SP BR-14040-903, Brazil.,Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, USA
| |
Collapse
|
12
|
Garate JA, Bernardin A, Escalona Y, Yanez C, English NJ, Perez-Acle T. Orientational and Folding Thermodynamics via Electric Dipole Moment Restraining. J Phys Chem B 2019; 123:2599-2608. [PMID: 30831028 DOI: 10.1021/acs.jpcb.8b09374] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The projection of molecular processes onto a small set of relevant descriptors, the so-called reaction coordinates or collective variables (CVs), is a technique nowadays routinely employed by the biomolecular simulation community. In this work, we implemented two CVs to manipulate the orientation (i.e., angle) (μ⃗a) and magnitude (|μ⃗|) of the electric dipole moment. In doing so, we studied the thermodynamics of water orientation under the application of external voltages and the folding of two polypeptides at zero-field conditions. The projection of the free-energy [potential of mean force (PMF)] along water orientation defined an upper limit of around 0.3 V for irrelevant thermodynamic effects. On the other hand, sufficiently strong μ⃗a restraints applied on 12-alanine (Ala12) triggered structural effects because of the alignment of local dipoles; for lower restraints, a full-body rotation is achieved. The manipulation of |μ⃗| produced strong perturbations on the secondary structure of Ala12, promoting an enhanced sampling to its configurational space. Rigorous free-energy calculations in the form of 2-D PMFs for deca-alanine showed the utility of |μ⃗| as a reaction coordinate to study folding in small α helices. As a whole, we propose that the manipulation of both components of the dipole moment, μ⃗a and |μ⃗|, provides thermodynamics insights into the structural conformation and stability of biomolecules. These new CVs are implemented in the Colvars module, available for NAMD and LAMMPS.
Collapse
Affiliation(s)
- Jose Antonio Garate
- Centro Interdisciplinario de Neurociencia de Valparaiso , Universidad de Valparaiso , Pasaje Harrington 287 , Playa Ancha, Valparaiso 2381850 , Chile
| | - Alejandro Bernardin
- Centro Interdisciplinario de Neurociencia de Valparaiso , Universidad de Valparaiso , Pasaje Harrington 287 , Playa Ancha, Valparaiso 2381850 , Chile.,Computational Biology Lab , Fundacion Ciencia & Vida , Avenida Zanartu 1482, Nunoa , Santiago 7780272 , Chile
| | - Yerko Escalona
- Institute for Molecular Modeling and Simulation , Muthgasse 18 , Vienna 1190 , Austria
| | - Carlos Yanez
- Computational Biology Lab , Fundacion Ciencia & Vida , Avenida Zanartu 1482, Nunoa , Santiago 7780272 , Chile
| | - Niall J English
- School of Chemical and Bioprocess Engineering , University College Dublin , Belfield, Dublin 4 , Ireland
| | - Tomas Perez-Acle
- Centro Interdisciplinario de Neurociencia de Valparaiso , Universidad de Valparaiso , Pasaje Harrington 287 , Playa Ancha, Valparaiso 2381850 , Chile.,Computational Biology Lab , Fundacion Ciencia & Vida , Avenida Zanartu 1482, Nunoa , Santiago 7780272 , Chile
| |
Collapse
|
13
|
Lu Y, Shi XF, Nguyen PH, Sterpone F, Salsbury FR, Derreumaux P. Amyloid-β(29-42) Dimeric Conformations in Membranes Rich in Omega-3 and Omega-6 Polyunsaturated Fatty Acids. J Phys Chem B 2019; 123:2687-2696. [PMID: 30813725 DOI: 10.1021/acs.jpcb.9b00431] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The omega-3 and omega-6 polyunsaturated fatty acids are two important components of cell membranes in human brains. When incorporated into phospholipids, omega-3 slows the progression of Alzheimer's disease (AD), whereas omega-6 is linked to increased risk of AD. Little is known on the amyloid-β (Aβ) conformations in membranes rich in omega-3 and omega-6 phospholipids. Herein, the structural properties of the Aβ29-42 dimer embedded in both fatty acid membranes were comparatively studied to a 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine (POPC) bilayer using all-atom molecular dynamics (MD) simulations. Starting from α-helix, both omega-6 and omega-3 membranes promote new orientations and conformations of the dimer, in agreement with the observed dependence of Aβ production upon addition of these two fatty acids. This conformational result is corroborated by atomistic MD simulations of the dimer of the 99 amino acid C-terminal fragment of amyloid precursor protein spanning the residues 15-55. Starting from β-sheet, omega-6 membrane promotes helical and disordered structures of Aβ29-42 dimer, whereas omega-3 membrane preserves the β-sheet structures differing however from those observed in POPC. Remarkably, the mixture of the two fatty acids and POPC depicts another conformational ensemble of the Aβ29-42 dimer. This finding demonstrates that variation in the abundance of the molecular phospholipids, which changes with age, modulates membrane-embedded Aβ oligomerization.
Collapse
Affiliation(s)
- Yan Lu
- School of Physics and Optoelectronic Engineering , Xidian University , Xi'an 710071 , China
| | - Xiao-Feng Shi
- School of Physics and Optoelectronic Engineering , Xidian University , Xi'an 710071 , China
| | - Phuong H Nguyen
- Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique (IBPC), UPR9080 CNRS , Université Paris Diderot, Sorbonne Paris Cite , 13 rue Pierre et Marie Curie , 75005 Paris , France
| | - Fabio Sterpone
- Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique (IBPC), UPR9080 CNRS , Université Paris Diderot, Sorbonne Paris Cite , 13 rue Pierre et Marie Curie , 75005 Paris , France
| | - Freddie R Salsbury
- Department of Physics , Wake Forest University , Winston-Salem , North Carolina 27106 , United States
| | - Philippe Derreumaux
- Laboratory of Theoretical Chemistry , Ton Duc Thang University , Ho Chi Minh City , Vietnam.,Faculty of Pharmacy , Ton Duc Thang University , Ho Chi Minh City , Vietnam
| |
Collapse
|
14
|
Molecular simulation of peptides coming of age: Accurate prediction of folding, dynamics and structures. Arch Biochem Biophys 2019; 664:76-88. [DOI: 10.1016/j.abb.2019.01.033] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2018] [Revised: 01/23/2019] [Accepted: 01/28/2019] [Indexed: 12/24/2022]
|
15
|
Man VH, He X, Derreumaux P, Ji B, Xie XQ, Nguyen PH, Wang J. Effects of All-Atom Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of Aβ 16-22 Dimer. J Chem Theory Comput 2019; 15:1440-1452. [PMID: 30633867 PMCID: PMC6745714 DOI: 10.1021/acs.jctc.8b01107] [Citation(s) in RCA: 87] [Impact Index Per Article: 17.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
We investigated the effects of 17 widely used atomistic molecular mechanics force fields (MMFFs) on the structures and kinetics of amyloid peptide assembly. To this end, we performed large-scale all-atom molecular dynamics simulations in explicit water on the dimer of the seven-residue fragment of the Alzheimer's amyloid-β peptide, Aβ16-22, for a total time of 0.34 ms. We compared the effects of these MMFFs by analyzing various global reaction coordinates, secondary structure contents, the fibril population, the in-register and out-of-register architectures, and the fibril formation time at 310 K. While the AMBER94, AMBER99, and AMBER12SB force fields do not predict any β-sheets, the seven force fields, AMBER96, GROMOS45a3, GROMOS53a5, GROMOS53a6, GROMOS43a1, GROMOS43a2, and GROMOS54a7, form β-sheets rapidly. In contrast, the following five force fields, AMBER99-ILDN, AMBER14SB, CHARMM22*, CHARMM36, and CHARMM36m, are the best candidates for studying amyloid peptide assembly, as they provide good balances in terms of structures and kinetics. We also investigated the assembly mechanisms of dimeric Aβ16-22 and found that the fibril formation rate is predominantly controlled by the total β-strand content.
Collapse
Affiliation(s)
- Viet Hoang Man
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Xibing He
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique UPR 9080, CNRS, Université Denis Diderot, Sorbonne Paris Cité, IBPC, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Beihong Ji
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Xiang-Qun Xie
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Phuong H. Nguyen
- Laboratoire de Biochimie Théorique UPR 9080, CNRS, Université Denis Diderot, Sorbonne Paris Cité, IBPC, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Junmei Wang
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
- Corresponding Author:
| |
Collapse
|
16
|
Li X, Lei J, Qi R, Xie L, Wei G. Mechanistic insight into E22Q-mutation-induced antiparallel-to-parallel β-sheet transition of Aβ16−22fibrils: an all-atom simulation study. Phys Chem Chem Phys 2019; 21:15686-15694. [DOI: 10.1039/c9cp02561h] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
E22Q mutation of Aβ16−22fibrils facilitates parallel β-sheet formation by enhancing Q22–Q22 hydrogen-bonding interaction and A21–A21, F20–F20, F19–F19 and V18–V18 hydrophobic interaction.
Collapse
Affiliation(s)
- Xuhua Li
- State Key Laboratory of Surface Physics
- Key Laboratory for Computational Physical Sciences (Ministry of Education), and Department of Physics
- Fudan University
- Shanghai
- China
| | - Jiangtao Lei
- State Key Laboratory of Surface Physics
- Key Laboratory for Computational Physical Sciences (Ministry of Education), and Department of Physics
- Fudan University
- Shanghai
- China
| | - Ruxi Qi
- State Key Laboratory of Surface Physics
- Key Laboratory for Computational Physical Sciences (Ministry of Education), and Department of Physics
- Fudan University
- Shanghai
- China
| | - Luogang Xie
- College of Physics and Electronic Engineering
- Zhengzhou University of Light Industry
- Zhengzhou 453002
- People's Republic of China
| | - Guanghong Wei
- State Key Laboratory of Surface Physics
- Key Laboratory for Computational Physical Sciences (Ministry of Education), and Department of Physics
- Fudan University
- Shanghai
- China
| |
Collapse
|
17
|
Passarini I, Rossiter S, Malkinson J, Zloh M. In Silico Structural Evaluation of Short Cationic Antimicrobial Peptides. Pharmaceutics 2018; 10:E72. [PMID: 29933540 PMCID: PMC6160961 DOI: 10.3390/pharmaceutics10030072] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2018] [Revised: 06/11/2018] [Accepted: 06/12/2018] [Indexed: 12/26/2022] Open
Abstract
Cationic peptides with antimicrobial properties are ubiquitous in nature and have been studied for many years in an attempt to design novel antibiotics. However, very few molecules are used in the clinic so far, sometimes due to their complexity but, mostly, as a consequence of the unfavorable pharmacokinetic profile associated with peptides. The aim of this work is to investigate cationic peptides in order to identify common structural features which could be useful for the design of small peptides or peptido-mimetics with improved drug-like properties and activity against Gram negative bacteria. Two sets of cationic peptides (AMPs) with known antimicrobial activity have been investigated. The first reference set comprised molecules with experimentally-known conformations available in the protein databank (PDB), and the second one was composed of short peptides active against Gram negative bacteria but with no significant structural information available. The predicted structures of the peptides from the first set were in excellent agreement with those experimentally-observed, which allowed analysis of the structural features of the second group using computationally-derived conformations. The peptide conformations, either experimentally available or predicted, were clustered in an “all vs. all” fashion and the most populated clusters were then analyzed. It was confirmed that these peptides tend to assume an amphipathic conformation regardless of the environment. It was also observed that positively-charged amino acid residues can often be found next to aromatic residues. Finally, a protocol was evaluated for the investigation of the behavior of short cationic peptides in the presence of a membrane-like environment such as dodecylphosphocholine (DPC) micelles. The results presented herein introduce a promising approach to inform the design of novel short peptides with a potential antimicrobial activity.
Collapse
Affiliation(s)
- Ilaria Passarini
- School of Life and Medical Sciences, University of Hertfordshire, College Lane, Hatfield AL10 9AB, UK.
| | - Sharon Rossiter
- School of Life and Medical Sciences, University of Hertfordshire, College Lane, Hatfield AL10 9AB, UK.
| | - John Malkinson
- UCL School of Pharmacy, University College London, 29/39 Brunswick Square, London WC1N 1AX, UK.
| | - Mire Zloh
- School of Life and Medical Sciences, University of Hertfordshire, College Lane, Hatfield AL10 9AB, UK.
- Faculty of Pharmacy, University Business Academy, Trg mladenaca 5, 21000 Novi Sad, Serbia.
- NanoPuzzle Medicines Design, Business & Technology Centre, Bessemer Drive, Stevenage SG1 2DX, UK.
| |
Collapse
|
18
|
Nguyen PH, del Castillo-Frias MP, Berthoumieux O, Faller P, Doig AJ, Derreumaux P. Amyloid-β/Drug Interactions from Computer Simulations and Cell-Based Assays. J Alzheimers Dis 2018; 64:S659-S672. [DOI: 10.3233/jad-179902] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Phuong H. Nguyen
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, IBPC, Paris, France
| | - Maria P. del Castillo-Frias
- Manchester Institute of Biotechnology and Department of Chemistry, Faculty of Science and Engineering, The University of Manchester, Manchester, UK
| | - Olivia Berthoumieux
- CNRS, LCC (Laboratoire de Chimie de Coordination), Toulouse Cedex 4, France et Université de Toulouse, UPS, INPT, Toulouse Cedex 4, France
| | - Peter Faller
- Biometals and Biology Chemistry, Institut de Chimie (CNRS UMR7177), Université de Strasbourg, Strasbourg, France
| | - Andrew J. Doig
- Manchester Institute of Biotechnology and Department of Chemistry, Faculty of Science and Engineering, The University of Manchester, Manchester, UK
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, IBPC, Paris, France
| |
Collapse
|
19
|
Liu X, Shi D, Zhou S, Liu H, Liu H, Yao X. Molecular dynamics simulations and novel drug discovery. Expert Opin Drug Discov 2017; 13:23-37. [DOI: 10.1080/17460441.2018.1403419] [Citation(s) in RCA: 129] [Impact Index Per Article: 18.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Xuewei Liu
- State Key Laboratory of Applied Organic Chemistry and Department of Chemistry, Lanzhou University, Lanzhou, China
| | - Danfeng Shi
- State Key Laboratory of Applied Organic Chemistry and Department of Chemistry, Lanzhou University, Lanzhou, China
| | | | - Hongli Liu
- School of Pharmacy, Lanzhou University, Lanzhou, China
| | - Huanxiang Liu
- School of Pharmacy, Lanzhou University, Lanzhou, China
| | - Xiaojun Yao
- State Key Laboratory of Applied Organic Chemistry and Department of Chemistry, Lanzhou University, Lanzhou, China
| |
Collapse
|
20
|
Zhang C, Drake JA, Ma J, Pettitt BM. Optimal updating magnitude in adaptive flat-distribution sampling. J Chem Phys 2017; 147:174105. [PMID: 29117700 PMCID: PMC5669982 DOI: 10.1063/1.5008618] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2017] [Accepted: 10/08/2017] [Indexed: 11/14/2022] Open
Abstract
We present a study on the optimization of the updating magnitude for a class of free energy methods based on flat-distribution sampling, including the Wang-Landau (WL) algorithm and metadynamics. These methods rely on adaptive construction of a bias potential that offsets the potential of mean force by histogram-based updates. The convergence of the bias potential can be improved by decreasing the updating magnitude with an optimal schedule. We show that while the asymptotically optimal schedule for the single-bin updating scheme (commonly used in the WL algorithm) is given by the known inverse-time formula, that for the Gaussian updating scheme (commonly used in metadynamics) is often more complex. We further show that the single-bin updating scheme is optimal for very long simulations, and it can be generalized to a class of bandpass updating schemes that are similarly optimal. These bandpass updating schemes target only a few long-range distribution modes and their optimal schedule is also given by the inverse-time formula. Constructed from orthogonal polynomials, the bandpass updating schemes generalize the WL and Langfeld-Lucini-Rago algorithms as an automatic parameter tuning scheme for umbrella sampling.
Collapse
Affiliation(s)
- Cheng Zhang
- Sealy Center for Structural Biology and Molecular Biophysics, The University of Texas Medical Branch, Galveston, Texas 77555-0304, USA
| | - Justin A Drake
- Sealy Center for Structural Biology and Molecular Biophysics, The University of Texas Medical Branch, Galveston, Texas 77555-0304, USA
| | - Jianpeng Ma
- Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, Texas 77030, USA
| | - B Montgomery Pettitt
- Sealy Center for Structural Biology and Molecular Biophysics, The University of Texas Medical Branch, Galveston, Texas 77555-0304, USA
| |
Collapse
|
21
|
Palazzesi F, Valsson O, Parrinello M. Conformational Entropy as Collective Variable for Proteins. J Phys Chem Lett 2017; 8:4752-4756. [PMID: 28906117 DOI: 10.1021/acs.jpclett.7b01770] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Many enhanced sampling methods rely on the identification of appropriate collective variables. For proteins, even small ones, finding appropriate descriptors has proven challenging. Here we suggest that the NMR S2 order parameter can be used to this effect. We trace the validity of this statement to the suggested relation between S2 and conformational entropy. Using the S2 order parameter and a surrogate for the protein enthalpy in conjunction with metadynamics or variationally enhanced sampling, we are able to reversibly fold and unfold a small protein and draw its free energy at a fraction of the time that is needed in unbiased simulations. We also use S2 in combination with the free energy flooding method to compute the unfolding rate of this peptide. We repeat this calculation at different temperatures to obtain the unfolding activation energy.
Collapse
Affiliation(s)
- Ferruccio Palazzesi
- Department of Chemistry and Applied Biosciences, ETH Zurich c/o USI Campus , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
- Facoltà di Informatica, Instituto di Scienze Computationali, Università della Svizzera italiana , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
| | - Omar Valsson
- Department of Chemistry and Applied Biosciences, ETH Zurich c/o USI Campus , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
- Facoltà di Informatica, Instituto di Scienze Computationali, Università della Svizzera italiana , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
- National Center for Computational Design and Discovery of Novel Materials MARVEL, Università della Svizzera italiana , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich c/o USI Campus , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
- Facoltà di Informatica, Instituto di Scienze Computationali, Università della Svizzera italiana , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
- National Center for Computational Design and Discovery of Novel Materials MARVEL, Università della Svizzera italiana , Via Giuseppe Buffi 13, CH-6900, Lugano, Switzerland
| |
Collapse
|
22
|
Multi-scale simulations of biological systems using the OPEP coarse-grained model. Biochem Biophys Res Commun 2017; 498:296-304. [PMID: 28917842 DOI: 10.1016/j.bbrc.2017.08.165] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2017] [Accepted: 08/31/2017] [Indexed: 12/14/2022]
Abstract
Biomolecules are complex machines that are optimized by evolution to properly fulfill or contribute to a variety of biochemical tasks in the cellular environment. Computer simulations based on quantum mechanics and atomistic force fields have been proven to be a powerful microscope for obtaining valuable insights into many biological, physical, and chemical processes. Many interesting phenomena involve, however, a time scale and a number of degrees of freedom, notably if crowding is considered, that cannot be explored at an atomistic resolution. To bridge the gap between reality and simulation, many different advanced computational techniques and coarse-grained (CG) models have been developed. Here, we report some applications of the CG OPEP protein model to amyloid fibril formation, the response of catch-bond proteins to two types of fluid flow, and interactive simulations to fold peptides with well-defined 3D structures or with intrinsic disorder.
Collapse
|
23
|
Chiricotto M, Melchionna S, Derreumaux P, Sterpone F. Hydrodynamic effects on β-amyloid (16-22) peptide aggregation. J Chem Phys 2017; 145:035102. [PMID: 27448906 DOI: 10.1063/1.4958323] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Computer simulations based on simplified representations are routinely used to explore the early steps of amyloid aggregation. However, when protein models with implicit solvent are employed, these simulations miss the effect of solvent induced correlations on the aggregation kinetics and lifetimes of metastable states. In this work, we apply the multi-scale Lattice Boltzmann Molecular Dynamics technique (LBMD) to investigate the initial aggregation phases of the amyloid Aβ16-22 peptide. LBMD includes naturally hydrodynamic interactions (HIs) via a kinetic on-lattice representation of the fluid kinetics. The peptides are represented by the flexible OPEP coarse-grained force field. First, we have tuned the essential parameters that control the coupling between the molecular and fluid evolutions in order to reproduce the experimental diffusivity of elementary species. The method is then deployed to investigate the effect of HIs on the aggregation of 100 and 1000 Aβ16-22 peptides. We show that HIs clearly impact the aggregation process and the fluctuations of the oligomer sizes by favouring the fusion and exchange dynamics of oligomers between aggregates. HIs also guide the growth of the leading largest cluster. For the 100 Aβ16-22 peptide system, the simulation of ∼300 ns allowed us to observe the transition from ellipsoidal assemblies to an elongated and slightly twisted aggregate involving almost the totality of the peptides. For the 1000 Aβ16-22 peptides, a system of unprecedented size at quasi-atomistic resolution, we were able to explore a branched disordered fibril-like structure that has never been described by other computer simulations, but has been observed experimentally.
Collapse
Affiliation(s)
- Mara Chiricotto
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, University Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Simone Melchionna
- CNR-ISC, Institute for Complex System, Consiglio Nazionale delle Ricerche, Rome, Italy
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, University Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Fabio Sterpone
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, University Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
24
|
In silico structural characterization of protein targets for drug development against Trypanosoma cruzi. J Mol Model 2016; 22:244. [DOI: 10.1007/s00894-016-3115-9] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2016] [Accepted: 09/02/2016] [Indexed: 10/21/2022]
|
25
|
Kührová P, Best RB, Bottaro S, Bussi G, Šponer J, Otyepka M, Banáš P. Computer Folding of RNA Tetraloops: Identification of Key Force Field Deficiencies. J Chem Theory Comput 2016; 12:4534-48. [PMID: 27438572 PMCID: PMC6169534 DOI: 10.1021/acs.jctc.6b00300] [Citation(s) in RCA: 91] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The computer-aided folding of biomolecules, particularly RNAs, is one of the most difficult challenges in computational structural biology. RNA tetraloops are fundamental RNA motifs playing key roles in RNA folding and RNA-RNA and RNA-protein interactions. Although state-of-the-art Molecular Dynamics (MD) force fields correctly describe the native state of these tetraloops as a stable free-energy basin on the microsecond time scale, enhanced sampling techniques reveal that the native state is not the global free energy minimum, suggesting yet unidentified significant imbalances in the force fields. Here, we tested our ability to fold the RNA tetraloops in various force fields and simulation settings. We employed three different enhanced sampling techniques, namely, temperature replica exchange MD (T-REMD), replica exchange with solute tempering (REST2), and well-tempered metadynamics (WT-MetaD). We aimed to separate problems caused by limited sampling from those due to force-field inaccuracies. We found that none of the contemporary force fields is able to correctly describe folding of the 5'-GAGA-3' tetraloop over a range of simulation conditions. We thus aimed to identify which terms of the force field are responsible for this poor description of TL folding. We showed that at least two different imbalances contribute to this behavior, namely, overstabilization of base-phosphate and/or sugar-phosphate interactions and underestimated stability of the hydrogen bonding interaction in base pairing. The first artifact stabilizes the unfolded ensemble, while the second one destabilizes the folded state. The former problem might be partially alleviated by reparametrization of the van der Waals parameters of the phosphate oxygens suggested by Case et al., while in order to overcome the latter effect we suggest local potentials to better capture hydrogen bonding interactions.
Collapse
Affiliation(s)
- Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic,
| | - Robert B. Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892-0520
| | - Sandro Bottaro
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, 34136 Trieste, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, 34136 Trieste, Italy
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic,
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic,
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic,
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| |
Collapse
|
26
|
Chiricotto M, Tran TT, Nguyen PH, Melchionna S, Sterpone F, Derreumaux P. Coarse-grained and All-atom Simulations towards the Early and Late Steps of Amyloid Fibril Formation. Isr J Chem 2016. [DOI: 10.1002/ijch.201600048] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Mara Chiricotto
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS; Université Paris Diderot, Sorbonne Paris Cité, IBPC; 13 Rue Pierre et Marie Curie 75005 Paris France
| | - Thanh Thuy Tran
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS; Université Paris Diderot, Sorbonne Paris Cité, IBPC; 13 Rue Pierre et Marie Curie 75005 Paris France
| | - Phuong H. Nguyen
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS; Université Paris Diderot, Sorbonne Paris Cité, IBPC; 13 Rue Pierre et Marie Curie 75005 Paris France
| | - Simone Melchionna
- Istituto Sistemi Complessi; Consiglio Nazionale delle Ricerche; P. le A. Moro 2 00185 Rome Italy
| | - Fabio Sterpone
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS; Université Paris Diderot, Sorbonne Paris Cité, IBPC; 13 Rue Pierre et Marie Curie 75005 Paris France
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS; Université Paris Diderot, Sorbonne Paris Cité, IBPC; 13 Rue Pierre et Marie Curie 75005 Paris France
| |
Collapse
|
27
|
Calero-Rubio C, Paik B, Jia X, Kiick KL, Roberts CJ. Predicting unfolding thermodynamics and stable intermediates for alanine-rich helical peptides with the aid of coarse-grained molecular simulation. Biophys Chem 2016; 217:8-19. [PMID: 27486699 DOI: 10.1016/j.bpc.2016.07.002] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2016] [Revised: 07/14/2016] [Accepted: 07/17/2016] [Indexed: 10/21/2022]
Abstract
This report focuses on the molecular-level processes and thermodynamics of unfolding of a series of helical peptides using a coarse-grained (CG) molecular model. The CG model was refined to capture thermodynamics and structural changes as a function of temperature for a set of published peptide sequences. Circular dichroism spectroscopy (CD) was used to experimentally monitor the temperature-dependent conformational changes and stability of published peptides and new sequences introduced here. The model predictions were quantitatively or semi-quantitatively accurate in all cases. The simulations and CD results showed that, as expected, in most cases the unfolding of helical peptides is well described by a simply 2-state model, and conformational stability increased with increased length of the helices. A notable exception in a 19-residue helix was when two Ala residues were each replaced with Phe. This stabilized a partly unfolded intermediate state via hydrophobic contacts, and also promoted aggregates at higher peptide concentrations.
Collapse
Affiliation(s)
- Cesar Calero-Rubio
- Chemical & Biomolecular Engineering Department, University of Delaware, Newark, DE 19716, United States
| | - Bradford Paik
- Material Science & Engineering Department, University of Delaware, Newark, DE 19716, United States
| | - Xinqiao Jia
- Material Science & Engineering Department, University of Delaware, Newark, DE 19716, United States
| | - Kristi L Kiick
- Material Science & Engineering Department, University of Delaware, Newark, DE 19716, United States.
| | - Christopher J Roberts
- Chemical & Biomolecular Engineering Department, University of Delaware, Newark, DE 19716, United States.
| |
Collapse
|
28
|
Lamiable A, Thevenet P, Tufféry P. A critical assessment of hidden markov model sub-optimal sampling strategies applied to the generation of peptide 3D models. J Comput Chem 2016; 37:2006-16. [PMID: 27317417 DOI: 10.1002/jcc.24422] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2016] [Revised: 05/03/2016] [Accepted: 05/17/2016] [Indexed: 12/23/2022]
Abstract
Hidden Markov Model derived structural alphabets are a probabilistic framework in which the complete conformational space of a peptidic chain is described in terms of probability distributions that can be sampled to identify conformations of largest probabilities. Here, we assess how three strategies to sample sub-optimal conformations-Viterbi k-best, forward backtrack and a taboo sampling approach-can lead to the efficient generation of peptide conformations. We show that the diversity of sampling is essential to compensate biases introduced in the estimates of the probabilities, and we find that only the forward backtrack and a taboo sampling strategies can efficiently generate native or near-native models. Finally, we also find such approaches are as efficient as former protocols, while being one order of magnitude faster, opening the door to the large scale de novo modeling of peptides and mini-proteins. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- A Lamiable
- INSERM UMR-S 973, Université Paris Diderot, Sorbonne Paris Cité
| | - P Thevenet
- INSERM UMR-S 973, Université Paris Diderot, Sorbonne Paris Cité
| | - P Tufféry
- INSERM UMR-S 973, Université Paris Diderot, Sorbonne Paris Cité
| |
Collapse
|
29
|
Nguyen PH, Sterpone F, Campanera JM, Nasica-Labouze J, Derreumaux P. Impact of the A2V Mutation on the Heterozygous and Homozygous Aβ1-40 Dimer Structures from Atomistic Simulations. ACS Chem Neurosci 2016; 7:823-32. [PMID: 27007027 DOI: 10.1021/acschemneuro.6b00053] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
The A2V mutation was reported to protect from Alzheimer's disease in its heterozygous form and cause an early Alzheimer's disease type dementia in its homozygous form. Experiments showed that the aggregation rate follows the order A2V > WT (wild-type) > A2V-WT. To understand the impact of this mutation, we carried out replica exchange molecular dynamics simulations of Aβ1-40 WT-A2V and A2V-A2V dimers and compared to the WT dimer. Our atomistic simulations reveal that the mean secondary structure remains constant, but there are substantial differences in the intramolecular and intermolecular conformations upon single and double A2V mutation. Upon single mutation, the intrinsic disorder is reduced, the intermolecular potential energies are reduced, the population of intramolecular three-stranded β-sheets is increased, and the number of all α dimer topologies is decreased. Taken together, these results offer an explanation for the reduced aggregation rate of the Aβ1-40 A2V-WT peptides and the protective effect of A2V in heterozygotes.
Collapse
Affiliation(s)
- Phuong H. Nguyen
- Laboratoire de
Biochimie Théorique, UPR 9080 CNRS, Université Paris
Diderot, Sorbonne Paris Cité, IBPC, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Fabio Sterpone
- Laboratoire de
Biochimie Théorique, UPR 9080 CNRS, Université Paris
Diderot, Sorbonne Paris Cité, IBPC, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Josep M. Campanera
- Departament
de Fisicoquimica, Facultat de Farmacia, Universitat de Barcelona, 08028 Barcelona, Catalonia, Spain
| | - Jessica Nasica-Labouze
- Laboratoire de
Biochimie Théorique, UPR 9080 CNRS, Université Paris
Diderot, Sorbonne Paris Cité, IBPC, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Philippe Derreumaux
- Laboratoire de
Biochimie Théorique, UPR 9080 CNRS, Université Paris
Diderot, Sorbonne Paris Cité, IBPC, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
30
|
Lei J, Qi R, Wei G, Nussinov R, Ma B. Self-aggregation and coaggregation of the p53 core fragment with its aggregation gatekeeper variant. Phys Chem Chem Phys 2016; 18:8098-107. [PMID: 26923710 PMCID: PMC6456058 DOI: 10.1039/c5cp06538k] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
Recent studies suggested that p53 aggregation can lead to loss-of-function (LoF), dominant-negative (DN) and gain-of-function (GoF) effects, with adverse cancer consequences. The p53 aggregation-nucleating (251)ILTIITL(257) fragment is a key segment in wild-type p53 aggregation; however, an I254R mutation can prevent it. It was suggested that self-assembly of wild-type p53 and its cross-interaction with mutants differ from the classical amyloid nucleation-growth mechanism. Here, using replica exchange molecular dynamics (REMD) simulations, we studied the cross-interactions of this p53 core fragment and its aggregation rescue I254R mutant. We found that the core fragment displays strong aggregation propensity, whereas the gatekeeper I254R mutant tends to be disordered, consistent with experiments. Our cross-interaction results reveal that the wild-type p53 fragment promotes β-sheet formation of the I254R mutant by shifting the disordered mutant peptides into aggregating states. As a result, the system has similar oligomeric structures, inter-peptide interactions and free energy landscape as the wild type fragment does, revealing a prion-like process. We also found that in the cross-interaction system, the wild-type species has higher tendency to interact with the mutant than with itself. This phenomenon illustrates synergistic effects between the p53 (251)ILTIITL(257) fragment and the mutant resembling prion cross-species propagation, cautioning against exploiting it in drug discovery.
Collapse
Affiliation(s)
- Jiangtao Lei
- State Key Laboratory of Surface Physics, Key Laboratory for Computational Physical Sciences (MOE), and Department of Physics, Fudan University, Shanghai, P. R. China.
| | - Ruxi Qi
- State Key Laboratory of Surface Physics, Key Laboratory for Computational Physical Sciences (MOE), and Department of Physics, Fudan University, Shanghai, P. R. China.
| | - Guanghong Wei
- State Key Laboratory of Surface Physics, Key Laboratory for Computational Physical Sciences (MOE), and Department of Physics, Fudan University, Shanghai, P. R. China.
| | - Ruth Nussinov
- Basic Science Program, Leidos Biomedical Research, Inc., Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland 21702, USA. and Sackler Inst. of Molecular Medicine Department of Human Genetics and Molecular Medicine Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| | - Buyong Ma
- Basic Science Program, Leidos Biomedical Research, Inc., Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland 21702, USA.
| |
Collapse
|
31
|
Enhanced, targeted sampling of high-dimensional free-energy landscapes using variationally enhanced sampling, with an application to chignolin. Proc Natl Acad Sci U S A 2016; 113:1150-5. [PMID: 26787868 DOI: 10.1073/pnas.1519712113] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The capabilities of molecular simulations have been greatly extended by a number of widely used enhanced sampling methods that facilitate escaping from metastable states and crossing large barriers. Despite these developments there are still many problems which remain out of reach for these methods which has led to a vigorous effort in this area. One of the most important problems that remains unsolved is sampling high-dimensional free-energy landscapes and systems that are not easily described by a small number of collective variables. In this work we demonstrate a new way to compute free-energy landscapes of high dimensionality based on the previously introduced variationally enhanced sampling, and we apply it to the miniprotein chignolin.
Collapse
|
32
|
Stirnemann G, Sterpone F. Recovering Protein Thermal Stability Using All-Atom Hamiltonian Replica-Exchange Simulations in Explicit Solvent. J Chem Theory Comput 2015; 11:5573-7. [DOI: 10.1021/acs.jctc.5b00954] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Guillaume Stirnemann
- CNRS Laboratoire de Biochimie
Théorique, Institut de Biologie Physico-Chimique, Univ. Paris
Denis Diderot, Sorbonne Paris Cité, PSL Research University, 13 rue Pierre et Marie Curie, 75005, Paris, France
| | - Fabio Sterpone
- CNRS Laboratoire de Biochimie
Théorique, Institut de Biologie Physico-Chimique, Univ. Paris
Denis Diderot, Sorbonne Paris Cité, PSL Research University, 13 rue Pierre et Marie Curie, 75005, Paris, France
| |
Collapse
|
33
|
Tarus B, Tran TT, Nasica-Labouze J, Sterpone F, Nguyen PH, Derreumaux P. Structures of the Alzheimer's Wild-Type Aβ1-40 Dimer from Atomistic Simulations. J Phys Chem B 2015; 119:10478-87. [PMID: 26228450 DOI: 10.1021/acs.jpcb.5b05593] [Citation(s) in RCA: 74] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We have studied the dimer of amyloid beta peptide Aβ of 40 residues by means of all-atom replica exchange molecular dynamics. The Aβ-dimers have been found to be the smallest toxic species in Alzheimer's disease, but their inherent flexibilities have precluded structural characterization by experimental methods. Though the 24-μs-scale simulation reveals a mean secondary structure of 18% β-strand and 10% α helix, we find transient configurations with an unstructured N-terminus and multiple β-hairpins spanning residues 17-21 and 30-36, but the antiparallel and perpendicular peptide orientations are preferred over the parallel organization. Short-lived conformational states also consist of all α topologies, and one compact peptide with β-sheet structure stabilized by a rather extended peptide with α-helical content. Overall, this first all-atom study provides insights into the equilibrium structure of the Aβ1-40 dimer in aqueous solution, opening a new avenue for a comprehensive understanding of the impact of pathogenic and protective mutations in early-stage Alzheimer's disease on a molecular level.
Collapse
Affiliation(s)
- Bogdan Tarus
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, IBPC, Université Paris Diderot, Sorbonne Paris Cité, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Thanh T Tran
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, IBPC, Université Paris Diderot, Sorbonne Paris Cité, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Jessica Nasica-Labouze
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, IBPC, Université Paris Diderot, Sorbonne Paris Cité, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Fabio Sterpone
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, IBPC, Université Paris Diderot, Sorbonne Paris Cité, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Phuong H Nguyen
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, IBPC, Université Paris Diderot, Sorbonne Paris Cité, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, IBPC, Université Paris Diderot, Sorbonne Paris Cité, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|