Mesoscopic Hamiltonian Research Articles

End-to-end distance and contour length distribution functions of DNA helices.

I present a computational method to evaluate the end-to-end and the contour length distribution functions of short DNA molecules described by a mesoscopic Hamiltonian. The method generates a large statistical ensemble of possible configurations for each dimer in the sequence, selects the global equilibrium twist conformation for the molecule, and determines the average base pair distances along the molecule backbone. Integrating over the base pair radial and angular fluctuations, I derive the room temperature distribution functions as a function of the sequence length. The obtained values for the most probable end-to-end distance and contour length distance, providing a measure of the global molecule size, are used to examine the DNA flexibility at short length scales. It is found that, also in molecules with less than ∼60 base pairs, coiled configurations maintain a large statistical weight and, consistently, the persistence lengths may be much smaller than in kilo-base DNA.

Twisting short dsDNA with applied tension

The twisting deformation of mechanically stretched DNA molecules is studied by a coarse grained Hamiltonian model incorporating the fundamental interactions that stabilize the double helix and accounting for the radial and angular base pair fluctuations. The latter are all the more important at short length scales in which DNA fragments maintain an intrinsic flexibility. The presented computational method simulates a broad ensemble of possible molecule conformations characterized by a specific average twist and determines the energetically most convenient helical twist by free energy minimization. As this is done for any external load, the method yields the characteristic twist-stretch profile of the molecule and also computes the changes in the macroscopic helix parameters i.e. average diameter and rise distance. It is predicted that short molecules under stretching should first over-twist and then untwist by increasing the external load. Moreover, applying a constant load and simulating a torsional strain which over-twists the helix, it is found that the average helix diameter shrinks while the molecule elongates, in agreement with the experimental trend observed in kilo-base long sequences. The quantitative relation between percent relative elongation and superhelical density at fixed load is derived. The proposed theoretical model and computational method offer a general approach to characterize specific DNA fragments and predict their macroscopic elastic response as a function of the effective potential parameters of the mesoscopic Hamiltonian.

Twist-stretch profiles of DNA chains

Helical molecules change their twist number under the effect of a mechanical load. We study the twist-stretch relation for a set of short DNA molecules modeled by a mesoscopic Hamiltonian. Finite temperature path integral techniques are applied to generate a large ensemble of possible configurations for the base pairs of the sequence. The model also accounts for the bending and twisting fluctuations between adjacent base pairs along the molecules stack. Simulating a broad range of twisting conformation, we compute the helix structural parameters by averaging over the ensemble of base pairs configurations. The method selects, for any applied force, the average twist angle which minimizes the molecule’s free energy. It is found that the chains generally over-twist under an applied stretching and the over-twisting is physically associated to the contraction of the average helix diameter, i.e. to the damping of the base pair fluctuations. Instead, assuming that the maximum amplitude of the bending fluctuations may decrease against the external load, the DNA molecule first over-twists for weak applied forces and then untwists above a characteristic force value. Our results are discussed in relation to available experimental information albeit for kilo-base long molecules.

Mesoscopic Hamiltonian for the fluctuations of adsorbed Lennard-Jones liquid films.

We use Monte Carlo simulations of a Lennard-Jones fluid adsorbed on a short-range planar wall substrate to study the fluctuations in the thickness of the wetting layer, and we get a quantitative and consistent characterization of their mesoscopic Hamiltonian, H[ξ]. We have observed important finite-size effects, which were hampering the analysis of previous results obtained with smaller systems. The results presented here support an appealing simple functional form for H[ξ], close but not exactly equal to the theoretical nonlocal proposal made on the basis a generic density-functional analysis by Parry and coworkers. We have analyzed systems under different wetting conditions, as a proof of principle for a method that provides a quantitative bridge between the molecular interactions and the phenomenology of wetting films at mesoscopic scales.

Twisting and bending stress in DNA minicircles.

The interplay between bending of the molecule axis and appearance of disruptions in circular DNA molecules, with ∼100 base pairs, is addressed. Three minicircles with different radii and almost equal contents of AT and GC pairs are investigated. The DNA sequences are modeled by a mesoscopic Hamiltonian which describes the essential interactions in the helix at the level of the base pair and incorporates twisting and bending degrees of freedom. Helix unwinding and bubble formation patterns are consistently computed by a path integral method that sums over a large number of molecule configurations compatible with the model potential. The path ensembles are determined, as a function of temperature, by minimizing the free energy of the system. Fluctuational openings appear along the helix to release the stress due to the bending of the molecule backbone. In agreement with the experimental findings, base pair disruptions are found with larger probability in the smallest minicircle of 66 bps whose bending angle is ∼6°. For this minicircle, a sizeable untwisting is obtained with the helical repeat showing a step-like increase at T = 315 K. The method can be generalized to determine the bubble probability profiles of open ends linear sequences.

Modeling DNA Dynamics by Path Integrals

Complementary strands in DNA double helix show temporary fluctuational openings which are essential to biological functions such as transcription and replication of the genetic information. Such large amplitude fluctuations, known as the breathing of DNA, are generally localized and, microscopically, are due to the breaking of the hydrogen bonds linking the base pairs (bps). I apply imaginary time path integral techniques to a mesoscopic Hamiltonian which accounts for the helicoidal geometry of a short circular DNA molecule. The bps displacements with respect to the ground state are interpreted as time dependent paths whose amplitudes are consistent with the model potential for the hydrogen bonds. The portion of the paths configuration space contributing to the partition function is determined by selecting the ensemble of paths which fulfill the second law of thermodynamics. Computations of the thermodynamics in the denaturation range show the energetic advantage for the equilibrium helicoidal geometry peculiar of B-DNA. I discuss the interplay between twisting of the double helix and anharmonic stacking along the molecule backbone suggesting an interesting relation between intrinsic nonlinear character of the microscopic interactions and molecular topology.