Improvement of the treatment of loop structures in the UNRES force field by inclusion of coupling between backbone- and side-chain-local conformational states
Paweł Krupa , Adam Sieradzan , S. Rackovsky , Maciej Baranowski , Stanisław Ołdziej , Harold A. Scheraga , Józef Adam Liwo , Cezary Czaplewski
AbstractThe UNited RESidue (UNRES) coarse-grained model of polypeptide chains, developed in our laboratory, enables us to carry out millisecond-scale molecular-dynamics simulations of large proteins e ff ectively. It performs well in ab initio predictions of protein structure, as demonstrated in the last Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP10). However, the resolution of the simulated structure is too coarse, especially in loop regions, which results from insu ffi cient speci fi city of the model of local interactions. To improve the representation of local interactions, in this work, we introduced new side-chain-backbone correlation potentials, derived from a statistical analysis of loop regions of 4585 proteins. To obtain su ffi cient statistics, we reduced the set of amino- acid-residue types to fi ve groups, derived in our earlier work on structurally optimized reduced alphabets [Solis, A. D.; Rackovsky, S. Proteins: Struct., Func., Bioinf. , 2000 , 38, 149 − 164], based on a statistical analysis of the properties of amino-acid structures. The new correlation potentials are expressed as one-dimensional Fourier series in the virtual-bond-dihedral angles involving side- chain centroids. The weight of these new terms was determined by a trial-and-error method, in which Multiplexed Replica Exchange Molecular Dynamics (MREMD) simulations were run on selected test proteins. The best average root-mean-square deviations (RMSDs) of the calculated structures from the experimental structures below the folding-transition temperatures were obtained with the weight of the new side-chain-backbone correlation potentials equal to 0.57. The resulting conformational ensembles were analyzed in detail by using the Weighted Histogram Analysis Method (WHAM) and Ward ’ s minimum-variance clustering. This analysis showed that the RMSDs from the experimental structures dropped by 0.5 Å on average, compared to simulations without the new terms, and the deviation of individual residues in the loop region of the computed structures from their counterparts in the experimental structures (after optimum superposition of the calculated and experimental structure) decreased by up to 8 Å. Consequently, the new terms improve the representation of local structure.
|Journal series||Journal of Chemical Theory and Computation, ISSN 1549-9618, e-ISSN 1549-9626, (A 40 pkt)|
|Publication size in sheets||0.6|
|Score|| = 40.0, 29-11-2018, ArticleFromJournal|
= 40.0, 29-11-2018, ArticleFromJournal
|Publication indicators||: 2013 = 5.31 (2) - 2013=6.205 (5)|
|Citation count*||33 (2018-12-18)|
* presented citation count is obtained through Internet information analysis and it is close to the number calculated by the Publish or Perish system.