Direct simulation of magnetic resonance relaxation rates and line shapes from molecular trajectories.

We simulate spin relaxation processes, which may be measured by either continuous wave or pulsed magnetic resonance techniques, using trajectory-based simulation methodologies. The spin-lattice relaxation rates are extracted numerically from the relaxation simulations. The rates obtained from the numerical fitting of the relaxation curves are compared to those obtained by direct simulation from the relaxation Bloch-Wangsness-Abragam-Redfield theory (BWART).

Estimation of the persistence length of DNA from the torsion elastic constant and supercoiling free energy: effect of ethylene glycol.

Two different methods are proposed to estimate the persistence length ( P) of DNA from the measured torsion elastic constant (alpha) and the twist energy parameter ( E T ) that governs the supercoiling free energy. The first method involves Monte Carlo simulations and reversible-work calculations of E T for model DNAs that possess the measured alpha and selected trial values of P. Comparison of the computed E T values with the experimental value allows estimation of P (or equivalently the bending elastic constant (kappa beta)) by interpolation.

Theory for spin-lattice relaxation of spin probes on weakly deformable DNA.

The weakly bending rod (WBR) model of double-stranded DNA (dsDNA) is adapted to analyze the internal dynamics of dsDNA as observed in electron paramagnetic resonance (EPR) measurements of the spin-lattice relaxation rate, R(1e), for spin probes rigidly attached to nucleic acid-bases. The WBR theory developed in this work models dsDNA base-pairs as diffusing rigid cylindrical discs connected by bending and twisting springs whose elastic force constants are kappa and alpha, respectively. Angular correlation functions for both rotational displacement and velocity are developed in detail so as to compute values for R(1e) due to four relaxation mechanisms: the chemical shift anisotropy (CSA), the electron-nuclear dipolar (END), the spin rotation (SR), and the generalized spin diffusion (GSD) relaxation processes.

Simulation of nitroxide electron paramagnetic resonance spectra from brownian trajectories and molecular dynamics simulations.

A simulated continuous wave electron paramagnetic resonance spectrum of a nitroxide spin label can be obtained from the Fourier transform of a free induction decay. It has been previously shown that the free induction decay can be calculated by solving the time-dependent stochastic Liouville equation for a set of Brownian trajectories defining the rotational dynamics of the label. In this work, a quaternion-based Monte Carlo algorithm has been developed to generate Brownian trajectories describing the global rotational diffusion of a spin-labeled protein.

Effects of ethylene glycol on the torsion elastic constant and hydrodynamic radius of p30delta DNA.

Upon increasing the concentration of ethylene glycol (EG) at 37 degrees C, the twist energy parameter, E(T), which governs the supercoiling free energy, was recently found to undergo a decreasing (or reverse) sigmoidal transition with a midpoint near 20 w/v % EG. In this study, the effects of adding 20 w/v % EG on the torsion elastic constant (alpha) of linear p30delta DNA and on the hydrodynamic radius (R(H)) of a synthetic 24 bp duplex DNA were examined at both 40 and 20 degrees C. The time-resolved fluorescence intensity and fluorescence polarization anisotropy (FPA) of intercalated ethidium were measured in order to assess the effects of 20 w/v % EG on: (1) alpha; (2) R(H); (3) the lifetimes of intercalated and non-intercalated dye; (4) the amplitude of dye wobble in its binding site; and (5) the binding constant for intercalation.

A contribution to the theory of preferential interaction coefficients.

A simple and complete derivation of the relation between concentration-based preferential interaction coefficients and integrals over the relevant pair correlation functions is presented for the first time. Certain omissions from the original treatment of pair correlation functions in multicomponent thermodynamics are also addressed. Connections between these concentration-based quantities and the more common molality-based preferential interaction coefficients are also derived.

Effects of small neutral osmolytes on the supercoiling free energy and intrinsic twist of p30delta DNA.

Both theory and experiments are employed to investigate the effects of small neutral osmolytes on the average intrinsic twist (l0), the torsion and bending elastic constants, and the twist energy parameter (ET) that governs the supercoiling free energy. The experimental data for ethylene glycol and acetamide at 37 degrees C suggest, and are interpreted in terms of, a model wherein the DNA exhibits an equilibrium between two distinct conformational states that possess different numbers of bound water molecules and exhibit different intrinsic twists and torsion and bending elastic constants. Expressions are derived to relate the effective ET and l0 to the equilibrium constant, water activity (aw), and number (n) of bound water molecules released per cooperative domain undergoing the two-state transition.

Monte Carlo simulations of locally melted supercoiled DNAs in 20 mM ionic strength.

Mesoscopic models of unmelted and locally melted supercoiled DNAs in 20 mM ionic strength are simulated over a range of linking difference from deltal = 0 to -26 turns, or superhelix density from sigma = 0 to -0.062. A domain containing m = 0, 28, or 56 melted basepairs (out of 4349 total) is modeled simply by a region of suitable length with substantially reduced torsion and bending elastic constants.