SwS: a solvation web service for nucleic acids

Knowledge-based prediction of DNA hydration using hydrated dinucleotides as building blocks

Database-derived water probability densities around structurally and sequentially distinct DNA dinucleotide fragments reproduce the known hydration motifs, which thus can be used as building blocks to predict DNA hydration.

Water plays an important role in stabilizing the structure of DNA and mediating its interactions. Here, the hydration of DNA was analyzed in terms of dinucleotide fragments from an ensemble of 2727 nonredundant DNA chains containing 41 853 dinucleotides and 316 265 associated first-shell water molecules. The dinucleotides were classified into categories based on their 16 sequences and the previously determined structural classes known as nucleotide conformers (NtCs). The construction of hydrated dinucleotide building blocks allowed dinucleotide hydration to be calculated as the probability of water density distributions. Peaks in the water densities, known as hydration sites (HSs), uncovered the interplay between base and sugar-phosphate hydration in the context of sequence and structure. To demonstrate the predictive power of hydrated DNA building blocks, they were then used to predict hydration in an independent set of crystal and NMR structures. In ten tested crystal structures, the positions of predicted HSs and experimental waters were in good agreement (more than 40% were within 0.5 Å) and correctly reproduced the known features of DNA hydration, for example the ‘spine of hydration’ in B-DNA. Therefore, it is proposed that hydrated building blocks can be used to predict DNA hydration in structures solved by NMR and cryo-EM, thus providing a guide to the interpretation of experimental data and computer models. The data for the hydrated building blocks and the predictions are available for browsing and visualization at the website https://watlas.datmos.org/watna/.

Role of Water in Defining the Structure and Properties of B-Form DNA

DNA in the cell is rarely naked but normally protein-bound in nucleosomes. Of special interest is the DNA bound to other factors that control its key functions of transcription, replication, and repair. For these several transactions of DNA, the state of hydration plays an important role in its function, and therefore needs to be defined in as much detail as possible. High-resolution crystallography of short B-form duplexes shows that the mixed polar and apolar surface of the major groove binds water molecules over the broad polar floor of the groove in a sequence-dependent varied manner. In contrast, the narrower minor groove, particularly at AT-rich segments, binds water molecules to the polar groups of the bases in a regular double layer reminiscent of the structure of ice. This review is largely devoted to measurements made in solution, principally calorimetric, that are fully consistent with the location of water molecules seen in crystals, thereby emphasizing the substantial difference between the hydration patterns of the two grooves.

Structural Bases for the Fitness Cost of the Antibiotic-Resistance and Lethal Mutations at Position 1408 of 16S rRNA

To understand a structural basis for the fitness cost of the A1408G antibiotic-resistance mutation in the ribosomal A-site RNA, we have determined crystal structures of its A1408C and A1408U lethal mutants, and made comparison with previously solved structures of the wild type and the antibiotic-resistant mutant. The A-site RNA containing an asymmetric internal loop functions as a molecular switch to discriminate a single cognate tRNA from several near-cognate tRNAs by its conformational ON/OFF switching. Overall structures of the “off” states of the A1408C/U lethal mutants are very similar to those of the wild type and the A1408G antibiotic-resistant mutant. However, significant differences are found in local base stacking interactions including the functionally important A1492 and A1493 residues. In the wild type and the A1408G antibiotic-resistant mutant “off” states, both adenines are exposed to the solvent region. On the other hand, one of the corresponding adenines of the lethal A1408C/U mutants stay deeply inside their A-site helices by forming a purine-pyrimidine AoC or A-U base pair and is sandwiched between the upper and lower bases. Therefore, the ON/OFF switching might unfavorably occur in the lethal mutants compared to the wild type and the A1408G antibiotic-resistant mutant. It is probable that bacteria manage to acquire antibiotic resistance without losing the function of the A-site molecular switch by mutating the position 1408 only from A to G, but not to pyrimidine base C or U.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

WatAA: Atlas of Protein Hydration. Exploring synergies between data mining and ab initio calculations

Water molecules represent an integral part of proteins and a key determinant of protein structure, dynamics and function. WatAA is a newly developed, web-based atlas of amino-acid hydration in proteins. The atlas provides information about the ordered first hydration shell of the most populated amino-acid conformers in proteins. The data presented in the atlas are drawn from two sources: experimental data and ab initio quantum-mechanics calculations. The experimental part is based on a data-mining study of a large set of high-resolution protein crystal structures. The crystal-derived data include 3D maps of water distribution around amino-acids and probability of occurrence of each of the identified hydration sites. The quantum mechanics calculations validate and extend this primary description by optimizing the water position for each hydration site, by providing hydrogen atom positions and by quantifying the interaction energy that stabilizes the water molecule at the particular hydration site position. The calculations show that the majority of experimentally derived hydration sites are positioned near local energy minima for water, and the calculated interaction energies help to assess the preference of water for the individual hydration sites. We propose that the atlas can be used to validate water placement in electron density maps in crystallographic refinement, to locate water molecules mediating protein-ligand interactions in drug design, and to prepare and evaluate molecular dynamics simulations. WatAA: Atlas of Protein Hydration is freely available without login at .

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.

Structure of the ordered hydration of amino acids in proteins: analysis of crystal structures

Crystallography provides unique information about the arrangement of water molecules near protein surfaces. Using a nonredundant set of 2818 protein crystal structures with a resolution of better than 1.8 Å, the extent and structure of the hydration shell of all 20 standard amino-acid residues were analyzed as function of the residue conformation, secondary structure and solvent accessibility. The results show how hydration depends on the amino-acid conformation and the environment in which it occurs. After conformational clustering of individual residues, the density distribution of water molecules was compiled and the preferred hydration sites were determined as maxima in the pseudo-electron-density representation of water distributions. Many hydration sites interact with both main-chain and side-chain amino-acid atoms, and several occurrences of hydration sites with less canonical contacts, such as carbon-donor hydrogen bonds, OH-π interactions and off-plane interactions with aromatic heteroatoms, are also reported. Information about the location and relative importance of the empirically determined preferred hydration sites in proteins has applications in improving the current methods of hydration-site prediction in molecular replacement, ab initio protein structure prediction and the set-up of molecular-dynamics simulations.

Fusion RNAs in crystallographic studies of double-stranded RNA from trypanosome RNA editing

Head-to-head fusions of two identical double-stranded fragments of RNA can be designed to self-assemble from a single RNA species and form a double-stranded helix with a twofold rotation axis relating the two strands. These symmetrical RNA molecules are more likely to crystallize without end-on-end statistical packing disorder because the two halves of the molecule are identical. This approach can be used to study many fragments of double-stranded RNA or many isolated helical domains from large single-stranded RNAs that may not yet be amenable to high-resolution studies by crystallography or NMR. We used fusion RNAs to study one (the U-helix) of three functional domains formed when guide RNA binds to its cognate pre-edited mRNA from the trypanosome RNA editing system. The U-helix forms when the 3' oligo(U) tail of the guide RNA (gRNA) binds to the purine-rich, pre-edited mRNA upstream from the current RNA editing site. Fusion RNAs 16-and 32-base pairs in length formed crystals that gave diffraction to 1.37 and 1.05 Å respectively. We provide the composition of a fusion RNA crystallization screen and describe the X-ray data collection, structure determination, and refinement of the crystal structures of fusion RNAs.

On using magnesium and potassium ions in RNA experiments

RNA is a dynamic molecule that can adopt different conformations under different conditions. In the cell, the negatively charged RNA is expected to interact with potassium and magnesium ions. The role of magnesium ions has been extensively studied in RNA, as it is necessary for proper folding of RNA and acts as a nucleophile in many reactions catalyzed by RNA. A short overview of the role of magnesium and potassium in RNA structures is provided here in order to demonstrate the need for carefully choosing the experimental conditions for RNA studies.

Structural studies of CNG repeats

CNG repeats (where N denotes one of the four natural nucleotides) are abundant in the human genome. Their tendency to undergo expansion can lead to hereditary diseases known as TREDs (trinucleotide repeat expansion disorders). The toxic factor can be protein, if the abnormal gene is expressed, or the gene transcript, or both. The gene transcripts have attracted much attention in the biomedical community, but their molecular structures have only recently been investigated. Model RNA molecules comprising CNG repeats fold into long hairpins whose stems generally conform to an A-type helix, in which the non-canonical N-N pairs are flanked by C-G and G-C pairs. Each homobasic pair is accommodated in the helical context in a unique manner, with consequences for the local helical parameters, solvent structure, electrostatic potential and potential to interact with ligands. The detailed three-dimensional profiles of RNA CNG repeats can be used in screening of compound libraries for potential therapeutics and in structure-based drug design. Here is a brief survey of the CNG structures published to date.

Structure of the myotonic dystrophy type 2 RNA and designed small molecules that reduce toxicity

Myotonic dystrophy type 2 (DM2) is an incurable neuromuscular disorder caused by a r(CCUG) expansion (r(CCUG)(exp)) that folds into an extended hairpin with periodically repeating 2×2 nucleotide internal loops (5'CCUG/3'GUCC). We designed multivalent compounds that improve DM2-associated defects using information about RNA-small molecule interactions. We also report the first crystal structure of r(CCUG) repeats refined to 2.35 Å. Structural analysis of the three 5'CCUG/3'GUCC repeat internal loops (L) reveals that the CU pairs in L1 are each stabilized by one hydrogen bond and a water-mediated hydrogen bond, while CU pairs in L2 and L3 are stabilized by two hydrogen bonds. Molecular dynamics (MD) simulations reveal that the CU pairs are dynamic and stabilized by Na(+) and water molecules. MD simulations of the binding of the small molecule to r(CCUG) repeats reveal that the lowest free energy binding mode occurs via the major groove, in which one C residue is unstacked and the cross-strand nucleotides are displaced. Moreover, we modeled the binding of our dimeric compound to two 5'CCUG/3'GUCC motifs, which shows that the scaffold on which the RNA-binding modules are displayed provides an optimal distance to span two adjacent loops.

Structural insights into how 5-hydroxymethylation influences transcription factor binding

Transcription factor binding and high resolution crystallographic studies (1.3 Å) of Dickerson-Drew duplexes with cytosine, methylcytosine and hydroxymethylcytosine bases provide evidence that C-5 cytosine modifications could regulate transcription by context dependent effects on DNA transcription factor interactions.

Analyzing and building nucleic acid structures with 3DNA

The 3DNA software package is a popular and versatile bioinformatics tool with capabilities to analyze, construct, and visualize three-dimensional nucleic acid structures. This article presents detailed protocols for a subset of new and popular features available in 3DNA, applicable to both individual structures and ensembles of related structures. Protocol 1 lists the set of instructions needed to download and install the software. This is followed, in Protocol 2, by the analysis of a nucleic acid structure, including the assignment of base pairs and the determination of rigid-body parameters that describe the structure and, in Protocol 3, by a description of the reconstruction of an atomic model of a structure from its rigid-body parameters. The most recent version of 3DNA, version 2.1, has new features for the analysis and manipulation of ensembles of structures, such as those deduced from nuclear magnetic resonance (NMR) measurements and molecular dynamic (MD) simulations; these features are presented in Protocols 4 and 5. In addition to the 3DNA stand-alone software package, the w3DNA web server, located at http://w3dna.rutgers.edu, provides a user-friendly interface to selected features of the software. Protocol 6 demonstrates a novel feature of the site for building models of long DNA molecules decorated with bound proteins at user-specified locations.

Studies of base pair sequence effects on DNA solvation based on all-atom molecular dynamics simulations

Detailed analyses of the sequence-dependent solvation and ion atmosphere of DNA are presented based on molecular dynamics (MD) simulations on all the 136 unique tetranucleotide steps obtained by the ABC consortium using the AMBER suite of programs. Significant sequence effects on solvation and ion localization were observed in these simulations. The results were compared to essentially all known experimental data on the subject. Proximity analysis was employed to highlight the sequence dependent differences in solvation and ion localization properties in the grooves of DNA. Comparison of the MD-calculated DNA structure with canonical A- and B-forms supports the idea that the G/C-rich sequences are closer to canonical A- than B-form structures, while the reverse is true for the poly A sequences, with the exception of the alternating ATAT sequence. Analysis of hydration density maps reveals that the flexibility of solute molecule has a significant effect on the nature of observed hydration. Energetic analysis of solute-solvent interactions based on proximity analysis of solvent reveals that the GC or CG base pairs interact more strongly with water molecules in the minor groove of DNA that the AT or TA base pairs, while the interactions of the AT or TA pairs in the major groove are stronger than those of the GC or CG pairs. Computation of solvent-accessible surface area of the nucleotide units in the simulated trajectories reveals that the similarity with results derived from analysis of a database of crystallographic structures is excellent. The MD trajectories tend to follow Manning's counterion condensation theory, presenting a region of condensed counterions within a radius of about 17 A from the DNA surface independent of sequence. The GC and CG pairs tend to associate with cations in the major groove of the DNA structure to a greater extent than the AT and TA pairs. Cation association is more frequent in the minor groove of AT than the GC pairs. In general, the observed water and ion atmosphere around the DNA sequences is the MD simulation is in good agreement with experimental observations.

Integrity of the P-site is probed during maturation of the 60S ribosomal subunit

The P-site of the 60S ribosomal subunit signals to Tif6 via Elf1 during ribosomal maturation, suggesting a quasifunctional check of the integrity of the 60S subunit before the first round of translation.

Eukaryotic ribosomes are preassembled in the nucleus and mature in the cytoplasm. Release of the antiassociation factor Tif6 by the translocase-like guanosine triphosphatase Efl1 is a critical late maturation step. In this paper, we show that a loop of Rpl10 that embraces the P-site transfer ribonucleic acid was required for release of Tif6, 90 Å away. Mutations in this P-site loop blocked 60S maturation but were suppressed by mutations in Tif6 or Efl1. Molecular dynamics simulations of the mutant Efl1 proteins suggest that they promote a conformation change in Efl1 equivalent to changes that elongation factor G and eEF2 undergo during translocation. These results identify molecular signaling from the P-site to Tif6 via Efl1, suggesting that the integrity of the P-site is interrogated during maturation. We propose that Efl1 promotes a functional check of the integrity of the 60S subunit before its first round of translation.

MINAS--a database of Metal Ions in Nucleic AcidS

Correctly folded into the respective native 3D structure, RNA and DNA are responsible for uncountable key functions in any viable organism. In order to exert their function, metal ion cofactors are closely involved in folding, structure formation and, e.g. in ribozymes, also the catalytic mechanism. The database MINAS, Metal Ions in Nucleic AcidS (http://www.minas.uzh.ch), compiles the detailed information on innersphere, outersphere and larger coordination environment of >70 000 metal ions of 36 elements found in >2000 structures of nucleic acids contained today in the PDB and NDB. MINAS is updated monthly with new structures and offers a multitude of search functions, e.g. the kind of metal ion, metal-ligand distance, innersphere and outersphere ligands defined by element or functional group, residue, experimental method, as well as PDB entry-related information. The results of each search can be saved individually for later use with so-called miniPDB files containing the respective metal ion together with the coordination environment within a 15 Å radius. MINAS thus offers a unique way to explore the coordination geometries and ligands of metal ions together with the respective binding pockets in nucleic acids.

Hydration sites of unpaired RNA bases: a statistical analysis of the PDB structures

BACKGROUND

Hydration is crucial for RNA structure and function. X-ray crystallography is the most commonly used method to determine RNA structures and hydration and, therefore, statistical surveys are based on crystallographic results, the number of which is quickly increasing.

RESULTS

A statistical analysis of the water molecule distribution in high-resolution X-ray structures of unpaired RNA nucleotides showed that: different bases have the same penchant to be surrounded by water molecules; clusters of water molecules indicate possible hydration sites, which, in some cases, match those of the major and minor grooves of RNA and DNA double helices; complex hydrogen bond networks characterize the solvation of the nucleotides, resulting in a significant rigidity of the base and its surrounding water molecules. Interestingly, the hydration sites around unpaired RNA bases do not match, in general, the positions that are occupied by the second nucleotide when the base-pair is formed.

CONCLUSIONS

The hydration sites around unpaired RNA bases were found. They do not replicate the atom positions of complementary bases in the Watson-Crick pairs.

Experimental and computational determination of tRNA dynamics

As the molecular representation of the genetic code, tRNA plays a central role in the translational machinery where it interacts with several proteins and other RNAs during the course of protein synthesis. These interactions exploit the dynamic flexibility of tRNA. In this minireview, we discuss the effects of modified bases, ions, and proteins on tRNA structure and dynamics and the challenges of observing its motions over the cycle of translation.

Hydrogen-bonding patterns of minor groove-binder-DNA complexes reveal criteria for discovery of new scaffolds

Minor groove-binding ligands are able to control gene expression and are of great interest for therapeutic applications. We extracted hydrogen-bonding geometries from all available structures of minor groove-binder-DNA complexes of two noncovalent binding modes, namely 1:1 (including hairpin and cyclic ligands) and 2:1 ligand/DNA binding. Positions of the ligand atoms involved in hydrogen bonding deviate from idealized hydrogen bond geometries and do not exploit the possibilities indicated by water molecules. Therefore, we suggest the inclusion of shape-based descriptors rather than hydrogen-bond patterns in virtual screening protocols for the identification of innovative minor groove-binding scaffolds.

Structural insights into CUG repeats containing the 'stretched U-U wobble': implications for myotonic dystrophy

Tracks containing CUG repeats are abundant in human gene transcripts. Their biological role includes modulation of pre-mRNA splicing, mRNA transport and regulation of translation. Expanded forms of CUG runs are associated with pathogenesis of several neurodegenerative diseases, including myotonic dystrophy type 1. We have analysed two crystal structures of RNA duplexes containing the CUG repeats: G(CUG)₂C and (CUG)₆. The first of the structures, analysed at 1.23 Å resolution, is of an oligomer designed by us. The second model was obtained after ‘detwinning’ the 1.58 Å X-ray data previously deposited in the PDB. The RNA duplexes are in the A-form in which all the C–G pairs form Watson–Crick interactions while all the uridine pairs can be described as U•U cis wobble having only one hydrogen bond between the bases. The residue, which accepts the H-bond, is inclined towards the minor groove. This previously unreported base pairing can be described as ‘stretched U–U wobble’. The regular hydrogen-bonding pattern of interactions with the solvent, the electrostatic charge distribution and surface features indicate the ligand binding potential of the CUG tracks.

Noncanonical G(syn)-G(anti) base pairs stabilized by sulphate anions in two X-ray structures of the (GUGGUCUGAUGAGGCC) RNA duplex

The structures of two crystal forms of the RNA 16-mer with the sequence GUGGUCUGAUGAGGCC, grown in the presence of a high concentration of sulphate ions, have been determined using synchrotron radiation at 1.4- and 2.0-A resolution. RNA with this sequence is known as one of the two strands of the noncleavable form of the hammerhead ribozyme. In both crystal structures, two G(syn)-G(anti) noncanonical base pairs are observed in the middle of a 14 base-pair (bp) duplex having 5'-dangling GU residues. Both structures contain sulphate anions interacting with the G-G bp stabilizing G in its syn conformation and bridging the two RNA strands. In both cases the interactions take place in the major groove, although the anions are accommodated within different helix geometries, most pronounced in the changing width of the major groove. In one structure, where a single sulphate spans both G-G pairs, the major groove is closed around the anion, while in the other structure, where each of the two G-G pairs is associated with a separate sulphate, the groove is open. This work provides the first examples of a G-G pair in syn-anti conformation, which minimizes the purine-purine clash in the center of the duplex, while utilizing its residual hydrogen bonding potential in specific interactions with sulphate anions.

3DNA: a versatile, integrated software system for the analysis, rebuilding and visualization of three-dimensional nucleic-acid structures

We present a set of protocols showing how to use the 3DNA suite of programs to analyze, rebuild and visualize three-dimensional nucleic-acid structures. The software determines a wide range of conformational parameters, including the identities and rigid-body parameters of interacting bases and base-pair steps, the nucleotides comprising helical fragments, the area of overlap of stacked bases and so on. The reconstruction of three-dimensional structure takes advantage of rigorously defined rigid-body parameters, producing rectangular block representations of the nucleic-acid bases and base pairs and all-atom models with approximate sugar-phosphate backbones. The visualization components create vector-based drawings and scenes that can be rendered as raster-graphics images, allowing for easy generation of publication-quality figures. The utility programs use geometric variables to control the view and scale of an object, for comparison of related structures. The commands run in seconds even for large structures. The software and related information are available at http://3dna.rutgers.edu/.

Tertiary structure and function of an RNA motif required for plant vascular entry to initiate systemic trafficking

Vascular entry is a decisive step for the initiation of long-distance movement of infectious and endogenous RNAs, silencing signals and developmental/defense signals in plants. However, the mechanisms remain poorly understood. We used Potato spindle tuber viroid (PSTVd) as a model to investigate the direct role of the RNA itself in vascular entry. We report here the identification of an RNA motif that is required for PSTVd to traffic from nonvascular into the vascular tissue phloem to initiate systemic infection. This motif consists of nucleotides U/C that form a water-inserted cis Watson-Crick/Watson-Crick base pair flanked by short helices that comprise canonical Watson-Crick/Watson-Crick base pairs. This tertiary structural model was inferred by comparison with X-ray crystal structures of similar motifs in rRNAs and is supported by combined mutagenesis and covariation analyses. Hydration pattern analysis suggests that water insertion induces a widened minor groove conducive to protein and/or RNA interactions. Our model and approaches have broad implications to investigate the RNA structural motifs in other RNAs for vascular entry and to study the basic principles of RNA structure-function relationships.

Nucleic acid solvation: from outside to insight

Nucleic acids are polyanionic molecules that were historically considered to be solely surrounded by a shell of water molecules and a neutralizing cloud of monovalent and divalent cations. In this respect, recent experimental and theoretical reports demonstrate that water molecules within complex nucleic acid structures can display very long residency times, and assist drug binding and catalytic reactions. Finally, anions can also bind to these polyanionic systems. Many of these recent insights are provided by state-of-the-art molecular dynamics simulations of nucleic acid systems, which will be described together with relevant methodological issues.