Dynamics of the Active Sites of Dimeric Seryl tRNA Synthetase from Methanopyrus kandleri

A Multiple Proton Transfer Mechanism for the Charging Step of the Aminoacylation Reaction at the Active Site of Aspartyl tRNA Synthetase

Aspartyl-tRNA synthetase catalyzes the attachment of aspartic acid to its cognate tRNA by the aminoacylation reaction during the initiation of the protein biosynthesis process. In the second step of the aminoacylation reaction, known as the charging step, the aspartate moiety is transferred from aspartyl-adenylate to the 3'-OH of A76 of tRNA through a proton transfer process. We have investigated different pathways for the charging step through three separate QM/MM simulations combined with the enhanced sampling method of well-sliced metadynamics and found out the most feasible pathway for the reaction at the active site of the enzyme. In the charging reaction, both the phosphate group and the ammonium group after deprotonation can potentially act as a base for proton transfer in the substrate-assisted mechanism. We have considered three possible mechanisms involving different pathways of proton transfer, and only one of them is determined to be enzymatically feasible. The free energy landscape along reaction coordinates where the phosphate group acts as the general base showed that, in the absence of water, the barrier height is 52.6 kcal/mol. The free energy barrier is reduced to 39.7 kcal/mol when the active site water molecules are also treated quantum mechanically, thus allowing a water mediated proton transfer. The charging reaction involving the ammonium group of the aspartyl adenylate is found to follow a path where first a proton from the ammonium group moves to a water in the vicinity forming a hydronium ion (H₃O⁺) and NH₂ group. The hydronium ion subsequently passes the proton to the Asp233 residue, thus minimizing the chance of back proton transfer from hydronium to the NH₂ group. The neutral NH₂ group subsequently takes the proton from the O3' of A76 with a free energy barrier of 10.7 kcal/mol. In the next step, the deprotonated O3' makes a nucleophilic attack to the carbonyl carbon forming a tetrahedral transition state with a free energy barrier of 24.8 kcal/mol. Thus, the present work shows that the charging step proceeds through a multiple proton transfer mechanism where the amino group formed after deprotonation acts as the base to capture a proton from O3' of A76 rather than the phosphate group. The current study also shows the important role played by Asp233 in the proton transfer process.

Computational Methods for Molecular Understanding of the Antibiotic-Aminoacyl tRNA Synthetase Interaction

Developing an understanding of the interactions between an antibiotic and its binding site in a pathogen cell is the key to antibiotic design-an important cost-saving methodology compared to the costly and time-consuming random trial-and-error approach. The rapid development of antibiotic resistance provides an impetus for such studies. Recent years have witnessed the beginning of the use of combined computational techniques, including computer simulations and quantum mechanical computations, to understand how antibiotics bind at the active site of aminoacyl tRNA synthetases (aaRSs) from pathogens. Such computational protocols assist the knowledge-based design of antibiotics targeting aaRSs, which are their validated targets. After the ideas behind the protocols and their strategic planning are discussed, the protocols are described along with their major outcomes. This is followed by an integration of results from the different basic protocols. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Analysis of active-site residues from primary sequence of synthetase and transfer RNAs Basic Protocol 2: Molecular dynamics simulation-based protocol to study the structure and dynamics of the aaRS active site:antibiotic complex Basic Protocol 3: Quantum mechanical method-based protocol to study the structure and dynamics of the aaRS active site:antibiotic complex.

Free Energy Landscape of the Adenylation Reaction of the Aminoacylation Process at the Active Site of Aspartyl tRNA Synthetase

The process of protein biosynthesis is initiated by the aminoacylation process where a transfer ribonucleic acid (tRNA) is charged by the attachment of its cognate amino acid at the active site of the corresponding aminoacyl tRNA synthetase enzyme. The first step of the aminoacylation process, known as the adenylation reaction, involves activation of the cognate amino acid where it reacts with a molecule of adenosine triphosphate (ATP) at the active site of the enzyme to form the aminoacyl adenylate and inorganic pyrophosphate. In the current work, we have investigated the adenylation reaction between aspartic acid and ATP at the active site of the fully solvated aspartyl tRNA synthetase (AspRS) from Escherichia coli in aqueous medium at room temperature through hybrid quantum mechanical/molecular mechanical (QM/MM) simulations combined with enhanced sampling methods of well-tempered and well-sliced metadynamics. The objective of the present work is to study the associated free energy landscape and reaction barrier and also to explore the effects of active site mutation on the free energy surface of the reaction. The current calculations include finite temperature effects on free energy profiles. In particular, apart from contributions of interaction energies, the current calculations also include contributions of conformational, vibrational, and translational entropy of active site residues, substrates, and also the rest of the solvated protein and surrounding water into the free energy calculations. The present QM/MM metadynamics simulations predict a free energy barrier of 23.35 and 23.5 kcal/mol for two different metadynamics methods used to perform the reaction at the active site of the wild type enzyme. The free energy barrier increases to 30.6 kcal/mol when Arg217, which is an important conserved residue of the wild type enzyme at its active site, is mutated by alanine. These free energy results including the effect of mutation compare reasonably well with those of kinetic experiments that are available in the literature. The current work also provides molecular details of structural changes of the reactants and surroundings as the system dynamically evolves along the reaction pathway from reactant to the product state through QM/MM metadynamics simulations.

Dynamics of the Catalytic Active Site of Isoleucyl tRNA Synthetase from Staphylococcus aureus bound with Adenylate and Mupirocin

The development of new antimicrobial drugs is critically needed due to the alarming increase in antibiotic resistance in bacterial pathogens. The active sites of different bacterial aminoacyl tRNA synthetases (aaRS) are validated targets of antibiotics. At present, the only aaRS inhibitor approved is mupirocin (MRC) which targets bacterial isoleucyl tRNA synthetase (IleRS). The present work is aimed at understanding the lacunae of knowledge concerning the active site conformational dynamics in IleRS in the presence of inhibitor mupirocin. With this end in view, we have carried out classical molecular dynamics simulation and metadynamics simulations of the open state of IleRS from Staphylococcus aureus (^SaIleRS), the closed state tripartite complex bound with cognate adenylate (Ile-AMP) and tRNA, the closed state tripartite complex bound with noncognate MRC, and the closed state tripartite complex bound with tRNA and MRC with mutated ^SaIleRS (V588F). The present simulation established a dynamic picture of ^SaIleRS complexed with cognate and the noncognate substrates which are completely consistent with crystallographic and biochemical studies and explain the existing lacunae of knowledge. The active site is significantly more compact in the Ile-AMP bound complex compared to the open state due to the closure of the KMSKS and HMGH loops and clamping down of the tRNA acceptor end near the active site gate. The present result shows that the unusual open conformational state of the KMSKS loop observed in the cognate substrate-bound complex in the crystal is due to crystallographic constraints. Although the mupirocin tightly fits the catalytic active site in the MRC-bound complex, the nonanoic acid moiety is partly exposed to the water. The KMSKS loop is pushed open in the MRC-bound complex to accommodate the noncognate MRC. New tunnels open up, extending to the editing site in the complex. Out of its three broad segments, the C12 to C17 segment, the conjugated segment, and the nonanoic moiety, the conjugated part of MRC binds most effectively with the active site of the MRC-bound complex. The aromatic residues packing around the C12 to C17 segment of MRC stabilize the tRNA hairpin conformation in a similar way as observed in the TrpRS. The V588F mutation is weakening the interaction between this region of the active site and weakens the binding of MRC in the active site. This result explains why the V588F mutation is responsible for low-level mupirocin resistance. The free energy of unbinding the conjugated segment (and C12 to C17 segment, as well) largely contributes to the total free energy of unbinding the MRC. The active site organization of IleRS from eukaryotic Candida albicans is compared with the bacterial IleRS active site to understand the low binding affinity in eukaryotic IleRS. The present study could be a starting point of future studies related to the development of effective drug binding in the ^SaIleRS.

Role of the transfer ribonucleic acid (tRNA) bound magnesium ions in the charging step of aminoacylation reaction in the glutamyl tRNA synthetase and the seryl tRNA synthetase bound with cognate tRNA

Aminoacylation reaction is the first step of protein biosynthesis. Transfer RNA (tRNA) is charged with an amino acid in this reaction and the reaction is catalyzed by aminoacyl tRNA synthetase enzyme (aaRS). In the present work, we use classical molecular dynamics simulation to show that the tRNA bound Mg²⁺ ions significantly influence the charging step of class I ^TtGluRS: Glu-AMP: tRNA^Glu and class II dimeric ^TtSerRS: Ser-AMP: tRNA^Ser. The CCA end of the acceptor terminal is disordered in the absence of coordinated Mg²⁺ ions and the CCA end can freely explore beyond the specific conformational space of the tRNA in its precharging state. A balance between the conformational disorder of the tRNA and the restriction imposed on the CCA terminal via coordination with the Mg²⁺ ions is needed for the placement of the CCA terminal in a precharging state organization. This result provides a molecular-level explanation of the experimental observation that the presence of Mg²⁺ ions is a necessary condition for a successful aminoacylation reaction.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulations of cognate and non-cognate AspRS-tRNAAsp complexes

Aspartyl tRNA synthetase (AspRS), one of the 20 aminoacyl-tRNA synthetases, plays an important role in protein synthesis by catalyzing the aminoacylation reaction and synthesises Aspartyl-tRNA (tRNA^Asp). A typical three-dimensional structure of AspRS comprises three distinct domains for the recognition of cognate tRNA and catalysis, namely, anti-codon binding domain/N-terminal domain, hinge domain and catalytic domain through their interactions with anti-codon loop, D-stem and acceptor arm of cognate tRNA, respectively. In this work, we have studied the structural characteristics of each domain of AspRS to understand the recognition mechanism of tRNA^Asp using molecular dynamics simulations. The dynamics of AspRS-tRNA^Asp complexes from E.coli (cognate and non-cognate), S.cerevisiae (cognate) and T.thermophilus (non-cognate) were compared to understand the differences in recognition of cognate and non-cognate tRNAs. Our results explain that the conformational changes associated with the recognition of tRNA occur only in the cognate complexes. Among the cognate complexes, the conformational changes in yeast AspRS are highly controlled during tRNA^Asp recognition than that of in the E. coli AspRS. Moreover, the functional motions required for the tRNA recognition are observed only in the cognate complexes, and the conformational changes in AspRS and their recognition of tRNA^Asp are organism specific.Communicated by Ramaswamy H. Sarma.

Dynamics of the active site loops in catalyzing aminoacylation reaction in seryl and histidyl tRNA synthetases

Aminoacylation reaction is the first step of protein biosynthesis. The catalytic reorganization at the active site of aminoacyl tRNA synthetases (aaRSs) is driven by the loop motions. There remain lacunae of understanding concerning the catalytic loop dynamics in aaRSs. We analyzed the functional loop dynamics in seryl tRNA synthetase from Methanopyrus kandleri (^mkSerRS) and histidyl tRNA synthetases from Thermus thermophilus (^ttHisRS), respectively, using molecular dynamics. Results confirm that the motif 2 loop and other active site loops are flexible spots within the catalytic domain. Catalytic residues of the loops form a network of interaction with the substrates to form a reactive state. The loops undergo transitions between closed state and open state and the relaxation of the constituent residues occurs in femtosecond to nanosecond time scale. Order parameters are higher for constituent catalytic residues which form a specific network of interaction with the substrates to form a reactive state compared to the Gly residues within the loop. The development of interaction is supported from mutation studies where the catalytic domain with mutated loop exhibits unfavorable binding energy with the substrates. During the open-close motion of the loops, the catalytic residues make relaxation by ultrafast librational motion as well as fast diffusive motion and subsequently relax rather slowly via slower diffusive motion. The Gly residues act as a hinge to facilitate the loop closing and opening by their faster relaxation behavior. The role of bound water is analyzed by comparing implicit solvent-based and explicit solvent-based simulations. Loops fail to form catalytically competent geometry in absence of water. The present result, for the first time reveals the nature of the active site loop dynamics in aaRS and their influence on catalysis.

MD Simulations of tRNA and Aminoacyl-tRNA Synthetases: Dynamics, Folding, Binding, and Allostery

While tRNA and aminoacyl-tRNA synthetases are classes of biomolecules that have been extensively studied for decades, the finer details of how they carry out their fundamental biological functions in protein synthesis remain a challenge. Recent molecular dynamics (MD) simulations are verifying experimental observations and providing new insight that cannot be addressed from experiments alone. Throughout the review, we briefly discuss important historical events to provide a context for how far the field has progressed over the past few decades. We then review the background of tRNA molecules, aminoacyl-tRNA synthetases, and current state of the art MD simulation techniques for those who may be unfamiliar with any of those fields. Recent MD simulations of tRNA dynamics and folding and of aminoacyl-tRNA synthetase dynamics and mechanistic characterizations are discussed. We highlight the recent successes and discuss how important questions can be addressed using current MD simulations techniques. We also outline several natural next steps for computational studies of AARS:tRNA complexes.