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.

What controls DNA looping?

The looping of DNA provides a means of communication between sequentially distant genomic sites that operate in tandem to express, copy, and repair the information encoded in the DNA base sequence. The short loops implicated in the expression of bacterial genes suggest that molecular factors other than the naturally stiff double helix are involved in bringing the interacting sites into close spatial proximity. New computational techniques that take direct account of the three-dimensional structures and fluctuations of protein and DNA allow us to examine the likely means of enhancing such communication. Here, we describe the application of these approaches to the looping of a 92 base-pair DNA segment between the headpieces of the tetrameric Escherichia coli Lac repressor protein. The distortions of the double helix induced by a second protein--the nonspecific nucleoid protein HU--increase the computed likelihood of looping by several orders of magnitude over that of DNA alone. Large-scale deformations of the repressor, sequence-dependent features in the DNA loop, and deformability of the DNA operators also enhance looping, although to lesser degrees. The correspondence between the predicted looping propensities and the ease of looping derived from gene-expression and single-molecule measurements lends credence to the derived structural picture.

Internucleosomal interactions mediated by histone tails allow distant communication in chromatin.

Action across long distances on chromatin is a hallmark of eukaryotic transcriptional regulation. Although chromatin structure per se can support long-range interactions, the mechanisms of efficient communication between widely spaced DNA modules in chromatin remain a mystery. The molecular simulations described herein suggest that transient binary internucleosomal interactions can mediate distant communication in chromatin. Electrostatic interactions between the N-terminal tails of the core histones and DNA enhance the computed probability of juxtaposition of sites that lie far apart along the DNA sequence. Experimental analysis of the rates of communication in chromatin constructs confirms that long-distance communication occurs efficiently and independently of distance on tail-containing, but not on tailless, chromatin. Taken together, our data suggest that internucleosomal interactions involving the histone tails are essential for highly efficient, long-range communication between regulatory elements and their targets in eukaryotic genomes.

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.

Long-range effects of histone point mutations on DNA remodeling revealed from computational analyses of SIN-mutant nucleosome structures.

The packaging of DNA into nucleosomes impedes the binding and access of molecules involved in its processing. The SWI/SNF multi-protein assembly, found in yeast, is one of many regulatory factors that stimulate the remodeling of DNA required for its transcription. Amino-acid point mutations in histones H3 or H4 partially bypass the requirement of the SWI/SNF complex in this system. The mechanisms underlying the observed remodeling, however, are difficult to discern from the crystal structures of nucleosomes bearing these so-called SIN (SWI/SNF INdependent) mutations. Here, we report detailed analyses of the conformations and interactions of the histones and DNA in these assemblies. We find that the loss of direct protein-DNA contacts near point-mutation sites, reported previously, is coupled to unexpected additional long-range effects, i.e. loss of intermolecular contacts and accompanying DNA conformational changes at sequentially and spatially distant sites. The SIN mutations seemingly transmit information relevant to DNA binding across the nucleosome. The energetic cost of deforming the DNA to the states found in the SIN-mutant structures helps to distinguish the mutants that show phenotypes in yeast from those that do not. Models incorporating these deformed dimer steps suggest ways that nucleosomal DNA may be remodeled during its biological processing.

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.

A novel roll-and-slide mechanism of DNA folding in chromatin: implications for nucleosome positioning.

How eukaryotic genomes encode the folding of DNA into nucleosomes and how this intrinsic organization of chromatin guides biological function are questions of wide interest. The physical basis of nucleosome positioning lies in the sequence-dependent propensity of DNA to adopt the tightly bent configuration imposed by the binding of the histone proteins. Traditionally, only DNA bending and twisting deformations are considered, while the effects of the lateral displacements of adjacent base pairs are neglected. We demonstrate, however, that these displacements have a much more important structural role than ever imagined. Specifically, the lateral Slide deformations observed at sites of local anisotropic bending of DNA define its superhelical trajectory in chromatin. Furthermore, the computed cost of deforming DNA on the nucleosome is sequence-specific: in optimally positioned sequences the most easily deformed base-pair steps (CA:TG and TA) occur at sites of large positive Slide and negative Roll (where the DNA bends into the minor groove). These conclusions rest upon a treatment of DNA that goes beyond the conventional ribbon model, incorporating all essential degrees of freedom of "real" duplexes in the estimation of DNA deformation energies. Indeed, only after lateral Slide displacements are considered are we able to account for the sequence-specific folding of DNA found in nucleosome structures. The close correspondence between the predicted and observed nucleosome locations demonstrates the potential advantage of our "structural" approach in the computer mapping of nucleosome positioning.

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.

Insights into Gene Expression and Packaging from Computer Simulations.

Within the nucleus of each cell lies DNA - an unfathomably long, twisted, and intricately coiled molecule - segments of which make up the genes that provide the instructions that a cell needs to operate. As we near the 60(th) anniversary of the discovery of the DNA double helix, crucial questions remain about how the physical arrangement of the DNA in cells affects how genes work. For example, how a cell stores the genetic information inside the nucleus is complicated by the necessity of maintaining accessibility to DNA for genetic processing. In order to gain insight into the roles played by various proteins in reading and compacting the genome, we have developed new methodologies to simulate the dynamic, three-dimensional structures of long, fluctuating, protein-decorated strands of DNA. Our a priori approach to the problem allows us to determine the effects of individual proteins and their chemical modifications on overall DNA structure and function. Here we present our recent treatment of the communication between regulatory proteins attached to precisely constructed stretches of chromatin. Our simulations account for the enhancement in communication detected experimentally on chromatin compared to protein-free DNA of the same chain length as well as the critical roles played by the cationic 'tails' of the histone proteins in this signaling. The states of chromatin captured in the simulations offer new insights into the ways that the DNA, histones, and regulatory proteins contribute to long-range communication along the genome.

Properties of the Nucleic-acid Bases in Free and Watson-Crick Hydrogen-bonded States: Computational Insights into the Sequence-dependent Features of Double-helical DNA.

The nucleic-acid bases carry structural and energetic signatures that contribute to the unique features of genetic sequences. Here we review the connection between the chemical structure of the constituent nucleotides and the polymeric properties of DNA. The sequence-dependent accumulation of charge on the major- and minor-groove edges of the Watson-Crick base pairs, obtained from ab initio calculations, presents unique motifs for direct sequence recognition. The optimization of base interactions generates a propellering of base-pair planes of the same handedness as that found in high-resolution double-helical structures. The optimized base pairs also deform along conformational pathways, i.e., normal modes, of the same type induced by the binding of proteins. Empirical energy computations that incorporate the properties of the base pairs account satisfactorily for general features of the next level of double-helical structure, but miss key sequence-dependent differences in dimeric structure and deformability. The latter discrepancies appear to reflect factors other than intrinsic base-pair structure.