1
|
Leśniewski M, Pyrka M, Czaplewski C, Co NT, Jiang Y, Gong Z, Tang C, Liwo A. Assessment of Two Restraint Potentials for Coarse-Grained Chemical-Cross-Link-Assisted Modeling of Protein Structures. J Chem Inf Model 2024; 64:1377-1393. [PMID: 38345917 PMCID: PMC10900291 DOI: 10.1021/acs.jcim.3c01890] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2023] [Revised: 01/20/2024] [Accepted: 01/22/2024] [Indexed: 02/27/2024]
Abstract
The influence of distance restraints from chemical cross-link mass spectroscopy (XL-MS) on the quality of protein structures modeled with the coarse-grained UNRES force field was assessed by using a protocol based on multiplexed replica exchange molecular dynamics, in which both simulated and experimental cross-link restraints were employed, for 23 small proteins. Six cross-links with upper distance boundaries from 4 Å to 12 Å (azido benzoic acid succinimide (ABAS), triazidotriazine (TATA), succinimidyldiazirine (SDA), disuccinimidyl adipate (DSA), disuccinimidyl glutarate (DSG), and disuccinimidyl suberate (BS3)) and two types of restraining potentials ((i) simple flat-bottom Lorentz-like potentials dependent on side chain distance (all cross-links) and (ii) distance- and orientation-dependent potentials determined based on molecular dynamics simulations of model systems (DSA, DSG, BS3, and SDA)) were considered. The Lorentz-like potentials with properly set parameters were found to produce a greater number of higher-quality models compared to unrestrained simulations than the MD-based potentials, because the latter can force too long distances between side chains. Therefore, the flat-bottom Lorentz-like potentials are recommended to represent cross-link restraints. It was also found that significant improvement of model quality upon the introduction of cross-link restraints is obtained when the sum of differences of indices of cross-linked residues exceeds 150.
Collapse
Affiliation(s)
- Mateusz Leśniewski
- Faculty
of Chemistry, University of Gdańsk, Fahrenheit Union of Universities, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Maciej Pyrka
- Faculty
of Chemistry, University of Gdańsk, Fahrenheit Union of Universities, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
- Department
of Physics and Biophysics, University of
Warmia and Mazury, ul. Oczapowskiego 4, 10-719 Olsztyn, Poland
| | - Cezary Czaplewski
- Faculty
of Chemistry, University of Gdańsk, Fahrenheit Union of Universities, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Nguyen Truong Co
- Faculty
of Chemistry, University of Gdańsk, Fahrenheit Union of Universities, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Yida Jiang
- College
of Chemistry and Molecular Engineering & Center for Quantitative
Biology & PKU-Tsinghua Center for Life Sciences & Beijing
National Laboratory for Molecular Sciences, Peking University, Beijing 100871, China
| | - Zhou Gong
- Innovation
Academy of Precision Measurement Science and Technology, Chinese Academy of Sciences, 30 W. Xiao Hong Shan, Wuhan 430071, China
| | - Chun Tang
- College
of Chemistry and Molecular Engineering & Center for Quantitative
Biology & PKU-Tsinghua Center for Life Sciences & Beijing
National Laboratory for Molecular Sciences, Peking University, Beijing 100871, China
| | - Adam Liwo
- Faculty
of Chemistry, University of Gdańsk, Fahrenheit Union of Universities, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
2
|
Borges-Araújo L, Patmanidis I, Singh AP, Santos LHS, Sieradzan AK, Vanni S, Czaplewski C, Pantano S, Shinoda W, Monticelli L, Liwo A, Marrink SJ, Souza PCT. Pragmatic Coarse-Graining of Proteins: Models and Applications. J Chem Theory Comput 2023; 19:7112-7135. [PMID: 37788237 DOI: 10.1021/acs.jctc.3c00733] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.
Collapse
Affiliation(s)
- Luís Borges-Araújo
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Ilias Patmanidis
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Akhil P Singh
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
| | - Lucianna H S Santos
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
- Institut de Pharmacologie Moléculaire et Cellulaire, Université Côte d'Azur, Inserm, CNRS, 06560 Valbonne, France
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Sergio Pantano
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita, Okayama 700-8530, Japan
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Adam Liwo
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| |
Collapse
|
3
|
Danielsson A, Samsonov SA, Liwo A, Sieradzan AK. Extension of the SUGRES-1P Coarse-Grained Model of Polysaccharides to Heparin. J Chem Theory Comput 2023; 19:6023-6036. [PMID: 37587433 PMCID: PMC10500997 DOI: 10.1021/acs.jctc.3c00511] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2023] [Indexed: 08/18/2023]
Abstract
Heparin is an unbranched periodic polysaccharide composed of negatively charged monomers and involved in key biological processes, including anticoagulation, angiogenesis, and inflammation. Its structure and dynamics have been studied extensively using experimental as well as theoretical approaches. The conventional approach of computational chemistry applied to the analysis of biomolecules is all-atom molecular dynamics, which captures the interactions of individual atoms by solving Newton's equation of motion. An alternative is molecular dynamics simulations using coarse-grained models of biomacromolecules, which offer a reduction of the representation and consequently enable us to extend the time and size scale of simulations by orders of magnitude. In this work, we extend the UNIfied COarse-gRaiNed (UNICORN) model of biological macromolecules developed in our laboratory to heparin. We carried out extensive tests to estimate the optimal weights of energy terms of the effective energy function as well as the optimal Debye-Hückel screening factor for electrostatic interactions. We applied the model to study unbound heparin molecules of polymerization degree ranging from 6 to 68 residues. We compare the obtained coarse-grained heparin conformations with models obtained from X-ray diffraction studies of heparin. The SUGRES-1P force field was able to accurately predict the general shape and global characteristics of heparin molecules.
Collapse
Affiliation(s)
- Annemarie Danielsson
- Faculty of Chemistry, University
of Gdansk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Sergey A. Samsonov
- Faculty of Chemistry, University
of Gdansk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam Liwo
- Faculty of Chemistry, University
of Gdansk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam K. Sieradzan
- Faculty of Chemistry, University
of Gdansk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
4
|
Liwo A, Pyrka M, Czaplewski C, Peng X, Niemi AJ. Long-Time Dynamics of Selected Molecular-Motor Components Using a Physics-Based Coarse-Grained Approach. Biomolecules 2023; 13:941. [PMID: 37371521 DOI: 10.3390/biom13060941] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2023] [Revised: 05/26/2023] [Accepted: 05/30/2023] [Indexed: 06/29/2023] Open
Abstract
Molecular motors are essential for the movement and transportation of macromolecules in living organisms. Among them, rotatory motors are particularly efficient. In this study, we investigated the long-term dynamics of the designed left-handed alpha/alpha toroid (PDB: 4YY2), the RBM2 flagellum protein ring from Salmonella (PDB: 6SD5), and the V-type Na+-ATPase rotor in Enterococcus hirae (PDB: 2BL2) using microcanonical and canonical molecular dynamics simulations with the coarse-grained UNRES force field, including a lipid-membrane model, on a millisecond laboratory time scale. Our results demonstrate that rotational motion can occur with zero total angular momentum in the microcanonical regime and that thermal motions can be converted into net rotation in the canonical regime, as previously observed in simulations of smaller cyclic molecules. For 6SD5 and 2BL2, net rotation (with a ratcheting pattern) occurring only about the pivot of the respective system was observed in canonical simulations. The extent and direction of the rotation depended on the initial conditions. This result suggests that rotatory molecular motors can convert thermal oscillations into net rotational motion. The energy from ATP hydrolysis is required probably to set the direction and extent of rotation. Our findings highlight the importance of molecular-motor structures in facilitating movement and transportation within living organisms.
Collapse
Affiliation(s)
- Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Maciej Pyrka
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities, Wita Stwosza 63, 80-308 Gdańsk, Poland
- Department of Physics and Biophysics, University of Warmia and Mazury, ul. Oczapowskiego 4, 10-719 Olsztyn, Poland
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Xubiao Peng
- Center for Quantum Technology Research, Key Laboratory of Advanced Optoelectronic Quantum Architecture and Measurements (MOE), School of Physics, Beijing Institute of Technology, Beijing 100081, China
| | - Antti J Niemi
- Nordita, Stockholm University and Uppsala University, Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden
| |
Collapse
|
5
|
Lipska AG, Sieradzan AK, Czaplewski C, Lipińska AD, Ocetkiewicz KM, Proficz J, Czarnul P, Krawczyk H, Liwo A. Long-time scale simulations of virus-like particles from three human-norovirus strains. J Comput Chem 2023; 44:1470-1483. [PMID: 36799410 DOI: 10.1002/jcc.27087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2022] [Revised: 12/22/2022] [Accepted: 01/29/2023] [Indexed: 02/18/2023]
Abstract
The dynamics of the virus like particles (VLPs) corresponding to the GII.4 Houston, GII.2 SMV, and GI.1 Norwalk strains of human noroviruses (HuNoV) that cause gastroenteritis was investigated by means of long-time (about 30 μs in the laboratory timescale) molecular dynamics simulations with the coarse-grained UNRES force field. The main motion of VLP units turned out to be the bending at the junction between the P1 subdomain (that sits in the VLP shell) and the P2 subdomain (that protrudes outside) of the major VP1 protein, this resulting in a correlated wagging motion of the P2 subdomains with respect to the VLP surface. The fluctuations of the P2 subdomain were found to be more pronounced and the P2 domain made a greater angle with the normal to the VLP surface for the GII.2 strain, which could explain the inability of this strain to bind the histo-blood group antigens (HBGAs).
Collapse
Affiliation(s)
- Agnieszka G Lipska
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Adam K Sieradzan
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Cezary Czaplewski
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Andrea D Lipińska
- Laboratory of Virus Molecular Biology, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Krzysztof M Ocetkiewicz
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Jerzy Proficz
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Paweł Czarnul
- Faculty of Electronics, Telecommunications and Informatics, Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Henryk Krawczyk
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Faculty of Electronics, Telecommunications and Informatics, Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Adam Liwo
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| |
Collapse
|
6
|
Sieradzan AK, Sans-Duñó J, Lubecka EA, Czaplewski C, Lipska AG, Leszczyński H, Ocetkiewicz KM, Proficz J, Czarnul P, Krawczyk H, Liwo A. Optimization of parallel implementation of UNRES package for coarse-grained simulations to treat large proteins. J Comput Chem 2023; 44:602-625. [PMID: 36378078 DOI: 10.1002/jcc.27026] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2022] [Revised: 09/19/2022] [Accepted: 10/17/2022] [Indexed: 11/16/2022]
Abstract
We report major algorithmic improvements of the UNRES package for physics-based coarse-grained simulations of proteins. These include (i) introduction of interaction lists to optimize computations, (ii) transforming the inertia matrix to a pentadiagonal form to reduce computing and memory requirements, (iii) removing explicit angles and dihedral angles from energy expressions and recoding the most time-consuming energy/force terms to minimize the number of operations and to improve numerical stability, (iv) using OpenMP to parallelize those sections of the code for which distributed-memory parallelization involves unfavorable computing/communication time ratio, and (v) careful memory management to minimize simultaneous access of distant memory sections. The new code enables us to run molecular dynamics simulations of protein systems with size exceeding 100,000 amino-acid residues, reaching over 1 ns/day (1 μs/day in all-atom timescale) with 24 cores for proteins of this size. Parallel performance of the code and comparison of its performance with that of AMBER, GROMACS and MARTINI 3 is presented.
Collapse
Affiliation(s)
- Adam K Sieradzan
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Jordi Sans-Duñó
- Department of Chemistry, University of Lleida, Lleida, Spain
| | - Emilia A Lubecka
- Faculty of Electronics, Telecommunication and Informatics, Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Agnieszka G Lipska
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Henryk Leszczyński
- Faculty of Mathematics, Physics and Informatics, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Krzysztof M Ocetkiewicz
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Jerzy Proficz
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Paweł Czarnul
- Faculty of Electronics, Telecommunication and Informatics, Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Henryk Krawczyk
- Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Faculty of Electronics, Telecommunication and Informatics, Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland.,Centre of Informatics Tri-city Academic Supercomputer and Network (CI TASK), Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| |
Collapse
|
7
|
Ślusarz R, Lubecka EA, Czaplewski C, Liwo A. Improvements and new functionalities of UNRES server for coarse-grained modeling of protein structure, dynamics, and interactions. Front Mol Biosci 2022; 9:1071428. [PMID: 36589235 PMCID: PMC9794589 DOI: 10.3389/fmolb.2022.1071428] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2022] [Accepted: 11/29/2022] [Indexed: 12/15/2022] Open
Abstract
In this paper we report the improvements and extensions of the UNRES server (https://unres-server.chem.ug.edu.pl) for physics-based simulations with the coarse-grained UNRES model of polypeptide chains. The improvements include the replacement of the old code with the recently optimized one and adding the recent scale-consistent variant of the UNRES force field, which performs better in the modeling of proteins with the β and the α+β structures. The scope of applications of the package was extended to data-assisted simulations with restraints from nuclear magnetic resonance (NMR) and chemical crosslink mass-spectroscopy (XL-MS) measurements. NMR restraints can be input in the NMR Exchange Format (NEF), which has become a standard. Ambiguous NMR restraints are handled without expert intervention owing to a specially designed penalty function. The server can be used to run smaller jobs directly or to prepare input data to run larger production jobs by using standalone installations of UNRES.
Collapse
Affiliation(s)
- Rafał Ślusarz
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Emilia A. Lubecka
- Faculty of Electronics, Telecommunication and Informatics, Gdańsk University of Technology, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Fahrenheit Union of Universities in Gdańsk, Gdańsk, Poland,*Correspondence: Adam Liwo,
| |
Collapse
|
8
|
Prediction of Aggregation of Biologically-Active Peptides with the UNRES Coarse-Grained Model. Biomolecules 2022; 12:biom12081140. [PMID: 36009034 PMCID: PMC9406146 DOI: 10.3390/biom12081140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2022] [Revised: 08/13/2022] [Accepted: 08/16/2022] [Indexed: 11/17/2022] Open
Abstract
The UNited RESidue (UNRES) model of polypeptide chains was applied to study the association of 20 peptides with sizes ranging from 6 to 32 amino-acid residues. Twelve of those were potentially aggregating hexa- or heptapeptides excised from larger proteins, while the remaining eight contained potentially aggregating sequences, functionalized by attaching larger ends rich in charged residues. For 13 peptides, the experimental data of aggregation were used. The remaining seven were synthesized, and their properties were measured in this work. Multiplexed replica-exchange simulations of eight-chain systems were conducted at 12 temperatures from 260 to 370 K at concentrations from 0.421 to 5.78 mM, corresponding to the experimental conditions. The temperature profiles of the fractions of monomers and octamers showed a clear transition corresponding to aggregate dissociation. Low simulated transition temperatures were obtained for the peptides, which did not precipitate after incubation, as well as for the H-GNNQQNY-NH2 prion–protein fragment, which forms small fibrils. A substantial amount of inter-strand β-sheets was found in most of the systems. The results suggest that UNRES simulations can be used to assess peptide aggregation except for glutamine- and asparagine-rich peptides, for which a revision of the UNRES sidechain–sidechain interaction potentials appears necessary.
Collapse
|
9
|
Smardz P, Sieradzan AK, Krupa P. Mechanical Stability of Ribonuclease A Heavily Depends on the Redox Environment. J Phys Chem B 2022; 126:6240-6249. [PMID: 35975925 PMCID: PMC9421896 DOI: 10.1021/acs.jpcb.2c04718] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Disulfide bonds are covalent bonds that connect nonlocal fragments of proteins, and they are unique post-translational modifications of proteins. They require the oxidizing environment to be stable, which occurs for example during oxidative stress; however, in a cell the reductive environment is maintained, lowering their stability. Despite many years of research on disulfide bonds, their role in the protein life cycle is not fully understood and seems to strictly depend on a system or process in which they are involved. In this article, coarse-grained UNited RESidue (UNRES), and all-atom Assisted Model Building with Energy Refinement (AMBER) force fields were applied to run a series of steered molecular dynamics (SMD) simulations of one of the most studied, but still not fully understood, proteins─ribonuclease A (RNase A). SMD simulations were performed to study the mechanical stability of RNase A in different oxidative-reductive environments. As disulfide bonds (and any other covalent bonds) cannot break/form in any classical all-atom force field, we applied additional restraints between sulfur atoms of reduced cysteines which were able to mimic the breaking of the disulfide bonds. On the other hand, the coarse-grained UNRES force field enables us to study the breaking/formation of the disulfide bonds and control the reducing/oxidizing environment owing to the presence of the designed distance/orientation-dependent potential. This study reveals that disulfide bonds have a strong influence on the mechanical stability of RNase A only in a highly oxidative environment. However, the local stability of the secondary structure seems to play a major factor in the overall stability of the protein. Both our thermal unfolding and mechanical stretching studies show that the most stable disulfide bond is Cys65-Cys72. The breaking of disulfide bonds Cys26-Cys84 and Cys58-Cys110 is associated with large force peaks. They are structural bridges, which are mostly responsible for stabilizing the RNase A conformation, while the presence of the remaining two bonds (Cys65-Cys72 and Cys40-Cys95) is most likely connected with the enzymatic activity rather than the structural stability of RNase A in the cytoplasm. Our results prove that disulfide bonds are indeed stabilizing fragments of the proteins, but their role is strongly redox environment-dependent.
Collapse
Affiliation(s)
- Pamela Smardz
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Paweł Krupa
- Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| |
Collapse
|
10
|
Rojas AV, Maisuradze GG, Scheraga HA, Liwo A. Probing Protein Aggregation Using the Coarse-Grained UNRES Force Field. Methods Mol Biol 2022; 2340:79-104. [PMID: 35167071 DOI: 10.1007/978-1-0716-1546-1_5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Protein aggregation is the cause of many, often lethal, diseases, including the Alzheimer's, Parkinson's, and Huntington's diseases, and familial amyloidosis. Theoretical investigation of the mechanism of this process, including the structures of the oligomeric intermediates which are the most toxic, is difficult because of long time scale of aggregation. Coarse-grained models, which enable us to extend the simulation time scale by three or more orders of magnitude, are, therefore, of great advantage in such studies. In this chapter, we describe the application of the physics-based UNited RESidue (UNRES) force field developed in our laboratory to study protein aggregation, in both free simulations and simulations of aggregation propagation from an existing template (seed), and illustrate it with the examples of Aβ-peptide aggregation and Aβ-peptide-assisted aggregation of the peptides derived from the repeat domains of tau (TauRD).
Collapse
Affiliation(s)
- Ana V Rojas
- Schrodinger Inc., 120 West 45th Street New York, New York, 10036, NY, USA
| | - Gia G Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, 14853-1301, NY, USA
| | - Harold A Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, 14853-1301, NY, USA
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk, 80-308, Poland.
| |
Collapse
|
11
|
Abstract
The biological significance of proteins attracted the scientific community in exploring their characteristics. The studies shed light on the interaction patterns and functions of proteins in a living body. Due to their practical difficulties, reliable experimental techniques pave the way for introducing computational methods in the interaction prediction. Automated methods reduced the difficulties but could not yet replace experimental studies as the field is still evolving. Interaction prediction problem being critical needs highly accurate results, but none of the existing methods could offer reliable performance that can parallel with experimental results yet. This article aims to assess the existing computational docking algorithms, their challenges, and future scope. Blind docking techniques are quite helpful when no information other than the individual structures are available. As more and more complex structures are being added to different databases, information-driven approaches can be a good alternative. Artificial intelligence, ruling over the major fields, is expected to take over this domain very shortly.
Collapse
|
12
|
Liwo A, Czaplewski C, Sieradzan AK, Lipska AG, Samsonov SA, Murarka RK. Theory and Practice of Coarse-Grained Molecular Dynamics of Biologically Important Systems. Biomolecules 2021; 11:1347. [PMID: 34572559 PMCID: PMC8465211 DOI: 10.3390/biom11091347] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Revised: 09/03/2021] [Accepted: 09/09/2021] [Indexed: 12/16/2022] Open
Abstract
Molecular dynamics with coarse-grained models is nowadays extensively used to simulate biomolecular systems at large time and size scales, compared to those accessible to all-atom molecular dynamics. In this review article, we describe the physical basis of coarse-grained molecular dynamics, the coarse-grained force fields, the equations of motion and the respective numerical integration algorithms, and selected practical applications of coarse-grained molecular dynamics. We demonstrate that the motion of coarse-grained sites is governed by the potential of mean force and the friction and stochastic forces, resulting from integrating out the secondary degrees of freedom. Consequently, Langevin dynamics is a natural means of describing the motion of a system at the coarse-grained level and the potential of mean force is the physical basis of the coarse-grained force fields. Moreover, the choice of coarse-grained variables and the fact that coarse-grained sites often do not have spherical symmetry implies a non-diagonal inertia tensor. We describe selected coarse-grained models used in molecular dynamics simulations, including the most popular MARTINI model developed by Marrink's group and the UNICORN model of biological macromolecules developed in our laboratory. We conclude by discussing examples of the application of coarse-grained molecular dynamics to study biologically important processes.
Collapse
Affiliation(s)
- Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Adam K. Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Agnieszka G. Lipska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Sergey A. Samsonov
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Rajesh K. Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India;
| |
Collapse
|
13
|
Antoniak A, Biskupek I, Bojarski KK, Czaplewski C, Giełdoń A, Kogut M, Kogut MM, Krupa P, Lipska AG, Liwo A, Lubecka EA, Marcisz M, Maszota-Zieleniak M, Samsonov SA, Sieradzan AK, Ślusarz MJ, Ślusarz R, Wesołowski PA, Ziȩba K. Modeling protein structures with the coarse-grained UNRES force field in the CASP14 experiment. J Mol Graph Model 2021; 108:108008. [PMID: 34419932 DOI: 10.1016/j.jmgm.2021.108008] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2021] [Revised: 08/12/2021] [Accepted: 08/13/2021] [Indexed: 12/31/2022]
Abstract
The UNited RESidue (UNRES) force field was tested in the 14th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP14), in which larger oligomeric and multimeric targets were present compared to previous editions. Three prediction modes were tested (i) ab initio (the UNRES group), (ii) contact-assisted (the UNRES-contact group), and (iii) template-assisted (the UNRES-template group). For most of the targets, the contact restraints were derived from the server models top-ranked by the DeepQA method, while the DNCON2 method was used for 11 targets. Our consensus-fragment procedure was used to run template-assisted predictions. Each group also processed the Nuclear Magnetic Resonance (NMR)- and Small Angle X-Ray Scattering (SAXS)-data assisted targets. The average Global Distance Test Total Score (GDT_TS) of the 'Model 1' predictions were 29.17, 39.32, and 56.37 for the UNRES, UNRES-contact, and UNRES-template predictions, respectively, increasing by 0.53, 2.24, and 3.76, respectively, compared to CASP13. It was also found that the GDT_TS of the UNRES models obtained in ab initio mode and in the contact-assisted mode decreases with the square root of chain length, while the exponent in this relationship is 0.20 for the UNRES-template group models and 0.11 for the best performing AlphaFold2 models, which suggests that incorporation of database information, which stems from protein evolution, brings in long-range correlations, thus enabling the correction of force-field inaccuracies.
Collapse
Affiliation(s)
- Anna Antoniak
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Iga Biskupek
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Krzysztof K Bojarski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Artur Giełdoń
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Mateusz Kogut
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Małgorzata M Kogut
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Paweł Krupa
- Institute of Physics, Polish Academy of Sciences, Aleja Lotników 32/46, Warsaw, PL-02668, Poland
| | - Agnieszka G Lipska
- 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; School of Computational Sciences, Korea Institute for Advanced Study, 87 Hoegiro, Dongdaemun-gu, 130-722, Seoul, Republic of Korea.
| | - Emilia A Lubecka
- Faculty of Electronics, Telecommunications and Informatics, Gdańsk University of Technology, G. Narutowicza 11/12, 80-233, Gdańsk, Poland
| | - Mateusz Marcisz
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland; Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, ul. Abrahama 58, 80-307, Gdańsk, Poland
| | | | - Sergey A Samsonov
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Magdalena J Ślusarz
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Rafał Ślusarz
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Patryk A Wesołowski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland; Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, ul. Abrahama 58, 80-307, Gdańsk, Poland
| | - Karolina Ziȩba
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| |
Collapse
|
14
|
Kogut M, Gong Z, Tang C, Liwo A. Pseudopotentials for coarse-grained cross-link-assisted modeling of protein structures. J Comput Chem 2021; 42:2054-2067. [PMID: 34402552 DOI: 10.1002/jcc.26736] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2021] [Revised: 07/09/2021] [Accepted: 08/03/2021] [Indexed: 11/08/2022]
Abstract
Pseudopotentials for the chemical cross-links comprising the glutamic- and aspartic-acid side chains bridged with adipic- (ADH) or pimelic-acid hydrazide (PDH), and the lysine side chains bridged with glutaric (BS2 G) or suberic acid (BS3 ) for coarse-grained cross-link-assisted simulations were determined by canonical molecular dynamics with the Amber14sb force field. The potentials depend on the distance between side-chain ends and on side-chain orientation, this preventing from making cross-link contacts across the globule in simulations. The potentials were implemented in the UNRES coarse-grained force field and their effect on the quality of models was assessed with 11 monomeric and 1 dimeric proteins, using synthetic or experimental cross-link data. Simulations with the new potentials resulted in improvement of the generated models compared to unrestrained simulations in more instances compared to those with the statistical potentials.
Collapse
Affiliation(s)
- Mateusz Kogut
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland
| | - Zhou Gong
- Innovation Academy of Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan, China
| | - Chun Tang
- College of Chemistry and Molecular Engineering, PKU-Tsinghua Center for Life Sciences, Beijing National Laboratory for Molecular Sciences, Peking University, Beijing, China
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland
| |
Collapse
|
15
|
Lubecka EA, Liwo A. ESCASA: Analytical estimation of atomic coordinates from coarse-grained geometry for nuclear-magnetic-resonance-assisted protein structure modeling. I. Backbone and H β protons. J Comput Chem 2021; 42:1579-1589. [PMID: 34048074 DOI: 10.1002/jcc.26695] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2021] [Revised: 05/06/2021] [Accepted: 05/11/2021] [Indexed: 12/13/2022]
Abstract
A method for the estimation of coordinates of atoms in proteins from coarse-grained geometry by simple analytical formulas (ESCASA), for use in nuclear-magnetic-resonance (NMR) data-assisted coarse-grained simulations of proteins is proposed. In this paper, the formulas for the backbone Hα and amide (HN ) protons, and the side-chain Hβ protons, given the Cα -trace, have been derived and parameterized, by using the interproton distances calculated from a set of 140 high-resolution non-homologous protein structures. The mean standard deviation over all types of proton pairs in the set was 0.44 Å after fitting. Validation against a set of 41 proteins with NMR-determined structures, which were not considered in parameterization, resulted in average standard deviation from average proton-proton distances of the NMR-determined structures of 0.25 Å, compared to 0.21 Å obtained with the PULCHRA all-atom-chain reconstruction algorithm and to the 0.12 Å standard deviation of the average-structure proton-proton distance of NMR-determined ensembles. The formulas provide analytical forces and can, therefore, be used in coarse-grained molecular dynamics.
Collapse
Affiliation(s)
- Emilia A Lubecka
- Faculty of Electronics, Telecommunications and Informatics, Gdańsk University of Technology, Gdańsk, Poland
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland
| |
Collapse
|
16
|
Sieradzan AK, Korneev A, Begun A, Kachlishvili K, Scheraga HA, Molochkov A, Senet P, Niemi AJ, Maisuradze GG. Investigation of Phosphorylation-Induced Folding of an Intrinsically Disordered Protein by Coarse-Grained Molecular Dynamics. J Chem Theory Comput 2021; 17:3203-3220. [PMID: 33909430 DOI: 10.1021/acs.jctc.1c00155] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
Apart from being the most common mechanism of regulating protein function and transmitting signals throughout the cell, phosphorylation has an ability to induce disorder-to-order transition in an intrinsically disordered protein. In particular, it was shown that folding of the intrinsically disordered protein, eIF4E-binding protein isoform 2 (4E-BP2), can be induced by multisite phosphorylation. Here, the principles that govern the folding of phosphorylated 4E-BP2 (pT37pT46 4E-BP218-62) are investigated by analyzing canonical and replica exchange molecular dynamics trajectories, generated with the coarse-grained united-residue force field, in terms of local and global motions and the time dependence of formation of contacts between Cαs of selected pairs of residues. The key residues involved in the folding of the pT37pT46 4E-BP218-62 are elucidated by this analysis. The correlations between local and global motions are identified. Moreover, for a better understanding of the physics of the formation of the folded state, the experimental structure of the pT37pT46 4E-BP218-62 is analyzed in terms of a kink (heteroclinic standing wave solution) of a generalized discrete nonlinear Schrödinger equation. It is shown that without molecular dynamics simulations the kinks are able to identify not only the phosphorylated sites of protein, the key players in folding, but also the reasons for the weak stability of the pT37pT46 4E-BP218-62.
Collapse
Affiliation(s)
- Adam K Sieradzan
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Anatolii Korneev
- Pacific Quantum Center, Far Eastern Federal University, 10 Ajax Bay, 690922 Russky Island, Vladivostok, Russia
| | - Alexander Begun
- Pacific Quantum Center, Far Eastern Federal University, 10 Ajax Bay, 690922 Russky Island, Vladivostok, Russia
| | - Khatuna Kachlishvili
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
| | - Harold A Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
| | - Alexander Molochkov
- Pacific Quantum Center, Far Eastern Federal University, 10 Ajax Bay, 690922 Russky Island, Vladivostok, Russia
| | - Patrick Senet
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States.,Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 6303 CNRS-Université de Bourgogne Franche-Comté, 9 Av. A. Savary, BP 47 870, F-21078 Dijon Cedex, France
| | - Antti J Niemi
- Pacific Quantum Center, Far Eastern Federal University, 10 Ajax Bay, 690922 Russky Island, Vladivostok, Russia.,Laboratoire de Mathematiques et Physique Theorique, CNRS UMR 6083, Fédération Denis Poisson, Université de Tours, Parc de Grandmont, F37200 Tours, France.,Nordita, Stockholm University and Uppsala University, Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden.,School of Physics, Beijing Institute of Technology, Beijing 100081, P.R. China
| | - Gia G Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
| |
Collapse
|
17
|
Wang Y, Song E, Pei Y, Li C, Liu J. Theoretical Study of Fast Calculation of Damping Loss Factors for Rubber Polymers. J Phys Chem Lett 2020; 11:6025-6031. [PMID: 32639736 DOI: 10.1021/acs.jpclett.0c01661] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Rubber polymers as acoustic damping materials are widely used in industrial sectors. However, there is still no theoretical method to rapidly predict the damping properties of rubber polymers. Here we propose a theoretical method for the fast calculation of damping loss factors in rubber polymers. Molecular dynamics simulations were employed to construct the optimal configuration of internal polymer chains with density as the descriptor in order to calculate the mechanical properties and thus the damping loss factors. The acceptable difference between the calculated and experimental damping loss factors shows that the proposed method is efficient to predict the intrinsic damping properties of rubber polymers. Through studying the relationship among composition, microstructure, and the damping property, we found that the carbon-carbon skeletal chains with lateral methyl groups produce high damping properties.
Collapse
Affiliation(s)
- Youwei Wang
- State Key Laboratory of High Performance Ceramics and Superfine Microstructures, Shanghai Institute of Ceramics, Chinese Academy of Sciences, Shanghai 200050, China
- College of Materials Science and Optoelectronic Technology, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Erhong Song
- State Key Laboratory of High Performance Ceramics and Superfine Microstructures, Shanghai Institute of Ceramics, Chinese Academy of Sciences, Shanghai 200050, China
- College of Materials Science and Optoelectronic Technology, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yongwei Pei
- School of Materials and Physics, China University of Mining and Technology, Daxue Road, No. 1, Xuzhou, Jiangsu 221116, China
| | - Chong Li
- School of Materials and Physics, China University of Mining and Technology, Daxue Road, No. 1, Xuzhou, Jiangsu 221116, China
| | - Jianjun Liu
- State Key Laboratory of High Performance Ceramics and Superfine Microstructures, Shanghai Institute of Ceramics, Chinese Academy of Sciences, Shanghai 200050, China
- College of Materials Science and Optoelectronic Technology, University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
18
|
Kachlishvili K, Korneev A, Maisuradze L, Liu J, Scheraga HA, Molochkov A, Senet P, Niemi AJ, Maisuradze GG. New Insights into Folding, Misfolding, and Nonfolding Dynamics of a WW Domain. J Phys Chem B 2020; 124:3855-3872. [PMID: 32271570 DOI: 10.1021/acs.jpcb.0c00628] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Intermediate states in protein folding are associated with formation of amyloid fibrils, which are responsible for a number of neurodegenerative diseases. Therefore, prevention of the aggregation of folding intermediates is one of the most important problems to overcome. Recently, we studied the origins and prevention of formation of intermediate states with the example of the Formin binding protein 28 (FBP28) WW domain. We demonstrated that the replacement of Leu26 by Asp26 or Trp26 (in ∼15% of the folding trajectories) can alter the folding scenario from three-state folding, a major folding scenario for the FBP28 WW domain (WT) and its mutants, toward two-state or downhill folding at temperatures below the melting point. Here, for a better understanding of the physics of the formation/elimination of intermediates, (i) the dynamics and energetics of formation of β-strands in folding, misfolding, and nonfolding trajectories of these mutants (L26D and L26W) is investigated; (ii) the experimental structures of WT, L26D, and L26W are analyzed in terms of a kink (heteroclinic standing wave solution) of a generalized discrete nonlinear Schrödinger equation. We show that the formation of each β-strand in folding trajectories is accompanied by the emergence of kinks in internal coordinate space as well as a decrease in local free energy. In particular, the decrease in downhill folding trajectory is ∼7 kcal/mol, while it varies between 31 and 48 kcal/mol for the three-state folding trajectory. The kink analyses of the experimental structures give new insights into formation of intermediates, which may become a useful tool for preventing aggregation.
Collapse
Affiliation(s)
- Khatuna Kachlishvili
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca 14853-1301, New York, United States
| | - Anatolii Korneev
- Laboratory of Physics of Living Matter, Far Eastern Federal University, Sukhanova 8, Vladivostok 690950, Russia
| | - Luka Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca 14853-1301, New York, United States.,Biochemistry, Quantitative Biology, Biophysics, and Structural Biology (BQBS) Track, Yale University, New Haven 06520-8084, ConnecticutUnited States
| | - Jiaojiao Liu
- School of Physics, Beijing Institute of Technology, Beijing 100081, P. R. China
| | - Harold A Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca 14853-1301, New York, United States
| | - Alexander Molochkov
- Laboratory of Physics of Living Matter, Far Eastern Federal University, Sukhanova 8, Vladivostok 690950, Russia
| | - Patrick Senet
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca 14853-1301, New York, United States.,Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR 6303 CNRS-Univ. de Bourgogne Franche-Comté, 9 Av. A. Savary, BP 47 870, Dijon Cedex F-21078, France
| | - Antti J Niemi
- Laboratory of Physics of Living Matter, Far Eastern Federal University, Sukhanova 8, Vladivostok 690950, Russia.,School of Physics, Beijing Institute of Technology, Beijing 100081, P. R. China.,Laboratoire de Mathematiques et Physique Theorique CNRS UMR 6083, Fédération Denis Poisson, Université de Tours, Parc de Grandmont, Tours F37200, France.,Nordita, Stockholm University, Roslagstullsbacken 23, Stockholm SE-106 91, Sweden
| | - Gia G Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca 14853-1301, New York, United States
| |
Collapse
|
19
|
Liwo A, Czaplewski C, Sieradzan AK, Lubecka EA, Lipska AG, Golon Ł, Karczyńska A, Krupa P, Mozolewska MA, Makowski M, Ganzynkowicz R, Giełdoń A, Maciejczyk M. Scale-consistent approach to the derivation of coarse-grained force fields for simulating structure, dynamics, and thermodynamics of biopolymers. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:73-122. [PMID: 32145953 DOI: 10.1016/bs.pmbts.2019.12.004] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
In this chapter the scale-consistent approach to the derivation of coarse-grained force fields developed in our laboratory is presented, in which the effective energy function originates from the potential of mean force of the system under consideration and embeds atomistically detailed interactions in the resulting energy terms through use of Kubo's cluster-cumulant expansion, appropriate selection of the major degrees of freedom to be averaged out in the derivation of analytical approximations to the energy terms, and appropriate expression of the interaction energies at the all-atom level in these degrees of freedom. Our approach enables the developers to find correct functional forms of the effective coarse-grained energy terms, without having to import them from all-atom force fields or deriving them on a heuristic basis. In particular, the energy terms derived in such a way exhibit correct dependence on coarse-grained geometry, in particular on site orientation. Moreover, analytical formulas for the multibody (correlation) terms, which appear to be crucial for coarse-grained modeling of many of the regular structures such as, e.g., protein α-helices and β-sheets, can be derived in a systematic way. Implementation of the developed theory to the UNIfied COarse-gRaiNed (UNICORN) model of biological macromolecules, which consists of the UNRES (for proteins), NARES-2P (for nucleic acids), and SUGRES-1P (for polysaccharides) components, and is being developed in our laboratory is described. Successful applications of UNICORN to the prediction of protein structure, simulating the folding and stability of proteins and nucleic acids, and solving biological problems are discussed.
Collapse
Affiliation(s)
- Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland; School of Computational Sciences, Korea Institute for Advanced Study, Seoul, Republic of Korea.
| | | | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland; School of Computational Sciences, Korea Institute for Advanced Study, Seoul, Republic of Korea
| | - Emilia A Lubecka
- Institute of Informatics, Faculty of Mathematics, Physics, and Informatics, University of Gdańsk, Gdańsk, Poland
| | | | - Łukasz Golon
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland
| | | | - Paweł Krupa
- Institute of Physics, Polish Academy of Sciences, Warsaw, Poland
| | | | | | | | - Artur Giełdoń
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland
| | - Maciej Maciejczyk
- Department of Physics and Biophysics, Faculty of Food Science, University of Warmia and Mazury in Olsztyn, Olsztyn, Poland
| |
Collapse
|
20
|
Karczyńska AS, Ziȩba K, Uciechowska U, Mozolewska MA, Krupa P, Lubecka EA, Lipska AG, Sikorska C, Samsonov SA, Sieradzan AK, Giełdoń A, Liwo A, Ślusarz R, Ślusarz M, Lee J, Joo K, Czaplewski C. Improved Consensus-Fragment Selection in Template-Assisted Prediction of Protein Structures with the UNRES Force Field in CASP13. J Chem Inf Model 2020; 60:1844-1864. [PMID: 31999919 PMCID: PMC7588044 DOI: 10.1021/acs.jcim.9b00864] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The method for protein-structure
prediction, which combines the
physics-based coarse-grained UNRES force field with knowledge-based
modeling, has been developed further and tested in the 13th Community
Wide Experiment on the Critical Assessment of Techniques for Protein
Structure Prediction (CASP13). The method implements restraints from
the consensus fragments common to server models. In this work, the
server models to derive fragments have been chosen on the basis of
quality assessment; a fully automatic fragment-selection procedure
has been introduced, and Dynamic Fragment Assembly pseudopotentials
have been fully implemented. The Global Distance Test Score (GDT_TS),
averaged over our “Model 1” predictions, increased by
over 10 units with respect to CASP12 for the free-modeling category
to reach 40.82. Our “Model 1” predictions ranked 20
and 14 for all and free-modeling targets, respectively (upper 20.2%
and 14.3% of all models submitted to CASP13 in these categories, respectively),
compared to 27 (upper 21.1%) and 24 (upper 18.9%) in CASP12, respectively.
For oligomeric targets, the Interface Patch Similarity (IPS) and Interface
Contact Similarity (ICS) averaged over our best oligomer models increased
from 0.28 to 0.36 and from 12.4 to 17.8, respectively, from CASP12
to CASP13, and top-ranking models of 2 targets (H0968 and T0997o)
were obtained (none in CASP12). The improvement of our method in CASP13
over CASP12 was ascribed to the combined effect of the overall enhancement
of server-model quality, our success in selecting server models and
fragments to derive restraints, and improvements of the restraint
and potential-energy functions.
Collapse
Affiliation(s)
| | - Karolina Ziȩba
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Urszula Uciechowska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Magdalena A Mozolewska
- Institute of Computer Science, Polish Academy of Sciences, ul. Jana Kazimierza 5, Warsaw PL-02668, Poland
| | - Paweł Krupa
- Institute of Physics, Polish Academy of Sciences, Aleja Lotników 32/46, Warsaw PL-02668, Poland
| | - Emilia A Lubecka
- Institute of Informatics, Faculty of Mathematics, Physics, and Informatics, University of Gdańsk, Wita Stwosza 57, Gdańsk 80-308, Poland
| | - Agnieszka G Lipska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Celina Sikorska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Sergey A Samsonov
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland.,School of Computational Sciences, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu, Seoul 130-722, Republic of Korea
| | - Artur Giełdoń
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland.,School of Computational Sciences, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu, Seoul 130-722, Republic of Korea
| | - Rafał Ślusarz
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Magdalena Ślusarz
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Jooyoung Lee
- 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, Gdańsk 80-308, Poland
| |
Collapse
|
21
|
Liwo A, Czaplewski C. Extension of the force-matching method to coarse-grained models with axially symmetric sites to produce transferable force fields: Application to the UNRES model of proteins. J Chem Phys 2020; 152:054902. [DOI: 10.1063/1.5138991] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- 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, 87 Hoegiro, Dongdaemun-gu, 130-722 Seoul, South Korea
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
22
|
Ziȩba K, Ślusarz M, Ślusarz R, Liwo A, Czaplewski C, Sieradzan AK. Extension of the UNRES Coarse-Grained Force Field to Membrane Proteins in the Lipid Bilayer. J Phys Chem B 2019; 123:7829-7839. [PMID: 31454484 DOI: 10.1021/acs.jpcb.9b06700] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The physics-based UNRES coarse-grained force field for the simulations of protein structure and dynamics has been extended to treat membrane proteins. The lipid bilayer has been modeled by introducing a continuous nonpolar phase with the water-interface region of appropriate thickness. The potentials for average electrostatic and correlation interactions of the peptide groups have been rescaled to account for the reduction of the dielectric permittivity compared to the water phase and new potentials for protein side-chain-side-chain interactions inside and across the lipid phase have been introduced. The model was implemented in the UNRES package for coarse-grained simulations of proteins, and the package with the new functionality was tested for total energy conservation and thermostat behavior in microcanonical and canonical molecular dynamics simulations runs, respectively. The method was validated by running unrestricted ab initio blind-prediction tests of 10 short α-helical membrane proteins, all runs started from the extended structures. The modified UNRES force field was able to predict correctly the overall folds of the membrane proteins studied.
Collapse
Affiliation(s)
- Karolina Ziȩba
- Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Magdalena Ślusarz
- Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Rafał Ślusarz
- 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
| | - Cezary Czaplewski
- Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63 , 80-308 Gdańsk , Poland
| | - Adam K Sieradzan
- Faculty of Chemistry , University of Gdańsk , Wita Stwosza 63 , 80-308 Gdańsk , Poland
| |
Collapse
|
23
|
Lubecka EA, Karczyńska AS, Lipska AG, Sieradzan AK, Ziȩba K, Sikorska C, Uciechowska U, Samsonov SA, Krupa P, Mozolewska MA, Golon Ł, Giełdoń A, Czaplewski C, Ślusarz R, Ślusarz M, Crivelli SN, Liwo A. Evaluation of the scale-consistent UNRES force field in template-free prediction of protein structures in the CASP13 experiment. J Mol Graph Model 2019; 92:154-166. [PMID: 31376733 DOI: 10.1016/j.jmgm.2019.07.013] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2019] [Revised: 07/20/2019] [Accepted: 07/23/2019] [Indexed: 01/17/2023]
Abstract
The recent NEWCT-9P version of the coarse-grained UNRES force field for proteins, with scale-consistent formulas for the local and correlation terms, has been tested in the CASP13 experiment of the blind-prediction of protein structure, in the ab initio, contact-assisted, and data-assisted modes. Significant improvement of the performance has been observed with respect to the CASP11 and CASP12 experiments (by over 10 GDT_TS units for the ab initio mode predictions and by over 15 GDT_TS units for the contact-assisted prediction, respectively), which is a result of introducing scale-consistent terms and improved handling of contact-distance restraints. As in previous CASP exercises, UNRES ranked higher in the free modeling category than in the general category that included template based modeling targets. Use of distance restraints from the predicted contacts, albeit many of them were wrong, resulted in the increase of GDT_TS by over 8 units on average and introducing sparse restraints from small-angle X-ray/neutron scattering and chemical cross-link-mass-spectrometry experiments, and ambiguous restraints from nuclear magnetic resonance experiments has also improved the predictions by 8.6, 9.7, and 10.7 GDT_TS units on average, respectively.
Collapse
Affiliation(s)
- Emilia A Lubecka
- Institute of Informatics, Faculty of Mathematics, Physics, and 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
| | | | - Agnieszka G Lipska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Karolina Ziȩba
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Celina Sikorska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Urszula Uciechowska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Sergey A Samsonov
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Paweł Krupa
- Institute of Physics, Polish Academy of Sciences, Aleja Lotników 32/46, Warsaw, PL, 02668, Poland
| | - Magdalena A Mozolewska
- Institute of Computer Science, Polish Academy of Sciences, ul. Jana Kazimierza 5, Warsaw, 01-248, Poland
| | - Łukasz Golon
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Artur Giełdoń
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Rafał Ślusarz
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Magdalena Ślusarz
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Silvia N Crivelli
- Department of Computer Science, UC Davis, One Shields Ave., Davis, CA, 95616, USA
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland; School of Computational Sciences, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu, 130-722, Seoul, Republic of Korea.
| |
Collapse
|
24
|
Lubecka EA, Liwo A. Introduction of a bounded penalty function in contact-assisted simulations of protein structures to omit false restraints. J Comput Chem 2019; 40:2164-2178. [PMID: 31037754 DOI: 10.1002/jcc.25847] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2019] [Revised: 03/29/2019] [Accepted: 04/14/2019] [Indexed: 12/26/2022]
Abstract
Contact-assisted simulations, the contacts being predicted or determined experimentally, have become very important in the determination of the structures of proteins and other biological macromolecules. In this work, the effect of contact-distance restraints on the simulated structures was investigated with the use of multiplexed replica exchange simulations with the coarse-grained UNRES force field. A modified bounded flat-bottom restraint function that does not generate a gradient when a restraint cannot be satisfied was implemented. Calculations were run with (i) a set of four small proteins, with contact restraints derived from experimental structures, and (ii) selected CASP11 and CASP12 targets, with restraints as used at prediction time. The bounded penalty function largely omitted false contacts, which were usually inconsistent. It was found that at least 20% of correct contacts must be present in the restraint set to improve model quality with respect to unrestrained simulations. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Emilia A Lubecka
- Institute of Informatics, Faculty of Mathematics, Physics and 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
|
25
|
Liwo A, Sieradzan AK, Lipska AG, Czaplewski C, Joung I, Żmudzińska W, Hałabis A, Ołdziej S. A general method for the derivation of the functional forms of the effective energy terms in coarse-grained energy functions of polymers. III. Determination of scale-consistent backbone-local and correlation potentials in the UNRES force field and force-field calibration and validation. J Chem Phys 2019; 150:155104. [PMID: 31005069 DOI: 10.1063/1.5093015] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
The general theory of the construction of scale-consistent energy terms in the coarse-grained force fields presented in Paper I of this series has been applied to the revision of the UNRES force field for physics-based simulations of proteins. The potentials of mean force corresponding to backbone-local and backbone-correlation energy terms were calculated from the ab initio energy surfaces of terminally blocked glycine, alanine, and proline, and the respective analytical expressions, derived by using the scale-consistent formalism, were fitted to them. The parameters of all these potentials depend on single-residue types, thus reducing their number and preventing over-fitting. The UNRES force field with the revised backbone-local and backbone-correlation terms was calibrated with a set of four small proteins with basic folds: tryptophan cage variant (TRP1; α), Full Sequence Design (FSD; α + β), villin headpiece (villin; α), and a truncated FBP-28 WW-domain variant (2MWD; β) (the NEWCT-4P force field) and, subsequently, with an enhanced set of 9 proteins composed of TRP1, FSD, villin, 1BDC (α), 2I18 (α), 1QHK (α + β), 2N9L (α + β), 1E0L (β), and 2LX7 (β) (the NEWCT-9P force field). The NEWCT-9P force field performed better than NEWCT-4P in a blind-prediction-like test with a set of 26 proteins not used in calibration and outperformed, in a test with 76 proteins, the most advanced OPT-WTFSA-2 version of UNRES with former backbone-local and backbone-correlation terms that contained more energy terms and more optimizable parameters. The NEWCT-9P force field reproduced the bimodal distribution of backbone-virtual-bond angles in the simulated structures, as observed in experimental protein structures.
Collapse
Affiliation(s)
- Adam Liwo
- Faculty of Chemistry, University of Gdańsk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdańsk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - Agnieszka G Lipska
- Faculty of Chemistry, University of Gdańsk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - InSuk Joung
- School of Computational Sciences, Korea Institute for Advanced Study, 87 Hoegiro, Dongdaemun-gu, 130-722 Seoul, South Korea
| | - Wioletta Żmudzińska
- Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, ul. Abrahama 58, 80-307 Gdańsk, Poland
| | - Anna Hałabis
- Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, ul. Abrahama 58, 80-307 Gdańsk, Poland
| | - Stanisław Ołdziej
- Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, ul. Abrahama 58, 80-307 Gdańsk, Poland
| |
Collapse
|
26
|
A structural model of the immune checkpoint CD160-HVEM complex derived from HDX-mass spectrometry and molecular modeling. Oncotarget 2019; 10:536-550. [PMID: 30728903 PMCID: PMC6355189 DOI: 10.18632/oncotarget.26570] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2018] [Accepted: 12/10/2018] [Indexed: 11/25/2022] Open
Abstract
CD160 is a T cell coinhibitory molecule that interacts with the herpes virus entry mediator (HVEM) on antigen-presenting cells to provide an inhibitory signal to T cells. To date, the structure of CD160 and its complex with HVEM are unknown. Here, we have identified the fragments of CD160 interacting with HVEM using ELISA tests, hydrogen/deuterium studies, affinity chromatography and mass spectrometry (MS). By combining hydrogen/deuterium exchange and mass spectrometry (HDX-MS) we obtained key information about the tertiary structure of CD160, predicting the 3D structure of the CD160–HVEM complex. Our results provide insights into the molecular architecture of this complex, serving as a useful basis for designing inhibitors for future immunotherapies.
Collapse
|
27
|
Faraggi E, Krupa P, Mozolewska MA, Liwo A, Kloczkowski A. Reoptimized UNRES Potential for Protein Model Quality Assessment. Genes (Basel) 2018; 9:genes9120601. [PMID: 30513992 PMCID: PMC6315818 DOI: 10.3390/genes9120601] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2018] [Revised: 11/25/2018] [Accepted: 11/27/2018] [Indexed: 11/16/2022] Open
Abstract
Ranking protein structure models is an elusive problem in bioinformatics. These models are evaluated on both the degree of similarity to the native structure and the folding pathway. Here, we simulated the use of the coarse-grained UNited RESidue (UNRES) force field as a tool to choose the best protein structure models for a given protein sequence among a pool of candidate models, using server data from the CASP11 experiment. Because the original UNRES was optimized for Molecular Dynamics simulations, we reoptimized UNRES using a deep feed-forward neural network, and we show that introducing additional descriptive features can produce better results. Overall, we found that the reoptimized UNRES performs better in selecting the best structures and tracking protein unwinding from its native state. We also found a relatively poor correlation between UNRES values and the model’s Template Modeling Score (TMS). This is remedied by reoptimization. We discuss some cases where our reoptimization procedure is useful.
Collapse
Affiliation(s)
- Eshel Faraggi
- Research and Information Systems, LLC, Indianapolis, IN 46240, USA.
- Department of Physics, Indiana University Purdue University Indianapolis, Indianapolis, IN 46202, USA.
- Battelle Center for Mathematical Medicine, Nationwide Children's Hospital, Columbus, OH 43215, USA.
| | - Pawel Krupa
- Battelle Center for Mathematical Medicine, Nationwide Children's Hospital, Columbus, OH 43215, USA.
- Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, PL-02-668 Warsaw, Poland.
| | - Magdalena A Mozolewska
- Battelle Center for Mathematical Medicine, Nationwide Children's Hospital, Columbus, OH 43215, USA.
- Institute of Computer Science, Polish Academy of Sciences, ul. Jana Kazimierza 5, 01-248 Warszawa, 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, Korea.
| | - Andrzej Kloczkowski
- Battelle Center for Mathematical Medicine, Nationwide Children's Hospital, Columbus, OH 43215, USA.
- Department of Pediatrics, The Ohio State University, Columbus, OH 43215, USA.
- Kavli Institute for Theoretical Physics China, Chinese Academy of Sciences, Beijing 100190, China.
| |
Collapse
|
28
|
Sieradzan AK, Giełdoń A, Yin Y, He Y, Scheraga HA, Liwo A. A new protein nucleic-acid coarse-grained force field based on the UNRES and NARES-2P force fields. J Comput Chem 2018; 39:2360-2370. [PMID: 30306573 DOI: 10.1002/jcc.25571] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2018] [Revised: 06/29/2018] [Accepted: 08/08/2018] [Indexed: 12/22/2022]
Abstract
Based on the coarse-grained UNRES and NARES-2P models of proteins and nucleic acids, respectively, developed in our laboratory, in this work we have developed a coarse-grained model of systems containing proteins and nucleic acids. The UNRES and NARES-2P effective energy functions have been applied to the protein and nucleic-acid components of a system, respectively, while protein-nucleic-acid interactions have been described by the respective coarse-grained potentials developed in our recent work (Yin et al., J. Chem Theory Comput. 2015, 11, 1792). The Debye-Hückel screening has been applied to the electrostatic-interaction energy between the phosphate groups and charged amino-acid side chains. The model has been integrated into the UNRES package for coarse-grained molecular dynamics simulations of proteins and the implementation has been tested for energy conservation in microcanonical molecular dynamics runs and for temperature conservation in canonical molecular dynamics runs. Two case studies were performed: (i) the dynamics of the Ku protein heterodimer bound to DNA, for which it was found that the Ku70/Ku80 protein complex plays an active role in DNA repairing and (ii) conformational changes of the multiple antibiotic resistance (MarA) protein occurring during DNA binding, for which the functionally important motions occurring during this process were identified. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Adam K Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Artur Giełdoń
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| | - Yanping Yin
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, 14853-1301
| | - Yi He
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, 14853-1301.,School of Engineering, University of California Merced, 5200 N. Lake Road, Merced, California, 95343
| | - Harold A Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, 14853-1301
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308, Gdańsk, Poland
| |
Collapse
|
29
|
Rojas AV, Maisuradze GG, Scheraga HA. Dependence of the Formation of Tau and Aβ Peptide Mixed Aggregates on the Secondary Structure of the N-Terminal Region of Aβ. J Phys Chem B 2018; 122:7049-7056. [PMID: 29940109 DOI: 10.1021/acs.jpcb.8b04647] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
One of the hallmarks of Alzheimer's disease is the formation of aggregates of the tau protein, a process that can be facilitated by the presence of fibrils formed by the amyloid β peptide (Aβ). However, the mechanism that triggers tau aggregation is still a matter of debate. The effect of Aβ40 fibrils on the aggregation of the repeat domain of tau (TauRD) is investigated here by employing coarse-grained molecular dynamics simulations. The results indicate that the repeat domain of tau has a high affinity for Aβ40 fibrils, with the 261GSTENLK267 fragment of tau driving TauRD toward the 16KLVFFA21 fragment in Aβ40. Monomeric Aβ40, in which the 16KLVFFA21 fragment is rarely found in an extended conformation (as in the fibril), has a low affinity for the TauRD, indicating that the ability of Aβ40 fibrils to bind to the TauRD depends on the 16KLVFFA21 fragment of Aβ adopting an extended conformation.
Collapse
Affiliation(s)
- Ana V Rojas
- Schrödinger, Inc. , 120 West 45th Street , New York , New York 10036 , United States
| | - Gia G Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology , Cornell University , Ithaca , New York 14853 , United States
| | - Harold A Scheraga
- Baker Laboratory of Chemistry and Chemical Biology , Cornell University , Ithaca , New York 14853 , United States
| |
Collapse
|
30
|
Czaplewski C, Karczyńska A, Sieradzan AK, Liwo A. UNRES server for physics-based coarse-grained simulations and prediction of protein structure, dynamics and thermodynamics. Nucleic Acids Res 2018; 46:W304-W309. [PMID: 29718313 PMCID: PMC6031057 DOI: 10.1093/nar/gky328] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2018] [Revised: 04/12/2018] [Accepted: 04/18/2018] [Indexed: 12/24/2022] Open
Abstract
A server implementation of the UNRES package (http://www.unres.pl) for coarse-grained simulations of protein structures with the physics-based UNRES model, coined a name UNRES server, is presented. In contrast to most of the protein coarse-grained models, owing to its physics-based origin, the UNRES force field can be used in simulations, including those aimed at protein-structure prediction, without ancillary information from structural databases; however, the implementation includes the possibility of using restraints. Local energy minimization, canonical molecular dynamics simulations, replica exchange and multiplexed replica exchange molecular dynamics simulations can be run with the current UNRES server; the latter are suitable for protein-structure prediction. The user-supplied input includes protein sequence and, optionally, restraints from secondary-structure prediction or small x-ray scattering data, and simulation type and parameters which are selected or typed in. Oligomeric proteins, as well as those containing D-amino-acid residues and disulfide links can be treated. The output is displayed graphically (minimized structures, trajectories, final models, analysis of trajectory/ensembles); however, all output files can be downloaded by the user. The UNRES server can be freely accessed at http://unres-server.chem.ug.edu.pl.
Collapse
Affiliation(s)
- Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Agnieszka Karczyńska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam K Sieradzan
- 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
|
31
|
Sieradzan AK, Mozolewska MA. Extension of coarse-grained UNRES force field to treat carbon nanotubes. J Mol Model 2018; 24:121. [PMID: 29700628 PMCID: PMC5920012 DOI: 10.1007/s00894-018-3656-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2017] [Accepted: 04/13/2018] [Indexed: 11/25/2022]
Abstract
Carbon nanotubes (CNTs) have recently received considerable attention because of their possible applications in various branches of nanotechnology. For their cogent application, knowledge of their interactions with biological macromolecules, especially proteins, is essential and computer simulations are very useful for such studies. Classical all-atom force fields limit simulation time scale and size of the systems significantly. Therefore, in this work, we implemented CNTs into the coarse-grained UNited RESidue (UNRES) force field. A CNT is represented as a rigid infinite-length cylinder which interacts with a protein through the Kihara potential. Energy conservation in microcanonical coarse-grained molecular dynamics simulations and temperature conservation in canonical simulations with UNRES containing the CNT component have been verified. Subsequently, studies of three proteins, bovine serum albumin (BSA), soybean peroxidase (SBP), and α-chymotrypsin (CT), with and without CNTs, were performed to examine the influence of CNTs on the structure and dynamics of these proteins. It was found that nanotubes bind to these proteins and influence their structure. Our results show that the UNRES force field can be used for further studies of CNT-protein systems with 3-4 order of magnitude larger timescale than using regular all-atom force fields. Graphical abstract Bovine serum albumin (BSA), soybean peroxidase (SBP), and α-chymotrypsin (CT), with and without CNTsᅟ.
Collapse
Affiliation(s)
- Adam K Sieradzan
- Faculty of Chemistry, University of Gdansk, ul. Wita Stwosza 63, 80-308, Gdansk, Poland
| | - Magdalena A Mozolewska
- Faculty of Chemistry, University of Gdansk, ul. Wita Stwosza 63, 80-308, Gdansk, Poland.
- Institute of Computer Science, Polish Academy of Sciences, ul. Jana Kazimierza 5, 01-248, Warsaw, Poland.
| |
Collapse
|
32
|
Lysosomal enzyme tripeptidyl peptidase 1 destabilizes fibrillar Aβ by multiple endoproteolytic cleavages within the β-sheet domain. Proc Natl Acad Sci U S A 2018; 115:1493-1498. [PMID: 29378960 DOI: 10.1073/pnas.1719808115] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
Accumulation of amyloid-beta (Aβ), which is associated with Alzheimer's disease, can be caused by excess production or insufficient clearance. Because of its β-sheet structure, fibrillar Aβ is resistant to proteolysis, which would contribute to slow degradation of Aβ plaques in vivo. Fibrillar Aβ can be internalized by microglia, which are the scavenger cells of the brain, but the fibrils are degraded only slowly in microglial lysosomes. Cathepsin B is a lysosomal protease that has been shown to proteolyze fibrillar Aβ. Tripeptidyl peptidase 1 (TPP1), a lysosomal serine protease, possesses endopeptidase activity and has been shown to cleave peptides between hydrophobic residues. Herein, we demonstrate that TPP1 is able to proteolyze fibrillar Aβ efficiently. Mass spectrometry analysis of peptides released from fibrillar Aβ digested with TPP1 reveals several endoproteolytic cleavages including some within β-sheet regions that are important for fibril formation. Using molecular dynamics simulations, we demonstrate that these cleavages destabilize fibrillar β-sheet structure. The demonstration that TPP1 can degrade fibrillar forms of Aβ provides insight into the turnover of fibrillar Aβ and may lead to new therapeutic methods to increase degradation of Aβ plaques.
Collapse
|
33
|
Karczyńska AS, Mozolewska MA, Krupa P, Giełdoń A, Liwo A, Czaplewski C. Prediction of protein structure with the coarse-grained UNRES force field assisted by small X-ray scattering data and knowledge-based information. Proteins 2017; 86 Suppl 1:228-239. [PMID: 29134679 DOI: 10.1002/prot.25421] [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: 07/20/2017] [Revised: 11/09/2017] [Accepted: 11/10/2017] [Indexed: 11/09/2022]
Abstract
A new approach to assisted protein-structure prediction has been proposed, which is based on running multiplexed replica exchange molecular dynamics simulations with the coarse-grained UNRES force field with restraints derived from knowledge-based models and distance distribution from small angle X-ray scattering (SAXS) measurements. The latter restraints are incorporated into the target function as a maximum-likelihood term that guides the shape of the simulated structures towards that defined by SAXS. The approach was first verified with the 1KOY protein, for which the distance distribution was calculated from the experimental structure, and subsequently used to predict the structures of 11 data-assisted targets in the CASP12 experiment. Major improvement of the GDT_TS was obtained for 2 targets, minor improvement for other 2 while, for 6 target GDT_TS deteriorated compared with that calculated for predictions without the SAXS data, partly because of assuming a wrong multimeric state (for Ts866) or because the crystal conformation was more compact than the solution conformation (for Ts942). Particularly good results were obtained for Ts909, in which use of SAXS data resulted in the selection of a correctly packed trimer and, subsequently, increased the GDT_TS of monomer prediction. It was found that running simulations with correct oligomeric state is essential for the success in SAXS-data-assisted prediction.
Collapse
Affiliation(s)
| | - Magdalena A Mozolewska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk, 80-308, Poland.,Institute of Computer Science, Polish Academy of Sciences, ul. Jana Kazimierza 5, Warsaw, 01-248, Poland
| | - Paweł Krupa
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk, 80-308, Poland.,Institute of Physics, Polish Academy of Sciences, Aleja Lotników 32/46, Warsaw, PL-02668, Poland
| | - Artur Giełdoń
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk, 80-308, Poland
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk, 80-308, Poland.,School of Computational Sciences, 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, Gdańsk, 80-308, Poland
| |
Collapse
|
34
|
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
|
35
|
Krupa P, Sieradzan AK, Mozolewska MA, Li H, Liwo A, Scheraga HA. Dynamics of Disulfide-Bond Disruption and Formation in the Thermal Unfolding of Ribonuclease A. J Chem Theory Comput 2017; 13:5721-5730. [DOI: 10.1021/acs.jctc.7b00724] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Paweł Krupa
- Baker
Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
| | - Adam K. Sieradzan
- Faculty
of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Magdalena A. Mozolewska
- Baker
Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
| | - Huiyu Li
- 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, 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
|
36
|
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
|
37
|
Ramezanghorbani F, Dalgicdir C, Sayar M. A multi-state coarse grained modeling approach for an intrinsically disordered peptide. J Chem Phys 2017; 147:094103. [DOI: 10.1063/1.5001087] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
38
|
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
|
39
|
Rojas A, Maisuradze N, Kachlishvili K, Scheraga HA, Maisuradze GG. Elucidating Important Sites and the Mechanism for Amyloid Fibril Formation by Coarse-Grained Molecular Dynamics. ACS Chem Neurosci 2017; 8:201-209. [PMID: 28095675 DOI: 10.1021/acschemneuro.6b00331] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
Fibrils formed by the β-amyloid (Aβ) peptide play a central role in the development of Alzheimer's disease. In this study, the principles governing their growth and stability are investigated by analyzing canonical and replica-exchange molecular dynamics trajectories of Aβ(9-40) fibrils. In particular, an unstructured monomer was allowed to interact freely with an Aβ fibril template. Trajectories were generated with the coarse-grained united-residue force field, and one- and two-dimensional free-energy landscapes (FELs) along the backbone virtual-bond angle θ and backbone virtual-bond-dihedral angle γ of each residue and principal components, respectively, were analyzed. Also, thermal unbinding (unfolding) of an Aβ peptide from the fibril template was investigated. These analyses enable us to illustrate the entire process of Aβ fibril elongation and to elucidate the key residues involved in it. Several different pathways were identified during the search for the fibril conformation by the monomer, which finally follows a dock-lock mechanism with two distinct locking stages. However, it was found that the correct binding, with native hydrogen bonds, of the free monomer to the fibril template at both stages is crucial for fibril elongation. In other words, if the monomer is incorrectly bound (with nonnative hydrogen bonds) to the fibril template during the first "docking" stage, it can remain attached to it for a long time before it dissociates and either attempts a different binding or allows another monomer to bind. This finding is consistent with an experimentally observed "stop-and-go" mechanism of fibril growth.
Collapse
Affiliation(s)
- Ana Rojas
- Baker Laboratory
of Chemistry
and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
| | - Nika Maisuradze
- Baker Laboratory
of Chemistry
and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
| | - Khatuna Kachlishvili
- Baker Laboratory
of Chemistry
and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
| | - Harold A. Scheraga
- Baker Laboratory
of Chemistry
and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
| | - Gia G. Maisuradze
- Baker Laboratory
of Chemistry
and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, United States
| |
Collapse
|
40
|
Sieradzan AK, Jakubowski R. Introduction of steered molecular dynamics into UNRES coarse-grained simulations package. J Comput Chem 2017; 38:553-562. [DOI: 10.1002/jcc.24685] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2016] [Revised: 11/08/2016] [Accepted: 11/09/2016] [Indexed: 01/21/2023]
Affiliation(s)
- Adam K. Sieradzan
- Faculty of Chemistry; University of Gdańsk; Wita Stwosza 63 Gdańsk 80-308 Poland
| | - Rafał Jakubowski
- Faculty of Physics, Astronomy and Informatics, Institute of Physics, Nicolaus Copernicus University; Grudziadzka 5 Torun 87-100 Poland
| |
Collapse
|
41
|
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
|
42
|
He M, Wang Z, Wang R, Zhang L, Jia Q. Preparation of Bio-Based Polyamide Elastomer by Using Green Plasticizers. Polymers (Basel) 2016; 8:E257. [PMID: 30974535 PMCID: PMC6432398 DOI: 10.3390/polym8070257] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2016] [Revised: 06/29/2016] [Accepted: 07/05/2016] [Indexed: 11/30/2022] Open
Abstract
The purpose of this work was to study the effects of three green plasticizers H₂O, glycerol, and soybean oil, on the properties of bio-based BDIS polyamides. The BDIS polyamides synthesized from the following biomass monomers: 1,4-butanediamine (BD), 1,10-decanediamine (DD), itaconic acid (IA), and sebacic acid (SA). It is interesting to note that the amorphous BDIS (IA-80%) polyamide was changed from the glassy state to the rubbery state after water soaking and induced crystallization at the same time. The H₂O-plasticized non-crosslinked BDIS (IA-80%) polyamides can be very useful for the preparation of physical water gel. The glycerol- and soybean oil-plasticized BDIS (IA-80%) polyamides displayed excellent toughness. The plasticized BDIS (IA-80%) polyamides were characterized by Fouriertransform infrared spectroscopy (FTIR), differential scanning calorimetry (DSC), thermogravimetric analysis (TGA), mechanical testing, and X-ray diffraction (XRD).
Collapse
Affiliation(s)
- Miaomiao He
- Key Laboratory of Beijing City on Preparation and Processing of Novel Polymer Materials, Beijing University of Chemical Technology, Beijing 100029, China.
| | - Zhao Wang
- Department of Chemistry, South Dakota School of Mines and Technology, Rapid, SD 57701, USA.
| | - Runguo Wang
- Key Laboratory of Beijing City on Preparation and Processing of Novel Polymer Materials, Beijing University of Chemical Technology, Beijing 100029, China.
| | - Liqun Zhang
- Key Laboratory of Beijing City on Preparation and Processing of Novel Polymer Materials, Beijing University of Chemical Technology, Beijing 100029, China.
| | - Qingxiu Jia
- College of Materials Science and Engineering, Beijing Institute of Fashion Technology, Beijing 100029, China.
| |
Collapse
|
43
|
Dai J, Niemi AJ, He J, Sieradzan A, Ilieva N. Bloch spin waves and emergent structure in protein folding with HIV envelope glycoprotein as an example. Phys Rev E 2016; 93:032409. [PMID: 27078392 DOI: 10.1103/physreve.93.032409] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2015] [Indexed: 06/05/2023]
Abstract
We inquire how structure emerges during the process of protein folding. For this we scrutinize collective many-atom motions during all-atom molecular dynamics simulations. We introduce, develop, and employ various topological techniques, in combination with analytic tools that we deduce from the concept of integrable models and structure of discrete nonlinear Schrödinger equation. The example we consider is an α-helical subunit of the HIV envelope glycoprotein gp41. The helical structure is stable when the subunit is part of the biological oligomer. But in isolation, the helix becomes unstable, and the monomer starts deforming. We follow the process computationally. We interpret the evolving structure both in terms of a backbone based Heisenberg spin chain and in terms of a side chain based XY spin chain. We find that in both cases the formation of protein supersecondary structure is akin the formation of a topological Bloch domain wall along a spin chain. During the process we identify three individual Bloch walls and we show that each of them can be modelled with a precision of tenths to several angstroms in terms of a soliton solution to a discrete nonlinear Schrödinger equation.
Collapse
Affiliation(s)
- Jin Dai
- School of Physics, Beijing Institute of Technology, Beijing 100081, People's Republic of China
| | - Antti J Niemi
- School of Physics, Beijing Institute of Technology, Beijing 100081, People's Republic of China
- 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
| | - Jianfeng He
- School of Physics, Beijing Institute of Technology, Beijing 100081, People's Republic of China
| | - Adam Sieradzan
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Nevena Ilieva
- Institute of Information and Communication Technologies, Bulgarian Academy of Sciences, 25A, Acad. G. Bonchev Str., Sofia 1113, Bulgaria
| |
Collapse
|
44
|
Abstract
The origins of formation of an intermediate state involved in amyloid formation and ways to prevent it are illustrated with the example of the Formin binding protein 28 (FBP28) WW domain, which folds with biphasic kinetics. Molecular dynamics of protein folding trajectories are used to examine local and global motions and the time dependence of formation of contacts between C(α)s and C(β)s of selected pairs of residues. Focus is placed on the WT FBP28 WW domain and its six mutants (L26D, L26E, L26W, E27Y, T29D, and T29Y), which have structures that are determined by high-resolution NMR spectroscopy. The origins of formation of an intermediate state are elucidated, viz. as formation of hairpin 1 by a hydrophobic collapse mechanism causing significant delay of formation of both hairpins, especially hairpin 2, which facilitates the emergence of an intermediate state. It seems that three-state folding is a major folding scenario for all six mutants and WT. Additionally, two-state and downhill folding scenarios were identified in ∼ 15% of the folding trajectories for L26D and L26W, in which both hairpins are formed by the Matheson-Scheraga mechanism much faster than in three-state folding. These results indicate that formation of hairpins connecting two antiparallel β-strands determines overall folding. The correlations between the local and global motions identified for all folding trajectories lead to the identification of the residues making the main contributions in the formation of the intermediate state. The presented findings may provide an understanding of protein folding intermediates in general and lead to a procedure for their prevention.
Collapse
|
45
|
Zaborowski B, Jagieła D, Czaplewski C, Hałabis A, Lewandowska A, Żmudzińska W, Ołdziej S, Karczyńska A, Omieczynski C, Wirecki T, Liwo A. A Maximum-Likelihood Approach to Force-Field Calibration. J Chem Inf Model 2015; 55:2050-70. [DOI: 10.1021/acs.jcim.5b00395] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Bartłomiej Zaborowski
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Dawid Jagieła
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Cezary Czaplewski
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Anna Hałabis
- Laboratory
of Biopolymer Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Kładki
24, 80-922 Gdańsk, Poland
| | - Agnieszka Lewandowska
- Laboratory
of Biopolymer Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Kładki
24, 80-922 Gdańsk, Poland
| | - Wioletta Żmudzińska
- Laboratory
of Biopolymer Structure, Intercollegiate Faculty of Biotechnology, University of Gdańsk and Medical University of Gdańsk, Kładki
24, 80-922 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, Kładki
24, 80-922 Gdańsk, Poland
| | - Agnieszka Karczyńska
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Christian Omieczynski
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Tomasz Wirecki
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
| | - Adam Liwo
- Faculty
of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, 80-308 Gdańsk, Poland
- Center
for In Silico Protein Structure and School of Computational Sciences, Korea Institute for Advanced Study, 87 Hoegiro, Dongdaemun-gu, Seoul 130-722, Republic of Korea
| |
Collapse
|
46
|
Krupa P, Mozolewska MA, Joo K, Lee J, Czaplewski C, Liwo A. Prediction of Protein Structure by Template-Based Modeling Combined with the UNRES Force Field. J Chem Inf Model 2015; 55:1271-81. [DOI: 10.1021/acs.jcim.5b00117] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Paweł Krupa
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
| | | | | | | | - Cezary Czaplewski
- 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
|
47
|
Abstract
This is a tour of a physical chemist through 65 years of protein chemistry from the time when emphasis was placed on the determination of the size and shape of the protein molecule as a colloidal particle, with an early breakthrough by James Sumner, followed by Linus Pauling and Fred Sanger, that a protein was a real molecule, albeit a macromolecule. It deals with the recognition of the nature and importance of hydrogen bonds and hydrophobic interactions in determining the structure, properties, and biological function of proteins until the present acquisition of an understanding of the structure, thermodynamics, and folding pathways from a linear array of amino acids to a biological entity. Along the way, with a combination of experiment and theoretical interpretation, a mechanism was elucidated for the thrombin-induced conversion of fibrinogen to a fibrin blood clot and for the oxidative-folding pathways of ribonuclease A. Before the atomic structure of a protein molecule was determined by x-ray diffraction or nuclear magnetic resonance spectroscopy, experimental studies of the fundamental interactions underlying protein structure led to several distance constraints which motivated the theoretical approach to determine protein structure, and culminated in the Empirical Conformational Energy Program for Peptides (ECEPP), an all-atom force field, with which the structures of fibrous collagen-like proteins and the 46-residue globular staphylococcal protein A were determined. To undertake the study of larger globular proteins, a physics-based coarse-grained UNited-RESidue (UNRES) force field was developed, and applied to the protein-folding problem in terms of structure, thermodynamics, dynamics, and folding pathways. Initially, single-chain and, ultimately, multiple-chain proteins were examined, and the methodology was extended to protein-protein interactions and to nucleic acids and to protein-nucleic acid interactions. The ultimate results led to an understanding of a variety of biological processes underlying natural and disease phenomena.
Collapse
|
48
|
Sieradzan AK. Introduction of periodic boundary conditions into UNRES force field. J Comput Chem 2015; 36:940-6. [DOI: 10.1002/jcc.23864] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2014] [Revised: 01/15/2015] [Accepted: 01/17/2015] [Indexed: 11/10/2022]
Affiliation(s)
- Adam K. Sieradzan
- Department of Physics and Astronomy; Uppsala University, Ångströmlaboratoriet, Lägerhyddsvägen; 1,751 20 Uppsala, Sweden; Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk Poland
| |
Collapse
|
49
|
Lipska AG, Sieradzan AK, Krupa P, Mozolewska MA, D’Auria S, Liwo A. Studies of conformational changes of an arginine-binding protein from Thermotoga maritima in the presence and absence of ligand via molecular dynamics simulations with the coarse-grained UNRES force field. J Mol Model 2015; 21:64. [DOI: 10.1007/s00894-015-2609-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2014] [Accepted: 02/08/2015] [Indexed: 11/30/2022]
|
50
|
Sieradzan AK, Krupa P, Scheraga HA, Liwo A, Czaplewski C. Physics-based potentials for the coupling between backbone- and side-chain-local conformational states in the UNited RESidue (UNRES) force field for protein simulations. J Chem Theory Comput 2015; 11:817-31. [PMID: 25691834 PMCID: PMC4327884 DOI: 10.1021/ct500736a] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
The UNited RESidue (UNRES) model of polypeptide chains is a coarse-grained model in which each amino-acid residue is reduced to two interaction sites, namely, a united peptide group (p) located halfway between the two neighboring α-carbon atoms (Cαs), which serve only as geometrical points, and a united side chain (SC) attached to the respective Cα. Owing to this simplification, millisecond molecular dynamics simulations of large systems can be performed. While UNRES predicts overall folds well, it reproduces the details of local chain conformation with lower accuracy. Recently, we implemented new knowledge-based torsional potentials (Krupa et al. J. Chem. Theory Comput. 2013, 9, 4620–4632) that depend on the virtual-bond dihedral angles involving side chains: Cα···Cα···Cα···SC (τ(1)), SC···Cα···Cα···Cα (τ(2)), and SC···Cα···Cα···SC (τ(3)) in the UNRES force field. These potentials resulted in significant improvement of the simulated structures, especially in the loop regions. In this work, we introduce the physics-based counterparts of these potentials, which we derived from the all-atom energy surfaces of terminally blocked amino-acid residues by Boltzmann integration over the angles λ(1) and λ(2) for rotation about the Cα···Cα virtual-bond angles and over the side-chain angles χ. The energy surfaces were, in turn, calculated by using the semiempirical AM1 method of molecular quantum mechanics. Entropy contribution was evaluated with use of the harmonic approximation from Hessian matrices. One-dimensional Fourier series in the respective virtual-bond-dihedral angles were fitted to the calculated potentials, and these expressions have been implemented in the UNRES force field. Basic calibration of the UNRES force field with the new potentials was carried out with eight training proteins, by selecting the optimal weight of the new energy terms and reducing the weight of the regular torsional terms. The force field was subsequently benchmarked with a set of 22 proteins not used in the calibration. The new potentials result in a decrease of the root-mean-square deviation of the average conformation from the respective experimental structure by 0.86 Å on average; however, improvement of up to 5 Å was observed for some proteins.
Collapse
Affiliation(s)
- Adam K. Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-180 Gdańsk, Poland
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, N.Y., 14853-1301, U.S.A
| | - Paweł Krupa
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-180 Gdańsk, Poland
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, N.Y., 14853-1301, U.S.A
| | - Harold A. Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, N.Y., 14853-1301, U.S.A
| | - Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-180 Gdańsk, Poland
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-180 Gdańsk, Poland
| |
Collapse
|