DNA supercoiling plays a major role in many cellular functions. The global DNA conformation is however intimately linked to local DNA-DNA interactions influencing both the physical properties and the biological functions of the supercoiled molecule. Juxtaposition of DNA double helices in ubiquitous crossover arrangements participates in multiple functions such as recombination, gene regulation and DNA packaging. However, little is currently known about how the structure and stability of direct DNA-DNA interactions influence the topological state of DNA. Here, a crystallographic analysis shows that due to the intrinsic helical chirality of DNA, crossovers of opposite handedness exhibit markedly different geometries. While right-handed crossovers are self-fitted by sequence-specific groove-backbone interaction and bridging Mg(2+) sites, left-handed crossovers are juxtaposed by groove-groove interaction. Our previous calculations have shown that the different geometries result in differential stabilisation in solution, in the presence of divalent cations. The present study reveals that the various topological states of the cell are associated with different inter-segmental interactions. While the unstable left-handed crossovers are exclusively formed in negatively supercoiled DNA, stable right-handed crossovers constitute the local signature of an unusual topological state in the cell, such as the positively supercoiled or relaxed DNA. These findings not only provide a simple mechanism for locally sensing the DNA topology but also lead to the prediction that, due to their different tertiary intra-molecular interactions, supercoiled molecules of opposite signs must display markedly different physical properties. Sticky inter-segmental interactions in positively supercoiled or relaxed DNA are expected to greatly slow down the slithering dynamics of DNA. We therefore suggest that the intrinsic helical chirality of DNA may have oriented the early evolutionary choices for DNA topology.

The three-dimensional crystal structures of HindIII bound to its cognate DNA with and without divalent cations were solved at 2.17 and 2.00 A resolution, respectively. HindIII forms a dimer. The structures showed that HindIII belongs to the EcoRI-like (alpha-class) subfamily of type II restriction endonucleases. The cognate DNA-complex structures revealed the specific DNA-recognition mechanism of HindIII by which it recognizes the palindromic sequence A/AGCTT. In the Mg(2+) ion-soaked structure the DNA was cleaved and two ions were bound at each active site, corresponding to the two-metal-ion mechanism.

For B-DNA, the strong linear correlation observed by nuclear magnetic resonance (NMR) between the (31)P chemical shifts (deltaP) and three recurrent internucleotide distances demonstrates the tight coupling between phosphate motions and helicoidal parameters. It allows to translate deltaP into distance restraints directly exploitable in structural refinement. It even provides a new method for refining DNA oligomers with restraints exclusively inferred from deltaP. Combined with molecular dynamics in explicit solvent, these restraints lead to a structural and dynamical view of the DNA as detailed as that obtained with conventional and more extensive restraints. Tests with the Jun-Fos oligomer show that this deltaP-based strategy can provide a simple and straightforward method to capture DNA properties in solution, from routine NMR experiments on unlabeled samples.

Experimental X-ray crystal structures and a database of calculated structural parameters of DNA octamers were used in combination to analyse the mechanics of DNA bending in the nucleosome core complex. The 1kx5 X-ray crystal structure of the nucleosome core complex was used to determine the relationship between local structure at the base step level and the global superhelical conformation observed for nucleosome-bound DNA. The superhelix is characterised by a large curvature (597 degrees ) in one plane and very little curvature (10 degrees ) in the orthogonal plane. Analysis of the curvature at the level of 10-step segments shows that there is a uniform curvature of 30 degrees per helical turn throughout most of the structure, but that there are two sharper kinks of 50 degrees at +/-2 helical turns from the central dyad base pair. The curvature is due almost entirely to the base step parameter roll. There are large periodic variations in roll, which are in phase with the helical twist, and account for 500 degrees of the total curvature. Although variations in the other base step parameters perturb the local path of the DNA, they make minimal contributions to the total curvature. This implies that DNA bending in the nucleosome is achieved using the roll-slide-twist degree of freedom previously identified as the major degree of freedom in naked DNA oligomers. The energetics of bending into a nucleosome-bound conformation were therefore analysed using a database of structural parameters that we have previously developed for naked DNA oligomers. The minimum energy roll, the roll flexibility force constant and the maximum and minimum accessible roll values were obtained for each base step in the relevant octanucleotide context to account for the effects of conformational coupling that vary with sequence context. The distribution of base step roll values and corresponding strain energy required to bend DNA into the nucleosome-bound conformation defined by the 1kx5 structure were obtained by applying a constant bending moment. When a single bending moment was applied to the entire sequence, the local details of the calculated structure did not match the experiment. However, when local 10-step bending moments were applied separately, the calculated structure showed excellent agreement with experiment. This implies that the protein applies variable bending forces along the DNA to maintain the superhelical path required for nucleosome wrapping. In particular, the 50 degrees kinks are constraints imposed by the protein rather than a feature of the 1kx5 DNA sequence. The kinks coincide with a relatively flexible region of the sequence, and this is probably a prerequisite for high affinity nucleosome binding, but the bending strain energy is significantly higher at these points than for the rest of the sequence. In the most rigid regions of the sequence, a higher strain energy is also required to achieve the standard 30 degrees curvature per helical turn. We conclude that matching of the DNA sequence to the local roll periodicity required to achieve bending, together with the increased flexibility required at the kinks, determines the sequence-selectivity of DNA wrapping in the nucleosome.

B-DNA flexibility, crucial for DNA-protein recognition, is sequence dependent. Free DNA in solution would in principle be the best reference state to uncover the relation between base sequences and their intrinsic flexibility; however, this has long been hampered by a lack of suitable experimental data. We investigated this relationship by compiling and analyzing a large dataset of NMR (31)P chemical shifts in solution. These measurements reflect the BI <--> BII equilibrium in DNA, intimately correlated to helicoidal descriptors of the curvature, winding and groove dimensions. Comparing the ten complementary DNA dinucleotide steps indicates that some steps are much more flexible than others. This malleability is primarily controlled at the dinucleotide level, modulated by the tetranucleotide environment. Our analyses provide an experimental scale called TRX that quantifies the intrinsic flexibility of the ten dinucleotide steps in terms of Twist, Roll, and X-disp (base pair displacement). Applying the TRX scale to DNA sequences optimized for nucleosome formation reveals a 10 base-pair periodic alternation of stiff and flexible regions. Thus, DNA flexibility captured by the TRX scale is relevant to nucleosome formation, suggesting that this scale may be of general interest to better understand protein-DNA recognition.

The 5R thymine glycol (5R-Tg) DNA lesion exists as a mixture of cis-(5R,6S) and trans-(5R,6R) epimers; these modulate base excision repair. We examine the 7:3 cis-(5R,6S):trans-(5R,6R) mixture of epimers paired opposite adenine in the 5'-GTgG-3' sequence with regard to nucleotide excision repair. Human XPA recognizes the lesion comparably to the C8-dG acetylaminoflourene (AAF) adduct, whereas XPC/HR23B recognition of Tg is superior. 5R-Tg is processed by the Escherichia coli UvrA and UvrABC proteins less efficiently than the C8-dG AAF adduct. For the cis-(5R, 6S) epimer Tg and A are inserted into the helix, remaining in the Watson-Crick alignment. The Tg N3H imine and A N(6) amine protons undergo increased solvent exchange. Stacking between Tg and the 3'-neighbor G*C base pair is disrupted. The solvent accessible surface and T(2) relaxation of Tg increases. Molecular dynamics calculations predict that the axial conformation of the Tg CH(3) group is favored; propeller twisting of the Tg*A pair and hydrogen bonding between Tg OH6 and the N7 atom of the 3'-neighbor guanine alleviate steric clash with the 5'-neighbor base pair. Tg also destabilizes the 5'-neighbor G*C base pair. This may facilitate flipping both base pairs from the helix, enabling XPC/HR23B recognition prior to recruitment of XPA.

Understanding the preference of non-specific proteins for certain DNA structural features requires an accurate description of the properties of free DNA, especially regarding their possible predisposition to adopt a conformation that favors the formation of a complex. Exploiting previous exhaustive NMR studies done on free DNA oligomers, we investigated the molecular basis of DNase I sensitivity, in conditions where DNase I binding limits the probability of cleavage. We showed that cleavage intensity was correlated with the adjacent 3' phosphate linkage flexibility, monitored by (31)P chemical shifts. Examining NMR-refined DNA structures highlighted that sequence-dependent flexible phosphates were associated with large minor groove variations that may promote the affinity of DNase I, according to the relevant DNA-protein complexes. In sum, this work demonstrates that specificity in DNA-DNase I interactions is mediated by DNA flexibility that influences the induced-fit transitions required to form productive complexes.

SUMMARY: The DNA inter base pair step parameters (Tilt, Roll, Twist, Shift, Slide, Rise) are a standard internal coordinate representation of DNA. In the absence of bend and shear, it is relatively easy to mentally visualize how Twist and Rise generate the familiar double helix. More complex structures do not readily yield to such intuition. For this reason we developed a plug-in for VMD that accepts a set of mathematical expressions as input and generates a coarse-grained model of DNA as output. This feature of VDNA appears to provide a unique approach to DNA modeling. Predefined expressions include: linear, sheared, bent, and circular DNA and models of the nucleosome superhelix, chromatin, thermal motion, and nucleosome unwrapping. AVAILABILITY: VDNA is preinstalled in VMD, http://www.ks.uiuc.edu/Research/vmd. Updates are at http://dna.ccs.tulane.edu. CONTACT: bishop@tulane.edu.

In a functioning genetic system, the information-encoding molecule must form a regular self-complementary complex (for example, the base-paired double helix of DNA) and it must be able to encode information and pass it on to new generations. Here we study a benzo-widened DNA-like molecule (yDNA) as a candidate for an alternative genetic set, and we explicitly test these two structural and functional requirements. The solution structure of a 10 bp yDNA duplex is measured by using 2D-NMR methods for a simple sequence composed of T-yA/yA-T pairs. The data confirm an antiparallel, right-handed, hydrogen-bonded helix resembling B-DNA but with a wider diameter and enlarged base-pair size. In addition to this, the abilities of two different polymerase enzymes (Klenow fragment of DNA pol I (Kf) and the repair enzyme Dpo4) to synthesize and extend the yDNA pairs T-yA, A-yT, and G-yC are measured by steady-state kinetics studies. Not surprisingly, insertion of complementary bases opposite yDNA bases is inefficient due to the larger base-pair size. We find that correct pairing occurs in several cases by both enzymes, but that common and relatively efficient mispairing involving T-yT and T-yC pairs interferes with fully correct formation and extension of pairs by these polymerases. Interestingly, the data show that extension of the large pairs is considerably more efficient with the flexible repair enzyme (Dpo4) than with the more rigid Kf enzyme. The results shed light on the properties of yDNA as a candidate for an alternative genetic information-encoding molecule and as a tool for application in basic science and biomedicine.

We describe Curves+, a new nucleic acid conformational analysis program which is applicable to a wide range of nucleic acid structures, including those with up to four strands and with either canonical or modified bases and backbones. The program is algorithmically simpler and computationally much faster than the earlier Curves approach, although it still provides both helical and backbone parameters, including a curvilinear axis and parameters relating the position of the bases to this axis. It additionally provides a full analysis of groove widths and depths. Curves+ can also be used to analyse molecular dynamics trajectories. With the help of the accompanying program Canal, it is possible to produce a variety of graphical output including parameter variations along a given structure and time series or histograms of parameter variations during dynamics.

Results concerning the flexibility of mononucleotide repeat DNA are presented using a novel methodology, denoted "SIR", to describe continuous changes in the structure of the nucleic acid. This methodology, combined with internal energy calculations and analytical energy gradients allows us to determine optimal conformations of poly(dG).poly(dC) and poly (dA).poly(dT) in both the A and B forms, taking into account the influence of the solvent medium and explicit counterions. Subsequently, several different types of distorsion of these optimal structures are investigated. It is shown that excellent correlation with experimental results concerning coupled changes in structural variables is obtained and several new correlations are also detected.

A novel and powerful methodology is developed which allows the alteration of molecular structures subjected to constraints and its application to polynucleotides with mononucleotide repeat symmetry, including the treatment of the flexible sugar rings is described. In contrast to procedures proposed by other authors, the constraints are formulated as differential equations which are linear with respect to the differentials of the geometrical variables. These equations can be solved easily by stepwise numerical integration involving sucessive infinitesimal rotations (SIR). Moreover, these equations define a set of independent curvilinear coordinates which can be used directly as the parameters of the energy functional in an energy minimisation procedure. This methodology allows the scanning of the full configurational space of a complex macromolecule, with direct access to the helicoidal variables in the case of periodic systems. Through this approach many problems involving biomacromolecular conformation, which would otherwise be intractable, may be studied with considerable ease.

A theoretical modelling of the interaction of putrescine (H3+N-(CH2)4-(+NH3) with DNA is carried out, introducing two new features which make the simulation of this interaction considerably more realistic. Firstly, the DNA to which putrescine is bound is fully flexible and thus able to respond to the distorting influence of the ligand. Secondly, the effect of changing the ratio of DNA base pairs per bound ligand is explicitly modelled. In this way, we have been able to confirm the experimentally known preference of putrescine binding with AT base pairs in B-DNA, but we also show, through the new features introduced, that the nature of the binding site of the ligand and the resulting impact on DNA conformation is strongly modified by the ligand binding density.

A new procedure for the analysis of the structure and molecular dynamics of duplex DNA is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the values of the conformational and helicoidal parameters for each structure entering the analysis using the method "Curves" developed by Lavery and Sklenar, J. Biomol. Str. Dyn. 6, 63 (1988), followed by a novel computer graphic display of the results. The graphic display is organized systematically using conformation wheels, or "dials", for each IUPAC torsional parameter and "windows" on the range of values assumed by the linear and angular helicoidal parameters, and is presented in a form isomorphous with the structure per se. The complete time evolution of the conformational and helicoidal parameters of a DNA double helix can then be depicted in a set of six composite figures. Dynamical aspects of helix bending are also subsumed in this analysis. The procedure is illustrated with an analysis of the structures of canonical A and B forms of DNA and the 300 degrees K native dodecamer duplex d(CGCGAATTCGCG). The "dials and windows" are then used for a comprehensive analysis of 30 psec of molecular dynamics on the dodecamer in the vicinity of a canonical B-DNA energy minimum. This involves presentation of the time evolution of 206 conformational and 230 helicoidal parameters for the dodecamer. A number of interesting structural features can be recognized in the analysis, including crankshaft motions, BI - BII transitions, sugar repuckerings, and a description of spontaneous helix bending at what corresponds to the 1 degrees and 2 degrees "hinge points" indicated in the crystal structure. Our approach is expected to be directly useful for critical analysis of the effects of various assumptions about force field parameters, hydration and electrostatic effects and thus contribute to the development of reliable simulation protocols for nucleic acid systems. Extension of the method to present differential changes in conformational and helicoidal parameters is expected to be valuable for the analysis of structural and molecular dynamics studies of the reorganization and adaptation of DNA on complexation with various drugs and regulatory proteins.

The algorithm "Curves", that we have recently presented in this journal (J. Biolmol. Str. Dynam. 6, 63-91 (1988], is updated to take into account the conventions developed at the Cambridge meeting on DNA curvature (September 1988) and extended to the calculation of local parameters. In addition, the principles which govern the choices made in establishing the Curves algorithm are compared with the approaches adopted by other authors.

We present a general and mathematically rigorous algorithm which allows the helicoidal structure of a protein to be calculated starting from the atomic coordinates of its peptide backbone. This algorithm yields a unique curved axis which quantifies the folding of the backbone and a full set of helicoidal parameters describing the location of each peptide unit. The parameters obtained form a complete and independent set and can therefore be used for analyzing, comparing, or reconstructing protein backbone geometry. This algorithm has been implemented in a computer program named P-Curve. Several examples of its possible applications are discussed.

A new procedure for the graphic analysis of molecular dynamics (MD) simulations on proteins is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the conformational and helicoidal parameters for each structure entering the analysis via the method "Curves," developed for proteins by Sklenar, Etchebest, and Lavery (Proteins: Structure, Function Genet. 6:46-60, 1989) followed by a novel computer graphic display of the results. The graphic display is organized systematically using conformation wheels ("dials") for each torsional parameter and "windows" on the range values assumed by the linear and angular helicoidal parameters, and is present in a form isomorphous with the primary structure per se. The complete time evolution of dynamic structure can then be depicted in a set of four composite figures. Dynamic aspects of secondary and tertiary structure are also provided. The procedure is illustrated with an analysis of a 50 psec in vacuo simulation on the 58 residue protein, bovine pancreatic trypsin inhibitor (BPTI), in the vicinity of the local minimum on the energy surface corresponding to a high resolution crystal structure. The time evolution of 272 conformational and 788 helicoidal parameters for BPTI is analyzed. A number of interesting features can be discerned in the analysis, including the dynamic range of conformational and helicoidal motions, the dynamic extent of 2 degrees structure motifs, and the calculated fluctuations in the helix axis. This approach is expected to be useful for a critical analysis of the effects of various assumptions about force field parameters, truncation of potentials, solvation, and electrostatic effects, and can thus contribute to the development of more reliable simulation protocols for proteins. Extensions of the analysis to present differential changes in conformational and helicoidal parameters is expected to be valuable in MD studies of protein complexes with substrates, inhibitors, and effectors and in determining the nature of structural changes in protein-protein interactions.

We describe Curves+, a new nucleic acid conformational analysis program which is applicable to a wide range of nucleic acid structures, including those with up to four strands and with either canonical or modified bases and backbones. The program is algorithmically simpler and computationally much faster than the earlier Curves approach, although it still provides both helical and backbone parameters, including a curvilinear axis and parameters relating the position of the bases to this axis. It additionally provides a full analysis of groove widths and depths. Curves+ can also be used to analyse molecular dynamics trajectories. With the help of the accompanying program Canal, it is possible to produce a variety of graphical output including parameter variations along a given structure and time series or histograms of parameter variations during dynamics.

In mammals, the sex of the embryo is determined during development by its commitment either to the male or female genetic program regulating testicular or ovarian organogenesis. Major steps towards unraveling sex determination in mammals are achieved by the identification of key genes involved in human pathologies and the application of mouse genetics to analyze their function. While the expression of Sry and Sox9 is sufficient to induce the male developmental program, the molecular pathways that specify ovarian differentiation were unclear before the recent demonstration that mutations in the RSPO1 gene induce female-to-male sex reversal in XX patients. By generating the corresponding mouse model, we have shown that Rspo1 is so far the earliest known gene controlling the female genetic developmental program. Rspo1 activates the canonical beta-catenin signaling pathway required for female somatic cell differentiation and germ cell commitment into meiosis. The aim of this review is to describe the roles of R-spondins (Rspo)in developmental processes and disorders and the current knowledge obtained from murine models. A particular focus will be on Rspo1 and its crucial function in sex determination.

We describe Curves+, a new nucleic acid conformational analysis program which is applicable to a wide range of nucleic acid structures, including those with up to four strands and with either canonical or modified bases and backbones. The program is algorithmically simpler and computationally much faster than the earlier Curves approach, although it still provides both helical and backbone parameters, including a curvilinear axis and parameters relating the position of the bases to this axis. It additionally provides a full analysis of groove widths and depths. Curves+ can also be used to analyse molecular dynamics trajectories. With the help of the accompanying program Canal, it is possible to produce a variety of graphical output including parameter variations along a given structure and time series or histograms of parameter variations during dynamics.

Regression equations of the Padé type (PADE W, PADE 1 /W) and equations with Simms' constants involved (SIMMS 1/W) were applied for the determination of internal parameters (i.e. concentrations, equilibrium constants, activity coefficient of hydrogen ions) related to a titrand-titrant system with one- or two-component sample tested with the species having acid-base properties. The parameters were searched according to an iteration procedure (MINUITS computer program applied).

Six rigid-body parameters (Shift, Slide, Rise, Tilt, Roll, Twist) are commonly used to describe the relative displacement and orientation of successive base pairs in a nucleic acid structure. The present work adapts this approach to describe the relative displacement and orientation of any two planes in an arbitrary molecule-specifically, planes which contain important pharmacophore elements. Relevant code from the 3DNA software package (Nucleic Acids Res. 2003, 31, 5108-5121) was generalized to treat molecular fragments other than DNA bases as input for the calculation of the corresponding rigid-body (or "planes") parameters. These parameters were used to construct feature vectors for a fuzzy relational clustering study of over 700 conformations of a flexible analogue of the dopamine reuptake inhibitor, GBR 12909. Several cluster validity measures were used to determine the optimal number of clusters. Translational (Shift, Slide, Rise) rather than rotational (Tilt, Roll, Twist) features dominate clustering based on planes that are relatively far apart, whereas both types of features are important to clustering when the pair of planes are close by. This approach was able to classify the data set of molecular conformations into groups and to identify representative conformers for use as template conformers in future Comparative Molecular Field Analysis studies of GBR 12909 analogues. The advantage of using the planes parameters, rather than the combination of atomic coordinates and angles between molecular planes used in our previous fuzzy relational clustering of the same data set (J. Chem. Inf. Model. 2005, 45, 610-623), is that the present clustering results are independent of molecular superposition and the technique is able to identify clusters in the molecule considered as a whole. This approach is easily generalizable to any two planes in any molecule.

Sodie and alluolminous orthoamphiboles from metamorphosed volcanic rocks develop helicoidal fractures upon cooling, due to elastic strain gradients radical to the c-axis. This condition may arise, fromradial composition gradients involving subastitution of tetrahedral silicon by aluminum in the double chain structure or from intergranular boundary constraints with adjacent quartz.

We develop a theoretical approach to hairpin-loop formation of single-stranded (ss) DNA by treating the strand as a two-state system in which bases are either "stacked" or "unstacked." The looping kinetics of ssDNA is shown to be intrinsically different from that of a wormlike chain; it is mainly controlled by stacking-breakage probability, not by the mean curvature of loops, and highly sensitive to the composition of the loop as seen in recent experiments. Our estimate of a stacking energy for poly ( dA ),-3.9 kcal/mol, is consistent with known results.