Conformational model of the Holliday junction transition deduced from molecular dynamics simulations

Homologous recombination plays a key role in the restart of stalled replication forks and in the generation of genetic diversity. During this process, two homologous DNA molecules undergo strand exchange to form a four-way DNA (Holliday) junction. In the presence of metal ions, the Holliday junction folds into the stacked-X structure that has two alternative conformers. Experiments have revealed the spontaneous transitions between these conformers, but their detailed pathways are not known. Here, we report a series of molecular dynamics simulations of the Holliday junction at physiological and elevated (400 K) temperatures. The simulations reveal new tetrahedral intermediates and suggest a schematic framework for conformer transitions. The tetrahedral intermediates bear resemblance to the junction conformation in complex with a junction-resolving enzyme, T7 endonuclease I, and indeed, one intermediate forms a stable complex with the enzyme as demonstrated in one simulation. We also describe free energy minima for various states of the Holliday junction system, which arise during conformer transitions. The results show that magnesium ions stabilize the stacked-X form and destabilize the open and tetrahedral intermediates. Overall, our study provides a detailed dynamic model of the Holliday junction undergoing a conformer transition.     

The solution structure of the four-way DNA junction at low-salt conditions: a fluorescence resonance energy transfer analysis.

The four-way DNA (Holliday) junction is an important postulated intermediate in the process of genetic recombination. Earlier studies have suggested that the junction exists in two alternative conformations, depending upon the salt concentration present. At high salt concentrations the junction folds into a stacked X structure, while at low salt concentrations the data indicate an extended unstacked conformation. The stereochemical conformation of the four-way DNA junction at low salt (low alkali ion concentration and no alkaline earth ions) was established by comparing the efficiency of fluorescence resonance energy transfer (FRET) between donor and acceptor molecules attached pairwise in three permutations to the 5' termini of the duplex arms. A new variation of FRET was implemented based upon a systematic variation of the fraction of donor labeled single strands. The FRET results indicate that the structure of the four-way DNA junction at low salt exists as an unstacked, extended, square arrangement of the four duplex arms. The donor titration measurements made in the presence of magnesium ions clearly show the folding of the junction into the X stacked structure. In addition, the FRET efficiency can be measured. The fluorescence anisotropy of the acceptor in the presence of Mg2+ during donor titrations was also measured; the FRET efficiency can be calculated from the anisotropy data and the results are consistent with the folded, stacked X structure.     

Interconversion between Parallel and Antiparallel Conformations of a 4H RNA junction in Domain 3 of Foot-and-Mouth Disease Virus IRES Captured by Dynamics Simulations

RNA junctions are common secondary structural elements present in a wide range of RNA species. They play crucial roles in directing the overall folding of RNA molecules as well as in a variety of biological functions. In particular, there has been great interest in the dynamics of RNA junctions, including conformational pathways of fully base-paired 4-way (4H) RNA junctions. In such constructs, all nucleotides participate in one of the four double-stranded stem regions, with no connecting loops. Dynamical aspects of these 4H RNAs are interesting because frequent interchanges between parallel and antiparallel conformations are thought to occur without binding of other factors. Gel electrophoresis and single-molecule fluorescence resonance energy transfer experiments have suggested two possible pathways: one involves a helical rearrangement via disruption of coaxial stacking, and the other occurs by a rotation between the helical axes of coaxially stacked conformers. Employing molecular dynamics simulations, we explore this conformational variability in a 4H junction derived from domain 3 of the foot-and-mouth disease virus internal ribosome entry site (IRES); this junction contains highly conserved motifs for RNA-RNA and RNA-protein interactions, important for IRES activity. Our simulations capture transitions of the 4H junction between parallel and antiparallel conformations. The interconversion is virtually barrier-free and occurs via a rotation between the axes of coaxially stacked helices with a transient perpendicular intermediate. We characterize this transition, with various interhelical orientations, by pseudodihedral angle and interhelical distance measures. The high flexibility of the junction, as also demonstrated experimentally, is suitable for IRES activity. Because foot-and-mouth disease virus IRES structure depends on long-range interactions involving domain 3, the perpendicular intermediate, which maintains coaxial stacking of helices and thereby consensus primary and secondary structure information, may be beneficial for guiding the overall organization of the RNA system in domain 3.     

Structural dynamics of individual Holliday junctions

The four-way DNA (Holliday) junction is the central intermediate of genetic recombination, but the dynamic aspects of this important structure are presently unclear. Although transitions between alternative stacking conformers have been predicted, conventional kinetic studies are precluded by the inability to synchronize the junction in a single conformer in bulk solution. Using single-molecule fluorescence methodology we have been able to detect these transitions. The sequence dependence, the influence of counterions and measured energetic barriers indicate that the conformer transition and branch migration processes share the unstacked, open structure as the common intermediate but have different rate-limiting steps. Relative rates indicate that multiple conformer transitions occur at each intermediate step of branch migration, allowing the junction to reach conformational equilibrium. This provides a mechanism whereby the sequence-dependent conformational bias could determine the extent of genetic exchange upon junction resolution.     

Conformations of an adenine bulge in a DNA octamer and its influence on DNA structure from molecular dynamics simulations.

Molecular dynamics simulations have been applied to the DNA octamer d(GCGCA-GAAC). d(GTTCGCGC), which has an adenine bulge at the center to determine the pathway for interconversion between the stacked and extended forms. These forms are known to be important in the molecular recognition of bulges. From a total of ~35 ns of simulation time with the most recent CHARMM27 force field a variety of distinct conformations and subconformations are found. Stacked and fully looped-out forms are in excellent agreement with experimental data from NMR and x-ray crystallography. Furthermore, in a number of conformations the bulge base associates with the minor groove to varying degrees. Transitions between many of the conformations are observed in the simulations and used to propose a complete transition pathway between the stacked and fully extended conformations. The effect on the surrounding DNA sequence is investigated and biological implications of the accessible conformational space and the suggested transition pathway are discussed, in particular for the interaction of the MS2 replicase operator RNA with its coat protein.     

Kinase conformations: a computational study of the effect of ligand binding.

Protein function is often controlled by ligand-induced conformational transitions. Yet, in spite of the increasing number of three-dimensional crystal structures of proteins in different conformations, not much is known about the driving forces of these transitions. As an initial step toward exploring the conformational and energetic landscape of protein kinases by computational methods, intramolecular energies and hydration free energies were calculated for different conformations of the catalytic domain of cAMP-dependent protein kinase (cAPK) with a continuum (Poisson) model for the electrostatics. Three protein kinase crystal structures for ternary complexes of cAPK with the peptide inhibitor PKI(5-24) and ATP or AMP-PNP were modeled into idealized intermediate and open conformations. Concordant with experimental observation, we find that the binding of PKI(5-24) is more effective in stabilizing the closed and intermediate forms of cAPK than ATP. PKI(5-24) seems to drive the final closure of the active site cleft from intermediate to closed state because ATP does not distinguish between these two states. Binding of PKI(5-24) and ATP is energetically additive.     

A comparison of the properties and the solution structure for RNA and DNA quadruplexes which contain two GGAGG sequences joined with a tetranucleotide linker

We have determined solution structure of r(GGAGGUUUUGGAGG) (R14) by NMR; the RNA 14-mer forms an intra-strand parallel quadruplex with a G-tetrad and a hexad, in which a G-tetrad core is augmented by association of two A residues. The quadruplex further forms a dimer through stacking interaction between the hexads. In order to obtain insight into the difference between RNA and DNA quadruplexes, we synthesized the corresponding DNA 14-mer, d(GGAGGTTTTGGAGG) (D14), and examined its properties and structure by CD, gel electrophoresis, and NMR. K+ ions increased the thermal stability of both R14 and D14 structures. The binding affinity of K+ ions to R14 was much higher than that to D14. The CD and gel electrophoretic studies suggest that D14 forms a quadruplex entirely different from that of R14 in the presence of K+ ions; two molecules of D14 form a quadruplex with both antiparallel and parallel strand alignments and with diagonal loops at both ends of the stacked G-tetrads. The NMR study also gave results that are consistent with such structure: alternate glycosidic conformation, 5'G(syn)-G(anti)3', and characteristic chemical shift data observed for many quadruplexes containing diagonal TTTT loops.     

Structural factors that determine selectivity of a high fidelity DNA polymerase for deoxy-, dideoxy-, and ribonucleotides

In addition to discriminating against base pair mismatches, DNA polymerases exhibit a high degree of selectivity for deoxyribonucleotides over ribo- or dideoxynucleotides. It has been proposed that a single active site residue (steric gate) blocks productive binding of nucleotides containing 2'-hydroxyls. Although this steric gate plays a role in sugar moiety discrimination, its interactions do not account fully for the observed behavior of mutants. Here we present 10 high resolution crystal structures and enzyme kinetic analyses of Bacillus DNA polymerase I large fragment variants complexed with deoxy-, ribo-, and dideoxynucleotides and a DNA substrate. Taken together, these data present a more nuanced and general mechanism for nucleotide discrimination in which ensembles of intermediate conformations in the active site trap non-cognate substrates. It is known that the active site O-helix transitions from an open state in the absence of nucleotide substrates to a ternary complex closed state in which the reactive groups are aligned for catalysis. Substrate misalignment in the closed state plays a fundamental part in preventing non-cognate nucleotide misincorpation. The structures presented here show that additional O-helix conformations intermediate between the open and closed state extremes create an ensemble of binding sites that trap and misalign non-cognate nucleotides. Water-mediated interactions, absent in the fully closed state, play an important role in formation of these binding sites and can be remodeled to accommodate different non-cognate substrates. This mechanism may extend also to base pair discrimination.     

Mismatch-induced conformational distortions in polymerase beta support an induced-fit mechanism for fidelity

Molecular dynamics simulations of DNA polymerase (pol) beta complexed with different incorrect incoming nucleotides (G x G, G x T, and T x T template base x incoming nucleotide combinations) at the template-primer terminus are analyzed to delineate structure-function relationships for aberrant base pairs in a polymerase active site. Comparisons, made to pol beta structure and motions in the presence of a correct base pair, are designed to gain atomically detailed insights into the process of nucleotide selection and discrimination. In the presence of an incorrect incoming nucleotide, alpha-helix N of the thumb subdomain believed to be required for pol beta's catalytic cycling moves toward the open conformation rather than the closed conformation as observed for the correct base pair (G x C) before the chemical reaction. Correspondingly, active-site residues in the microenvironment of the incoming base are in intermediate conformations for non-Watson-Crick pairs. The incorrect incoming nucleotide and the corresponding template residue assume distorted conformations and do not form Watson-Crick bonds. Furthermore, the coordination number and the arrangement of ligands observed around the catalytic and nucleotide binding magnesium ions are mismatch specific. Significantly, the crucial nucleotidyl transferase reaction distance (P(alpha)-O3') for the mismatches between the incoming nucleotide and the primer terminus is not ideally compatible with the chemical reaction of primer extension that follows these conformational changes. Moreover, the extent of active-site distortion can be related to experimentally determined rates of nucleotide misincorporation and to the overall energy barrier associated with polymerase activity. Together, our studies provide structure-function insights into the DNA polymerase-induced constraints (i.e., alpha-helix N conformation, DNA base pair bonding, conformation of protein residues in the vicinity of dNTP, and magnesium ions coordination) during nucleotide discrimination and pol beta-nucleotide interactions specific to each mispair and how they may regulate fidelity. They also lend further support to our recent hypothesis that additional conformational energy barriers are involved following nucleotide binding but prior to the chemical reaction.     

Conformational flexibility of four-way junctions in RNA

Helical junctions are common architectural features in RNA. They are particularly important in autonomously folding molecules, as exemplified by the hairpin ribozyme. We have used single-molecule fluorescence spectroscopy to study the dynamic properties of the perfect (4H) four-way helical junction derived from the hairpin ribozyme. In the presence of Mg(2+), the junction samples parallel and antiparallel conformations and both stacking conformers, with a bias towards one antiparallel stacking conformer. There is continual interconversion between the forms, such that there are several transitions per second under physiological conditions. Our data suggest that interconversion proceeds via an open intermediate with reduced cation binding in which coaxial stacking between helices is disrupted. The rate of interconversion becomes slower at higher Mg(2+) concentrations, yet the activation barrier decreases under these conditions, indicating that entropic effects are important. Transitions also occur in the presence of Na(+) only; however, the coaxial stacking appears incomplete under these conditions. The polymorphic and dynamic character of the four-way RNA junction provides a source of structural diversity, from which particular conformations required for biological function might be stabilised by additional RNA interactions or protein binding.     

The stereochemistry of a four-way DNA junction: a theoretical study.

The stereochemical conformation of the four-way helical junction in DNA (the Holliday junction; the postulated central intermediate of genetic recombination) has been analysed, using molecular mechanical computer modelling. A version of the AMBER program package was employed, that had been modified to include the influence of counterions and a global optimisation procedure. Starting from an extended planar structure, the conformation was varied in order to minimise the energy, and we discuss three structures obtained by this procedure. One structure is closely related to a square-planar cross, in which there is no stacking interaction between the four double helical stems. This structure is probably closely similar to that observed experimentally in the absence of cations. The remaining two structures are based on related, yet distinct, conformations, in which there is pairwise coaxial stacking of neighbouring stems. In these structures, the four DNA stems adopt the form of two quasi-continuous helices, in which base stacking is very similar to that found in standard B-DNA geometry. The two stacked helices so formed are not aligned parallel to each other, but subtend an angle of approximately 60 degrees. The strands that exchange between one stacked helix and the other are disposed about the smaller angle of the cross (i.e. 60 degrees rather than 120 degrees), generating an approximately antiparallel alignment of DNA sequences. This structure is precisely the stacked X-structure proposed on the basis of experimental data. The calculations indicate distortions from standard B-DNA conformation that are required to adopt the stacked X-structure; a widening of the minor groove at the junction, and reorientation of the central phosphate groups of the exchanging strands. An important feature of the stacked X-structure is that it presents two structurally distinct sides. These may be recognised differently by enzymes, providing a rationalisation for the points of cleavage by Holliday resolvases.     

Complementary strand relocation may play vital roles in RecA-based homology recognition

RecA-family proteins mediate homologous recombination and recombinational DNA repair through homology search and strand exchange. Initially, the protein forms a filament with the incoming single-stranded DNA (ssDNA) bound in site I. The RecA–ssDNA filament then binds double-stranded DNA (dsDNA) in site II. Non-homologous dsDNA rapidly unbinds, whereas homologous dsDNA undergoes strand exchange yielding heteroduplex dsDNA in site I and the leftover outgoing strand in site II. We show that applying force to the ends of the complementary strand significantly retards strand exchange, whereas applying the same force to the outgoing strand does not. We also show that crystallographically determined binding site locations require an intermediate structure in addition to the initial and final structures. Furthermore, we demonstrate that the characteristic dsDNA extension rates due to strand exchange and free RecA binding are the same, suggesting that relocation of the complementary strand from its position in the intermediate structure to its position in the final structure limits both rates. Finally, we propose that homology recognition is governed by transitions to and from the intermediate structure, where the transitions depend on differential extension in the dsDNA. This differential extension drives strand exchange forward for homologs and increases the free energy penalty for strand exchange of non-homologs.     

Kinetic Dissection of the Pre-existing Conformational Equilibrium in the Trypsin Fold

Structural biology has recently documented the conformational plasticity of the trypsin fold for both the protease and zymogen in terms of a pre-existing equilibrium between closed (E*) and open (E) forms of the active site region. How such plasticity is manifested in solution and affects ligand recognition by the protease and zymogen is poorly understood in quantitative terms. Here we dissect the E*-E equilibrium with stopped-flow kinetics in the presence of excess ligand or macromolecule. Using the clotting protease thrombin and its zymogen precursor prethrombin-2 as relevant models we resolve the relative distribution of the E* and E forms and the underlying kinetic rates for their interconversion. In the case of thrombin, the E* and E forms are distributed in a 1:4 ratio and interconvert on a time scale of 45 ms. In the case of prethrombin-2, the equilibrium is shifted strongly (10:1 ratio) in favor of the closed E* form and unfolds over a faster time scale of 4.5 ms. The distribution of E* and E forms observed for thrombin and prethrombin-2 indicates that zymogen activation is linked to a significant shift in the pre-existing equilibrium between closed and open conformations that facilitates ligand binding to the active site. These findings broaden our mechanistic understanding of how conformational transitions control ligand recognition by thrombin and its zymogen precursor prethrombin-2 and have direct relevance to other members of the trypsin fold.     

Uranyl photoprobing of a four-way DNA junction: evidence for specific metal ion binding.

Metal ions are very important in mediating the folding of nucleic acids, as exemplified by the folding of the four-way DNA junction into the stacked X-conformation. Uranyl ion-mediated photocleavage provides a method for the localization of high-affinity ion binding sites in nucleic acids, and we have applied this to the four-way DNA junction. We have made the following observations. (i) Uranyl ions (UO2(2+)) suppressed the reactivity of junction thymine bases against attack by osmium tetroxide, indicating that the uranyl ion induces folding of the junction into a stacked conformation. (ii) DNA located immediately at the point of strand exchange on the two exchanging strands was hypersensitive to uranyl photocleavage. The relative hypersensitivity was considerably accentuated when the photocleavage was carried out in the presence of citrate ions. This suggests the presence of a tight binding site for the uranyl ion in the junction. (iii) The same positions were significantly protected from uranyl cleavage by the presence of hexamminecobalt (III) or spermidine. These ions are known to induce the folded conformation of the four-way junction with high efficiency, suggesting a direct competition between the ions. By contrast, magnesium ions failed to generate a similar protection against photocleavage. These results suggest that the uranyl, hexamminecobalt (III) and spermidine ions compete for the same high affinity binding site on the junction. This site is located at the centre of the junction, at the point where the exchanging strands pass between the stacked helices. We believe that we have observed the first known example of a metal ion 'footprint' on a folded nucleic acid structure.     

Molecular Dynamics Study of the Opening Mechanism for DNA Polymerase I

During DNA replication, DNA polymerases follow an induced fit mechanism in order to rapidly distinguish between correct and incorrect dNTP substrates. The dynamics of this process are crucial to the overall effectiveness of catalysis. Although X-ray crystal structures of DNA polymerase I with substrate dNTPs have revealed key structural states along the catalytic pathway, solution fluorescence studies indicate that those key states are populated in the absence of substrate. Herein, we report the first atomistic simulations showing the conformational changes between the closed, open, and ajar conformations of DNA polymerase I in the binary (enzyme∶DNA) state to better understand its dynamics. We have applied long time-scale, unbiased molecular dynamics to investigate the opening process of the fingers domain in the absence of substrate for B. stearothermophilis DNA polymerase in silico. These simulations are biologically and/or physiologically relevant as they shed light on the transitions between states in this important enzyme. All closed and ajar simulations successfully transitioned into the fully open conformation, which is known to be the dominant binary enzyme-DNA conformation from solution and crystallographic studies. Furthermore, we have detailed the key stages in the opening process starting from the open and ajar crystal structures, including the observation of a previously unknown key intermediate structure. Four backbone dihedrals were identified as important during the opening process, and their movements provide insight into the recognition of dNTP substrate molecules by the polymerase binary state. In addition to revealing the opening mechanism, this study also demonstrates our ability to study biological events of DNA polymerase using current computational methods without biasing the dynamics.     

RNA approaches the B‐form in stacked single strand dinucleotide contexts

Duplex RNA adopts an A‐form structure, while duplex DNA interconverts between the A‐ and B‐forms depending on the environment. The C2′‐endo sugar pucker seen in B‐form DNA can occur infrequently in ribose sugars as well, but RNA is not understood to assume B‐form conformations. Through analysis of over 45,000 stacked single strand dinucleotide (SSD) crystal structure conformations, this study demonstrates that RNA is capable of adopting a wide conformational range between the canonical A‐ and B‐forms at the localized SSD level, including many B‐form‐like conformations. It does so through C2′‐endo ribose conformations in one or both nucleotides, and B‐form‐like neighboring base stacking patterns. As chemical reactions on nucleic acids involve localized changes in chemical bonds, the understanding of how enzymes distinguish between DNA and RNA nucleotides is altered by the energetic accessibility of these rare B‐form‐like RNA SSD conformations. The existence of these conformations also has direct implications in parametrization of molecular mechanics energy functions used extensively to model nucleic acid behavior., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 65–82, 2016     

Integrin αIIbβ3 intermediates: From molecular dynamics to adhesion assembly

The platelet integrin α_IIbβ₃ undergoes long-range conformational transitions associated with its functional conversion from inactive (low-affinity) to active (high-affinity) during hemostasis. Although new conformations that are intermediate between the well-characterized bent and extended states have been identified, their molecular dynamic properties and functions in the assembly of adhesions remain largely unexplored. In this study, we evaluated the properties of intermediate conformations of integrin α_IIbβ₃ and characterized their effects on the assembly of adhesions by combining all-atom simulations, principal component analysis, and mesoscale modeling. Our results show that in the low-affinity, bent conformation, the integrin ectodomain tends to pivot around the legs; in intermediate conformations, the headpiece becomes partially extended, away from the lower legs. In the fully open, active state, α_IIbβ₃ is flexible, and the motions between headpiece and lower legs are accompanied by fluctuations of the transmembrane helices. At the mesoscale, bent integrins form only unstable adhesions, but intermediate or open conformations stabilize the adhesions. These studies reveal a mechanism by which small variations in ligand binding affinity and enhancement of the ligand-bound lifetime in the presence of actin retrograde flow stabilize α_IIbβ₃ integrin adhesions.     

Crystal structure analysis of the exocytosis-sensitive phosphoprotein, pp63/parafusin (phosphoglucomutase), from Paramecium reveals significant conformational variability.

During exocytosis of dense-core secretory vesicles (trichocysts) in Paramecium, the protein pp63/parafusin (pp63/pf) is transiently dephosphorylated. We report here the structures of two crystal forms of one isoform of this protein which has a high degree of homology with rabbit phosphoglucomutase, whose structure has been reported. As expected, both proteins possess highly similar structures, showing the same four domains forming two lobes with an active-site crevice in between. The two X-ray structures that we report here were determined after crystallization in the presence of sulfate and tartrate, and show the lobes arranged as a closed and an open conformation, respectively. While both conformations possess a bound divalent cation, only the closed (sulfate-bound) conformation shows bound sulfate ions in the "phosphate-transfer site" near the catalytic serine residue and in the "phosphate-binding site". Comparison with the open form shows that the latter dianion is placed in the centre of three arginine residues, one contributed by subunit II and two by subunit IV, suggesting that it causes a contraction of the arginine triangle, which establishes the observed conformational closure of the lobes. It is therefore likely that the closed conformation forms only when a phosphoryl group is bound to the phosphate-binding site. The previously published structure of rabbit phosphoglucomutase is intermediate between these two conformers. Several of the known reversible phosphorylation sites of pp63/pf-1 are at positions critical for transition between the conformations and for binding of the ligands and thus give hints as to possible roles of pp63/pf-1 in the course of exocytosis.     

Influence of copper(II) and magnesium(II) ions on the ciprofloxacin binding to DNA

The influence of magnesium(II) and copper(II) ions on the binding of ciprofloxacin to double stranded calf thymus DNA was studied by fluorescence emission spectroscopy, ultraviolet- and circular dichroism (CD) spectroscopy. The interaction of ciprofloxacin and copper(II) ions was followed by strong fluorescence quenching which was almost unaffected by the presence of DNA. On the other hand, only a slight decrease in fluorescence emission intensity, which was enhanced in the presence of DNA, was observed for ciprofloxacin interaction with magnesium(II) ions. Furthermore, magnesium(II) ions increase the thermal stability of the DNA, while, in the presence of ciprofloxacin, the degree of stabilisation is smaller. In contrast, copper(II) ions destabilise double helical DNA to heat, while ciprofloxacin slightly affects only the second transition of the biphasic melting curve of calf thymus DNA. Magnesium(II) ions at 25 degrees C induce conformational transitions of DNA at concentrations of 1.5 mM and 2.5 M, as monitored by CD. On the other hand copper(II) ions induce only one conformational transition, at a concentration of 12.7 microM. At higher concentrations of copper(II) ions (c>700 microM) DNA starts to precipitate. Significant changes in the CD spectra of DNA were observed after addition of ciprofloxacin to a solution containing DNA and copper(II) ions, but not to DNA and magnesium(II) ions. Based on our spectroscopic results, we propose that copper(II) ions are not directly involved into ciprofloxacin binding to DNA via phosphate groups as it has been suggested for magnesium(II) ions.     

Cation-specific structural accommodation within a catalytic RNA

Metal ions facilitate the folding of the hairpin ribozyme but do not participate directly in catalysis. The metal complex cobalt(III) hexaammine supports folding and activity of the ribozyme and also mediates specific internucleotide photocrosslinks, several of which retain catalytic ability. These crosslinks imply that the active core structure organized by [Co(NH3)6]3+ is different from that organized by Mg2+ and that revealed in the crystal structure [Rupert, P. B., and Ferre-D'Amare, A. R. (2001) Nature 410, 780-786] (1). Residues U+2 and C+3 of the substrate, in particular, adopt different conformations in [Co(NH3)6]3+. U+2 is bulged out of loop A and stacked on residue G36, whereas the nucleotide at position +3 is stacked on G8, a nucleobase crucial for catalysis. Cleavage kinetics performed with +2 variants and a C+3 U variant correlate with the crosslinking observations. Variants that decreased cleavage rates in magnesium up to 70-fold showed only subtle decreases or even increases in observed rates when assayed in [Co(NH3)6]3+. Here, we propose a model of the [Co(NH3)6]3+-mediated catalytic core generated by MC-SYM that is consistent with these data.