Characterization of the Reaction Path and Transition States for RNA Transphosphorylation Models from Theory and Experiment

Structural modifications as tools in mechanistic studies of the cleavage of RNA phosphodiester linkages

The cleavage of RNA phosphodiester bonds by RNase A and hammerhead ribozyme at neutral pH fundamentally differs from the spontaneous reactions of these bonds under the same conditions. While the predominant spontaneous reaction is isomerization of the 3',5'-phosphodiester linkages to their 2',5'-counterparts, this reaction has never been reported to compete with the enzymatic cleavage reaction, not even as a minor side reaction. Comparative kinetic measurements with structurally modified di-nucleoside monophosphates and oligomeric phosphodiesters have played an important role in clarification of mechanistic details of the buffer-independent and buffer-catalyzed reactions. More recently, heavy atom isotope effects and theoretical calculations have refined the picture. The primary aim of all these studies has been to form a solid basis for mechanistic analyses of the action of more complicated catalytic machineries. In other words, to contribute to conception of a plausible unified picture of RNA cleavage by biocatalysts, such as RNAse A, hammerhead ribozyme and DNAzymes. In addition, structurally modified trinucleoside monophosphates as transition state models for Group I and II introns have clarified some features of the action of large ribozymes.

Phosphoryl-Transfer Reaction in RNA under Alkaline Conditions

The phosphoryl-transfer reaction in RNA under alkaline conditions by exploring the influence of several solvents theoretically was studied. The calculations were carried out by using the M06-2X functional and the solvents were taken as a continuum by using the solvent model density (SMD) method. The main findings show that the O2'-P-O5' angle in the reactants, the free activation energies, and the reaction mechanism are clearly dependent on the dielectric constant of the environment, thus showing that the electrostatic term is the determining factor for this chemical system with two negative charges. Our study seems to indicate that water, the solvent with the greatest dielectric constant, would be the solvent that increases the reaction rate the most. As this outcome was not the case in enzymatic catalysis, one has to conclude that, in the case of proteins as well as for ribozymes, the enzymatic catalysis is not mainly due to the solvent reaction field, but to local electrical fields as a result of enzyme preorganization.

Transition State Features in the Hepatitis Delta Virus Ribozyme Reaction Revealed by Atomic Perturbations

Endonucleolytic ribozymes constitute a class of non-coding RNAs that catalyze single-strand RNA scission. With crystal structures available for all of the known ribozymes, a major challenge involves relating functional data to the physically observed RNA architecture. In the case of the hepatitis delta virus (HDV) ribozyme, there are three high-resolution crystal structures, the product state of the reaction and two precursor variants, with distinct mechanistic implications. Here, we develop new strategies to probe the structure and catalytic mechanism of a ribozyme. First, we use double-mutant cycles to distinguish differences in functional group proximity implicated by the crystal structures. Second, we use a corrected form of the Brønsted equation to assess the functional significance of general acid catalysis in the system. Our results delineate the functional relevance of atomic interactions inferred from structure, and suggest that the HDV ribozyme transition state resembles the cleavage product in the degree of proton transfer to the leaving group.

Mechanistic insights into RNA transphosphorylation from kinetic isotope effects and linear free energy relationships of model reactions

Phosphoryl transfer reactions are ubiquitous in biology and the understanding of the mechanisms whereby these reactions are catalyzed by protein and RNA enzymes is central to reveal design principles for new therapeutics. Two of the most powerful experimental probes of chemical mechanism involve the analysis of linear free energy relations (LFERs) and the measurement of kinetic isotope effects (KIEs). These experimental data report directly on differences in bonding between the ground state and the rate-controlling transition state, which is the most critical point along the reaction free energy pathway. However, interpretation of LFER and KIE data in terms of transition-state structure and bonding optimally requires the use of theoretical models. In this work, we apply density-functional calculations to determine KIEs for a series of phosphoryl transfer reactions of direct relevance to the 2'-O-transphosphorylation that leads to cleavage of the phosphodiester backbone of RNA. We first examine a well-studied series of phosphate and phosphorothioate mono-, di- and triesters that are useful as mechanistic probes and for which KIEs have been measured. Close agreement is demonstrated between the calculated and measured KIEs, establishing the reliability of our quantum model calculations. Next, we examine a series of RNA transesterification model reactions with a wide range of leaving groups in order to provide a direct connection between observed Brønsted coefficients and KIEs with the structure and bonding in the transition state. These relations can be used for prediction or to aid in the interpretation of experimental data for similar non-enzymatic and enzymatic reactions. Finally, we apply these relations to RNA phosphoryl transfer catalyzed by ribonuclease A, and demonstrate the reaction coordinate-KIE correlation is reasonably preserved. A prediction of the secondary deuterium KIE in this reaction is also provided. These results demonstrate the utility of building up knowledge of mechanism through the systematic study of model systems to provide insight into more complex biological systems such as phosphoryl transfer enzymes and ribozymes.