1
|
Andrews B, Schweitzer-Stenner R, Urbanc B. Intrinsic Conformational Dynamics of Glycine and Alanine in Polarizable Molecular Dynamics Force Fields: Comparison to Spectroscopic Data. J Phys Chem B 2024; 128:6217-6231. [PMID: 38877893 PMCID: PMC11215781 DOI: 10.1021/acs.jpcb.4c02278] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Revised: 05/29/2024] [Accepted: 05/30/2024] [Indexed: 06/28/2024]
Abstract
Molecular dynamics (MD) is a great tool for elucidating conformational dynamics of proteins and peptides in water at the atomistic level that often surpasses the level of detail available experimentally. Structure predictions, however, are limited by the accuracy of the underlying MD force field. This limitation is particularly stark in the case of intrinsically disordered peptides and proteins, which are characterized by solvent-accessible and disordered peptide regions and domains. Recent studies show that most additive MD force fields, including CHARMM36m, do not reproduce the intrinsic conformational distributions of guest amino acid residues x in cationic GxG peptides in water in line with experimental data. Positing that a lack of polarizability in additive MD force fields may be the culprit for the reported discrepancies, we here examine the conformational dynamics of guest glycine and alanine residues in cationic GxG peptides in water using two polarizable MD force fields, CHARMM Drude and AMOEBA. Our results indicate that while AMOEBA captures the experimental data better than CHARMM Drude, neither of the two polarizable force fields offers an improvement of the Ramachandran distributions of glycine and alanine residues in cationic GGG and GAG peptides, respectively, over CHARMM36m.
Collapse
Affiliation(s)
- Brian Andrews
- Department
of Physics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| | | | - Brigita Urbanc
- Department
of Physics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
2
|
Shi C, Kaffy J, Ha-Duong T, Gallard JF, Pruvost A, Mabondzo A, Ciccone L, Ongeri S, Tonali N. Proteolytically Stable Diaza-Peptide Foldamers Mimic Helical Hot Spots of Protein-Protein Interactions and Act as Natural Chaperones. J Med Chem 2023; 66:12005-12017. [PMID: 37632446 DOI: 10.1021/acs.jmedchem.3c00611] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/28/2023]
Abstract
A novel class of peptidomimetic foldamers based on diaza-peptide units are reported. Circular dichroism, attenuated total reflection -Fourier transform infrared, NMR, and molecular dynamics studies demonstrate that unlike the natural parent nonapeptide, the specific incorporation of one diaza-peptide unit at the N-terminus allows helical folding in water, which is further reinforced by the introduction of a second unit at the C-terminus. The ability of these foldamers to resist proteolysis, to mimic the small helical hot spot of transthyretin-amyloid β (Aβ) cross-interaction, and to decrease pathological Aβ aggregation demonstrates that the introduction of diaza-peptide units is a valid approach for designing mimics or inhibitors of protein-protein interaction and other therapeutic peptidomimetics. This study also reveals that small peptide foldamers can play the same role as physiological chaperone proteins and opens a new way to design inhibitors of amyloid protein aggregation, a hallmark of more than 20 serious human diseases such as Alzheimer's disease.
Collapse
Affiliation(s)
- Chenghui Shi
- Université Paris-Saclay, CNRS, BioCIS, Bat. Henri Moissan, 17 av. des Sciences, 91400 Orsay, France
| | - Julia Kaffy
- Université Paris-Saclay, CNRS, BioCIS, Bat. Henri Moissan, 17 av. des Sciences, 91400 Orsay, France
| | - Tâp Ha-Duong
- Université Paris-Saclay, CNRS, BioCIS, Bat. Henri Moissan, 17 av. des Sciences, 91400 Orsay, France
| | - Jean-François Gallard
- Equipe Biologie et Chimie Structurales, Dept Chimie et Biologie Structurales et Analytiques, ICSN CNRS, Université Paris Saclay, 1 avenue de la terrasse, 91190 Gif sur Yvette, France
| | - Alain Pruvost
- CEA, INRAE, Département Médicaments et Technologies pour La Santé, Université Paris-Saclay, SPI 91191 Gif-sur-Yvette, France
| | - Aloise Mabondzo
- CEA, INRAE, Département Médicaments et Technologies pour La Santé, Université Paris-Saclay, SPI 91191 Gif-sur-Yvette, France
| | - Lidia Ciccone
- Department of Pharmacy, University of Pisa, 56126 Pisa, Italy
| | - Sandrine Ongeri
- Université Paris-Saclay, CNRS, BioCIS, Bat. Henri Moissan, 17 av. des Sciences, 91400 Orsay, France
| | - Nicolo Tonali
- Université Paris-Saclay, CNRS, BioCIS, Bat. Henri Moissan, 17 av. des Sciences, 91400 Orsay, France
| |
Collapse
|
3
|
Schweitzer-Stenner R, Kurbaj R, O'Neill N, Andrews B, Shah R, Urbanc B. Conformational Manifold Sampled by Two Short Linear Motif Segments Probed by Circular Dichroism, Vibrational, and Nuclear Magnetic Resonance Spectroscopy. Biochemistry 2023; 62:2571-2586. [PMID: 37595285 DOI: 10.1021/acs.biochem.3c00212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/20/2023]
Abstract
Disordered protein segments called short linear motifs (SLiM) serve as recognition sites for a variety of biological processes and act as targeting signals, modification, and ligand binding sites. While SLiMs do not adopt one of the known regular secondary structures, the conformational distribution might still reflect the structural propensities of their amino acid residues and possible interactions between them. In the past, conformational analyses of short peptides provided compelling evidence for the notion that individual residues are less conformationally flexible than locally expected for a random coil. Here, we combined various spectroscopies (NMR, IR, vibrational, and UV circular dichroism) to determine the Ramachandran plots of two SLiM motifs, i.e., GRRDSG and GRRTSG. They are two representatives of RxxS motifs that are capable of being phosphorylated by protein kinase A, an enzyme that plays a fundamental role in a variety of biological processes. Our results reveal that the nearest and non-nearest interactions between residues cause redistributions between polyproline II and β-strand basins while concomitantly stabilizing extended relative to turn-forming and helical structures. They also cause shifts in basin positions. With increasing temperature, β-strand populations become more populated at the expense of polyproline II. While molecular dynamics simulations with Amber ff14SB and CHARMM 36m force fields indicate residue-residue interactions, they do not account for the observed structural changes.
Collapse
Affiliation(s)
| | - Raghed Kurbaj
- Department of Chemistry, Drexel University, Philadelphia, PA19104Pennsylvania,United States
| | - Nichole O'Neill
- Department of Chemistry, Drexel University, Philadelphia, PA19104Pennsylvania,United States
| | - Brian Andrews
- Department of Physics, Drexel University, Philadelphia,PA19104Pennsylvania,United States
| | - Riya Shah
- Department of Physics, Drexel University, Philadelphia,PA19104Pennsylvania,United States
| | - Brigita Urbanc
- Department of Physics, Drexel University, Philadelphia,PA19104Pennsylvania,United States
| |
Collapse
|
4
|
Schweitzer-Stenner R. The relevance of short peptides for an understanding of unfolded and intrinsically disordered proteins. Phys Chem Chem Phys 2023; 25:11908-11933. [PMID: 37096579 DOI: 10.1039/d3cp00483j] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023]
Abstract
Over the last thirty years the unfolded state of proteins has attracted considerable interest owing to the discovery of intrinsically disordered proteins which perform a plethora of functions despite resembling unfolded proteins to a significant extent. Research on both, unfolded and disordered proteins has revealed that their conformational properties can deviate locally from random coil behavior. In this context results from work on short oligopeptides suggest that individual amino acid residues sample the sterically allowed fraction of the Ramachandran plot to a different extent. Alanine has been found to exhibit a peculiarity in that it has a very high propensity for adopting polyproline II like conformations. This Perspectives article reviews work on short peptides aimed at exploring the Ramachandran distributions of amino acid residues in different contexts with experimental and computational means. Based on the thus provided overview the article discussed to what extent short peptides can serve as tools for exploring unfolded and disordered proteins and as benchmarks for the development of a molecular dynamics force field.
Collapse
|
5
|
Kovács D, Bodor A. The influence of random-coil chemical shifts on the assessment of structural propensities in folded proteins and IDPs. RSC Adv 2023; 13:10182-10203. [PMID: 37006359 PMCID: PMC10065145 DOI: 10.1039/d3ra00977g] [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/13/2023] [Accepted: 03/15/2023] [Indexed: 04/03/2023] Open
Abstract
In studying secondary structural propensities of proteins by nuclear magnetic resonance (NMR) spectroscopy, secondary chemical shifts (SCSs) serve as the primary atomic scale observables. For SCS calculation, the selection of an appropriate random coil chemical shift (RCCS) dataset is a crucial step, especially when investigating intrinsically disordered proteins (IDPs). The scientific literature is abundant in such datasets, however, the effect of choosing one over all the others in a concrete application has not yet been studied thoroughly and systematically. Hereby, we review the available RCCS prediction methods and to compare them, we conduct statistical inference by means of the nonparametric sum of ranking differences and comparison of ranks to random numbers (SRD-CRRN) method. We try to find the RCCS predictors best representing the general consensus regarding secondary structural propensities. The existence and the magnitude of resulting differences on secondary structure determination under varying sample conditions (temperature, pH) are demonstrated and discussed for globular proteins and especially IDPs.
Collapse
Affiliation(s)
- Dániel Kovács
- ELTE, Eötvös Loránd University, Institute of Chemistry, Analytical and BioNMR Laboratory Pázmány Péter sétány 1/A Budapest 1117 Hungary
- Eötvös Loránd University, Hevesy György PhD School of Chemistry Pázmány Péter sétány 1/A Budapest 1117 Hungary
| | - Andrea Bodor
- ELTE, Eötvös Loránd University, Institute of Chemistry, Analytical and BioNMR Laboratory Pázmány Péter sétány 1/A Budapest 1117 Hungary
| |
Collapse
|
6
|
Mortensen JC, Damjanovic J, Miao J, Hui T, Lin Y. A backbone-dependent rotamer library with high (ϕ, ψ) coverage using metadynamics simulations. Protein Sci 2022; 31:e4491. [PMID: 36327064 PMCID: PMC9679973 DOI: 10.1002/pro.4491] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2022] [Revised: 10/26/2022] [Accepted: 10/28/2022] [Indexed: 12/06/2023]
Abstract
Backbone-dependent rotamer libraries are commonly used to assign the side chain dihedral angles of amino acids when modeling protein structures. Most rotamer libraries are created by curating protein crystal structure data and using various methods to extrapolate the existing data to cover all possible backbone conformations. However, these rotamer libraries may not be suitable for modeling the structures of cyclic peptides and other constrained peptides because these molecules frequently sample backbone conformations rarely seen in the crystal structures of linear proteins. To provide backbone-dependent side chain information beyond the α-helix, β-sheet, and PPII regions, we used explicit-solvent metadynamics simulations of model dipeptides to create a new rotamer library that has high coverage in the (ϕ, ψ) space. Furthermore, this approach can be applied to build high-coverage rotamer libraries for noncanonical amino acids. The resulting Metadynamics of Dipeptides for Rotamer Distribution (MEDFORD) rotamer library predicts the side chain conformations of high-resolution protein crystal structures with similar accuracy (~80%) to a state-of-the-art rotamer library. Our ability to test the accuracy of MEDFORD at predicting the side chain dihedral angles of amino acids in noncanonical backbone conformation is restricted by the limited structural data available for cyclic peptides. For the cyclic peptide data that are currently available, MEDFORD and the state-of-the-art rotamer library perform comparably. However, the two rotamer libraries indeed make different rotamer predictions in noncanonical (ϕ, ψ) regions. For noncanonical amino acids, the MEDFORD rotamer library predicts the χ1 values with approximately 75% accuracy.
Collapse
Affiliation(s)
| | | | - Jiayuan Miao
- Department of ChemistryTufts UniversityMedfordMassachusettsUSA
| | - Tiffani Hui
- Department of ChemistryTufts UniversityMedfordMassachusettsUSA
| | - Yu‐Shan Lin
- Department of ChemistryTufts UniversityMedfordMassachusettsUSA
| |
Collapse
|
7
|
Cui X, Liu H, Chen HF. Polarizable Force Field of Intrinsically Disordered Proteins with CMAP and Reweighting Optimization. J Chem Inf Model 2022; 62:4970-4982. [PMID: 36178373 DOI: 10.1021/acs.jcim.2c00835] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Intrinsically disordered proteins (IDPs) are highly structurally heterogeneous without a specific tertiary structure under physiology conditions and play key roles in the development of human diseases. Due to the characteristics of diverse conformations, as one of the important methods, molecular dynamics simulation can complement information for experimental methods. Because of the enrichment for charged amino acids for IDPs, polarizable force fields should be a good choice for the simulation of IDPs. However, current polarizable force fields are limited in sampling conformer features of IDPs. Therefore, a polarizable force field was released and named Drude2019IDP based on Drude2019 with reweighting and grid-based potential energy correction map optimization. In order to evaluate the performance of Drude2019IDP, 16 dipeptides, 18 short peptides, 3 representative IDPs, and 5 structural proteins were simulated. The results show that the NMR observables driven by Drude2019IDP are in better agreement with the experiment data than those by Drude2019 on short peptides and IDPs. Drude2019IDP can sample more diverse conformations than Drude2019. Furthermore, the performances of the two force fields are similar to the sample ordered proteins. These results confirm that the developed Drude2019IDP can improve the reproduction of conformers for intrinsically disordered proteins and can be used to gain insight into the paradigm of sequence-disorder for IDPs.
Collapse
Affiliation(s)
- Xiaochen Cui
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai200240, China
| | - Hao Liu
- Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai200240, China
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai200240, China.,Shanghai Center for Bioinformation Technology, Shanghai200235, China
| |
Collapse
|
8
|
Schweitzer-Stenner R. Exploring Nearest Neighbor Interactions and Their Influence on the Gibbs Energy Landscape of Unfolded Proteins and Peptides. Int J Mol Sci 2022; 23:ijms23105643. [PMID: 35628453 PMCID: PMC9147007 DOI: 10.3390/ijms23105643] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2022] [Revised: 05/12/2022] [Accepted: 05/13/2022] [Indexed: 11/17/2022] Open
Abstract
The Flory isolated pair hypothesis (IPH) is one of the corner stones of the random coil model, which is generally invoked to describe the conformational dynamics of unfolded and intrinsically disordered proteins (IDPs). It stipulates, that individual residues sample the entire sterically allowed space of the Ramachandran plot without exhibiting any correlations with the conformational dynamics of its neighbors. However, multiple lines of computational, bioinformatic and experimental evidence suggest that nearest neighbors have a significant influence on the conformational sampling of amino acid residues. This implies that the conformational entropy of unfolded polypeptides and proteins is much less than one would expect based on the Ramachandran plots of individual residues. A further implication is that the Gibbs energies of residues in unfolded proteins or polypeptides are not additive. This review provides an overview of what is currently known and what has yet to be explored regarding nearest neighbor interactions in unfolded proteins.
Collapse
|
9
|
Stoppelman JP, Ng TT, Nerenberg PS, Wang LP. Development and Validation of AMBER-FB15-Compatible Force Field Parameters for Phosphorylated Amino Acids. J Phys Chem B 2021; 125:11927-11942. [PMID: 34668708 DOI: 10.1021/acs.jpcb.1c07547] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
Phosphorylation of select amino acid residues is one of the most common biological mechanisms for regulating protein structures and functions. While computational modeling can be used to explore the detailed structural changes associated with phosphorylation, most molecular mechanics force fields developed for the simulation of phosphoproteins have been noted to be inconsistent with experimental data. In this work, we parameterize force fields for the phosphorylated forms of the amino acids serine, threonine, and tyrosine using the ForceBalance software package with the goal of improving agreement with experiments for these residues. Our optimized force field, denoted as FB18, is parameterized using high-quality ab initio potential energy scans and is designed to be fully compatible with the AMBER-FB15 protein force field. When utilized in MD simulations together with the TIP3P-FB water model, we find that FB18 consistently enhances the prediction of experimental quantities such as 3J NMR couplings and intramolecular hydrogen-bonding propensities in comparison to previously published models. As was reported with AMBER-FB15, we also see improved agreement with the reference QM calculations in regions at and away from local minima. We thus believe that the FB18 parameter set provides a promising route for the further investigation of the varied effects of protein phosphorylation.
Collapse
Affiliation(s)
- John P Stoppelman
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Tracey T Ng
- Department of Physics & Astronomy, California State University, Los Angeles, California 90032, United States
| | - Paul S Nerenberg
- Department of Physics & Astronomy, California State University, Los Angeles, California 90032, United States.,Department of Biological Sciences, California State University, Los Angeles, California 90032, United States
| | - Lee-Ping Wang
- Department of Chemistry, University of California, Davis, California 95616, United States
| |
Collapse
|
10
|
Yuan Y, Wang F. A comparison of three DFT exchange-correlation functionals and two basis sets for the prediction of the conformation distribution of hydrated polyglycine. J Chem Phys 2021; 155:094104. [PMID: 34496578 DOI: 10.1063/5.0059669] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The performance of three density functional theory (DFT) exchange-correlation functionals, namely, Perdew-Burke-Ernzerhof (PBE), BP86, and B3LYP, in predicting conformational distributions of a hydrated glycine peptide is tested with two different basis sets in the framework of adaptive force matching (AFM). The conformational distributions yielded the free energy profiles of the DFT functional and basis set combinations. Unlike traditional validations of potential energy and structural parameters, our approach allows the free energy of DFT to be validated. When compared to experimental distributions, the def2-TZVP basis set provides better agreement than a slightly trimmed aug-cc-pVDZ basis set. B3LYP is shown to be better than BP86 and PBE. The glycine model fitted against B3LYP-D3(BJ) with the def2-TZVP basis set is the most accurate and named the AFM2021 model for glycine. The AFM2021 glycine model provides better agreement with experimental J-coupling constants than C36m and ff14SB, although the margin is very small when compared to C36m. Our previously published alanine model is also refitted with the slightly simplified AFM2021 energy expression. This work shows good promise of AFM for developing force fields for a range of proteinogenic peptides using only DFT as reference.
Collapse
Affiliation(s)
- Ying Yuan
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Feng Wang
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| |
Collapse
|
11
|
Zhang S, Andrews B, Schweitzer-Stenner R, Urbanc B. Intrinsic Conformational Dynamics of Alanine in Water/Ethanol Mixtures: An Experiment-Driven Molecular Dynamics Study. J Phys Chem B 2020; 124:11600-11616. [PMID: 33300341 DOI: 10.1021/acs.jpcb.0c08245] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In vitro, cationic glycylalanylglycine (GAG) forms a hydrogel in binary mixtures of water and ethanol. In water, alanine residue is known for its high polyproline II (pPII) content. Spectroscopic data, including three J-coupling constants and amide I' profiles, indicate that addition of 42% ethanol to water significantly reduces the pPII content of alanine residue in GAG. Here, experiment-based Gaussian Ramachandran distributions of alanine in GAG at different ethanol fractions are examined and three MD force fields are evaluated with respect to their ability to capture these ethanol-induced conformational changes. MD simulations on monomeric GAG in eight different water/ethanol mixtures within Amber ff14SB, OPLS-AA/M, and CHARMM36m reveal that only Amber ff14SB partially captures the ethanol-induced conformational changes of alanine residue in monomeric GAG when 42% ethanol is added to water. MD simulations of 200 mM GAG ensembles in pure water and in the aqueous solution with 42% ethanol showcase the ability of CHARMM36m to capture the effect of ethanol on the average pPII content of alanine in GAG and provide a plausible explanation for this effect, which may stem from an increased propensity of GAG to form oligomers in the presence of ethanol.
Collapse
Affiliation(s)
- Shuting Zhang
- Department of Physics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| | - Brian Andrews
- Department of Physics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| | | | - Brigita Urbanc
- Department of Physics, Drexel University, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
12
|
Glycine in Water Favors the Polyproline II State. Biomolecules 2020; 10:biom10081121. [PMID: 32751224 PMCID: PMC7463814 DOI: 10.3390/biom10081121] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2020] [Revised: 07/22/2020] [Accepted: 07/27/2020] [Indexed: 01/09/2023] Open
Abstract
Conformational preferences of amino acid residues in water are determined by the backbone and side-chain properties. Alanine is known for its high polyproline II (pPII) propensity. The question of relative contributions of the backbone and side chain to the conformational preferences of alanine and other amino acid residues in water is not fully resolved. Because glycine lacks a heavy-atom side chain, glycine-based peptides can be used to examine to which extent the backbone properties affect the conformational space. Here, we use published spectroscopic data for the central glycine residue of cationic triglycine in water to demonstrate that its conformational space is dominated by the pPII state. We assess three commonly used molecular dynamics (MD) force fields with respect to their ability to capture the conformational preferences of the central glycine residue in triglycine. We show that pPII is the mesostate that enables the functional backbone groups of the central residue to form the most hydrogen bonds with water. Our results indicate that the pPII propensity of the central glycine in GGG is comparable to that of alanine in GAG, implying that the water-backbone hydrogen bonding is responsible for the high pPII content of these residues.
Collapse
|
13
|
Song D, Liu H, Luo R, Chen HF. Environment-Specific Force Field for Intrinsically Disordered and Ordered Proteins. J Chem Inf Model 2020; 60:2257-2267. [PMID: 32227937 PMCID: PMC10449432 DOI: 10.1021/acs.jcim.0c00059] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The need for accurate and efficient force fields for modeling 3D structures of macrobiomolecules and in particular intrinsically disordered proteins (IDPs) has increased with recent findings to associate IDPs and human diseases. However, most conventional protein force fields and recent IDP-specific force fields are limited in reproducing accurate structural features of IDPs. Here, we present an environmental specific precise force field (ESFF1) based on CMAP corrections of 71 different sequence environments to improve the accuracy and efficiency of MD simulation for both IDPs and folded proteins. MD simulations of 84 different short peptides, IDPs, and structured proteins show that ESFF1 can accurately reproduce spectroscopic properties for different peptides and proteins whether they are disordered or ordered. The successful ab initio folding of five fast-folding proteins further supports the reliability of ESFF1. The extensive analysis documented here shows that ESFF1 is able to achieve a reasonable balance between ordered and disordered states in protein simulations.
Collapse
Affiliation(s)
- Dong Song
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Hao Liu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Ray Luo
- Departments of Molecular Biology and Biochemistry, Chemical and Molecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine, California 92697-3900, United States
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
- Shanghai Center for Bioinformation Technology, Shanghai 200235, China
| |
Collapse
|
14
|
Zhai L, Otani Y, Hori Y, Tomita T, Ohwada T. Peptide-based short single β-strand mimics without hydrogen bonding or aggregation. Chem Commun (Camb) 2020; 56:1573-1576. [PMID: 31930273 DOI: 10.1039/c9cc08378b] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Generation of short peptides with a single β-strand structure in solution is difficult. Herein, we design a new class of single β-strand peptidic mimics that are stable without self-aggregation in protic and non-protic solvents. Introduction of one present β-strand mimic can induce and propagate the β-strand structure for at least a penta-peptide sequence.
Collapse
Affiliation(s)
- Luhan Zhai
- Laboratory of Organic and Medicinal Chemistry, Graduate School of Pharmaceutical Sciences, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-0033, Japan.
| | | | | | | | | |
Collapse
|
15
|
Diem M, Oostenbrink C. Hamiltonian Reweighing To Refine Protein Backbone Dihedral Angle Parameters in the GROMOS Force Field. J Chem Inf Model 2020; 60:279-288. [PMID: 31873012 DOI: 10.1021/acs.jcim.9b01034] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics simulations of proteins depend critically on the underlying force field, which may be parameterized against experimental data or high-quality quantum calculations. Here, we develop search algorithms based on Monte Carlo and steepest descent calculations to optimize the backbone dihedral angle parameters from a single reference simulation. We apply these tools to improve the agreement between simulations of single, capped amino acids and experimentally determined J values and secondary structure propensities of these molecules. The parameters are further refined based on simulations of a set of seven proteins and finally validated in simulations on a large set of 52 protein structures. Improvements in the dihedral angle distributions are observed, and structural propensities of the proteins are reproduced very well. Overall, the GROMOS 54A8_bb parameter set forms an improvement to previous parameter sets, both for small molecules and for protein simulations.
Collapse
Affiliation(s)
- Matthias Diem
- Institute for Molecular Modeling and Simulation , University of Natural Resources and Life Sciences , Muthgasse 18 , 1190 Vienna , Austria
| | - Chris Oostenbrink
- Institute for Molecular Modeling and Simulation , University of Natural Resources and Life Sciences , Muthgasse 18 , 1190 Vienna , Austria
| |
Collapse
|
16
|
Yu L, Li DW, Brüschweiler R. Balanced Amino-Acid-Specific Molecular Dynamics Force Field for the Realistic Simulation of Both Folded and Disordered Proteins. J Chem Theory Comput 2019; 16:1311-1318. [DOI: 10.1021/acs.jctc.9b01062] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
|
17
|
Tian C, Kasavajhala K, Belfon KAA, Raguette L, Huang H, Migues AN, Bickel J, Wang Y, Pincay J, Wu Q, Simmerling C. ff19SB: Amino-Acid-Specific Protein Backbone Parameters Trained against Quantum Mechanics Energy Surfaces in Solution. J Chem Theory Comput 2019; 16:528-552. [PMID: 31714766 DOI: 10.1021/acs.jctc.9b00591] [Citation(s) in RCA: 848] [Impact Index Per Article: 169.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Molecular dynamics (MD) simulations have become increasingly popular in studying the motions and functions of biomolecules. The accuracy of the simulation, however, is highly determined by the molecular mechanics (MM) force field (FF), a set of functions with adjustable parameters to compute the potential energies from atomic positions. However, the overall quality of the FF, such as our previously published ff99SB and ff14SB, can be limited by assumptions that were made years ago. In the updated model presented here (ff19SB), we have significantly improved the backbone profiles for all 20 amino acids. We fit coupled φ/ψ parameters using 2D φ/ψ conformational scans for multiple amino acids, using as reference data the entire 2D quantum mechanics (QM) energy surface. We address the polarization inconsistency during dihedral parameter fitting by using both QM and MM in aqueous solution. Finally, we examine possible dependency of the backbone fitting on side chain rotamer. To extensively validate ff19SB parameters, and to compare to results using other Amber models, we have performed a total of ∼5 ms MD simulations in explicit solvent. Our results show that after amino-acid-specific training against QM data with solvent polarization, ff19SB not only reproduces the differences in amino-acid-specific Protein Data Bank (PDB) Ramachandran maps better but also shows significantly improved capability to differentiate amino-acid-dependent properties such as helical propensities. We also conclude that an inherent underestimation of helicity is present in ff14SB, which is (inexactly) compensated for by an increase in helical content driven by the TIP3P bias toward overly compact structures. In summary, ff19SB, when combined with a more accurate water model such as OPC, should have better predictive power for modeling sequence-specific behavior, protein mutations, and also rational protein design. Of the explicit water models tested here, we recommend use of OPC with ff19SB.
Collapse
Affiliation(s)
- Chuan Tian
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Koushik Kasavajhala
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Kellon A A Belfon
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Lauren Raguette
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - He Huang
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Angela N Migues
- Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - John Bickel
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Yuzhang Wang
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Jorge Pincay
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Qin Wu
- Center for Functional Nanomaterials , Brookhaven National Laboratory , Upton , New York 11973 , United States
| | - Carlos Simmerling
- Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States.,Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States
| |
Collapse
|
18
|
Wang S, Otani Y, Zhai L, Su A, Nara M, Kawahata M, Yamaguchi K, Sada A, Ohki R, Ohwada T. Overall Shape Constraint of Alternating α/β-Hybrid Peptides Containing Bicyclic β-Proline. Org Lett 2019; 21:7813-7817. [PMID: 31518151 DOI: 10.1021/acs.orglett.9b02799] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Our NMR, IR/Raman, CD spectroscopic, and X-ray crystallographic studies, as well as accelerated molecular dynamics simulations, showed that alternating hybrid α/β-peptides containing a bicyclic β-proline surrogate form unique extended curved folds, regardless of the peptide length and solvent environment. It is suggested that extended β/PPII structures are preferred in the insulating α-alanine moieties between the rigid bicyclic β-proline structures. These hybrid peptides inhibit p53-MDM2 and p53-MDMX protein-protein interactions.
Collapse
Affiliation(s)
- Siyuan Wang
- Graduate School of Pharmaceutical Sciences , University of Tokyo , 7-3-1 Hongo , Bunkyo-ku , Tokyo 113-0033 , Japan.,Research Foundation Itsuu Laboratory , C1232 Kanagawa Science Park R&D Building, 3-2-1 Sakado, Takatsu-ku , Kawasaki , Kanagawa 213-0012 , Japan
| | - Yuko Otani
- Graduate School of Pharmaceutical Sciences , University of Tokyo , 7-3-1 Hongo , Bunkyo-ku , Tokyo 113-0033 , Japan
| | - Luhan Zhai
- Graduate School of Pharmaceutical Sciences , University of Tokyo , 7-3-1 Hongo , Bunkyo-ku , Tokyo 113-0033 , Japan
| | - Aoze Su
- Graduate School of Pharmaceutical Sciences , University of Tokyo , 7-3-1 Hongo , Bunkyo-ku , Tokyo 113-0033 , Japan
| | - Masayuki Nara
- Department of Chemistry, College of Liberal Arts and Sciences , Tokyo Medical and Dental University , 2-8-30 Kohnodai , Ichikawa , Chiba 272-0827 , Japan
| | - Masatoshi Kawahata
- Department of Pharmaceutical Sciences at Kagawa Campus , Tokushima Bunri University , 1314-1 Shido , Sanuki , Kagawa 769-2193 , Japan
| | - Kentaro Yamaguchi
- Department of Pharmaceutical Sciences at Kagawa Campus , Tokushima Bunri University , 1314-1 Shido , Sanuki , Kagawa 769-2193 , Japan
| | - Akane Sada
- Laboratory of Fundamental Oncology , National Cancer Center Research Institute , Tsukiji 5-1-1 , Chuo-ku , Tokyo 104-0045 , Japan
| | - Rieko Ohki
- Laboratory of Fundamental Oncology , National Cancer Center Research Institute , Tsukiji 5-1-1 , Chuo-ku , Tokyo 104-0045 , Japan
| | - Tomohiko Ohwada
- Graduate School of Pharmaceutical Sciences , University of Tokyo , 7-3-1 Hongo , Bunkyo-ku , Tokyo 113-0033 , Japan
| |
Collapse
|
19
|
Allen AA, Robertson MJ, Payne MC, Cole DJ. Development and Validation of the Quantum Mechanical Bespoke Protein Force Field. ACS OMEGA 2019; 4:14537-14550. [PMID: 31528808 PMCID: PMC6740169 DOI: 10.1021/acsomega.9b01769] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/14/2019] [Accepted: 07/30/2019] [Indexed: 06/10/2023]
Abstract
Molecular mechanics force field parameters for macromolecules, such as proteins, are traditionally fit to reproduce experimental properties of small molecules, and thus, they neglect system-specific polarization. In this paper, we introduce a complete protein force field that is designed to be compatible with the quantum mechanical bespoke (QUBE) force field by deriving nonbonded parameters directly from the electron density of the specific protein under study. The main backbone and sidechain protein torsional parameters are rederived in this work by fitting to quantum mechanical dihedral scans for compatibility with QUBE nonbonded parameters. Software is provided for the preparation of QUBE input files. The accuracy of the new force field, and the derived torsional parameters, is tested by comparing the conformational preferences of a range of peptides and proteins with experimental measurements. Accurate backbone and sidechain conformations are obtained in molecular dynamics simulations of dipeptides, with NMR J coupling errors comparable to the widely used OPLS force field. In simulations of five folded proteins, the secondary structure is generally retained, and the NMR J coupling errors are similar to standard transferable force fields, although some loss of the experimental structure is observed in certain regions of the proteins. With several avenues for further development, the use of system-specific nonbonded force field parameters is a promising approach for next-generation simulations of biological molecules.
Collapse
Affiliation(s)
- Alice
E. A. Allen
- TCM
Group, Cavendish Laboratory, 19 JJ Thomson Ave, Cambridge CB3 0HE, United Kingdom
| | - Michael J. Robertson
- Department of Molecular and Cellular Physiology and Department of Structural Biology Stanford University School of Medicine, 279 Campus Drive, Stanford, California 94305, United States
| | - Michael C. Payne
- TCM
Group, Cavendish Laboratory, 19 JJ Thomson Ave, Cambridge CB3 0HE, United Kingdom
| | - Daniel J. Cole
- School
of Natural and Environmental Sciences, Newcastle
University, Newcastle
upon Tyne NE1 7RU, United
Kingdom
| |
Collapse
|
20
|
Denèfle T, Pramil E, Gómez-Morales L, Levasseur MD, Lardé E, Newton C, Herry K, Herbi L, Lamotte Y, Odile E, Ancellin N, Grondin P, Martinez-Torres AC, Viviani F, Merle-Beral H, Lequin O, Susin SA, Karoyan P. Homotrimerization Approach in the Design of Thrombospondin-1 Mimetic Peptides with Improved Potency in Triggering Regulated Cell Death of Cancer Cells. J Med Chem 2019; 62:7656-7668. [PMID: 31403795 DOI: 10.1021/acs.jmedchem.9b00024] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
In order to optimize the potency of the first serum-stable peptide agonist of CD47 (PKHB1) in triggering regulated cell death of cancer cells, we designed a maturation process aimed to mimic the trimeric structure of the thrombospondin-1/CD47 binding epitope. For that purpose, an N-methylation scan of the PKHB1 sequence was realized to prevent peptide aggregation. Structural and pharmacological analyses were conducted in order to assess the conformational impact of these chemical modifications on the backbone structure and the biological activity. This structure-activity relationship study led to the discovery of a highly soluble N-methylated peptide that we termed PKT16. Afterward, this monomer was used for the design of a homotrimeric peptide mimic that we termed [PKT16]3, which proved to be 10-fold more potent than its monomeric counterpart. A pharmacological evaluation of [PKT16]3 in inducing cell death of adherent (A549) and nonadherent (MEC-1) cancer cell lines was also performed.
Collapse
Affiliation(s)
- Thomas Denèfle
- Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM , 75005 Paris , France.,Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM, Site OncoDesign , 25-27 Avenue du Québec , 91140 Les Ulis , France
| | - Elodie Pramil
- Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM , 75005 Paris , France.,Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM, Site OncoDesign , 25-27 Avenue du Québec , 91140 Les Ulis , France.,Cell Death and Drug Resistance in Lymphoproliferative Disorders Team , Centre de Recherche des Cordeliers, INSERM UMRS 1138 , 75006 Paris , France
| | - Luis Gómez-Morales
- Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM , 75005 Paris , France.,Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM, Site OncoDesign , 25-27 Avenue du Québec , 91140 Les Ulis , France.,Laboratory of Immunology and Virology , Autonomous University of Nuevo Leon , 66451 San Nicolas de los Garza , NL , Mexico
| | - Mikail D Levasseur
- Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM , 75005 Paris , France.,Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM, Site OncoDesign , 25-27 Avenue du Québec , 91140 Les Ulis , France
| | - Eva Lardé
- Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM , 75005 Paris , France.,Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM, Site OncoDesign , 25-27 Avenue du Québec , 91140 Les Ulis , France
| | - Clara Newton
- Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM , 75005 Paris , France.,Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM, Site OncoDesign , 25-27 Avenue du Québec , 91140 Les Ulis , France
| | - Kenny Herry
- OncoDesign , 25 Avenue du Québec , 91140 Les Ulis , France
| | - Linda Herbi
- Cell Death and Drug Resistance in Lymphoproliferative Disorders Team , Centre de Recherche des Cordeliers, INSERM UMRS 1138 , 75006 Paris , France
| | - Yann Lamotte
- OncoDesign , 25 Avenue du Québec , 91140 Les Ulis , France
| | - Estelle Odile
- Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM , 75005 Paris , France.,Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM, Site OncoDesign , 25-27 Avenue du Québec , 91140 Les Ulis , France
| | | | - Pascal Grondin
- OncoDesign , 25 Avenue du Québec , 91140 Les Ulis , France
| | - Ana-Carolina Martinez-Torres
- Laboratory of Immunology and Virology , Autonomous University of Nuevo Leon , 66451 San Nicolas de los Garza , NL , Mexico
| | | | - Hélène Merle-Beral
- Cell Death and Drug Resistance in Lymphoproliferative Disorders Team , Centre de Recherche des Cordeliers, INSERM UMRS 1138 , 75006 Paris , France
| | - Olivier Lequin
- Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM , 75005 Paris , France
| | - Santos A Susin
- Cell Death and Drug Resistance in Lymphoproliferative Disorders Team , Centre de Recherche des Cordeliers, INSERM UMRS 1138 , 75006 Paris , France
| | - Philippe Karoyan
- Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM , 75005 Paris , France.,Sorbonne Université , Ecole Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM, Site OncoDesign , 25-27 Avenue du Québec , 91140 Les Ulis , France.,SiRIC CURAMUS (CANCER UNITED RESEARCH ASSOCIATING MEDICINE, UNIVERSITY & SOCIETY, Site de Recherche Intégrée sur le Cancer) IUC, AP-HP.6 , Sorbonne Université 75005 Paris , France.,Kayvisa AG , Industriestrasse, 44 , 6300 Zug , Switzerland.,Kaybiotix GmbH , Zugerstrasse 32 , 6340 Baar , Switzerland
| |
Collapse
|
21
|
Turupcu A, Diem M, Smith LJ, Oostenbrink C. Structural Aspects of the O-glycosylation Linkage in Glycopeptides via MD Simulations and Comparison with NMR Experiments. Chemphyschem 2019; 20:1527-1537. [PMID: 30920077 PMCID: PMC6563056 DOI: 10.1002/cphc.201900079] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2019] [Revised: 03/27/2019] [Indexed: 11/26/2022]
Abstract
A powerful conformational searching and enhanced sampling simulation method, and unbiased molecular dynamics simulations have been used along with NMR spectroscopic observables to provide a detailed structural view of O-glycosylation. For four model systems, the force-field parameters can accurately predict experimental NMR observables (J couplings and NOE's). This enables us to derive conclusions based on the generated ensembles, in which O-glycosylation affects the peptide backbone conformation by forcing it towards to an extended conformation. An exception is described for β-GalNAc-Thr where the α content is increased and stabilized via hydrogen bonding between the sugar and the peptide backbone, which was not observed in the rest of the studied systems. These observations might offer an explanation for the evolutionary preference of α-linked GalNAc glycosylation instead of a β link.
Collapse
Affiliation(s)
- Aysegül Turupcu
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and SimulationUniversity of Natural Resources and Life Sciences ViennaViennaAustria
| | - Matthias Diem
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and SimulationUniversity of Natural Resources and Life Sciences ViennaViennaAustria
| | | | - Chris Oostenbrink
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and SimulationUniversity of Natural Resources and Life Sciences ViennaViennaAustria
| |
Collapse
|
22
|
Amino acid conformations control the morphological and chiral features of the self-assembled peptide nanostructures: Young investigators perspective. J Colloid Interface Sci 2019; 548:244-254. [PMID: 31004957 DOI: 10.1016/j.jcis.2019.04.019] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2018] [Revised: 04/03/2019] [Accepted: 04/05/2019] [Indexed: 01/11/2023]
Abstract
HYPOTHESIS A variety of nanostructures with different chiral features can be self-assembled from short peptides with highly similar sequences. We hypothesize that these supramolecular nanostructures are ruled by the constituent amino acid residues which adopt their conformations under the influence of intra-/inter-molecular interactions during peptide self-assembly. APPROACH Through reviewing recent advances in the self-assembly of short peptides and focusing on the relationship between amino acid conformations, peptide secondary structures and intra-/inter-molecular interactions within the supramolecular architectures, we aim to rationalize the complex interactive processes involved in the self-assembly of short, designed peptides. RESULTS Given the highly complexing interactive processes, the adoption of amino acid conformations and their control over peptide self-assembly consist of 4 main steps: (1) Each amino acid residue adopts its unique conformation in a specific sequence; (2) The sequence exhibits its own main chain geometry and determines the propensity of the intermolecular alignment within the building block; (3) The structural propensity of the building block and the packing mode between them determine the self-assembled structural features such as twisting, growth and chirality; (4) In addition to intra-/inter-molecular interactions, inter-sheet and inter-building block interactions could also affect the residue conformations and nanostructures, causing structural readjustment.
Collapse
|
23
|
Choi JM, Pappu RV. Experimentally Derived and Computationally Optimized Backbone Conformational Statistics for Blocked Amino Acids. J Chem Theory Comput 2019; 15:1355-1366. [PMID: 30516982 PMCID: PMC10846683 DOI: 10.1021/acs.jctc.8b00572] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Experimentally derived, amino acid specific backbone dihedral angle distributions are invaluable for modeling data-driven conformational equilibria of proteins and for enabling quantitative assessments of the accuracies of molecular mechanics force fields. The protein coil library that is extracted from analysis of high-resolution structures of proteins has served as a useful proxy for quantifying intrinsic and context-dependent conformational distributions of amino acids. However, data that go into coil libraries will have hidden biases, and ad hoc procedures must be used to remove these biases. Here, we combine high-resolution biased information from protein structural databases with unbiased low-resolution information from spectroscopic measurements of blocked amino acids to obtain experimentally derived and computationally optimized coil-library landscapes for each of the 20 naturally occurring amino acids. Quantitative descriptions of conformational distributions require parsing of data into conformational basins with defined envelopes, centers, and statistical weights. We develop and deploy a numerical method to extract conformational basins. The weights of conformational basins are optimized to reproduce quantitative inferences drawn from spectroscopic experiments for blocked amino acids. The optimized distributions serve as touchstones for assessments of intrinsic conformational preferences and for quantitative comparisons of molecular mechanics force fields.
Collapse
Affiliation(s)
- Jeong-Mo Choi
- Department of Biomedical Engineering and Center for Biological Systems Engineering, Washington University in St. Louis, One Brookings Drive, Campus Box 1097, St. Louis, Missouri 63130
| | - Rohit V. Pappu
- Department of Biomedical Engineering and Center for Biological Systems Engineering, Washington University in St. Louis, One Brookings Drive, Campus Box 1097, St. Louis, Missouri 63130
| |
Collapse
|
24
|
Choi JM, Pappu RV. Improvements to the ABSINTH Force Field for Proteins Based on Experimentally Derived Amino Acid Specific Backbone Conformational Statistics. J Chem Theory Comput 2019; 15:1367-1382. [PMID: 30633502 PMCID: PMC10749164 DOI: 10.1021/acs.jctc.8b00573] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
We present an improved version of the ABSINTH implicit solvation model and force field paradigm (termed ABSINTH-C) by incorporating a grid-based term that bootstraps against experimentally derived and computationally optimized conformational statistics for blocked amino acids. These statistics provide high-resolution descriptions of the intrinsic backbone dihedral angle preferences for all 20 amino acids. The original ABSINTH model generates Ramachandran plots that are too shallow in terms of the basin structures and too permissive in terms of dihedral angle preferences. We bootstrap against the reference optimized landscapes and incorporate CMAP-like residue-specific terms that help us reproduce the intrinsic dihedral angle preferences of individual amino acids. These corrections that lead to ABSINTH-C are achieved by balancing the incorporation of the new residue-specific terms with the accuracies inherent to the original ABSINTH model. We demonstrate the robustness of ABSINTH-C through a series of examples to highlight the preservation of accuracies as well as examples that demonstrate the improvements. Our efforts show how the recent experimentally derived and computationally optimized coil-library landscapes can be used as a touchstone for quantifying errors and making improvements to molecular mechanics force fields.
Collapse
Affiliation(s)
- Jeong-Mo Choi
- Department of Biomedical Engineering and Center for Biological Systems Engineering, Washington University in St. Louis, One Brookings Drive, Campus Box 1097, St. Louis, Missouri 63130
| | - Rohit V. Pappu
- Department of Biomedical Engineering and Center for Biological Systems Engineering, Washington University in St. Louis, One Brookings Drive, Campus Box 1097, St. Louis, Missouri 63130
| |
Collapse
|
25
|
Vymětal J, Jurásková V, Vondrášek J. AMBER and CHARMM Force Fields Inconsistently Portray the Microscopic Details of Phosphorylation. J Chem Theory Comput 2018; 15:665-679. [PMID: 30468703 DOI: 10.1021/acs.jctc.8b00715] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Phosphorylation of serine, threonine, and tyrosine is one of the most frequently occurring and crucial post-translational modifications of proteins often associated with important structural and functional changes. We investigated the direct effect of phosphorylation on the intrinsic conformational preferences of amino acids as a potential trigger of larger structural events. We conducted a comparative study of force fields on terminally capped amino acids (dipeptides) as the simplest model for phosphorylation. Our bias-exchange metadynamics simulations revealed that all model dipeptides sampled a great heterogeneity of ensembles affected by introduction of mono- and dianionic phosphate groups. However, the detected changes in populations of backbone conformers and side-chain rotamers did not reveal a strong discriminatory shift in preferences, as could be anticipated for the bulky, charged phosphate group. Furthermore, the AMBER and CHARMM force fields provided inconsistent populations of individual conformers as well as net structural trends upon phosphorylation. Detailed analysis of ensembles revealed competition between hydration and formation of internal hydrogen bonds involving amide hydrogens and the phosphate group. The observed difference in hydration free energy and potential for hydrogen bonding in individual force fields could be attributed to the different partial atomic charges used in each force field and, hence, the different parametrization strategies. Nevertheless, conformational propensities and net structural changes upon phosphorylation are difficult to extract from experimental measurements, and existing experimental data provide limited guidance for force field assessment and further development.
Collapse
Affiliation(s)
- Jiří Vymětal
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences , Flemingovo náměstí 542/2 , 166 10 Praha 6 , Czech Republic
| | - Veronika Jurásková
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences , Flemingovo náměstí 542/2 , 166 10 Praha 6 , Czech Republic
| | - Jiří Vondrášek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences , Flemingovo náměstí 542/2 , 166 10 Praha 6 , Czech Republic
| |
Collapse
|
26
|
Zhai L, Wang S, Nara M, Takeuchi K, Shimada I, Otani Y, Ohwada T. Application of C-Terminal 7-Azabicyclo[2.2.1]heptane to Stabilize β-Strand-like Extended Conformation of a Neighboring α-Amino Acid. J Org Chem 2018; 83:13063-13079. [DOI: 10.1021/acs.joc.8b01756] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Luhan Zhai
- Laboratory of Organic and Medicinal Chemistry, Graduate School of Pharmaceutical Sciences, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Siyuan Wang
- Laboratory of Organic and Medicinal Chemistry, Graduate School of Pharmaceutical Sciences, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Masayuki Nara
- Department of Chemistry, College of Liberal Arts and Sciences, Tokyo Medical and Dental University, Ichikawa, Chiba 272-0827, Japan
| | - Koh Takeuchi
- Molecular Profiling Research Center for Drug Discovery, National Institute of Advanced Industrial Science and Technology (AIST), Aomi, Koto-ku, Tokyo 135-0064, Japan
| | - Ichio Shimada
- Molecular Profiling Research Center for Drug Discovery, National Institute of Advanced Industrial Science and Technology (AIST), Aomi, Koto-ku, Tokyo 135-0064, Japan
- Laboratory of Physical Chemistry, Graduate School of Pharmaceutical Sciences, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Yuko Otani
- Laboratory of Organic and Medicinal Chemistry, Graduate School of Pharmaceutical Sciences, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Tomohiko Ohwada
- Laboratory of Organic and Medicinal Chemistry, Graduate School of Pharmaceutical Sciences, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| |
Collapse
|
27
|
Margreitter C, Oostenbrink C. Correction to Optimization of Protein Backbone Dihedral Angles by Means of Hamiltonian Reweighting. J Chem Inf Model 2018; 58:1716-1720. [PMID: 30062876 PMCID: PMC8154260 DOI: 10.1021/acs.jcim.8b00470] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Christian Margreitter
- Institute for molecular modeling and simulation , University of natural resources and life sciences , Muthgasse 18 , 1190 Vienna , Austria
| | - Chris Oostenbrink
- Institute for molecular modeling and simulation , University of natural resources and life sciences , Muthgasse 18 , 1190 Vienna , Austria
| |
Collapse
|
28
|
Kaminker* R, Anastasaki A, Gutekunst WR, Luo Y, Lee S, Hawker* CJ. Tuning of protease resistance in oligopeptides through N-alkylation. Chem Commun (Camb) 2018; 54:9631-9634. [PMID: 30095837 PMCID: PMC6141190 DOI: 10.1039/c8cc04407d] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
N-Methylation of amino acids is an effective way to create protease resistance in both natural and synthetic peptides. However, alkyl substituents other than N-methyl have not been extensively studied. Here, we prepare and examine a series of N-substituted peptides in which the size and length of the alkyl group is modulated. These design insights provide a unique and modular handle for tuning proteolysis in oligopeptides.
Collapse
Affiliation(s)
- R. Kaminker*
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States. ;
| | - A. Anastasaki
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States. ;
| | - W. R. Gutekunst
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States. ;
| | - Y. Luo
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States. ;
| | - S. Lee
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States. ;
| | - C. J. Hawker*
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States. ;
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United Sates
| |
Collapse
|
29
|
Zhang Y, Zhou Y, He L, Fu Y, Zhang W, Hu J, Shi Z. Hydration effects on Leu's polyproline II population in AcLXPNH 2. Chem Commun (Camb) 2018; 54:5764-5767. [PMID: 29781018 DOI: 10.1039/c8cc02402b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Hydration is important in many fundamental processes. To investigate hydration effects on peptide conformations, we examined neighboring-residue and side-chain blocking effects in AcLXPNH2. A correlation between two effects suggests that hydration stabilizes PII more than β-structures. Our results are important for understanding the hydration effects on peptide conformations and hydration-forces in general.
Collapse
Affiliation(s)
- Yan Zhang
- School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 430074, P. R. China.
| | | | | | | | | | | | | |
Collapse
|
30
|
Structure-activity relationships of β-hairpin mimics as modulators of amyloid β-peptide aggregation. Eur J Med Chem 2018; 154:280-293. [DOI: 10.1016/j.ejmech.2018.05.018] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2018] [Revised: 05/11/2018] [Accepted: 05/12/2018] [Indexed: 12/23/2022]
|
31
|
Xu Y, Correia I, Ha-Duong T, Kihal N, Soulier JL, Kaffy J, Crousse B, Lequin O, Ongeri S. The use of 4,4,4-trifluorothreonine to stabilize extended peptide structures and mimic β-strands. Beilstein J Org Chem 2018; 13:2842-2853. [PMID: 29564012 PMCID: PMC5753055 DOI: 10.3762/bjoc.13.276] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2017] [Accepted: 12/13/2017] [Indexed: 01/06/2023] Open
Abstract
Pentapeptides having the sequence R-HN-Ala-Val-X-Val-Leu-OMe, where the central residue X is L-serine, L-threonine, (2S,3R)-L-CF3-threonine and (2S,3S)-L-CF3-threonine were prepared. The capacity of (2S,3S)- and (2S,3R)-CF3-threonine analogues to stabilize an extended structure when introduced in the central position of pentapeptides is demonstrated by NMR conformational studies and molecular dynamics simulations. CF3-threonine containing pentapeptides are more prone to mimic β-strands than their natural Ser and Thr pentapeptide analogues. The proof of concept that these fluorinated β-strand mimics are able to disrupt protein–protein interactions involving β-sheet structures is provided. The CF3-threonine containing pentapeptides interact with the amyloid peptide Aβ1-42 in order to reduce the protein–protein interactions mediating its aggregation process.
Collapse
Affiliation(s)
- Yaochun Xu
- Molécules Fluorées et Chimie Médicinale, BioCIS, Univ. Paris-Sud, CNRS, Université Paris Saclay, 5 rue Jean-Baptiste Clément, 92296 Châtenay-Malabry Cedex, France
| | - Isabelle Correia
- Sorbonne Universités, UPMC Univ Paris 06, Ecole Normale Supérieure, PSL Research University, CNRS, Laboratoire des Biomolécules, 4 place Jussieu, 75252 Paris Cedex 05, France
| | - Tap Ha-Duong
- Molécules Fluorées et Chimie Médicinale, BioCIS, Univ. Paris-Sud, CNRS, Université Paris Saclay, 5 rue Jean-Baptiste Clément, 92296 Châtenay-Malabry Cedex, France
| | - Nadjib Kihal
- Molécules Fluorées et Chimie Médicinale, BioCIS, Univ. Paris-Sud, CNRS, Université Paris Saclay, 5 rue Jean-Baptiste Clément, 92296 Châtenay-Malabry Cedex, France
| | - Jean-Louis Soulier
- Molécules Fluorées et Chimie Médicinale, BioCIS, Univ. Paris-Sud, CNRS, Université Paris Saclay, 5 rue Jean-Baptiste Clément, 92296 Châtenay-Malabry Cedex, France
| | - Julia Kaffy
- Molécules Fluorées et Chimie Médicinale, BioCIS, Univ. Paris-Sud, CNRS, Université Paris Saclay, 5 rue Jean-Baptiste Clément, 92296 Châtenay-Malabry Cedex, France
| | - Benoît Crousse
- Molécules Fluorées et Chimie Médicinale, BioCIS, Univ. Paris-Sud, CNRS, Université Paris Saclay, 5 rue Jean-Baptiste Clément, 92296 Châtenay-Malabry Cedex, France
| | - Olivier Lequin
- Sorbonne Universités, UPMC Univ Paris 06, Ecole Normale Supérieure, PSL Research University, CNRS, Laboratoire des Biomolécules, 4 place Jussieu, 75252 Paris Cedex 05, France
| | - Sandrine Ongeri
- Molécules Fluorées et Chimie Médicinale, BioCIS, Univ. Paris-Sud, CNRS, Université Paris Saclay, 5 rue Jean-Baptiste Clément, 92296 Châtenay-Malabry Cedex, France
| |
Collapse
|
32
|
Rani L, Mallajosyula SS. Phosphorylation versus O-GlcNAcylation: Computational Insights into the Differential Influences of the Two Competitive Post-Translational Modifications. J Phys Chem B 2017; 121:10618-10638. [DOI: 10.1021/acs.jpcb.7b08790] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Lata Rani
- Department of Chemistry, Indian Institute of Technology Gandhinagar, Palaj, Gandhinagar, Gujarat, India - 382355
| | - Sairam S. Mallajosyula
- Department of Chemistry, Indian Institute of Technology Gandhinagar, Palaj, Gandhinagar, Gujarat, India - 382355
| |
Collapse
|
33
|
Rubio-Martinez J, Tomas MS, Perez JJ. Effect of the solvent on the conformational behavior of the alanine dipeptide deduced from MD simulations. J Mol Graph Model 2017; 78:118-128. [PMID: 29055185 DOI: 10.1016/j.jmgm.2017.10.005] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2017] [Revised: 10/07/2017] [Accepted: 10/09/2017] [Indexed: 11/19/2022]
Abstract
In general, peptides do not exhibit a well-defined conformational profile in solution. However, despite the experimental blurred picture associated with their structure, compelling spectroscopic evidence shows that peptides exhibit local order. The conformational profile of a peptide is the result of a balance between intramolecular interactions between different atoms of the molecule and intermolecular interactions between atoms of the molecule and the solvent. Accordingly, the conformational profile of a peptide will change upon the properties of the solvent it is soaked. To get insight into the balance between intra- and intermolecular interactions on the conformational preferences of the peptide backbone we have studied the conformational profile of the alanine dipeptide in diverse solvents using molecular dynamics as sampling technique. Solvents studied include chloroform, methanol, dimethyl sulfoxide, water and N-methylacetamide. Different treatments of the solvent have been studied in the present work including explicit solvent molecules, a generalized Born model and using the bulk dielectric constant of the solvent. The diverse calculations identify four major conformations with different populations in the diverse solvents: the C7eq only sampled in chloroform; the C5 or extended conformation; the polyproline (PII) conformation and the right-handed α-helix conformation (αR). The results of present calculations permit to analyze how the balance between intra- and intermolecular interactions explains the populations of the diverse conformations observed.
Collapse
Affiliation(s)
- Jaime Rubio-Martinez
- Dept. of Physical Chemistry, Faculty of Chemistry, Universitat de Barcelona and the Institut de Recerca en Quimica Teorica i Computacional (IQTCUB), Mati i Franques 1-3, 08028 Barcelona, Spain
| | - M Santos Tomas
- Department of Architecture Technology, Universitat Politecnica de Catalunya, Av. Diagonal, 649, E-08028 Barcelona, Spain
| | - Juan J Perez
- Department of Chemical Engineering, Universitat Politecnica de Catalunya- Barcelona Tech, Av. Diagonal, 647, 08028 Barcelona, Spain.
| |
Collapse
|
34
|
Denèfle T, Boullet H, Herbi L, Newton C, Martinez-Torres AC, Guez A, Pramil E, Quiney C, Pourcelot M, Levasseur MD, Lardé E, Moumné R, Ogi FX, Grondin P, Merle-Beral H, Lequin O, Susin SA, Karoyan P. Thrombospondin-1 Mimetic Agonist Peptides Induce Selective Death in Tumor Cells: Design, Synthesis, and Structure-Activity Relationship Studies. J Med Chem 2016; 59:8412-21. [PMID: 27526615 DOI: 10.1021/acs.jmedchem.6b00781] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Thrombospondin-1 (TSP-1) is a glycoprotein considered as a key actor within the tumor microenvironment. Its binding to CD47, a cell surface receptor, triggers programmed cell death. Previous studies allowed the identification of 4N1K decapeptide derived from the TSP-1/CD47 binding epitope. Here, we demonstrate that this peptide is able to induce selective apoptosis of various cancer cell lines while sparing normal cells. A structure-activity relationship study led to the design of the first serum stable TSP-1 mimetic agonist peptide able to trigger selective programmed cell death (PCD) of at least lung, breast, and colorectal cancer cells. Altogether, these results will be of valuable interest for further investigation in the design of potent CD47 agonist peptides, opening new perspectives for the development of original anticancer therapies.
Collapse
Affiliation(s)
- Thomas Denèfle
- Sorbonne Universités, UPMC Université Paris 06, Ecole Normale Supérieure, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,Département de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,UPMC Université Paris 06, Laboratoire des BioMolécules, Site GSK, 25-27 Avenue du Québec, 91140 Les Ulis, France
| | - Héloise Boullet
- Sorbonne Universités, UPMC Université Paris 06, Ecole Normale Supérieure, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,Département de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules, 75005 Paris, France
| | - Linda Herbi
- Cell Death and Drug Resistance in Lymphoproliferative Disorders Team, Centre de Recherche des Cordeliers, INSERM UMRS, 1138, Paris, France.,Sorbonne Universités, UPMC Université Paris 06, UMRS 1138, Paris, France.,Université Paris Descartes, Sorbonne Paris Cité, UMRS 1138, Paris, France
| | - Clara Newton
- Sorbonne Universités, UPMC Université Paris 06, Ecole Normale Supérieure, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,Département de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,UPMC Université Paris 06, Laboratoire des BioMolécules, Site GSK, 25-27 Avenue du Québec, 91140 Les Ulis, France
| | - Ana-Carolina Martinez-Torres
- Cell Death and Drug Resistance in Lymphoproliferative Disorders Team, Centre de Recherche des Cordeliers, INSERM UMRS, 1138, Paris, France.,Sorbonne Universités, UPMC Université Paris 06, UMRS 1138, Paris, France.,Université Paris Descartes, Sorbonne Paris Cité, UMRS 1138, Paris, France
| | - Alexandre Guez
- Sorbonne Universités, UPMC Université Paris 06, Ecole Normale Supérieure, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,Département de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules, 75005 Paris, France
| | - Elodie Pramil
- Sorbonne Universités, UPMC Université Paris 06, Ecole Normale Supérieure, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,Département de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,Cell Death and Drug Resistance in Lymphoproliferative Disorders Team, Centre de Recherche des Cordeliers, INSERM UMRS, 1138, Paris, France.,Sorbonne Universités, UPMC Université Paris 06, UMRS 1138, Paris, France.,Université Paris Descartes, Sorbonne Paris Cité, UMRS 1138, Paris, France
| | - Claire Quiney
- Cell Death and Drug Resistance in Lymphoproliferative Disorders Team, Centre de Recherche des Cordeliers, INSERM UMRS, 1138, Paris, France.,Sorbonne Universités, UPMC Université Paris 06, UMRS 1138, Paris, France.,Université Paris Descartes, Sorbonne Paris Cité, UMRS 1138, Paris, France
| | - Marilyne Pourcelot
- Sorbonne Universités, UPMC Université Paris 06, Ecole Normale Supérieure, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,Département de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,UPMC Université Paris 06, Laboratoire des BioMolécules, Site GSK, 25-27 Avenue du Québec, 91140 Les Ulis, France
| | - Mikail D Levasseur
- Sorbonne Universités, UPMC Université Paris 06, Ecole Normale Supérieure, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,Département de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,UPMC Université Paris 06, Laboratoire des BioMolécules, Site GSK, 25-27 Avenue du Québec, 91140 Les Ulis, France
| | - Eva Lardé
- Sorbonne Universités, UPMC Université Paris 06, Ecole Normale Supérieure, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,Département de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,UPMC Université Paris 06, Laboratoire des BioMolécules, Site GSK, 25-27 Avenue du Québec, 91140 Les Ulis, France
| | - Roba Moumné
- Sorbonne Universités, UPMC Université Paris 06, Ecole Normale Supérieure, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,Département de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules, 75005 Paris, France
| | | | | | - Hélène Merle-Beral
- Cell Death and Drug Resistance in Lymphoproliferative Disorders Team, Centre de Recherche des Cordeliers, INSERM UMRS, 1138, Paris, France.,Sorbonne Universités, UPMC Université Paris 06, UMRS 1138, Paris, France.,AP-HP, GH Pitié-Salpêtrière, Service d'Hématologie Biologique, 75013 Paris, France
| | - Olivier Lequin
- Sorbonne Universités, UPMC Université Paris 06, Ecole Normale Supérieure, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,Département de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules, 75005 Paris, France
| | - Santos A Susin
- Cell Death and Drug Resistance in Lymphoproliferative Disorders Team, Centre de Recherche des Cordeliers, INSERM UMRS, 1138, Paris, France.,Sorbonne Universités, UPMC Université Paris 06, UMRS 1138, Paris, France.,Université Paris Descartes, Sorbonne Paris Cité, UMRS 1138, Paris, France
| | - Philippe Karoyan
- Sorbonne Universités, UPMC Université Paris 06, Ecole Normale Supérieure, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,Département de Chimie, Ecole Normale Supérieure, PSL Research University, UPMC Université Paris 06, CNRS, Laboratoire des Biomolécules, 75005 Paris, France.,UPMC Université Paris 06, Laboratoire des BioMolécules, Site GSK, 25-27 Avenue du Québec, 91140 Les Ulis, France
| |
Collapse
|
35
|
Margreitter C, Oostenbrink C. Optimization of Protein Backbone Dihedral Angles by Means of Hamiltonian Reweighting. J Chem Inf Model 2016; 56:1823-34. [PMID: 27559757 PMCID: PMC5039763 DOI: 10.1021/acs.jcim.6b00399] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Molecular dynamics
simulations depend critically on the accuracy
of the underlying force fields in properly representing biomolecules.
Hence, it is crucial to validate the force-field parameter sets in
this respect. In the context of the GROMOS force field, this is usually
achieved by comparing simulation data to experimental observables
for small molecules. In this study, we develop new amino acid backbone
dihedral angle potential energy parameters based on the widely used
54A7 parameter set by matching to experimental J values
and secondary structure propensity scales. In order to find the most
appropriate backbone parameters, close to 100 000 different
combinations of parameters have been screened. However, since the
sheer number of combinations considered prohibits actual molecular
dynamics simulations for each of them, we instead predicted the values
for every combination using Hamiltonian reweighting. While the original
54A7 parameter set fails to reproduce the experimental data, we are
able to provide parameters that match significantly better. However,
to ensure applicability in the context of larger peptides and full
proteins, further studies have to be undertaken.
Collapse
Affiliation(s)
- Christian Margreitter
- Institute for Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Muthgasse 18, 1190 Vienna, Austria
| | - Chris Oostenbrink
- Institute for Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Muthgasse 18, 1190 Vienna, Austria
| |
Collapse
|
36
|
Bünnemann K, Merten C. Solvation of N,C-Protected Valine: Interactions with DMSO and a Chiral Solvating Agent. J Phys Chem B 2016; 120:9434-42. [DOI: 10.1021/acs.jpcb.6b05897] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Christian Merten
- Organische
Chemie 2, Ruhr Universität Bochum, 44801 Bochum, Germany
| |
Collapse
|
37
|
Ilawe NV, Raeber AE, Schweitzer-Stenner R, Toal SE, Wong BM. Assessing backbone solvation effects in the conformational propensities of amino acid residues in unfolded peptides. Phys Chem Chem Phys 2016; 17:24917-24. [PMID: 26343224 DOI: 10.1039/c5cp03646a] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Conformational ensembles of individual amino acid residues within model GxG peptides (x representing different amino acid residues) are dominated by a mixture of polyproline II (pPII) and β-strand like conformations. We recently discovered rather substantial differences between the enthalpic and entropic contributions to this equilibrium for different amino acid residues. Isoleucine and valine exceed all other amino acid residues in terms of their rather large enthalpic stabilization and entropic destabilization of polyproline II. In order to shed light on these underlying physical mechanisms, we performed high-level DFT calculations to explore the energetics of four representative GxG peptides where x = alanine (A), leucine (L), valine (V), and isoleucine (I) in explicit water (10 H2O molecules with a polarizable continuum water model) and in vacuo. We found that the large energetic contributions to the stabilization of pPII result, to a major extent, from peptide-water, water-water interactions, and changes of the solvent self-energy. Differences between the peptide-solvent interaction energies of hydration in pPII and β-strand peptides are particularly important for the pPII ⇌ β equilibria of the more aliphatic peptides GIG and GLG. Furthermore, we performed a vibrational analysis of the four peptides in both conformations and discovered a rather substantial mixing between water motions and peptide vibrations below 700 cm(-1). We found that the respective vibrational entropies are substantially different for the considered conformations, and their contributions to the Gibbs/Helmholtz energy stabilize β-strand conformations. Taken together, our results underscore the notion of the solvent being the predominant determinant of peptide (and protein) conformations in the unfolded state.
Collapse
Affiliation(s)
- Niranjan V Ilawe
- Department of Chemical & Environmental Engineering, and Materials Science & Engineering Program, University of California, Riverside, 900 University Avenue, Riverside, CA 92521, USA.
| | | | | | | | | |
Collapse
|
38
|
Li S, Andrews CT, Frembgen-Kesner T, Miller MS, Siemonsma SL, Collingsworth TD, Rockafellow IT, Ngo NA, Campbell BA, Brown RF, Guo C, Schrodt M, Liu YT, Elcock AH. Molecular Dynamics Simulations of 441 Two-Residue Peptides in Aqueous Solution: Conformational Preferences and Neighboring Residue Effects with the Amber ff99SB-ildn-NMR Force Field. J Chem Theory Comput 2016; 11:1315-29. [PMID: 26579777 DOI: 10.1021/ct5010966] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Understanding the intrinsic conformational preferences of amino acids and the extent to which they are modulated by neighboring residues is a key issue for developing predictive models of protein folding and stability. Here we present the results of 441 independent explicit-solvent MD simulations of all possible two-residue peptides that contain the 20 standard amino acids with histidine modeled in both its neutral and protonated states. (3)J(HNHα) coupling constants and δ(Hα) chemical shifts calculated from the MD simulations correlate quite well with recently published experimental measurements for a corresponding set of two-residue peptides. Neighboring residue effects (NREs) on the average (3)J(HNHα) and δ(Hα) values of adjacent residues are also reasonably well reproduced, with the large NREs exerted experimentally by aromatic residues, in particular, being accurately captured. NREs on the secondary structure preferences of adjacent amino acids have been computed and compared with corresponding effects observed in a coil library and the average β-turn preferences of all amino acid types have been determined. Finally, the intrinsic conformational preferences of histidine, and its NREs on the conformational preferences of adjacent residues, are both shown to be strongly affected by the protonation state of the imidazole ring.
Collapse
Affiliation(s)
- Shuxiang Li
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Casey T Andrews
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | | | - Mark S Miller
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Stephen L Siemonsma
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | | | - Isaac T Rockafellow
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Nguyet Anh Ngo
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Brady A Campbell
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Reid F Brown
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Chengxuan Guo
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Michael Schrodt
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Yu-Tsan Liu
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Adrian H Elcock
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| |
Collapse
|
39
|
Childers MC, Towse CL, Daggett V. The effect of chirality and steric hindrance on intrinsic backbone conformational propensities: tools for protein design. Protein Eng Des Sel 2016; 29:271-80. [PMID: 27284086 DOI: 10.1093/protein/gzw023] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2016] [Accepted: 05/11/2016] [Indexed: 01/30/2023] Open
Abstract
The conformational propensities of amino acids are an amalgamation of sequence effects, environmental effects and underlying intrinsic behavior. Many have attempted to investigate neighboring residue effects to aid in our understanding of protein folding and improve structure prediction efforts, especially with respect to difficult to characterize states, such as disordered or unfolded states. Host-guest peptide series are a useful tool in examining the propensities of the amino acids free from the surrounding protein structure. Here, we compare the distributions of the backbone dihedral angles (φ/ψ) of the 20 proteogenic amino acids in two different sequence contexts using the AAXAA and GGXGG host-guest pentapeptide series. We further examine their intrinsic behaviors across three environmental contexts: water at 298 K, water at 498 K, and 8 M urea at 298 K. The GGXGG systems provide the intrinsic amino acid propensities devoid of any conformational context. The alanine residues in the AAXAA series enforce backbone chirality, thereby providing a model of the intrinsic behavior of amino acids in a protein chain. Our results show modest differences in φ/ψ distributions due to the steric constraints of the Ala side chains, the magnitudes of which are dependent on the denaturing conditions. One of the strongest factors modulating φ/ψ distributions was the protonation of titratable side chains, and the largest differences observed were in the amino acid propensities for the rarely sampled αL region.
Collapse
Affiliation(s)
| | - Clare-Louise Towse
- Department of Bioengineering, University of Washington, Seattle, WA 98195-5013, USA
| | - Valerie Daggett
- Department of Bioengineering, University of Washington, Seattle, WA 98195-5013, USA
| |
Collapse
|
40
|
Zhou Y, He L, Zhang W, Hu J, Shi Z. Populations of the Minor α-Conformation in AcGXGNH2 and the α-Helical Nucleation Propensities. Sci Rep 2016; 6:27197. [PMID: 27256621 PMCID: PMC4891685 DOI: 10.1038/srep27197] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2016] [Accepted: 05/16/2016] [Indexed: 01/25/2023] Open
Abstract
Intrinsic backbone conformational preferences of different amino acids are important for understanding the local structure of unfolded protein chains. Recent evidence suggests α-structure is relatively minor among three major backbone conformations for unfolded proteins. The α-helices are the dominant structures in many proteins. For these proteins, how could the α-structures occur from the least in unfolded to the most in folded states? Populations of the minor α-conformation in model peptides provide vital information. Reliable determination of populations of the α-conformers in these peptides that exist in multiple equilibriums of different conformations remains a challenge. Combined analyses on data from AcGXPNH2 and AcGXGNH2 peptides allow us to derive the populations of PII, β and α in AcGXGNH2. Our results show that on average residue X in AcGXGNH2 adopt PII, β, and α 44.7%, 44.5% and 10.8% of time, respectively. The contents of α-conformations for different amino acids define an α-helix nucleation propensity scale. With derived PII, β and α-contents, we can construct a free energy-conformation diagram on each AcGXGNH2 in aqueous solution for the three major backbone conformations. Our results would have broad implications on early-stage events of protein folding.
Collapse
Affiliation(s)
- Yanjun Zhou
- School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 430074, P.R. China
| | - Liu He
- School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 430074, P.R. China
| | - Wenwen Zhang
- School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 430074, P.R. China
| | - Jingjing Hu
- School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 430074, P.R. China
| | - Zhengshuang Shi
- School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 430074, P.R. China
| |
Collapse
|
41
|
Mirtič A, Merzel F, Grdadolnik J. The amide III vibrational circular dichroism band as a probe to detect conformational preferences of alanine dipeptide in water. Biopolymers 2016; 101:814-8. [PMID: 24436080 DOI: 10.1002/bip.22460] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2013] [Revised: 01/08/2014] [Accepted: 01/08/2014] [Indexed: 01/24/2023]
Abstract
The conformational preferences of blocked alanine dipeptide (ADP), Ac-Ala-NHMe, in aqueous solution were studied using vibrational circular dichroism (VCD) together with density functional theory (DFT) calculations. DFT calculations of three most representative conformations of ADP surrounded by six explicit water molecules immersed in a dielectric continuum have proven high sensitivity of amide III VCD band shape that is characteristic for each conformation of the peptide backbone. The polyproline II (PII ) and αR conformation of ADP are associated with a positive VCD band while β conformation has a negative VCD band in amide III region. Knowing this spectral characteristic of each conformation allows us to assign the experimental amide III VCD spectrum of ADP. Moreover, the amide III region of the VCD spectrum was used to determine the relative populations of conformations of ADP in water. Based on the interpretation of the amide III region of VCD spectrum we have shown that dominant conformation of ADP in water is PII which is stabilized by hydrogen bonded water molecules between CO and NH groups on the peptide backbone.
Collapse
Affiliation(s)
- Andreja Mirtič
- National Institute of Chemistry, Hajdrihova 19,, SI-1000, Ljubljana, Slovenia
| | | | | |
Collapse
|
42
|
Christensen A, Kubař T, Cui Q, Elstner M. Semiempirical Quantum Mechanical Methods for Noncovalent Interactions for Chemical and Biochemical Applications. Chem Rev 2016; 116:5301-37. [PMID: 27074247 PMCID: PMC4867870 DOI: 10.1021/acs.chemrev.5b00584] [Citation(s) in RCA: 246] [Impact Index Per Article: 30.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2015] [Indexed: 12/28/2022]
Abstract
Semiempirical (SE) methods can be derived from either Hartree-Fock or density functional theory by applying systematic approximations, leading to efficient computational schemes that are several orders of magnitude faster than ab initio calculations. Such numerical efficiency, in combination with modern computational facilities and linear scaling algorithms, allows application of SE methods to very large molecular systems with extensive conformational sampling. To reliably model the structure, dynamics, and reactivity of biological and other soft matter systems, however, good accuracy for the description of noncovalent interactions is required. In this review, we analyze popular SE approaches in terms of their ability to model noncovalent interactions, especially in the context of describing biomolecules, water solution, and organic materials. We discuss the most significant errors and proposed correction schemes, and we review their performance using standard test sets of molecular systems for quantum chemical methods and several recent applications. The general goal is to highlight both the value and limitations of SE methods and stimulate further developments that allow them to effectively complement ab initio methods in the analysis of complex molecular systems.
Collapse
Affiliation(s)
- Anders
S. Christensen
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Tomáš Kubař
- Institute of Physical
Chemistry & Center for Functional Nanostructures and Institute of Physical
Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Qiang Cui
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Marcus Elstner
- Institute of Physical
Chemistry & Center for Functional Nanostructures and Institute of Physical
Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| |
Collapse
|
43
|
Almeida GG, Cordeiro JMM, Martín ME, Aguilar MA. Conformational Changes of the Alanine Dipeptide in Water–Ethanol Binary Mixtures. J Chem Theory Comput 2016; 12:1514-24. [DOI: 10.1021/acs.jctc.5b00952] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Glauco G. Almeida
- Faculdade de Engenharia,
Universidade Estadual Paulista “Julio de Mesquita Filho″, Ilha Solteira 15385-000, Brasil
| | - João M. M. Cordeiro
- Faculdade de Engenharia,
Universidade Estadual Paulista “Julio de Mesquita Filho″, Ilha Solteira 15385-000, Brasil
| | - M. Elena Martín
- Área
de Química Física, University of Extremadura, Avda.
Elvas s/n, Edif. José Ma Viguera
Lobo, 3a planta, Badajoz 06006, Spain
| | - Manuel A. Aguilar
- Área
de Química Física, University of Extremadura, Avda.
Elvas s/n, Edif. José Ma Viguera
Lobo, 3a planta, Badajoz 06006, Spain
| |
Collapse
|
44
|
Appadurai R, Senapati S. Dynamical Network of HIV-1 Protease Mutants Reveals the Mechanism of Drug Resistance and Unhindered Activity. Biochemistry 2016; 55:1529-40. [PMID: 26892689 DOI: 10.1021/acs.biochem.5b00946] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
HIV-1 protease variants resist drugs by active and non-active-site mutations. The active-site mutations, which are the primary or first set of mutations, hamper the stability of the enzyme and resist the drugs minimally. As a result, secondary mutations that not only increase protein stability for unhindered catalytic activity but also resist drugs very effectively arise. While the mechanism of drug resistance of the active-site mutations is through modulating the active-site pocket volume, the mechanism of drug resistance of the non-active-site mutations is unclear. Moreover, how these allosteric mutations, which are 8-21 Å distant, communicate to the active site for drug efflux is completely unexplored. Results from molecular dynamics simulations suggest that the primary mechanism of drug resistance of the secondary mutations involves opening of the flexible protease flaps. Results from both residue- and community-based network analyses reveal that this precise action of protease is accomplished by the presence of robust communication paths between the mutational sites and the functionally relevant regions: active site and flaps. While the communication is more direct in the wild type, it traverses across multiple intermediate residues in mutants, leading to weak signaling and unregulated motions of flaps. The global integrity of the protease network is, however, maintained through the neighboring residues, which exhibit high degrees of conservation, consistent with clinical data and mutagenesis studies.
Collapse
Affiliation(s)
- Rajeswari Appadurai
- BJM School of Biosciences and Department of Biotechnology, Indian Institution of Technology Madras , Chennai 600 036, India
| | - Sanjib Senapati
- BJM School of Biosciences and Department of Biotechnology, Indian Institution of Technology Madras , Chennai 600 036, India
| |
Collapse
|
45
|
Kubař T, Welke K, Groenhof G. New QM/MM implementation of the DFTB3 method in the gromacs package. J Comput Chem 2015; 36:1978-89. [PMID: 26238364 DOI: 10.1002/jcc.24029] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2015] [Accepted: 07/09/2015] [Indexed: 01/07/2023]
Abstract
The approximate density-functional tight-binding theory method DFTB3 has been implemented in the quantum mechanics/molecular mechanics (QM/MM) framework of the Gromacs molecular simulation package. We show that the efficient smooth particle-mesh Ewald implementation of Gromacs extends to the calculation of QM/MM electrostatic interactions. Further, we make use of the various free-energy functionalities provided by Gromacs and the PLUMED plugin. We exploit the versatility and performance of the current framework in three typical applications of QM/MM methods to solve biophysical problems: (i) ultrafast proton transfer in malonaldehyde, (ii) conformation of the alanine dipeptide, and (iii) electron-induced repair of a DNA lesion. Also discussed is the further development of the framework, regarding mostly the options for parallelization.
Collapse
Affiliation(s)
- Tomáš Kubař
- Institute of Physical Chemistry and Center for Functional Nanostructures, Karlsruhe Institute of Technology, Karlsruhe, 76131, Germany
| | - Kai Welke
- Department of Chemistry, Nagoya University, Nagoya, 464-8602, Japan
| | - Gerrit Groenhof
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, Jyväskylä, 40014, Finland
| |
Collapse
|
46
|
Robertson M, Tirado-Rives J, Jorgensen WL. Improved Peptide and Protein Torsional Energetics with the OPLSAA Force Field. J Chem Theory Comput 2015; 11:3499-509. [PMID: 26190950 PMCID: PMC4504185 DOI: 10.1021/acs.jctc.5b00356] [Citation(s) in RCA: 482] [Impact Index Per Article: 53.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2015] [Indexed: 11/28/2022]
Abstract
The development and validation of new peptide dihedral parameters are reported for the OPLS-AA force field. High accuracy quantum chemical methods were used to scan φ, ψ, χ1, and χ2 potential energy surfaces for blocked dipeptides. New Fourier coefficients for the dihedral angle terms of the OPLS-AA force field were fit to these surfaces, utilizing a Boltzmann-weighted error function and systematically examining the effects of weighting temperature. To prevent overfitting to the available data, a minimal number of new residue-specific and peptide-specific torsion terms were developed. Extensive experimental solution-phase and quantum chemical gas-phase benchmarks were used to assess the quality of the new parameters, named OPLS-AA/M, demonstrating significant improvement over previous OPLS-AA force fields. A Boltzmann weighting temperature of 2000 K was determined to be optimal for fitting the new Fourier coefficients for dihedral angle parameters. Conclusions are drawn from the results for best practices for developing new torsion parameters for protein force fields.
Collapse
Affiliation(s)
- Michael
J. Robertson
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - Julian Tirado-Rives
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - William L. Jorgensen
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| |
Collapse
|
47
|
Correct folding of an α-helix and a β-hairpin using a polarized 2D torsional potential. Sci Rep 2015; 5:10359. [PMID: 26039188 PMCID: PMC5380191 DOI: 10.1038/srep10359] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2014] [Accepted: 04/09/2015] [Indexed: 01/12/2023] Open
Abstract
A new modification to the AMBER force field that incorporates the coupled two-dimensional main chain torsion energy has been evaluated for the balanced representation of secondary structures. In this modified AMBER force field (AMBER03(2D)), the main chain torsion energy is represented by 2-dimensional Fourier expansions with parameters fitted to the potential energy surface generated by high-level quantum mechanical calculations of small peptides in solution. Molecular dynamics simulations are performed to study the folding of two model peptides adopting either α-helix or β-hairpin structures. Both peptides are successfully folded into their native structures using an AMBER03(2D) force field with the implementation of a polarization scheme (AMBER03(2D)p). For comparison, simulations using a standard AMBER03 force field with and without polarization, as well as AMBER03(2D) without polarization, fail to fold both peptides successfully. The correction to secondary structure propensity in the AMBER03 force field and the polarization effect are critical to folding Trpzip2; without these factors, a helical structure is obtained. This study strongly suggests that this new force field is capable of providing a more balanced preference for helical and extended conformations. The electrostatic polarization effect is shown to be indispensable to the growth of secondary structures.
Collapse
|
48
|
Toal SE, Kubatova N, Richter C, Linhard V, Schwalbe H, Schweitzer-Stenner R. Randomizing the unfolded state of peptides (and proteins) by nearest neighbor interactions between unlike residues. Chemistry 2015; 21:5173-92. [PMID: 25728043 DOI: 10.1002/chem.201406539] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2014] [Indexed: 12/29/2022]
Abstract
To explore the influence of nearest neighbors on conformational biases in unfolded peptides, we combined vibrational and 2D NMR spectroscopy to obtain the conformational distributions of selected "GxyG" host-guest peptides in aqueous solution: GDyG, GSyG, GxLG, GxVG, where x/y=A, K, L, V. Large changes of conformational propensities were observed due to nearest-neighbor interactions, at variance with the isolated pair hypothesis. We found that protonated aspartic acid and serine lose their above-the-average preference for turn-like structures in favor of polyproline II (pPII) populations in the presence of neighbors with bulky side chains. Such residues also decrease the above-the-average pPII preference of alanine. These observations suggest that the underlying mechanism involves a disruption of the hydration shell. Thermodynamic analysis of (3) J(H(N) ,H(α) ) (T) data for each x,y residue reveals that modest changes in the conformational ensemble masks larger changes of enthalpy and entropy governing the pPII↔β equilibrium indicating a significant residue dependent temperature dependence of the peptides' conformational ensembles. These results suggest that nearest-neighbor interactions between unlike residues act as conformational randomizers close to the enthalpy-entropy compensation temperature, eliminating intrinsic biases in favor of largely balanced pPII/β dominated ensembles at physiological temperatures.
Collapse
Affiliation(s)
- Siobhan E Toal
- Department of Chemistry, Drexel University, 3141 Chestnut Street, Philadelphia, PA 10104 (USA); Present address: Department of Biophysics and Biochemistry, Yale University, New Haven, CT 06250 (USA)
| | | | | | | | | | | |
Collapse
|
49
|
Beauchamp KA, Pande VS, Das R. Bayesian energy landscape tilting: towards concordant models of molecular ensembles. Biophys J 2014; 106:1381-90. [PMID: 24655513 DOI: 10.1016/j.bpj.2014.02.009] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2013] [Revised: 01/19/2014] [Accepted: 02/06/2014] [Indexed: 11/28/2022] Open
Abstract
Predicting biological structure has remained challenging for systems such as disordered proteins that take on myriad conformations. Hybrid simulation/experiment strategies have been undermined by difficulties in evaluating errors from computational model inaccuracies and data uncertainties. Building on recent proposals from maximum entropy theory and nonequilibrium thermodynamics, we address these issues through a Bayesian energy landscape tilting (BELT) scheme for computing Bayesian hyperensembles over conformational ensembles. BELT uses Markov chain Monte Carlo to directly sample maximum-entropy conformational ensembles consistent with a set of input experimental observables. To test this framework, we apply BELT to model trialanine, starting from disagreeing simulations with the force fields ff96, ff99, ff99sbnmr-ildn, CHARMM27, and OPLS-AA. BELT incorporation of limited chemical shift and (3)J measurements gives convergent values of the peptide's α, β, and PPII conformational populations in all cases. As a test of predictive power, all five BELT hyperensembles recover set-aside measurements not used in the fitting and report accurate errors, even when starting from highly inaccurate simulations. BELT's principled framework thus enables practical predictions for complex biomolecular systems from discordant simulations and sparse data.
Collapse
Affiliation(s)
- Kyle A Beauchamp
- Computational Biology Center, Memorial Sloan-Kettering Cancer Center, New York, New York
| | - Vijay S Pande
- Departments of Chemistry, Computer Science, and Structural Biology and Biophysics Program, Stanford University, Stanford, California.
| | - Rhiju Das
- Departments of Biochemistry and Physics and Biophysics Program, Stanford University, Stanford, California.
| |
Collapse
|
50
|
Zhou CY, Jiang F, Wu YD. Residue-Specific Force Field Based on Protein Coil Library. RSFF2: Modification of AMBER ff99SB. J Phys Chem B 2014; 119:1035-47. [DOI: 10.1021/jp5064676] [Citation(s) in RCA: 78] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Chen-Yang Zhou
- College
of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China
| | - Fan Jiang
- Laboratory of Computational Chemistry and Drug Design, Laboratory of Chemical Genomics, Peking University Shenzhen Graduate School, Shenzhen, 518055, China
| | - Yun-Dong Wu
- College
of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China
- Laboratory of Computational Chemistry and Drug Design, Laboratory of Chemical Genomics, Peking University Shenzhen Graduate School, Shenzhen, 518055, China
| |
Collapse
|