Modeling superhelical DNA: recent analytical and dynamic approaches

Multiscale integration of environmental stimuli in plant tropism produces complex behaviors

Plant tropism refers to the directed movement of an organ or organism in response to external stimuli. Typically, these stimuli induce hormone transport that triggers cell growth or deformation. In turn, these local cellular changes create mechanical forces on the plant tissue that are balanced by an overall deformation of the organ, hence changing its orientation with respect to the stimuli. This complex feedback mechanism takes place in a three-dimensional growing plant with varying stimuli depending on the environment. We model this multiscale process in filamentary organs for an arbitrary stimulus by explicitly linking hormone transport to local tissue deformation leading to the generation of mechanical forces and the deformation of the organ in three dimensions. We show, as examples, that the gravitropic, phototropic, nutational, and thigmotropic dynamic responses can be easily captured by this framework. Further, the integration of evolving stimuli and/or multiple contradictory stimuli can lead to complex behavior such as sun following, canopy escape, and plant twining.

Traveling Salesman Finds Random Walk: A Curve Reconstruction Algorithm for Supercoiled DNA

DNA supercoiling plays an important role in a variety of cellular processes, including transcription, replication, and DNA compaction. To fully understand these processes, we must uncover and characterize the dynamics of supercoiled DNA. However, supercoil dynamics are difficult to access because of the wide range of relevant length and timescales. In this work, we present an algorithm to reconstruct the arrangement of identical fluorescent particles distributed around a circular DNA molecule, given their three-dimensional trajectories through time. We find that this curve reconstruction problem is analogous to solving the traveling salesman problem. We demonstrate that our approach converges to the correct arrangement with a sufficiently long observation time. In addition, we show that the time required to accurately reconstruct the fluorophore arrangement is reduced by increasing the fluorophore density or reducing the level of supercoiling. This curve reconstruction algorithm, when paired with next-generation super-resolution imaging systems, could be used to access and thereby advance our understanding of supercoil dynamics.

Influence of Nucleoid-Associated Proteins on DNA Supercoiling

DNA supercoiling, where the DNA strand forms a writhe to relieve torsional stress, plays a vital role in packaging the genetic material in cells. Experiment, simulation, and theory have all demonstrated how supercoiling emerges due to the over- or underwinding of the DNA strand. Nucleoid-associated proteins (NAPs) help structure DNA in prokaryotes, yet the role that they play in the supercoiling process has not been as thoroughly investigated. We develop a coarse-grained simulation to model DNA supercoiling in the presence of proteins, providing a rigorous physical understanding of how NAPs affect supercoiling behavior. Specifically, we demonstrate how the force and torque necessary to form supercoils are affected by the presence of NAPs. NAPs that bend DNA stabilize the supercoil, thus shifting the transition between extended and supercoiled DNAs. We develop a theory to explain how NAP binding affects DNA supercoiling. This provides insight into how NAPs modulate DNA compaction via a combination of supercoiling and local protein-dependent deformations.

B7 Liquid Crystal Filament Growth in Presence of Carbon Nanotubes

Liquid crystal phases formed from bent-core mesogens have attracted much interest of the liquid crystal research community, due to the manifestation of chirality effects from achiral molecules. One of the most elusive of the bent-core phases is the B7 phase, which at its early stage often forms in a helical filament fashion. We investigate the growth of such filaments in the presence of single-walled nanotubes to elucidate possible effects on the growth dynamics and helicity of B7 helical filaments. It is found that the filament width slightly decreases in comparison to the neat B7 material, suggesting a more tightly bound structure around the nanotubes, with the nanotubes likely acting as the core of the helical filament. No effects on pitch or periodicity of the helical superstructure is observed. The filament growth velocity quickly decreases as nanotubes are added to the B7 phase, indicating that a more tightly bound structure needs a longer time of formation. An observed buckling instability is of interest as a microscopic example for the study of nonlinear dynamics theories of filaments. The present investigation is thus of general importance for nanoparticle directed growth of filaments, which has applications in biomolecular growth and high tensile strength fibres.

Hyperplectonemes: A Higher Order Compact and Dynamic DNA Self-Organization

Bacterial chromosome has a compact structure that dynamically changes its shape in response to bacterial growth rate and growth phase. Determining how chromatin remains accessible to DNA binding proteins, and transcription machinery is crucial to understand the link between genetic regulation, DNA structure, and topology. Here, we study very large supercoiled dsDNA using high-resolution characterization, theoretical modeling, and molecular dynamics calculations. We unveil a new type of highly ordered DNA organization forming in the presence of attractive DNA-DNA interactions, which we call hyperplectonemes. We demonstrate that their formation depends on DNA size, supercoiling, and bacterial physiology. We compare structural, nanomechanic, and dynamic properties of hyperplectonemes bound by three highly abundant nucleoid-associated proteins (FIS, H-NS, and HU). In all these cases, the negative supercoiling of DNA determines molecular dynamics, modulating their 3D shape. Overall, our findings provide a mechanistic insight into the critical role of DNA topology in genetic regulation.

Structural Basis for Elastic Mechanical Properties of the DNA Double Helix

In this article, we investigate the principal structural features of the DNA double helix and their effects on its elastic mechanical properties. We develop, in the pursuit of this purpose, a helical continuum model consisting of a soft helical core and two stiff ribbons wrapping around it. The proposed model can reproduce the negative twist-stretch coupling of the helix successfully as well as its global stretching, bending, and torsional rigidities measured experimentally. Our parametric study of the model using the finite element method further reveals that the stiffness of phosphate backbones is a crucial factor for the counterintuitive overwinding behavior of the duplex and its extraordinarily high torsional rigidity, the major-minor grooves augment the twist-stretch coupling, and the change of the helicity might be responsible for the transition from a negative to a positive twist-stretching coupling when a tensile force is applied to the duplex.

DNA buckling in bacteriophage cavities as a mechanism to aid virus assembly

While relatively simple biologically, bacteriophages are sophisticated biochemical machines that execute a precise sequence of events during virus assembly, DNA packaging, and ejection. These stages of the viral life cycle require intricate coordination of viral components whose structures are being revealed by single molecule experiments and high resolution (cryo-electron microscopy) reconstructions. For example, during packaging, bacteriophages employ some of the strongest known molecular motors to package DNA against increasing pressure within the viral capsid shell. Located upstream of the motor is an elaborate portal system through which DNA is threaded. A high resolution reconstruction of the portal system for bacteriophage ϕ29 reveals that DNA buckles inside a small cavity under large compressive forces. In this study, we demonstrate that DNA can also buckle in other bacteriophages including T7 and P22. Using a computational rod model for DNA, we demonstrate that a DNA buckle can initiate and grow within the small confines of a cavity under biologically-attainable force levels. The forces of DNA-cavity contact and DNA-DNA electrostatic repulsion ultimately limit cavity filling. Despite conforming to very different cavity geometries, the buckled DNA within T7 and P22 exhibits near equal volumetric energy density (∼1kT/nm(3)) and energetic cost of packaging (∼22kT). We hypothesize that a DNA buckle creates large forces on the cavity interior to signal the conformational changes to end packaging. In addition, a DNA buckle may help retain the genome prior to tail assembly through significantly increased contact area with the portal.

Buckling of a soap film spanning a flexible loop resistant to bending and twisting

A generalization of the Euler–Plateau problem to account for the energy contribution due to twisting of the bounding loop is proposed. Euler–Lagrange equations are derived in a parametrized setting and a buckling analysis is performed. A pair of dimensionless parameters govern buckling from a flat, circular ground state. While one of these is familiar from the Euler–Plateau problem, the other encompasses information about the ratio of the torsional rigidity to the bending rigidity, the twist density and the size of the loop. For sufficiently small values of the latter parameter, two separate groups of buckling modes are identified. However, for values of that parameter exceeding the critical twist density arising in Michell's study of the stability of a twisted elastic ring, only one group of buckling modes exists. Buckling diagrams indicate that a loop with greater torsional rigidity shows more resistance to transverse buckling. Additionally, a twisted flexible loop spanned by a soap film buckles at a value of the twist density less that the value at which buckling would occur if the soap film were absent.

Role of microscopic flexibility in tightly curved DNA

The genetic material in living cells is organized into complex structures in which DNA is subjected to substantial contortions. Here we investigate the difference in structure, dynamics, and flexibility between two topological states of a short (107 base pair) DNA sequence in a linear form and a covalently closed, tightly curved circular DNA form. By employing a combination of all-atom molecular dynamics (MD) simulations and elastic rod modeling of DNA, which allows capturing microscopic details while monitoring the global dynamics, we demonstrate that in the highly curved regime the microscopic flexibility of the DNA drastically increases due to the local mobility of the duplex. By analyzing vibrational entropy and Lipari–Szabo NMR order parameters from the simulation data, we propose a novel model for the thermodynamic stability of high-curvature DNA states based on vibrational untightening of the duplex. This novel view of DNA bending provides a fundamental explanation that bridges the gap between classical models of DNA and experimental studies on DNA cyclization, which so far have been in substantial disagreement.

An elastic rod model to evaluate effects of ionic concentration on equilibrium configuration of DNA in salt solution

As a coarse-gained model, a super-thin elastic rod subjected to interfacial interactions is used to investigate the condensation of DNA in a multivalent salt solution. The interfacial traction between the rod and the solution environment is determined in terms of the Young-Laplace equation. Kirchhoff's theory of elastic rod is used to analyze the equilibrium configuration of a DNA chain under the action of the interfacial traction. Two models are established to characterize the change of the interfacial traction and elastic modulus of DNA with the ionic concentration of the salt solution, respectively. The influences of the ionic concentration on the equilibrium configuration of DNA are discussed. The results show that the condensation of DNA is mainly determined by competition between the interfacial energy and elastic strain energy of the DNA itself, and the interfacial traction is one of forces that drive DNA condensation. With the change of concentration, the DNA segments will undergo a series of alteration from the original configuration to the condensed configuration, and the spiral-shape appearing in the condensed configuration of DNA is independent of the original configuration.

Perspective: Coarse-grained models for biomolecular systems

By focusing on essential features, while averaging over less important details, coarse-grained (CG) models provide significant computational and conceptual advantages with respect to more detailed models. Consequently, despite dramatic advances in computational methodologies and resources, CG models enjoy surging popularity and are becoming increasingly equal partners to atomically detailed models. This perspective surveys the rapidly developing landscape of CG models for biomolecular systems. In particular, this review seeks to provide a balanced, coherent, and unified presentation of several distinct approaches for developing CG models, including top-down, network-based, native-centric, knowledge-based, and bottom-up modeling strategies. The review summarizes their basic philosophies, theoretical foundations, typical applications, and recent developments. Additionally, the review identifies fundamental inter-relationships among the diverse approaches and discusses outstanding challenges in the field. When carefully applied and assessed, current CG models provide highly efficient means for investigating the biological consequences of basic physicochemical principles. Moreover, rigorous bottom-up approaches hold great promise for further improving the accuracy and scope of CG models for biomolecular systems.

Bacterial promoter repression by DNA looping without protein-protein binding competition

The Escherichia coli lactose operon provides a paradigm for understanding gene control by DNA looping where the lac repressor (LacI) protein competes with RNA polymerase for DNA binding. Not all promoter loops involve direct competition between repressor and RNA polymerase. This raises the possibility that positioning a promoter within a tightly constrained DNA loop is repressive per se, an idea that has previously only been considered in vitro. Here, we engineer living E. coli bacteria to measure repression due to promoter positioning within such a tightly constrained DNA loop in the absence of protein–protein binding competition. We show that promoters held within such DNA loops are repressed ∼100-fold, with up to an additional ∼10-fold repression (∼1000-fold total) dependent on topological positioning of the promoter on the inner or outer face of the DNA loop. Chromatin immunoprecipitation data suggest that repression involves inhibition of both RNA polymerase initiation and elongation. These in vivo results show that gene repression can result from tightly looping promoter DNA even in the absence of direct competition between repressor and RNA polymerase binding.

Structural ensemble and dynamics of toroidal-like DNA shapes in bacteriophage ϕ29 exit cavity

In the bacteriophage ϕ29, DNA is packed into a preassembled capsid from which it ejects under high pressure. A recent cryo-EM reconstruction of ϕ29 revealed a compact toroidal DNA structure (30-40 basepairs) lodged within the exit cavity formed by the connector-lower collar protein complex. Using multiscale models, we compute a detailed structural ensemble of intriguing DNA toroids of various lengths, all highly compatible with experimental observations. In particular, coarse-grained (elastic rod) and atomistic (molecular dynamics) models predict the formation of DNA toroids under significant compression, a largely unexplored state of DNA. Model predictions confirm that a biologically attainable compressive force of 25 pN sustains the toroid and yields DNA electron density maps highly consistent with the experimental reconstruction. The subsequent simulation of dynamic toroid ejection reveals large reactions on the connector that may signal genome release.

Molecular modeling of nucleic acid structure

This unit is the first in a series of four units covering the analysis of nucleic acid structure by molecular modeling. The unit provides an overview of the computer simulation of nucleic acids. Topics include the static structure model, computational graphics and energy models, the generation of an initial model, and characterization of the overall three-dimensional structure.

Exploring polymorphisms in B-DNA helical conformations

The traditional mesoscopic paradigm represents DNA as a series of base-pair steps whose energy response to equilibrium perturbations is elastic, with harmonic oscillations (defining local stiffness) around a single equilibrium conformation. In addition, base sequence effects are often analysed as a succession of independent XpY base-pair steps (i.e. a nearest-neighbour (NN) model with only 10 unique cases). Unfortunately, recent massive simulations carried out by the ABC consortium suggest that the real picture of DNA flexibility may be much more complex. The paradigm of DNA flexibility therefore needs to be revisited. In this article, we explore in detail one of the most obvious violations of the elastic NN model of flexibility: the bimodal distributions of some helical parameters. We perform here an in-depth statistical analysis of a very large set of MD trajectories and also of experimental structures, which lead to very solid evidence of bimodality. We then suggest ways to improve mesoscopic models to account for this deviation from the elastic regime.

DNA modeling reveals an extended lac repressor conformation in classic in vitro binding assays

Protein-mediated DNA looping, such as that induced by the lactose repressor (LacI) of Escherichia coli, is a well-known gene regulation mechanism. Although researchers have given considerable attention to DNA looping by LacI, many unanswered questions about this mechanism, including the role of protein flexibility, remain. Recent single-molecule observations suggest that the two DNA-binding domains of LacI are capable of splaying open about the tetramerization domain into an extended conformation. We hypothesized that if recent experiments were able to reveal the extended conformation, it is possible that such structures occurred in previous studies as well. In this study, we tested our hypothesis by reevaluating two classic in vitro binding assays using a computational rod model of DNA. The experiments and computations evaluate the looping of both linear DNA and supercoiled DNA minicircles over a broad range of DNA interoperator lengths. The computed energetic minima align well with the experimentally observed interoperator length for optimal loop stability. Of equal importance, the model reveals that the most stable loops for linear DNA occur when LacI adopts the extended conformation. In contrast, for DNA minicircles, optimal stability may arise from either the closed or the extended protein conformation depending on the degree of supercoiling and the interoperator length.

A multiscale dynamic model of DNA supercoil relaxation by topoisomerase IB

In this study, we report what we believe to be the first multiscale simulation of the dynamic relaxation of DNA supercoils by human topoisomerase IB (topo IB). We leverage our previous molecular dynamics calculations of the free energy landscape describing the interaction between a short DNA fragment and topo IB. Herein, this landscape is used to prescribe boundary conditions for a computational, elastodynamic continuum rod model of a long length of supercoiled DNA. The rod model, which accounts for the nonlinear bending, twisting, and electrostatic interaction of the (negatively charged) DNA backbone, is extended to include the hydrodynamic drag induced by the surrounding physiological buffer. Simulations for a 200-bp-long DNA supercoil in complex with topo IB reveal a relaxation timescale of ∼0.1-1.0 μs. The relaxation follows a sequence of cascading reductions in the supercoil linking number (Lk), twist (Tw), and writhe (Wr) that follow companion cascading reductions in the supercoil elastic and electrostatic energies. The novel (to our knowledge) multiscale modeling method may enable simulations of the entire experimental setup that measures DNA supercoiling and relaxation via single molecule magnetic trapping.

Gene repression by minimal lac loops in vivo

The inflexibility of double-stranded DNA with respect to bending and twisting is well established in vitro. Understanding apparent DNA physical properties in vivo is a greater challenge. Here, we exploit repression looping with components of the Escherichia coli lac operon to monitor DNA flexibility in living cells. We create a minimal system for testing the shortest possible DNA repression loops that contain an E. coli promoter, and compare the results to prior experiments. Our data reveal that loop-independent repression occurs for certain tight operator/promoter spacings. When only loop-dependent repression is considered, fits to a thermodynamic model show that DNA twisting limits looping in vivo, although the apparent DNA twist flexibility is 2- to 4-fold higher than in vitro. In contrast, length-dependent resistance to DNA bending is not observed in these experiments, even for the shortest loops constraining <0.4 persistence lengths of DNA. As observed previously for other looping configurations, loss of the nucleoid protein heat unstable (HU) markedly disables DNA looping in vivo. Length-independent DNA bending energy may reflect the activities of architectural proteins and the structure of the DNA topological domain. We suggest that the shortest loops are formed in apical loops rather than along the DNA plectonemic superhelix.

Cosserat curve model for superelasticity of helices

Superelasticity behavior of helices has been the focus of recent research in micro-/nano-engineering, while the traditional Kirchhoff rod model restricts itself in the bending and torsion conditions. With the aid of the concept of a Cosserat curve, a novel theoretical basis has been established for statics and dynamics of helices with essential extension and shear, which is able to quantitatively analyze the superelastic mechanical properties. Except for a good agreement with the experimental observation, numerical solutions have shown that we cannot only predict two important properties of the superelasticity characteristics: the breaking force and the stretch of the coil wire under the axial loading, but also precisely describe and explain the Hooke's constant and torque in the entire stretching and breaking processes. The present work has provided useful information for the future experimental investigation on superelasticity as well as its application in meta-/quantum devices.

Computational analysis of looping of a large family of highly bent DNA by LacI

Sequence-dependent intrinsic curvature of DNA influences looping by regulatory proteins such as LacI and NtrC. Curvature can enhance stability and control shape, as observed in LacI loops formed with three designed sequences with operators bracketing an A-tract bend. We explore geometric, topological, and energetic effects of curvature with an analysis of a family of highly bent sequences, using the elastic rod model from previous work. A unifying straight-helical-straight representation uses two phasing parameters to describe sequences composed of two straight segments that flank a common helically supercoiled segment. We exercise the rod model over this two-dimensional space of phasing parameters to evaluate looping behaviors. This design space is found to comprise two subspaces that prefer parallel versus anti-parallel binding topologies. The energetic cost of looping varies from 4 to 12 kT. Molecules can be designed to yield distinct binding topologies as well as hyperstable or hypostable loops and potentially loops that can switch conformations. Loop switching could be a mechanism for control of gene expression. Model predictions for linking numbers and sizes of LacI-DNA loops can be tested using multiple experimental approaches, which coupled with theory could address whether proteins or DNA provide the observed flexibility of protein-DNA loops.

Helical coherence of DNA in crystals and solution

The twist, rise, slide, shift, tilt and roll between adjoining base pairs in DNA depend on the identity of the bases. The resulting dependence of the double helix conformation on the nucleotide sequence is important for DNA recognition by proteins, packaging and maintenance of genetic material, and other interactions involving DNA. This dependence, however, is obscured by poorly understood variations in the stacking geometry of the same adjoining base pairs within different sequence contexts. In this article, we approach the problem of sequence-dependent DNA conformation by statistical analysis of X-ray and NMR structures of DNA oligomers. We evaluate the corresponding helical coherence length--a cumulative parameter quantifying sequence-dependent deviations from the ideal double helix geometry. We find, e.g. that the solution structure of synthetic oligomers is characterized by 100-200 A coherence length, which is similar to approximately 150 A coherence length of natural, salmon-sperm DNA. Packing of oligomers in crystals dramatically alters their helical coherence. The coherence length increases to 800-1200 A, consistent with its theoretically predicted role in interactions between DNA at close separations.

Molecular modeling of nucleic acid structure

This unit is the first in a series of four units covering the analysis of nucleic acid structure by molecular modeling. This unit provides an overview of computer simulation of nucleic acids. Topics include the static structure model, computational graphics and energy models, generation of an initial model, and characterization of the overall three-dimensional structure.

Intrinsic curvature of DNA influences LacR-mediated looping

Protein-mediated DNA looping is a common mechanism for regulating gene expression. Loops occur when a protein binds to two operators on the same DNA molecule. The probability of looping is controlled, in part, by the basepair sequence of inter-operator DNA, which influences its structural properties. One structural property is the intrinsic or stress-free curvature. In this article, we explore the influence of sequence-dependent intrinsic curvature by exercising a computational rod model for the inter-operator DNA as applied to looping of the LacR-DNA complex. Starting with known sequences for the inter-operator DNA, we first compute the intrinsic curvature of the helical axis as input to the rod model. The crystal structure of the LacR (with bound operators) then defines the requisite boundary conditions needed for the dynamic rod model that predicts the energetics and topology of the intervening DNA loop. A major contribution of this model is its ability to predict a broad range of published experimental data for highly bent (designed) sequences. The model successfully predicts the loop topologies known from fluorescence resonance energy transfer measurements, the linking number distribution known from cyclization assays with the LacR-DNA complex, the relative loop stability known from competition assays, and the relative loop size known from gel mobility assays. In addition, the computations reveal that highly curved sequences tend to lower the energetic cost of loop formation, widen the energy distribution among stable and meta-stable looped states, and substantially alter loop topology. The inclusion of sequence-dependent intrinsic curvature also leads to nonuniform twist and necessitates consideration of eight distinct binding topologies from the known crystal structure of the LacR-DNA complex.

Modeling DNA loops using the theory of elasticity

An elastic rod model of a protein-bound DNA loop is adapted for application in multi-scale simulations of protein-DNA complexes. The classical Kirchhoff system of equations which describes the equilibrium structure of the elastic loop is modified to account for the intrinsic twist and curvature, anisotropic bending properties, and electrostatic charge of DNA. The effects of bending anisotropy and electrostatics are studied for the DNA loop clamped by the lac repressor protein. For two possible lengths of the loop, several topologically different conformations are predicted and extensively analyzed over the broad range of model parameters describing DNA bending and electrostatic properties. The scope and applications of the model in already accomplished and in future multi-scale studies of protein-DNA complexes are discussed.

Structural dynamics of the lac repressor-DNA complex revealed by a multiscale simulation

A multiscale simulation of a complex between the lac repressor protein (LacI) and a 107-bp-long DNA segment is reported. The complex between the repressor and two operator DNA segments is described by all-atom molecular dynamics; the size of the simulated system comprises either 226,000 or 314,000 atoms. The DNA loop connecting the operators is modeled as a continuous elastic ribbon, described mathematically by the nonlinear Kirchhoff differential equations with boundary conditions obtained from the coordinates of the terminal base pairs of each operator. The forces stemming from the looped DNA are included in the molecular dynamics simulations; the loop structure and the forces are continuously recomputed because the protein motions during the simulations shift the operators and the presumed termini of the loop. The simulations reveal the structural dynamics of the LacI-DNA complex in unprecedented detail. The multiple domains of LacI exhibit remarkable structural stability during the simulation, moving much like rigid bodies. LacI is shown to absorb the strain from the looped DNA mainly through its mobile DNA-binding head groups. Even with large fluctuating forces applied, the head groups tilt strongly and keep their grip on the operator DNA, while the remainder of the protein retains its V-shaped structure. A simulated opening of the cleft of LacI by 500-pN forces revealed the interactions responsible for locking LacI in the V-conformation.

Stability results for constrained calculus of variations problems: an analysis of the twisted elastic loop

Problems with a variational structure are ubiquitous throughout the physical sciences and have a distinguished scientific history. Constrained variational problems have been much less studied, particularly the theory of stability, which determines which solutions are physically realizable. In this paper, we develop stability exchange results appropriate for parameter-dependent calculus of variations problems with two particular features: either the parameter appears in the boundary conditions, or there are isoperimetric constraints. In particular, we identify an associated distinguished bifurcation diagram, which encodes the direction of stability exchange at folds. We apply the theory to a twisted elastic loop, which can naturally be formulated as a calculus of variations problem with both isoperimetric constraints and parameter-dependent boundary conditions. In combination with a perturbation expansion that classifies certain pitchfork bifurcations as sub- or super-critical, the distinguished diagram for the twisted loop provides a classification of the stability properties of all equilibria. In particular, an unanticipated sensitive dependence of stability properties on the ratio of twisting to bending stiffness is revealed.

Structural basis for cooperative DNA binding by CAP and lac repressor

Catabolite gene activator protein (CAP) and lac repressor (LR) are celebrated transcription-regulating proteins that bind to DNA cooperatively forming a ternary complex with the promoter loop. Here we present a multiscale model of the ternary complex derived from crystal structures of the proteins and a continuous structure of the DNA loop built using the theory of elasticity. We predict that the loop is underwound in the binary complex with the LR, whereas in the ternary complex with the LR and CAP, the loop is overwound and extended due to an upstream relocation of a DNA binding hand of LR. The computed relocation distance matches the experimental observations and the energy balance of the system explains the cooperativity effect. Using the multiscale approach, we build an all-atom model of the ternary complex that suggests a series of further experimental investigations.

Methods for determining stability in continuum elastic-rod models of DNA

Elastic-rod models of DNA have offered an alternative method for studying the macroscopic properties of the molecule. An essential component of the modelling effort is to identify the biologically accessible, or stable, solutions. The underlying variational structure of the elastic-rod model can be exploited to derive methods that identify stable equilibrium configurations. We present two methods for determining the stability of the equilibria of elastic-rod models: the conjugate-point method and the distinguished-diagram method. Additionally, we apply these methods to two intrinsically curved DNA molecules: a DNA filament with an A-tract bend and a DNA minicircle with a catabolite gene activator protein binding site. The stable solutions of these models provide visual insight into the three-dimensional structure of the DNA molecules.

Symmetric equilibria of a thin elastic rod with self-contacts

The thin elastic rod is a conventional model for the mesoscale structure of DNA. The symmetric multi-leafed closed equilibria with one multiple self-contact point are examined. The contact force is assumed to be pointwise and frictionless. The bifurcation diagram is presented for first branches of two- and three-leafed families including trefoil knots.

Modelling DNA loops using continuum and statistical mechanics

The classical Kirchhoff elastic-rod model applied to DNA is extended to account for sequence-dependent intrinsic twist and curvature, anisotropic bending rigidity, electrostatic force interactions, and overdamped Brownian motion in a solvent. The zero-temperature equilibrium rod model is then applied to study the structural basis of the function of the lac repressor protein in the lac operon of Escherichia coli. The structure of a DNA loop induced by the clamping of two distant DNA operator sites by lac repressor is investigated and the optimal geometries for the loop of length 76 bp are predicted. Further, the mimicked binding of catabolite gene activator protein (CAP) inside the loop provides solutions that might explain the experimentally observed synergy in DNA binding between the two proteins. Finally, a combined Monte Carlo and Brownian dynamics solver for a worm-like chain model is described and a preliminary analysis of DNA loop-formation kinetics is presented.

DNA basepair step deformability inferred from molecular dynamics simulations

The sequence-dependent DNA deformability at the basepair step level was investigated using large-scale atomic resolution molecular dynamics simulation of two 18-bp DNA oligomers: d(GCCTATAAACGCCTATAA) and d(CTAGGTGGATGACTCATT). From an analysis of the structural fluctuations, the harmonic potential energy functions for all 10 unique steps with respect to the six step parameters have been evaluated. In the case of roll, three distinct groups of steps have been identified: the flexible pyrimidine-purine (YR) steps, intermediate purine-purine (RR), and stiff purine-pyrimidine (RY). The YR steps appear to be the most flexible in tilt and partially in twist. Increasing stiffness from YR through RR to RY was observed for rise, whereas shift and slide lack simple trends. A proposed measure of the relative importance of couplings identifies the slide-rise, twist-roll, and twist-slide couplings to play a major role. The force constants obtained are of similar magnitudes to those based on a crystallographic ensemble. However, the current data have a less complicated and less pronounced sequence dependence. A correlation analysis reveals concerted motions of neighboring steps and thus exposes limitations in the dinucleotide model. The comparison of DNA deformability from this and other studies with recent quantum-chemical stacking energy calculations suggests poor correlation between the stacking and flexibility.

Mechanical properties of nanosprings

Nanostructures (nanotubes, nanowires, etc.) have been the object of intense theoretical and experimental investigations in recent years. Among these structures, helical nanosprings or nanocoils have attracted particular interest due to their special mechanical properties. In this work, we investigated structural properties of nanosprings in the Kirchhoff rod model. We derived expressions that can be used experimentally to obtain nanospring Young's modulus and Poisson's ratio values. Our results also might explain why the presence of catalytic particles is so important in nanostructure growth.

Local sequential minimization of double stranded B-DNA using Monte Carlo annealing

A software algorithm has been developed to investigate the folding process in B-DNA structures in vacuum under a simple and accurate force field. This algorithm models linear double stranded B-DNA sequences based on a local, sequential minimization procedure. The original B-DNA structures were modeled using initial nucleotide structures taken from the Brookhaven database. The models contain information at the atomic level allowing one to investigate as accurately as possible the structure and characteristics of the resulting DNA structures. A variety of DNA sequences and sizes were investigated containing coding and non-coding, random and real, homogeneous or heterogeneous sequences in the range of 2 to 40 base pairs. The force field contains terms such as angle bend, Lennard-Jones, electrostatic interactions and hydrogen bonding which are set up using the Dreiding II force field and defined to account for the helical parameters such as twist, tilt and rise. A close comparison was made between this local minimization algorithm and a global one (previously published) in order to find out advantages and disadvantages of the different methods. From the comparison, this algorithm gives better and faster results than the previous method, allowing one to minimize larger DNA segments. DNA segments with a length of 40 bases need approximately 4 h, while 2.5 weeks are needed with the previous method. After each minimization the angles between phosphate-oxygen-carbon A1, the oxygen-phosphate-oxygen A2 and the average helical twists were calculated. From the generated fragments it was found that the bond angles are A1=150 degrees +/-2 degrees and A2=130 degrees +/-10 degrees, while the helical twist is 36.6 degrees +/-2 degrees in the A strand and A1=150 degrees +/-6 degrees and A2=130+/-6 degrees with helical twist 39.6 degrees +/-2 degrees in the B strand for the DNA segment with the same sequence as the Dickerson dodecamer.

Link, twist, energy, and the stability of DNA minicircles

We describe how the stability properties of DNA minicircles can be directly read from plots of various biologically intuitive quantities along families of equilibrium configurations. Our conclusions follow from extensions of the mathematical theory of distinguished bifurcation diagrams that are applied within the specific context of an elastic rod model of minicircles. Families of equilibria arise as a twisting angle alpha is varied. This angle is intimately related to the continuously varying linking number Lk for nicked DNA configurations that is defined as the sum of Twist and Writhe. We present several examples of such distinguished bifurcation diagrams involving plots of the energy E, linking number Lk, and a twist moment m3, along families of cyclized equilibria of both intrinsically straight and intrinsically curved DNA fragments.

Effect of DNA superhelicity and bound proteins on mechanistic aspects of the Hin-mediated and Fis-enhanced inversion

Using a recently developed inhomogeneous, macroscopic model for long DNA bound to proteins, we examine topological and geometric aspects of DNA/protein structures and dynamics on various stages of the Hin inversion pathway. This biological reaction involves exchange of DNA in a synaptic complex that brings together several DNA sites bound to Hin dimers as well as Fis enhancers. Brownian dynamics simulations in the millisecond timescale allow us to follow and analyze the DNA/protein dynamics trajectories and to examine the effects of DNA superhelicity and protein binding on various reaction steps. Analysis of the generated kinetic pathways helps explain mechanistic aspects regarding the process by which two or three protein-bound DNA sites come to close spatial proximity and show that how topological selectivity (two trapped supercoils), enhancer binding, and properties of supercoiled DNA play critical roles in regulating the inversion reaction. Specifically, a critical amount of DNA superhelicity (e.g., |sigma| > 0.02) leads to an optimal interplay for the first reaction step-two-site juxtaposition-between large-scale random rearrangements of Hin-bound DNA and local slithering within branches of plectonemes. The three-site juxtaposition, the second step, is significantly accelerated by the presence of an enhancer protein that, due to severe local bending, also alters juxtaposition mechanisms, especially for superhelical density magnitude greater than around 0.04.

Modeling the dynamics of the solvated SL1 domain of HIV-1 genomic RNA

A mode-coupling solution of the Smoluchowski diffusion equation (MCD theory), designed to describe the dynamics of wobbling macromolecules in water, is applied to a macromolecular bead model including water beads in the nearest layers. The necessary statistical averages are evaluated by time averaging along a molecular dynamics (MD) trajectory where both solute and water are introduced as atomistic models. The cross peaks in (1)H nuclear Overhauser effect spectroscopy (NOESY) NMR spectra that are routinely measured to determine biological structures are here calculated for the mutated 23 nucleotides stem-loop fragment of the SL1 domain in the HIV-1(Lai) genomic RNA. The calculations are in acceptable agreement with experiments without requiring any screening of the hydrodynamic interactions. The screening of hydrodynamics was necessary in previous MCD calculations obtained by using the same full atomistic MD trajectory, but a nonsolvated frictional model.

Time-lapse imaging of conformational changes in supercoiled DNA by scanning force microscopy

Most of the scanning force microscopy (SFM) images of supercoiled DNA on untreated mica thus far reported have not shown tight plectonemic structure seen by electron microscopy, but instead less coiled molecules and sometimes a partly "condensed" state with intimate chain-chain interactions. By observing time-lapse images of conformational changes of DNA induced by decreasing ionic strength of imaging buffer in solution SFM, we could show that the process of water rinsing, an indispensable step for preparation of dried samples, may be responsible for some of the conformational anomalies in the images previously reported. We have studied several protocols to observe supercoiled DNA molecules by SFM and discuss the merits and the demerits. Images obtained following uranyl acetate treatment may be ideal for the detection of DNA damage, as the supercoiled and nicked forms are easily distinguishable.

DNA folding: structural and mechanical properties of the two-angle model for chromatin

We present a theoretical analysis of the structural and mechanical properties of the 30-nm chromatin fiber. Our study is based on the two-angle model introduced by Woodcock et al. (Woodcock, C. L., S. A. Grigoryev, R. A. Horowitz, and N. Whitaker. 1993. Proc. Natl. Acad. Sci. USA. 90:9021-9025) that describes the chromatin fiber geometry in terms of the entry-exit angle of the nucleosomal DNA and the rotational setting of the neighboring nucleosomes with respect to each other. We analytically explore the different structures that arise from this building principle, and demonstrate that the geometry with the highest density is close to the one found in native chromatin fibers under physiological conditions. On the basis of this model we calculate mechanical properties of the fiber under stretching. We obtain expressions for the stress-strain characteristics that show good agreement with the results of recent stretching experiments (Cui, Y., and C. Bustamante. 2000. Proc. Natl. Acad. Sci. USA. 97:127-132) and computer simulations (Katritch, V., C. Bustamante, and W. K. Olson. 2000. J. Mol. Biol. 295:29-40), and which provide simple physical insights into correlations between the structural and elastic properties of chromatin.

Computational modeling predicts the structure and dynamics of chromatin fiber

BACKGROUND

The compact form of the chromatin fiber is a critical regulator of fundamental processes such as transcription and replication. These reactions can occur only when the fiber is unraveled and the DNA strands contained within are exposed to interact with nuclear proteins. While progress on identifying the biochemical mechanisms that control localized folding and hence govern access to genetic information continues, the internal structure of the chromatin fiber, let alone the structural pathways for folding and unfolding, remain unknown.

RESULTS

To offer structural insights into how this nucleoprotein complex might be organized, we present a macroscopic computer model describing the mechanics of the chromatin fiber on the polymer level. We treat the core particles as electrostatically charged disks linked via charged elastic DNA segments and surrounded by a microionic hydrodynamic solution. Each nucleosome unit is represented by several hundred charges optimized so that the effective Debye-Hückel electrostatic field matches the field predicted by the nonlinear Poisson-Boltzmann equation. On the basis of Brownian dynamics simulations, we show that oligonucleosomes condense and unfold in a salt-dependent manner analogous to the chromatin fiber.

CONCLUSIONS

Our predicted chromatin model shows good agreement with experimental diffusion coefficients and small-angle X-ray scattering data. A fiber of width 30 nm, organized in a compact helical zigzag pattern with about 4 nucleosomes per 10 nm, naturally emerges from a repeating nucleosome folding motif. This fiber has a cross-sectional radius of gyration of R(c) = 8.66 nm, in close agreement with corresponding values for rat thymus and chicken erythrocyte chromatin (8.82 and 8.5 nm, respectively).

Near equilibrium dynamics of nonhomogeneous Kirchhoff filaments in viscous media

We study the near equilibrium dynamics of nonhomogeneous elastic filaments in viscous media using the Kirchhoff model of rods. Viscosity is incorporated in the model as an external force, which we approximate by the resistance felt by an infinite cylinder immersed in a slowly moving fluid. We use the recently developed method of Goriely and Tabor [Phys. Rev. Lett. 77, 3537 (1996); Physica D 105, 20 (1997); 105, 45 (1997)] to study the dynamics in the vicinity of the simplest equilibrium solution for a closed rod with nonhomogeneous distribution of mass, namely, the planar ring configuration. We show that small variations of the mass density along the rod are sufficient to couple the symmetric modes of the homogeneous rod problem, producing asymmetric deformations that modify substantially the dynamical coiling, even at quite low Reynolds number. The higher-density segments of the rod tend to become more rigid and less coiled. We comment on possible applications to DNA.

DNA rings with multiple energy minima

Within the context of DNA rings, we analyze the relationship between intrinsic shape and the existence of multiple stable equilibria, either nicked or cyclized with the same link. A simple test, based on a perturbation expansion of symmetry breaking within a continuum elastic rod model, provides good predictions of the occurrence of such multiple equilibria. The reliability of these predictions is verified by direct computation of nicked and cyclized equilibria for several thousand DNA minicircles with lengths of 200 and 900 bp. Furthermore, our computations of equilibria for nicked rings predict properties of the equilibrium distribution of link, as calculated by much more computationally intensive Monte Carlo simulations.

Modeling DNA deformations

Recent developments have been made in modeling double-helical DNA at four levels of three-dimensional structure: the all-atom level, whereby an oligonucleotide duplex is surrounded by a shroud of solvent molecules; the base-pair level, with explicit backbone atoms; the mesoscopic level, that is, a few hundred base pairs, with the local duplex conformation described by knowledge-based harmonic energy functions; and the scale of several thousand nucleotides, with the duplex described as an ideal elastic rod. Predictions of the sequence-dependent bending and twisting of the double helix, as well as solvent- and force-induced B-->A and over-stretching conformational transitions, are compared with experimental data. These subtle conformational changes are critical to the functioning of the double helix, including its packaging in the close confines of the cell, the mutual fit of DNA and protein in nucleoprotein complexes, and the effective recognition of base pairs in recombination and transcription.

Mathematical modelling of interwound DNA supercoils

Small loops of DNA are affected by a variety of enzymes which remove turns of twist relative to the underlying double-helical structure. The molecule adopts a complex three-dimensional shape known as a supercoil in order to relieve the resulting internal stresses. This article describes an approach to modelling the overall shape of the supercoiled structure using elastic rod theory, which leads to simple expressions for predicting the shape of the structure. Predictions on the number of crossings in the balanced ply and the length of the end loops are compared to data in the literature and show reasonable agreement. The effect of the charged phosphate groups along the backbone of the DNA on the resulting supercoiled shape are also examined, and it is shown that this shape is very sensitive to the ionic concentration of the solution.