551
|
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
|
552
|
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
|
553
|
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
|
554
|
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
|
555
|
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
|
556
|
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]
|
557
|
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
|
558
|
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
|
559
|
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
|
560
|
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
|
561
|
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
|
562
|
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
|
563
|
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
|
564
|
Carpency TN, Gunton JD, Rickman JM. Phase behavior of patchy spheroidal fluids. J Chem Phys 2016; 145:214904. [DOI: 10.1063/1.4969074] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Affiliation(s)
- T. N. Carpency
- Department of Physics, Lehigh University, Bethlehem, Pennsylvania 18015, USA
| | - J. D. Gunton
- Department of Physics, Lehigh University, Bethlehem, Pennsylvania 18015, USA
| | - J. M. Rickman
- Department of Physics, Lehigh University, Bethlehem, Pennsylvania 18015, USA
- Department of Materials Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, USA
| |
Collapse
|
565
|
Ciemny MP, Debinski A, Paczkowska M, Kolinski A, Kurcinski M, Kmiecik S. Protein-peptide molecular docking with large-scale conformational changes: the p53-MDM2 interaction. Sci Rep 2016; 6:37532. [PMID: 27905468 PMCID: PMC5131342 DOI: 10.1038/srep37532] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2016] [Accepted: 10/27/2016] [Indexed: 12/27/2022] Open
Abstract
Protein-peptide interactions are often associated with large-scale conformational changes that are difficult to study either by classical molecular modeling or by experiment. Recently, we have developed the CABS-dock method for flexible protein-peptide docking that enables large-scale rearrangements of the protein chain. In this study, we use CABS-dock to investigate the binding of the p53-MDM2 complex, an element of the cell cycle regulation system crucial for anti-cancer drug design. Experimental data suggest that p53-MDM2 binding is affected by significant rearrangements of a lid region - the N-terminal highly flexible MDM2 fragment; however, the details are not clear. The large size of the highly flexible MDM2 fragments makes p53-MDM2 intractable for exhaustive binding dynamics studies using atomistic models. We performed extensive dynamics simulations using the CABS-dock method, including large-scale structural rearrangements of MDM2 flexible regions. Without a priori knowledge of the p53 peptide structure or its binding site, we obtained near-native models of the p53-MDM2 complex. The simulation results match well the experimental data and provide new insights into the possible role of the lid fragment in p53 binding. The presented case study demonstrates that CABS-dock methodology opens up new opportunities for protein-peptide docking with large-scale changes of the protein receptor structure.
Collapse
Affiliation(s)
- Maciej Pawel Ciemny
- University of Warsaw, Faculty of Chemistry, Warsaw 02-093, Poland
- University of Warsaw, Faculty of Physics, Warsaw, 02-093, Poland
| | | | - Marta Paczkowska
- University of Warsaw, Faculty of Chemistry, Warsaw 02-093, Poland
| | - Andrzej Kolinski
- University of Warsaw, Faculty of Chemistry, Warsaw 02-093, Poland
| | | | | |
Collapse
|
566
|
Tang Y, Zhang H. Theoretical understanding of bio-interfaces/bio-surfaces by simulation: A mini review. BIOSURFACE AND BIOTRIBOLOGY 2016. [DOI: 10.1016/j.bsbt.2016.11.003] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022] Open
|
567
|
Pulawski W, Jamroz M, Kolinski M, Kolinski A, Kmiecik S. Coarse-Grained Simulations of Membrane Insertion and Folding of Small Helical Proteins Using the CABS Model. J Chem Inf Model 2016; 56:2207-2215. [PMID: 27775349 DOI: 10.1021/acs.jcim.6b00350] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
The CABS coarse-grained model is a well-established tool for modeling globular proteins (predicting their structure, dynamics, and interactions). Here we introduce an extension of the CABS representation and force field (CABS-membrane) to the modeling of the effect of the biological membrane environment on the structure of membrane proteins. We validate the CABS-membrane model in folding simulations of 10 short helical membrane proteins not using any knowledge about their structure. The simulations start from random protein conformations placed outside the membrane environment and allow for full flexibility of the modeled proteins during their spontaneous insertion into the membrane. In the resulting trajectories, we have found models close to the experimental membrane structures. We also attempted to select the correctly folded models using simple filtering followed by structural clustering combined with reconstruction to the all-atom representation and all-atom scoring. The CABS-membrane model is a promising approach for further development toward modeling of large protein-membrane systems.
Collapse
Affiliation(s)
- Wojciech Pulawski
- Faculty of Chemistry, University of Warsaw , Pasteura 1, 02-093 Warsaw, Poland
| | - Michal Jamroz
- Faculty of Chemistry, University of Warsaw , Pasteura 1, 02-093 Warsaw, Poland
| | - Michal Kolinski
- Bioinformatics Laboratory, Mossakowski Medical Research Center of the Polish Academy of Sciences , Pawinskiego 5, 02-106 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
|
568
|
Mozolewska MA, Krupa P, Zaborowski B, Liwo A, Lee J, Joo K, Czaplewski C. Use of Restraints from Consensus Fragments of Multiple Server Models To Enhance Protein-Structure Prediction Capability of the UNRES Force Field. J Chem Inf Model 2016; 56:2263-2279. [DOI: 10.1021/acs.jcim.6b00189] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
| | - Paweł Krupa
- 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
- Center
for In Silico Protein Structure and School of Computational Sciences, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu, Seoul 130-722, Republic of Korea
| | - Jooyoung Lee
- Center
for In Silico Protein Structure and School of Computational Sciences, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu, Seoul 130-722, Republic of Korea
| | - Keehyoung Joo
- Center
for Advanced Computation, Korea Institute for Advanced Study, 85
Hoegiro, Dongdaemun-gu, Seoul 130-722, Republic of Korea
| | - Cezary Czaplewski
- Faculty
of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
569
|
Multiscale design of coarse-grained elastic network-based potentials for the μ opioid receptor. J Mol Model 2016; 22:227. [DOI: 10.1007/s00894-016-3092-z] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2016] [Accepted: 08/08/2016] [Indexed: 01/10/2023]
|