1
|
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
|
2
|
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
|
3
|
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
|
4
|
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
|
5
|
Lee J, Lee J, Sasaki TN, Sasai M, Seok C, Lee J. De novo
protein structure prediction by dynamic fragment assembly and conformational space annealing. Proteins 2011; 79:2403-17. [DOI: 10.1002/prot.23059] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2010] [Revised: 03/24/2011] [Accepted: 04/12/2011] [Indexed: 12/25/2022]
|
6
|
Lee J, Kim SY, Lee J. Protein structure prediction based on fragment assembly and parameter optimization. Biophys Chem 2005; 115:209-14. [PMID: 15752606 DOI: 10.1016/j.bpc.2004.12.046] [Citation(s) in RCA: 35] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2004] [Revised: 11/09/2004] [Accepted: 12/10/2004] [Indexed: 11/28/2022]
Abstract
We propose a novel method for ab-initio prediction of protein tertiary structures based on the fragment assembly and global optimization. Fifteen residue long fragment libraries are constructed using the secondary structure prediction method PREDICT, and fragments in these libraries are assembled to generate full-length chains of a query protein. Tertiary structures of 50 to 100 conformations are obtained by minimizing an energy function for proteins, using the conformational space annealing method that enables one to sample diverse low-lying local minima of the energy. Then in order to enhance the performance of the prediction method, we optimize the linear parameters of the energy function, so that the native-like conformations become energetically more favorable than the non-native ones for proteins with known structures. We test the feasibility of the parameter optimization procedure by applying it to the training set consisting of three proteins: the 10-55 residue fragment of staphylococcal protein A (PDB ID 1bdd), a designed protein betanova, and 1fsd.
Collapse
Affiliation(s)
- Julian Lee
- Department of Bioinformatics and Life Science, Computer Aided Molecular Design Research Center, Bioinformatics and Molecular Design Technology Innovation Center, Soongsil University, Seoul 156-743, South Korea.
| | | | | |
Collapse
|
7
|
Lee J, Kim SY, Joo K, Kim I, Lee J. Prediction of protein tertiary structure using PROFESY, a novel method based on fragment assembly and conformational space annealing. Proteins 2004; 56:704-14. [PMID: 15281124 DOI: 10.1002/prot.20150] [Citation(s) in RCA: 66] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
A novel method for ab initio prediction of protein tertiary structures, PROFESY (PROFile Enumerating SYstem), is proposed. This method utilizes the secondary structure prediction information of a query sequence and the fragment assembly procedure based on global optimization. Fifteen-residue-long fragment libraries are constructed using the secondary structure prediction method PREDICT, and fragments in these libraries are assembled to generate full-length chains of a query protein. Tertiary structures of 50 to 100 conformations are obtained by minimizing an energy function for proteins, using the conformational space annealing method that enables one to sample diverse low-lying local minima of the energy. We apply PROFESY for benchmark tests to proteins with known structures to demonstrate its feasibility. In addition, we participated in CASP5 and applied PROFESY to four new-fold targets for blind prediction. The results are quite promising, despite the fact that PROFESY was in its early stages of development. In particular, PROFESY successfully provided us the best model-one structure for the target T0161.
Collapse
Affiliation(s)
- Julian Lee
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul, Korea
| | | | | | | | | |
Collapse
|
8
|
Abstract
Empirical force field-based studies of biological macromolecules are becoming a common tool for investigating their structure-activity relationships at an atomic level of detail. Such studies facilitate interpretation of experimental data and allow for information not readily accessible to experimental methods to be obtained. A large part of the success of empirical force field-based methods is the quality of the force fields combined with the algorithmic advances that allow for more accurate reproduction of experimental observables. Presented is an overview of the issues associated with the development and application of empirical force fields to biomolecular systems. This is followed by a summary of the force fields commonly applied to the different classes of biomolecules; proteins, nucleic acids, lipids, and carbohydrates. In addition, issues associated with computational studies on "heterogeneous" biomolecular systems and the transferability of force fields to a wide range of organic molecules of pharmacological interest are discussed.
Collapse
Affiliation(s)
- Alexander D Mackerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland 21201, USA.
| |
Collapse
|
9
|
Abstract
Extensive Monte Carlo folding simulations for four proteins of various structural classes are carried out, using a single continuous potential (united-residue force field). In all cases, collapse occurs at a very early stage, and proteins fold into their nativelike conformations at appropriate temperatures. We also observe that glassy transitions occur at low temperatures. The simulation results demonstrate that the folding mechanism is controlled not only by thermodynamic factors but also by kinetic factors: The way a protein folds into its native structure is also determined by the convergence point of early folding trajectories, which cannot be obtained by the free energy surface.
Collapse
Affiliation(s)
- Seung-Yeon Kim
- School of Computational Sciences, Korea Institute for Advanced Study, 207-43 Cheongryangri-dong, Dongdaemun-gu, Seoul 130-722, Korea
| | | | | |
Collapse
|
10
|
Lee J, Kim SY, Lee J. Design of a Protein Potential Energy Landscape by Parameter Optimization. J Phys Chem B 2004. [DOI: 10.1021/jp037076c] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Julian Lee
- Department of Bioinformatics and Life Sciences and Bioinformatics and Molecular Design Technology Innovation Center, Soongsil University, Seoul 156-743, Korea, and School of Computational Sciences, Korea Institute for Advanced Study, Seoul 130-722, Korea
| | - Seung-Yeon Kim
- Department of Bioinformatics and Life Sciences and Bioinformatics and Molecular Design Technology Innovation Center, Soongsil University, Seoul 156-743, Korea, and School of Computational Sciences, Korea Institute for Advanced Study, Seoul 130-722, Korea
| | - Jooyoung Lee
- Department of Bioinformatics and Life Sciences and Bioinformatics and Molecular Design Technology Innovation Center, Soongsil University, Seoul 156-743, Korea, and School of Computational Sciences, Korea Institute for Advanced Study, Seoul 130-722, Korea
| |
Collapse
|