Recognition of the spliceosomal branch site RNA helix on the basis of surface and electrostatic features

Interaction between the Spliceosomal Pre-mRNA Branch Site and U2 snRNP Protein p14

We have probed the molecular basis of recognition between human spliceosomal U2 snRNP protein p14 and RNA targets representing the intron branch site region. Interaction of an RNA duplex representing the branch site helix perturbed at least 10 nuclear magnetic resonance cross-peaks of (15)N-labeled p14. However, similar chemical shift changes were observed upon interaction with a duplex without the bulged branch site residue, suggesting that binding of p14 to RNA is nonspecific and does not recognize the branch site. We propose that the p14-RNA interaction screens charges on the backbone of the branch site during spliceosome assembly.

Features of CPB: a Poisson-Boltzmann solver that uses an adaptive Cartesian grid

The capabilities of an adaptive Cartesian grid (ACG)-based Poisson-Boltzmann (PB) solver (CPB) are demonstrated. CPB solves various PB equations with an ACG, built from a hierarchical octree decomposition of the computational domain. This procedure decreases the number of points required, thereby reducing computational demands. Inside the molecule, CPB solves for the reaction-field component (ϕrf ) of the electrostatic potential (ϕ), eliminating the charge-induced singularities in ϕ. CPB can also use a least-squares reconstruction method to improve estimates of ϕ at the molecular surface. All surfaces, which include solvent excluded, Gaussians, and others, are created analytically, eliminating errors associated with triangulated surfaces. These features allow CPB to produce detailed surface maps of ϕ and compute polar solvation and binding free energies for large biomolecular assemblies, such as ribosomes and viruses, with reduced computational demands compared to other Poisson-Boltzmann equation solvers. The reader is referred to http://www.continuum-dynamics.com/solution-mm.html for how to obtain the CPB software.

Progress in developing Poisson-Boltzmann equation solvers

This review outlines the recent progress made in developing more accurate and efficient solutions to model electrostatics in systems comprised of bio-macromolecules and nano-objects, the last one referring to objects that do not have biological function themselves but nowadays are frequently used in biophysical and medical approaches in conjunction with bio-macromolecules. The problem of modeling macromolecular electrostatics is reviewed from two different angles: as a mathematical task provided the specific definition of the system to be modeled and as a physical problem aiming to better capture the phenomena occurring in the real experiments. In addition, specific attention is paid to methods to extend the capabilities of the existing solvers to model large systems toward applications of calculations of the electrostatic potential and energies in molecular motors, mitochondria complex, photosynthetic machinery and systems involving large nano-objects.

Impact of base pair identity 5' to the spliceosomal branch site adenosine on branch site conformation

The branch site helix from Saccharomyces cerevisiae with pseudouridine (ψ) incorporated in a phylogenetically conserved position of U2 snRNA features an extrahelical branch site adenosine (A) that forms a base triple interaction with the minor groove edge of a widely conserved purine(U2 strand)-pyrimidine(intron strand) (R(U2)-Y(intron)) base pair two positions upstream. In these studies, NMR spectra of a duplex in which 2-aminopurine (2ap), a fluorescent analog of adenine lacking the proposed hydrogen bond donor, was substituted for the branch site A, indicated that the substitution does not alter the extrahelical position of the branch site residue; thus, it appears that a hydrogen bond between the adenine amino group and the R-Y pair is not obligatory for stabilization of the extrahelical conformation. In contrast, reversal of the orientation of A(U2)-U(intron) to U(U2)-A(intron) resulted in an intrahelical position for the branch site A or 2ap. Fluorescence intensity of 2ap substituted for the branch site A with the original R(U2)-Y(intron) orientation (AU or GC) was high, consistent with an extrahelical position, whereas fluorescence in helices with the reversed R-Y orientation, or with a mismatched pair (A-U → G•A or U•C), was markedly quenched, implying that the residue was stacked in the helix. The A 5' to the branch site residue was not extrahelical in any of the duplexes. These findings suggest that the R(U2)-Y(intron) base pair orientation in the ψ-dependent branch site helix plays an important role in positioning the branch site A for recognition and/or function.

The ionic atmosphere around A-RNA: Poisson-Boltzmann and molecular dynamics simulations

The distributions of different cations around A-RNA are computed by Poisson-Boltzmann (PB) equation and replica exchange molecular dynamics (MD). Both the nonlinear PB and size-modified PB theories are considered. The number of ions bound to A-RNA, which can be measured experimentally, is well reproduced in all methods. On the other hand, the radial ion distribution profiles show differences between MD and PB. We showed that PB results are sensitive to ion size and functional form of the solvent dielectric region but not the solvent dielectric boundary definition. Size-modified PB agrees with replica exchange molecular dynamics much better than nonlinear PB when the ion sizes are chosen from atomistic simulations. The distribution of ions 14 Å away from the RNA central axis are reasonably well reproduced by size-modified PB for all ion types with a uniform solvent dielectric model and a sharp dielectric boundary between solvent and RNA. However, this model does not agree with MD for shorter distances from the A-RNA. A distance-dependent solvent dielectric function proposed by another research group improves the agreement for sodium and strontium ions, even for shorter distances from the A-RNA. However, Mg(2+) distributions are still at significant variances for shorter distances.

A Fast and Robust Poisson-Boltzmann Solver Based on Adaptive Cartesian Grids

An adaptive Cartesian grid (ACG) concept is presented for the fast and robust numerical solution of the 3D Poisson-Boltzmann Equation (PBE) governing the electrostatic interactions of large-scale biomolecules and highly charged multi-biomolecular assemblies such as ribosomes and viruses. The ACG offers numerous advantages over competing grid topologies such as regular 3D lattices and unstructured grids. For very large biological molecules and multi-biomolecule assemblies, the total number of grid-points is several orders of magnitude less than that required in a conventional lattice grid used in the current PBE solvers thus allowing the end user to obtain accurate and stable nonlinear PBE solutions on a desktop computer. Compared to tetrahedral-based unstructured grids, ACG offers a simpler hierarchical grid structure, which is naturally suited to multigrid, relieves indirect addressing requirements and uses fewer neighboring nodes in the finite difference stencils. Construction of the ACG and determination of the dielectric/ionic maps are straightforward, fast and require minimal user intervention. Charge singularities are eliminated by reformulating the problem to produce the reaction field potential in the molecular interior and the total electrostatic potential in the exterior ionic solvent region. This approach minimizes grid-dependency and alleviates the need for fine grid spacing near atomic charge sites. The technical portion of this paper contains three parts. First, the ACG and its construction for general biomolecular geometries are described. Next, a discrete approximation to the PBE upon this mesh is derived. Finally, the overall solution procedure and multigrid implementation are summarized. Results obtained with the ACG-based PBE solver are presented for: (i) a low dielectric spherical cavity, containing interior point charges, embedded in a high dielectric ionic solvent - analytical solutions are available for this case, thus allowing rigorous assessment of the solution accuracy; (ii) a pair of low dielectric charged spheres embedded in a ionic solvent to compute electrostatic interaction free energies as a function of the distance between sphere centers; (iii) surface potentials of proteins, nucleic acids and their larger-scale assemblies such as ribosomes; and (iv) electrostatic solvation free energies and their salt sensitivities - obtained with both linear and nonlinear Poisson-Boltzmann equation - for a large set of proteins. These latter results along with timings can serve as benchmarks for comparing the performance of different PBE solvers.

Revisiting the association of cationic groove-binding drugs to DNA using a Poisson-Boltzmann approach

Proper modeling of nonspecific salt-mediated electrostatic interactions is essential to understanding the binding of charged ligands to nucleic acids. Because the linear Poisson-Boltzmann equation (PBE) and the more approximate generalized Born approach are applied routinely to nucleic acids and their interactions with charged ligands, the reliability of these methods is examined vis-à-vis an efficient nonlinear PBE method. For moderate salt concentrations, the negative derivative, SK(pred), of the electrostatic binding free energy, DeltaG(el), with respect to the logarithm of the 1:1 salt concentration, [M(+)], for 33 cationic minor groove drugs binding to AT-rich DNA sequences is shown to be consistently negative and virtually constant over the salt range considered (0.1-0.4 M NaCl). The magnitude of SK(pred) is approximately equal to the charge on the drug, as predicted by counterion condensation theory (CCT) and observed in thermodynamic binding studies. The linear PBE is shown to overestimate the magnitude of SK(pred), whereas the nonlinear PBE closely matches the experimental results. The PBE predictions of SK(pred) were not correlated with DeltaG(el) in the presence of a dielectric discontinuity, as would be expected from the CCT. Because this correlation does not hold, parameterizing the PBE predictions of DeltaG(el) against the reported experimental data is not possible. Moreover, the common practice of extracting the electrostatic and nonelectrostatic contributions to the binding of charged ligands to biopolyelectrolytes based on the simple relation between experimental SK values and the electrostatic binding free energy that is based on CCT is called into question by the results presented here. Although the rigid-docking nonlinear PB calculations provide reliable predictions of SK(pred), at least for the charged ligand-nucleic acid complexes studied here, accurate estimates of DeltaG(el) will require further development in theoretical and experimental approaches.

Simulations of RNA interactions with monovalent ions

RNA folding and binding reactions are mediated by interactions with ions that make up the surrounding aqueous electrolytic milieu. Although Mg(2+) ions are often implicated as being crucial for RNA folding, it is known that folding is feasible in high concentrations of monovalent alkali-halide salts. Experiments have yielded important information regarding the salt dependence of RNA stability. Recent work has shown that molecular simulations based on explicit representations of solvent molecules and monovalent ions can also provide useful insights regarding the ionic atmospheres around model RNA systems. These insights can help rationalize intriguing observations regarding the dependence of RNA stability on cation type providing one pays attention to important considerations that go into the proper design of molecular simulations. These issues are discussed in detail and the methods are applied to an A-form RNA and B-form DNA sequence. The results of these simulations are compared to previous work on a kissing-loop system with analogous sequence. In particular, ionic atmospheres obtained from molecular simulations are compared to those obtained using the nonlinear Poisson Boltzmann model for continuum electrostatics for these three different nucleic acid systems. The comparisons indicate reasonable agreement in terms of coarse-grained observables such as the numbers of counterions accumulated around the solutes. However, details of the ionic atmospheres, captured in terms of spatial free energy density profiles, are quite different between the two approaches. These comparisons suggest the need for improvements in continuum models to capture sequence-specific effects, ion-ion correlation, and the effects of partial dehydration of ions.

Different mechanisms for pseudouridine formation in yeast 5S and 5.8S rRNAs

The selection of sites for pseudouridylation in eukaryotic cytoplasmic rRNA occurs by the base pairing of the rRNA with specific guide sequences within the RNA components of box H/ACA small nucleolar ribonucleoproteins (snoRNPs). Forty-four of the 46 pseudouridines (Psis) in the cytoplasmic rRNA of Saccharomyces cerevisiae have been assigned to guide snoRNAs. Here, we examine the mechanism of Psi formation in 5S and 5.8S rRNA in which the unassigned Psis occur. We show that while the formation of the Psi in 5.8S rRNA is associated with snoRNP activity, the pseudouridylation of 5S rRNA is not. The position of the Psi in 5.8S rRNA is guided by snoRNA snR43 by using conserved sequence elements that also function to guide pseudouridylation elsewhere in the large-subunit rRNA; an internal stem-loop that is not part of typical yeast snoRNAs also is conserved in snR43. The multisubstrate synthase Pus7 catalyzes the formation of the Psi in 5S rRNA at a site that conforms to the 7-nucleotide consensus sequence present in other substrates of Pus7. The different mechanisms involved in 5S and 5.8S rRNA pseudouridylation, as well as the multiple specificities of the individual trans factors concerned, suggest possible roles in linking ribosome production to other processes, such as splicing and tRNA synthesis.

The electrostatic characteristics of G.U wobble base pairs

G.U wobble base pairs are the most common and highly conserved non-Watson-Crick base pairs in RNA. Previous surface maps imply uniformly negative electrostatic potential at the major groove of G.U wobble base pairs embedded in RNA helices, suitable for entrapment of cationic ligands. In this work, we have used a Poisson-Boltzmann approach to gain a more detailed and accurate characterization of the electrostatic profile. We found that the major groove edge of an isolated G.U wobble displays distinctly enhanced negativity compared with standard GC or AU base pairs; however, in the context of different helical motifs, the electrostatic pattern varies. G.U wobbles with distinct widening have similar major groove electrostatic potentials to their canonical counterparts, whereas those with minimal widening exhibit significantly enhanced electronegativity, ranging from 0.8 to 2.5 kT/e, depending upon structural features. We propose that the negativity at the major groove of G.U wobble base pairs is determined by the combined effect of the base atoms and the sugar-phosphate backbone, which is impacted by stacking pattern and groove width as a result of base sequence. These findings are significant in that they provide predictive power with respect to which G.U sites in RNA are most likely to bind cationic ligands.

Use of a novel Förster resonance energy transfer method to identify locations of site-bound metal ions in the U2-U6 snRNA complex

U2 and U6 snRNAs pair to form a phylogenetically conserved complex at the catalytic core of the spliceosome. Interactions with divalent metal ions, particularly Mg(II), at specific sites are essential for its folding and catalytic activity. We used a novel Förster resonance energy transfer (FRET) method between site-bound luminescent lanthanide ions and a covalently attached fluorescent dye, combined with supporting stoichiometric and mutational studies, to determine locations of site-bound Tb(III) within the human U2-U6 complex. At pH 7.2, we detected three metal-ion-binding sites in: (1) the consensus ACACAGA sequence, which forms the internal loop between helices I and III; (2) the four-way junction, which contains the conserved AGC triad; and (3) the internal loop of the U6 intra-molecular stem loop (ISL). Binding at each of these sites is supported by previous phosphorothioate substitution studies and, in the case of the ISL site, by NMR. Binding of Tb(III) at the four-way junction and the ISL sites was found to be pH-dependent, with no ion binding observed below pH 6 and 7, respectively. This pH dependence of metal ion binding suggests that the local environment may play a role in the binding of metal ions, which may impact on splicing activity.

A new RNA branching activity: the GIR1 ribozyme

The formation of lariat intermediates during the first step of splicing of group II introns and spliceosomal introns is a well-studied fundamental reaction in molecular biology. Apart from this prominent example, there are surprisingly few occurrences of branched nucleotides or even 2',5'-phosphodiester bonds in biology. We recently described a new ribozyme, the GIR1 branching ribozyme, which catalyzes the formation of a tiny lariat that caps an mRNA. This new example together with work on artificial branching ribozymes and deoxyribozymes shows that branching is facile and points to the possibility that branching reactions could be more prevalent than previously recognized.

Pseudouridine modification in Caenorhabditis elegans spliceosomal snRNAs: unique modifications are found in regions involved in snRNA-snRNA interactions

Background

Pseudouridine (Ψ) is an abundant modified nucleoside in RNA and a number of studies have shown that the presence of Ψ affects RNA structure and function. The positions of Ψ in spliceosomal small nuclear RNAs (snRNAs) have been determined for a number of species but not for the snRNAs from Caenorhabditis elegans (C. elegans), a popular experimental model system of development.

Results

As a prelude to determining the function of or requirement for this modification in snRNAs, we have mapped the positions of Ψ in U1, U2, U4, U5, and U6 snRNAs from worms using a specific primer extension method. As with other species, C. elegans U2 snRNA has the greatest number of Ψ residues, with nine, located in the 5' half of the U2 snRNA. U5 snRNA has three Ψs, in or near the loop of the large stem-loop that dominates the structure of this RNA. U6 and U1 snRNAs each have one Ψ, and two Ψ residues were found in U4 snRNA.

Conclusion

The total number of Ψs found in the snRNAs of C. elegans is significantly higher than the minimal amount found in yeasts but it is lower than that seen in sequenced vertebrate snRNAs. When the actual sites of modification on C. elegans snRNAs are compared with other sequenced snRNAs most of the positions correspond to modifications found in other species. However, two of the positions modified on C. elegans snRNAs are unique, one at position 28 on U2 snRNA and one at position 62 on U4 snRNA. Both of these modifications are in regions of these snRNAs that interact with U6 snRNA either in the spliceosome or in the U4/U6 small nuclear ribonucleoprotein particle (snRNP) and the presence of Ψ may be involved in strengthening the intermolecular association of the snRNAs.