1
|
Kumar H, Kim P. Artificial intelligence in fusion protein three-dimensional structure prediction: Review and perspective. Clin Transl Med 2024; 14:e1789. [PMID: 39090739 PMCID: PMC11294035 DOI: 10.1002/ctm2.1789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2024] [Revised: 07/16/2024] [Accepted: 07/19/2024] [Indexed: 08/04/2024] Open
Abstract
Recent advancements in artificial intelligence (AI) have accelerated the prediction of unknown protein structures. However, accurately predicting the three-dimensional (3D) structures of fusion proteins remains a difficult task because the current AI-based protein structure predictions are focused on the WT proteins rather than on the newly fused proteins in nature. Following the central dogma of biology, fusion proteins are translated from fusion transcripts, which are made by transcribing the fusion genes between two different loci through the chromosomal rearrangements in cancer. Accurately predicting the 3D structures of fusion proteins is important for understanding the functional roles and mechanisms of action of new chimeric proteins. However, predicting their 3D structure using a template-based model is challenging because known template structures are often unavailable in databases. Deep learning (DL) models that utilize multi-level protein information have revolutionized the prediction of protein 3D structures. In this review paper, we highlighted the latest advancements and ongoing challenges in predicting the 3D structure of fusion proteins using DL models. We aim to explore both the advantages and challenges of employing AlphaFold2, RoseTTAFold, tr-Rosetta and D-I-TASSER for modelling the 3D structures. HIGHLIGHTS: This review provides the overall pipeline and landscape of the prediction of the 3D structure of fusion protein. This review provides the factors that should be considered in predicting the 3D structures of fusion proteins using AI approaches in each step. This review highlights the latest advancements and ongoing challenges in predicting the 3D structure of fusion proteins using deep learning models. This review explores the advantages and challenges of employing AlphaFold2, RoseTTAFold, tr-Rosetta, and D-I-TASSER to model 3D structures.
Collapse
Affiliation(s)
- Himansu Kumar
- Department of Bioinformatics and Systems MedicineMcWilliams School of Biomedical InformaticsThe University of Texas Health Science Center at HoustonHoustonTexasUSA
| | - Pora Kim
- Department of Bioinformatics and Systems MedicineMcWilliams School of Biomedical InformaticsThe University of Texas Health Science Center at HoustonHoustonTexasUSA
| |
Collapse
|
2
|
Setiadi J, Boothroyd S, Slochower DR, Dotson DL, Thompson MW, Wagner JR, Wang LP, Gilson MK. Tuning Potential Functions to Host-Guest Binding Data. J Chem Theory Comput 2024; 20:239-252. [PMID: 38147689 PMCID: PMC10838530 DOI: 10.1021/acs.jctc.3c01050] [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] [Indexed: 12/28/2023]
Abstract
Software to more rapidly and accurately predict protein-ligand binding affinities is of high interest for early-stage drug discovery, and physics-based methods are among the most widely used technologies for this purpose. The accuracy of these methods depends critically on the accuracy of the potential functions that they use. Potential functions are typically trained against a combination of quantum chemical and experimental data. However, although binding affinities are among the most important quantities to predict, experimental binding affinities have not to date been integrated into the experimental data set used to train potential functions. In recent years, the use of host-guest complexes as simple and tractable models of binding thermodynamics has gained popularity due to their small size and simplicity, relative to protein-ligand systems. Host-guest complexes can also avoid ambiguities that arise in protein-ligand systems such as uncertain protonation states. Thus, experimental host-guest binding data are an appealing additional data type to integrate into the experimental data set used to optimize potential functions. Here, we report the extension of the Open Force Field Evaluator framework to enable the systematic calculation of host-guest binding free energies and their gradients with respect to force field parameters, coupled with the curation of 126 host-guest complexes with available experimental binding free energies. As an initial application of this novel infrastructure, we optimized generalized Born (GB) cavity radii for the OBC2 GB implicit solvent model against experimental data for 36 host-guest systems. This refitting led to a dramatic improvement in accuracy for both the training set and a separate test set with 90 additional host-guest systems. The optimized radii also showed encouraging transferability from host-guest systems to 59 protein-ligand systems. However, the new radii are significantly smaller than the baseline radii and lead to excessively favorable hydration free energies (HFEs). Thus, users of the OBC2 GB model currently may choose between GB cavity radii that yield more accurate binding affinities and GB cavity radii that yield more accurate HFEs. We suspect that achieving good accuracy on both will require more far-reaching adjustments to the GB model. We note that binding free-energy calculations using the OBC2 model in OpenMM gain about a 10× speedup relative to corresponding explicit solvent calculations, suggesting a future role for implicit solvent absolute binding free-energy (ABFE) calculations in virtual compound screening. This study proves the principle of using host-guest systems to train potential functions that are transferrable to protein-ligand systems and provides an infrastructure that enables a range of applications.
Collapse
Affiliation(s)
- Jeffry Setiadi
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, 9255 Pharmacy Lane, La Jolla, California 92093, United States
| | - Simon Boothroyd
- Boothroyd Scientific Consulting Ltd., London WC2H 9JQ, U.K
- Psivant Therapeutics, Boston, Massachusetts 02210, United States
| | | | - David L Dotson
- Datryllic LLC, Phoenix, Arizona 85003, United States
- The Open Force Field Consortium, Open Molecular Software Foundation, Davis, California 95616, United States
| | - Matthew W Thompson
- The Open Force Field Consortium, Open Molecular Software Foundation, Davis, California 95616, United States
| | - Jeffrey R Wagner
- The Open Force Field Consortium, Open Molecular Software Foundation, Davis, California 95616, United States
| | - Lee-Ping Wang
- Chemistry Department, University of California Davis, Davis, California 95616, United States
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, 9255 Pharmacy Lane, La Jolla, California 92093, United States
| |
Collapse
|
3
|
Kim S, Hyun DG, Nam Y, Shin SJ, Im D, Kim HS, Leem SH, Park HH, Kim BH, Park YH, Cho E, Goddard WA, Kim DH, Kim HI, Moon M. Genipin and pyrogallol: Two natural small molecules targeting the modulation of disordered proteins in Alzheimer's disease. Biomed Pharmacother 2023; 168:115770. [PMID: 37865990 DOI: 10.1016/j.biopha.2023.115770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Revised: 10/03/2023] [Accepted: 10/18/2023] [Indexed: 10/24/2023] Open
Abstract
Alzheimer's disease (AD) is characterized by the aggregation of disordered proteins, such as amyloid beta (Aβ) and tau, leading to neurotoxicity and disease progression. Despite numerous efforts, effective inhibitors of Aβ and tau aggregates have not been developed. Thus, we aimed to screen natural small molecules from crude extracts that target various pathologies and are prescribed for patients with neurological diseases. In this study, we screened 162 natural small molecules prescribed for neurological diseases and identified genipin and pyrogallol as hit compounds capable of simultaneously regulating the aggregation of Aβ and tau K18. Moreover, we confirmed the dual modulatory effects of these compounds on the reduction of amyloid-mediated neurotoxicity in vitro and the disassembly of preformed Aβ42 and tau K18 fibrils. Furthermore, we observed the alleviatory effects of genipin and pyrogallol against AD-related pathologies in triple transgenic AD mice. Molecular dynamics and docking simulations revealed the molecular interaction dynamics of genipin and pyrogallol with Aβ42 and tau K18, providing insights into their suppression of aggregation. Our findings suggest the therapeutic potential of genipin and pyrogallol as dual modulators for the treatment of AD by inhibiting aggregation or promoting dissociation of Aβ and tau.
Collapse
Affiliation(s)
- Sujin Kim
- Department of Biochemistry, College of Medicine, Konyang University, Daejeon 35365, the Republic of Korea; Research Institute for Dementia Science, Konyang University, Daejeon 35365, the Republic of Korea
| | - Da Gyeong Hyun
- Department of Chemistry, Korea University, Seoul 02841, the Republic of Korea; Center for Proteogenome Research, Korea University, Seoul 02841, the Republic of Korea; Single Cell Analysis Laboratory, Korea University, Seoul 02841, the Republic of Korea; Division of Chemistry and Chemical Engineering and Materials Process and Simulation Center, California Institute of Technology, Pasadena, CA 91125, United States
| | - Yunkwon Nam
- Department of Biochemistry, College of Medicine, Konyang University, Daejeon 35365, the Republic of Korea
| | - Soo Jung Shin
- Department of Biochemistry, College of Medicine, Konyang University, Daejeon 35365, the Republic of Korea
| | - Dongjoon Im
- Department of Chemistry, Korea University, Seoul 02841, the Republic of Korea; Center for Proteogenome Research, Korea University, Seoul 02841, the Republic of Korea; Single Cell Analysis Laboratory, Korea University, Seoul 02841, the Republic of Korea; Division of Chemistry and Chemical Engineering and Materials Process and Simulation Center, California Institute of Technology, Pasadena, CA 91125, United States
| | - Hyeon Soo Kim
- Department of Biochemistry, College of Medicine, Konyang University, Daejeon 35365, the Republic of Korea
| | - Seol Hwa Leem
- Department of Biochemistry, College of Medicine, Konyang University, Daejeon 35365, the Republic of Korea
| | - Hyun Ha Park
- Department of Biochemistry, College of Medicine, Konyang University, Daejeon 35365, the Republic of Korea
| | - Byeong-Hyeon Kim
- Department of Biochemistry, College of Medicine, Konyang University, Daejeon 35365, the Republic of Korea
| | - Yong Ho Park
- Department of Biochemistry, College of Medicine, Konyang University, Daejeon 35365, the Republic of Korea
| | - Eunbi Cho
- Department of Pharmacology, School of Medicine, Konkuk University, Seoul 05029, the Republic of Korea
| | - William A Goddard
- Division of Chemistry and Chemical Engineering and Materials Process and Simulation Center, California Institute of Technology, Pasadena, CA 91125, United States
| | - Dong Hyun Kim
- Department of Pharmacology, School of Medicine, Konkuk University, Seoul 05029, the Republic of Korea.
| | - Hugh I Kim
- Department of Chemistry, Korea University, Seoul 02841, the Republic of Korea; Center for Proteogenome Research, Korea University, Seoul 02841, the Republic of Korea; Single Cell Analysis Laboratory, Korea University, Seoul 02841, the Republic of Korea; Division of Chemistry and Chemical Engineering and Materials Process and Simulation Center, California Institute of Technology, Pasadena, CA 91125, United States.
| | - Minho Moon
- Department of Biochemistry, College of Medicine, Konyang University, Daejeon 35365, the Republic of Korea; Research Institute for Dementia Science, Konyang University, Daejeon 35365, the Republic of Korea.
| |
Collapse
|
4
|
Tsirigoni AM, Goktas M, Atris Z, Valleriani A, Vila Verde A, Blank KG. Chain Sliding versus β-Sheet Formation upon Shearing Single α-Helical Coiled Coils. Macromol Biosci 2023; 23:e2200563. [PMID: 36861255 DOI: 10.1002/mabi.202200563] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Revised: 02/20/2023] [Indexed: 03/03/2023]
Abstract
Coiled coils (CCs) are key building blocks of biogenic materials and determine their mechanical response to large deformations. Of particular interest is the observation that CC-based materials display a force-induced transition from α-helices to mechanically stronger β-sheets (αβT). Steered molecular dynamics simulations predict that this αβT requires a minimum, pulling speed-dependent CC length. Here, de novo designed CCs with a length between four to seven heptads are utilized to probe if the transition found in natural CCs can be mimicked with synthetic sequences. Using single-molecule force spectroscopy and molecular dynamics simulations, these CCs are mechanically loaded in shear geometry and their rupture forces and structural responses to the applied load are determined. Simulations at the highest pulling speed (0.01 nm ns-1 ) show the appearance of β-sheet structures for the five- and six-heptad CCs and a concomitant increase in mechanical strength. The αβT is less probable at a lower pulling speed of 0.001 nm ns-1 and is not observed in force spectroscopy experiments. For CCs loaded in shear geometry, the formation of β-sheets competes with interchain sliding. β-sheet formation is only possible in higher-order CC assemblies or in tensile-loading geometries where chain sliding and dissociation are prohibited.
Collapse
Affiliation(s)
- Anna-Maria Tsirigoni
- Max Planck Institute of Colloids and Interfaces, Mechano(bio)chemistry, Am Mühlenberg 1, 14476, Potsdam, Germany.,Max Planck Institute of Colloids and Interfaces, Department of Biomaterials, Am Mühlenberg 1, 14476, Potsdam, Germany
| | - Melis Goktas
- Max Planck Institute of Colloids and Interfaces, Mechano(bio)chemistry, Am Mühlenberg 1, 14476, Potsdam, Germany
| | - Zeynep Atris
- Max Planck Institute of Colloids and Interfaces, Mechano(bio)chemistry, Am Mühlenberg 1, 14476, Potsdam, Germany.,Max Planck Institute of Colloids and Interfaces, Department of Biomaterials, Am Mühlenberg 1, 14476, Potsdam, Germany
| | - Angelo Valleriani
- Max Planck Institute of Colloids and Interfaces, Department of Biomaterials, Am Mühlenberg 1, 14476, Potsdam, Germany
| | - Ana Vila Verde
- University of Duisburg-Essen, Faculty of Physics, Lotharstrasse 1, 47057, Duisburg, Germany
| | - Kerstin G Blank
- Max Planck Institute of Colloids and Interfaces, Mechano(bio)chemistry, Am Mühlenberg 1, 14476, Potsdam, Germany.,Johannes Kepler University Linz, Institute of Experimental Physics, Department of Biomolecular & Selforganizing Matter, Altenberger Strasse 69, Linz, 4040, Austria
| |
Collapse
|
5
|
Kunche L, Natarajan U. Conformations and Solvation of Synthetic Polymers in Water by Generalized Born Implicit-Solvent Molecular Dynamics Simulations: Stereoisomers of Poly(acrylic acid) and Poly(methacrylic acid). J Phys Chem B 2023; 127:1244-1253. [PMID: 36705523 DOI: 10.1021/acs.jpcb.2c06658] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
We present the GB-OBC model as an approach for implicit-solvent MD simulations of a synthetic macromolecule in water. The model is tested and found to be successful in reproducing the chain dimensions and predicting the free energy of solvation of carboxylic acid vinyl polymers. The influence of stereochemistry and the hydrophobic nature of the polymer was investigated as a function of chain length (20 < N < 600) for poly(acrylic acid) (PAA) and poly(methacrylic acid) (PMA). The dimensionless parameters of the GB-OBC model were parameterized to be applicable to PAA and PMA. Scaling relations for chain dimensions obtained using implicit-solvent MD simulations in this study are in good agreement with those from experiments, theory of solvated chains in good solvents, and all-atom MD simulations in explicit water. Results show that ⟨Rg2⟩/NL2 is greater for the atactic chain as compared to the isotactic chain, for PAA as well as PMA. ⟨Rg2⟩/NL2 values of chains attain constancy in water for N = 200, with the values being greater for PMA. The PMA chain is conformationally more perturbed than the PAA chain, for both isotactic and atactic stereochemistry. The solvation free energy ΔGhyd of PAA and PMA in water is negative for all chain lengths (N = 20-600) and becomes more favorable with an increase in molecular weight. The ΔGhyd values for isotactic and atactic chains are identical at lower values of N but differ slightly for N > 300. Irrespective of the hydrophobic nature of the polymer, the atactic chain is thermodynamically more soluble in water as compared to the isotactic chain. The isotactic chain is less hydrophilic as compared to the atactic chain due to the closer proximity of the COOH groups along the backbone. This implicit solvent method is an effective way to accurately simulate the configurational properties and solvation of synthetic polymers in water.
Collapse
Affiliation(s)
- Lakshmikumar Kunche
- Macromolecular Modeling and Simulation Lab, Department of Chemical Engineering, Indian Institute of Technology (IIT) Madras, Chennai600036, India
| | - Upendra Natarajan
- Macromolecular Modeling and Simulation Lab, Department of Chemical Engineering, Indian Institute of Technology (IIT) Madras, Chennai600036, India
| |
Collapse
|
6
|
Bo S, Sun Q, Ning P, Yuan N, Weng Y, Liang Y, Wang H, Lu Z, Li Z, Zhao X. A novel approach to analyze the association characteristics between post-spliced introns and their corresponding mRNA. Front Genet 2023; 14:1151172. [PMID: 36923795 PMCID: PMC10008863 DOI: 10.3389/fgene.2023.1151172] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2023] [Accepted: 02/15/2023] [Indexed: 03/03/2023] Open
Abstract
Studies have shown that post-spliced introns promote cell survival when nutrients are scarce, and intron loss/gain can influence many stages of mRNA metabolism. However, few approaches are currently available to study the correlation between intron sequences and their corresponding mature mRNA sequences. Here, based on the results of the improved Smith-Waterman local alignment-based algorithm method (SW method) and binding free energy weighted local alignment algorithm method (BFE method), the optimal matched segments between introns and their corresponding mature mRNAs in Caenorhabditis elegans (C.elegans) and their relative matching frequency (RF) distributions were obtained. The results showed that although the distributions of relative matching frequencies on mRNAs obtained by the BFE method were similar to the SW method, the interaction intensity in 5'and 3'untranslated regions (UTRs) regions was weaker than the SW method. The RF distributions in the exon-exon junction regions were comparable, the effects of long and short introns on mRNA and on the five functional sites with BFE method were similar to the SW method. However, the interaction intensity in 5'and 3'UTR regions with BFE method was weaker than with SW method. Although the matching rate and length distribution shape of the optimal matched fragment were consistent with the SW method, an increase in length was observed. The matching rates and the length of the optimal matched fragments were mainly in the range of 60%-80% and 20-30bp, respectively. Although we found that there were still matching preferences in the 5'and 3'UTR regions of the mRNAs with BFE, the matching intensities were significantly lower than the matching intensities between introns and their corresponding mRNAs with SW method. Overall, our findings suggest that the interaction between introns and mRNAs results from synergism among different types of sequences during the evolutionary process.
Collapse
Affiliation(s)
- Suling Bo
- College of Computer Information, Inner Mongolia Medical University, Hohhot, China
| | - Qiuying Sun
- Department of Oncology, Inner Mongolia Cancer Hospital and The Affiliated People's Hospital of Inner Mongolia Medical University, Hohhot, China
| | - Pengfei Ning
- College of Computer Information, Inner Mongolia Medical University, Hohhot, China
| | - Ningping Yuan
- College of Computer Information, Inner Mongolia Medical University, Hohhot, China
| | - Yujie Weng
- College of Computer Information, Inner Mongolia Medical University, Hohhot, China
| | - Ying Liang
- College of Computer Information, Inner Mongolia Medical University, Hohhot, China
| | - Huitao Wang
- College of Computer Information, Inner Mongolia Medical University, Hohhot, China
| | - Zhanyuan Lu
- Inner Mongolia Academy of Agricultural and Animal Husbandry Sciences, Hohhot, China.,School of Life Science, Inner Mongolia University, Hohhot, China.,Key Laboratory of Black Soil Protection And Utilization (Hohhot), Ministry of Agriculture and Rural Affairs, Hohhot, China.,6 Inner Mongolia Key Laboratory of Degradation Farmland Ecological Restoration and Pollution Control, Hohhot, China
| | - Zhongxian Li
- College of Computer Information, Inner Mongolia Medical University, Hohhot, China
| | - Xiaoqing Zhao
- Inner Mongolia Academy of Agricultural and Animal Husbandry Sciences, Hohhot, China.,School of Life Science, Inner Mongolia University, Hohhot, China.,Key Laboratory of Black Soil Protection And Utilization (Hohhot), Ministry of Agriculture and Rural Affairs, Hohhot, China.,6 Inner Mongolia Key Laboratory of Degradation Farmland Ecological Restoration and Pollution Control, Hohhot, China
| |
Collapse
|
7
|
Arsiccio A, Ganguly P, Shea JE. A Transfer Free Energy Based Implicit Solvent Model for Protein Simulations in Solvent Mixtures: Urea-Induced Denaturation as a Case Study. J Phys Chem B 2022; 126:4472-4482. [PMID: 35679169 DOI: 10.1021/acs.jpcb.2c00889] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
We developed a method for implicit solvent molecular dynamics simulations of proteins in solvent mixtures (model with implicit solvation thermodynamics, MIST). The MIST method introduces experimental group transfer free energies to the generalized Born formulation for generating molecular trajectories without the need for developing rigorous explicit-solvent force fields for multicomponent solutions. As a test case, we studied the urea-induced denaturation of the Trp-cage miniprotein in water. We demonstrate that our method allows efficient exploration of the conformational space of the protein in only a few hundreds of nanoseconds of all-atom unbiased simulations. Furthermore, selective implementation of the transfer free energies of specific peptide groups, backbone, and side chains enables us to decouple their specific energetic contributions to the conformational changes of the protein. The approach herein developed can readily be extended to the investigation of complex matrices as well as to the characterization of protein aggregation. The MIST method is implemented in Plumed (ver. 2.8) as a separate module called SASA.
Collapse
Affiliation(s)
- Andrea Arsiccio
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, Santa Barbara, California 93106, United States
| | - Pritam Ganguly
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, Santa Barbara, California 93106, United States
| | - Joan-Emma Shea
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, Santa Barbara, California 93106, United States.,Department of Physics, University of California, Santa Barbara, Santa Barbara, California 93106, United States
| |
Collapse
|
8
|
Dingfelder F, Macocco I, Benke S, Nettels D, Faccioli P, Schuler B. Slow Escape from a Helical Misfolded State of the Pore-Forming Toxin Cytolysin A. JACS AU 2021; 1:1217-1230. [PMID: 34467360 PMCID: PMC8397351 DOI: 10.1021/jacsau.1c00175] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/20/2021] [Indexed: 05/12/2023]
Abstract
The pore-forming toxin cytolysin A (ClyA) is expressed as a large α-helical monomer that, upon interaction with membranes, undergoes a major conformational rearrangement into the protomer conformation, which then assembles into a cytolytic pore. Here, we investigate the folding kinetics of the ClyA monomer with single-molecule Förster resonance energy transfer spectroscopy in combination with microfluidic mixing, stopped-flow circular dichroism experiments, and molecular simulations. The complex folding process occurs over a broad range of time scales, from hundreds of nanoseconds to minutes. The very slow formation of the native state occurs from a rapidly formed and highly collapsed intermediate with large helical content and nonnative topology. Molecular dynamics simulations suggest pronounced non-native interactions as the origin of the slow escape from this deep trap in the free-energy surface, and a variational enhanced path-sampling approach enables a glimpse of the folding process that is supported by the experimental data.
Collapse
Affiliation(s)
- Fabian Dingfelder
- Department
of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Iuri Macocco
- Department
of Physics, Trento University, Via Sommarive 14, 38123 Povo (Trento), Italy
- SISSA, Via Bonomea 265, 34136 Trieste, Italy
| | - Stephan Benke
- Department
of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Daniel Nettels
- Department
of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Pietro Faccioli
- Department
of Physics, Trento University, Via Sommarive 14, 38123 Povo (Trento), Italy
- INFN-TIFPA, Via Sommarive 14, 38123 Povo (Trento), Italy
| | - Benjamin Schuler
- Department
of Biochemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
- Department
of Physics, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| |
Collapse
|
9
|
Kulik M, Mori T, Sugita Y. Multi-Scale Flexible Fitting of Proteins to Cryo-EM Density Maps at Medium Resolution. Front Mol Biosci 2021; 8:631854. [PMID: 33842541 PMCID: PMC8025875 DOI: 10.3389/fmolb.2021.631854] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2020] [Accepted: 01/26/2021] [Indexed: 11/13/2022] Open
Abstract
Structure determination using cryo-electron microscopy (cryo-EM) medium-resolution density maps is often facilitated by flexible fitting. Avoiding overfitting, adjusting force constants driving the structure to the density map, and emulating complex conformational transitions are major concerns in the fitting. To address them, we develop a new method based on a three-step multi-scale protocol. First, flexible fitting molecular dynamics (MD) simulations with coarse-grained structure-based force field and replica-exchange scheme between different force constants replicas are performed. Second, fitted Cα atom positions guide the all-atom structure in targeted MD. Finally, the all-atom flexible fitting refinement in implicit solvent adjusts the positions of the side chains in the density map. Final models obtained via the multi-scale protocol are significantly better resolved and more reliable in comparison with long all-atom flexible fitting simulations. The protocol is useful for multi-domain systems with intricate structural transitions as it preserves the secondary structure of single domains.
Collapse
Affiliation(s)
- Marta Kulik
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako-shi, Japan
| | - Takaharu Mori
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako-shi, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako-shi, Japan.,RIKEN Center for Computational Science, Kobe, Japan.,RIKEN Center for Biosystems Dynamics Research, Kobe, Japan
| |
Collapse
|
10
|
Damjanovic J, Miao J, Huang H, Lin YS. Elucidating Solution Structures of Cyclic Peptides Using Molecular Dynamics Simulations. Chem Rev 2021; 121:2292-2324. [PMID: 33426882 DOI: 10.1021/acs.chemrev.0c01087] [Citation(s) in RCA: 44] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Protein-protein interactions are vital to biological processes, but the shape and size of their interfaces make them hard to target using small molecules. Cyclic peptides have shown promise as protein-protein interaction modulators, as they can bind protein surfaces with high affinity and specificity. Dozens of cyclic peptides are already FDA approved, and many more are in various stages of development as immunosuppressants, antibiotics, antivirals, or anticancer drugs. However, most cyclic peptide drugs so far have been natural products or derivatives thereof, with de novo design having proven challenging. A key obstacle is structural characterization: cyclic peptides frequently adopt multiple conformations in solution, which are difficult to resolve using techniques like NMR spectroscopy. The lack of solution structural information prevents a thorough understanding of cyclic peptides' sequence-structure-function relationship. Here we review recent development and application of molecular dynamics simulations with enhanced sampling to studying the solution structures of cyclic peptides. We describe novel computational methods capable of sampling cyclic peptides' conformational space and provide examples of computational studies that relate peptides' sequence and structure to biological activity. We demonstrate that molecular dynamics simulations have grown from an explanatory technique to a full-fledged tool for systematic studies at the forefront of cyclic peptide therapeutic design.
Collapse
Affiliation(s)
- Jovan Damjanovic
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| | - Jiayuan Miao
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| | - He Huang
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| | - Yu-Shan Lin
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| |
Collapse
|
11
|
Zhou P, Wang H, Chen Z, Liu Q. Context contribution to the intermolecular recognition of human ACE2-derived peptides by SARS-CoV-2 spike protein: implications for improving the peptide affinity but not altering the peptide specificity by optimizing indirect readout. Mol Omics 2020; 17:86-94. [PMID: 33174576 DOI: 10.1039/d0mo00103a] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an etiological agent of the current rapidly growing outbreak of coronavirus disease (COVID-19), which is straining health systems around the world. Disrupting the intermolecular association of SARS-CoV-2 spike glycoprotein (S protein) with its cell surface receptor human angiotensin-converting enzyme 2 (hACE2) has been recognized as a promising therapeutic strategy against COVID-19. The association is a typical peptide-mediated interaction, where the hACE adopts an α1-helix, which can form a two-helix bundle with the α2-helix, to pack against a flat pocket on the S protein surface. Here, we demonstrate that the protein context of full-length hACE plays an essential role in supporting the hACE2 α1-helix recognition by viral S protein. Energetic analysis reveals that the α1-helical peptide (αHP) and also the two-helix bundle peptide (tBP) cannot bind effectively to S protein when they are split from the hACE protein. The context contributes moderately and considerably to the direct readout (DR) and indirect readout (IR) of peptide recognition, respectively. Dynamics simulation suggests that the two free peptides exhibit a large intrinsic disorder without the support of protein context, which would incur a considerable entropy penalty upon binding to S protein. To restore the IR effect lost by splitting peptides from hACE, we herein propose employing hydrocarbon stapling and cyclization strategies to constrain the free αHP and tBP peptides into their native ordered conformations, respectively. The stapling and cyclization are carefully designed in order to avoid influencing the peptide DR effect, which has been demonstrated to improve the peptide binding affinity (but not specificity) to S protein. The stapling/cyclization-imposed conformational constraint can effectively minimize the unfavorable IR effect (i) by reducing the peptide flexibility and entropy cost upon their binding to S protein, and (ii) by helping peptide pre-folding into their native state to facilitate the conformational selection by S protein.
Collapse
Affiliation(s)
- Peng Zhou
- Center for Informational Biology, University of Electronic Science and Technology of China (UESTC) at Qingshuihe Campus, No. 2006 Xiyuan Ave West Hi-Tech Zone, Chengdu 611731, China.
| | | | | | | |
Collapse
|
12
|
Litman J, Thiel AC, Schnieders MJ. Scalable Indirect Free Energy Method Applied to Divalent Cation-Metalloprotein Binding. J Chem Theory Comput 2019; 15:4602-4614. [DOI: 10.1021/acs.jctc.9b00147] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jacob Litman
- University of Iowa, Department of Biochemistry, 51 Newton Road, 4-403 Bowen Science Building, Iowa City, Iowa 52242, United States
| | - Andrew C. Thiel
- University of Iowa, Department of Biomedical Engineering, 103 South Capitol Street, 5601 Seamans Center for the Engineering Arts and Sciences, Iowa City, Iowa 52242, United States
| | - Michael J. Schnieders
- University of Iowa, Department of Biochemistry, 51 Newton Road, 4-403 Bowen Science Building, Iowa City, Iowa 52242, United States
- University of Iowa, Department of Biomedical Engineering, 103 South Capitol Street, 5601 Seamans Center for the Engineering Arts and Sciences, Iowa City, Iowa 52242, United States
| |
Collapse
|
13
|
Nunes-Alves A, Zuckerman DM, Arantes GM. Escape of a Small Molecule from Inside T4 Lysozyme by Multiple Pathways. Biophys J 2019. [PMID: 29539393 DOI: 10.1016/j.bpj.2018.01.014] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
The T4 lysozyme L99A mutant is often used as a model system to study small-molecule binding to proteins, but pathways for ligand entry and exit from the buried binding site and the associated protein conformational changes have not been fully resolved. Here, molecular dynamics simulations were employed to model benzene exit from its binding cavity using the weighted ensemble (WE) approach to enhance sampling of low-probability unbinding trajectories. Independent WE simulations revealed four pathways for benzene exit, which correspond to transient tunnels spontaneously formed in previous simulations of apo T4 lysozyme. Thus, benzene unbinding occurs through multiple pathways partially created by intrinsic protein structural fluctuations. Motions of several α-helices and side chains were involved in ligand escape from metastable microstates. WE simulations also provided preliminary estimates of rate constants for each exit pathway. These results complement previous works and provide a semiquantitative characterization of pathway heterogeneity for binding of small molecules to proteins.
Collapse
Affiliation(s)
- Ariane Nunes-Alves
- Department of Biochemistry, Instituto de Química, Universidade de São Paulo, São Paulo, Brazil
| | - Daniel M Zuckerman
- Department of Biomedical Engineering, School of Medicine, Oregon Health & Science University, Portland, Oregon.
| | | |
Collapse
|
14
|
Gilabert JF, Lecina D, Estrada J, Guallar V. Monte Carlo Techniques for Drug Design: The Success Case of PELE. BIOMOLECULAR SIMULATIONS IN STRUCTURE-BASED DRUG DISCOVERY 2018. [DOI: 10.1002/9783527806836.ch5] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Joan F. Gilabert
- Barcelona Supercomputing Center (BSC); Life Science Department; Jordi Girona 29 08034 Barcelona Spain
| | - Daniel Lecina
- Barcelona Supercomputing Center (BSC); Life Science Department; Jordi Girona 29 08034 Barcelona Spain
| | - Jorge Estrada
- Barcelona Supercomputing Center (BSC); Life Science Department; Jordi Girona 29 08034 Barcelona Spain
| | - Victor Guallar
- Barcelona Supercomputing Center (BSC); Life Science Department; Jordi Girona 29 08034 Barcelona Spain
- ICREA; Passeig Lluís Companys 23 08010 Barcelona Spain
| |
Collapse
|
15
|
Huang H, Simmerling C. Fast Pairwise Approximation of Solvent Accessible Surface Area for Implicit Solvent Simulations of Proteins on CPUs and GPUs. J Chem Theory Comput 2018; 14:5797-5814. [PMID: 30303377 DOI: 10.1021/acs.jctc.8b00413] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We propose a pairwise and readily parallelizable SASA-based nonpolar solvation approach for protein simulations, inspired by our previous pairwise GB polar solvation model development. In this work, we developed a novel function to estimate the atomic and molecular SASAs of proteins, which results in comparable accuracy as the LCPO algorithm in reproducing numerical icosahedral-based SASA values. Implemented in Amber software and tested on consumer GPUs, our pwSASA method reasonably reproduces LCPO simulation results, but accelerates MD simulations up to 30 times compared to the LCPO implementation, which is greatly desirable for protein simulations facing sampling challenges. The value of incorporating the nonpolar term in implicit solvent simulations is explored on a peptide fragment containing the hydrophobic core of HP36 and evaluating thermal stability profiles of four small proteins.
Collapse
Affiliation(s)
- He Huang
- Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States.,Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Carlos Simmerling
- Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States.,Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States
| |
Collapse
|
16
|
Dongmo Foumthuim CJ, Corazza A, Esposito G, Fogolari F. Molecular dynamics simulations of β2-microglobulin interaction with hydrophobic surfaces. MOLECULAR BIOSYSTEMS 2018; 13:2625-2637. [PMID: 29051937 DOI: 10.1039/c7mb00464h] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Hydrophobic surfaces are known to adsorb and unfold proteins, a process that has been studied only for a few proteins. Here we address the interaction of β2-microglobulin, a paradigmatic protein for the study of amyloidogenesis, with hydrophobic surfaces. A system with 27 copies of the protein surrounded by a model cubic hydrophobic box is studied by implicit solvent molecular dynamics simulations. Most proteins adsorb on the walls of the box without major distortions in local geometry, whereas free molecules maintain proper structures and fluctuations as observed in explicit solvent molecular dynamics simulations. The major conclusions from the simulations are as follows: (i) the adopted implicit solvent model is adequate to describe protein dynamics and thermodynamics; (ii) adsorption occurs readily and is irreversible on the simulated timescale; (iii) the regions most involved in molecular encounters and stable interactions with the walls are the same as those that are important in protein-protein and protein-nanoparticle interactions; (iv) unfolding following adsorption occurs at regions found to be flexible by both experiments and simulations; (v) thermodynamic analysis suggests a very large contribution from van der Waals interactions, whereas unfavorable electrostatic interactions are not found to contribute much to adsorption energy. Surfaces with different degrees of hydrophobicity may occur in vivo. Our simulations show that adsorption is a fast and irreversible process which is accompanied by partial unfolding. The results and the thermodynamic analysis presented here are consistent with and rationalize previous experimental work.
Collapse
|
17
|
Valverde LR, Thurston BA, Ferguson AL, Wilson WL. Evidence for Prenucleated Fibrilogenesis of Acid-Mediated Self-Assembling Oligopeptides via Molecular Simulation and Fluorescence Correlation Spectroscopy. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2018; 34:7346-7354. [PMID: 29842783 DOI: 10.1021/acs.langmuir.8b00312] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
An important step in controlling biomimetic amyloid systems is understanding the self-assembly reaction kinetics. We are interested in a family of such materials characterized by symmetric sequences of amino acids flanking a π-conjugated functional core. Many of these materials rapidly self-assemble into long fibers upon protonation in an acidic environment. Despite extensive investigation of these materials' properties, little is yet understood regarding their reaction kinetics. Based on previous studies, we have chosen DFAG-4T-GAFD as a representative system and conducted molecular dynamics simulations to show that although large-scale assembly is induced by lowering pH, some degree of assembly is thermodynamically favorable in high-pH nonprotonating environments. These results are consistent with findings for other systems such as DFAG-OPV-GAFD. The nonprotonated aggregation also appears to be concentration dependent, occurring at concentrations of 100 nM and above. Single molecule measurements using fluorescence correlation spectroscopy provide experimental support for these computational predictions. We find evidence of spontaneous aggregation in aqueous solutions of peptides with concentrations as low as 100 nM; however, 10 nM solutions appear to be largely homogeneous solutions of unassembled monomer. These results indicate that the simplest explanations for kinetics of acid-mediated assembly-protonation-induced nucleation by monomeric addition followed by subsequent stages of aggregation and elongation-are inappropriate in this system. In fact, the system only exists as pure monomer in very low concentrations, nucleation actually occurs in the absence of protonating elements at concentrations typically used for experiments, and pH triggered assembly proceeds from these preassembled aggregates. Accordingly, triggered assembly must be considered to operate outside the domain of nucleation-dependent models.
Collapse
Affiliation(s)
- Lawrence R Valverde
- Department of Materials Science and Engineering , University of Illinois at Urbana-Champaign , 1304 West Green Street , Urbana , Illinois 61801 , United States
- Frederick Seitz Materials Research Laboratory , University of Illinois at Urbana-Champaign , Urbana , Illinois 61801 , United States
| | - Bryce A Thurston
- Department of Physics , University of Illinois at Urbana-Champaign , 1110 West Green Street , Urbana , Illinois 61801 , United States
| | - Andrew L Ferguson
- Department of Materials Science and Engineering , University of Illinois at Urbana-Champaign , 1304 West Green Street , Urbana , Illinois 61801 , United States
- Frederick Seitz Materials Research Laboratory , University of Illinois at Urbana-Champaign , Urbana , Illinois 61801 , United States
- Department of Chemical and Biomolecular Engineering , University of Illinois at Urbana-Champaign , 600 South Mathews Avenue , Urbana , Illinois 61801 , United States
- Department of Physics , University of Illinois at Urbana-Champaign , 1110 West Green Street , Urbana , Illinois 61801 , United States
| | - William L Wilson
- Department of Materials Science and Engineering , University of Illinois at Urbana-Champaign , 1304 West Green Street , Urbana , Illinois 61801 , United States
- Frederick Seitz Materials Research Laboratory , University of Illinois at Urbana-Champaign , Urbana , Illinois 61801 , United States
- Center for Nanoscale Systems, Faculty of Arts and Sciences , Harvard University , Cambridge , Massachusetts 02138 , United States
| |
Collapse
|
18
|
Orioli S, A Beccara S, Faccioli P. Self-consistent calculation of protein folding pathways. J Chem Phys 2018; 147:064108. [PMID: 28810783 DOI: 10.1063/1.4997197] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
We introduce an iterative algorithm to efficiently simulate protein folding and other conformational transitions, using state-of-the-art all-atom force fields. Starting from the Langevin equation, we obtain a self-consistent stochastic equation of motion, which directly yields the reaction pathways. From the solution of this set of equations we derive a stochastic estimate of the reaction coordinate. We validate this approach against the results of plain MD simulations of the folding of a small protein, which were performed on the Anton supercomputer. In order to explore the computational efficiency of this algorithm, we apply it to generate a folding pathway of a protein that consists of 130 amino acids and has a folding rate of the order of s-1.
Collapse
Affiliation(s)
- S Orioli
- Dipartimento di Fisica, Università degli Studi di Trento, Via Sommarive 14, Povo, Trento I-38123, Italy and INFN-TIFPA, Via Sommarive 14, Povo, Trento I-38123, Italy
| | - S A Beccara
- Dipartimento di Fisica, Università degli Studi di Trento, Via Sommarive 14, Povo, Trento I-38123, Italy and INFN-TIFPA, Via Sommarive 14, Povo, Trento I-38123, Italy
| | - P Faccioli
- Dipartimento di Fisica, Università degli Studi di Trento, Via Sommarive 14, Povo, Trento I-38123, Italy and INFN-TIFPA, Via Sommarive 14, Povo, Trento I-38123, Italy
| |
Collapse
|
19
|
Thurston BA, Ferguson AL. Machine learning and molecular design of self-assembling -conjugated oligopeptides. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1469754] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Affiliation(s)
- Bryce A. Thurston
- Department of Physics, University of Illinois at Urbana-Champaign , Urbana, IL, USA
| | - Andrew L. Ferguson
- Department of Physics, University of Illinois at Urbana-Champaign , Urbana, IL, USA
- Department of Materials Science and Engineering, University of Illinois at Urbana-Champaign , Urbana, IL, USA
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign , Urbana, IL, USA
| |
Collapse
|
20
|
Folding a viral peptide in different membrane environments: pathway and sampling analyses. J Biol Phys 2018; 44:195-209. [PMID: 29644513 DOI: 10.1007/s10867-018-9490-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Accepted: 03/16/2018] [Indexed: 10/17/2022] Open
Abstract
Flock House virus (FHV) is a well-characterized model system to study infection mechanisms in non-enveloped viruses. A key stage of the infection cycle is the disruption of the endosomal membrane by a component of the FHV capsid, the membrane active γ peptide. In this study, we perform all-atom molecular dynamics simulations of the 21 N-terminal residues of the γ peptide interacting with membranes of differing compositions. We carry out umbrella sampling calculations to study the folding of the peptide to a helical state in homogenous and heterogeneous membranes consisting of neutral and anionic lipids. From the trajectory data, we evaluate folding energetics and dissect the mechanism of folding in the different membrane environments. We conclude the study by analyzing the extent of configurational sampling by performing time-lagged independent component analysis.
Collapse
|
21
|
Ianeselli A, Orioli S, Spagnolli G, Faccioli P, Cupellini L, Jurinovich S, Mennucci B. Atomic Detail of Protein Folding Revealed by an Ab Initio Reappraisal of Circular Dichroism. J Am Chem Soc 2018; 140:3674-3682. [DOI: 10.1021/jacs.7b12399] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Alan Ianeselli
- Centre for Integrative Biology, Trento University, Via Sommarive 9, 38128 Povo, Trento, Italy
| | - Simone Orioli
- Physics Department, Trento University, Via Sommarive 14, 38128 Povo, Trento, Italy
- INFN-TIFPA, Via Sommarive 14, 38128 Povo, Trento, Italy
| | - Giovanni Spagnolli
- Centre for Integrative Biology, Trento University, Via Sommarive 9, 38128 Povo, Trento, Italy
| | - Pietro Faccioli
- Physics Department, Trento University, Via Sommarive 14, 38128 Povo, Trento, Italy
- INFN-TIFPA, Via Sommarive 14, 38128 Povo, Trento, Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124, Pisa, Italy
| | - Sandro Jurinovich
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124, Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124, Pisa, Italy
| |
Collapse
|
22
|
Zhou HX, Pang X. Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation. Chem Rev 2018; 118:1691-1741. [PMID: 29319301 DOI: 10.1021/acs.chemrev.7b00305] [Citation(s) in RCA: 477] [Impact Index Per Article: 79.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.
Collapse
Affiliation(s)
- Huan-Xiang Zhou
- Department of Chemistry and Department of Physics, University of Illinois at Chicago , Chicago, Illinois 60607, United States.,Department of Physics and Institute of Molecular Biophysics, Florida State University , Tallahassee, Florida 32306, United States
| | - Xiaodong Pang
- Department of Physics and Institute of Molecular Biophysics, Florida State University , Tallahassee, Florida 32306, United States
| |
Collapse
|
23
|
Watts CR, Gregory A, Frisbie C, Lovas S. Effects of force fields on the conformational and dynamic properties of amyloid β(1-40) dimer explored by replica exchange molecular dynamics simulations. Proteins 2017; 86:279-300. [PMID: 29235155 DOI: 10.1002/prot.25439] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2017] [Accepted: 12/10/2017] [Indexed: 12/14/2022]
Abstract
The conformational space and structural ensembles of amyloid beta (Aβ) peptides and their oligomers in solution are inherently disordered and proven to be challenging to study. Optimum force field selection for molecular dynamics (MD) simulations and the biophysical relevance of results are still unknown. We compared the conformational space of the Aβ(1-40) dimers by 300 ns replica exchange MD simulations at physiological temperature (310 K) using: the AMBER-ff99sb-ILDN, AMBER-ff99sb*-ILDN, AMBER-ff99sb-NMR, and CHARMM22* force fields. Statistical comparisons of simulation results to experimental data and previously published simulations utilizing the CHARMM22* and CHARMM36 force fields were performed. All force fields yield sampled ensembles of conformations with collision cross sectional areas for the dimer that are statistically significantly larger than experimental results. All force fields, with the exception of AMBER-ff99sb-ILDN (8.8 ± 6.4%) and CHARMM36 (2.7 ± 4.2%), tend to overestimate the α-helical content compared to experimental CD (5.3 ± 5.2%). Using the AMBER-ff99sb-NMR force field resulted in the greatest degree of variance (41.3 ± 12.9%). Except for the AMBER-ff99sb-NMR force field, the others tended to under estimate the expected amount of β-sheet and over estimate the amount of turn/bend/random coil conformations. All force fields, with the exception AMBER-ff99sb-NMR, reproduce a theoretically expected β-sheet-turn-β-sheet conformational motif, however, only the CHARMM22* and CHARMM36 force fields yield results compatible with collapse of the central and C-terminal hydrophobic cores from residues 17-21 and 30-36. Although analyses of essential subspace sampling showed only minor variations between force fields, secondary structures of lowest energy conformers are different.
Collapse
Affiliation(s)
- Charles R Watts
- Department of Neurosurgery, Mayo Clinic, College of Medicine, Rochester, Minnesota.,Department of Neurosurgery, Mayo Clinic Health System, La Crosse, Wisconsin
| | - Andrew Gregory
- Department of Neurosurgery, Mayo Clinic Health System, La Crosse, Wisconsin
| | - Cole Frisbie
- Department of Neurosurgery, Mayo Clinic Health System, La Crosse, Wisconsin.,Department of Biomedical Sciences, Creighton University, Omaha, Nebraska
| | - Sándor Lovas
- Department of Biomedical Sciences, Creighton University, Omaha, Nebraska
| |
Collapse
|
24
|
Nunes-Alves A, Arantes GM. Mechanical Unfolding of Macromolecules Coupled to Bond Dissociation. J Chem Theory Comput 2017; 14:282-290. [DOI: 10.1021/acs.jctc.7b00805] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Ariane Nunes-Alves
- Department of Biochemistry,
Instituto de Quı́mica, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-900 São Paulo, SP, Brazil
| | - Guilherme Menegon Arantes
- Department of Biochemistry,
Instituto de Quı́mica, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-900 São Paulo, SP, Brazil
| |
Collapse
|
25
|
Zhang H, Jiang Y, Yan H, Cui Z, Yin C. Comparative Assessment of Computational Methods for Free Energy Calculations of Ionic Hydration. J Chem Inf Model 2017; 57:2763-2775. [DOI: 10.1021/acs.jcim.7b00485] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Haiyang Zhang
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Beijing
Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
| | - Hai Yan
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Ziheng Cui
- Beijing
Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
| | - Chunhua Yin
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
26
|
Carpinteri A, Lacidogna G, Piana G, Bassani A. Terahertz mechanical vibrations in lysozyme: Raman spectroscopy vs modal analysis. J Mol Struct 2017. [DOI: 10.1016/j.molstruc.2017.02.099] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
27
|
Mizuhara Y, Parkin D, Umezawa K, Ohnuki J, Takano M. Over-Destabilization of Protein-Protein Interaction in Generalized Born Model and Utility of Energy Density Integration Cutoff. J Phys Chem B 2017; 121:4669-4677. [PMID: 28426223 DOI: 10.1021/acs.jpcb.7b01438] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The generalize Born (GB) model is frequently used in MD simulations of biomolecular systems in aqueous solution. The GB model is usually based on the so-called Coulomb field approximation (CFA) for the energy density integration. In this study, we report that the GB model with CFA overdestabilizes the long-range electrostatic attraction between oppositely charged molecules (ionic bond forming two-helix system and kinesin-tubulin system) when the energy density integration cutoff, rmax, which is used to calculate the Born energy, is set to a large value. We show that employing large rmax, which is usually expected to make simulation results more accurate, worsens the accuracy so that the attraction is changed into repulsion. It is demonstrated that the overdestabilization is caused by the overestimation of the desolvation penalty upon binding that originates from CFA. We point out that the overdestabilization can be corrected by employing a relatively small cutoff (rmax = 10-15 Å), affirming that the GB models, even with CFA, can be used as a powerful tool to theoretically study the protein-protein interaction, particularly on its dynamical aspect, such as binding and unbinding.
Collapse
Affiliation(s)
- Yukinobu Mizuhara
- Department of Pure and Applied Physics, Waseda University , Okubo 3-4-1, Shinjuku-Ku, Tokyo 169-8555, Japan
| | - Dan Parkin
- Department of Advanced Science and Engineering, Waseda University , Okubo 3-4-1, Shinjuku-Ku, Tokyo 169-8555, Japan
| | | | - Jun Ohnuki
- Department of Pure and Applied Physics, Waseda University , Okubo 3-4-1, Shinjuku-Ku, Tokyo 169-8555, Japan
| | - Mitsunori Takano
- Department of Pure and Applied Physics, Waseda University , Okubo 3-4-1, Shinjuku-Ku, Tokyo 169-8555, Japan
| |
Collapse
|
28
|
Zhang J, Zhang H, Wu T, Wang Q, van der Spoel D. Comparison of Implicit and Explicit Solvent Models for the Calculation of Solvation Free Energy in Organic Solvents. J Chem Theory Comput 2017; 13:1034-1043. [DOI: 10.1021/acs.jctc.7b00169] [Citation(s) in RCA: 109] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jin Zhang
- Department
of Chemistry and Soft Matter Research Center, Zhejiang University, Hangzhou 310027, China
| | - Haiyang Zhang
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Tao Wu
- Department
of Chemistry and Soft Matter Research Center, Zhejiang University, Hangzhou 310027, China
| | - Qi Wang
- Department
of Chemistry and Soft Matter Research Center, Zhejiang University, Hangzhou 310027, China
| | - David van der Spoel
- Uppsala
Center for Computational Chemistry, Science for Life Laboratory, Department
of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box
596, SE-75124 Uppsala, Sweden
| |
Collapse
|
29
|
Yoo J, Aksimentiev A. Refined Parameterization of Nonbonded Interactions Improves Conformational Sampling and Kinetics of Protein Folding Simulations. J Phys Chem Lett 2016; 7:3812-3818. [PMID: 27617340 DOI: 10.1021/acs.jpclett.6b01747] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
Recent advances in computational technology have enabled brute-force molecular dynamics (MD) simulations of protein folding using physics-based molecular force fields. The extensive sampling of protein conformations afforded by such simulations revealed, however, considerable compaction of the protein conformations in the unfolded state, which is inconsistent with experiment. Here, we show that a set of surgical corrections to nonbonded interactions between amine nitrogen-carboxylate oxygen and aliphatic carbon-carbon atom pairs can considerably improve the realism of protein folding simulations. Specifically, we show that employing our corrections in ∼500 μs all-atom replica-exchange MD simulations of the WW domain and villin head piece proteins increases the size of the denatured proteins' conformations and does not destabilize the native conformations of the proteins. In addition to making the folded conformations a global minimum of the respective free energy landscapes at room temperature, our corrections also make the free energy landscape smoother, considerably accelerating the folding kinetics and, hence, reducing the computational expense of a protein folding simulation.
Collapse
Affiliation(s)
- Jejoong Yoo
- Center for the Physics of Living Cells, Department of Physics, University of Illinois at Urbana-Champaign , 1110 West Green Street, Urbana, Illinois 61801, United States
| | - Aleksei Aksimentiev
- Center for the Physics of Living Cells, Department of Physics, University of Illinois at Urbana-Champaign , 1110 West Green Street, Urbana, Illinois 61801, United States
| |
Collapse
|
30
|
Zhang H, Yin C, Yan H, van der Spoel D. Evaluation of Generalized Born Models for Large Scale Affinity Prediction of Cyclodextrin Host–Guest Complexes. J Chem Inf Model 2016; 56:2080-2092. [DOI: 10.1021/acs.jcim.6b00418] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Haiyang Zhang
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Chunhua Yin
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Hai Yan
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - David van der Spoel
- Uppsala
Center for Computational Chemistry, Science for Life Laboratory, Department
of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box
596, SE-75124 Uppsala, Sweden
| |
Collapse
|
31
|
Zwier MC, Pratt AJ, Adelman JL, Kaus JW, Zuckerman DM, Chong LT. Efficient Atomistic Simulation of Pathways and Calculation of Rate Constants for a Protein-Peptide Binding Process: Application to the MDM2 Protein and an Intrinsically Disordered p53 Peptide. J Phys Chem Lett 2016; 7:3440-5. [PMID: 27532687 PMCID: PMC5008990 DOI: 10.1021/acs.jpclett.6b01502] [Citation(s) in RCA: 70] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/19/2023]
Abstract
The characterization of protein binding processes - with all of the key conformational changes - has been a grand challenge in the field of biophysics. Here, we have used the weighted ensemble path sampling strategy to orchestrate molecular dynamics simulations, yielding atomistic views of protein-peptide binding pathways involving the MDM2 oncoprotein and an intrinsically disordered p53 peptide. A total of 182 independent, continuous binding pathways were generated, yielding a kon that is in good agreement with experiment. These pathways were generated in 15 days using 3500 cores of a supercomputer, substantially faster than would be possible with "brute force" simulations. Many of these pathways involve the anchoring of p53 residue F19 into the MDM2 binding cleft when forming the metastable encounter complex, indicating that F19 may be a kinetically important residue. Our study demonstrates that it is now practical to generate pathways and calculate rate constants for protein binding processes using atomistic simulation on typical computing resources.
Collapse
Affiliation(s)
- Matthew C. Zwier
- Department of Chemistry, Drake University, Des Moines, Iowa 50311, United States
| | - Adam J. Pratt
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Joshua L. Adelman
- Department of Biological Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Joseph W. Kaus
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Daniel M. Zuckerman
- Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Lillian T. Chong
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
- Institute of Biochemistry and Biotechnology, Martin-Luther Universität Halle-Wittenberg, Halle 06120, Germany
- Corresponding Author:
| |
Collapse
|
32
|
Wang F, Cazzolli G, Wintrode P, Faccioli P. Folding Mechanism of Proteins Im7 and Im9: Insight from All-Atom Simulations in Implicit and Explicit Solvent. J Phys Chem B 2016; 120:9297-307. [PMID: 27532482 DOI: 10.1021/acs.jpcb.6b05819] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
Im7 and Im9 are evolutionary related proteins with almost identical native structures. In spite of their structural similarity, experiments show that Im7 folds through a long-lived on-pathway intermediate, while Im9 folds according to two-state kinetics. In this work, we use a recently developed enhanced path sampling method to generate many folding trajectories for these proteins, using realistic atomistic force fields, in both implicit and explicit solvent. Overall, our results are in good agreement with the experimental ϕ values and with the result of ϕ-value-restrained molecular dynamics (MD) simulations. However, our implicit solvent simulations fail to predict a qualitative difference in the folding pathways of Im7 and Im9. In contrast, our simulations in explicit solvent correctly reproduce the fact that only protein Im7 folds through a on-pathway intermediate. By analyzing our atomistic trajectories, we provide a physical picture which explains the observed difference in the folding kinetics of these chains.
Collapse
Affiliation(s)
- F Wang
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy , Baltimore, Maryland 21201, United States
| | - G Cazzolli
- Physics Department, University of Trento , via Sommarive 14 Povo, Trento 38128, Italy
| | - P Wintrode
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy , Baltimore, Maryland 21201, United States
| | - P Faccioli
- Physics Department, University of Trento , via Sommarive 14 Povo, Trento 38128, Italy
| |
Collapse
|
33
|
Thurston BA, Tovar JD, Ferguson AL. Thermodynamics, morphology, and kinetics of early-stage self-assembly of π-conjugated oligopeptides. MOLECULAR SIMULATION 2016. [DOI: 10.1080/08927022.2015.1125997] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Bryce A. Thurston
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - John D. Tovar
- Department of Chemistry, Johns Hopkins University, Baltimore, MD, USA
- Institute of NanoBioTechnology, Johns Hopkins University, Baltimore, MD, USA
- Department of Materials Science and Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Andrew L. Ferguson
- Department of Materials Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, USA
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| |
Collapse
|
34
|
Arthur EJ, Brooks CL. Parallelization and improvements of the generalized born model with a simple sWitching function for modern graphics processors. J Comput Chem 2016; 37:927-39. [PMID: 26786647 DOI: 10.1002/jcc.24280] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2015] [Revised: 11/24/2015] [Accepted: 11/24/2015] [Indexed: 12/12/2022]
Abstract
Two fundamental challenges of simulating biologically relevant systems are the rapid calculation of the energy of solvation and the trajectory length of a given simulation. The Generalized Born model with a Simple sWitching function (GBSW) addresses these issues by using an efficient approximation of Poisson-Boltzmann (PB) theory to calculate each solute atom's free energy of solvation, the gradient of this potential, and the subsequent forces of solvation without the need for explicit solvent molecules. This study presents a parallel refactoring of the original GBSW algorithm and its implementation on newly available, low cost graphics chips with thousands of processing cores. Depending on the system size and nonbonded force cutoffs, the new GBSW algorithm offers speed increases of between one and two orders of magnitude over previous implementations while maintaining similar levels of accuracy. We find that much of the algorithm scales linearly with an increase of system size, which makes this water model cost effective for solvating large systems. Additionally, we utilize our GPU-accelerated GBSW model to fold the model system chignolin, and in doing so we demonstrate that these speed enhancements now make accessible folding studies of peptides and potentially small proteins.
Collapse
Affiliation(s)
- Evan J Arthur
- Department of Chemistry, University of Michigan, 930 N. University Ave., Ann Arbor, Michigan, 48109
| | - Charles L Brooks
- Department of Chemistry, University of Michigan, 930 N. University Ave., Ann Arbor, Michigan, 48109.,LSA Biophysics, University of Michigan, 930 N. University Ave, Ann Arbor, Michigan, 48109
| |
Collapse
|
35
|
Ulmschneider JP, Ulmschneider MB. Folding Simulations of the Transmembrane Helix of Virus Protein U in an Implicit Membrane Model. J Chem Theory Comput 2015; 3:2335-46. [PMID: 26636223 DOI: 10.1021/ct700103k] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Vpu is an 81-amino-acid auxiliary membrane protein encoded by human immunodeficiency virus type 1 (HIV-1). One of its roles is to amplify viral release by self-assembling in homo-oligomers to form functional water-filled pores enabling the flux of ions across the membrane. Various NMR and CD studies have shown that the transmembrane domain of Vpu has a helical conformation. With a recently developed implicit membrane model and an efficient Monte Carlo (MC) algorithm using concerted backbone rotations, we simulate the folding of the transmembrane domain of Vpu at atomic resolution. The implicit membrane environment is based on the generalized Born theory and enables very long time scale events, such as folding to be observed using detailed all-atom representation of the protein. Such studies are currently computationally unfeasible with fully explicit lipid bilayer molecular dynamics simulations. The correct helical transmembrane structure of Vpu is predicted from extended conformations and remains stably inserted. Tilt and kink angles agree well with experimental estimates from NMR measurements. The experimentally observed change in tilt angle in membranes of varying hydrophobic width is accurately reproduced. The extensive simulation of a pentamer of the Vpu transmembrane domain in the implicit membrane gives results similar to the ones reported previously for fully explicit bilayer simulations.
Collapse
Affiliation(s)
- Jakob P Ulmschneider
- Department of Chemistry, University of Rome "La Sapienza", Rome, Italy, and Department of Biochemistry, University of Oxford, Oxford, U.K
| | - Martin B Ulmschneider
- Department of Chemistry, University of Rome "La Sapienza", Rome, Italy, and Department of Biochemistry, University of Oxford, Oxford, U.K
| |
Collapse
|
36
|
Zhang H, Tan T, van der Spoel D. Generalized Born and Explicit Solvent Models for Free Energy Calculations in Organic Solvents: Cyclodextrin Dimerization. J Chem Theory Comput 2015; 11:5103-13. [DOI: 10.1021/acs.jctc.5b00620] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Haiyang Zhang
- Department of Biological Science and Engineering,
School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
- Beijing Key Laboratory of Bioprocess, Department of Biochemical
Engineering, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
| | - Tianwei Tan
- Beijing Key Laboratory of Bioprocess, Department of Biochemical
Engineering, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
| | - David van der Spoel
- Uppsala Center for
Computational Chemistry, Science for Life Laboratory, Department of
Cell and Molecular Biology, Uppsala University, Husargatan 3, Box
596, SE-75124 Uppsala, Sweden
| |
Collapse
|
37
|
Zachmann M, Mathias G, Antes I. Parameterization of the Hamiltonian Dielectric Solvent (HADES) Reaction-Field Method for the Solvation Free Energies of Amino Acid Side-Chain Analogs. Chemphyschem 2015; 16:1739-49. [PMID: 25820235 DOI: 10.1002/cphc.201402861] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2014] [Revised: 02/02/2015] [Indexed: 11/10/2022]
Abstract
Optimization of the Hamiltonian dielectric solvent (HADES) method for biomolecular simulations in a dielectric continuum is presented with the goal of calculating accurate absolute solvation free energies while retaining the model's accuracy in predicting conformational free-energy differences. The solvation free energies of neutral and polar amino acid side-chain analogs calculated by using HADES, which may optionally include nonpolar contributions, were optimized against experimental data to reach a chemical accuracy of about 0.5 kcal mol(-1). The new parameters were evaluated for charged side-chain analogs. The HADES results were compared with explicit-solvent, generalized Born, Poisson-Boltzmann, and QM-based methods. The potentials of mean force (PMFs) between pairs of side-chain analogs obtained by using HADES and explicit-solvent simulations were used to evaluate the effects of the improved parameters optimized for solvation free energies on intermolecular potentials.
Collapse
Affiliation(s)
- Martin Zachmann
- Theoretical Chemical Biology and Protein Modelling Group, Technische Universiät München (Germany)
| | - Gerald Mathias
- Lehrstuhl für Biomolekulare Optik, Ludwig-Maximilians Universität München (Germany).
| | - Iris Antes
- Theoretical Chemical Biology and Protein Modelling Group, Technische Universiät München (Germany).
| |
Collapse
|
38
|
Lewis I, Schaefer M, Wagner T, Oberer L, Sager E, Wipfli P, Vorherr T. A Detailed Investigation on Conformation, Permeability and PK Properties of Two Related Cyclohexapeptides. Int J Pept Res Ther 2014. [DOI: 10.1007/s10989-014-9447-3] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
39
|
Mondal J, Friesner RA, Berne BJ. Role of Desolvation in Thermodynamics and Kinetics of Ligand Binding to a Kinase. J Chem Theory Comput 2014; 10:5696-5705. [PMID: 25516727 PMCID: PMC4263462 DOI: 10.1021/ct500584n] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2014] [Indexed: 12/29/2022]
Abstract
![]()
Computer
simulations are used to determine the free energy landscape
for the binding of the anticancer drug Dasatinib to its src kinase
receptor and show that before settling into a free energy basin the
ligand must surmount a free energy barrier. An analysis based on using
both the ligand-pocket separation and the pocket-water occupancy as
reaction coordinates shows that the free energy barrier is a result
of the free energy cost for almost complete desolvation of the binding
pocket. The simulations further show that the barrier is not a result
of the reorganization free energy of the binding pocket. Although
a continuum solvent model gives the location of free energy minima,
it is not able to reproduce the intermediate free energy barrier.
Finally, it is shown that a kinetic model for the on rate constant
in which the ligand diffuses up to a doorway state and then surmounts
the desolvation free energy barrier is consistent with published microsecond
time-scale simulations of the ligand binding kinetics for this system
[Shaw, D. E. et al. J. Am.
Chem. Soc.2011, 133, 9181−918321545110].
Collapse
Affiliation(s)
- Jagannath Mondal
- Department of Chemistry, Columbia University , 3000 Broadway, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University , 3000 Broadway, New York, New York 10027, United States
| | - B J Berne
- Department of Chemistry, Columbia University , 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
40
|
Hédin F, Plattner N, Doll JD, Meuwly M. Spatial Averaging: Sampling Enhancement for Exploring Configurational Space of Atomic Clusters and Biomolecules. J Chem Theory Comput 2014; 10:4284-96. [DOI: 10.1021/ct500529w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Florent Hédin
- Department
of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Nuria Plattner
- Department
of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, D-14195 Berlin, Germany
| | - J. D. Doll
- Department
of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Markus Meuwly
- Department
of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
- Department
of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
41
|
Razzokov J, Naderi S, van der Schoot P. Prediction of the structure of a silk-like protein in oligomeric states using explicit and implicit solvent models. SOFT MATTER 2014; 10:5362-5374. [PMID: 24937549 DOI: 10.1039/c4sm00384e] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
We perform Replica Exchange Molecular Dynamics (REMD) simulations on a silk-like protein design with amino-acid sequence [(Gly-Ala)3-Gly-Glu]5 to investigate the stability of a single protein, a dimer, a trimer and a tetramer made up of these proteins starting from β-roll and β-sheet structures in both explicit (TIP3P) and implicit (GBSA) solvent models. Our simulation results for the implicit solvent model agree with those for the explicit solvent model for simulation times up to the longest tested, being 30 ns per replica. From this we infer that the implicit solvent model that we use is reliable, allowing us to reach much longer time scales (up to 200 ns per replica). We find that the self-assembly of fibers of these proteins in solution must be a nucleated process, involving nuclei made up of at least three monomers. We also find that the conformation of the protein changes upon assembly, i.e., there is a transition from a disordered globular state to an ordered β-sheet structure in the self-assembled state of aggregates containing more than two monomers. This indicates that autosteric effects must be important in the polymerization of this protein, reminiscent of what is observed for β-amyloids. Our findings are consistent with recent experimental results on a protein with an amino acid sequence similar to that of the protein we study.
Collapse
Affiliation(s)
- Jamoliddin Razzokov
- Institute Ion-Plasma and Laser Technologies, Academy of Sciences of Uzbekistan, Dormon yoli Str. 33, 100125, Tashkent, Uzbekistan.
| | | | | |
Collapse
|
42
|
Li M, Duan M, Fan J, Han L, Huo S. Graph representation of protein free energy landscape. J Chem Phys 2014; 139:185101. [PMID: 24320303 DOI: 10.1063/1.4829768] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
The thermodynamics and kinetics of protein folding and protein conformational changes are governed by the underlying free energy landscape. However, the multidimensional nature of the free energy landscape makes it difficult to describe. We propose to use a weighted-graph approach to depict the free energy landscape with the nodes on the graph representing the conformational states and the edge weights reflecting the free energy barriers between the states. Our graph is constructed from a molecular dynamics trajectory and does not involve projecting the multi-dimensional free energy landscape onto a low-dimensional space defined by a few order parameters. The calculation of free energy barriers was based on transition-path theory using the MSMBuilder2 package. We compare our graph with the widely used transition disconnectivity graph (TRDG) which is constructed from the same trajectory and show that our approach gives more accurate description of the free energy landscape than the TRDG approach even though the latter can be organized into a simple tree representation. The weighted-graph is a general approach and can be used on any complex system.
Collapse
Affiliation(s)
- Minghai Li
- Gustaf H. Carlson School of Chemistry and Biochemistry, Clark University, 950 Main Street, Worcester, Massachusetts 01610, USA
| | | | | | | | | |
Collapse
|
43
|
Nunes-Alves A, Arantes GM. Ligand-receptor affinities computed by an adapted linear interaction model for continuum electrostatics and by protein conformational averaging. J Chem Inf Model 2014; 54:2309-19. [PMID: 25076043 DOI: 10.1021/ci500301s] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Accurate calculations of free energies involved in small-molecule binding to a receptor are challenging. Interactions between ligand, receptor, and solvent molecules have to be described precisely, and a large number of conformational microstates has to be sampled, particularly for ligand binding to a flexible protein. Linear interaction energy models are computationally efficient methods that have found considerable success in the prediction of binding free energies. Here, we parametrize a linear interaction model for implicit solvation with coefficients adapted by ligand and binding site relative polarities in order to predict ligand binding free energies. Results obtained for a diverse series of ligands suggest that the model has good predictive power and transferability. We also apply implicit ligand theory and propose approximations to average contributions of multiple ligand-receptor poses built from a protein conformational ensemble and find that exponential averages require proper energy discrimination between plausible binding poses and false-positives (i.e., decoys). The linear interaction model and the averaging procedures presented can be applied independently of each other and of the method used to obtain the receptor structural representation.
Collapse
Affiliation(s)
- Ariane Nunes-Alves
- Department of Biochemistry, Instituto de Química, Universidade de São Paulo , Av. Prof. Lineu Prestes 748, 05508-900, São Paulo, SP, Brazil
| | | |
Collapse
|
44
|
Yu H, Zhou P, Deng M, Shang Z. Indirect Readout in Protein-Peptide Recognition: A Different Story from Classical Biomolecular Recognition. J Chem Inf Model 2014; 54:2022-32. [DOI: 10.1021/ci5000246] [Citation(s) in RCA: 96] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
| | - Peng Zhou
- Center
of Bioinformatics (COBI), School of Life Science and Technology, University of Electronic Science and Technology of China (UESTC), Chengdu, Sichuan 610054, China
| | | | | |
Collapse
|
45
|
Duan M, Li M, Han L, Huo S. Euclidean sections of protein conformation space and their implications in dimensionality reduction. Proteins 2014; 82:2585-96. [PMID: 24913095 DOI: 10.1002/prot.24622] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2014] [Revised: 05/06/2014] [Accepted: 05/30/2014] [Indexed: 01/05/2023]
Abstract
Dimensionality reduction is widely used in searching for the intrinsic reaction coordinates for protein conformational changes. We find the dimensionality-reduction methods using the pairwise root-mean-square deviation (RMSD) as the local distance metric face a challenge. We use Isomap as an example to illustrate the problem. We believe that there is an implied assumption for the dimensionality-reduction approaches that aim to preserve the geometric relations between the objects: both the original space and the reduced space have the same kind of geometry, such as Euclidean geometry vs. Euclidean geometry or spherical geometry vs. spherical geometry. When the protein free energy landscape is mapped onto a 2D plane or 3D space, the reduced space is Euclidean, thus the original space should also be Euclidean. For a protein with N atoms, its conformation space is a subset of the 3N-dimensional Euclidean space R(3N). We formally define the protein conformation space as the quotient space of R(3N) by the equivalence relation of rigid motions. Whether the quotient space is Euclidean or not depends on how it is parameterized. When the pairwise RMSD is employed as the local distance metric, implicit representations are used for the protein conformation space, leading to no direct correspondence to a Euclidean set. We have demonstrated that an explicit Euclidean-based representation of protein conformation space and the local distance metric associated to it improve the quality of dimensionality reduction in the tetra-peptide and β-hairpin systems.
Collapse
Affiliation(s)
- Mojie Duan
- Gustaf H. Carlson School of Chemistry and Biochemistry, Clark University, Worcester, Massachusetts, 01610
| | | | | | | |
Collapse
|
46
|
Shen L, Hu H. Resolution-Adapted All-Atomic and Coarse-Grained Model for Biomolecular Simulations. J Chem Theory Comput 2014; 10:2528-36. [DOI: 10.1021/ct401029k] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Lin Shen
- Department
of Chemistry, University of Hong Kong, Pokfulam Road, Hong Kong, China
| | - Hao Hu
- Department
of Chemistry, University of Hong Kong, Pokfulam Road, Hong Kong, China
- The University of Hong Kong-Shenzhen Institute of Research and Innovation, Kejizhong Second Road, Shenzhen, China
| |
Collapse
|
47
|
Mereuta L, Roy M, Asandei A, Lee JK, Park Y, Andricioaei I, Luchian T. Slowing down single-molecule trafficking through a protein nanopore reveals intermediates for peptide translocation. Sci Rep 2014; 4:3885. [PMID: 24463372 PMCID: PMC3902492 DOI: 10.1038/srep03885] [Citation(s) in RCA: 93] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2013] [Accepted: 12/19/2013] [Indexed: 12/20/2022] Open
Abstract
The microscopic details of how peptides translocate one at a time through nanopores are crucial determinants for transport through membrane pores and important in developing nano-technologies. To date, the translocation process has been too fast relative to the resolution of the single molecule techniques that sought to detect its milestones. Using pH-tuned single-molecule electrophysiology and molecular dynamics simulations, we demonstrate how peptide passage through the α-hemolysin protein can be sufficiently slowed down to observe intermediate single-peptide sub-states associated to distinct structural milestones along the pore, and how to control residence time, direction and the sequence of spatio-temporal state-to-state dynamics of a single peptide. Molecular dynamics simulations of peptide translocation reveal the time- dependent ordering of intermediate structures of the translocating peptide inside the pore at atomic resolution. Calculations of the expected current ratios of the different pore-blocking microstates and their time sequencing are in accord with the recorded current traces.
Collapse
Affiliation(s)
- Loredana Mereuta
- Department of Physics, Alexandru I. Cuza University, Iasi, Romania
- These authors contributed equally to this work
| | - Mahua Roy
- Department of Chemistry, University of California, Irvine CA 92697, USA
- These authors contributed equally to this work
| | - Alina Asandei
- Department of Interdisciplinary Research, Alexandru I. Cuza University, Iasi, Romania
| | - Jong Kook Lee
- Research Center for Proteineous Materials, Chosun University, Gwangju, South Korea
| | - Yoonkyung Park
- Research Center for Proteineous Materials, Chosun University, Gwangju, South Korea
| | - Ioan Andricioaei
- Department of Chemistry, University of California, Irvine CA 92697, USA
| | - Tudor Luchian
- Department of Physics, Alexandru I. Cuza University, Iasi, Romania
| |
Collapse
|
48
|
Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration. ENTROPY 2013. [DOI: 10.3390/e16010163] [Citation(s) in RCA: 291] [Impact Index Per Article: 26.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
|
49
|
Co NT, Hu CK, Li MS. Dual effect of crowders on fibrillation kinetics of polypeptide chains revealed by lattice models. J Chem Phys 2013; 138:185101. [PMID: 23676073 DOI: 10.1063/1.4804299] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
We have developed the lattice model for describing polypeptide chains in the presence of crowders. The influence of crowding confinement on the fibrillation kinetics of polypeptide chains is studied using this model. We observed the non-trivial behavior of the fibril formation time τfib that it decreases with the concentration of crowders if crowder sizes are large enough, but the growth is observed for crowders of small sizes. This allows us to explain the recent experimental observation on the dual effect of crowding particles on fibril growth of proteins that for a fixed crowder concentration the fibrillation kinetics is fastest at intermediate values of total surface of crowders. It becomes slow at either small or large coverages of cosolutes. It is shown that due to competition between the energetics and entropic effects, the dependence of τfib on the size of confined space is described by a parabolic function.
Collapse
Affiliation(s)
- Nguyen Truong Co
- Institute for Computational Science and Technology, SBI Building, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam
| | | | | |
Collapse
|
50
|
Zheng G, Schaefer M, Karplus M. Hemoglobin Bohr Effects: Atomic Origin of the Histidine Residue Contributions. Biochemistry 2013; 52:8539-55. [DOI: 10.1021/bi401126z] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Guishan Zheng
- Department
of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Michael Schaefer
- Novartis Institutes for BioMedical Research, Klybeckstrasse 141, 4053 Basel, Switzerland
| | - Martin Karplus
- Department
of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
- Laboratoire
de Chimie Biophysique, ISIS Universite de Strasbourg, 67000 Strasbourg, France
| |
Collapse
|