Orientation and distance dependent chiral discrimination in the first step of the aminoacylation reaction: integrated molecular orbital and semi-empirical method (ONIOM) based calculation

Roles of the Active Site Zn(II) and Residues in Substrate Discrimination by Threonyl-tRNA Synthetase: An MD and QM/MM Investigation

Threonyl-tRNA synthetase (ThrRS) is a Zn(II) containing enzyme that catalyzes the activation of threonine and its subsequent transfer to the cognate tRNA. This process is accomplished with remarkable fidelity, with ThrRS being able to discriminate its cognate substrate from similar analogues such as serine and valine. Molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) methods have been used to elucidate the role of Zn(II) in the aminoacylation mechanism of ThrRS. More specifically, the role of Zn(II) and active site residues in ThrRS's ability to discriminate between its cognate substrate l-threonine and the noncognate l-serine, l-valine, and d-threonine has been examined. The present results suggest that a role of the Zn(II) ion, with its Lewis acidity, is to facilitate deprotonation of the side chain hydroxyl groups of the aminoacyl moieties of cognate Thr-AMP and noncognate Ser-AMP substrates. In their deprotonated forms, these substrates are able to adopt a conformation preferable for aminoacyl transfer from aa-AMP onto the Ado-3'OH of the tRNA^Thr cosubstrate. Relative to the neutral substrates, when the substrates are first deprotonated with the assistance of the Zn(II) ion, the barrier for the rate-limiting step is decreased significantly by 42.0 and 39.2 kJ mol^-1 for l-Thr-AMP and l-Ser-AMP, respectively. An active site arginyl also plays a key role in stabilizing the buildup of negative charge on the substrate's bridging phosphate oxygen during the mechanism. For the enantiomeric substrate analogue d-Thr-AMP, product formation is highly disfavored, and as a result, the reverse reaction has a very low barrier of 16.0 kJ mol^-1.

Natural Mechanism of Origination and Conservation of Monochirality of Amino Acids

The proteins belonging to different organisms as well as the natural amino acids as incorporated into proteins occurring free in tissues and lymphatic fluids are, as a rule, levorotatory. This article contains a simple explanation of this structural discrimination and of the extremely slow racemization of natural organic substances.

A three-layer ONIOM model for the outside binding of cationic porphyrins and nucleotide pair DNA

In this work we investigated the outside binding mode between a cationic porphyrin and a nucleotide pair of DNA, adenine-thymine and guanine-cytosine, in a supramolecular assembly. We used two structural models (semi-extended, extended) that differ in the size of porphyrin, two kinds of theoretical methods: a three layer ONIOM (B3LYP/6-31G(d)/PM3/UFF), and DFT B3LYP/6-31G(d,p), and three cationic porphyrins. ONIOM method was first tested on the semi-extended model that was calculated in four environments: gas phase, solution phase using an explicit solvent model (H(2)O), in the presence of a sodium cation (Na(+)) and in both (H(2)O + Na(+)). From interaction energy results, we found that the affinity of the cationic substituent by the adenine nucleotide is favored upon the thymine nucleotide. The extended model that considers the whole porphyrin was applied in the gas phase to the four nucleotides. All the cationic porphyrins showed affinity by the nucleotides in the order adenine > guanine > thymine > cytosine. The interaction energy values for outside binding showed a strong porphyrin-nucleotide interaction (≈-90 kcal mol(-1)), that slightly varies between the nucleotides suggesting that this kind of cationic porphyrin has a little selectivity for some of them. We also found that the effect of the nature of the cationic substituent (chain length) in the porphyrin on the outside binding is small (≈2-13 kcal mol(-1)). Coherence between the results showed that ONIOM is a useful tool to get a reasonable molecular geometry to be used as a starting point in calculations of density functional theory.

Mechanism of the activation step of the aminoacylation reaction: a significant difference between class I and class II synthetases

In the present work we report, for the first time, a novel difference in the molecular mechanism of the activation step of aminoacylation reaction between the class I and class II aminoacyl tRNA synthetases (aaRSs). The observed difference is in the mode of nucleophilic attack by the oxygen atom of the carboxylic group of the substrate amino acid (AA) to the αP atom of adenosine triphosphate (ATP). The syn oxygen atom of the carboxylic group attacks the α-phosphorous atom (αP) of ATP in all class I aaRSs (except TrpRS) investigated, while the anti oxygen atom attacks in the case of class II aaRSs. The class I aaRSs investigated are GluRS, GlnRS, TyrRS, TrpRS, LeuRS, ValRS, IleRS, CysRS, and MetRS and class II aaRSs investigated are HisRS, LysRS, ProRS, AspRS, AsnRS, AlaRS, GlyRS, PheRS, and ThrRS. The variation of the electron density at bond critical points as a function of the conformation of the attacking oxygen atom measured by the dihedral angle ψ (C(α)-C') conclusively proves this. The result shows that the strength of the interaction of syn oxygen and αP is stronger than the interaction with the anti oxygen for class I aaRSs. This indicates that the syn oxygen is the most probable candidate for the nucleophilic attack in class I aaRSs. The result is further supported by the computation of the variation of the nonbonded interaction energies between αP atom and anti oxygen as well as syn oxygen in class I and II aaRSs, respectively. The difference in mechanism is explained based on the analysis of the electrostatic potential of the AA and ATP which shows that the relative arrangement of the ATP with respect to the AA is opposite in class I and class II aaRSs, which is correlated with the organization of the active site in respective aaRSs. A comparative study of the reaction mechanisms of the activation step in a class I aaRS (Glutaminyl tRNA synthetase) and in a class II aaRS (Histidyl tRNA synthetase) is carried out by the transition state analysis. The atoms in molecule analysis of the interaction between active site residues or ions and substrates are carried out in the reactant state and the transition state. The result shows that the observed novel difference in the mechanism is correlated with the organizations of the active sites of the respective aaRSs. The result has implication in understanding the experimentally observed different modes of tRNA binding in the two classes of aaRSs.

Aminoacylation reaction in the histidyl-tRNA synthetase: fidelity mechanism of the activation step

Aminoacylation is a vital step of natural biosynthesis of peptide. Correct aminoacylation is a necessary prerequisite for the elimination of noncognate amino acids such as D-amino acids. In the present work, we studied the fidelity mechanism of histidine (His) activation (first step of aminoacylation reaction) using a combined quantum mechanical/semiempirical method based on a model of crystal structure of the oligomeric complex of histidyl-tRNA synthetase (HisRS) from Escherichia coli. The study of the variation in the energy during the mutual approach of the His and ATP to form adenylate shows that the surrounding nanospace of synthetase confines the reactants (L-His and ATP) and proximally places in a geometry suitable for the in-line nucleophilic attack. The significantly higher energy of the energy surface of the model containing D-His is due to unfavorable interaction of D-His with ATP and surrounding residues. This indicates that the network of interaction (principally electrostatic) is highly unfavorable when D-amino acid is incorporated. The reorganization of the surrounding nanospace can lower the unfavorable nature of the intermolecular energy surface of D-His and surrounding residues. However, such a rearrangement requires large-scale structural reorganization of the synthetase structure and is unfavorable. The variation in the bond angles and distances in going from the reactant state to the product state via transition state confirms the mechanism of nucleophilic attack and concomitant inversion of oxygen atoms around alpha-phosphorus (alpha-P). Calculation of the electrostatic potential indicates that in addition to the Mg(2+) the Arg residues in the active site facilitate the nucleophilic attack by reducing the negative charge distributed over the oxygen atoms attached to the alpha-P of ATP. Arg 259 residue has a role similar to that played by the two Mg(2+) cations as this residue is in close proximity of the alpha-P of ATP. Arg 113 also facilitates the reduction of the negative charge on the other side of the reaction center. The favorable electrostatic interaction of the Arg 259 with ATP and His is also concluded from the calculation of the binding energy. The Arg 259 anchors the carboxylic acid group of His and the oxygen atom of the alpha-phosphate group during the progress of reaction. Consequently, Arg 259 plays an important catalytic role in the activation step rather than merely reducing the negative charge density over the ATP.