Heat conductivity of DNA double helix

Thermal conductivity of isolated single molecule DNA fragments is of importance for nanotechnology, but has not yet been measured experimentally. Theoretical estimates based on simplified (1D) models predict anomalously high thermal conductivity. To investigate thermal properties of single molecule DNA we have developed a 3D coarse-grained (CG) model that retains the realism of the full all-atom description, but is significantly more efficient. Within the proposed model each nucleotide is represented by 6 particles or grains; the grains interact via effective potentials inferred from classical molecular dynamics (MD) trajectories based on a well-established all-atom potential function. Comparisons of 10 ns long MD trajectories between the CG and the corresponding all-atom model show similar root-mean-square deviations from the canonical B-form DNA, and similar structural fluctuations. At the same time, the CG model is 10 to 100 times faster depending on the length of the DNA fragment in the simulation. Analysis of dispersion curves derived from the CG model yields longitudinal sound velocity and torsional stiffness in close agreement with existing experiments. The computational efficiency of the CG model makes it possible to calculate thermal conductivity of a single DNA molecule not yet available experimentally. For a uniform (polyG-polyC) DNA, the estimated conductivity coefficient is 0.3 W/mK which is half the value of thermal conductivity for water. This result is in stark contrast with estimates of thermal conductivity for simplified, effectively 1D chains ("beads on a spring") that predict anomalous (infinite) thermal conductivity. Thus, full 3D character of DNA double-helix retained in the proposed model appears to be essential for describing its thermal properties at a single molecule level.Comment: 16 pages, 12 figure

Understanding Nucleic Acid Structural Changes by Comparing Wide-Angle X-ray Scattering (WAXS) Experiments to Molecular Dynamics Simulations

Wide-angle x-ray scattering (WAXS) is emerging as a powerful tool for increasing the resolution of solution structure measurements of biomolecules. Compared to its better known complement, small angle x-ray scattering (SAXS), WAXS targets higher scattering angles and can enhance structural studies of molecules by accessing finer details of solution structures. Although the extension from SAXS to WAXS is easy to implement experimentally, the computational tools required to fully harness the power of WAXS are still under development. Currently, WAXS is employed to study structural changes and ligand binding in proteins; however the methods are not as fully developed for nucleic acids. Here, we show how WAXS can qualitatively characterize nucleic acid structures as well as the small but significant structural changes driven by the addition of multivalent ions. We show the potential of WAXS to test all-atom molecular dynamics (MD) simulations and to provide insight in understanding how the trivalent ion cobalt(III) hexammine (CoHex) affects the structure of RNA and DNA helices. We find that MD simulations capture the RNA structural change that occurs due to addition of CoHex

SU(N) quantum spin models: A variational wavefunction study

The study of SU(N) quantum spin models is relevant to a variety of physical systems including ultracold atoms in optical lattices, and also leads to insights into novel quantum phases and phase transitions of SU(2) spin models. We use Gutzwiller projected fermionic variational wavefunctions to explore the phase diagram and correlation functions of SU(N) spin models in the self-conjugate representation, with Heisenberg bilinear and biquadratic interactions. In 1D, the variational phase diagram of the SU(4) spin chain is constructed by examining instabilities of the Gutzwiller projected free fermion ground state to various broken symmetries, and it agrees well with exact results.The spin and dimer correlations of the Gutzwiller projected free fermion state with N flavors of fermions are also in good agreement with exact and 1/N calculations for the critical points of SU(N) spin chains. In 2D, the variational phase diagram on the square lattice is obtained by studying instabilities of the Gutzwiller projected pi-flux state. The variational ground state of the pure Heisenberg model is found to exhibit long range Neel order for N=2,4 and spin Peierls order for N > 4. For N=4 and 6, biquadratic interactions lead to a complex phase diagram which includes an extended valence bond crystal in both cases, as well as a stable pi-flux phase for N=6. The spin correlations of the projected pi-flux state at N=4 are in good agreement with 1/N calculations. We find that this state also shows strongly enhanced dimer correlations, in qualitative accord with the large-N results. We compare our results with a recent QMC study of the SU(4) Heisenberg model.Comment: 22 pages, 7 figs, added references to arxiv versio

Two-phase stretching of molecular chains

While stretching of most polymer chains leads to rather featureless force-extension diagrams, some, notably DNA, exhibit non-trivial behavior with a distinct plateau region. Here we propose a unified theory that connects force-extension characteristics of the polymer chain with the convexity properties of the extension energy profile of its individual monomer subunits. Namely, if the effective monomer deformation energy as a function of its extension has a non-convex (concave up) region, the stretched polymer chain separates into two phases: the weakly and strongly stretched monomers. Simplified planar and 3D polymer models are used to illustrate the basic principles of the proposed model. Specifically, we show rigorously that when the secondary structure of a polymer is mostly due to weak non-covalent interactions, the stretching is two-phase, and the force-stretching diagram has the characteristic plateau. We then use realistic coarse-grained models to confirm the main findings and make direct connection to the microscopic structure of the monomers. We demostrate in detail how the two-phase scenario is realized in the \alpha-helix, and in DNA double helix. The predicted plateau parameters are consistent with single molecules experiments. Detailed analysis of DNA stretching demonstrates that breaking of Watson-Crick bonds is not necessary for the existence of the plateau, although some of the bonds do break as the double-helix extends at room temperature. The main strengths of the proposed theory are its generality and direct microscopic connection.Comment: 16 pges, 22 figure

Predicting Transcription Factor Specificity with All-Atom Models

The binding of a transcription factor (TF) to a DNA operator site can initiate or repress the expression of a gene. Computational prediction of sites recognized by a TF has traditionally relied upon knowledge of several cognate sites, rather than an ab initio approach. Here, we examine the possibility of using structure-based energy calculations that require no knowledge of bound sites but rather start with the structure of a protein-DNA complex. We study the PurR E. coli TF, and explore to which extent atomistic models of protein-DNA complexes can be used to distinguish between cognate and non-cognate DNA sites. Particular emphasis is placed on systematic evaluation of this approach by comparing its performance with bioinformatic methods, by testing it against random decoys and sites of homologous TFs. We also examine a set of experimental mutations in both DNA and the protein. Using our explicit estimates of energy, we show that the specificity for PurR is dominated by direct protein-DNA interactions, and weakly influenced by bending of DNA.Comment: 26 pages, 3 figure

Phase diagram and symmetry breaking of SU(4) spin-orbital chain in a generalized external field

The ground state phases of a one-dimensional SU(4) spin-orbital Hamiltonian in a generalized external field are studied on the basis of Bethe-ansatz solution. Introducing three Land\'e $g$ factors for spin, orbital and their products in the SU(4) Zeeman term, we discuss systematically the various symmetry breaking. The magnetization versus external field are obtained by solving Bethe-ansatz equations numerically. The phase diagrams corresponding to distinct residual symmetries are given by means of both numerical and analytical methods.Comment: Revtex4, 16 pages, 7 figure

Integrity of H1 helix in prion protein revealed by molecular dynamic simulations to be especially vulnerable to changes in the relative orientation of H1 and its S1 flank

In the template-assistance model, normal prion protein (PrPC), the pathogenic cause of prion diseases such as Creutzfeldt-Jakob (CJD) in human, Bovine Spongiform Encephalopathy (BSE) in cow, and scrapie in sheep, converts to infectious prion (PrPSc) through an autocatalytic process triggered by a transient interaction between PrPC and PrPSc. Conventional studies suggest the S1-H1-S2 region in PrPC to be the template of S1-S2 $\beta$-sheet in PrPSc, and the conformational conversion of PrPC into PrPSc may involve an unfolding of H1 in PrPC and its refolding into the $\beta$-sheet in PrPSc. Here we conduct a series of simulation experiments to test the idea of transient interaction of the template-assistance model. We find that the integrity of H1 in PrPC is vulnerable to a transient interaction that alters the native dihedral angles at residue Asn$^{143}$, which connects the S1 flank to H1, but not to interactions that alter the internal structure of the S1 flank, nor to those that alter the relative orientation between H1 and the S2 flank.Comment: A major revision on statistical analysis method has been made. The paper now has 23 pages, 11 figures. This work was presented at 2006 APS March meeting session K29.0004 at Baltimore, MD, USA 3/13-17, 2006. This paper has been accepted for pubcliation in European Biophysical Journal on Feb 2, 200