Base sequence effects in double-helical DNA. II. Configurational statistics of rodlike chains

Biophysical Reviews' "Meet the Editors Series"-a profile of Wilma K. Olson

As one of the five Executive Editors of Biophysical Reviews I have been asked to provide this short biographical sketch for the readers of the journal. I have been a member of the Editorial Board since the inception of the journal in 2008 and an Executive Editor since 2014. I hold a B.S. degree in Chemistry from the University of Delaware and a Ph.D. in Physical Chemistry from Stanford University. Except for a year as a Damon Runyon Postdoctoral Fellow at Columbia University, I have spent my entire professional career at Rutgers, the State University of New Jersey, where I am currently the Mary I. Bunting Professor of Chemistry and Chemical Biology. I served for many years as Founding Director of the Rutgers University Center for Molecular Biophysics and Biophysical Chemistry, and have trained undergraduate, graduate, and postdoctoral students from a variety of academic disciplines.

RNA 3D structure prediction by using a coarse-grained model and experimental data

RNAs form complex secondary and three-dimensional structures, and their biological functions highly rely on their structures and dynamics. Here we developed a general coarse-grained framework for RNA 3D structure prediction. A new, hybrid coarse-grained model that explicitly describes the electrostatics and hydrogen-bond interactions has been constructed based on experimental structural statistics. With the simulated annealing simulation protocol, several RNAs of less than 30-nt were folded to within 4.0 Å of the native structures. In addition, with limited restraints on Watson-Crick basepairing based on the data from NMR spectroscopy and small-angle X-ray scattering (SAXS) information, the current model was able to characterize the complex tertiary structures of large size RNAs, such as 5S ribosome and U2/U6 snRNA. We also demonstrated that the pseudoknot structure was better captured when the coordinating Mg(2+) cations and limited basepairing restraints were included. The accuracy of our model has been compared favorably with other RNA structure prediction methods presented in the previous study of RNA-Puzzles. Therefore the coarse-grained model presented here offers a unique approach for accurate prediction and modeling of RNA structures.

Coarse-grained model for simulation of RNA three-dimensional structures

The accurate prediction of an RNA's three-dimensional structure from its "primary structure" will have a tremendous influence on the experimental design and its interpretation and ultimately our understanding of the many functions of RNA. This paper presents a general coarse-grained (CG) potential for modeling RNA 3-D structures. Each nucleotide is represented by five pseudo atoms, two for the backbone (one for the phosphate and another for the sugar) and three for the base to represent base-stacking interactions. The CG potential has been parametrized from statistical analysis of 688 RNA experimental structures. Molecular dynamic simulations of 15 RNA molecules with the length of 12-27 nucleotides have been performed using the CG potential, with performance comparable to that from all-atom simulations. For ~75% of systems tested, simulated annealing led to native-like structures at least once out of multiple repeated runs. Furthermore, with weak distance restraints based on the knowledge of three to five canonical Watson-Crick pairs, all 15 RNAs tested are successfully folded to within 6.5 Å of native structures using the CG potential and simulated annealing. The results reveal that with a limited secondary structure model the current CG potential can reliably predict the 3-D structures for small RNA molecules. We also explored an all-atom force field to construct atomic structures from the CG simulations.

Group theory and biomolecular conformation: I. Mathematical and computational models

Biological macromolecules, and the complexes that they form, can be described in a variety of ways ranging from quantum mechanical and atomic chemical models, to coarser grained models of secondary structure and domains, to continuum models. At each of these levels, group theory can be used to describe both geometric symmetries and conformational motion. In this survey, a detailed account is provided of how group theory has been applied across computational structural biology to analyze the conformational shape and motion of macromolecules and complexes.

The Stochastic Elastica and Excluded-Volume Perturbations of DNA Conformational Ensembles

A coordinate-free Lie-group formulation for generating ensembles of DNA conformations in solution is presented. In this formulation, stochastic differential equations define sample paths on the Euclidean motion group. The ensemble of these paths exhibits the same behavior as solutions of the Fokker-Planck equation for the stochastically forced elastica. Longer chains for which the effects of excluded volume become important are handled by piecing together shorter chains and modeling their interactions. It is assumed that the final chain lengths of interest are long enough for excluded volume effects to become important, but not so long that the semi-flexible nature of the chain is lost. The effect of excluded volume is then taken into account by grouping short self-avoiding conformations into 'bundles' with common end constraints and computing average interaction effects between bundles. The accuracy of this approximation is shown to be good when using a numerically generated ensemble of self-avoiding sample paths as the baseline for comparison.

Effects of the nucleoid protein HU on the structure, flexibility, and ring-closure properties of DNA deduced from Monte Carlo simulations

The histone-like HU (heat unstable) protein plays a key role in the organization and regulation of the Escherichia coli genome. The nonspecific nature of HU binding to DNA complicates analysis of the mechanism by which the protein contributes to the looping of DNA. Conventional models of the looping of HU-bound duplexes attribute the changes in biophysical properties of DNA brought about by the random binding of protein to changes in the effective parameters of an ideal helical wormlike chain. Here, we introduce a novel Monte Carlo approach to study the effects of nonspecific HU binding on the configurational properties of DNA directly. We randomly decorated segments of an ideal double-helical DNA with HU molecules that induce the bends and other structural distortions of the double helix find in currently available X-ray structures. We find that the presence of HU at levels approximating those found in the cell reduces the persistence length by roughly threefold compared with that of naked DNA. The binding of protein has particularly striking effects on the cyclization properties of short duplexes, altering the dependence of ring closure on chain length in a way that cannot be mimicked by a simple wormlike model and accumulating at higher-than-expected levels on successfully closed chains. Moreover, the uptake of protein on small minicircles depends on chain length, taking advantage of the HU-induced deformations of DNA structure to facilitate ligation. Circular duplexes with bound HU show much greater propensity than protein-free DNA to exist as negatively supercoiled topoisomers, suggesting a potential role of HU in organizing the bacterial nucleoid. The local bending and undertwisting of DNA by HU, in combination with the number of bound proteins, provide a structural rationale for the condensation of DNA and the observed expression levels of reporter genes in vivo.

A coarse grain model for DNA

Understanding the behavior of DNA at the molecular level is of considerable fundamental and engineering importance. While adequate representations of DNA exist at the atomic and continuum level, there is a relative lack of models capable of describing the behavior of DNA at mesoscopic length scales. We present a mesoscale model of DNA that reduces the complexity of a nucleotide to three interactions sites, one each for the phosphate, sugar, and base, thereby rendering the investigation of DNA up to a few microns in length computationally tractable. The charges on these sites are considered explicitly. The model is parametrized using thermal denaturation experimental data at a fixed salt concentration. The validity of the model is established by its ability to predict several aspects of DNA behavior, including salt-dependent melting, bubble formation and rehybridization, and the mechanical properties of the molecule as a function of salt concentration.

Conformational Statistics of Semi-Flexible Macromolecular Chains with Internal Joints

Fluctuations in the bending angles at internal irregularities of DNA and RNA (such as symmetric loops, bulges, and nicks/gaps) have been observed from various experiments. However, little effort has been made to computationally predict and explain the statistical behavior of semi-flexible chains with internal defects. In this paper, we describe the general structure of these macromolecular chains as inextensible elastic chains with one or more internal joints which have limited ranges of rotation, and propose a method to compute the probability density functions of the end-to-end pose of these macromolecular chains. Our method takes advantage of the operational properties of the non-commutative Fourier transform for the group of rigid-body motions in three-dimensional space, SE(3). Two representative types of joints, the hinge for planar rotation and the ball joint for spatial rotation, are discussed in detail. The proposed method applies to various stiffness models of semi-flexible chain-like macromolecules. Examples are calculated using the Kratky-Porod model with specified stiffness, angular fluctuation, and joint locations. Entropic effects associated with internal angular fluctuations of semi-flexible macromolecular chains with internal joints can be computed using this formulation. Our method also provides a potential tool to detect the existence of internal irregularities.

Conformational Analysis of Stiff Chiral Polymers with End-Constraints

We present a Lie-group-theoretic method for the kinematic and dynamic analysis of chiral semi-flexible polymers with end constraints. The first is to determine the minimum energy conformations of semi-flexible polymers with end constraints, and the second is to perform normal mode analysis based on the determined minimum energy conformations. In this paper, we use concepts from the theory of Lie groups and principles of variational calculus to model such polymers as inextensible or extensible chiral elastic rods with coupling between twisting and bending stiffnesses, and/or between twisting and extension stiffnesses. This method is general enough to include any stiffness and chirality parameters in the context of elastic filament models with the quadratic elastic potential energy function. As an application of this formulation, the analysis of DNA conformations is discussed. We demonstrate our method with examples of DNA conformations in which topological properties such as writhe, twist, and linking number are calculated from the results of the proposed method. Given these minimum energy conformations, we describe how to perform the normal mode analysis. The results presented here build both on recent experimental work in which DNA mechanical properties have been measured, and theoretical work in which the mechanics of non-chiral elastic rods has been studied.

Sequence-dependent motions of DNA: a normal mode analysis at the base-pair level

Computer simulation of the dynamic structure of DNA can be carried out at various levels of resolution. Detailed high resolution information about the motions of DNA is typically collected for the atoms in a few turns of double helix. At low resolution, by contrast, the sequence-dependence features of DNA are usually neglected and molecules with thousands of base pairs are treated as ideal elastic rods. The present normal mode analysis of DNA in terms of six base-pair "step" parameters per chain residue addresses the dynamic structure of the double helix at intermediate resolution, i.e., the mesoscopic level of a few hundred base pairs. Sequence-dependent effects are incorporated into the calculations by taking advantage of "knowledge-based" harmonic energy functions deduced from the mean values and dispersion of the base-pair "step" parameters in high-resolution DNA crystal structures. Spatial arrangements sampled along the dominant low frequency modes have end-to-end distances comparable to those of exact polymer models which incorporate all possible chain configurations. The normal mode analysis accounts for the overall bending, i.e., persistence length, of the double helix and shows how known discrepancies in the measured twisting constants of long DNA molecules could originate in the deformability of neighboring base-pair steps. The calculations also reveal how the natural coupling of local conformational variables affects the global motions of DNA. Successful correspondence of the computed stretching modulus with experimental data requires that the DNA base pairs be inclined with respect to the direction of stretching, with chain extension effected by low energy transverse motions that preserve the strong van der Waals' attractions of neighboring base-pair planes. The calculations further show how one can "engineer" the macroscopic properties of DNA in terms of dimer deformability so that polymers which are intrinsically straight in the equilibrium state exhibit the mesoscopic bending anisotropy essential to DNA curvature and loop formation.

A refined prediction method for gel retardation of DNA oligonucleotides from dinucleotide step parameters: reconciliation of DNA bending models with crystal structure data

The development and assessment of a prediction method for gel retardation and sequence dependent curvature of DNA based on dinulcleotide step parameters are described. The method is formulated using the Babcock-Olson equations for base pair step geometry (1) and employs Monte Carlo simulated annealing for parameter optimization against experimental data. The refined base pair step parameters define a stuctural construct which, when the width of observed parameter distributions is taken into account, is consistent with the results of DNA oligonucleotide crystal structures. The predictive power of the method is demonstrated and tested via comparisons with DNA bending data on sets of sequences not included in the training set, including A-tracts with and without periodic helix phasing, phased A4T4 and T4A4 motifs, a sequence with a phased GGGCCC motif, some "unconventional" helix phasing sequences, and three short fragments of kinetoplast DNA from Crithidia fasiculata that exhibit significantly different behavior on non-denaturing polyacrylamide gels. The nature of the structural construct produced by the methodology is discussed with respect to static and dynamic models of structure and representations of bending and bendability. An independent theoretical account of sequence dependent chemical footprinting results is provided. Detailed analysis of sequences with A-tract induced axis bending forms the basis for a critical discussion of the applicability of wedge models,junction models and non A-tract, general sequence models for understanding the origin of DNA curvature at the molecular level.

Nucleic acids: theory and computer simulation, Y2K

Molecular dynamics simulations on DNA and RNA that include solvent are now being performed under realistic environmental conditions of water activity and salt. Improvements to force-fields and treatments of long-range interactions have significantly increased the reliability of simulations. New studies of sequence effects, axis bending, solvation and conformational transitions have appeared.

Factors Contributing to Aromatic Stacking in Water: Evaluation in the Context of DNA

We report the use of thermodynamic measurements in a self-complementary DNA duplex (5'-dXCGCGCG)(2), where X is an unpaired natural or nonnatural deoxynucleoside, to study the forces that stabilize aqueous aromatic stacking in the context of DNA. Thermal denaturation experiments show that the core duplex (lacking X) is formed with a free energy (37 °C) of -8.1 kcal·mol(-1) in a pH 7.0 buffer containing 1 M Na(+). We studied the effects of adding single dangling nucleosides (X) where the aromatic "base" is adenine, guanine, thymine, cytosine, pyrrole, benzene, 4-methylindole, 5-nitroindole, trimethylbenzene, difluorotoluene, naphthalene, phenanthrene, and pyrene. Adding these dangling residues is found to stabilize the duplex by an additional -0.8 to -3.4 kcal·mol(-1). At 5 μM DNA concentration, T(m) values range from 41.7 °C (core sequence) to 64.1 °C (with dangling pyrene residues). For the four natural bases, the order of stacking ability is A > G ≥ T = C. The nonpolar analogues stack more strongly in general than the more polar natural bases. The stacking geometry was confirmed in two cases (X = adenine and pyrene) by 2-D NOESY experiments. Also studied is the effect of ethanol cosolvent on the stacking of natural bases and pyrene. Stacking abilities were compared to calculated values for hydrophobicity, dipole moment, polarizability, and surface area. In general, hydrophobic effects are found to be larger than other effects stabilizing stacking (electrostatic effects, dispersion forces); however, the natural DNA bases are found to be less dependent on hydrophobic effects than are the more nonpolar compounds. The results also point out strategies for the design nucleoside analogues that stack considerably more strongly than the natural bases; such compounds may be useful in stabilizing designed DNA structures and complexes.

Molecular dynamics studies of axis bending in d(G5-(GA4T4C)2-C5) and d(G5-(GT4A4C)2-C5): effects of sequence polarity on DNA curvature

Gel retardation studies and other experiments indicate that DNA sequences containing the d(GA4T4C)n motif are curved, whereas those of identical composition but with a reverse sequence polarity, the d(GT4A4C)n motif, are straight. Hydroxyl radical cleavage experiments show that d(GA4T4C)n shows a unique signature, whereas d(GT4A4C)n behaves normally. To explain these results at a molecular level, molecular dynamics (MD) simulations were performed on the DNA duplexes d(G5-(GA4T4C)2-C5) and d(G5-(GT4A4C)2-C5) to 3.0 and 2.5 ns, respectively. The MD simulations are based on the Cornell force field implemented in the AMBER 4.1 modeling package and performed in a neutral solution of anionic DNA with K+, Cl- and Mg2+ at concentrations roughly comparable to a ligase buffer. Long range interactions were treated by the particle mesh Ewald method. Analysis of the results shows that the calculated dynamical structure of d(G5-(GA4T4C)2-C5) exhibits strong gross curvature, consistent with the observed behavior. The most significant locus of curvature in the MD structure is found at the central C15-G16 step, with an average roll angle of 12.8(+/-6.40)deg. The d(G5-(GT4A4C)2-C5) MD structure exhibited significantly less gross curvature. Analysis of results indicates that the reduction in gross curvature in the d(G5-(GT4A4C)2-C5) trajectory originates from the effect of the T10-A11 and T20-A21 steps, which showed average roll angles of 12.5(+/-5)deg. These three steps, T10-A11, C15-G16 and T20-A21, are half-helix turns away from one another, and their contributions to concerted bending cancel out. The A-tracts in the MD structure are essentially straight. The dynamical structure of d(G5-(GA4T4C)2-C5) exhibited minor groove deformation comprised of expansion at the 5' end of A-tracts and progressive narrowing towards the 3' end, consistent with and elaborating the interpretation of hydroxyl radical chemical probing results.

Measurement of the DNA bend angle induced by the catabolite activator protein using Monte Carlo simulation of cyclization kinetics

A Monte Carlo simulation method for studying DNA cyclization (or ring-closure) has been extended to the case of protein-induced bending, and its application to experimental data has been demonstrated. Estimates for the geometric parameters describing the DNA bend induced by the catabolite activator protein (CAP or CRP) were obtained which correctly predict experimental DNA cyclization probabilities (J factors), determined for a set of 11 150 to 166 bp DNA restriction fragments bearing A tracts phased against CAP binding sites. We find that simulation of out-of-phase molecules is difficult and time consuming, requiring the geometric parameters to be optimized individually rather than globally. A wedge angle model for DNA bending was found to make reasonable predictions for the free DNA. The bend angle in the CAP-DNA complex is estimated to be 85 to 90 degrees, in agreement with estimates from gel electrophoresis and X-ray co-crystal structures. Since the DNA is found to have a pre-existing bend of 15 degrees, the change in bend angle induced by CAP is 70 to 75 degrees, in a agreement with an estimate from topological measurements. We find evidence for slight (approximately 10 degrees) unwinding by CAP. The persistence length and helical repeat of the unbound portion of the DNA are in accord with literature-cited values, but the best-fit DNA torsional modulus C is found to be 1.7 (+/- 0.2) x 10(-19) erg. cm, versus literature estimates and best-fit values for the free DNA of 2.0 x 10(-19) to 3.4 x 10(-19) erg.com. Simulations using this low value of C predict that cyclization of molecules with out-of-phase bends proceeds via undertwisting or overtwisting of the DNA between the bends, so as to align the bends, rather than through conformations with substantial writhe. We present experiments on the topoisomers formed by cyclization with CAP which support this conclusion, and thereby rationalize the surprising result that cyclization can actually be enhanced by out-of-phase bends if the twist required to align the bends improves the torsional alignment of the ends. The relationship between the present work and previous studies on DNA bending by CAP is discussed, and recommendations are given for the efficient application of the cyclization/simulation approach to DNA bending.

Configurational transitions in Fourier series-represented DNA supercoils

A new Fourier series representation of supercoiled DNA is employed in Langevin dynamics simulations to study large-scale configurational motions of intermediate-length chains. The polymer is modeled as an ideal elastic rod subject to long-range van der Waals' interactions. The van der Waals' term prevents the self-contact of distant chain segments and also mimics attractive forces thought to stabilize the association of closely spaced charged rods. The finite Fourier series-derived polymer formulation is an alternative to the piecewise B-spline curves used in past work to describe the motion of smoothly deformed supercoiled DNA in terms of a limited number of independent variables. This study focuses on two large-scale configurational events: the interconversion between circular and figure-8 forms at a relatively low level of supercoiling, and the transformation between branched and interwound structures at a higher superhelical density.

Molecular dynamics simulations of small DNA plasmids: effects of sequence and supercoiling on intramolecular motions

Small (600 base pair) DNA plasmids were modeled with a simplified representation (3DNA) and the intramolecular motions were studied using molecular mechanics and molecular dynamics techniques. The model is detailed enough to incorporate sequence effects. At the same time, it is simple enough to allow long molecular dynamics simulations. The simulations revealed that large-scale slithering occurs in a homogeneous sequence. In a heterogeneous sequence, containing numerous small intrinsic curves, the centers of the curves are preferentially positioned at the tips of loops. With more curves than loop tips (two in unbranched supercoiled DNA), the heterogeneous sequence plasmid slithers short distances to reposition other curves into the loop tips. However, the DNA is immobilized most of the time, with the loop tips positioned over a few favored curve centers. Branching or looping also appears in the heterogeneous sequence as a new method of repositioning the loop tips. Instead of a smooth progression of increasing writhing with increasing linking difference, theoretical studies have predicted that there is a threshold between unwrithed and writhed DNA at a linking difference between one and two. This has previously been observed in simulations of static structures and is demonstrated here for dynamic homogeneous closed DNA. Such an abrupt transition is not found in the heterogeneous sequence in both the static and dynamic cases.

Analysis of local helix bending in crystal structures of DNA oligonucleotides and DNA-protein complexes

Sequence-dependent bending of the helical axes in 112 oligonucleotide duplex crystal structures resident in the Nucleic Acid Database have been analyzed and compared with the use of bending dials, a computer graphics tool. Our analysis includes structures of both A and B forms of DNA and considers both uncomplexed forms of the double helix as well as those bound to drugs and proteins. The patterns in bending preferences in the crystal structures are analyzed by base pair steps, and emerging trends are noted. Analysis of the 66 B-form structures in the Nucleic Acid Database indicates that uniform trends within all pyrimidine-purine and purine-pyrimidine steps are not necessarily observed but are found particularly at CG and GC steps of dodecamers. The results support the idea that AA steps are relatively straight and that larger roll bends occur at or near the junctions of these A-tracts with their flanking sequences. The data on 16 available crystal structures of protein-DNA complexes indicate that the majority of the DNA bends induced via protein binding are sharp localized kinks. The analysis of the 30 available A-form DNA structures indicates that these structures are also bent and show a definitive preference for bending into the deep major groove over the shallow minor groove.

Flexing and folding double helical DNA

DNA base sequence, once thought to be interesting only as a carrier of the genetic blueprint, is now recognized as playing a structural role in modulating the biological activity of genes. Primary sequences of nucleic acid bases describe real three-dimensional structures with properties reflecting those structures. Moreover, the structures are base sequence dependent with individual residues adopting characteristic spatial forms. As a consequence, the double helix can fold into tertiary arrangements, although the deformation is much more gradual and spread over a larger molecular scale than in proteins. As part of an effort to understand how local structural irregularities are translated at the macromolecular level in DNA and recognized by proteins, a series of calculations probing the structure and properties of the double helix have been performed. By combining several computational techniques, complementary information as well as a series of built-in checks and balances for assessing the significance of the findings are obtained. The known sequence dependent bending, twisting, and translation of simple dimeric fragments have been incorporated into computer models of long open DNAs of varying length and chemical composition as well as in closed double helical circles and loops. The extent to which the double helix can be forced to bend and twist is monitored with newly parameterized base sequence dependent elastic energy potentials based on the observed configurations of adjacent base pairs in the B-DNA crystallographic literature.

Static contributions to the persistence length of DNA and dynamic contributions to DNA curvature

Long molecules of DNA have the statistical properties of a worm-like coil. Deviations from linearity occur both because of small dynamic bends induced by thermal motion and from a random distribution of static bends. The latter originate in the different conformations of each of the possible base pair sequences. In this paper a statistical theory of the persistence length of DNA is developed which includes both static and dynamic effects for each base pair sequence, as well as the sequence-dependent correlations of bending angles. The result applies to a generic DNA, i.e., the average over an ensemble of all possible sequences. The theory is also applied to the generation of the average properties of curved DNAs by an analytic method that includes dynamic averaging as well as correlated bends. These results provide information which supplements that obtained by others using Monte Carlo methods. The additivity relation 1/P = 1/P(S) + 1/P(d) proposed by Trifonov et al., where P is the persistence length and P(S) and P(d) are the persistence lengths arising from purely static and dynamic effects, respectively, has been verified to be accurate to better than 0.5%. This is true for both a simplified model and one that includes a complete set of static bends at all base pair sequences.

Hydrophobic, Non-Hydrogen-Bonding Bases and Base Pairs in DNA

We report the properties of hydrophobic isosteres of pyrimidines and purines in synthetic DNA duplexes. Phenyl nucleosides 1 and 2 are nonpolar isosteres of the natural thymidine nucleoside, and indole nucleoside 3 is an analog of the complementary purine 2-aminodeoxyadenosine. The nucleosides were incorporated into synthetic oligodeoxynucleotides and were paired against each other and against the natural bases. Thermal denaturation experiments were used to measure the stabilities of the duplexes at neutral pH. It is found that the hydrophobic base analogs are nonselective in pairing with the four natural bases but selective for pairing with each other rather than with the natural bases. For example, compound 2 selectively pairs with itself rather than with A, T, G, or C; the magnitude of this selectivity is found to be 6.5-9.3 °C in Tm or 1.5-1.8 kcal/mol in free energy (25 °C). All possible hydrophobic pairing combinations of 1, 2, and 3 were examined. Results show that the pairing affinity depends on the nature of the pairs and on position in the duplex. The highest affinity pairs are found to be the 1-1 and 2-2 self-pairs and the 1-2 heteropair. The best stabilization occurs when the pairs are placed at the ends of duplexes rather than internally; the internal pairs may be destabilized by imperfect steric mimicry which leads to non-ideal duplex structure. In some cases the hydrophobic pairs are significantly stabilizing to the DNA duplex; for example, when situated at the end of a duplex, the 1-1 pair is more stabilizing than a T-A pair. When situated internally, the affinity of the 1-1 pair is the same as, or slightly better than, the analogous T-T mismatch pair, which is known to have two hydrogen bonds. The studies raise the possibility that hydrogen bonds may not always be required for the formation of stable duplex DNA-like structure. In addition, the results point out the importance of solvation and desolvation in natural base pairing, and lend new support to the idea that hydrogen bonds in DNA may be more important for specificity of pairing than for affinity. Finally, the study raises the possibility of using these or related base pairs to expand the genetic code beyond the natural A-T and G-C pairs.

Configurational statistics of the DNA duplex: extended generator matrices to treat the rotations and translations of adjacent residues

The base-to-base virtual bond treatment of nucleic acids used in statistical mechanical calculations of polynucleotide chain properties has been refined by incorporating the six parameters that relate the positions and orientations of sequential rigid bodies. The scheme allows for the sequence-dependent bending, twisting, and displacement of base pairs as well as for asymmetry in the angular and translational fluctuations of individual residues. Expressions are developed for the generator matrices required for the computation, as a function of chain length, of various parameters measuring the overall mean extension and shape of the DNA. Quantities of interest include the end-to-end vector r, the square of the end-to-end distance r2, the square radius of gyration s2, the center-of-gravity vector g, the second moments of inertia Sx2, and the higher moments of r and g. The matrix expressions introduced in the 1960s by Flory and co-workers for the determination of configuration-dependent polymer chain averages are decomposed into their translational and orientational contributions so that the methods can be extended to the rigid body analysis of chemical moieties. The new expressions permit, for the first time, examination of the effects of sequence-dependent translations, such as the lateral sliding of residues in A- and B-helices and the vertical opening of base pairs in drug-DNA complexes, on the average extension and shape of the long flexible double helix. The approach is in the following paper using conformational energy estimates of the base sequence-dependent flexibility of successive B-DNA base pairs.

Spatial translational motions of base pairs in DNA molecules: application of the extended matrix generator method

We have used the elementary generator matrices outlined in the preceding paper to examine the conformational plasticity of the nucleic acid double helix. Here we investigate kinked DNA structures made up of alternating B- and A-type helices and intrinsically curved duplexes perturbed by the intercalation of ligands. We model the B-to-A transition by the lateral translation of adjacent base pairs, and the intercalation of ligands by the vertical displacement of neighboring residues. We report a complete set of average configuration-dependent parameters, ranging from scalars (i.e., persistence lengths) to first- and second-order tensor parameters (i.e., average second moments of inertia), as well as approximations of the associated spatial distributions of the DNA and their angular correlations. The average structures of short chains (of lengths less than 100 base pairs) with local kinks or intrinsically curved sequences are essentially rigid rods. At the smallest chain lengths (10 base pairs), the kinked and curved chains exhibit similar average properties, although they are structurally perturbed compared to the standard B-DNA duplex. In contrast, at lengths of 200 base pairs, the curved and kinked chains are more compact on average and are located in a different space from the standard B- or A-DNA helix. While A-DNA is shorter and thicker than B-DNA in x-ray models, the long flexible A-DNA helix is thinner and more extended on average than its B-DNA counterpart because of more limited fluctuations in local structure. Curved polymers of 50 base pairs or longer also show significantly greater asymmetry than other DNAs (in terms of the distribution of base pairs with respect to the center of gravity of the chain). The intercalation of drugs in the curved DNA straightens and extends the smoothly deformed template. The dimensions of the average ellipsoidal boundaries defining the configurations of the intercalated polymers are roughly double those of the intrinsically curved chain. The altered proportions and orientations of these density functions reflect the changing shape and flexibility of the double helix. The calculations shed new light on the possible structural role of short A-DNA fragments in long B-type duplexes and also offer a model for understanding how GC-specific intercalative ligands can straighten naturally curved DNA. The mechanism is not immediately obvious from current models of DNA curvature, which attribute the bending of the chain to a perturbed structure in repeating tracts of A.T base pairs.

Strong sequence patterns in eukaryotic promoter regions: potential implications for DNA structure

1. Analysis of eukaryotic sequences reveals recurring trends in upstream regions. Oligomers composed of (G/C)n and (A/T)m blocks are preferentially flanked by (G/C)2 doublets on their 3' rather than on their 5' ends, that is (G/C)n(A/T)m(G/C)2 > (G/C)n+2(A/T)m. 2. These trends are stronger for larger n and smaller m. Additional trends are outlined below. 3. The trends are correlated with DNA structural parameters, in particular with twist and roll angles. 4. Generally, the trends hold if the base pair step joining the 5' (G/C)2 doublet to the (G/C)n (A/T)m oligomer is not undertwisted and is not strongly rolled into the major groove. 5. Other DNA parameters crucial for DNA-protein interactions are discussed as well.

Persistence analysis of the static and dynamical helix deformations of DNA oligonucleotides: application to the crystal structure and molecular dynamics simulation of d(CGCGAATTCGCG)2

A theory and graphical presentation for the analysis of helix structure and deformations in oligonucleotides is presented. The parameters "persistence" and "flexibility" as defined in the configurational statistics of polymers of infinite length are reformulated at the oligonucleotide level in an extension of J. A. Schellman's method [(1974) Biopolymers, Vol. 17, pp. 217-226], and used as a basis for a systematic "Persistence Analysis" of the helix deformation properties for all possible subsequences in the structure. The basis for the analysis is a set of link vectors referenced to individual base pairs, and is limited to sequences exhibiting only perturbed rod-like behavior, i.e., below the threshold for supercoiling. The present application of the method is concerned with a physical model for the angular component of bending, so the link vectors are defined as the unit components of a global helix axis obtained by the procedure "Curves" of R. Lavery and H. Sklenar [(1988) J. Biomol. Struct. Dynam., Vol. 6, pp. 63-91; (1989) ibid., Vol. 6, pp. 655-667]. A discussion of the relationship between global bending and relative orientation of base pairs is provided. Our approach is illustrated by analysis of some model oligonucleotide structures with intrinsic kinks, the crystal structure of the dodecamer d(CGCGAATTCGCG)2, and the results of two molecular dynamics simulations on this dodecamer using two variations of the GROMOS force field. The results indicate that essentially all aspects of curvature in short oligonucleotides can be determined, such as the position and orientation of each bend, the sharpness or smoothness, and the location and linearity of subsequences. In the case of molecular dynamics simulations, where a Boltzmann ensemble of structures is analyzed, the spatial extent of the deformations (flexibility) is also considered.

The electrostatic contribution to DNA base-stacking interactions

Base-stacking and phosphate-phosphate interactions in B-DNA are studied using the finite difference Poisson-Boltzmann equation. Interaction energies and dielectric constants are calculated and compared to the predictions of simple dielectric models. No extant simple dielectric model adequately describes phosphate-phosphate interactions. Electrostatic effects contribute negligibly to the sequence and conformational dependence of base-stacking interactions. Electrostatic base-stacking interactions can be adequately modeled using the Hingerty screening function. The repulsive and dispersive Lennard-Jones interactions dominate the dependence of the stacking interactions on roll, tilt, twist, and propellor. The Lennard-Jones stacking energy in ideal B-DNA is found to be essentially independent of sequence.

Analysis of the creA gene, a regulator of carbon catabolite repression in Aspergillus nidulans

The complete nucleotide sequence derived from a genomic clone and two cDNA clones of the creA gene of Aspergillus nidulans is presented. The gene contains no introns. The derived polypeptide of 415 amino acids contains two zinc fingers of the C2H2 class, frequent S(T)PXX motifs, and an alanine-rich region indicative of a DNA-binding repressor protein. The amino acid sequence of the zinc finger region has 84% similarity to the zinc finger region of Mig1, a protein involved in carbon catabolite repression in yeast cells, and it is related both to the mammalian Egr1 and Egr2 proteins and to the Wilms' tumor protein. A deletion removing the creA gene was obtained, by using in vitro techniques, in both a heterokaryon and a diploid strain but was unobtainable in a pure haploid condition. Evidence is presented suggesting that the phenotype of such a deletion, when not complemented by another creA allele, is leaky lethality allowing limited germination of the spore but not colony formation. This phenotype is far more extreme than that of any of the in vivo-generated mutations, and thus either the gene product may have an activator activity as well as a repressor function or some residual repressor function may be required for full viability.

Analysis of the creA gene, a regulator of carbon catabolite repression in Aspergillus nidulans

The complete nucleotide sequence derived from a genomic clone and two cDNA clones of the creA gene of Aspergillus nidulans is presented. The gene contains no introns. The derived polypeptide of 415 amino acids contains two zinc fingers of the C2H2 class, frequent S(T)PXX motifs, and an alanine-rich region indicative of a DNA-binding repressor protein. The amino acid sequence of the zinc finger region has 84% similarity to the zinc finger region of Mig1, a protein involved in carbon catabolite repression in yeast cells, and it is related both to the mammalian Egr1 and Egr2 proteins and to the Wilms' tumor protein. A deletion removing the creA gene was obtained, by using in vitro techniques, in both a heterokaryon and a diploid strain but was unobtainable in a pure haploid condition. Evidence is presented suggesting that the phenotype of such a deletion, when not complemented by another creA allele, is leaky lethality allowing limited germination of the spore but not colony formation. This phenotype is far more extreme than that of any of the in vivo-generated mutations, and thus either the gene product may have an activator activity as well as a repressor function or some residual repressor function may be required for full viability.

Pulse-field electrophoresis indicates full-length Mycoplasma chromosomes range widely in size

Full-size linear chromosomes were prepared from mycoplasmas by using gamma-irradiation to introduce one (on average) double-strand break in their circular chromosomes. Chromosome sizes were estimated by pulsed-field gel electrophoresis (PFGE) from the mobilities of these full-length molecules relative to DNA size references. Sizes estimated for Ureaplasma urealyticum T960 and 16 Mycoplasma species ranged from 684 kbp (M. hominis) to 1315 kbp (M. iowae). Using this sample, we found no correlation between the mobility of the full-size linear chromosomes and their G + C content. Sizes for A. laidlawii and A. hippikon were within the range expected from renaturation kinetics. PFGE size estimates are in good agreement with sizes determined by other methods, including electron microscopy, an ordered clone library, and summation of restriction fragments. Our estimates also agree with those from renaturation kinetics for both the largest and some of the smallest chromosomes, but in the intermediate size range, renaturation kinetics consistently provides lower values than PFGE or electron microscopy. Our PFGE estimates show that mycoplasma chromosomes span a continual range of sizes, with several intermediate values falling between the previously recognized large and small chromosome size clusters.

Sequence signals in eukaryotic upstream regions

Two DNA sequence elements are known to recur frequently upstream of eukaryotic polymerase II-transcribed genes. The TATAAA, at position -40, specifies the transcription initiation site. The GGCCAATCT is less frequent around -80. Sequence analysis of upstream regions reveals that the underlined yeast UAS2 consensus sequence, TGATTGGT, is also very frequent at -80 in higher polymerase II-transcribed animal sequences. The underlined CCAAT box and yeast UAS sequences are complementary. Structural analysis suggests some symmetry in their DNA structures. Upstream of the TATAAT-rich region there is an abundance of GC sequences. Analysis of nucleotide tracts indicates that these are preferentially flanked by their complementary nucleotides with a pyrimidine-purine junction, i.e., TTAN, CCGn, CnGG, TnAA. Here, I discuss DNA structural consideration in upstream regions along with protein readout of the major and minor groove information content. These sequence-structure aspects are put in the general context of protein (factors)-DNA (elements) recognition and regulation.

Self-attraction and natural curvature in null DNA

Forces of self-attraction inherent in DNA are unmasked when its ionic charge is neutralized. On the global level, self-attraction operates between segments to condense null (charge-neutralized) DNA into a segment-rich particle. Locally, self-attraction tends to contract an individual segment along its axis. If certain conditions are satisfied, the compressed segment buckles outward from the original line of the axis. Its most stable shape is then curved, or, as an extreme case, even completely folded. Buckling conditions are derived and shown to be met by DNA, thus explaining the high degree of ordered curvature and folding in the observed morphologies of condensed null DNA. The central concept employed is the buckling persistence length. It is evaluated for null DNA (40-50 bp) and agrees with experimental data (less than 60 bp). It helps in understanding the observed cooperative unit in the condensation/decondensation equilibrium (about 60 bp) and the observed size of digestion fragments unstable in the condensed phase (about 80 bp). The root-mean-square thermal compression/extension fluctuation in DNA is estimated at about 0.1 A/bp.