Conformations of closed DNA

Bistability of a helical filament confined on a cylinder

The natural configuration of an intrinsically curved and twisted filament is uniquely a helix so that it can be referred to as a helical filament. We find that confining a helical filament on a cylinder can create a bistable state. When c_{0}R=0.5, where c_{0} is the intrinsic curvature of filament and R is the radius of cylinder, the phase diagram for the stability of a helix contains three regimes. Regime I has a small intrinsic twisting rate (ITR) and exhibits a bistable state which consists of two isoenergic helices. In regime II, the filament has a moderate ITR and the bistable state consists of a metastable low-pitch helix and a stable nonhelix. In regime III, the helix is unstable, owing to a large ITR. A similar phenomenon occurs when c_{0}R∼0.5. Monte Carlo simulation confirms these conclusions and indicates further that there are bistable nonhelices in regime III. This bistable system offers a prospective green material since the wide range of parameters and distinctive configurations for bistable states favor its realization and application.

Bistability induced by a spontaneous twisting rate for a two-dimensional intrinsically curved filament

We find that a moderate intrinsic twisting rate (ITR) can induce a bistable state for a force-free two-dimensional intrinsically curved filament. There are two different configurations of equal energy in a bistable state so that the filament is clearly different from its three-dimensional counterpart. The smaller the ITR or the larger the intrinsic curvature (IC), the clearer the distinction between two isoenergetic configurations and the longer the filament. In bistable states, the relationship between length and ITR is approximately a hyperbola and relationship between IC and critical ITR is approximately linear. Thermal fluctuation can result in a shift between two isoenergetic configurations, but large bending and twisting rigidities can prevent the shift and maintain the filament in one of these two configurations. Moreover, a filament can have a metastable state and at a finite temperature such a filament has the similar property as that of a filament with bistable state.

Random motion theory of an optical vortex in nonlinear birefringent media

A theoretical study is presented for the random aspect of an optical vortex inherent in the nonlinear birefringent Kerr effect, which is called the optical spin vortex. We start with the two-component nonlinear Schrödinger equation. The vortex is inherent in the spin texture caused by an anisotropy of the dielectric tensor, for which the role of spin is played by the Stokes vector (or pseudospin). The evolutional equation is derived for the vortex center coordinate using the effective Lagrangian of the pseudospin field. This is converted to the Langevin equation in the presence of the fluctuation together with the dissipation. The corresponding Fokker-Planck equation is derived and analytically solved for a particular form of the birefringence inspired from the Faraday effect. The main consequence is that the relaxation distance for the distribution function is expressed by the universal constant in the Faraday effect and the size of optical vortex. The result would provide a possible clue for future experimental study in polarization optics from a stochastic aspect.

Mechanical property of the helical configuration for a twisted intrinsically straight biopolymer

We explore the effects of two typical torques on the mechanical property of the helical configuration for an intrinsically straight filament or biopolymer either in three-dimensional space or on a cylinder. One torque is parallel to the direction of a uniaxial applied force, and is coupled to the cross section of the filament. We obtain some algebraic equations for the helical configuration and find that the boundary conditions are crucial. In three-dimensional space, we show that the extension is always a monotonic function of the applied force. On the other hand, for a filament confined on a cylinder, the twisting rigidity and torque coupled to the cross section are irrelevant in forming a helix if the filament is isotropic and under free boundary condition. However, the twisting rigidity and the torque coupled to the cross section become crucial when the Euler angle at two ends of the filament are fixed. Particularly, the extension of a helix can subject to a first-order transition so that in such a condition a biopolymer can act as a switch or sensor in some biological processes. We also present several phase diagrams to provide the conditions to form a helix.

Knot soliton in DNA and geometric structure of its free-energy density

In general, the geometric structure of DNA is characterized using an elastic rod model. The Landau model provides us a new theory to study the geometric structure of DNA. By using the decomposition of the arc unit in the helical axis of DNA, we find that the free-energy density of DNA is similar to the free-energy density of a two-condensate superconductor. By using the φ-mapping topological current theory, the torus knot soliton hidden in DNA is demonstrated. We show the relation between the geometric structure and free-energy density of DNA and the Frenet equations in differential geometry theory are considered. Therefore, the free-energy density of DNA can be expressed by the curvature and torsion of the helical axis.

DNA as a one-dimensional chiral material. II. Dynamics of the structural transition between B form and Z form

We analyze the dynamics of structural transitions between normal right-handed B form and unusual left-handed Z form for a linear DNA molecule. The dynamics under the external torque in physiological buffer is modeled by a Langevin equation, with the potential term given by the authors previously [Phys. Rev. E 84, 021926 (2011)]. With this model, we first simulate the relaxation processes around B-form structure after sudden changes of the external torques, where slow relaxation ~t(-1/2) as a function of the elapsed time t is observed. Then, the dynamics of structural transition from Z form to B form is computed under various external torque strength. For small external torques, the transition proceeds via nucleation and the growth, while for higher torques, Z-form structure becomes unstable, and the transition mechanism is switched to a spinodal-like process. These numerical results are qualitatively understood by simple phenomenological arguments.

Biopolymer elasticity: Mechanics and thermal fluctuations

We present an analytical study of the role of thermal fluctuations in shaping molecular elastic properties of semiflexible polymers. Our study interpolates between mechanics and statistical mechanics in a controlled way and shows how thermal fluctuations modify the elastic properties of biopolymers. We present a study of the minimum-energy configurations with explicit expressions for their energy and writhe and plots of the extension versus link for these configurations and a study of fluctuations around the local minima of energy and approximate analytical formulas for the free energy of stretched twisted polymers. The central result of our study is a closed-form expression for the leading thermal fluctuation correction to the free energy around the nonperturbative writhing family solution for the configuration of a biopolymer. From the derived formulas, the predictions of the wormlike chain model for molecular elasticity can be worked out for a comparison against numerical simulations and experiments.

DNA as a one-dimensional chiral material: application to the structural transition between B form and Z form

A dynamical model is presented for chiral change in DNA molecules. The model is an extension of the conventional elastic model, which incorporates the structure of base pairs and uses a spinor representation for the DNA configuration together with a gauge principle. Motivated by a recent experiment reporting chiral transitions between right-handed B-DNA and left-handed Z-DNA [Lee et al., Proc. Natl. Acad. Sci. (USA) 107, 4985 (2010)], we analyze the free energy for the particular case of linear DNA with an externally applied torque. The model shows that there exists, at low temperature, a rapid structural change depending on the torque exerted on the DNA, which causes switching in B and Z domain sizes. This can explain the frequent switches of DNA extension observed in experiments.

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.

Effect of bending anisotropy on the 3D conformation of short DNA loops

The equilibrium three dimensional shape of relatively short loops of DNA is studied using an elastic model that takes into account anisotropy in bending rigidities. Using a reasonable estimate for the anisotropy, it is found that cyclized DNA with lengths that are not integer multiples of the pitch take on nontrivial shapes that involve bending out of planes and formation of kinks. The effect of sequence inhomogeneity on the shape of DNA is addressed, and shown to enhance the geometrical features. These findings could shed some light on the role of DNA conformation in protein-DNA interactions.

Writhe formulas and antipodal points in plectonemic DNA configurations

The linking and writhing numbers are key quantities when characterizing the structure of a piece of supercoiled DNA. Defined as double integrals over the shape of the double helix, these numbers are not always straightforward to compute, though a simplified formula was established in a theorem by Fuller [Proc. Natl. Acad. Sci. U.S.A. 75, 3557 (1978)]. We examine the range of applicability of this widely used simplified formula, and show that it cannot be employed for plectonemic DNA. We show that inapplicability is due to a hypothesis of Fuller theorem that is not met. The hypothesis seems to have been overlooked in many works.

Quantification of the shape of fiber tracts

The fiber tracts generated using diffusion MRI are usually simply displayed and assessed visually for a specific clinical or medical research purpose. This paper proposes computational techniques that can be used to study the shape of the tracts and make interindividual comparisons. These methods make use of fundamental geometric invariants, such as curvatures and torsions, or Fourier descriptors, together with the link of a pair of curves. Intersubject comparisons only require that the starting and ending points of the tracts can be defined and do not require point-by-point correspondences such as obtained using image registration. Principal component analysis-based shape analysis is also investigated. The invariants are tested on simulations and in vivo datasets, and the scale dependence and noise sensitivity of the measures are assessed. The potential for these techniques to be used in neuroscience research and clinical applications is demonstrated.

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.

Elastic correlations in nucleosomal DNA structure

The structure of DNA in the nucleosome core particle is studied using an elastic model that incorporates anisotropy in the bending energetics and twist-bend coupling. Using the experimentally determined structure of nucleosomal DNA [T. J. Richmond and C. A. Davey, Nature (London) 423, 145 (2003)], it is shown that elastic correlations exist between twist, roll, tilt, and stretching of DNA, as well as the distance between phosphate groups. The twist-bend coupling term is shown to be able to capture these correlations to a large extent, and a fit to the experimental data yields a new estimate of G = 25 nm for the value of the twist-bend coupling constant.

Elasticity and stability of a helical filament

We derive the general shape equations in terms of Euler angles for a uniform elastic rod with spontaneous torsion and curvatures and subjected to external force and torque. Our results based on an analytic formalism show that the extension of a helical rod may undergo a one-step discontinuous transition with increasing stretching force. This agrees quantitatively with experimental observations for a helix in a chemically defined lipid concentrate. The larger the twisting rigidity, the larger the jump in the extension. The effect of torque on the jump is, however, dependent on the value of the spontaneous torsion. In contrast, increasing the spontaneous torsion encourages the continuous variation of the extension. An "over-collapse" behavior is observed for the rod with asymmetric bending rigidity, and an intrinsic asymmetric elasticity under twisting force is found.

Nucleosome repositioning via loop formation

Active (catalyzed) and passive (intrinsic) nucleosome repositioning is known to be a crucial event during the transcriptional activation of certain eukaryotic genes. Here we consider theoretically the intrinsic mechanism and study in detail the energetics and dynamics of DNA-loop-mediated nucleosome repositioning, as previously proposed by earlier works. The surprising outcome of the present study is the inherent nonlocality of nucleosome motion within this model-being a direct physical consequence of the loop mechanism. On long enough DNA templates the longer jumps dominate over the previously predicted local motion, a fact that contrasts simple diffusive mechanisms considered before. The possible experimental outcome resulting from the considered mechanism is predicted, discussed, and compared to existing experimental findings.

Single, double, and multiple double strand breaks induced in DNA by 3-100 eV electrons

Nonthermal secondary electrons with initial kinetic energies below 100 eV are an abundant transient species created in irradiated cells and thermalize within picoseconds through successive multiple energy loss events. Here we show that below 15 eV such low-energy electrons induce single (SSB) and double (DSB) strand breaks in plasmid DNA exclusively via formation and decay of molecular resonances involving DNA components (base, sugar, hydration water, etc.). Furthermore, the strand break quantum yields (per incident electron) due to resonances occur with intensities similar to those that appear between 25 and 100 eV electron energy, where nonresonant mechanisms related to excitation/ionizations/dissociations are shown to dominate the yields, although with some contribution from multiple scattering electron energy loss events. We also present the first measurements of the electron energy dependence of multiple double strand breaks (MDSB) induced in DNA by electrons with energies below 100 eV. Unlike the SSB and DSB yields, which remain relatively constant above 25 eV, the MDSB yields show a strong monotonic increase above 30 eV, however with intensities at least 1 order of magnitude smaller than the combined SSB and DSB yields. The observation of MDSB above 30 eV is attributed to strand break clusters (nano-tracks) involving multiple successive interactions of one single electron at sites that are distant in primary sequence along the DNA double strand, but are in close contact; such regions exist in supercoiled DNA (as well as cellular DNA) where the double helix crosses itself or is in close proximity to another part of the same DNA molecule.

Automatic classification of protein structure by using Gauss integrals

We introduce a method of looking at, analyzing, and comparing protein structures. The topology of a protein is captured by 30 numbers inspired by Vassiliev knot invariants. To illustrate the simplicity and power of this topological approach, we construct a measure (scaled Gauss metric, SGM) of similarity of protein shapes. Under this metric, protein chains naturally separate into fold clusters. We use SGM to construct an automatic classification procedure for the CATH2.4 database. The method is very fast because it requires neither alignment of the chains nor any chain-chain comparison. It also has only one adjustable parameter. We assign 95.51% of the chains into the proper C (class), A (architecture), T (topology), and H (homologous superfamily) fold, find all new folds, and detect no false geometric positives. Using the SGM, we display a "map" of the space of folds projected onto two dimensions, show the relative locations of the major structural classes, and "zoom into" the space of proteins to show architecture, topology, and fold clusters. The existence of a simple measure of a protein fold computed from the chain path will have a major impact on automatic fold classification.

Monte Carlo simulations of supercoiled DNAs confined to a plane

Recent advances in atomic force microscopy (AFM) have enabled researchers to obtain images of supercoiled DNAs deposited on mica surfaces in buffered aqueous milieux. Confining a supercoiled DNA to a plane greatly restricts its configurational freedom, and could conceivably alter certain structural properties, such as its twist and writhe. A program that was originally written to perform Monte Carlo simulations of supercoiled DNAs in solution was modified to include a surface potential. This potential flattens the DNAs to simulate the effect of deposition on a surface. We have simulated transfers of a 3760-basepair supercoiled DNA from solution to a surface in both 161 and 10 mM ionic strength. In both cases, the geometric and thermodynamic properties of the supercoiled DNAs on the surface differ significantly from the corresponding quantities in solution. At 161 mM ionic strength, the writhe/twist ratio is 1.20-1.33 times larger for DNAs on the surface than for DNAs in solution and significant differences in the radii of gyration are also observed. Simulated surface structures in 161 mM ionic strength closely resemble those observed by AFM. Simulated surface structures in 10 mM ionic strength are similar to a minority of the structures observed by AFM, but differ from the majority of such structures for unknown reasons. In 161 mM ionic strength, the internal energy (excluding the surface potential) decreases substantially as the DNA is confined to the surface. Evidently, supercoiled DNAs in solution are typically deformed farther from the minimum energy configuration than are the corresponding surface-confined DNAs. Nevertheless, the work (Delta A(int)) done on the internal coordinates, which include uniform rotations at constant configuration, during the transfer is positive and 2.6-fold larger than the decrease in internal energy. The corresponding entropy change is negative, and its contribution to Delta A(int) is positive and exceeds the decrease in internal energy by 3.6 fold. The work done on the internal coordinates during the solution-to-surface transfer is directed primarily toward reducing their entropy. Evidently, the number of configurations available to the more deformed solution DNA is vastly greater than for the less deformed surface-confined DNA.

Constraints, histones, and the 30-nm spiral

We investigate the mechanical stability of a segment of DNA wrapped around a histone in the nucleosome configuration, under the assumption that the proper model for this packaging arrangement is that of an elastic rod that is free to twist and that writhes subject to mechanical constraints. We find that the number of constraints required to stabilize the nuclesome configuration is determined by the length of the segment, the number of times the DNA wraps around the histone spool, and the specific constraints utilized. While it can be shown that four constraints suffice, in principle, to insure stability of the nucleosome, a proper choice must be made to guarantee the effectiveness of this minimal number. The optimal choice of constraints appears to bear a relation to the existence of a spiral ridge on the surface of the histone octamer. The particular configuration that we investigate is related to the 30-nm spiral, a higher-order organization of DNA in chromatin.