Cation-induced kinetic heterogeneity of the intron-exon recognition in single group II introns.

RNA is commonly believed to undergo a number of sequential folding steps before reaching its functional fold, i.e., the global minimum in the free energy landscape. However, there is accumulating evidence that several functional conformations are often in coexistence, corresponding to multiple (local) minima in the folding landscape. Here we use the 5'-exon-intron recognition duplex of a self-splicing ribozyme as a model system to study the influence of Mg(2+) and Ca(2+) on RNA tertiary structure formation. Bulk and single-molecule spectroscopy reveal that near-physiological M(2+) concentrations strongly promote interstrand association. Moreover, the presence of M(2+) leads to pronounced kinetic heterogeneity, suggesting the coexistence of multiple docked and undocked RNA conformations. Heterogeneity is found to decrease at saturating M(2+) concentrations. Using NMR, we locate specific Mg(2+) binding pockets and quantify their affinity toward Mg(2+). Mg(2+) pulse experiments show that M(2+) exchange occurs on the timescale of seconds. This unprecedented combination of NMR and single-molecule Förster resonance energy transfer demonstrates for the first time to our knowledge that a rugged free energy landscape coincides with incomplete occupation of specific M(2+) binding sites at near-physiological M(2+) concentrations. Unconventional kinetics in nucleic acid folding frequently encountered in single-molecule experiments are therefore likely to originate from a spectrum of conformations that differ in the occupation of M(2+) binding sites.

NMR structure of the 5' splice site in the group IIB intron Sc.ai5γ--conformational requirements for exon-intron recognition.

A crucial step of the self-splicing reaction of group II intron ribozymes is the recognition of the 5' exon by the intron. This recognition is achieved by two regions in domain 1 of the intron, the exon-binding sites EBS1 and EBS2 forming base pairs with the intron-binding sites IBS1 and IBS2 located at the end of the 5' exon. The complementarity of the EBS1•IBS1 contact is most important for ensuring site-specific cleavage of the phosphodiester bond between the 5' exon and the intron. Here, we present the NMR solution structures of the d3' hairpin including EBS1 free in solution and bound to the IBS1 7-mer. In the unbound state, EBS1 is part of a flexible 11-nucleotide (nt) loop. Binding of IBS1 restructures and freezes the entire loop region. Mg(2+) ions are bound near the termini of the EBS1•IBS1 helix, stabilizing the interaction. Formation of the 7-bp EBS1•IBS1 helix within a loop of only 11 nt forces the loop backbone to form a sharp turn opposite of the splice site, thereby presenting the scissile phosphate in a position that is structurally unique.

Dynamics may significantly influence the estimation of interatomic distances in biomolecular X-ray structures.

Atomic positions obtained by X-ray crystallography are time and space averages over many molecules in the crystal. Importantly, interatomic distances, calculated between such average positions and frequently used in structural and mechanistic analyses, can be substantially different from the more appropriate time-average and ensemble-average interatomic distances. Using crystallographic B-factors, one can deduce corrections, which have so far been applied exclusively to small molecules, to obtain correct average distances as a function of the type of atomic motion. Here, using 4774 high-quality protein X-ray structures, we study the significance of such corrections for different types of atomic motion. Importantly, we show that for distances shorter than 5 Å, corrections greater than 0.5 Å may apply, especially for noncorrelated or anticorrelated motion. For example, 14% of the studied structures have at least one pair of atoms with a correction of ≥0.5 Å in the case of noncorrelated motion. Using molecular dynamics simulations of villin headpiece, ubiquitin, and SH3 domain unit cells, we demonstrate that the majority of average interatomic distances in these proteins agree with noncorrelated corrections, suggesting that such deviations may be truly relevant. Importantly, we demonstrate that the corrections do not significantly affect stereochemistry and the overall quality of final refined X-ray structures, but can provide marked improvements in starting unrefined models obtained from low-resolution X-ray data. Finally, we illustrate the potential mechanistic and biological significance of the calculated corrections for KcsA ion channel and show that they provide indirect evidence that motions in its selectivity filter are highly correlated.

Conformational averaging in structural biology: issues, challenges and computational solutions.

Most experimental methods in structural biology provide time- and ensemble-averaged signals and, consequently, molecular structures based on such signals often exhibit only idealized, average features. Second, most experimental signals are only indirectly related to real, molecular geometries, and solving a structure typically involves a complicated procedure, which may not always result in a unique solution. To what extent do such conformationally-averaged, non-linear experimental signals and structural models derived from them accurately represent the underlying microscopic reality? Are there some structural motifs that are actually artificially more likely to be "seen" in an experiment simply due to the averaging artifact? Finally, what are the practical consequences of ignoring the averaging effects when it comes to functional and mechanistic implications that we try to glean from experimentally-based structural models? In this review, we critically address the work that has been aimed at studying such questions. We summarize the details of experimental methods typically used in structural biology (most notably nuclear magnetic resonance, X-ray crystallography and different types of spectroscopy), discuss their individual susceptibility to conformational (motional) averaging, and review several theoretical approaches, most importantly molecular dynamics simulations that are increasingly being used to aid experimentalists in interpreting structural biology experiments.

Divalent metal ions promote the formation of the 5'-splice site recognition complex in a self-splicing group II intron.

Group II introns are ribozymes occurring in genes of plants, fungi, lower eukaryotes, and bacteria. These large RNA molecular machines, ranging in length from 400 to 2500 nucleotides, are able to catalyze their own excision from pre-mRNA, as well as to reinsert themselves into RNA or sometimes even DNA. The intronic domain 1 contains two sequences (exon binding sites 1 and 2, EBS1 and EBS2) that pair with their complementary regions at the 3'-end of the 5'-exon (intron binding sites 1 and 2, IBS1 and IBS2) such defining the 5'-splice site. The correct recognition of the 5'-splice site stands at the beginning of the two steps of splicing and is thus crucial for catalysis. It is known that metal ions play an important role in folding and catalysis of ribozymes in general. Here, we characterize the specific metal ion requirements for the formation of the 5'-splice site recognition complex from the mitochondrial yeast group II intron Sc.ai5gamma. Circular dichroism studies reveal that the formation of the EBS1.IBS1 duplex does not necessarily require divalent metal ions, as large amounts of monovalent metal ions also promote the duplex, albeit at a 5000 times higher concentration. Nevertheless, micromolar amounts of divalent metal ions, e.g. Mg2+ or Cd2+, strongly promote the formation of the 5'-splice site. These observations illustrate that a high charge density independent of the nature of the ion is needed for binding EBS1 to IBS1, but divalent metal ions are presumably the better players.