The tetrahymena ribozyme cleaves a 5'-methylene phosphonate monoester approximately 10(2)-fold faster than a normal phosphate diester: implications for enzyme catalysis of phosphoryl transfer reactions

Unveiling the Cleavage Mechanism of an RNA Model Compound on the whole pH Scale: Computations Meet Experiments in the Determination of Reaction Rates

The knowledge of the mechanism of reactions occurring in solution is a primary research line both in the context of theoretical-computational chemistry and in the field of organic and bio-organic chemistry. Given the importance of the hydrolysis of nucleic acids in life-related phenomena, here we present a combined experimental and computational study on the cleavage of an RNA model compound. This phosphodiester features a cleavage rate strictly dependent on the pH with three different dependence domains. Such experimental evidence, highlighted by an in-depth kinetic investigation, unequivocally suggests a change in the reaction mechanism along the pH scale. In order to interpret the data and to explain the experimental behavior, we have applied a theoretical-computational procedure, involving a hybrid quantum/classical approach, able to model chemical reactions in complex environments, i. e. in solution. This study turns out to quantitatively reproduce the experimental data with accuracy and, in addition, provides useful mechanistic insight into the transesterification process of the investigated compound. The study indicates that the cleavage can occur through anA N D N ${A_N D_N }$ , anA N + D N ${A_N + D_N }$ , and aD N A N ${D_N A_N }$ mechanism depending on the pH values.

A kinetic and thermodynamic framework for the Azoarcus group I ribozyme reaction

Determination of quantitative thermodynamic and kinetic frameworks for ribozymes derived from the Azoarcus group I intron and comparisons to their well-studied analogs from the Tetrahymena group I intron reveal similarities and differences between these RNAs. The guanosine (G) substrate binds to the Azoarcus and Tetrahymena ribozymes with similar equilibrium binding constants and similar very slow association rate constants. These and additional literature observations support a model in which the free ribozyme is not conformationally competent to bind G and in which the probability of assuming the binding-competent state is determined by tertiary interactions of peripheral elements. As proposed previously, the slow binding of guanosine may play a role in the specificity of group I intron self-splicing, and slow binding may be used analogously in other biological processes. The internal equilibrium between ribozyme-bound substrates and products is similar for these ribozymes, but the Azoarcus ribozyme does not display the coupling in the binding of substrates that is observed with the Tetrahymena ribozyme, suggesting that local preorganization of the active site and rearrangements within the active site upon substrate binding are different for these ribozymes. Our results also confirm the much greater tertiary binding energy of the 5'-splice site analog with the Azoarcus ribozyme, binding energy that presumably compensates for the fewer base-pairing interactions to allow the 5'-exon intermediate in self splicing to remain bound subsequent to 5'-exon cleavage and prior to exon ligation. Most generally, these frameworks provide a foundation for design and interpretation of experiments investigating fundamental properties of these and other structured RNAs.

Synthesis and anti-acetylcholinesterase properties of novel β- and γ-substituted alkoxy organophosphonates

Activated organophosphate (OP) insecticides and chemical agents inhibit acetylcholinesterase (AChE) to form OP-AChE adducts. Whereas the structure of the OP correlates with the rate of inhibition, the structure of the OP-AChE adduct influences the rate at which post-inhibitory reactivation or aging phenomena occurs. In this report, we prepared a panel of β-substituted ethoxy and γ-substituted propoxy phosphonoesters of the type p-NO(2)PhO-P(X)(R)[(O(CH(2))(n)Z] (R=Me, Et; X=O, S; n=2, 3; Z=halogen, OTs) and examined the inhibition of three AChEs by select structures in the panel. The β-fluoroethoxy methylphosphonate analog (R=Me, Z=F, n=2) was the most potent anti-AChE compound comparable (ki ∼6 × 10(6)M(-1)min(-1)) to paraoxon against EEAChE. Analogs with Z=Br, I, or OTs were weak inhibitors of the AChEs, and methyl phosphonates (R=Me) were more potent than the corresponding ethyl phosphonates (R=Et). As expected, analogs with a thionate linkage (PS) were poor inhibitors of the AChEs.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

Frequency and effect of the binding of Mg2+, Mn2+, and Co2+ ions on the guanine base in Watson-Crick and reverse Watson-Crick base pairs

We performed MP2 calculations to elucidate the structure and energetics of the Mg(2+), Mn(2+), and Co(2+) hexahydrated aquaions, and the effect of the metal binding to the N7 atom of (i) a single guanine, (ii) a guanine involved in a Watson-Crick pair, and (iii) a guanine involved in a reverse Watson-Crick base pair. Our comparative analysis of the three aquaions indicates a clear inverse correlation between the radius of the cation and the binding energy, that indeed increases in the order Mn(2+) < Co(2+) < Mg(2+). The trend in the binding energies of the pentahydrated cations to the N7 atom of the guanine is instead Mg(2+) < Mn(2+) < Co(2+), suggesting a rather different bonding scheme that, for the two transition metals, involves back-donation from the aromatic ring of the guanine to their empty d orbitals. In the gas phase, the three hydrated metals significantly stabilize both G-C base pair geometries, Watson-Crick and reverse Watson-Crick, we investigated. Inclusion of a continuous solvent model, however, remarkably reduces this additional stabilization, which becomes almost negligible in the case of the Mg(2+) cation coordinated to the guanine in the standard Watson-Crick geometry. Conversely, all three metal ions sensibly stabilize the reverse Watson-Crick geometry, also in water. Our results are supported by a screening of the structures available in the Protein Data Bank, which clearly indicates that the two transition metals we investigated have a tendency greater than Mg(2+) to coordinate to the N7 atom of guanines, and that there is no clear correlation between the number of guanines in experimental structures with a metal bound to N7 atom and their involvement in Watson-Crick base pairs.

Use of phosphorothioates to identify sites of metal-ion binding in RNA

Single atom substitutions provide an exceptional opportunity to investigate RNA structure and function. Replacing a phosphoryl oxygen with a sulfur represents one of the most common and powerful single atom substitutions and can be used to determine the sites of metal-ion binding. Using functional assays of ribozyme catalysis, based on pre-steady-state kinetics, it is possible to extend this analysis to the transition state, capturing ligands for catalytic metal ions in this fleeting state. In conjunction with data determined from X-ray crystallography, this technique can provide a picture of the environment surrounding catalytic metal ions in both the ground state and the transition state at atomic resolution. Here, we describe the principles of such analysis, explain limitations of the method, and provide a practical example based on our results with the Tetrahymena group I ribozyme.

Probing the role of a secondary structure element at the 5'- and 3'-splice sites in group I intron self-splicing: the tetrahymena L-16 ScaI ribozyme reveals a new role of the G.U pair in self-splicing

Several ribozyme constructs have been used to dissect aspects of the group I self-splicing reaction. The Tetrahymena L-21 ScaI ribozyme, the best studied of these intron analogues, catalyzes a reaction analogous to the first step of self-splicing, in which a 5'-splice site analogue (S) and guanosine (G) are converted into a 5'-exon analogue (P) and GA. This ribozyme preserves the active site but lacks a short 5'-terminal segment (called the IGS extension herein) that forms dynamic helices, called the P1 extension and P10 helix. The P1 extension forms at the 5'-splice site in the first step of self-splicing, and P10 forms at the 3'-splice site in the second step of self-splicing. To dissect the contributions from the IGS extension and the helices it forms, we have investigated the effects of each of these elements at each reaction step. These experiments were performed with the L-16 ScaI ribozyme, which retains the IGS extension, and with 5'- and 3'-splice site analogues that differ in their ability to form the helices. The presence of the IGS extension strengthens binding of P by 40-fold, even when no new base pairs are formed. This large effect was especially surprising, as binding of S is essentially unaffected for S analogues that do not form additional base pairs with the IGS extension. Analysis of a U.U pair immediately 3' to the cleavage site suggests that a previously identified deleterious effect from a dangling U residue on the L-21 ScaI ribozyme arises from a fortuitous active site interaction and has implications for RNA tertiary structure specificity. Comparisons of the affinities of 5'-splice site analogues that form only a subset of base pairs reveal that inclusion of the conserved G.U base pair at the cleavage site of group I introns destabilizes the P1 extension >100-fold relative to the stability of a helix with all Watson-Crick base pairs. Previous structural data with model duplexes and the recent intron structures suggest that this effect can be attributed to partial unstacking of the P1 extension at the G.U step. These results suggest a previously unrecognized role of the G.U wobble pair in self-splicing: breaking cooperativity in base pair formation between P1 and the P1 extensions. This effect may facilitate replacement of the P1 extension with P10 after the first chemical step of self-splicing and release of the ligated exons after the second step of self-splicing.

Novel isosteric, isopolar phosphonate analogs of oligonucleotides: preparation and properties

A synthetic approach leading to novel-type modified oligothymidylates containing an isosteric, isopolar, enzyme-stable C3'-O-P-CH(2)-O-C4'' phosphonate alternative to phosphodiester internucleotide bond was elaborated. The suitable monomers were prepared from 4'-phosphonomethoxy derivatives of alpha-L-threo and beta-D-erythro-2',5'-dideoxythymidine, which were considered interesting as structurally related to nucleoside 5'-monophosphates. The phosphotriester method was applied to the automated synthesis of both homooligomeric phosphonate 15-mer chains and alternating phosphonate-phosphate constructs. The fully modified homooligomers did not hybridize while homooligomers with alternating sequences containing alpha-L-threo-configured units (but not beta-D-erythro-) showed a significant decrease in T(m) values in comparison with natural dT(15). For a comparative study, phosphodiester 4'-CH(3)-substituted oligothymidylate was synthesized and physical studies (NMR, CD, MDS modeling) were undertaken to shed more light on the changes in conformational behavior arising from the chosen structural alterations.

Alkaline phosphatase mono- and diesterase reactions: comparative transition state analysis

Enzyme-catalyzed phosphoryl transfer reactions have frequently been suggested to proceed through transition states that are altered from their solution counterparts. Previous work with Escherichia coli alkaline phosphatase (AP), however, suggests that this enzyme catalyzes the hydrolysis of phosphate monoesters through a loose, dissociative transition state, similar to that in solution. AP also exhibits catalytic promiscuity, with a low level of phosphodiesterase activity, despite the tighter, more associative transition state for phosphate diester hydrolysis in solution. Because AP is evolutionarily optimized for phosphate monoester hydrolysis, it is possible that the active site environment alters the transition state for diester hydrolysis to be looser in its bonding to the incoming and outgoing groups. To test this possibility, we have measured the nonenzymatic and AP-catalyzed rate of reaction for a series of substituted methyl phenyl phosphate diesters. The values of beta(lg) and additional data suggest that the transition state for AP-catalyzed phosphate diester hydrolysis is indistinguishable from that in solution. Instead of altering transition state structure, AP catalyzes phosphoryl transfer reactions by recognizing and stabilizing transition states similar to those in aqueous solution. The AP active site therefore has the ability to recognize different transition states, a property that could assist in the evolutionary optimization of promiscuous activities.

Understanding the transition states of phosphodiester bond cleavage: insights from heavy atom isotope effects

The nucleotides of DNA and RNA are joined by phosphodiester linkages whose synthesis and hydrolysis are catalyzed by numerous essential enzymes. Two prominent mechanisms have been proposed for RNA and protein enzyme catalyzed cleavage of phosphodiester bonds in RNA: (a) intramolecular nucleophilic attack by the 2'-hydroxyl group adjacent to the reactive phosphate; and (b) intermolecular nucleophilic attack by hydroxide, or other oxyanion. The general features of these two mechanisms have been established by physical organic chemical analyses; however, a more detailed understanding of the transition states of these reactions is emerging from recent kinetic isotope effect (KIE) studies. The recent data show interesting differences between the chemical mechanisms and transition state structures of the inter- and intramolecular reactions, as well as provide information on the impact of metal ion, acid, and base catalysis on these mechanisms. Importantly, recent nonenzymatic model studies show that interactions with divalent metal ions, an important feature of many phosphodiesterase active sites, can influence both the mechanism and transition state structure of nonenzymatic phosphodiester cleavage. Such detailed investigations are important because they mimic catalytic strategies employed by both RNA and protein phosphodiesterases, and so set the stage for explorations of enzyme-catalyzed transition states. Application of KIE analyses for this class of enzymes is just beginning, and several important technical challenges remain to be overcome. Nonetheless, such studies hold great promise since they will provide novel insights into the role of metal ions and other active site interactions.

EPR spectroscopic analysis of TAR RNA-metal ion interactions

Metal ion-induced changes in HIV-1 TAR RNA internal dynamics were determined by the changes in EPR spectral width for TAR RNAs containing spin-labeled nucleotides (U23, U25, U38, and U40). This gave a dynamic signature for each of 10 metal ions studied, which fell into one of three distinct groups. While Li(+) and K(+) had little effect on TAR RNA internal dynamics, Na(+) unexpectedly had a dynamic signature that was similar to Ca(2+) and Sr(2+), with a decrease in mobility at U23 and U38, little or no change at U25, and an increase in mobility at U40. Mg(2+), Co(2+), Ni(2+), Zn(2+), and Ba(2+) had similar effects on U23, U38, and U40, but the mobility of U25 was markedly increased. Our results show that RNA dynamics change upon metal binding to the TAR RNA bulge, indicating that RNA structure adapts to accommodate metal ions of different size and coordination properties.

Extraordinarily slow binding of guanosine to the Tetrahymena group I ribozyme: implications for RNA preorganization and function

The Tetrahymena ribozyme derived from the self-splicing group I intron binds a 5'-splice site analog (S) and guanosine (G), catalyzing their conversion to a 5'-exon analog (P) and GA. Herein, we show that binding of guanosine is exceptionally slow, limiting the reaction at near neutral pH. Our results implicate a conformational rearrangement on guanosine binding, likely because the binding site is not prearranged in the absence of ligand. The fast accommodation of guanosine (10(2) to 10(3) x s(-1)) and prior structural data suggest local rather than global rearrangements, raising the possibility that folding of this and perhaps other large RNAs is not fully cooperative. Guanosine binding is accelerated by addition of residues that form helices, referred to as P9.0 and P10, immediately 5' and 3' to the guanosine. These rate enhancements provide evidence for binding intermediates that have the adjacent helices formed before accommodation of guanosine into its binding site. Because the ability to form the P9.0 and P10 helices distinguishes the guanosine at the correct 3'-splice site from other guanosine residues, the faster binding of the correct guanosine can enhance specificity of 3'-splice site selection. Thus, paradoxically, the absence of a preformed binding site and the resulting slow guanosine binding can contribute to splicing specificity by providing an opportunity for the adjacent helices to increase the rate of binding of the guanosine specifying the 3'-splice site.