501
|
Wang L, Murphy-Ullrich JE, Song Y. Multiscale simulation of the interaction of calreticulin-thrombospondin-1 complex with a model membrane microdomain. J Biomol Struct Dyn 2018; 37:811-822. [PMID: 29380675 DOI: 10.1080/07391102.2018.1433065] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Cell surface calreticulin (CRT) binding to thrombospondin-1 (TSP1), regulates cell adhesion, migration, anoikis resistance, and collagen production. Due to the essential role of membrane microdomains in CRT-mediated focal adhesion disassembly, we previously studied the effect of raft-like bilayers on TSP1-CRT interactions with all-atom molecular dynamics (AAMD) simulations. However, the simulated systems of protein on the surface of the bilayer(s) in the explicit solvent are too large for long timescale AAMD simulations due to computational expense. In this study, we adopted a multiscale modeling approach of combining AAMD, coarse-grained molecule dynamics (CGMD), and reversed AAMD (REV AAMD) simulations to investigate the interactions of single CRT or of the TSP1-CRT complex with a membrane microdomain at microsecond timescale. Results showed that CRT conformational stabilization by binding of TSP1 in AAMD simulation was undetectable in CGMD simulation, but it was recovered in REV AAMD simulation. Similarly, interactions of the CRT N-domain and TSP1 with the membrane microdomain were lost in CGMD simulations but they were re-gained in the REV AAMD simulations. There was the higher coordination of the CRT P-domain in the TSP1-CRT complex with the lipid components of membrane microdomain compared to that of single CRT, which could directly affect the conformation of CRT and further mediate CRT recruitment of LDL receptor-related protein for signaling events. This study provides structural and molecular insights into TSP1-CRT interactions in a membrane microdomain environment and demonstrates the feasibility of using multiscale simulations to investigate the interactions between protein and membrane microdomains at a long timescale.
Collapse
Affiliation(s)
- Lingyun Wang
- a Department of Biomedical Engineering , The University of Alabama at Birmingham , Birmingham 35294 , AL , USA
| | - Joanne E Murphy-Ullrich
- b Department of Pathology , The University of Alabama at Birmingham , Birmingham 35294 , AL , USA
| | - Yuhua Song
- a Department of Biomedical Engineering , The University of Alabama at Birmingham , Birmingham 35294 , AL , USA
| |
Collapse
|
502
|
Li B, Fooksa M, Heinze S, Meiler J. Finding the needle in the haystack: towards solving the protein-folding problem computationally. Crit Rev Biochem Mol Biol 2018; 53:1-28. [PMID: 28976219 PMCID: PMC6790072 DOI: 10.1080/10409238.2017.1380596] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2017] [Revised: 08/22/2017] [Accepted: 09/13/2017] [Indexed: 12/22/2022]
Abstract
Prediction of protein tertiary structures from amino acid sequence and understanding the mechanisms of how proteins fold, collectively known as "the protein folding problem," has been a grand challenge in molecular biology for over half a century. Theories have been developed that provide us with an unprecedented understanding of protein folding mechanisms. However, computational simulation of protein folding is still difficult, and prediction of protein tertiary structure from amino acid sequence is an unsolved problem. Progress toward a satisfying solution has been slow due to challenges in sampling the vast conformational space and deriving sufficiently accurate energy functions. Nevertheless, several techniques and algorithms have been adopted to overcome these challenges, and the last two decades have seen exciting advances in enhanced sampling algorithms, computational power and tertiary structure prediction methodologies. This review aims at summarizing these computational techniques, specifically conformational sampling algorithms and energy approximations that have been frequently used to study protein-folding mechanisms or to de novo predict protein tertiary structures. We hope that this review can serve as an overview on how the protein-folding problem can be studied computationally and, in cases where experimental approaches are prohibitive, help the researcher choose the most relevant computational approach for the problem at hand. We conclude with a summary of current challenges faced and an outlook on potential future directions.
Collapse
Affiliation(s)
- Bian Li
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
| | - Michaela Fooksa
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
- Chemical and Physical Biology Graduate Program, Vanderbilt University, Nashville, TN, USA
| | - Sten Heinze
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
| | - Jens Meiler
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
| |
Collapse
|
503
|
Bore SL, Milano G, Cascella M. Hybrid Particle-Field Model for Conformational Dynamics of Peptide Chains. J Chem Theory Comput 2018; 14:1120-1130. [DOI: 10.1021/acs.jctc.7b01160] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Sigbjørn Løland Bore
- Department
of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O.
Box 1033 Blindern, 0315 Oslo, Norway
| | - Giuseppe Milano
- Department
of Organic Materials Science, University of Yamagata, 4-3-16 Jonan
Yonezawa, Yamagata-ken 992-8510, Japan
| | - Michele Cascella
- Department
of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O.
Box 1033 Blindern, 0315 Oslo, Norway
| |
Collapse
|
504
|
Nasedkin A, Davidsson J, Niemi AJ, Peng X. Solution x-ray scattering and structure formation in protein dynamics. Phys Rev E 2018; 96:062405. [PMID: 29347365 DOI: 10.1103/physreve.96.062405] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2017] [Indexed: 11/07/2022]
Abstract
We propose a computationally effective approach that builds on Landau mean-field theory in combination with modern nonequilibrium statistical mechanics to model and interpret protein dynamics and structure formation in small- to wide-angle x-ray scattering (S/WAXS) experiments. We develop the methodology by analyzing experimental data in the case of Engrailed homeodomain protein as an example. We demonstrate how to interpret S/WAXS data qualitatively with a good precision and over an extended temperature range. We explain experimental observations in terms of protein phase structure, and we make predictions for future experiments and for how to analyze data at different ambient temperature values. We conclude that the approach we propose has the potential to become a highly accurate, computationally effective, and predictive tool for analyzing S/WAXS data. For this, we compare our results with those obtained previously in an all-atom molecular dynamics simulation.
Collapse
Affiliation(s)
- Alexandr Nasedkin
- Department of Physics, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden
| | - Jan Davidsson
- Department of Chemistry, Uppsala University, P. O. Box 803, S-75108, Uppsala, Sweden
| | - Antti J Niemi
- Department of Physics, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden.,Nordita, Stockholm University, Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden.,Department of Physics and Astronomy, Uppsala University, P. O. Box 803, S-75108, Uppsala, Sweden.,Laboratoire de Mathematiques et Physique Theorique CNRS UMR 6083, Fédération Denis Poisson, Université de Tours, Parc de Grandmont, F37200, Tours, France.,School of Physics, Beijing Institute of Technology, Beijing 100081, P.R. China.,Laboratory of Physics of Living Matter, School of Biomedicine, Far Eastern Federal University, Vladivostok 690090, Russia¶
| | - Xubiao Peng
- Department of Physics and Astronomy, University of British Columbia, Vancouver, British Columbia V6T1Z4, Canada
| |
Collapse
|
505
|
Hills RD. Refining amino acid hydrophobicity for dynamics simulation of membrane proteins. PeerJ 2018; 6:e4230. [PMID: 29340240 PMCID: PMC5767086 DOI: 10.7717/peerj.4230] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2017] [Accepted: 12/14/2017] [Indexed: 11/20/2022] Open
Abstract
Coarse-grained (CG) models have been successful in simulating the chemical properties of lipid bilayers, but accurate treatment of membrane proteins and lipid-protein molecular interactions remains a challenge. The CgProt force field, original developed with the multiscale coarse graining method, is assessed by comparing the potentials of mean force for sidechain insertion in a DOPC bilayer to results reported for atomistic molecular dynamics simulations. Reassignment of select CG sidechain sites from the apolar to polar site type was found to improve the attractive interfacial behavior of tyrosine, phenylalanine and asparagine as well as charged lysine and arginine residues. The solvation energy at membrane depths of 0, 1.3 and 1.7 nm correlates with experimental partition coefficients in aqueous mixtures of cyclohexane, octanol and POPC, respectively, for sidechain analogs and Wimley-White peptides. These experimental values serve as important anchor points in choosing between alternate CG models based on their observed permeation profiles, particularly for Arg, Lys and Gln residues where the all-atom OPLS solvation energy does not agree well with experiment. Available partitioning data was also used to reparameterize the representation of the peptide backbone, which needed to be made less attractive for the bilayer hydrophobic core region. The newly developed force field, CgProt 2.4, correctly predicts the global energy minimum in the potentials of mean force for insertion of the uncharged membrane-associated peptides LS3 and WALP23. CgProt will find application in studies of lipid-protein interactions and the conformational properties of diverse membrane protein systems.
Collapse
Affiliation(s)
- Ronald D Hills
- Department of Pharmaceutical Sciences, College of Pharmacy, University of New England, Portland, ME, United States of America
| |
Collapse
|
506
|
Otero-Mato JM, Montes-Campos H, Calvelo M, García-Fandiño R, Gallego LJ, Piñeiro Á, Varela LM. GADDLE Maps: General Algorithm for Discrete Object Deformations Based on Local Exchange Maps. J Chem Theory Comput 2018; 14:466-478. [DOI: 10.1021/acs.jctc.7b00861] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Affiliation(s)
- J. Manuel Otero-Mato
- Nanomaterials,
Photonics and Soft Matter Group, Departamento de Física de
Partículas y Departamento de Física Aplicada, Facultade
de Física, Universidade de Santiago de Compostela, Campus Vida s/n, E-15782 Santiago de Compostela, Spain
| | - Hadrián Montes-Campos
- Nanomaterials,
Photonics and Soft Matter Group, Departamento de Física de
Partículas y Departamento de Física Aplicada, Facultade
de Física, Universidade de Santiago de Compostela, Campus Vida s/n, E-15782 Santiago de Compostela, Spain
| | - Martín Calvelo
- Department
of Organic Chemistry, Center for Research in Biological Chemistry
and Molecular Materials, University of Santiago de Compostela, Campus
Vida s/n, E-15782 Santiago de Compostela, Spain
| | - Rebeca García-Fandiño
- CIQUP,
Department of Chemistry and Biochemistry, Faculty of Sciences, University of Porto, R. Campo alegre, 687, P-4169-007 Porto, Portugal
| | - Luis J. Gallego
- Nanomaterials,
Photonics and Soft Matter Group, Departamento de Física de
Partículas y Departamento de Física Aplicada, Facultade
de Física, Universidade de Santiago de Compostela, Campus Vida s/n, E-15782 Santiago de Compostela, Spain
| | - Ángel Piñeiro
- Soft
Matter and Molecular Biophysics Group, Departamento de Física
Aplicada, Facultade de Física, Universidade de Santiago de Compostela, Campus Vida s/n, E-15782 Santiago de Compostela, Spain
| | - Luis M. Varela
- Nanomaterials,
Photonics and Soft Matter Group, Departamento de Física de
Partículas y Departamento de Física Aplicada, Facultade
de Física, Universidade de Santiago de Compostela, Campus Vida s/n, E-15782 Santiago de Compostela, Spain
| |
Collapse
|
507
|
Du P, Rick SW, Kumar R. Towards a coarse-grained model of the peptoid backbone: the case of N,N-dimethylacetamide. Phys Chem Chem Phys 2018; 20:23386-23396. [DOI: 10.1039/c8cp03283a] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Coarse-grained model of DMA, containing the basic motif of the peptoid backbone, based on short ranged many-body ranged interactions.
Collapse
Affiliation(s)
- Pu Du
- Department of Chemistry
- Louisiana State University
- Baton Rouge
- USA
| | - Steven W. Rick
- Department of Chemistry
- University of New Orleans
- New Orleans
- USA
| | - Revati Kumar
- Department of Chemistry
- Louisiana State University
- Baton Rouge
- USA
| |
Collapse
|
508
|
Wang L, Yuan Y, Chen X, Chen J, Guo Y, Li M, Li C, Pu X. Probing the cooperative mechanism of the μ–δ opioid receptor heterodimer by multiscale simulation. Phys Chem Chem Phys 2018; 20:29969-29982. [DOI: 10.1039/c8cp06652c] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The activation-cooperativity of the μ–δ opioid receptor heterodimer was probed by multiscale simulation coupled with a protein structure network.
Collapse
Affiliation(s)
- Longrong Wang
- Faculty of Chemistry
- Sichuan University
- Chengdu
- People's Republic of China
| | - Yuan Yuan
- College of Management
- Southwest University for Nationalities
- Chengdu 610041
- P. R. China
| | - Xin Chen
- Faculty of Chemistry
- Sichuan University
- Chengdu
- People's Republic of China
| | - Jiangfan Chen
- Faculty of Chemistry
- Sichuan University
- Chengdu
- People's Republic of China
| | - Yanzhi Guo
- Faculty of Chemistry
- Sichuan University
- Chengdu
- People's Republic of China
| | - Menglong Li
- Faculty of Chemistry
- Sichuan University
- Chengdu
- People's Republic of China
| | - Chuan Li
- College of Computer Science
- Sichuan University
- Chengdu
- P. R. China
| | - Xuemei Pu
- Faculty of Chemistry
- Sichuan University
- Chengdu
- People's Republic of China
| |
Collapse
|
509
|
Salmon L, Ahlstrom LS, Bardwell JCA, Horowitz S. Selecting Conformational Ensembles Using Residual Electron and Anomalous Density (READ). Methods Mol Biol 2018; 1764:491-504. [PMID: 29605935 PMCID: PMC6148353 DOI: 10.1007/978-1-4939-7759-8_31] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Heterogeneous and dynamic biomolecular complexes play a central role in many cellular processes but are poorly understood due to experimental challenges in characterizing their structural ensembles. To address these difficulties, we developed a hybrid methodology that combines X-ray crystallography with ensemble selections typically used in NMR studies to determine structural ensembles of heterogeneous biomolecular complexes. The method, termed READ, for residual electron and anomalous density, enables the visualization of heterogeneous conformational ensembles of complexes within crystals. Here we present a detailed protocol for performing the ensemble selections to construct READ ensembles. From a diverse pool of binding poses, a selection scheme is used to determine a subset of conformations that maximizes agreement with the X-ray data. Overall, READ is a general approach for obtaining a high-resolution view of dynamic protein-protein complexes.
Collapse
Affiliation(s)
- Loïc Salmon
- Centre de RMN à Très Hauts Champs, Institut des Sciences Analytiques, UMR 5280, CNRS, ENS Lyon, UCB Lyon 1, Université de Lyon, Villeurbanne, France.
| | - Logan S Ahlstrom
- Department of Molecular, Cellular, and Developmental Biology, Howard Hughes Medical Institute, University of Michigan, Ann Arbor, MI, USA
- Department of Molecular, Cellular, and Developmental Biology, University of Michigan, Ann Arbor, MI, USA
| | - James C A Bardwell
- Department of Molecular, Cellular, and Developmental Biology, Howard Hughes Medical Institute, University of Michigan, Ann Arbor, MI, USA
- Department of Molecular, Cellular, and Developmental Biology, University of Michigan, Ann Arbor, MI, USA
| | - Scott Horowitz
- Department of Chemistry and Biochemistry, Knoebel Institute for Healthy Aging, University of Denver, Denver, CO, USA.
| |
Collapse
|
510
|
Computational Methods for Efficient Sampling of Protein Landscapes and Disclosing Allosteric Regions. COMPUTATIONAL MOLECULAR MODELLING IN STRUCTURAL BIOLOGY 2018; 113:33-63. [DOI: 10.1016/bs.apcsb.2018.06.001] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
511
|
Liu X, Chen J. HyRes: a coarse-grained model for multi-scale enhanced sampling of disordered protein conformations. Phys Chem Chem Phys 2017; 19:32421-32432. [PMID: 29186229 PMCID: PMC5729119 DOI: 10.1039/c7cp06736d] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Efficient coarse-grained (CG) models can be coupled with atomistic force fields to accelerate the sampling of atomistic energy landscapes in the multi-scale enhanced sampling (MSES) framework. This approach may be particularly suitable for generating atomistic conformational ensembles of intrinsically disordered proteins (IDPs). While MSES is relatively robust to inherent CG artifacts, achieving optimal sampling efficiency requires CG modeling to generate the local and long-range fluctuations that are largely consistent with those at the atomistic level. Here, we describe a new hybrid resolution CG model (HyRes) for MSES simulations of disordered protein states, which is specifically designed to provide semi-quantitative secondary structure propensities together with a qualitative description of long-range nonspecific interactions. The HyRes model contains an atomistic description of the backbone with intermediate resolution side chains. The secondary structure propensities are tuned by adjusting the backbone hydrogen-bonding strength and the ϕ/ψ torsion profile. The sizes and covalent geometries of the side chains are parameterized to reproduce distributions derived from atomistic simulations. Lennard-Jones parameters for sidechain beads are assigned to reproduce statistical potentials derived from the protein structural database, and then globally parameterized with nonspecific electrostatic interactions to reproduce the free energy profiles of pair wise interactions and the key conformational properties of model peptides. Application of HyRes to MSES simulations of small IDPs suggests that it is capable of driving faster structural transitions at the atomistic level and increasing the convergence rate compared to the Cα-only Gō-like models previously utilized. With further optimization, we believe that the new CG model could greatly improve the efficiency of MSES simulations of the larger and more complex IDPs frequently involved in cellular signalling and regulation.
Collapse
Affiliation(s)
- Xiaorong Liu
- Department of Chemistry, University of Massachusetts Amherst, Amherst, MA 01003, USA.
| | | |
Collapse
|
512
|
Meyer H, Voigtmann T, Schilling T. On the non-stationary generalized Langevin equation. J Chem Phys 2017; 147:214110. [DOI: 10.1063/1.5006980] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Hugues Meyer
- Physikalisches Institut, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
- Research Unit in Engineering Science, Université du Luxembourg, L-4364 Esch-sur-Alzette, Luxembourg
| | - Thomas Voigtmann
- Institut für Materialphysik im Weltraum, Deutsches Zentrum für Luft- und Raumfahrt (DLR), 51170 Köln, Germany
- Department of Physics, Heinrich Heine University, Universitätsstraße 1, 40225 Düsseldorf, Germany
| | - Tanja Schilling
- Physikalisches Institut, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
| |
Collapse
|
513
|
Trovato F, Fumagalli G. Molecular simulations of cellular processes. Biophys Rev 2017; 9:941-958. [PMID: 29185136 DOI: 10.1007/s12551-017-0363-6] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2017] [Accepted: 11/19/2017] [Indexed: 12/12/2022] Open
Abstract
It is, nowadays, possible to simulate biological processes in conditions that mimic the different cellular compartments. Several groups have performed these calculations using molecular models that vary in performance and accuracy. In many cases, the atomistic degrees of freedom have been eliminated, sacrificing both structural complexity and chemical specificity to be able to explore slow processes. In this review, we will discuss the insights gained from computer simulations on macromolecule diffusion, nuclear body formation, and processes involving the genetic material inside cell-mimicking spaces. We will also discuss the challenges to generate new models suitable for the simulations of biological processes on a cell scale and for cell-cycle-long times, including non-equilibrium events such as the co-translational folding, misfolding, and aggregation of proteins. A prominent role will be played by the wise choice of the structural simplifications and, simultaneously, of a relatively complex energetic description. These challenging tasks will rely on the integration of experimental and computational methods, achieved through the application of efficient algorithms. Graphical abstract.
Collapse
Affiliation(s)
- Fabio Trovato
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195, Berlin, Germany.
| | - Giordano Fumagalli
- Nephrology and Dialysis Unit, USL Toscana Nord Ovest, 55041, Lido di Camaiore, Lucca, Italy
| |
Collapse
|
514
|
Dosztányi Z. Prediction of protein disorder based on IUPred. Protein Sci 2017; 27:331-340. [PMID: 29076577 DOI: 10.1002/pro.3334] [Citation(s) in RCA: 117] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2017] [Revised: 10/25/2017] [Accepted: 10/25/2017] [Indexed: 12/19/2022]
Abstract
Many proteins contain intrinsically disordered regions (IDRs), functional polypeptide segments that in isolation adopt a highly flexible conformational ensemble instead of a single, well-defined structure. Disorder prediction methods, which can discriminate ordered and disordered regions from the amino acid sequence, have contributed significantly to our current understanding of the distinct properties of intrinsically disordered proteins by enabling the characterization of individual examples as well as large-scale analyses of these protein regions. One popular method, IUPred provides a robust prediction of protein disorder based on an energy estimation approach that captures the fundamental difference between the biophysical properties of ordered and disordered regions. This paper reviews the energy estimation method underlying IUPred and the basic properties of the web server. Through an example, it also illustrates how the prediction output can be interpreted in a more complex case by taking into account the heterogeneous nature of IDRs. Various applications that benefited from IUPred to provide improved disorder predictions, complementing domain annotations and aiding the identification of functional short linear motifs are also described here. IUPred is freely available for noncommercial users through the web server (http://iupred.enzim.hu and http://iupred.elte.hu) . The program can also be downloaded and installed locally for large-scale analyses.
Collapse
Affiliation(s)
- Zsuzsanna Dosztányi
- MTA-ELTE Lendület Bioinformatics Research Group, Department of Biochemistry, Eötvös Loránd University, Budapest, H-1117, Hungary
| |
Collapse
|
515
|
Onufriev AV, Izadi S. Water models for biomolecular simulations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1347] [Citation(s) in RCA: 94] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Alexey V. Onufriev
- Department of Physics; Virginia Tech; Blacksburg VA USA
- Department of Computer Science; Virginia Tech; Blacksburg VA USA
- Center for Soft Matter and Biological Physics; Virginia Tech; Blacksburg VA USA
| | - Saeed Izadi
- Early Stage Pharmaceutical Development; Genentech Inc.; South San Francisco, CA USA
| |
Collapse
|
516
|
Xu G, Ma T, Zang T, Wang Q, Ma J. OPUS-CSF: A C-atom-based scoring function for ranking protein structural models. Protein Sci 2017; 27:286-292. [PMID: 29047165 PMCID: PMC5734313 DOI: 10.1002/pro.3327] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2017] [Revised: 10/14/2017] [Accepted: 10/16/2017] [Indexed: 12/12/2022]
Abstract
We report a C‐atom‐based scoring function, named OPUS‐CSF, for ranking protein structural models. Rather than using traditional Boltzmann formula, we built a scoring function (CSF score) based on the native distributions (derived from the entire PDB) of coordinate components of mainchain C (carbonyl) atoms on selected residues of peptide segments of 5, 7, 9, and 11 residues in length. In testing OPUS‐CSF on decoy recognition, it maximally recognized 257 native structures out of 278 targets in 11 commonly used decoy sets, significantly outperforming other popular all‐atom empirical potentials. The average correlation coefficient with TM‐score was also comparable with those of other potentials. OPUS‐CSF is a highly coarse‐grained scoring function, which only requires input of partial mainchain information, and very fast. Thus, it is suitable for applications at early stage of structural building.
Collapse
Affiliation(s)
- Gang Xu
- School of Life Sciences, Tsinghua University, Beijing, China
| | - Tianqi Ma
- Applied Physics Program, Rice University, Houston, Texas.,Department of Bioengineering, Rice University, Houston, Texas
| | - Tianwu Zang
- Applied Physics Program, Rice University, Houston, Texas.,Department of Bioengineering, Rice University, Houston, Texas
| | - Qinghua Wang
- Verna and Marrs McLean Department of Biochemistry and Molecular Biology, Baylor College of Medicine, One Baylor Plaza, BCM-125, Houston, Texas
| | - Jianpeng Ma
- School of Life Sciences, Tsinghua University, Beijing, China.,Applied Physics Program, Rice University, Houston, Texas.,Department of Bioengineering, Rice University, Houston, Texas.,Verna and Marrs McLean Department of Biochemistry and Molecular Biology, Baylor College of Medicine, One Baylor Plaza, BCM-125, Houston, Texas
| |
Collapse
|
517
|
Sieradzan AK, Lipska AG, Lubecka EA. Shielding effect in protein folding. J Mol Graph Model 2017; 79:118-132. [PMID: 29161634 DOI: 10.1016/j.jmgm.2017.10.018] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2017] [Revised: 10/19/2017] [Accepted: 10/20/2017] [Indexed: 01/01/2023]
Abstract
One of the most important interactions responsible for protein folding and stability are hydrogen bonds between peptide groups. There is a constant competition between the water molecules and peptide groups in a hydrogen bond formation. Also side-chains take part in this process by reducing hydration of peptide group (shielding effect) that promotes the protein folding. In this paper, a new approach to take into account a shielding effect is presented. A modification of the energy function is derived and incorporated into the UNited RESidue (UNRES) force field. Canonical Molecular Dynamics and Replica Exchange Molecular Dynamics with UNRES force field is applied to study the influence of this effect on protein structure, folding kinetics and free energy landscapes. The results of test calculations suggest that even small contribution of this effect into energy function changes force field behavior as well as speeds up the folding process significantly.
Collapse
Affiliation(s)
- Adam K Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland.
| | - Agnieszka G Lipska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Emilia A Lubecka
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; Institute of Informatics, University of Gdańsk, Wita Stwosza 57, 80-308 Gdańsk, Poland
| |
Collapse
|
518
|
Dawid AE, Gront D, Kolinski A. SURPASS Low-Resolution Coarse-Grained Protein Modeling. J Chem Theory Comput 2017; 13:5766-5779. [PMID: 28992694 DOI: 10.1021/acs.jctc.7b00642] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Coarse-grained modeling of biomolecules has a very important role in molecular biology. In this work we present a novel SURPASS (Single United Residue per Pre-Averaged Secondary Structure fragment) model of proteins that can be an interesting alternative for existing coarse-grained models. The design of the model is unique and strongly supported by the statistical analysis of structural regularities characteristic for protein systems. Coarse-graining of protein chain structures assumes a single center of interactions per residue and accounts for preaveraged effects of four adjacent residue fragments. Knowledge-based statistical potentials encode complex interaction patterns of these fragments. Using the Replica Exchange Monte Carlo sampling scheme and a generic version of the SURPASS force field we performed test simulations of a representative set of single-domain globular proteins. The method samples a significant part of conformational space and reproduces protein structures, including native-like, with surprisingly good accuracy. Future extension of the SURPASS model on large biomacromolecular systems is briefly discussed.
Collapse
Affiliation(s)
- Aleksandra E Dawid
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw , Pasteura 1, 02-093 Warsaw, Poland
| | - Dominik Gront
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw , Pasteura 1, 02-093 Warsaw, Poland
| | - Andrzej Kolinski
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw , Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
519
|
Okumura H, Higashi M, Yoshida Y, Sato H, Akiyama R. Theoretical approaches for dynamical ordering of biomolecular systems. Biochim Biophys Acta Gen Subj 2017; 1862:212-228. [PMID: 28988931 DOI: 10.1016/j.bbagen.2017.10.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2017] [Revised: 09/30/2017] [Accepted: 10/04/2017] [Indexed: 01/21/2023]
Abstract
BACKGROUND Living systems are characterized by the dynamic assembly and disassembly of biomolecules. The dynamical ordering mechanism of these biomolecules has been investigated both experimentally and theoretically. The main theoretical approaches include quantum mechanical (QM) calculation, all-atom (AA) modeling, and coarse-grained (CG) modeling. The selected approach depends on the size of the target system (which differs among electrons, atoms, molecules, and molecular assemblies). These hierarchal approaches can be combined with molecular dynamics (MD) simulation and/or integral equation theories for liquids, which cover all size hierarchies. SCOPE OF REVIEW We review the framework of quantum mechanical/molecular mechanical (QM/MM) calculations, AA MD simulations, CG modeling, and integral equation theories. Applications of these methods to the dynamical ordering of biomolecular systems are also exemplified. MAJOR CONCLUSIONS The QM/MM calculation enables the study of chemical reactions. The AA MD simulation, which omits the QM calculation, can follow longer time-scale phenomena. By reducing the number of degrees of freedom and the computational cost, CG modeling can follow much longer time-scale phenomena than AA modeling. Integral equation theories for liquids elucidate the liquid structure, for example, whether the liquid follows a radial distribution function. GENERAL SIGNIFICANCE These theoretical approaches can analyze the dynamic behaviors of biomolecular systems. They also provide useful tools for exploring the dynamic ordering systems of biomolecules, such as self-assembly. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.
Collapse
Affiliation(s)
- Hisashi Okumura
- Research Center for Computational Science, Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan; Department of Structural Molecular Science, The Graduate University for Advanced Studies, Okazaki, Aichi 444-8585, Japan.
| | - Masahiro Higashi
- Department of Chemistry, Biology and Marine Science, University of the Ryukyus, Nishihara, Okinawa 903-0213, Japan
| | - Yuichiro Yoshida
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan; Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Japan
| | - Ryo Akiyama
- Department of Chemistry, Kyushu University, Fukuoka 819-0395, Japan
| |
Collapse
|
520
|
Alemasov NA, Ivanisenko NV, Ivanisenko VA. Regression model for predicting pathogenic properties of SOD1 mutants based on the analysis of conformational stability and conservation of hydrogen bonds. J Mol Graph Model 2017; 77:378-385. [DOI: 10.1016/j.jmgm.2017.09.014] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2017] [Revised: 09/07/2017] [Accepted: 09/15/2017] [Indexed: 10/18/2022]
|
521
|
Karczyńska AS, Czaplewski C, Krupa P, Mozolewska MA, Joo K, Lee J, Liwo A. Ergodicity and model quality in template-restrained canonical and temperature/Hamiltonian replica exchange coarse-grained molecular dynamics simulations of proteins. J Comput Chem 2017; 38:2730-2746. [DOI: 10.1002/jcc.25070] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2017] [Revised: 07/10/2017] [Accepted: 09/01/2017] [Indexed: 01/22/2023]
Affiliation(s)
- Agnieszka S. Karczyńska
- Faculty of Chemistry; University of Gdańsk, ul. Wita Stwosza 63; Gdańsk 80-308 Poland
- Center for In Silico Protein Science; Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu; Seoul 02455 Republic of Korea
- School of Computational Sciences; Korea Institute for Advanced Study, 85 Hoegiro Dongdaemun-gu; Seoul 02455 Republic of Korea
| | - Cezary Czaplewski
- Faculty of Chemistry; University of Gdańsk, ul. Wita Stwosza 63; Gdańsk 80-308 Poland
| | - Paweł Krupa
- Faculty of Chemistry; University of Gdańsk, ul. Wita Stwosza 63; Gdańsk 80-308 Poland
- Institute of Physics, Polish Academy of Sciences, Aleja Lotników 32/46; Warsaw PL 02668 Poland
| | - Magdalena A. Mozolewska
- Faculty of Chemistry; University of Gdańsk, ul. Wita Stwosza 63; Gdańsk 80-308 Poland
- Institute of Computer Science, Polish Academy of Sciences, ul. Jana Kazimierza 5; Warsaw 01-248 Poland
| | - Keehyoung Joo
- School of Computational Sciences; Korea Institute for Advanced Study, 85 Hoegiro Dongdaemun-gu; Seoul 02455 Republic of Korea
- Center for Advanced Computation, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu; Seoul 02455 Republic of Korea
| | - Jooyoung Lee
- Center for In Silico Protein Science; Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu; Seoul 02455 Republic of Korea
- School of Computational Sciences; Korea Institute for Advanced Study, 85 Hoegiro Dongdaemun-gu; Seoul 02455 Republic of Korea
- Center for Advanced Computation, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu; Seoul 02455 Republic of Korea
| | - Adam Liwo
- Faculty of Chemistry; University of Gdańsk, ul. Wita Stwosza 63; Gdańsk 80-308 Poland
- Center for In Silico Protein Science; Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu; Seoul 02455 Republic of Korea
- School of Computational Sciences; Korea Institute for Advanced Study, 85 Hoegiro Dongdaemun-gu; Seoul 02455 Republic of Korea
| |
Collapse
|
522
|
Lubecka EA, Liwo A. A general method for the derivation of the functional forms of the effective energy terms in coarse-grained energy functions of polymers. II. Backbone-local potentials of coarse-grained O1→4-bonded polyglucose chains. J Chem Phys 2017; 147:115101. [DOI: 10.1063/1.4994130] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Affiliation(s)
- Emilia A. Lubecka
- Institute of Informatics, University of Gdańsk, Wita Stwosza 57, 80-308 Gdańsk, Poland
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
523
|
Machado MR, González HC, Pantano S. MD Simulations of Viruslike Particles with Supra CG Solvation Affordable to Desktop Computers. J Chem Theory Comput 2017; 13:5106-5116. [DOI: 10.1021/acs.jctc.7b00659] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Matı́as R. Machado
- Biomolecular Simulations
Group, Institut Pasteur de Montevideo, Mataojo 2020, Montevideo CP 11400, Uruguay
| | - Humberto C. González
- Biomolecular Simulations
Group, Institut Pasteur de Montevideo, Mataojo 2020, Montevideo CP 11400, Uruguay
| | - Sergio Pantano
- Biomolecular Simulations
Group, Institut Pasteur de Montevideo, Mataojo 2020, Montevideo CP 11400, Uruguay
| |
Collapse
|
524
|
Schiffrin B, Calabrese AN, Higgins AJ, Humes JR, Ashcroft AE, Kalli AC, Brockwell DJ, Radford SE. Effects of Periplasmic Chaperones and Membrane Thickness on BamA-Catalyzed Outer-Membrane Protein Folding. J Mol Biol 2017; 429:3776-3792. [PMID: 28919234 PMCID: PMC5692476 DOI: 10.1016/j.jmb.2017.09.008] [Citation(s) in RCA: 51] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2017] [Revised: 09/08/2017] [Accepted: 09/09/2017] [Indexed: 11/18/2022]
Abstract
The biogenesis of outer-membrane proteins (OMPs) in gram-negative bacteria involves delivery by periplasmic chaperones to the β-barrel assembly machinery (BAM), which catalyzes OMP insertion into the outer membrane. Here, we examine the effects of membrane thickness, the Escherichia coli periplasmic chaperones Skp and SurA, and BamA, the central subunit of the BAM complex, on the folding kinetics of a model OMP (tOmpA) using fluorescence spectroscopy, native mass spectrometry, and molecular dynamics simulations. We show that prefolded BamA promotes the release of tOmpA from Skp despite the nM affinity of the Skp:tOmpA complex. This activity is located in the BamA β-barrel domain, but is greater when full-length BamA is present, indicating that both the β-barrel and polypeptide transport-associated (POTRA) domains are required for maximal activity. By contrast, SurA is unable to release tOmpA from Skp, providing direct evidence against a sequential chaperone model. By varying lipid acyl chain length in synthetic liposomes we show that BamA has a greater catalytic effect on tOmpA folding in thicker bilayers, suggesting that BAM catalysis involves lowering of the kinetic barrier imposed by the hydrophobic thickness of the membrane. Consistent with this, molecular dynamics simulations reveal that increases in membrane thinning/disorder by the transmembrane domain of BamA is greatest in thicker bilayers. Finally, we demonstrate that cross-linking of the BamA barrel does not affect tOmpA folding kinetics in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) liposomes, suggesting that lateral gating of the BamA barrel and/or hybrid barrel formation is not required, at least for the assembly of a small 8-stranded OMP in vitro. Mechanisms of OMP periplasmic transport and folding by BAM are poorly understood. BamA catalyzes folding by reducing the kinetic barrier imposed by membrane thickness. BamA proteoliposomes promote folding of Skp-bound tOmpA. Lateral gating is not required for BamA-catalyzed folding of tOmpA in DMPC bilayers.
Collapse
Affiliation(s)
- Bob Schiffrin
- Astbury Centre for Structural Molecular Biology, School of Molecular and Cellular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, UK
| | - Antonio N Calabrese
- Astbury Centre for Structural Molecular Biology, School of Molecular and Cellular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, UK
| | - Anna J Higgins
- Astbury Centre for Structural Molecular Biology, School of Molecular and Cellular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, UK
| | - Julia R Humes
- Astbury Centre for Structural Molecular Biology, School of Molecular and Cellular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, UK
| | - Alison E Ashcroft
- Astbury Centre for Structural Molecular Biology, School of Molecular and Cellular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, UK
| | - Antreas C Kalli
- Astbury Centre for Structural Molecular Biology, School of Molecular and Cellular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, UK; Leeds Institute of Cancer and Pathology, University of Leeds, St. James's University Hospital, Wellcome Trust Brenner Building, Leeds LS9 7TF, UK
| | - David J Brockwell
- Astbury Centre for Structural Molecular Biology, School of Molecular and Cellular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, UK.
| | - Sheena E Radford
- Astbury Centre for Structural Molecular Biology, School of Molecular and Cellular Biology, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, UK.
| |
Collapse
|
525
|
Multi-scale simulations of biological systems using the OPEP coarse-grained model. Biochem Biophys Res Commun 2017; 498:296-304. [PMID: 28917842 DOI: 10.1016/j.bbrc.2017.08.165] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2017] [Accepted: 08/31/2017] [Indexed: 12/14/2022]
Abstract
Biomolecules are complex machines that are optimized by evolution to properly fulfill or contribute to a variety of biochemical tasks in the cellular environment. Computer simulations based on quantum mechanics and atomistic force fields have been proven to be a powerful microscope for obtaining valuable insights into many biological, physical, and chemical processes. Many interesting phenomena involve, however, a time scale and a number of degrees of freedom, notably if crowding is considered, that cannot be explored at an atomistic resolution. To bridge the gap between reality and simulation, many different advanced computational techniques and coarse-grained (CG) models have been developed. Here, we report some applications of the CG OPEP protein model to amyloid fibril formation, the response of catch-bond proteins to two types of fluid flow, and interactive simulations to fold peptides with well-defined 3D structures or with intrinsic disorder.
Collapse
|
526
|
Krupa P, Hałabis A, Żmudzińska W, Ołdziej S, Scheraga HA, Liwo A. Maximum Likelihood Calibration of the UNRES Force Field for Simulation of Protein Structure and Dynamics. J Chem Inf Model 2017; 57:2364-2377. [DOI: 10.1021/acs.jcim.7b00254] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Paweł Krupa
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
- Institute of Physics, Polish Academy of Sciences, Aleja Lotników 32/46, PL-02668 Warsaw, Poland
| | - Anna Hałabis
- Laboratory of Biopolymer
Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Abrahama 58, 80-307 Gdańsk, Poland
| | - Wioletta Żmudzińska
- Laboratory of Biopolymer
Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Abrahama 58, 80-307 Gdańsk, Poland
| | - Stanisław Ołdziej
- Laboratory of Biopolymer
Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Abrahama 58, 80-307 Gdańsk, Poland
| | - Harold A. Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
- School of Computational Sciences, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu, Seoul 130-722, Republic of Korea
| |
Collapse
|
527
|
Xu G, Ma T, Zang T, Sun W, Wang Q, Ma J. OPUS-DOSP: A Distance- and Orientation-Dependent All-Atom Potential Derived from Side-Chain Packing. J Mol Biol 2017; 429:3113-3120. [PMID: 28864201 DOI: 10.1016/j.jmb.2017.08.013] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2017] [Revised: 07/27/2017] [Accepted: 08/22/2017] [Indexed: 01/18/2023]
Abstract
We report a new distance- and orientation-dependent, all-atom statistical potential derived from side-chain packing, named OPUS-DOSP, for protein structure modeling. The framework of OPUS-DOSP is based on OPUS-PSP, previously developed by us [JMB (2008), 376, 288-301], with refinement and new features. In particular, distance or orientation contribution is considered depending on the range of contact distance. A new auxiliary function in energy function is also introduced, in addition to the traditional Boltzmann term, in order to adjust the contributions of extreme cases. OPUS-DOSP was tested on 11 decoy sets commonly used for statistical potential benchmarking. Among 278 native structures, 239 and 249 native structures were recognized by OPUS-DOSP without and with the auxiliary function, respectively. The results show that OPUS-DOSP has an increased decoy recognition capability comparing with those of other relevant potentials to date.
Collapse
Affiliation(s)
- Gang Xu
- School of Life Sciences, Tsinghua University, Beijing 100084, China
| | - Tianqi Ma
- Applied Physics Program, Rice University, Houston, TX 77005, United States; Department of Bioengineering, Rice University, Houston, TX 77005, United States
| | - Tianwu Zang
- Applied Physics Program, Rice University, Houston, TX 77005, United States; Department of Bioengineering, Rice University, Houston, TX 77005, United States
| | - Weitao Sun
- Zhou Pei-Yuan Center for Applied Mathematics, Tsinghua University, Beijing 100084, China
| | - Qinghua Wang
- Verna and Marrs Mclean Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, TX 77030, United States
| | - Jianpeng Ma
- School of Life Sciences, Tsinghua University, Beijing 100084, China; Applied Physics Program, Rice University, Houston, TX 77005, United States; Department of Bioengineering, Rice University, Houston, TX 77005, United States; Verna and Marrs Mclean Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, TX 77030, United States.
| |
Collapse
|
528
|
Kurcinski M, Blaszczyk M, Ciemny MP, Kolinski A, Kmiecik S. A protocol for CABS-dock protein-peptide docking driven by side-chain contact information. Biomed Eng Online 2017; 16:73. [PMID: 28830545 PMCID: PMC5568604 DOI: 10.1186/s12938-017-0363-6] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
Background The characterization of protein–peptide interactions is a challenge for computational molecular docking. Protein–peptide docking tools face at least two major difficulties: (1) efficient sampling of large-scale conformational changes induced by binding and (2) selection of the best models from a large set of predicted structures. In this paper, we merge an efficient sampling technique with external information about side-chain contacts to sample and select the best possible models. Methods In this paper we test a new protocol that uses information about side-chain contacts in CABS-dock protein–peptide docking. As shown in our recent studies, CABS-dock enables efficient modeling of large-scale conformational changes without knowledge about the binding site. However, the resulting set of binding sites and poses is in many cases highly diverse and difficult to score. Results As we demonstrate here, information about a single side-chain contact can significantly improve the prediction accuracy. Importantly, the imposed constraints for side-chain contacts are quite soft. Therefore, the developed protocol does not require precise contact information and ensures large-scale peptide flexibility in the broad contact area. Conclusions The demonstrated protocol provides the extension of the CABS-dock method that can be practically used in the structure prediction of protein–peptide complexes guided by the knowledge of the binding interface.
Collapse
Affiliation(s)
- Mateusz Kurcinski
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Żwirki i Wigury 101, 02-089, Warsaw, Poland
| | - Maciej Blaszczyk
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Żwirki i Wigury 101, 02-089, Warsaw, Poland
| | - Maciej Pawel Ciemny
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Żwirki i Wigury 101, 02-089, Warsaw, Poland.,Faculty of Physics, University of Warsaw, Pasteura 5, Warsaw, 02-093, Poland
| | - Andrzej Kolinski
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Żwirki i Wigury 101, 02-089, Warsaw, Poland
| | - Sebastian Kmiecik
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Żwirki i Wigury 101, 02-089, Warsaw, Poland.
| |
Collapse
|
529
|
Mishra NK, Jain A, Peter C, Verma S. Combining Experimental and Simulation Techniques to Understand Morphology Control in Pentapeptide Nanostructures. J Phys Chem B 2017; 121:8155-8161. [DOI: 10.1021/acs.jpcb.7b06005] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Narendra Kumar Mishra
- Department
of Chemistry, Center for Nanoscience and Soft Nanotechnology, Indian Institute of Technology Kanpur, Kanpur, 208016 (UP), India
| | - Alok Jain
- Fachbereich
Chemie, Universität Konstanz, D-78457 Konstanz, Germany
| | - Christine Peter
- Fachbereich
Chemie, Universität Konstanz, D-78457 Konstanz, Germany
| | - Sandeep Verma
- Department
of Chemistry, Center for Nanoscience and Soft Nanotechnology, Indian Institute of Technology Kanpur, Kanpur, 208016 (UP), India
| |
Collapse
|
530
|
Ciemny MP, Kurcinski M, Blaszczyk M, Kolinski A, Kmiecik S. Modeling EphB4-EphrinB2 protein-protein interaction using flexible docking of a short linear motif. Biomed Eng Online 2017; 16:71. [PMID: 28830442 PMCID: PMC5568603 DOI: 10.1186/s12938-017-0362-7] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Background Many protein–protein interactions are mediated by a short linear motif. Usually, amino acid sequences of those motifs are known or can be predicted. It is much harder to experimentally characterize or predict their structure in the bound form. In this work, we test a possibility of using flexible docking of a short linear motif to predict the interaction interface of the EphB4-EphrinB2 complex (a system extensively studied for its significance in tumor progression). Methods In the modeling, we only use knowledge about the motif sequence and experimental structures of EphB4-EphrinB2 complex partners. The proposed protocol enables efficient modeling of significant conformational changes in the short linear motif fragment during molecular docking simulation. For the docking simulations, we use the CABS-dock method for docking fully flexible peptides to flexible protein receptors (available as a server at http://biocomp.chem.uw.edu.pl/CABSdock/). Based on the docking result, the protein–protein complex is reconstructed and refined. Results Using this novel protocol, we obtained an accurate EphB4-EphrinB2 interaction model. Conclusions The results show that the CABS-dock method may be useful as the primary docking tool in specific protein–protein docking cases similar to EphB4-EphrinB2 complex—that is, where a short linear motif fragment can be identified.
Collapse
Affiliation(s)
- Maciej Pawel Ciemny
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Żwirki i Wigury 101, 02-089, Warsaw, Poland.,Faculty of Physics, University of Warsaw, Pasteura 5, Warsaw, Poland
| | - Mateusz Kurcinski
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Żwirki i Wigury 101, 02-089, Warsaw, Poland
| | - Maciej Blaszczyk
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Żwirki i Wigury 101, 02-089, Warsaw, Poland
| | - Andrzej Kolinski
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Żwirki i Wigury 101, 02-089, Warsaw, Poland
| | - Sebastian Kmiecik
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Żwirki i Wigury 101, 02-089, Warsaw, Poland.
| |
Collapse
|
531
|
Li X, Jia J, Mei Y, Latour RA. Molecular modeling to predict peptide accessibility for peptide-functionalized hydrogels. Biointerphases 2017; 12:031008. [PMID: 28821213 PMCID: PMC5562506 DOI: 10.1116/1.4992101] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2017] [Revised: 07/28/2017] [Accepted: 08/07/2017] [Indexed: 12/16/2022] Open
Abstract
Peptide-functionalized (PF) hydrogels are being widely investigated by the tissue engineering and regenerative medicine communities for a broad range of applications because of their unique potential to mimic the natural extracellular matrix and promote tissue regeneration. In order for these complex material systems to perform their intended bioactive function (e.g., cell signaling), the peptides that are tethered to the hydrogel matrix must be accessible at the hydrogel surface for cell-receptor binding. The factors influencing the surface accessibility of the tethered peptide mainly include the length of the tethers, the loading (i.e., concentration) of the peptide, and the association between the tethered peptide and the hydrogel matrix. In the present work, the authors developed coarse-grained molecular models based on the all-atom polymer consistent force field for a type of poly(ethylene glycol)-based PF hydrogel and conducted molecular simulations to investigate the distribution of the peptide within the hydrogel and its surface accessibility as a function of tether length and peptide concentration. The calculated results of the effects of these design parameters on the surface accessibility of the peptide agree very well with corresponding experimental measurements in which peptide accessibility was quantified by the number of cells attached to the hydrogel surface per unit area. The developed modeling methods are able to provide unique insights into the molecular behavior of PF hydrogels and the distribution of the tethered peptides, which can serve as a guide for hydrogel design optimization.
Collapse
Affiliation(s)
- Xianfeng Li
- Department of Bioengineering, Clemson University, Clemson, South Carolina 29634
| | - Jia Jia
- Department of Bioengineering, Clemson University, Clemson, South Carolina 29634
| | - Ying Mei
- Department of Bioengineering, Clemson University, Clemson, South Carolina 29634
| | - Robert A Latour
- Department of Bioengineering, Clemson University, Clemson, South Carolina 29634
| |
Collapse
|
532
|
Toward more efficient simulations of slow processes in large biomolecular systems: Comment on "Ligand diffusion in proteins via enhanced sampling in molecular dynamics" by Jakub Rydzewski and Wieslaw Nowak. Phys Life Rev 2017; 22-23:75-76. [PMID: 28781239 DOI: 10.1016/j.plrev.2017.07.003] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2017] [Accepted: 07/28/2017] [Indexed: 12/31/2022]
|
533
|
Culka M, Gisdon FJ, Ullmann GM. Computational Biochemistry-Enzyme Mechanisms Explored. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2017; 109:77-112. [PMID: 28683923 DOI: 10.1016/bs.apcsb.2017.04.004] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Understanding enzyme mechanisms is a major task to achieve in order to comprehend how living cells work. Recent advances in biomolecular research provide huge amount of data on enzyme kinetics and structure. The analysis of diverse experimental results and their combination into an overall picture is, however, often challenging. Microscopic details of the enzymatic processes are often anticipated based on several hints from macroscopic experimental data. Computational biochemistry aims at creation of a computational model of an enzyme in order to explain microscopic details of the catalytic process and reproduce or predict macroscopic experimental findings. Results of such computations are in part complementary to experimental data and provide an explanation of a biochemical process at the microscopic level. In order to evaluate the mechanism of an enzyme, a structural model is constructed which can be analyzed by several theoretical approaches. Several simulation methods can and should be combined to get a reliable picture of the process of interest. Furthermore, abstract models of biological systems can be constructed combining computational and experimental data. In this review, we discuss structural computational models of enzymatic systems. We first discuss various models to simulate enzyme catalysis. Furthermore, we review various approaches how to characterize the enzyme mechanism both qualitatively and quantitatively using different modeling approaches.
Collapse
Affiliation(s)
- Martin Culka
- Computational Biochemistry, University of Bayreuth, Bayreuth, Germany
| | - Florian J Gisdon
- Computational Biochemistry, University of Bayreuth, Bayreuth, Germany
| | | |
Collapse
|
534
|
Sequence statistics of tertiary structural motifs reflect protein stability. PLoS One 2017; 12:e0178272. [PMID: 28552940 PMCID: PMC5446159 DOI: 10.1371/journal.pone.0178272] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2017] [Accepted: 05/10/2017] [Indexed: 11/19/2022] Open
Abstract
The Protein Data Bank (PDB) has been a key resource for learning general rules of sequence-structure relationships in proteins. Quantitative insights have been gained by defining geometric descriptors of structure (e.g., distances, dihedral angles, solvent exposure, etc.) and observing their distributions and sequence preferences. Here we argue that as the PDB continues to grow, it may become unnecessary to reduce structure into a set of elementary descriptors. Instead, it could be possible to deduce quantitative sequence-structure relationships in the context of precisely-defined complex structural motifs by mining the PDB for closely matching backbone geometries. To validate this idea, we turned to the the task of predicting changes in protein stability upon amino-acid substitution—a difficult problem of broad significance. We defined non-contiguous tertiary motifs (TERMs) around a protein site of interest and extracted sequence preferences from ensembles of closely-matching substructures in the PDB to predict mutational stability changes at the site, ΔΔGm. We demonstrate that these ensemble statistics predict ΔΔGm on par with state-of-the-art statistical and machine-learning methods on large thermodynamic datasets, and outperform these, along with a leading structure-based modeling approach, when tested in the context of unbiased diverse mutations. Further, we show that the performance of the TERM-based method is directly related to the amount of available relevant structural data, automatically improving with the growing PDB. This enables a means of estimating prediction accuracy. Our results clearly demonstrate that: 1) statistics of non-contiguous structural motifs in the PDB encode fundamental sequence-structure relationships related to protein thermodynamic stability, and 2) the PDB is now large enough that such statistics are already useful in practice, with their accuracy expected to continue increasing as the database grows. These observations suggest new ways of using structural data towards addressing problems of computational structural biology.
Collapse
|
535
|
Peterson LX, Roy A, Christoffer C, Terashi G, Kihara D. Modeling disordered protein interactions from biophysical principles. PLoS Comput Biol 2017; 13:e1005485. [PMID: 28394890 PMCID: PMC5402988 DOI: 10.1371/journal.pcbi.1005485] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2016] [Revised: 04/24/2017] [Accepted: 03/29/2017] [Indexed: 12/12/2022] Open
Abstract
Disordered protein-protein interactions (PPIs), those involving a folded protein and an intrinsically disordered protein (IDP), are prevalent in the cell, including important signaling and regulatory pathways. IDPs do not adopt a single dominant structure in isolation but often become ordered upon binding. To aid understanding of the molecular mechanisms of disordered PPIs, it is crucial to obtain the tertiary structure of the PPIs. However, experimental methods have difficulty in solving disordered PPIs and existing protein-protein and protein-peptide docking methods are not able to model them. Here we present a novel computational method, IDP-LZerD, which models the conformation of a disordered PPI by considering the biophysical binding mechanism of an IDP to a structured protein, whereby a local segment of the IDP initiates the interaction and subsequently the remaining IDP regions explore and coalesce around the initial binding site. On a dataset of 22 disordered PPIs with IDPs up to 69 amino acids, successful predictions were made for 21 bound and 18 unbound receptors. The successful modeling provides additional support for biophysical principles. Moreover, the new technique significantly expands the capability of protein structure modeling and provides crucial insights into the molecular mechanisms of disordered PPIs.
Collapse
Affiliation(s)
- Lenna X. Peterson
- Department of Biological Sciences, Purdue University, West Lafayette, Indiana, United States of America
| | - Amitava Roy
- Department of Biological Sciences, Purdue University, West Lafayette, Indiana, United States of America
- Department of Medicinal Chemistry and Molecular Pharmacology, Purdue University, West Lafayette, Indiana, United States of America
- Bioinformatics and Computational Biosciences Branch, Rocky Mountain Laboratories, NIAID, National Institutes of Health, Hamilton, Montana, United States of America
| | - Charles Christoffer
- Department of Computer Science, Purdue University, West Lafayette, Indiana, United States of America
| | - Genki Terashi
- Department of Biological Sciences, Purdue University, West Lafayette, Indiana, United States of America
- School of Pharmacy, Kitasato University, Tokyo, Japan
| | - Daisuke Kihara
- Department of Biological Sciences, Purdue University, West Lafayette, Indiana, United States of America
- Department of Computer Science, Purdue University, West Lafayette, Indiana, United States of America
| |
Collapse
|
536
|
Affiliation(s)
- Adam K. Sieradzan
- Chemistry
Department, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Paweł Krupa
- Chemistry
Department, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
- Institute
of Physics, Polish Academy of Sciences, Aleja Lotnikow 32/46, PL-02668 Warsaw, Poland
| | - David J. Wales
- Department
of Chemistry, Cambridge University, Lensfield Road, Cambridge CB2 1EW, U.K
| |
Collapse
|
537
|
Di Pasquale N, Carbone P. Local and global dynamics of multi-resolved polymer chains: Effects of the interactions atoms-beads on the dynamic of the chains. J Chem Phys 2017; 146:084905. [PMID: 28249431 DOI: 10.1063/1.4976645] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- N. Di Pasquale
- School of Chemistry, The University of Manchester, Manchester M13 9PL, United Kingdom
| | - P. Carbone
- School of Chemical Engineering and Analytical Science, The University of Manchester, Manchester M13 9PL, United Kingdom
| |
Collapse
|
538
|
Koehl P. Minimum action transition paths connecting minima on an energy surface. J Chem Phys 2017; 145:184111. [PMID: 27846680 DOI: 10.1063/1.4966974] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Dynamics is essential to the biological functions of many bio-molecules, yet our knowledge of dynamics remains fragmented. Experimental techniques for studying bio-molecules either provide high resolution information on static conformations of the molecule or provide low-resolution, ensemble information that does not shed light on single molecule dynamics. In parallel, bio-molecular dynamics occur at time scale that are not yet attainable through detailed simulation methods. These limitations are especially noticeable when studying transition paths. To address this issue, we report in this paper two methods that derive meaningful trajectories for proteins between two of their conformations. The first method, MinActionPath, uses approximations of the potential energy surface for the molecule to derive an analytical solution of the equations of motion related to the concept of minimum action path. The second method, RelaxPath, follows the same principle of minimum action path but implements a more sophisticated potential, including a mixed elastic potential and a collision term to alleviate steric clashes. Using this new potential, the equations of motion cannot be solved analytically. We have introduced a relaxation method for solving those equations. We describe both the theories behind the two methods and their implementations, focusing on the specific techniques we have used that make those implementations amenable to study large molecular systems. We have illustrated the performance of RelaxPath on simple 2D systems. We have also compared MinActionPath and RelaxPath to other methods for generating transition paths on a well suited test set of large proteins, for which the end points of the trajectories as well as an intermediate conformation between those end points are known. We have shown that RelaxPath outperforms those other methods, including MinActionPath, in its ability to generate trajectories that get close to the known intermediates. We have also shown that the structures along the RelaxPath trajectories remain protein-like. Open source versions of the two programs MinActionPath and RelaxPath are available by request.
Collapse
Affiliation(s)
- Patrice Koehl
- Department of Computer Science and Genome Center, University of California, Davis, California 95616, USA
| |
Collapse
|
539
|
Koehl P, Poitevin F, Navaza R, Delarue M. The Renormalization Group and Its Applications to Generating Coarse-Grained Models of Large Biological Molecular Systems. J Chem Theory Comput 2017; 13:1424-1438. [PMID: 28170254 DOI: 10.1021/acs.jctc.6b01136] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Understanding the dynamics of biomolecules is the key to understanding their biological activities. Computational methods ranging from all-atom molecular dynamics simulations to coarse-grained normal-mode analyses based on simplified elastic networks provide a general framework to studying these dynamics. Despite recent successes in studying very large systems with up to a 100,000,000 atoms, those methods are currently limited to studying small- to medium-sized molecular systems due to computational limitations. One solution to circumvent these limitations is to reduce the size of the system under study. In this paper, we argue that coarse-graining, the standard approach to such size reduction, must define a hierarchy of models of decreasing sizes that are consistent with each other, i.e., that each model contains the information of the dynamics of its predecessor. We propose a new method, Decimate, for generating such a hierarchy within the context of elastic networks for normal-mode analysis. This method is based on the concept of the renormalization group developed in statistical physics. We highlight the details of its implementation, with a special focus on its scalability to large systems of up to millions of atoms. We illustrate its application on two large systems, the capsid of a virus and the ribosome translation complex. We show that highly decimated representations of those systems, containing down to 1% of their original number of atoms, still capture qualitatively and quantitatively their dynamics. Decimate is available as an OpenSource resource.
Collapse
Affiliation(s)
- Patrice Koehl
- Department of Computer Sciences and Genome Center, University of California, Davis , Davis, California 95616, United States
| | - Frédéric Poitevin
- Department of Structural Biology, Stanford University , Stanford, California 94305, United States.,Stanford PULSE Institute, SLAC National Accelerator Laboratory, Standford University , Menlo Park, California 94025, United States
| | - Rafael Navaza
- Platform of Crystallogenesis and Crystallography, CiTech, Institut Pasteur , 75015 Paris, France
| | - Marc Delarue
- Unité de Dynamique Structurale des Macromolécules, UMR 3528 du CNRS, Institut Pasteur , 75015 Paris, France
| |
Collapse
|
540
|
Poma AB, Cieplak M, Theodorakis PE. Combining the MARTINI and Structure-Based Coarse-Grained Approaches for the Molecular Dynamics Studies of Conformational Transitions in Proteins. J Chem Theory Comput 2017; 13:1366-1374. [PMID: 28195464 DOI: 10.1021/acs.jctc.6b00986] [Citation(s) in RCA: 112] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The application of coarse-grained (CG) models in biology is essential to access large length and time scales required for the description of many biological processes. The ELNEDIN protein model is based on the well-known MARTINI CG force-field and incorporates additionally harmonic bonds of a certain spring constant within a defined cutoff distance between pairs of residues, in order to preserve the native structure of the protein. In this case, the use of unbreakable harmonic bonds hinders the study of unfolding and folding processes. To overcome this barrier we have replaced the harmonic bonds with Lennard-Jones interactions based on the contact map of the native protein structure as is done in Go̅-like models. This model exhibits very good agreement with all-atom simulations and the ELNEDIN. Moreover, it can capture the structural motion linked to particular catalytic activity in the Man5B protein, in agreement with all-atom simulations. In addition, our model is based on the van der Waals radii, instead of a cutoff distance, which results in a smaller contact map. In conclusion, we anticipate that our model will provide further possibilities for studying biological systems based on the MARTINI CG force-field by using advanced-sampling methods, such as parallel tempering and metadynamics.
Collapse
Affiliation(s)
- Adolfo B Poma
- Institute of Physics, Polish Academy of Sciences , Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences , Al. Lotników 32/46, 02-668 Warsaw, Poland
| | | |
Collapse
|
541
|
Zhang Y, Cao Z, Zhang JZ, Xia F. Performance Comparison of Systematic Methods for Rigorous Definition of Coarse-Grained Sites of Large Biomolecules. J Chem Inf Model 2017; 57:214-222. [PMID: 28128949 DOI: 10.1021/acs.jcim.6b00683] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Construction of coarse-grained (CG) models for large biomolecules used for multiscale simulations demands a rigorous definition of CG sites for them. Several coarse-graining methods such as the simulated annealing and steepest descent (SASD) based on the essential dynamics coarse-graining (ED-CG) or the stepwise local iterative optimization (SLIO) based on the fluctuation maximization coarse-graining (FM-CG), were developed to do it. However, the practical applications of these methods such as SASD based on ED-CG are subject to limitations because they are too expensive. In this work, we extend the applicability of ED-CG by combining it with the SLIO algorithm. A comprehensive comparison of optimized results and accuracy of various algorithms based on ED-CG show that SLIO is the fastest as well as the most accurate algorithm among them. ED-CG combined with SLIO could give converged results as the number of CG sites increases, which demonstrates that it is another efficient method for coarse-graining large biomolecules. The construction of CG sites for Ras protein by using MD fluctuations demonstrates that the CG sites derived from FM-CG can reflect the fluctuation properties of secondary structures in Ras accurately.
Collapse
Affiliation(s)
- Yuwei Zhang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemistry Engineering, Xiamen University , Xiamen 361005, China.,School of Chemistry and Molecular Engineering, East China Normal University , Shanghai 200062, China
| | - Zexing Cao
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemistry Engineering, Xiamen University , Xiamen 361005, China
| | - John Zenghui Zhang
- School of Chemistry and Molecular Engineering, East China Normal University , Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry, NYU Shanghai , Shanghai 200062, China
| | - Fei Xia
- School of Chemistry and Molecular Engineering, East China Normal University , Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry, NYU Shanghai , Shanghai 200062, China
| |
Collapse
|
542
|
Abstract
Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.
Collapse
Affiliation(s)
| | - Kenneth M. Merz
- Department of Chemistry, Department of Biochemistry and Molecular Biology, and Institute of Cyber-Enabled Research, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
543
|
Zheng W, Wen H. A survey of coarse-grained methods for modeling protein conformational transitions. Curr Opin Struct Biol 2017; 42:24-30. [DOI: 10.1016/j.sbi.2016.10.008] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2016] [Revised: 10/07/2016] [Accepted: 10/10/2016] [Indexed: 01/28/2023]
|
544
|
Meng Q, Chen J. Ring-polymer molecular dynamics studies on the rate coefficient of the abstraction channel of hydrogen plus ethane, propane, and dimethyl ether. J Chem Phys 2017; 146:024108. [PMID: 28088140 DOI: 10.1063/1.4973831] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Qingyong Meng
- Department of Applied Chemistry, Northwestern Polytechnical University, Youyi West Road 127, 710072 Xi’an, China
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, 116023 Dalian, China
| | - Jun Chen
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, 116023 Dalian, China
- iChEM, College of Chemistry and Chemical Engineering, Xiamen University, Siming South Road 422, 361005 Xiamen, China
| |
Collapse
|
545
|
Makowski M, Liwo A, Scheraga HA. Simple Physics-Based Analytical Formulas for the Potentials of Mean Force of the Interaction of Amino Acid Side Chains in Water. VII. Charged-Hydrophobic/Polar and Polar-Hydrophobic/Polar Side Chains. J Phys Chem B 2017; 121:379-390. [PMID: 28000446 DOI: 10.1021/acs.jpcb.6b08541] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The physics-based potentials of side-chain-side-chain interactions corresponding to pairs composed of charged and polar, polar and polar, charged and hydrophobic, and hydrophobic and hydrophobic side chains have been determined. A total of 144 four-dimensional potentials of mean force (PMFs) of all possible pairs of molecules modeling these pairs were determined by umbrella-sampling molecular dynamics simulations in explicit water as functions of distance and orientation, and the analytical expressions were then fitted to the PMFs. Depending on the type of interacting sites, the analytical approximation to the PMF is a sum of terms corresponding to van der Waals interactions and cavity-creation involving the nonpolar sections of the side chains and van der Waals, cavity-creation, and electrostatic (charge-dipole or dipole-dipole) interaction energies and polarization energies involving the charged or polar sections of the side chains. The model used in this work reproduces all features of the interacting pairs. The UNited RESidue force field with the new side-chain-side-chain interaction potentials was preliminarily tested with the N-terminal part of the B-domain of staphylococcal protein A (PDBL 1BDD ; a three-α-helix bundle) and UPF0291 protein YnzC from Bacillus subtilis (PDB: 2HEP ; an α-helical hairpin).
Collapse
Affiliation(s)
- Mariusz Makowski
- Faculty of Chemistry, University of Gdańsk , Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk , Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Harold A Scheraga
- Baker laboratory of Chemistry and Chemical Biology, Cornell University , Ithaca, New York 14853-1301, United States
| |
Collapse
|
546
|
Fosso-Tande J, Black C, G. Aller S, Lu L, D. Hills Jr R. Simulation of lipid-protein interactions with the CgProt force field. AIMS MOLECULAR SCIENCE 2017. [DOI: 10.3934/molsci.2017.3.352] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
|
547
|
Kmiecik S, Kolinski A. One-Dimensional Structural Properties of Proteins in the Coarse-Grained CABS Model. Methods Mol Biol 2017; 1484:83-113. [PMID: 27787822 DOI: 10.1007/978-1-4939-6406-2_8] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Despite the significant increase in computational power, molecular modeling of protein structure using classical all-atom approaches remains inefficient, at least for most of the protein targets in the focus of biomedical research. Perhaps the most successful strategy to overcome the inefficiency problem is multiscale modeling to merge all-atom and coarse-grained models. This chapter describes a well-established CABS coarse-grained protein model. The CABS (C-Alpha, C-Beta, and Side chains) model assumes a 2-4 united-atom representation of amino acids, knowledge-based force field (derived from the statistical regularities seen in known protein sequences and structures) and efficient Monte Carlo sampling schemes (MC dynamics, MC replica-exchange, and combinations). A particular emphasis is given to the unique design of the CABS force-field, which is largely defined using one-dimensional structural properties of proteins, including protein secondary structure. This chapter also presents CABS-based modeling methods, including multiscale tools for de novo structure prediction, modeling of protein dynamics and prediction of protein-peptide complexes. CABS-based tools are freely available at http://biocomp.chem.uw.edu.pl/tools.
Collapse
Affiliation(s)
- Sebastian Kmiecik
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warszawa, 02-093, Poland
| | - Andrzej Kolinski
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warszawa, 02-093, Poland.
| |
Collapse
|
548
|
Ciemny MP, Kurcinski M, Kozak KJ, Kolinski A, Kmiecik S. Highly Flexible Protein-Peptide Docking Using CABS-Dock. Methods Mol Biol 2017; 1561:69-94. [PMID: 28236234 DOI: 10.1007/978-1-4939-6798-8_6] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Protein-peptide molecular docking is a difficult modeling problem. It is even more challenging when significant conformational changes that may occur during the binding process need to be predicted. In this chapter, we demonstrate the capabilities and features of the CABS-dock server for flexible protein-peptide docking. CABS-dock allows highly efficient modeling of full peptide flexibility and significant flexibility of a protein receptor. During CABS-dock docking, the peptide folding and binding process is explicitly simulated and no information about the peptide binding site or its structure is used. This chapter presents a successful CABS-dock use for docking a potentially therapeutic peptide to a protein target. Moreover, simulation contact maps, a new CABS-dock feature, are described and applied to the docking test case. Finally, a tutorial for running CABS-dock from the command line or command line scripts is provided. The CABS-dock web server is available from http://biocomp.chem.uw.edu.pl/CABSdock/ .
Collapse
Affiliation(s)
- Maciej Paweł Ciemny
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093, Warsaw, Poland
| | - Mateusz Kurcinski
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093, Warsaw, Poland
| | - Konrad Jakub Kozak
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093, Warsaw, Poland
| | - Andrzej Kolinski
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093, Warsaw, Poland
| | - Sebastian Kmiecik
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093, Warsaw, Poland.
| |
Collapse
|
549
|
Hsieh YC, Poitevin F, Delarue M, Koehl P. Comparative Normal Mode Analysis of the Dynamics of DENV and ZIKV Capsids. Front Mol Biosci 2016; 3:85. [PMID: 28083537 PMCID: PMC5187361 DOI: 10.3389/fmolb.2016.00085] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2016] [Accepted: 12/12/2016] [Indexed: 11/13/2022] Open
Abstract
Key steps in the life cycle of a virus, such as the fusion event as the virus infects a host cell and its maturation process, relate to an intricate interplay between the structure and the dynamics of its constituent proteins, especially those that define its capsid, much akin to an envelope that protects its genomic material. We present a comprehensive, comparative analysis of such interplay for the capsids of two viruses from the flaviviridae family, Dengue (DENV) and Zika (ZIKV). We use for that purpose our own software suite, DD-NMA, which is based on normal mode analysis. We describe the elements of DD-NMA that are relevant to the analysis of large systems, such as virus capsids. In particular, we introduce our implementation of simplified elastic networks and justify their parametrization. Using DD-NMA, we illustrate the importance of packing interactions within the virus capsids on the dynamics of the E proteins of DENV and ZIKV. We identify differences between the computed atomic fluctuations of the E proteins in DENV and ZIKV and relate those differences to changes observed in their high resolution structures. We conclude with a discussion on additional analyses that are needed to fully characterize the dynamics of the two viruses.
Collapse
Affiliation(s)
- Yin-Chen Hsieh
- Department of Computer Science and Genome Center, University of California, Davis Davis, CA, USA
| | - Frédéric Poitevin
- Department of Structural Biology, Stanford UniversityStanford, CA, USA; SLAC National Accelerator Laboratory, Stanford PULSE InstituteMenlo Park, CA, USA
| | - Marc Delarue
- Unit of Structural Dynamics of Macromolecules, UMR 3528 du Centre National de la Recherche Scientifique, Institut Pasteur Paris, France
| | - Patrice Koehl
- Department of Computer Science and Genome Center, University of California, Davis Davis, CA, USA
| |
Collapse
|
550
|
Fraaije JGEM, van Male J, Becherer P, Serral Gracià R. Coarse-Grained Models for Automated Fragmentation and Parametrization of Molecular Databases. J Chem Inf Model 2016; 56:2361-2377. [DOI: 10.1021/acs.jcim.6b00003] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Johannes G. E. M. Fraaije
- Leiden
Institute of Chemistry, Leiden University, Einsteinweg 55, 2300 RA Leiden, The Netherlands
- Culgi BV, Galileiweg 8, 2333 BD Leiden, The Netherlands
| | - Jan van Male
- Culgi BV, Galileiweg 8, 2333 BD Leiden, The Netherlands
| | - Paul Becherer
- Culgi BV, Galileiweg 8, 2333 BD Leiden, The Netherlands
| | | |
Collapse
|