Probing the stereospecificity of tyrosyl- and glutaminyl-tRNA synthetase with molecular dynamics

The Dual Role of the 2'-OH Group of A76 tRNATyr in the Prevention of d-tyrosine Mistranslation

Aminoacyl-tRNA-synthetases are crucial enzymes for initiation step of translation. Possessing editing activity, they protect living cells from misincorporation of non-cognate and non-proteinogenic amino acids into proteins. Tyrosyl-tRNA synthetase (TyrRS) does not have such editing properties, but it shares weak stereospecificity in recognition of d-/l-tyrosine (Tyr). Nevertheless, an additional enzyme, d-aminoacyl-tRNA-deacylase (DTD), exists to overcome these deficiencies. The precise catalytic role of hydroxyl groups of the tRNA^Tyr A76 in the catalysis by TyrRS and DTD remained unknown. To address this issue, [³²P]-labeled tRNA^Tyr substrates have been tested in aminoacylation and deacylation assays. TyrRS demonstrates similar activity in charging the 2' and 3'-OH groups of A76 with l-Tyr. This synthetase can effectively use both OH groups as primary sites for aminoacylation with l-Tyr, but demonstrates severe preference toward 2'-OH, in charging with d-Tyr. In both cases, the catalysis is not substrate-assisted: neither the 2'-OH nor the 3'-OH group assists catalysis. In contrast, DTD catalyzes deacylation of d-Tyr-tRNA^Tyr specifically from the 3'-OH group, while the 2'-OH assists in this hydrolysis.

Comparing pairwise-additive and many-body generalized Born models for acid/base calculations and protein design

Generalized Born (GB) solvent models are common in acid/base calculations and protein design. With GB, the interaction between a pair of solute atoms depends on the shape of the protein/solvent boundary and, therefore, the positions of all solute atoms, so that GB is a many-body potential. For compute-intensive applications, the model is often simplified further, by introducing a mean, native-like protein/solvent boundary, which removes the many-body property. We investigate a method for both acid/base calculations and protein design that uses Monte Carlo simulations in which side chains can explore rotamers, bind/release protons, or mutate. The fluctuating protein/solvent dielectric boundary is treated in a way that is numerically exact (within the GB framework), in contrast to a mean boundary. Its originality is that it captures the many-body character while retaining the residue-pairwise complexity given by a fixed boundary. The method is implemented in the Proteus protein design software. It yields a slight but systematic improvement for acid/base constants in nine proteins and a significant improvement for the computational design of three PDZ domains. It eliminates a source of model uncertainty, which will facilitate the analysis of other model limitations. © 2017 Wiley Periodicals, Inc.