Analysis of predictions for the catalytic mechanism of ribosomal peptidyl transfer

Reversing the relative time courses of the peptide bond reaction with oligopeptides of different lengths and charged amino acid distributions in the ribosome exit tunnel

The kinetics of the protein elongation cycle by the ribosome depends on intertwined factors. One of these factors is the electrostatic interaction of the nascent protein with the ribosome exit tunnel. In this computational biology theoretical study, we focus on the rate of the peptide bond formation and its dependence on the ribosome exit tunnel electrostatic potential profile. We quantitatively predict how oligopeptides of variable lengths can affect the peptide bond formation rate. We applied the Michaelis-Menten model as previously extended to incorporate the mechano-biochemical effects of forces on the rate of reaction at the catalytic site of the ribosome. For a given pair of carboxy-terminal amino acid substrate at the P- and an aminoacyl-tRNA at the A-sites, the relative time courses of the peptide bond formation reaction can be reversed depending on the oligopeptide sequence embedded in the tunnel and their variable lengths from the P-site. The reversal is predicted to occur from a shift in positions of charged amino acids upstream in the oligopeptidyl-tRNA at the P-site. The position shift must be adjusted by clever design of the oligopeptide probes using the electrostatic potential profile along the exit tunnel axial path. These predicted quantitative results bring strong evidence of the importance and relative contribution of the electrostatic interaction of the ribosome exit tunnel with the nascent peptide chain during elongation.

The energy landscape of the ribosome

The ribosome is a prototypical assembly that can be used to establish general principles and techniques for the study of biological molecular machines. Motivated by the fact that the dynamics of every biomolecule is governed by an underlying energy landscape, there has been great interest to understand and quantify ribosome energetics. In the present review, we will focus on theoretical and computational strategies for probing the interactions that shape the energy landscape of the ribosome, with an emphasis on more recent studies of the elongation cycle. These efforts include the application of quantum mechanical methods for describing chemical kinetics, as well as classical descriptions to characterize slower (microsecond to millisecond) large-scale (10-100 Å) rearrangements, where motion is described in terms of diffusion across an energy landscape. Together, these studies provide broad insights into the factors that control a diverse range of dynamical processes in this assembly.

Dynamics and Function of sRNA/mRNAs Under the Scrutiny of Computational Simulation Methods

Molecular dynamics simulations have proved extremely useful in investigating the functioning of proteins with atomic-scale resolution. Many applications to the study of RNA also exist, and their number increases by the day. However, implementing MD simulations for RNA molecules in solution faces challenges that the MD practitioner must be aware of for the appropriate use of this tool. In this chapter, we present the fundamentals of MD simulations, in general, and the peculiarities of RNA simulations, in particular. We discuss the strengths and limitations of the technique and provide examples of its application to elucidate small RNA's performance.

Caracal: A Versatile Ring Polymer Molecular Dynamics Simulation Package

A new open-source program package named Caracal covering simulations of molecular systems with ring polymer molecular dynamics (RPMD) is presented. It combines a powerful RPMD implementation including chemical reaction rate calculations and biased periodic and nonperiodic samplings with a collection of easy to set up potential energy surface (PES) methodologies, thus delivering an all-inclusive approach. Most implemented PESs are based on the QMDFF and EVB-QMDFF methods. Where the quantum mechanically derived force field (QMDFF) can be set up for an arbitrary molecular system in a black-box fashion, the empirical valence bond (EVB)-QMDFF connects two QMDFFs and is able to represent the PES of a chemical reaction. With our previously published flavors of this composite method, PESs for almost arbitrary gas phase thermal ground state reactions can be set up. Given an optimized reaction path, the mechanism of the reaction can be classified and RPMD rate constants can be obtained via umbrella sampling and recrossing calculations on an EVB-QMDFF PES. Further, QMDFFs can be polymerized for the description of liquid systems. In this paper, the internal structure as well as the handling philosophy of Caracal are outlined. Further, examples of the different possible kinds of calculations are given.

The Peptidyl Transferase Center: a Window to the Past

In his 2001 article, "Translation: in retrospect and prospect," the late Carl Woese made a prescient observation that there was a need for the then-current view of translation to be "reformulated to become an all-embracing perspective about which 21st century Biology can develop" (RNA 7:1055-1067, 2001, https://doi.org/10.1017/s1355838201010615). The quest to decipher the origins of life and the road to the genetic code are both inextricably linked with the history of the ribosome. After over 60 years of research, significant progress in our understanding of how ribosomes work has been made. Particularly attractive is a model in which the ribosome may facilitate an ∼180° rotation of the CCA end of the tRNA from the A-site to the P-site while the acceptor stem of the tRNA would then undergo a translation from the A-site to the P-site. However, the central question of how the ribosome originated remains unresolved. Along the path from a primitive RNA world or an RNA-peptide world to a proto-ribosome world, the advent of the peptidyl transferase activity would have been a seminal event. This functionality is now housed within a local region of the large-subunit (LSU) rRNA, namely, the peptidyl transferase center (PTC). The PTC is responsible for peptide bond formation during protein synthesis and is usually considered to be the oldest part of the modern ribosome. What is frequently overlooked is that by examining the origins of the PTC itself, one is likely going back even further in time. In this regard, it has been proposed that the modern PTC originated from the association of two smaller RNAs that were once independent and now comprise a pseudosymmetric region in the modern PTC. Could such an association have survived? Recent studies have shown that the extant PTC is largely depleted of ribosomal protein interactions. It is other elements like metallic ion coordination and nonstandard base/base interactions that would have had to stabilize the association of RNAs. Here, we present a detailed review of the literature focused on the nature of the extant PTC and its proposed ancestor, the proto-ribosome.

Structural and functional studies revealed key mechanisms underlying elongation step of protein translation

The ribosome is an ancient and universally conserved macromolecular machine that synthesizes proteins in all organisms. Since the discovery of the ribosome by electron microscopy in the mid-1950s, rapid progress has been made in research on it, regarding its architecture and functions. As a machine that synthesizes polypeptides, the sequential addition of amino acids to a growing polypeptide chain occurs during a phase called the elongation cycle. This is the core step of protein translation and is highly conserved between bacteria and eukarya. The elongation cycle involves codon recognition by aminoacyl tRNAs, catalysis of peptide bond formation, and the most complex operation of translation-translocation. In this review, we discuss the fundamental results from structural and functional studies over the past decades that have led to understanding of the three key questions underlying translation.

Mechanistic alternatives for peptide bond formation on the ribosome

The peptidyl transfer reaction on the large ribosomal subunit depends on the protonation state of the amine nucleophile and exhibits a large kinetic solvent isotope effect (KSIE ∼8). In contrast, the related peptidyl-tRNA hydrolysis reaction involved in termination shows a KSIE of ∼4 and a pH-rate profile indicative of base catalysis. It is, however, unclear why these reactions should proceed with different mechanisms, as the experimental data suggests. One explanation is that two competing mechanisms may be operational in the peptidyl transferase center (PTC). Herein, we explored this possibility by re-examining the previously proposed proton shuttle mechanism and testing the feasibility of general base catalysis also for peptide bond formation. We employed a large cluster model of the active site and different reaction mechanisms were evaluated by density functional theory calculations. In these calculations, the proton shuttle and general base mechanisms both yield activation energies comparable to the experimental values. However, only the proton shuttle mechanism is found to be consistent with the experimentally observed pH-rate profile and the KSIE. This suggests that the PTC promotes the proton shuttle mechanism for peptide bond formation, while prohibiting general base catalysis, although the detailed mechanism by which general base catalysis is excluded remains unclear.

Substrate-Induced Formation of Ribosomal Decoding Center for Accurate and Rapid Genetic Code Translation

Accurate translation of genetic information is crucial for synthesis of functional proteins in all organisms. We use recent experimental data to discuss how induced fit affects accuracy of initial codon selection on the ribosome by aminoacyl transfer RNA in ternary complex ( T₃) with elongation factor Tu (EF-Tu) and guanosine-5'-triphosphate (GTP). We define actual accuracy ([Formula: see text]) of a particular protein synthesis system as its current accuracy and the effective selectivity ([Formula: see text]) as [Formula: see text] in the limit of zero ribosomal binding affinity for T₃. Intrinsic selectivity ([Formula: see text]), defined as the upper thermodynamic limit of [Formula: see text], is determined by the free energy difference between near-cognate and cognate T₃ in the pre-GTP hydrolysis state on the ribosome. [Formula: see text] is much larger than [Formula: see text], suggesting the possibility of a considerable increase in [Formula: see text] and [Formula: see text] at negligible kinetic cost. Induced fit increases [Formula: see text] and [Formula: see text] without affecting [Formula: see text], and aminoglycoside antibiotics reduce [Formula: see text] and [Formula: see text] at unaltered [Formula: see text].

Mechanistic insights into the slow peptide bond formation with D-amino acids in the ribosomal active site

During protein synthesis, ribosomes discriminate chirality of amino acids and prevent incorporation of D-amino acids into nascent proteins by slowing down the rate of peptide bond formation. Despite this phenomenon being known for nearly forty years, no structures have ever been reported that would explain the poor reactivity of D-amino acids. Here we report a 3.7Å-resolution crystal structure of a bacterial ribosome in complex with a D-aminoacyl-tRNA analog bound to the A site. Although at this resolution we could not observe individual chemical groups, we could unambiguously define the positions of the D-amino acid side chain and the amino group based on chemical restraints. The structure reveals that similarly to L-amino acids, the D-amino acid binds the ribosome by inserting its side chain into the ribosomal A-site cleft. This binding mode does not allow optimal nucleophilic attack of the peptidyl-tRNA by the reactive α-amino group of a D-amino acid. Also, our structure suggests that the D-amino acid cannot participate in hydrogen-bonding with the P-site tRNA that is required for the efficient proton transfer during peptide bond formation. Overall, our work provides the first mechanistic insight into the ancient mechanism that helps living cells ensure the stereochemistry of protein synthesis.

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.

Molecular simulations of the ribosome and associated translation factors

The ribosome is a macromolecular complex which is responsible for protein synthesis in all living cells according to their transcribed genetic information. Using X-ray crystallography and, more recently, cryo-electron microscopy (cryo-EM), the structure of the ribosome was resolved at atomic resolution in many functional and conformational states. Molecular dynamics simulations have added information on dynamics and energetics to the available structural information, thereby have bridged the gap to the kinetics obtained from single-molecule and bulk experiments. Here, we review recent computational studies that brought notable insights into ribosomal structure and function.

Mechanistic Insights Into Catalytic RNA-Protein Complexes Involved in Translation of the Genetic Code

The contemporary world is an "RNA-protein world" rather than a "protein world" and tracing its evolutionary origins is of great interest and importance. The different RNAs that function in close collaboration with proteins are involved in several key physiological processes, including catalysis. Ribosome-the complex megadalton cellular machinery that translates genetic information encoded in nucleotide sequence to amino acid sequence-epitomizes such an association between RNA and protein. RNAs that can catalyze biochemical reactions are known as ribozymes. They usually employ general acid-base catalytic mechanism, often involving the 2'-OH of RNA that activates and/or stabilizes a nucleophile during the reaction pathway. The protein component of such RNA-protein complexes (RNPCs) mostly serves as a scaffold which provides an environment conducive for the RNA to function, or as a mediator for other interacting partners. In this review, we describe those RNPCs that are involved at different stages of protein biosynthesis and in which RNA performs the catalytic function; the focus of the account is on highlighting mechanistic aspects of these complexes. We also provide a perspective on such associations in the context of proofreading during translation of the genetic code. The latter aspect is not much appreciated and recent works suggest that this is an avenue worth exploring, since an understanding of the subject can provide useful insights into how RNAs collaborate with proteins to ensure fidelity during these essential cellular processes. It may also aid in comprehending evolutionary aspects of such associations.

Induced fit of the peptidyl-transferase center of the ribosome and conformational freedom of the esterified amino acids

The catalytic site of most enzymes can efficiently handle only one substrate. In contrast, the ribosome is capable of polymerizing at a similar rate at least 20 different kinds of amino acids from aminoacyl-tRNA carriers while using just one catalytic site, the peptidyl-transferase center (PTC). An induced-fit mechanism has been uncovered in the PTC, but a possible connection between this mechanism and the uniform handling of the substrates has not been investigated. We present an analysis of published ribosome structures supporting the hypothesis that the induced fit eliminates unreactive rotamers predominantly populated for some A-site aminoacyl esters before induction. We show that this hypothesis is fully consistent with the wealth of kinetic data obtained with these substrates. Our analysis reveals that induction constrains the amino acids into a reactive conformation in a side-chain independent manner. It allows us to highlight the rationale of the PTC structural organization, which confers to the ribosome the very unusual ability to handle large as well as small substrates.

CADEE: Computer-Aided Directed Evolution of Enzymes

The tremendous interest in enzymes as biocatalysts has led to extensive work in enzyme engineering, as well as associated methodology development. Here, a new framework for computer-aided directed evolution of enzymes (CADEE) is presented which allows a drastic reduction in the time necessary to prepare and analyze in silico semi-automated directed evolution of enzymes. A pedagogical example of the application of CADEE to a real biological system is also presented in order to illustrate the CADEE workflow.

The role of caffeine as an inhibitor in the aggregation of amyloid forming peptides: a unified molecular dynamics simulation and experimental study

10 : 1 or more caffeine to SwP ratio can prevent the formation of β-sheet conformation.

Exceptionally large entropy contributions enable the high rates of GTP hydrolysis on the ribosome

Protein synthesis on the ribosome involves hydrolysis of GTP in several key steps of the mRNA translation cycle. These steps are catalyzed by the translational GTPases of which elongation factor Tu (EF-Tu) is the fastest GTPase known. Here, we use extensive computer simulations to explore the origin of its remarkably high catalytic rate on the ribosome and show that it is made possible by a very large positive activation entropy. This entropy term (TΔS(‡)) amounts to more than 7 kcal/mol at 25 °C. It is further found to be characteristic of the reaction mechanism utilized by the translational, but not other, GTPases and it enables these enzymes to attain hydrolysis rates exceeding 500 s(-1). This entropy driven mechanism likely reflects the very high selection pressure on the speed of protein synthesis, which drives the rate of each individual GTPase towards maximal turnover rate of the whole translation cycle.

The P-site A76 2′-OH acts as a peptidyl shuttle in a stepwise peptidyl transfer mechanism

The P-site-A76-2′OH transfers the polypeptide chain to the A-site α-amine and A2451 facilitates this transfer by acting as proton shuttle.

Heavy water reduces GFP expression in prokaryotic cell-free assays at the translation level while stimulating its transcription

The in vitro proliferation of prokaryotic and eukaryotic cells is remarkably hampered in the presence of heavy water (D2O). Impairment of gene expression at the transcription or translation level can be the base for this effect. However, insights into the underlying mechanisms are lacking. Here, we employ a cell-free expression system for the quantitative analysis of the effect of increasing percentages of D2O on the kinetics of in-vitro GFP expression. Experiments are designed to discriminate the rates of transcription, translation, and protein folding using pDNA and mRNA vectors, respectively. We find that D2O significantly stimulates GFP expression at the transcription level but acts as a suppressor at translation and maturation (folding) in a linear dose-dependent manner. At a D2O concentration of 60%, the GFP expression rate was reduced to 40% of an undisturbed sample. We observed a similar inhibition of GFP expression by D2O in a recombinant Escherichia coli strain, although the inhibitory effect is less pronounced. These results demonstrate the suitability of cell-free systems for quantifying the impact of heavy water on gene expression and establish a platform to further assess the potential therapeutic use of heavy water as antiproliferative agent.

The ribosome as a versatile catalyst: reactions at the peptidyl transferase center

In all contemporary organisms, the active site of the ribosome--the peptidyl transferase center--catalyzes two distinct reactions, peptide bond formation between peptidyl-tRNA and aminoacyl-tRNA as well as the hydrolysis of peptidyl-tRNA with the help of a release factor. However, when provided with appropriate substrates, ribosomes can also catalyze a broad range of other chemical reaction, which provides the basis for orthogonal translation and synthesis of alloproteins from unnatural building blocks. Advances in understanding the mechanisms of the two ubiquitous reactions, the peptide bond formation and peptide release, provide insights into the versatility of the active site of the ribosome. Release factors 1 and 2 and elongation factor P are auxiliary factors that augment the intrinsic catalytic activity of the ribosome in special cases.

ONIOM and ab-initio calculations on the mechanism of uncatalyzed peptide bond formation

Finding a proper transition structure for the peptide bond formation process can lead one to a better understanding of the role of ribosome in catalyzing this reaction. Using computer simulations, we performed the potential energy surface scan on the ester bond dissociation of P-site aminoacyl-tRNA and the peptide bond formation of P-site and A-site amino acids. The full fragments of initiator tRNA(i)(met) and elongator tRNA(phe) are attached to both cognate and non-cognate amino acids as the P-site substrate. The A-site amino acid for all four calculations is methionine. We used ONIOM calculations to reduce the computational cost. Our study illustrates the reduced rate of peptide bond formation for misacylated tRNA(i)(met) in the absence of ribosomal bases. The misacylated elongator tRNA(phe), however, did not show any difference in its PES compared with that for the phe-tRNA(phe). This demonstrates the structural specification of initiator tRNA(i)(met) for the amino acids side chain.

Bridging the gap between ribosome structure and biochemistry by mechanistic computations

The wealth of structural and biochemical data now available for protein synthesis on the ribosome presents major new challenges for computational biochemistry. Apart from technical difficulties in modeling ribosome systems, the complexity of the overall translation cycle with a multitude of different kinetic steps presents a formidable problem for computational efforts where we have only seen the beginning. However, a range of methodologies including molecular dynamics simulations, free energy calculations, molecular docking and quantum chemical approaches have already been put to work with promising results. In particular, the combined efforts of structural biology, biochemistry, kinetics and computational modeling can lead towards a quantitative structure-based description of translation.

Analysis of enzymatic transacylase Brønsted studies with application to the ribosome

Preferential binding of an enzyme to the transition state relative to the ground state is a key strategy for enzyme catalysis. When there is a difference between the ground and transition state charge distributions, enzymes maximize electrostatic interactions to achieve this enhanced transition state binding. Although the transition state is difficult to observe directly by structural methods, the chemical details of this transient species can be characterized by studies of substituent effects (Brønsted, Hammett, Swain-Scott, etc.) and isotope effects. Brønsted analysis can provide an estimate of transition state charges for the nucleophile and leaving group of a reaction. This Account will discuss the theoretical basis of Brønsted analysis and describe its practical application to the study of transacylase enzyme systems including the peptidyl transferase reaction of the ribosome. The Brønsted coefficient is derived from the linear free energy relationship (LFER) that correlates the acidity (pK(a)) of a reactive atom to the log of its rate constant. The Brønsted coefficient establishes the change in atomic charge as the reaction proceeds from the ground state to the transition state. Bonding events alter the electrostatics of atoms and the extent of bonding can be extrapolated from transition state charges. Therefore, well-defined nucleophile and leaving group transition state charges limit the number of mechanisms that are consistent with a particular transition state. Brønsted results are most informative when interpreted in the context of other mechanistic data, especially for enzymatic studies where an active site may promote a transition state that differs significantly from a prediction based on uncatalyzed solution reactions. Here we review Brønsted analyses performed on transacylases to illustrate how these data enhanced the enzymatic mechanistic studies. Through a systematic comparison of five enzymes, we reveal a wide spectrum of Brønsted values that are possible for what otherwise appear to be similar chemical reactions. The variations in the Brønsted coefficients predict different transition states for the various enzymes. This Account explores an overriding theme in the enzymatic mechanisms that catalysis enhances commensurate bond formation and proton abstraction events. The extent of the two bonding events in relationship to each other can be inferred from the Brønsted coefficient. When viewed in the context of recent ribosomal studies, this interpretation provides mechanistic insights into peptide bond formation.

Different substrate-dependent transition states in the active site of the ribosome

The active site of the ribosome, the peptidyl transferase centre, catalyses two reactions, namely, peptide bond formation between peptidyl-tRNA and aminoacyl-tRNA as well as the release-factor-dependent hydrolysis of peptidyl-tRNA. Unlike peptide bond formation, peptide release is strongly impaired by mutations of nucleotides within the active site, in particular by base exchanges at position A2602 (refs 1, 2). The 2'-OH group of A76 of the peptidyl-tRNA substrate seems to have a key role in peptide release. According to computational analysis, the 2'-OH may take part in a concerted 'proton shuttle' by which the leaving group is protonated, in analogy to similar current models of peptide bond formation. Here we report kinetic solvent isotope effects and proton inventories (reaction rates measured in buffers with increasing content of deuterated water, D(2)O) of the two reactions catalysed by the active site of the Escherichia coli ribosome. The transition state of the release factor 2 (RF2)-dependent hydrolysis reaction is characterized by the rate-limiting formation of a single strong hydrogen bond. This finding argues against a concerted proton shuttle in the transition state of the hydrolysis reaction. In comparison, the proton inventory for peptide bond formation indicates the rate-limiting formation of three hydrogen bonds with about equal contributions, consistent with a concerted eight-membered proton shuttle in the transition state. Thus, the ribosome supports different rate-limiting transition states for the two reactions that take place in the peptidyl transferase centre.

pH-sensitivity of the ribosomal peptidyl transfer reaction dependent on the identity of the A-site aminoacyl-tRNA

We studied the pH-dependence of ribosome catalyzed peptidyl transfer from fMet-tRNA(fMet) to the aa-tRNAs Phe-tRNA(Phe), Ala-tRNA(Ala), Gly-tRNA(Gly), Pro-tRNA(Pro), Asn-tRNA(Asn), and Ile-tRNA(Ile), selected to cover a large range of intrinsic pK(a)-values for the α-amino group of their amino acids. The peptidyl transfer rates were different at pH 7.5 and displayed different pH-dependence, quantified as the pH-value, pK(a)(obs), at which the rate was half maximal. The pK(a)(obs)-values were downshifted relative to the intrinsic pK(a)-value of aa-tRNAs in bulk solution. Gly-tRNA(Gly) had the smallest downshift, while Ile-tRNA(Ile) and Ala-tRNA(Ala) had the largest downshifts. These downshifts correlate strongly with molecular dynamics (MD) estimates of the downshifts in pK(a)-values of these aa-tRNAs upon A-site binding. Our data show the chemistry of peptide bond formation to be rate limiting for peptidyl transfer at pH 7.5 in the Gly and Pro cases and indicate rate limiting chemistry for all six aa-tRNAs.

Editing mechanism of aminoacyl-tRNA synthetases operates by a hybrid ribozyme/protein catalyst

Aminoacyl-tRNA synthetases (aaRSs) are critical for the translational process, catalyzing the attachment of specific amino acids to their cognate tRNAs. To ensure formation of the correct aminoacyl-tRNA, and thereby enhance the reliability of translation, several aaRSs have an editing capability that hinders formation of misaminoacylated tRNAs. We investigated theoretically the mechanism of the editing reaction for a class I enzyme, leucyl-tRNA synthetase (LeuRS), complexed with a misaminoacylated tRNA(Leu), employing ab initio hybrid quantum mechanical/molecular mechanical potentials in conjunction with molecular dynamics simulations. It is shown that the water molecule that acts as the nucleophile in the editing reaction is activated by a 3'-hydroxyl group at the 3'-end of tRNA(Leu) and that the O2' atom of the leaving group of the substrate is capped by one of the water's hydrogen atoms. Thus, it is shown that editing is a self-cleavage reaction of the tRNA and so it is the tRNA, and not the protein, that drives the reaction. The protein does, however, have an important stabilizing effect on some high-energy intermediates along the reaction path, which is more efficient than the ribozyme would be alone. This indicates that editing is achieved by a novel "hybrid ribozyme/protein catalyst". Analysis of existing experimental data and additional modeling shows that this ribozymal mechanism appears to be widespread, occurring in the ribosome as well as in other aaRSs. It also suggests transitional forms that could have played an important role in the RNA world hypothesis for the origin of life.

Transition states of uncatalyzed hydrolysis and aminolysis reactions of a ribosomal P-site substrate determined by kinetic isotope effects

The ester bond of peptidyl-tRNA undergoes nucleophilic attack in solution and when catalyzed by the ribosome. To characterize the uncatalyzed hydrolysis reaction, a model of peptide release, the transition state structure for hydrolysis of a peptidyl-tRNA mimic was determined. Kinetic isotope effects were measured at several atoms that potentially undergo a change in bonding in the transition state. Large kinetic isotope effects of carbonyl (18)O and alpha-deuterium substitutions on uncatalyzed hydrolysis indicate the transition state is nearly tetrahedral. Kinetic isotope effects were also measured for aminolysis by hydroxylamine to study a reaction similar to the formation of a peptide bond. In contrast to hydrolysis, the large leaving group (18)O isotope effect indicates the C-O3' bond has undergone significant scission in the transition state. The smaller carbonyl (18)O and alpha-deuterium effects are consistent with a later transition state. The assay developed here can also be used to measure isotope effects for the ribosome-catalyzed reactions. These uncatalyzed reactions serve as a basis for determining what aspects of the transition states are stabilized by the ribosome to achieve a rate enhancement.

The transition state for peptide bond formation reveals the ribosome as a water trap

Recent progress in elucidating the peptide bond formation process on the ribosome has led to notion of a proton shuttle mechanism where the 2'-hydroxyl group of the P-site tRNA plays a key role in mediating proton transfer between the nucleophile and leaving group, whereas ribosomal groups do not actively participate in the reaction. Despite these advances, the detailed nature of the transition state for peptidyl transfer and the role of several trapped water molecules in the peptidyl transferase center remain major open questions. Here, we employ high-level quantum chemical ab initio calculations to locate and characterize global transition states for the reaction, described by a molecular model encompassing all the key elements of the reaction center. The calculated activation enthalpy as well as structures are in excellent agreement with experimental data and point to feasibility of an eight-membered "double proton shuttle" mechanism in which an auxiliary water molecule, observed both in computer simulations and crystal structures, actively participates. A second conserved water molecule is found to be of key importance for stabilizing developing negative charge on the substrate oxyanion and its presence is catalytically favorable both in terms of activation enthalpy and entropy. Transition states calculated both for six- and eight-membered mechanisms are invariably late and do not involve significant charge development on the attacking amino group. Predicted kinetic isotope effects consistent with this picture are similar to those observed for uncatalyzed ester aminolysis reactions in solution.

tRNA as an active chemical scaffold for diverse chemical transformations

During protein synthesis, tRNA serves as the intermediary between cognate amino acids and their corresponding RNA trinucleotide codons. Aminoacyl-tRNA is also a biosynthetic precursor and amino acid donor for other macromolecules. AA-tRNAs allow transformations of acidic amino acids into their amide-containing counterparts, and seryl-tRNA(Ser) donates serine for antibiotic synthesis. Aminoacyl-tRNA is also used to cross-link peptidoglycan, to lysinylate the lipid bilayer, and to allow proteolytic turnover via the N-end rule. These alternative functions may signal the use of RNA in early evolution as both a biological scaffold and a catalyst to achieve a wide variety of chemical transformations.

Theoretical studies of RNA catalysis: hybrid QM/MM methods and their comparison with MD and QM

Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome. We compare and correlate QM/MM results with those from QM and/or molecular dynamics (MD) simulations, and discuss scope and limitations with a critical eye on current shortcomings in available methodologies and computer resources. We thus hope to foster mutual appreciation and facilitate collaboration between experimentalists and theorists to jointly advance our understanding of RNA catalysis at an atomistic level.

Total syntheses of a conformationally locked North-type methanocarba puromycin analogue and a dinucleotide derivative

An original synthetic approach for the first synthesis of an enantiopure methanocarba puromycin (3'-alpha-aminoacylamino-3'-deoxyadenosine) analogue and its cytidine dinucleotide derivative is described. Each compound is conformationally locked in a North-type pucker and exhibits both a pseudoaxial hydroxy group and a pseudoequatorial aminoacyl group. The syntheses were accomplished from D-ribose in 18 and 19 steps, respectively, with key steps being a ring-closing metathesis, a Luche reduction, a Simmons-Smith cyclopropanation, a Mitsunobu coupling, a Mattocks bromoacetylation, a regioselective and a stereoselective nucleophilic substitution, a chemoselective phosphoramidite coupling and a Staudinger-Vilarrasa coupling. Both molecules are being tested for peptidyl transfer efficiency in ribosomes for comparison with the peptidyl transfer kinetics of natural puromycin and other natural and synthetic ribosomal A site substrates.