Molecular Dynamics Investigation into Substrate Binding and Identity of the Catalytic Base in the Mechanism of Threonyl-tRNA Synthetase

QM/MM investigation of the discriminatory pre-transfer editing mechanism operated by Lysyl-tRNA synthetase

Aminoacyl-tRNA synthetases (aaRSs) are essential enzymes that remarkable facilitate the aminoacylation process during translation. With a high fidelity, the mischarged tRNA is prevented through implementing pre- and post-transfer proofreading mechanisms. For instance, Lysine-tRNA synthetase charges the native substrate, lysine, to its cognate tRNA. In spite of the great structural similarity between lysine to the noncognate and toxic ornithine, with the side chain of lysine being only one methylene group longer, LysRS is able to achieve this discrimination with a high efficiency. In this work, the hybrid quantum mechanics/molecular mechanics (QM/MM) investigation was applied to probe the pre-transfer editing mechanism catalyzed by lysyl-tRNA synthetase to reject the noncognte aminoacyl, L-ornityl (Orn), compared to the cognate substrate, L-lysyl. Particularly, the self-cyclization pre-transfer editing mechanism was explored for the two substrates. The substrate-assisted self-cyclization editing of Orn-AMP, where its phosphate moiety acts as the catalytic base, is found to be the rate-determining step with an energy barrier of 101.2 kJ mol-1. Meanwhile, the corresponding rate-limiting pathway for the native Lys-AMP lies at 140.2 kJ mol-1. This observation clearly indicated the infeasibility of this catalytic scenario in the presence of the native substrate. Interestingly, a thermodynamically favorable cyclic product of -92.9 kJ mol-1 with respect to the aminoacyl reactant complex demonstrated evidence of a successful pre-transfer editing. This reaction resulted in the discharge of the on-cognate -ornithine derivative from LysU's active site. These valuable mechanistic insights are valuable to enrich our knowledge of this extremely efficient and specific catalytic machinery of LysRS.Communicated by Ramaswamy H. Sarma.

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.

Fundamentals behind the specificity of Cysteinyl-tRNA synthetase: MD and QM/MM joint investigations

Cysteinyl-tRNA synthetase (CysRS) catalyzes the aminoacylation reaction of cysteine to its cognate tRNA^Cys in the first step of protein translation. It is found that CysRS is different from other aaRSs as it transfers cysteine without the need for an editing reaction, which is not applicable in the case of serine despite the similarity in their structures. Surprisingly, the reasons why CysRS has high amino acid specificity are not clear yet. In this research, the binding configurations of Cys-AMP and its near-cognate amino acid Ser-AMP with CysRS are compared by Molecular Dynamics (MD). The results reveal that CysRS screens the substrate Cys-AMP to a certain extent in the process of combination and recognition, thus providing a guarantee for the high selectivity of the next reaction. While Ser-AMP is in a folded state in CysRS. In the meanwhile, the interaction between Cys-AMP and Zn963 in CysRS is much stronger than Ser-AMP. The substrate-assisted aminoacylation mechanism in CysRS is also explored by Quantum Mechanics/Molecular Mechanics (QM/MM) modeling. According to the QM/MM potential energies, the energy barrier of TS_Cys-AMP is 91.75 kJ/mol, while that of TS_Ser-AMP is close to 150 kJ/mol. Based on thermochemistry calculations, it is found that the product of Cys-AMP is more stable than the reactant. In contrast, Ser-AMP has a reactant that is more stable than its product. As a result, it reflects that the specificity of CysRS originates from both the kinetic and thermodynamical perspectives of the reaction. Our investigations demonstrate comprehensively on how CysRS recognizes and catalyzes the substrate Cys-AMP, hoping to provide some guidance for researchers in this area.

Role of salt-bridging interactions in recognition of viral RNA by arginine-rich peptides

Interactions between RNA molecules and proteins are critical to many cellular processes and are implicated in various diseases. The RNA-peptide complexes are good model systems to probe the recognition mechanism of RNA by proteins. In this work, we report studies on the binding-unbinding process of a helical peptide from a viral RNA element using nonequilibrium molecular dynamics simulations. We explored the existence of various dissociation pathways with distinct free-energy profiles that reveal metastable states and distinct barriers to peptide dissociation. We also report the free-energy differences for each of the four pathways to be 96.47 ± 12.63, 96.1 ± 10.95, 91.83 ± 9.81, and 92 ± 11.32 kcal/mol. Based on the free-energy analysis, we further propose the preferred pathway and the mechanism of peptide dissociation. The preferred pathway is characterized by the formation of sequential hydrogen-bonding and salt-bridging interactions between several key arginine amino acids and the viral RNA nucleotides. Specifically, we identified one arginine amino acid (R8) of the peptide to play a significant role in the recognition mechanism of the peptide by the viral RNA molecule.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

A transition state “trapped”? QM-cluster models of engineered threonyl-tRNA synthetase

QM-cluster models demonstrate how protein bioengineering alters the local energy landscape of p-biphenylalanine to stabilize a transition state analogue.

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.

Insights from molecular dynamics on substrate binding and effects of active site mutations in Δ1-pyrroline-5-carboxylate dehydrogenase

The NAD+-dependent enzyme, Δ1-pyrroline-5-carboxylate dehydrogenase (P5CDH), has an important role in proline and hydroxyproline catabolism for humans. Specifically, this aldehyde dehydrogenase is responsible for the oxidation of both l-glutamate-γ-semialdehyde (GSA) and 4-erythro-hydroxy-l-glutamate-γ-semialdehyde (4-OH-GSA) to their respective l-glutamate product forms. We have performed a detailed molecular dynamics (MD) study of both the reactant and product complex structures of P5CDH to gain insights into ligand binding (i.e., GSA, 4-OH-GSA, NAD+, GLU) in the active site. Moreover, our investigations were further extended to examine the structural impact of S352L, S352A, and E314A mutations on the deficiency in the P5CDH enzymatic activity. Our in silico mutation analysis indicated that the conserved Glu447 has significantly shifted in both the S352L and E314A mutants, causing NAD+ to be displaced from its predictive orientation in the binding site and hence forming a catalytically inactive enzyme. However in the case of S352A, the catalytic site including the oxyanion hole and Cys348 remain virtually unchanged, and the coenzyme maintains its binding position.

Aminoacyl-tRNA Synthetases in the Bacterial World

Aminoacyl-tRNA synthetases (aaRSs) are modular enzymes globally conserved in the three kingdoms of life. All catalyze the same two-step reaction, i.e., the attachment of a proteinogenic amino acid on their cognate tRNAs, thereby mediating the correct expression of the genetic code. In addition, some aaRSs acquired other functions beyond this key role in translation. Genomics and X-ray crystallography have revealed great structural diversity in aaRSs (e.g., in oligomery and modularity, in ranking into two distinct groups each subdivided in 3 subgroups, by additional domains appended on the catalytic modules). AaRSs show huge structural plasticity related to function and limited idiosyncrasies that are kingdom or even species specific (e.g., the presence in many Bacteria of non discriminating aaRSs compensating for the absence of one or two specific aaRSs, notably AsnRS and/or GlnRS). Diversity, as well, occurs in the mechanisms of aaRS gene regulation that are not conserved in evolution, notably between distant groups such as Gram-positive and Gram-negative Bacteria. The review focuses on bacterial aaRSs (and their paralogs) and covers their structure, function, regulation, and evolution. Structure/function relationships are emphasized, notably the enzymology of tRNA aminoacylation and the editing mechanisms for correction of activation and charging errors. The huge amount of genomic and structural data that accumulated in last two decades is reviewed, showing how the field moved from essentially reductionist biology towards more global and integrated approaches. Likewise, the alternative functions of aaRSs and those of aaRS paralogs (e.g., during cell wall biogenesis and other metabolic processes in or outside protein synthesis) are reviewed. Since aaRS phylogenies present promiscuous bacterial, archaeal, and eukaryal features, similarities and differences in the properties of aaRSs from the three kingdoms of life are pinpointed throughout the review and distinctive characteristics of bacterium-like synthetases from organelles are outlined.

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is inactivated by S-sulfuration in vitro

Hydrogen sulfide (H2S) is produced enzymatically by cystathionine β-synthase (CBS) and cystathionine γ-lyase (CSE), as well as other enzymes in mammalian tissues. These discoveries have led to the crowning of H2S as yet another toxic gas that serves as a gasotransmitter like NO and CO. H2S is thought to exert its biological effects through its reaction with cysteine thiols in proteins, yielding sulfurated thiol (-SSH) derivatives. One of the first proteins shown to be modified by H2S was glyceraldehyde 3-phosphate dehydrogenase (GAPDH) [1] where the S-sulfuration of the active site cysteine (Cys 152) resulted in ~7-fold increase in the activity of the enzyme. In the present study we have attempted to reproduce this result with no success. GAPDH in its reduced, or hydrogen peroxide, or glutathione disulfide, or nitrosonium oxidized forms was reacted with sulfide or polysulfides. Sulfide had no effect on reduced GAPDH activity, while polysulfides inhibited GAPDH to ~42% of control. S-sulfuration of GAPDH occurred at Cys 247 after sulfide treatment, Cys 156 and Cys 247 after polysulfide treatment. No evidence of S-sulfuration at active site Cys 152 was discovered. Both sulfide and polysulfide was able to restore the activity of glutathione disulfide oxidized GAPDH, but not to control untreated levels. Treatment of glutathione disulfide oxidized GAPDH with polysulfide also produced S-sulfuration of Cys 156. Treatment of a C156S mutant of GAPDH with sulfide and polysulfide resulted in S-sulfuration of Cys 152, which also caused a decrease and not an increase in enzymatic activity. Computational chemistry shows S-sulfuration of Cys 156 may affect the position of catalytic Cys 152, raising its pKa by 0.5, which may affect the nucleophilicity of Cys 152. The current study raises significant questions about the reported ability of H2S to activate GAPDH by the sulfuration of its active site thiol, and indicates that polysulfide is a stronger protein S-sulfurating agent than sulfide.

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.

Computational investigations on the catalytic mechanism of maleate isomerase: the role of the active site cysteine residues

The maleate isomerase (MI) catalysed isomerization of maleate to fumarate has been investigated using a wide range of computational modelling techniques, including small model DFT calculations, QM-cluster approach, quantum mechanical/molecular mechanical approach (QM/MM in the ONIOM formalism) and molecular dynamics simulations. Several fundamental questions regarding the mechanism were answered in detail, such as the activation and stabilization of the catalytic Cys in a rather hydrophobic active site. The two previously proposed mechanisms were considered, where either enediolate or succinyl-Cys intermediate forms. Small model calculations as well as an ONIOM-based approach suggest that an enediolate intermediate is too unstable. Furthermore, the formation of succinyl-Cys intermediate via the nucleophilic attack of Cys76(-) on the substrate C2 (as proposed experimentally) was found to be energetically unfeasible in both QM-cluster and ONIOM approaches. Instead, our results show that Cys194, upon activation via the substrate, acts as a nucleophile and Cys76 acts as an acid/base catalyst, forming a succinyl-Cys intermediate in a concerted fashion. Indeed, the calculated PA of Cys76 is always higher than that of Cys194 before or upon substrate binding in the active site. Furthermore, the mechanism proceeds via multiple steps by substrate rotation around C2-C3 with the assistance of the now negatively charged Cys76, leading to the formation of fumarate. Finally, our calculated barrier is in good agreement with experiment. These findings represent a novel mechanism in the racemase superfamily.

A multi-scale computational study on the mechanism of Streptococcus pneumoniae Nicotinamidase (SpNic)

Nicotinamidase (Nic) is a key zinc-dependent enzyme in NAD metabolism that catalyzes the hydrolysis of nicotinamide to give nicotinic acid. A multi-scale computational approach has been used to investigate the catalytic mechanism, substrate binding and roles of active site residues of Nic from Streptococcus pneumoniae (SpNic). In particular, density functional theory (DFT), molecular dynamics (MD) and ONIOM quantum mechanics/molecular mechanics (QM/MM) methods have been employed. The overall mechanism occurs in two stages: (i) formation of a thioester enzyme-intermediate (IC2) and (ii) hydrolysis of the thioester bond to give the products. The polar protein environment has a significant effect in stabilizing reaction intermediates and in particular transition states. As a result, both stages effectively occur in one step with Stage 1, formation of IC2, being rate limiting barrier with a cost of 53.5 kJ·mol-1 with respect to the reactant complex, RC. The effects of dispersion interactions on the overall mechanism were also considered but were generally calculated to have less significant effects with the overall mechanism being unchanged. In addition, the active site lysyl (Lys103) is concluded to likely play a role in stabilizing the thiolate of Cys136 during the reaction.

Molecular recognition mechanism of peptide chain bound to the tRNA(Lys3) anticodon loop in silico

The mechanism by which proteins recognize and bind the post-transcriptional modifications of RNAs is unknown, yet these interactions play important functions in biology. Atomistic molecular dynamics simulations were performed to examine the folding of the model peptide chain -RVTHHAFLGAHRTVG- and the complex formed by the folded peptide with the native anticodon stem and loop of the human tRNA(Lys3) (hASL(Lys3)) in order to explore the binding mechanism. By analyzing and comparing two folded conformations of this peptide obtained from the folding simulation, we found that the van der Waals (VDW) energy is necessary for the thermal stability of the peptide, and the charge-charge (ELE + EGB) energy is crucial for determining the three-dimensional folded structure of the peptide backbone. Subsequently, two conformations of the peptide were employed to investigate their binding behaviors to hASL(Lys3). The metastable folded peptide was found to bind to hASL(Lys3) much easier than the stable folded peptide in the binding simulations. An energetic analysis reveals that the VDW energy favors the binding, whereas the ELE + EGB energies disfavor the binding. Arginines on the peptide preferentially attract the phosphate backbone via the inter-chain ELE + EGB interaction, significantly contributing to the binding affinity. The hydrophobic phenylalanine interacts with the anticodon loop of hASL(Lys3) via the inter-chain VDW interaction, significantly contributing to the binding specificity.

Multi-scale computational enzymology: enhancing our understanding of enzymatic catalysis

Elucidating the origin of enzymatic catalysis stands as one the great challenges of contemporary biochemistry and biophysics. The recent emergence of computational enzymology has enhanced our atomistic-level description of biocatalysis as well the kinetic and thermodynamic properties of their mechanisms. There exists a diversity of computational methods allowing the investigation of specific enzymatic properties. Small or large density functional theory models allow the comparison of a plethora of mechanistic reactive species and divergent catalytic pathways. Molecular docking can model different substrate conformations embedded within enzyme active sites and determine those with optimal binding affinities. Molecular dynamics simulations provide insights into the dynamics and roles of active site components as well as the interactions between substrate and enzymes. Hybrid quantum mechanical/molecular mechanical (QM/MM) can model reactions in active sites while considering steric and electrostatic contributions provided by the surrounding environment. Using previous studies done within our group, on OvoA, EgtB, ThrRS, LuxS and MsrA enzymatic systems, we will review how these methods can be used either independently or cooperatively to get insights into enzymatic catalysis.