Tertiary structure in transfer ribonucleic acids

MetBP: A Software Tool for Detection of Interaction between Metal Ion-RNA Base Pairs

MOTIVATION

The role of metals in shaping and functioning of RNA is a well established fact and the understanding of that through the analysis of structural data has biological relevance. Often metal ions bind to one or more atoms of the nucleobase of an RNA. This fact becomes more interesting when such bases form a base pair with any other base. Furthermore, when metal ions bind to any residue of an RNA, the secondary structural features of the residue (helix, loop, unpaired etc) are also biologically important. The available metal binding related software tools cannot address such type specific queries.

RESULTS

To fill this limitation, we have designed a software tool, called MetBP, that meets the goal. This tool is a stand-alone command line based tool and has no dependency on the other existing software. It accepts a structure file in mmCIF or PDB format and computes the base pairs and thereafter reports all metals that bind to one or more nucleotides that form pairs with another. It reports binding distance, angles along with base pair stability. It also reports several other important aspects, e.g. secondary structure of the residue in the RNA. MetBP can be used as a generalized metal binding site detection tool for Proteins and DNA as well.

AVAILABILITY

https://github.com/computational-biology/metbp.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

From static to dynamic: the need for structural ensembles and a predictive model of RNA folding and function

To understand RNA, it is necessary to move beyond a descriptive categorization towards quantitative predictions of its molecular conformations and functional behavior. An incisive approach to understanding the function and folding of biological RNA systems involves characterizing small, simple components that are largely responsible for the behavior of complex systems including helix-junction-helix elements and tertiary motifs. State-of-the-art methods have permitted unprecedented insight into the conformational ensembles of these elements revealing, for example, that conformations of helix-junction-helix elements are confined to a small region of the ensemble, that this region is highly dependent on the junction's topology, and that the correct alignment of tertiary motifs may be a rare conformation on the overall folding landscape. Further characterization of RNA components and continued development of experimental and computational methods with the goal of quantitatively predicting RNA folding and functional behavior will be critical to understanding biological RNA systems.

Transient RNA-protein interactions in RNA folding

The RNA folding trajectory features numerous off-pathway folding traps, which represent conformations that are often equally as stable as the native functional ones. Therefore, the conversion between these off-pathway structures and the native correctly folded ones is the critical step in RNA folding. This process, referred to as RNA refolding, is slow, and is represented by a transition state that has a characteristic high free energy. Because this kinetically limiting process occurs in vivo, proteins (called RNA chaperones) have evolved that facilitate the (re)folding of RNA molecules. Here, we present an overview of how proteins interact with RNA molecules in order to achieve properly folded states. In this respect, the discrimination between static and transient interactions is crucial, as different proteins have evolved a multitude of mechanisms for RNA remodeling. For RNA chaperones that act in a sequence-unspecific manner and without the use of external sources of energy, such as ATP, transient RNA–protein interactions represent the basis of the mode of action. By presenting stretches of positively charged amino acids that are positioned in defined spatial configurations, RNA chaperones enable the RNA backbone, via transient electrostatic interactions, to sample a wider conformational space that opens the route for efficient refolding reactions.

Interaction of the ruthenium red cation with nucleic acid double helices

The hexapositive complex cation ruthenium red very effectively stabilizes DNA and RNA double helices against thermal denaturation. In the presence of nucleic acid helices, this symmetric cation acquires an extrinsic CD spectrum near the wavelength of the dye's maximum absorbance. Competition experiments with single-stranded polyd(T) show this induced CD to be the result of selective binding to helical sites. The preferential affinity of ruthenium red for double helical binding sites is so great that it brings about biphasic absorbance- temperature profiles of polyd(A-T) at low [cation]: [polynucleotide phosphate]. The visible CD signal and fraction of helix melting at the upper transition increases with ruthenium red concentration until approximate charge neutrality is reached. These interactions, which have been studied in detail with the poly(U-U) helix as well as polyd(A-T), are likely largely electrostatic, since sufficient [NaCl] eliminates the bipliasic melting of polyd(A-T), renders the ultraviolet absorbance of poly (U) insensitive to ruthenium red, and abolishes the induced CD effects. The bipliasic melting of polyd(A-T) at intermediate [dye] is attributed to saturation of remaining double helical segments by cation migration from newly melted regions- Furthermore, virtually no change was observed in the induced CD upon melting through the first transition, whereas the effect is destroyed upon inciting through the second transition. A quantitative treatment of the data is used to obtain binding site size and association constant for the complex. The induced effect may prove useful in the exploration of exposed nucleic acid helical structure in such complex particles as nucleosomes or ribosomes.

MeRNA: a database of metal ion binding sites in RNA structures

Metal ions are essential for the folding of RNA into stable tertiary structures and for the catalytic activity of some RNA enzymes. To aid in the study of the roles of metal ions in RNA structural biology, we have created MeRNA (Metals in RNA), a comprehensive compilation of all metal binding sites identified in RNA 3D structures available from the PDB and Nucleic Acid Database. Currently, our database contains information relating to binding of 9764 metal ions corresponding to 23 distinct elements, in 256 RNA structures. The metal ion locations were confirmed and ligands characterized using original literature references. MeRNA includes eight manually identified metal-ion binding motifs, which are described in the literature. MeRNA is searchable by PDB identifier, metal ion, method of structure determination, resolution and R-values for X-ray structure and distance from metal to any RNA atom or to water. New structures with their respective binding motifs will be added to the database as they become available. The MeRNA database will further our understanding of the roles of metal ions in RNA folding and catalysis and have applications in structural and functional analysis, RNA design and engineering. The MeRNA database is accessible at .

A thermodynamic framework for the magnesium-dependent folding of RNA

The goal of this review is to present a unified picture of the relationship between ion binding and RNA folding based on recent theoretical and computational advances. In particular, we present a model describing how the association of magnesium ions is coupled to the tertiary structure folding of several well-characterized RNA molecules. This model is developed in terms of the nonlinear Poisson-Boltzmann (NLPB) equation, which provides a rigorous electrostatic description of the interaction between Mg(2+) and specific RNA structures. In our description, most of the ions surrounding an RNA behave as a thermally fluctuating ensemble distributed according to a Boltzmann weighted average of the mean electrostatic potential around the RNA. In some cases, however, individual ions near the RNA may shed some of their surrounding waters to optimize their Coulombic interactions with the negatively charged ligands on the RNA. These chelated ions are energetically distinct from the surrounding ensemble and must be treated explicitly. This model is used to explore several different RNA systems that interact differently with Mg(2+). In each case, the NLPB equation accurately describes the stoichiometric and energetic linkage between Mg(2+) binding and RNA folding without requiring any fitted parameters in the calculation. Based on this model, we present a physical description of how Mg(2+) binds and stabilizes specific RNA structures to promote the folding reaction.

Assembly of core helices and rapid tertiary folding of a small bacterial group I ribozyme

Compact but non-native intermediates have been implicated in the hierarchical folding of several large RNAs, but there is little information on their structure. In this article, ribonuclease and hydroxyl radical cleavage protection assays showed that base pairing of core helices stabilize a compact state of a small group I ribozyme from Azoarcus pre-tRNA(ile). Base pairing of the ribozyme core requires 10-fold less Mg(2+) than stable tertiary interactions, indicating that assembly of helices in the catalytic core represents a distinct phase that precedes the formation of native tertiary structure. Tertiary folding occurs in <100 ms at 37 degrees C. Such rapid folding is unprecedented among group I ribozymes and illustrates the association between structural complexity and folding time. A 3D model of the Azoarcus ribozyme was constructed by identifying homologous sequence motifs in rRNA. The model reveals distinct structural features, such as a large interface between the P4-P6 and P3-P9 domains, that may explain the unusual stability of the Azoarcus ribozyme and the cooperativity of folding.

A physical origin for functional domain structure in nucleic acids as evidenced by cross-linking entropy: II

In Part I, cross-linking entropy (CLE) was proposed as a mechanism that limits the size of functional domains of RNA. To test this hypothesis, the theory is developed into an RNA secondary structure prediction filter which is applied to nearest-neighbor secondary structure (NNSS) algorithms that utilize a free energy (FE) minimization strategy. (The NNSS strategies are also referred to as the dynamic programming algorithm in the literature.) The cross-linking entropy for RNA is derived from a generalized Gaussian polymer chain model where the entropic contributions caused by the formation of base pairs (stacking) in RNA are analysed globally. Local entropic contributions are associated with the freezing out of degrees of freedom in the links. Both global and local entropic effects are strongly influenced by the persistence length. The cross-linking entropy provides a physical origin for the size of functional domains in long nucleic acid sequences and may go further to explain as to why the majority of the domain regions in typical sequences tend to be less than 600 nucleotides in length. In addition, improvements were observed in the "best guess" predictive capacity over NNSS prediction strategies. The thermodynamic distribution is more representative of the expected structures and is strongly governed by such physical parameters as the persistence length and the excluded volume. The CLE appears to generalize the tabulated penalties used in NNSS algorithms. The principal parameter influencing this entropy is the persistence length. The model is shown to accomodate a variable persistence length and is capable of describing the folding dynamics of RNA. A two-state kinetic model based on the CLE principle is used to help elucidate the folding kinetics of functional domains in the group I introns.

Early events in RNA folding

We describe a conceptual framework for understanding the way large RNA molecules fold based on the notion that their free-energy landscape is rugged. A key prediction of our theory is that RNA folding can be described by the kinetic partitioning mechanism (KPM). According to KPM a small fraction of molecules folds rapidly to the native state whereas the remaining fraction is kinetically trapped in a low free-energy non-native state. This model provides a unified description of the way RNA and proteins fold. Single-molecule experiments on Tetrahymena ribozyme, which directly validate our theory, are analyzed using KPM. We also describe the earliest events that occur on microsecond time scales in RNA folding. These must involve collapse of RNA molecules that are mediated by counterion-condensation. Estimates of time scales for the initial events in RNA folding are provided for the Tetrahymena ribozyme.

Role of counterion condensation in folding of the Tetrahymena ribozyme. II. Counterion-dependence of folding kinetics

Condensed counterions contribute to the stability of compact structures in RNA, largely by reducing electrostatic repulsion among phosphate groups. Varieties of cations induce a collapsed state in the Tetrahymena ribozyme that is readily transformed to the catalytically active structure in the presence of Mg2+. Native gel electrophoresis was used to compare the effects of the valence and size of the counterion on the kinetics of this transition. The rate of folding was found to decrease with the charge of the counterion. Transitions in monovalent ions occur 20- to 40-fold faster than transitions induced by multivalent metal ions. These results suggest that multivalent cations yield stable compact structures, which are slower to reorganize to the native conformation than those induced by monovalent ions. The folding kinetics are 12-fold faster in the presence of spermidine3+ than [Co(NH3)6]3+, consistent with less effective stabilization of long-range RNA interactions by polyamines. Under most conditions, the observed folding rate decreases with increasing counterion concentration. In saturating amounts of counterion, folding is accelerated by addition of urea. These observations indicate that reorganization of compact intermediates involves partial unfolding of the RNA. We find that folding of the ribozyme is most efficient in a mixture of monovalent salt and Mg2+. This is attributed to competition among counterions for binding to the RNA. The counterion dependence of the folding kinetics is discussed in terms of the ability of condensed ions to stabilize compact structures in RNA.

Role of counterion condensation in folding of the Tetrahymena ribozyme. I. Equilibrium stabilization by cations

Folding of RNA into an ordered, compact structure requires substantial neutralization of the negatively charged backbone by positively charged counterions. Using a native gel electrophoresis assay, we have examined the effects of counterion condensation upon the equilibrium folding of the Tetrahymena ribozyme. Incubation of the ribozyme in the presence of mono-, di- and trivalent ions induces a conformational state that is capable of rapidly forming the native structure upon brief exposure to Mg2+. The cation concentration dependence of this transition is directly correlated with the charge of the counterion used to induce folding. Substrate cleavage assays confirm the rapid onset of catalytic activity under these conditions. These results are discussed in terms of classical counterion condensation theory. A model for folding is proposed which predicts effects of charge, ionic radius and temperature on counterion-induced RNA folding transitions.

Mg(2+) binding to tRNA revisited: the nonlinear Poisson-Boltzmann model

Our current understanding of Mg(2+) binding to RNA, in both thermodynamic and structural terms, is largely based on classical studies of transfer RNAs. Based on these studies, it is clear that magnesium ions are crucial for stabilizing the folded structure of tRNA. We present here a rigorous theoretical model based on the nonlinear Poisson-Boltzmann (NLPB) equation for understanding Mg(2+) binding to yeast tRNA(Phe). We use this model to interpret a variety of experimental Mg(2+) binding data. In particular, we find that the NLPB equation provides a remarkably accurate description of both the overall stoichiometry and the free energy of Mg(2+) binding to yeast tRNA(Phe) without any fitted parameters. In addition, the model accurately describes the interaction of Mg(2+) with localized regions of the RNA as determined by the pK(a) shift of differently bound fluorophores. In each case, we find that the model also reproduces the univalent salt-dependence and the anticooperativity of Mg(2+) binding. Our results lead us to a thermodynamic description of Mg(2+) binding to yeast tRNA(Phe) based on the NLPB equation. In this model, Mg(2+) binding is simply explained by an ensemble of ions distributed according to a Boltzmann weighted average of the mean electrostatic potential around the RNA. It appears that the entire ensemble of electrostatically bound ions superficially mimics a few strongly coordinated ions. In this regard, we find that Mg(2+) stabilizes the tertiary structure of yeast tRNA(Phe) in part by accumulating in regions of high negative electrostatic potential. These regions of Mg(2+) localization correspond to bound ions that are observed in the X-ray crystallographic structures of yeast tRNA(Phe). Based on our results and the available thermodynamic data, there is no evidence that specifically coordinated Mg ions have a significant role in stabilizing the native tertiary structure of yeast tRNA(Phe) in solution.

Correct folding of a ribozyme induced by nonspecific macromolecules

The 50-nucleotide hammerhead ribozyme HH-S was tested for self-cleavage. The self-cleavage was very inefficient, and only 13% of HH-S was transformed to its cleavage products. Surprisingly, the percentage of cleavage of HH-S was increased to 30% when 1 microg of tRNA was added to the reaction mixture (6 microL). Other macromolecules such as DNAs and proteins were examined to see if they also augmented cleavage of HH-S, and it was found that most of the macromolecules tested, except nucleotide monomers, did indeed enhance HH-S cleavage. The self-cleaving reaction was almost saturated in 30 min, and only 13% of HH-S was cleaved at 37 degrees C for a 70-min reaction, indicating that 87% of HH-S was in kinetically trapped inactive conformations. Time courses for the reaction of the HH-S self-cleavage were also measured in the presence of tRNA, an oligodeoxyribonucleotide, or BSA. Cleavage of HH-S, which had already reached a plateau of 13% cleaved, increased gradually after the addition of the effector molecules. The first-order rate constant for the self-cleavage reaction in the absence of an effector was comparable to that in the presence of BSA, indicating that the effector molecules do not affect the chemical step of self-cleavage. These results demonstrate that a variety of nonspecific macromolecules can induce conformational change of the hammerhead even in such a low concentration as 0.003% (w/v). This conformational change may occur by macromolecular collisions, or nonspecific weak interactions between HH-S and effectors. Alternatively, a molecular crowding effect may cause the conformational change.

Fast folding of a ribozyme by stabilizing core interactions: evidence for multiple folding pathways in RNA

Folding of the Tetrahymena ribozyme under physiological conditions in vitro is limited by slow conversion of long-lived intermediates to the active structure. These intermediates arise because the most stable domain of the ribozyme folds 10-50 times more rapidly than the core region containing helix P3. Native gel electrophoresis and time-resolved X-ray-dependent hydroxyl radical cleavage revealed that mutations that weaken peripheral interactions between domains accelerated folding fivefold, while a point mutation that stabilizes P3 enabled 80 % of the mutant RNA to reach the native conformation within 30 seconds at 22 degrees C. The P3 mutation increased the folding rate of the catalytic core as much as 50-fold, so that both domains of the ribozyme were formed at approximately the same rate. The results show that the ribozyme folds rapidly without significantly populating metastable intermediates when native interactions in the ribozyme core are stabilized relative to peripheral structural elements.

Stitching together RNA tertiary architectures

The powerful explanatory paradigm of molecular biology requiring form to co-evolve with function has again been proven successful when, over the recent two decades, a wealth of biological functions have been uncovered for RNA. Previously considered as a mere mediator of the genetic code, RNA is now acknowledged as a key player in a wide variety of cellular processes. Along with the discovery of novel biological functions of RNA molecules, a number of RNA three-dimensional structures have been solved which beautifully demonstrate the molecular adaptability which allows RNA to participate as a key player in these functions. A distinct repertoire of molecular motifs provides a basis for the assembly of complex RNA tertiary architectures.

Magnesium-dependent alternative foldings of active and inactive Escherichia coli tRNA(Glu) revealed by chemical probing

A stable conformer of Escherichia coli tRNA(Glu), obtained in the absence of Mg(2+), is inactive in the aminoacylation reaction. Probing it with diethylpyrocarbonate, dimethyl sulfate and ribonuclease V1 revealed that it has a hairpin structure with two internal loops; the helical segments at both extremities have the same structure as the acceptor stem and the anticodon arm of the native conformer of tRNA(Glu)and the middle helix is formed of nucleotides from the D-loop (G15-C20:2) and parts of the T-loop and stem (G51-C56), with G19 bulging out. This model is consistent with other known properties of this inactive conformer, including its capacity to dimerize. Therefore, this tRNA requires magnesium to acquire a conformation that can be aminoacylated, as others require a post-transcriptional modification to reach this active conformation.

A new method to monitor the rate of conformational transitions in RNA

Many RNAs need Mg2+to produce stable tertiary structures. Here we describe a simple method to measure the rate and activation parameters of tertiary structure unfolding that exploits this Mg2+dependence. Our approach is based on mixing an RNA solution with excess EDTA in a stopped-flow instrument equipped with an absorbance detector, under conditions of temperature and ionic strength where, after chelation of Mg2+, tertiary structure unfolds. We have demonstrated the utility of this method by studying phenylalanine-specific transfer RNA from yeast (tRNAPhe) because the unfolding rates and the corresponding activation parameters have been determined previously and provide a benchmark for our technique. We find that within error, our stopped-flow method reproduces both the rate and activation enthalpy for tertiary unfolding of yeast tRNAPhe measured previously by temperature-jump relaxation kinetics. Since many different RNAs require divalent magnesium for tertiary structure stabilization, this technique should be applicable to study the folding of other RNAs.

Hierarchy and dynamics of RNA folding

The evidence showing that the self-assembly of complex RNAs occurs in discrete transitions, each relating to the folding of sub-systems of increasing size and complexity starting from a state with most of the secondary structure, is reviewed. The reciprocal influence of the concentration of magnesium ions and nucleotide mutations on tertiary structure is analyzed. Several observations demonstrate that detrimental mutations can be rescued by high magnesium concentrations, while stabilizing mutations lead to a lesser dependence on magnesium ion concentration. Recent data point to the central controlling and monitoring roles of RNA-binding proteins that can bind to the different folding stages, either before full establishment of the secondary structure or at the molten globule state before the cooperative transition to the final three-dimensional structure.

RNA structure and the regulation of gene expression

RNA secondary and tertiary structure is involved in post-transcriptional regulation of gene expression either by exposing specific sequences or through the formation of specific structural motifs. An overview of RNA secondary and tertiary structures known from biophysical studies is followed by a review of examples of the elements of RNA processing, mRNA stability and translation of the messenger. These structural elements comprise sense-antisense double-stranded RNA, hairpin and stem-loop structures, and more complex structures such as bifurcations, pseudoknots and triple-helical elements. Metastable structures formed during RNA folding pathway are also discussed. The examples presented are mostly chosen from plant systems, plant viruses, and viroids. Examples from bacteria or fungi are discussed only when unique regulatory properties of RNA structures have been elucidated in these systems.

The first solvation shell of magnesium and calcium ions in a model nucleic acid environment: an ab initio study

The interaction of organophosphate anions with divalent metal ions is central to many biological catalytic events. While experimental structural studies can give insight into the likely geometries that can be adopted, quantum mechanics allows for a more complete exploration of the competing forms. Ab initio quantum mechanical calculations have been performed on a series of complexes comprised of dimethyl phosphate, a divalent metal ion (either Mg(II) or Ca(II)) and water of hydration. An additional series of complexes were studied that included a Cl(I) ion to provide for charge neutrality. The most stable orientation of the hydrated metal ion complexed with the phosphate anion occurs when the metal ion is in a unidentate, rather than bidentate, orientation. The question of whether the divalent metal ion is located in the phosphinyl (-PO2(-)-) plane depends on the identity of the divalent metal ion and on the charge state of the complex.

Modeling study on the cleavage step of the self-splicing reaction in group I introns

A three-dimensional model of the Tetrahymena thermophila group I intron is used to further explore the catalytic mechanism of the transphosphorylation reaction of the cleavage step. Based on the coordinates of the catalytic core model proposed by Michel and Westhof (Michel, F., Westhof, E. J. Mol. Biol. 216, 585-610 (1990)), we first converted their ligation step model into a model of the cleavage step by the substitution of several bases and the removal of helix P9. Next, an attempt to place a trigonal bipyramidal transition state model in the active site revealed that this modified model for the cleavage step could not accommodate the transition state due to insufficient space. A lowering of P1 helix relative to surrounding helices provided the additional space required. Simultaneously, it provided a better starting geometry to model the molecular contacts proposed by Pyle et al. (Pyle, A. M., Murphy, F. L., Cech, T. R. Nature 358, 123-128. (1992)), based on mutational studies involving the J8/7 segment. Two hydrated Mg2+ complexes were placed in the active site of the ribozyme model, using the crystal structure of the functionally similar Klenow fragment (Beese, L.S., Steitz, T.A. EMBO J. 10, 25-33 (1991)) as a guide. The presence of two metal ions in the active site of the intron differs from previous models, which incorporate one metal ion in the catalytic site to fulfill the postulated roles of Mg2+ in catalysis. The reaction profile is simulated based on a trigonal bipyramidal transition state, and the role of the hydrated Mg2+ complexes in catalysis is further explored using molecular orbital calculations.

Visualizing the higher order folding of a catalytic RNA molecule

The higher order folding process of the catalytic RNA derived from the self-splicing intron of Tetrahymena thermophila was monitored with the use of Fe(II)-EDTA-induced free radical chemistry. The overall tertiary structure of the RNA molecule forms cooperatively with the uptake of at least three magnesium ions. Local folding transitions display different metal ion dependencies, suggesting that the RNA tertiary structure assembles through a specific folding intermediate before the catalytic core is formed. Enzymatic activity, assayed with an RNA substrate that is complementary to the catalytic RNA active site, coincides with the cooperative structural transition. The higher order RNA foldings produced by Mg(II), Ca(II), and Sr(II) are similar; however, only the Mg(II)-stabilized RNA is catalytically active. Thus, these results directly demonstrate that divalent metal ions participate in general folding of the ribozyme tertiary structure, and further indicate a more specific involvement of Mg(II) in catalysis.

Structural biology of zinc

The biological function of zinc is governed by the composition of its tetrahedral coordination polyhedron in the metalloprotein, and each ligand group that coordinates to the metal ion does so with a well-defined stereochemical preference. Consequently, protein-zinc recognition and discrimination requires proper chemical composition and proper stereochemistry of the metal-ligand environment. However, it should be noted that the entire protein behaves as the "zinc ligand," since residues that are quite distant from the metal affect recognition and function by through-space (either solvent or the protein milieu) or through-hydrogen bond coulombic interactions. Additionally, long-range interactions across hydrogen bonds serve to orient ligands and therefore minimize the entropy loss incurred on metal binding. Since zinc is not subject to ligand field stabilization effects, it is easy for the tetrahedral protein-binding site to discriminate zinc from other first-row transition metal ions: It is only for Zn2+ that the change from an octahedral to a tetrahedral ligand field is not energetically disfavored. Structural considerations such as these must illuminate the engineering of de novo zinc-binding sites in proteins. Zinc serves chemical, structural, and regulatory roles in biological systems. In biological chemistry zinc serves as an electrophilic catalyst; that is, it stabilizes negative charges encountered during an enzyme-catalyzed reaction. The coordination polyhedron of catalytic zinc is usually dominated by histidine side chains. In biological structure zinc is typically sequestered from solvent, and its coordination polyhedron is almost exclusively dominated by cysteine thiolates. Structural or regulatory zinc is found as either a single metal ion or as part of a cluster of two or more metals. In multinuclear clusters cysteine thiolates either bridge two metal ions or serve as terminal ligands to a single metal ion. Even in complex multinuclear clusters, Zn2+ displays tetrahedral coordination. The structural biology of zinc continues to receive attention in catalytic and regulatory systems such as leucine aminopeptidase, alkaline phosphatase, transcription factors, and steroid receptors. For example, zinc-mediated hormone-receptor association has recently been demonstrated in the binding of human growth hormone to the extracellular binding domain of the human prolactin receptor (Cunningham et al., 1990). To be sure, structural studies of zinc in biology will continue to be a fruitful source of bioinorganic advances, as well as surprises, in the future.

Defining the inside and outside of a catalytic RNA molecule

Ribozymes are RNA molecules that catalyze biochemical reactions. Fe(II)-EDTA, a solvent-based reagent which cleaves both double- and single-stranded RNA, was used to investigate the structure of the Tetrahymena ribozyme. Regions of cleavage alternate with regions of substantial protection along the entire RNA molecule. In particular, most of the catalytic core shows greatly reduced cleavage. These data constitute experimental evidence that an RNA enzyme, like a protein enzyme, has an interior and an exterior. Determination of positions where the phosphodiester backbone of the RNA is on the inside or on the outside of the molecule provides major constraints for modeling the three-dimensional structure of the Tetrahymena ribozyme. This approach should be generally informative for structured RNA molecules.

Secondary structure of Bombyx mori and Dictyostelium discoideum 5S rRNA from S1 nuclease and cobra venom ribonuclease susceptibility, and computer assisted analysis

The 5S rRNAs from Bombyx mori and Dictyostelium discoideum were end-labeled with [32-P] at either the 5' or 3' end and sequenced using enzymatic digestion. The secondary structure of these molecules was studied using the single-strand specific S1 nuclease and the base-pair specific cobra venom ribonuclease. Computer analysis of these results was performed and was used to generate a consensus secondary structure for each molecule. A comparison of these results with those of other workers is presented.

Spin-labelled nucleic acids

The spin label method developed by McConnell 15 years ago is now widely used in studies of the structure and dynamic properties of a variety of the biological systems such as proteins and protein complexes, lipids and membranes, nucleic acids, nucleoproteins, etc.

The ESR spectrum of the nitroxide radcal – the spin label – is very sensitive to its microenvironment and permits easy registration of even subtle alterations in it. If spin labels are attached to different sites of a macromolecule the information can be gained about conformational properties of all these local regions and, as a result, about the dynamic behaviour of the object as a whole.

The interaction of polyamines with DNA: a 23Na NMR study

The interaction between a variety of polyamines, both naturally occurring and synthetic, and calf thymus DNA has been studied using 23Na NMR. The relaxation behaviour of 23Na reflects the extent of interaction of Na+ with DNA phosphate groups and therefore the extent of charge neutralisation of DNA phosphate groups (P) by polyamine amino and imino groups (N) in solutions of DNa, polyamine and Na+. The studies reveal that whereas spermine and spermidine are capable of expelling nearly all of the Na+ ions from DNA at N/P approximately 1, diamines such as putrescine and homologues of spermine and spermidine are capable of neutralising only roughly 50% of DNA phosphates. The results provide a challenge to current models of DNA-polyamine interactions.

Chemical probes for higher-order structure in RNA

Three chemical reactions can probe the secondary and tertiary interactions of RNA molecules in solution. Dimethyl sulfate monitors the N-7 of guanosines and senses tertiary interactions there, diethyl pyrocarbonate detects stacking of adenosines, and an alternate dimethyl sulfate reaction examines the N-3 of cytidines and thus probes base pairing. The reactions work between 0 degrees C and 90 degrees C and at pH 4.5--8.5 in a variety of buffers. As an example we follow the progressive denaturation of yeast tRNAPhe terminally labeled with 32P as the tertiary and secondary structures sequentially melt out. A single autoradiograph of a terminally labeled molecule locates regions of higher-order structure and identifies the bases involved.

Conformational states of cowpea chlorotic mottle virus ribonucleic acid components

The conditions determining conformational changes in the four ribonucleic acid components of cowpea chlorotic mottle virus have been studied. All four components have at least two electrophoretically separable conformers, the occurrence of which can be regulated by both monovalent and polyvalent cations. This phenomenon also occurs, in a much less striking way, in the ribonucleic acids of the two other members of the bromovirus group, brome mosaic virus and broad bean mottle virus. Although specific in some respects, these changes have much in common with effects which have been observed in tRNAs, 55 RNAs and rRNAs. A provisional interpretation of the conformational behaviour of the viral RNAs is given in terms which have been proposed for certain tRNAs which have been studied in great detail.

On the binding of Mg2+ and Mn2+ to tRNA

In this paper previous binding studies of Mg2+ and Mn2+ ions to tRNA's are reconsidered. Binding data of some representative examples are interpreted including interactions between charges located on the macroion. Both curved and bell-shaped Scatchard plots can be accounted for quantitatively if corrections are made for electrostatic interactions and, if necessary, for the effect of conformational changes on these interactions. It appears that there is no need to invoke more than one class of binding sites on tRNA's, meaning that the experimental binding data can be described using the same intrinsic pK value for all phosphate groups.

RNA structure

This review is concerned primarily with the physical structure and changes in the structure of RNA molecules. It will be evident that we have not attempted comprehensive coverage of what amounts to a vast literature. We have tried to stay away from particular areas that have been recently reviewed elsewhere. Citations to and information from them are included, however, so that access to the literature is available. Much of what we treat in depth deals with the crystal structures and solution behaviour of model RNA compounds, including synthetic polymers and molecular fragments such as dinucleoside phosphates. Sequence data on natural RNA are cited, but not in detail. Similarly, apart from tRNA, natural RNAs the structural determinations of which are presently not so far advanced, are not dwelt upon. We have tried to present in detail the available structural data with scaled drawings that permit facile comparisons of molecular geometries.

Studies on secondary structure of single-stranded RNA from bacteriophage MS2 by electron microscopy

A method allowing the demonstration and study by electron microscopy of secondary structure of viral RNA has been developed. Single-stranded RNA from the bacteriophage MS2 has been analyzed in the electron microscope in the presence of various concentrations of MgCl2. Depending on the salt concentration, the molecules display one to three large open loops which range in size from 10 to 20% of the total RNA length, and smaller closed loops which are approximately 3-5% of the total RNA length. Within one spreading, the conformation of the molecules is variable. However, the average complexity of the molecules increases with increasing salt, and individual loops which are infrequent at low salt increase in frequency with increasing salt. By analyzing the manner in which the individual loop appeared, it was possible to show that all molecules could be described by one basic pattern of secondary structure formation.

Proton nuclear magnetic resonance of spin-labeled Escherichia coli tRNAf1MET

Thiouridine at position 8 (s4U8) of tRNAf1Met was spin-labeled with the nitroxide free radical, N-(1-oxyl-2,2,5,5-Tetramethyl-3-pyrrolidinyl) bromacetamide, for proton nuclear magnetic resonance spectroscopic studies. The well-resolved methyl peak of ribothymidine is unperturbed, but the peak tentatively assigned to the C-5 methylene group of dihydrouridine is considerably broadened in spin-labeled tRNAf1Met. Of the approximately 27 slowly exchanging protons observed in the region between 11 and 15 ppm downfield from 4,4-dimethyl-4-silapentane-1-sulfonic acid, the equivalent of about five protons apparently disappeared in spin-labeled tRNAf1Met. The well-resolved single proton at 14.8 ppm was missing not only in the paramagnetic species, but also in the diamagnetic reduced form of spin-labeled tRNAf1Met, and was unequivocally identified as a hydrogen bond involving s4U8 by comparison of several forms of tRNAf1Met specifically modified at s4U. Evidence that the perturbation of a second single proton resonance at 14.6 ppm (shift and broadening) is coupled to the loss of a tertiary hydrogen bond involving residue 8, arises from the same modified forms. The resolved resonances in the methyl and N-H regions, particularly the resonance at 14.6 ppm as well as the four N-bonded proton resonances at higher field which are broadened solely due to their proximity to the unpaired electron of the spin label, provide specific indicators of the geometry of tRNAf1Met structure in solution. Their observability by nuclear magnetic resonance spectroscopy opens up the possibility of monitoring distance changes among the base residues of spin-labeled tRNAf1Met upon its interaction with aminoacyl-tRNA synthetase and other enzymes.

tRNA conformation and magnesium binding. A study of a yeast phenylalanine-specific tRNA by a fluorescent indicator and differential melting curves

The binding of Mg2+ to tRNAPhe (yeast) in three conformational states was studied at 10, 30, 45, and 70 degrees C by the fluorescence indicator 8-hydroxyquinoline 5-sulphonic acid in the presence of 0.032 M monovalent cations (Na+). At temperatures below those characteristic for early melting (completely folded tRNA) the Scatchard plots are biphasic. They were well fitted by two classes of noninteracting binding sites with stability constants independent of temperature (KA = 9X10(4), KB = 6X10(3) M-1). In partially unfolded tRNA the strong binding process is co-operative. A single class of weak sites was found in the statistically coiled conformation at 70 degrees C (KB = 3.3X103 M-1). The total number of binding sites is 23 +/- 5; differences for the folded and unfolded conformations are smaller than 1. The influence of Mg2+ on the stability of the conformational elements of tRNAPhe (yeast) and its CCA-half (i.e. nucleotides 38--76) was determined by differential ultraviolet absorbance and depolarisation melting curves using the fluorescence of the Y base. Tertiary structure corresponding to early melting is stabilized by strongly bound Mg2+, whereas all other melting transitions are only influenced by Mg2+ bound at weak sites. The stability constants of tertiary structure obtained from the melting experiments can quantitatively be described by assuming that 5 +/- 1 non-interacting strong sites as characterized by the fluorescence titrations are converted to weak sites upon unfolding of the tertiary structure. Co-operative interaction of Mg2+ with the 5 strong sites in the folded conformation of tRNA can be ruled out. Strong binding of Mg2+ to completely folded tRNA does not produce a conformational transition changing ultraviolet absorbance, circular dichroism and sedimentation coefficient.

Method for predicting RNA secondary structure

We report a method for predicting the most stable secondary structure of RNA from its primary sequence of nucleotides. The technique consists of a series of three computer programs interfaced to take the nucleotide sequence of any RNA and (a) list all possible helical regions, using modified Watson-Crick base-pairing rules; (b) create all possible secondary structures by forming permutations of compatible helical regions; and (c)evaluate each structure for total free energy of formation from a completely extended chain. A free energy distribution and the base-by-base bonding interactions of each possible structure are catalogued by the system and are readily available for examination. The method has been applied to 62 tRNA sequences. The total free-energy of the predicted most stable structures ranged from -19 to -41 kcal/mole (-22 to -49 kJ/mole). The number of structures created was also highly sequence-dependent and ranged from 200 to 13,000. In nearly all cases the cloverleaf is predicted to be the structure with the lowest free energy of formation.