Bubble statistics and positioning in superhelically stressed DNA

We present a general framework to study the thermodynamic denaturation of double-stranded DNA under superhelical stress. We report calculations of position- and size-dependent opening probabilities for bubbles along the sequence. Our results are obtained from transfer-matrix solutions of the Zimm-Bragg model for unconstrained DNA and of a self-consistent linearization of the Benham model for superhelical DNA. The numerical efficiency of our method allows for the analysis of entire genomes and of random sequences of corresponding length ($10^6-10^9$ base pairs). We show that, at physiological conditions, opening in superhelical DNA is strongly cooperative with average bubble sizes of $10^2-10^3$ base pairs (bp), and orders of magnitude higher than in unconstrained DNA. In heterogeneous sequences, the average degree of base-pair opening is self-averaging, while bubble localization and statistics are dominated by sequence disorder. Compared to random sequences with identical GC-content, genomic DNA has a significantly increased probability to open large bubbles under superhelical stress. These bubbles are frequently located directly upstream of transcription start sites.Comment: to be appeared in Physical Review

Thermodynamics of Twisted DNA with Solvent Interaction

The imaginary time path integral formalism is applied to a nonlinear Hamiltonian for a short fragment of heterogeneous DNA with a stabilizing solvent interaction term. Torsional effects are modeled by a twist angle between neighboring base pairs stacked along the molecule backbone. The base pair displacements are described by an ensemble of temperature dependent paths thus incorporating those fluctuational effects which shape the multisteps thermal denaturation. By summing over $\sim 10^7 - 10^8$ base pair paths, a large number of double helix configurations is taken into account consistently with the physical requirements of the model potential. The partition function is computed as a function of the twist. It is found that the equilibrium twist angle, peculiar of B-DNA at room temperature, yields the stablest helicoidal geometry against thermal disruption of the base pair hydrogen bonds. This result is corroborated by the computation of thermodynamical properties such as fractions of open base pairs and specific heat.Comment: The Journal of Chemical Physics (2011) in pres

Bubble generation in a twisted and bent DNA-like model

The DNA molecule is modeled by a parabola embedded chain with long-range interactions between twisted base pair dipoles. A mechanism for bubble generation is presented and investigated in two different configurations. Using random normally distributed initial conditions to simulate thermal fluctuations, a relationship between bubble generation, twist and curvature is established. An analytical approach supports the numerical results.Comment: 7 pages, 8 figures. Accepted for Phys. Rev. E (in press

Fluctuation-Facilitated Charge Migration along DNA

We propose a model Hamiltonian for charge transfer along the DNA double helix with temperature driven fluctuations in the base pair positions acting as the rate limiting factor for charge transfer between neighboring base pairs. We compare the predictions of the model with the recent work of J.K. Barton and A.H. Zewail (Proc.Natl.Acad.Sci.USA, {\bf 96}, 6014 (1999)) on the unusual two-stage charge transfer of DNA.Comment: 4 pages, 2 figure

Statics and Dynamics of the Wormlike Bundle Model

Bundles of filamentous polymers are primary structural components of a broad range of cytoskeletal structures, and their mechanical properties play key roles in cellular functions ranging from locomotion to mechanotransduction and fertilization. We give a detailed derivation of a wormlike bundle model as a generic description for the statics and dynamics of polymer bundles consisting of semiflexible polymers interconnected by crosslinking agents. The elastic degrees of freedom include bending as well as twist deformations of the filaments and shear deformation of the crosslinks. We show that a competition between the elastic properties of the filaments and those of the crosslinks leads to renormalized effective bend and twist rigidities that become mode-number dependent. The strength and character of this dependence is found to vary with bundle architecture, such as the arrangement of filaments in the cross section and pretwist. We discuss two paradigmatic cases of bundle architecture, a uniform arrangement of filaments as found in F-actin bundles and a shell-like architecture as characteristic for microtubules. Each architecture is found to have its own universal ratio of maximal to minimal bending rigidity, independent of the specific type of crosslink induced filament coupling; our predictions are in reasonable agreement with available experimental data for microtubules. Moreover, we analyze the predictions of the wormlike bundle model for experimental observables such as the tangent-tangent correlation function and dynamic response and correlation functions. Finally, we analyze the effect of pretwist (helicity) on the mechanical properties of bundles. We predict that microtubules with different number of protofilaments should have distinct variations in their effective bending rigidity

Cymantrene–Triazole "Click" Products: Structural Characterization and Electrochemical Properties

We report the first known examples of triazole-derivatized cymantrene complexes (η5-[4-substituted triazol-1-yl]cyclopentadienyl)tricarbonylmanganese(I), obtained via a “click” chemical synthesis, bearing a phenyl, 3-aminophenyl, or 4-aminophenyl moiety at the 4-position of the triazole ring. Structural characterization data using multinuclear NMR, UV–vis, ATR-IR, and mass spectrometric methods are provided, as well as crystallographic data for (η5-[4-phenyltriazol-1-yl]cyclopentadienyl)tricarbonylmanganese(I) and (η5-[4-(3-aminophenyl)triazol-1-yl]cyclopentadienyl)tricarbonylmanganese(I). Cyclic voltammetric characterization of the redox behavior of each of the three cymantrene–triazole complexes is presented together with digital simulations, in situ infrared spectroelectrochemistry, and DFT calculations to extract the associated kinetic and thermodynamic parameters. The trypanocidal activity of each cymantrene–triazole complex is also examined, and these complexes are found to be more active than cymantrene alone

Two-level system with a thermally fluctuating transfer matrix element: Application to the problem of DNA charge transfer

Charge transfer along the base-pair stack in DNA is modeled in terms of thermally-assisted tunneling between adjacent base pairs. Central to our approach is the notion that tunneling between fluctuating pairs is rate-limited by the requirement of their optimal alignment. We focus on this aspect of the process by modeling two adjacent base pairs in terms of a classical damped oscillator subject to thermal fluctuations as described by a Fokker-Planck equation. We find that the process is characterized by two time scales, a result that is in accord with experimental findings.Comment: original file is revtex4, 10 pages, three eps figure

Calculating singlet excited states: comparison with fast time-resolved infrared spectroscopy of coumarins

In contrast to the ground state, the calculation of the infrared (IR) spectroscopy of molecular singlet excited states represents a substantial challenge. Here we use the structural IR fingerprint of the singlet excited states of a range of coumarin dyes to assess the accuracy of density functional theory based methods for the calculation of excited state IR spectroscopy. It is shown that excited state Kohn-Sham density functional theory provides a high level of accuracy and represents an alternative approach to time-dependent density functional theory for simulating the IR spectroscopy of the singlet excited states

A stitch in time: Efficient computation of genomic DNA melting bubbles

Background: It is of biological interest to make genome-wide predictions of the locations of DNA melting bubbles using statistical mechanics models. Computationally, this poses the challenge that a generic search through all combinations of bubble starts and ends is quadratic. Results: An efficient algorithm is described, which shows that the time complexity of the task is O(NlogN) rather than quadratic. The algorithm exploits that bubble lengths may be limited, but without a prior assumption of a maximal bubble length. No approximations, such as windowing, have been introduced to reduce the time complexity. More than just finding the bubbles, the algorithm produces a stitch profile, which is a probabilistic graphical model of bubbles and helical regions. The algorithm applies a probability peak finding method based on a hierarchical analysis of the energy barriers in the Poland-Scheraga model. Conclusions: Exact and fast computation of genomic stitch profiles is thus feasible. Sequences of several megabases have been computed, only limited by computer memory. Possible applications are the genome-wide comparisons of bubbles with promotors, TSS, viral integration sites, and other melting-related regions.Comment: 16 pages, 10 figure

Statistical Mechanics of Torque Induced Denaturation of DNA

A unifying theory of the denaturation transition of DNA, driven by temperature T or induced by an external mechanical torque Gamma is presented. Our model couples the hydrogen-bond opening and the untwisting of the helicoidal molecular structure. We show that denaturation corresponds to a first-order phase transition from B-DNA to d-DNA phases and that the coexistence region is naturally parametrized by the degree of supercoiling sigma. The denaturation free energy, the temperature dependence of the twist angle, the phase diagram in the T,Gamma plane and isotherms in the sigma, Gamma plane are calculated and show a good agreement with experimental data.Comment: 5 pages, 3 figures, model improve