Biopolym. Cell. 2018; 34(3):239-248.
Computational approaches for parameterization of aminoacyl-tRNA synthetase substrates
1Rayevsky O. V., 1Tukalo M. A.
  1. Institute of Molecular Biology and Genetics, NAS of Ukraine
    150, Akademika Zabolotnoho Str., Kyiv, Ukraine, 03143


Aim. To parameterize a modified chained residue and use a newborn topology for molecular dynamics simulation. Methods. To deal with the problem, a series of ab initio and semi-empirical methods were combined. The RESP (Restrained ElectroStatic Potential) program, which fits molecular electrostatic potential (MEP) at molecular surfaces using an atom-centered point charge model. All parameters were quantum mechanically calculated and processed with R.E.D.III server. Results. The method of molecular dynamics has potential advantages, like its capability to explore large systems, and disadvantages, like not being feasible to run on fly without a preliminary prepared topologies for identification of each molecule. In an attempt to find a balance between both features speed and accuracy and apply the approach in a computational study, a functional mechanism of prolyl-tRNA synthetase from E. faecalis was investigated. In addition, a well validated protocol of topology preparation for non-canonical structure was developed. Conclusions. Computational approaches like molecular dynamics simulation and molecular docking had now become a strong in silico methods to study biological processes. The major benefits of this methods are expensiveness and speed. It also a strong competitor of quantum modeling approach. However, there is a need to include new structures that are not exist in the GROMACS library is associated with the growing use of different types of modified amino and nucleic acids (DNA derivatives and RNAs) both in fundamental studies in molecular biology, molecular biophysics, etc., and in applied research related to the search for new drugs. The obtained data well correlate with the experimental data.
Keywords: force field, RESP, R.E.D.III, molecular dynamics, aminoacyl-tRNA synthetase


[1] Yadavalli SS, Ibba M. Quality control in aminoacyl-tRNA synthesis its role in translational fidelity. Adv Protein Chem Struct Biol. 2012;86:1-43.
[2] Schüttelkopf AW, van Aalten DM. PRODRG: a tool for high-throughput crystallography of protein-ligand complexes. Acta Crystallogr D Biol Crystallogr. 2004;60(Pt 8):1355-63.
[3] Vanommeslaeghe K, Raman EP, MacKerell AD Jr. Automation of the CHARMM General Force Field (CGenFF) II: assignment of bonded parameters and partial atomic charges. J Chem Inf Model. 2012;52(12):3155-68.
[4] Case DA, Darden TA, Cheatham TE, Simmerling CL, Wang J, Duke RE, et al. AMBER 12. University of California, San Francisco, 2012.
[5] Hess B, Kutzner C, van der Spoel D, Lindahl E. GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J Chem Theory Comput. 2008;4(3):435-47.
[6] Naray-Szabo G, Ferenczy GG, Molecular Electrostatics. Chem Rev. 1995; 95(4):829-87.
[7] Dupradeau FY, Pigache A, Zaffran T, Savineau C, Lelong R, Grivel N, Lelong D, Rosanski W, Cieplak P. The R.E.D. tools: advances in RESP and ESP charge derivation and force field library building. Phys Chem Chem Phys. 2010;12(28):7821-39.
[8] Crepin T, Yaremchuk A, Tukalo M, Cusack S. Structures of two bacterial prolyl-tRNA synthetases with and without a cis-editing domain. Structure. 2006;14(10):1511-25.
[9] Benkert P, Biasini M, Schwede T. Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics. 2011;27(3):343-50.
[10] Boyarshin KS, Priss AE, Rayevskiy AV, Ilchenko MM, Dubey IY, Kriklivyi IA, Yaremchuk AD, Tukalo MA. A new mechanism of post-transfer editing by aminoacyl-tRNA synthetases: catalysis of hydrolytic reaction by bacterial-type prolyl-tRNA synthetase. J Biomol Struct Dyn. 2017;35(3):669-682.
[11] Jensen F. Introduction to Computational Chemistry. Chichester: John Wiley & Sons; 2006; 620 p.
[12] Wang J, Cieplak P, Kollman PA. How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? J Comp Chem 2000;21(12):1049-74.
[13] Schmidt MW, Baldridge KK, Boatz JA, Elbert ST, Gordon MS, Jensen JH, Koseki S, Matsunaga N, Nguyen KA, Su S, Windus TL, Dupuis M. Montgomery JA Jr. General Atomic and Molecular Electronic Structure System. J Comput Chem. 1993; 14(11):1347-63.
[14] Cornell WD, Cieplak P, Bayly CI, Kollman PA. Application of RESP charges to calculate conformational energies, hydrogen bond energies and free energies of solvation. J Am Chem Soc. 1993;115(21):9620–31.
[15] Bayly CI, Cieplak P, Cornell W, Kollman PA. A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges: the RESP model. J Phys Chem. 1993; 97(40):10269-80.
[16] Kumar S, Das M, Hadad CM, Musier-Forsyth K. Substrate specificity of bacterial prolyl-tRNA synthetase editing domain is controlled by a tunable hydrophobic pocket. J Biol Chem. 2012; 287(5):3175–84.
[17] Rayevsky AV, Sharifi M, Tukalo MA. Molecular modeling and molecular dynamics simulation study of archaeal leucyl-tRNA synthetase in complex with different mischarged tRNA in editing conformation. J Mol Graph Model. 2017;76:289-295.
[18] Rayevsky A, Sharifi M, Tukalo M. A molecular dynamics simulation study of amino acid selectivity of LeuRS editing domain from Thermus thermophilus. J Mol Graph Model. 2018;84:74-81.
[19] Bullwinkle T,Lazazzera B, Ibba M. Quality control and infiltration of translation by amino acids outside of the genetic code. Ann Rev Genet. 2014; 48(1):149-66.