DSSR-enabled innovative schematics of 3D nucleic acid structures with PyMOL.

Sophisticated analysis and simplified visualization are crucial for understanding complicated structures of biomacromolecules. DSSR (Dissecting the Spatial Structure of RNA) is an integrated computational tool that has streamlined the analysis and annotation of 3D nucleic acid structures. The program creates schematic block representations in diverse styles that can be seamlessly integrated into PyMOL and complement its other popular visualization options. In addition to portraying individual base blocks, DSSR can draw Watson-Crick pairs as long blocks and highlight the minor-groove edges. Notably, DSSR can dramatically simplify the depiction of G-quadruplexes by automatically detecting G-tetrads and treating them as large square blocks. The DSSR-enabled innovative schematics with PyMOL are aesthetically pleasing and highly informative: the base identity, pairing geometry, stacking interactions, double-helical stems, and G-quadruplexes are immediately obvious. These features can be accessed via four interfaces: the command-line interface, the DSSR plugin for PyMOL, the web application, and the web application programming interface. The supplemental PDF serves as a practical guide, with complete and reproducible examples. Thus, even beginners or occasional users can get started quickly, especially via the web application at http://skmatic.x3dna.org.

Web 3DNA 2.0 for the analysis, visualization, and modeling of 3D nucleic acid structures.

Web 3DNA (w3DNA) 2.0 is a significantly enhanced version of the widely used w3DNA server for the analysis, visualization, and modeling of 3D nucleic-acid-containing structures. Since its initial release in 2009, the w3DNA server has continuously served the community by making commonly-used features of the 3DNA suite of command-line programs readily accessible. However, due to the lack of updates, w3DNA has clearly shown its age in terms of modern web technologies and it has long lagged behind further developments of 3DNA per se. The w3DNA 2.0 server presented here overcomes all known shortcomings of w3DNA while maintaining its battle-tested characteristics. Technically, w3DNA 2.0 implements a simple and intuitive interface (with sensible defaults) for increased usability, and it complies with HTML5 web standards for broad accessibility. Featurewise, w3DNA 2.0 employs the most recent version of 3DNA, enhanced with many new functionalities, including: the automatic handling of modified nucleotides; a set of 'simple' base-pair and step parameters for qualitative characterization of non-Watson-Crick double-helical structures; new structural parameters that integrate the rigid base plane and the backbone phosphate group, the two nucleic acid components most reliably determined with X-ray crystallography; in silico base mutations that preserve the backbone geometry; and a notably improved module for building models of single-stranded RNA, double-helical DNA, Pauling triplex, G-quadruplex, or DNA structures 'decorated' with proteins. The w3DNA 2.0 server is freely available, without registration, at http://web.x3dna.org.

DNA Conformational Changes Play a Force-Generating Role during Bacteriophage Genome Packaging.

Motors that move DNA, or that move along DNA, play essential roles in DNA replication, transcription, recombination, and chromosome segregation. The mechanisms by which these DNA translocases operate remain largely unknown. Some double-stranded DNA (dsDNA) viruses use an ATP-dependent motor to drive DNA into preformed capsids. These include several human pathogens as well as dsDNA bacteriophages-viruses that infect bacteria. We previously proposed that DNA is not a passive substrate of bacteriophage packaging motors but is instead an active component of the machinery. We carried out computational studies on dsDNA in the channels of viral portal proteins, and they reveal DNA conformational changes consistent with that hypothesis. dsDNA becomes longer ("stretched") in regions of high negative electrostatic potential and shorter ("scrunched") in regions of high positive potential. These results suggest a mechanism that electrostatically couples the energy released by ATP hydrolysis to DNA translocation: The chemical cycle of ATP binding, hydrolysis, and product release drives a cycle of protein conformational changes. This produces changes in the electrostatic potential in the channel through the portal, and these drive cyclic changes in the length of dsDNA as the phosphate groups respond to the protein's electrostatic potential. The DNA motions are captured by a coordinated protein-DNA grip-and-release cycle to produce DNA translocation. In short, the ATPase, portal, and dsDNA work synergistically to promote genome packaging.

Effects of Noncanonical Base Pairing on RNA Folding: Structural Context and Spatial Arrangements of G·A Pairs.

Noncanonical base pairs play important roles in assembling the three-dimensional structures critical to the diverse functions of RNA. These associations contribute to the looped segments that intersperse the canonical double-helical elements within folded, globular RNA molecules. They stitch together various structural elements, serve as recognition elements for other molecules, and act as sites of intrinsic stiffness or deformability. This work takes advantage of new software (DSSR) designed to streamline the analysis and annotation of RNA three-dimensional structures. The multiscale structural information gathered for individual molecules, combined with the growing number of unique, well-resolved RNA structures, makes it possible to examine the collective features deeply and to uncover previously unrecognized patterns of chain organization. Here we focus on a subset of noncanonical base pairs involving guanine and adenine and the links between their modes of association, secondary structural context, and contributions to tertiary folding. The rigorous descriptions of base-pair geometry that we employ facilitate characterization of recurrent geometric motifs and the structural settings in which these arrangements occur. Moreover, the numerical parameters hint at the natural motions of the interacting bases and the pathways likely to connect different spatial forms. We draw attention to higher-order multiplexes involving two or more G·A pairs and the roles these associations appear to play in bridging different secondary structural units. The collective data reveal pairing propensities in base organization, secondary structural context, and deformability and serve as a starting point for further multiscale investigations and/or simulations of RNA folding.

DSSR-enhanced visualization of nucleic acid structures in Jmol.

Sophisticated and interactive visualizations are essential for making sense of the intricate 3D structures of macromolecules. For proteins, secondary structural components are routinely featured in molecular graphics visualizations. However, the field of RNA structural bioinformatics is still lagging behind; for example, current molecular graphics tools lack built-in support even for base pairs, double helices, or hairpin loops. DSSR (Dissecting the Spatial Structure of RNA) is an integrated and automated command-line tool for the analysis and annotation of RNA tertiary structures. It calculates a comprehensive and unique set of features for characterizing RNA, as well as DNA structures. Jmol is a widely used, open-source Java viewer for 3D structures, with a powerful scripting language. JSmol, its reincarnation based on native JavaScript, has a predominant position in the post Java-applet era for web-based visualization of molecular structures. The DSSR-Jmol integration presented here makes salient features of DSSR readily accessible, either via the Java-based Jmol application itself, or its HTML5-based equivalent, JSmol. The DSSR web service accepts 3D coordinate files (in mmCIF or PDB format) initiated from a Jmol or JSmol session and returns DSSR-derived structural features in JSON format. This seamless combination of DSSR and Jmol/JSmol brings the molecular graphics of 3D RNA structures to a similar level as that for proteins, and enables a much deeper analysis of structural characteristics. It fills a gap in RNA structural bioinformatics, and is freely accessible (via the Jmol application or the JSmol-based website http://jmol.x3dna.org).

DSSR: an integrated software tool for dissecting the spatial structure of RNA.

Insight into the three-dimensional architecture of RNA is essential for understanding its cellular functions. However, even the classic transfer RNA structure contains features that are overlooked by existing bioinformatics tools. Here we present DSSR (Dissecting the Spatial Structure of RNA), an integrated and automated tool for analyzing and annotating RNA tertiary structures. The software identifies canonical and noncanonical base pairs, including those with modified nucleotides, in any tautomeric or protonation state. DSSR detects higher-order coplanar base associations, termed multiplets. It finds arrays of stacked pairs, classifies them by base-pair identity and backbone connectivity, and distinguishes a stem of covalently connected canonical pairs from a helix of stacked pairs of arbitrary type/linkage. DSSR identifies coaxial stacking of multiple stems within a single helix and lists isolated canonical pairs that lie outside of a stem. The program characterizes 'closed' loops of various types (hairpin, bulge, internal, and junction loops) and pseudoknots of arbitrary complexity. Notably, DSSR employs isolated pairs and the ends of stems, whether pseudoknotted or not, to define junction loops. This new, inclusive definition provides a novel perspective on the spatial organization of RNA. Tests on all nucleic acid structures in the Protein Data Bank confirm the efficiency and robustness of the software, and applications to representative RNA molecules illustrate its unique features. DSSR and related materials are freely available at http://x3dna.org/.

Structure-based prediction and characterization of photo-crosslinking in native protein-RNA complexes.

UV-crosslinking of protein and RNA in direct contacts has been widely used to study protein-RNA complexes while our understanding of the photo-crosslinking mechanisms remains poor. This knowledge gap is due to the challenge of precisely mapping the crosslink sites in protein and RNA simultaneously in their native sequence and structural contexts. Here we systematically analyze protein-RNA interactions and photo-crosslinking by bridging crosslinked nucleotides and amino acids mapped using different assays with protein-RNA complex structures. We developed a computational method PxR3D-map which reliably predicts crosslink sites using structural information characterizing protein-RNA interaction interfaces. Analysis of the informative features revealed that photo-crosslinking is facilitated by base stacking with not only aromatic residues, but also dipeptide bonds that involve glycine, and distinct mechanisms are utilized by different RNA-binding domains. Our work suggests protein-RNA photo-crosslinking is highly selective in the cellular environment, which can guide data interpretation and further technology development for UV-crosslinking-based assays.

Weak operator binding enhances simulated Lac repressor-mediated DNA looping.

The 50th anniversary of Biopolymers coincides closely with the like celebration of the discovery of the Escherichia coli (lac) lactose operon, a classic genetic system long used to illustrate the influence of biomolecular structure on function. The looping of DNA induced by the binding of the Lac repressor protein to sequentially distant operator sites on DNA continues to serve as a paradigm for understanding long-range genomic communication. Advances in analyses of DNA structures and in incorporation of proteins in computer simulations of DNA looping allow us to address long-standing questions about the role of protein-mediated DNA loop formation in transcriptional control. Here we report insights gained from studies of the sequence-dependent contributions of the natural lac operators to Lac repressor-mediated DNA looping. Novel superposition of the ensembles of protein-bound operator structures derived from NMR measurements reveals variations in DNA folding missed in conventional structural alignments. The changes in folding affect the predicted ease with which the repressor induces loop formation and the ways that DNA closes between the protein headpieces. The peeling of the auxiliary operators away from the repressor enhances the formation of loops with the 92-bp wildtype spacing and hints of a structural reason behind their weak binding.

The RNA backbone plays a crucial role in mediating the intrinsic stability of the GpU dinucleotide platform and the GpUpA/GpA miniduplex.

The side-by-side interactions of nucleobases contribute to the organization of RNA, forming the planar building blocks of helices and mediating chain folding. Dinucleotide platforms, formed by side-by-side pairing of adjacent bases, frequently anchor helices against loops. Surprisingly, GpU steps account for over half of the dinucleotide platforms observed in RNA-containing structures. Why GpU should stand out from other dinucleotides in this respect is not clear from the single well-characterized H-bond found between the guanine N2 and the uracil O4 groups. Here, we describe how an RNA-specific H-bond between O2'(G) and O2P(U) adds to the stability of the GpU platform. Moreover, we show how this pair of oxygen atoms forms an out-of-plane backbone 'edge' that is specifically recognized by a non-adjacent guanine in over 90% of the cases, leading to the formation of an asymmetric miniduplex consisting of 'complementary' GpUpA and GpA subunits. Together, these five nucleotides constitute the conserved core of the well-known loop-E motif. The backbone-mediated intrinsic stabilities of the GpU dinucleotide platform and the GpUpA/GpA miniduplex plausibly underlie observed evolutionary constraints on base identity. We propose that they may also provide a reason for the extreme conservation of GpU observed at most 5'-splice sites.

3DNALandscapes: a database for exploring the conformational features of DNA.

3DNALandscapes, located at: http://3DNAscapes.rutgers.edu, is a new database for exploring the conformational features of DNA. In contrast to most structural databases, which archive the Cartesian coordinates and/or derived parameters and images for individual structures, 3DNALandscapes enables searches of conformational information across multiple structures. The database contains a wide variety of structural parameters and molecular images, computed with the 3DNA software package and known to be useful for characterizing and understanding the sequence-dependent spatial arrangements of the DNA sugar-phosphate backbone, sugar-base side groups, base pairs, base-pair steps, groove structure, etc. The data comprise all DNA-containing structures--both free and bound to proteins, drugs and other ligands--currently available in the Protein Data Bank. The web interface allows the user to link, report, plot and analyze this information from numerous perspectives and thereby gain insight into DNA conformation, deformability and interactions in different sequence and structural contexts. The data accumulated from known, well-resolved DNA structures can serve as useful benchmarks for the analysis and simulation of new structures. The collective data can also help to understand how DNA deforms in response to proteins and other molecules and undergoes conformational rearrangements.

[Measuring soil water content by using near infrared spectral characteristics of soil].

Different textured soils (sandy loam, silty clay loam and clay) from Manas County, Xinjiang were researched with indoor spectral reflectance, continuum removal was used to process soil spectra curve, and correlation analysis was made about normalized spectral reflectance and water content for modeling. The results show that different textured soils have reflectance in a order that clay>silty clay loam>sandy loam; the critical points of field capacity in sandy loam, silty clay loam and clay were 20.01%, 24.10% and 30.43% respectively, and water content was inversely proportional to spectral reflectance below such critical points while proportional above the points. Within 1390-1623 nm band, the negative correlation coefficients of soil water content and normalized spectra reflectance show better negative correlation and reaching significant levels, R2 of the model established for soil water content prediction exceeded mostly 0.8, and the average relative error of the model was 10%. The model could accurately reflect the soil moisture content, Its advantages such as accuracy, non-destruction and rapidness provide a new approach to measuring soil water content.

Web 3DNA--a web server for the analysis, reconstruction, and visualization of three-dimensional nucleic-acid structures.

The w3DNA (web 3DNA) server is a user-friendly web-based interface to the 3DNA suite of programs for the analysis, reconstruction, and visualization of three-dimensional (3D) nucleic-acid-containing structures, including their complexes with proteins and other ligands. The server allows the user to determine a wide variety of conformational parameters in a given structure--such as the identities and rigid-body parameters of interacting nucleic-acid bases and base-pair steps, the nucleotides comprising helical fragments, etc. It is also possible to build 3D models of arbitrary nucleotide sequences and helical types, customized single-stranded and double-helical structures with user-defined base-pair parameters and sequences, and models of DNA 'decorated' at user-defined sites with proteins and other molecules. The visualization component offers unique, publication-quality representations of nucleic-acid structures, such as 'block' images of bases and base pairs and stacking diagrams of interacting nucleotides. The w3DNA web server, located at http://w3dna.rutgers.edu, is free and open to all users with no login requirement.

New restraints and validation approaches for nucleic acid structures in PDB-REDO.

The quality of macromolecular structure models crucially depends on refinement and validation targets, which optimally describe the expected chemistry. Commonly used software for these two procedures has been designed and developed in a protein-centric manner, resulting in relatively few established features for the refinement and validation of nucleic acid-containing structure models. Here, new nucleic acid-specific approaches implemented in PDB-REDO are described, including a new restraint model using noncovalent geometries (base-pair hydrogen bonding and base-pair stacking) as refinement targets. New validation routines are also presented, including a metric for Watson-Crick base-pair geometry normality (ZbpG). Applying the PDB-REDO pipeline with the new restraint model to the whole Protein Data Bank (PDB) demonstrates an overall positive effect on the quality of nucleic acid-containing structure models. Finally, we discuss examples of improvements in the geometry of specific nucleic acid structures in the PDB. The new PDB-REDO models and pipeline are available at https://pdb-redo.eu/.

New information content in RNA base pairing deduced from quantitative analysis of high-resolution structures.

Non-canonical base pairs play important roles in organizing the complex three-dimensional folding of RNA. Here, we outline methodology developed both to analyze the spatial patterns of interacting base pairs in known RNA structures and to reconstruct models from the collective experimental information. We focus attention on the structural context and deformability of the seven pairing patterns found in greatest abundance in the helical segments in a set of well-resolved crystal structures, including (i-ii) the canonical A.U and G.C Watson-Crick base pairs, (iii) the G.U wobble pair, (iv) the sheared G.A pair, (v) the A.U Hoogsteen pair, (vi) the U.U wobble pair, and (vii) the G.A Watson-Crick-like pair. The non-canonical pairs stand out from the canonical associations in terms of apparent deformability, spanning a broader range of conformational states as measured by the six rigid-body parameters used to describe the spatial arrangements of the interacting bases, the root-mean-square deviations of the base-pair atoms, and the fluctuations in hydrogen-bonding geometry. The deformabilties, the modes of base-pair deformation, and the preferred sites of occurrence depend on sequence. We also characterize the positioning and overlap of the base pairs with respect to the base pairs that stack immediately above and below them in double-helical fragments. We incorporate the observed positions of the bases, base pairs, and intervening phosphorus atoms in models to predict the effects of the non-canonical interactions on overall helical structure.

Toward a Mechanistic Understanding of DNA Methylation Readout by Transcription Factors.

Epigenetic DNA modification impacts gene expression, but the underlying molecular mechanisms are only partly understood. Adding a methyl group to a cytosine base locally modifies the structural features of DNA in multiple ways, which may change the interaction with DNA-binding transcription factors (TFs) and trigger a cascade of downstream molecular events. Cells can be probed using various functional genomics assays, but it is difficult to disentangle the confounded effects of DNA modification on TF binding, chromatin accessibility, intranuclear variation in local TF concentration, and rate of transcription. Here we discuss how high-throughput in vitro profiling of protein-DNA interactions has enabled comprehensive characterization and quantification of the methylation sensitivity of TFs. Despite the limited structural data for DNA containing methylated cytosine, automated analysis of structural information in the Protein Data Bank (PDB) shows how 5-methylcytosine (5mC) can be recognized in various ways by amino acid side chains. We discuss how a context-dependent effect of methylation on DNA groove geometry can affect DNA binding by homeodomain proteins and how principled modeling of ChIP-seq data can overcome the confounding that makes the interpretation of in vivo data challenging. The emerging picture is that epigenetic modifications affect TF binding in a highly context-specific manner, with a direction and effect size that depend critically on their position within the TF binding site and the amino acid sequence of the TF. With this improved mechanistic knowledge, we have come closer to understanding how cells use DNA modification to acquire, retain, and change their identity.

Reduced Gja5 expression in arterial endothelial cells impairs arteriogenesis during acute ischemic cardiovascular disease.

The aim of the present study was to investigate the functional role of gap junction protein α 5 (Gja5) in arterial endothelial cells in the arteriogenesis that occurs during acute ischemic cardiovascular disease. Gja5 knockout mice and the femoral artery occlusion (FAO) model were used in the current study. Perfusions of both hindlimbs were obtained separately prior to FAO, immediately following FAO and 1, 3, 7, 14 and 21 days after FAO using a Laser Doppler Flow Imager. Genetic evidence concerning the gastrocnemicus (GC) muscle was collected by reverse transcription-quantitative polymerase chain reaction. There were significant reductions in the hindlimb perfusion of Gja5-/- mice compared with Gja5+/+ mice 1, 3, 7, 14 and 21 days following FAO. In Gja5+/- and in Gja5+/+ mice, the expression of Gja5 in the GC muscle was increased 4-fold in the ischemic hindlimb 3 days following FAO. Levels of Gja5 expression then returned to baseline values 7 days after FAO. The results of the present study demonstrated that arterial Gja5 expression serves a functional role in acute ischemic cardiovascular disease.

Parallels between DNA and collagen - comparing elastic models of the double and triple helix.

Multi-stranded helices are widespread in nature. The interplay of polymeric properties with biological function is seldom discussed. This study probes analogies between structural and mechanical properties of collagen and DNA. We modeled collagen with Eulerian rotational and translational parameters of adjacent rungs in the triple-helix ladder and developed statistical potentials by extracting the dispersion of the parameters from a database of atomic-resolution structures. The resulting elastic model provides a common quantitative way to describe collagen deformations upon interacting with integrins or matrix metalloproteinase and DNA deformations upon protein binding. On a larger scale, deformations in Type I collagen vary with a periodicity consistent with the D-periodic banding of higher-order fibers assemblies. This indicates that morphologies of natural higher-order collagen packing might be rooted in the characteristic deformation patterns.

3DNA: a software package for the analysis, rebuilding and visualization of three-dimensional nucleic acid structures.

We present a comprehensive software package, 3DNA, for the analysis, reconstruction and visualization of three-dimensional nucleic acid structures. Starting from a coordinate file in Protein Data Bank (PDB) format, 3DNA can handle antiparallel and parallel double helices, single-stranded structures, triplexes, quadruplexes and other complex tertiary folding motifs found in both DNA and RNA structures. The analysis routines identify and categorize all base interactions and classify the double helical character of appropriate base pair steps. The program makes use of a recently recommended reference frame for the description of nucleic acid base pair geometry and a rigorous matrix-based scheme to calculate local conformational parameters and rebuild the structure from these parameters. The rebuilding routines produce rectangular block representations of nucleic acids as well as full atomic models with the sugar-phosphate backbone and publication quality 'standardized' base stacking diagrams. Utilities are provided to locate the base pairs and helical regions in a structure and to reorient structures for effective visualization. Regular helical models based on X-ray diffraction measurements of various repeating sequences can also be generated within the program.

DNA Scrunching in the Packaging of Viral Genomes.

The motors that drive double-stranded DNA (dsDNA) genomes into viral capsids are among the strongest of all biological motors for which forces have been measured, but it is not known how they generate force. We previously proposed that the DNA is not a passive substrate but that it plays an active role in force generation. This "scrunchworm hypothesis" holds that the motor proteins repeatedly dehydrate and rehydrate the DNA, which then undergoes cyclic shortening and lengthening motions. These are captured by a coupled protein-DNA grip-and-release cycle to rectify the motion and translocate the DNA into the capsid. In this study, we examined the interactions of dsDNA with the dodecameric connector protein of bacteriophage ϕ29, using molecular dynamics simulations on four different DNA sequences, starting from two different conformations (A-DNA and B-DNA). In all four simulations starting with the protein equilibrated with A-DNA in the channel, we observed transitions to a common, metastable, highly scrunched conformation, designated A*. This conformation is very similar to one recently reported by Kumar and Grubmüller in much longer MD simulations on B-DNA docked into the ϕ29 connector. These results are significant for four reasons. First, the scrunched conformations occur spontaneously, without requiring lever-like protein motions often believed to be necessary for DNA translocation. Second, the transition takes place within the connector, providing the location of the putative "dehydrator". Third, the protein has more contacts with one strand of the DNA than with the other; the former was identified in single-molecule laser tweezer experiments as the "load-bearing strand". Finally, the spontaneity of the DNA-protein interaction suggests that it may play a role in the initial docking of DNA in motors like that of T4 that can load and package any sequence.

Transitions of Double-Stranded DNA Between the A- and B-Forms.

The structure of double-stranded DNA (dsDNA) is sensitive to solvent conditions. In solution, B-DNA is the favored conformation under physiological conditions, while A-DNA is the form found under low water activity. The A-form is induced locally in some protein-DNA complexes, and repeated transitions between the B- and A-forms have been proposed to generate the forces used to drive dsDNA into viral capsids during genome packaging. Here, we report analyses on previous molecular dynamics (MD) simulations on B-DNA, along with new MD simulations on the transition from A-DNA to B-DNA in solution. We introduce the A-B Index (ABI), a new metric along the A-B continuum, to quantify our results. When A-DNA is placed in an equilibrated solution at physiological ionic strength, there is no energy barrier to the transition to the B-form, which begins within about 1 ns. The transition is essentially complete within 5 ns, although occasionally a stretch of a few base pairs will remain A-like for up to ∼10 ns. A comparison of four sequences with a range of predicted A-phobicities shows that more A-phobic sequences make the transition more rapidly than less A-phobic sequences. Simulations on dsDNA with a region of roughly one turn locked in the A-form allow us to characterize the A/B junction, which has an average bend angle of 20-30°. Fluctuations in this angle occur with characteristic times of about 10 ns.