Richard Lavery

DNA: An Extensible Molecule

The force-displacement response of a single duplex DNA molecule was measured. The force saturates at a plateau around 70 piconewtons, which ends when the DNA has been stretched about 1.7 times its contour length. This behavior reveals a highly cooperative transition to a state here termed S-DNA. Addition of an intercalator suppresses this transition. Molecular modeling of the process also yields a force plateau and suggests a structure for the extended form. These results may shed light on biological processes involving DNA extension and open the route for mechanical studies on individual molecules in a previously unexplored range.

A New Approach to the Rapid Determination of Protein Side Chain Conformations

Two efficient algorithms have been developed which allow amino acid side chain conformations to be optimized rapidly for a given peptide backbone conformation. Both these approaches are based on the assumption that each side chain can be represented by a small number of rotameric states. These states have been obtained by a dynamic cluster analysis of a large data base of known crystallographic structures. Successful applications of these algorithms to the prediction of known protein conformations are presented.

A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA

It is well recognized that base sequence exerts a significant influence on the properties of DNA and plays a significant role in protein–DNA interactions vital for cellular processes. Understanding and predicting base sequence effects requires an extensive structural and dynamic dataset which is currently unavailable from experiment. A consortium of laboratories was consequently formed to obtain this information using molecular simulations. This article describes results providing information not only on all 10 unique base pair steps, but also on all possible nearest-neighbor effects on these steps. These results are derived from simulations of 50–100 ns on 39 different DNA oligomers in explicit solvent and using a physiological salt concentration. We demonstrate that the simulations are converged in terms of helical and backbone parameters. The results show that nearest-neighbor effects on base pair steps are very significant, implying that dinucleotide models are insufficient for predicting sequence-dependent behavior. Flanking base sequences can notably lead to base pair step parameters in dynamic equilibrium between two conformational sub-states. Although this study only provides limited data on next-nearest-neighbor effects, we suggest that such effects should be analyzed before attempting to predict the sequence-dependent behavior of DNA.

JUMNA (junction minimisation of nucleic acids)

The latest version of JUMNA program for the energy optimisation of nucleic acids and nucleic acid-ligand complexes is described. The force field used, input and output data, various options for symmetry, conformational constraints and energy mapping are discussed as well as recent combinatorial search techniques.

On the molecular mechanism of drug intercalation into DNA: a simulation study of the intercalation pathway, free energy, and DNA structural changes.

Intercalation into DNA (insertion between a pair of base pairs) is a critical step in the function of many anticancer drugs. Despite its importance, a detailed mechanistic understanding of this process at the molecular level is lacking. We have constructed, using extensive atomistic computer simulations and umbrella sampling techniques, a free energy landscape for the intercalation of the anticancer drug daunomycin into a twelve base pair B-DNA. A similar free energy landscape has been constructed for a probable intermediate DNA minor groove-bound state. These allow a molecular level understanding of aspects of the thermodynamics, DNA structural changes, and kinetic pathways of the intercalation process. Key DNA structural changes involve opening the future intercalation site base pairs toward the minor groove (positive roll), followed by an increase in the rise, accompanied by hydrogen bonding changes of the minor groove waters. The calculated intercalation free energy change is -12.3 kcal/mol, in reasonable agreement with the experimental estimate -9.4 kcal/mol. The results point to a mechanism in which the drug first binds to the minor groove and then intercalates into the DNA in an activated process, which is found to be in general agreement with experimental kinetic results.

Wringing Out DNA

The chiral nature of DNA plays a crucial role in cellular processes. Here we use magnetic tweezers to explore one of the signatures of this chirality, the coupling between stretch and twist deformations. We show that the extension of a stretched DNA molecule increases linearly by 0.42 nm per excess turn applied to the double helix. This result contradicts the intuition that DNA should lengthen as it is unwound and get shorter with overwinding. We then present numerical results of energy minimizations of torsionally restrained DNA that display a behavior similar to the experimental data and shed light on the molecular details of this surprising effect.

Unraveling Proteins: A Molecular Mechanics Study

An internal coordinate molecular mechanics study of unfolding peptide chains by external stretching has been carried out to predict the type of force spectra that may be expected from single-molecule manipulation experiments currently being prepared. Rather than modeling the stretching of a given protein, we have looked at the behavior of simple secondary structure elements (alpha-helix, beta-ribbon, and interacting alpha-helices) to estimate the magnitude of the forces involved in their unfolding or separation and the dependence of these forces on the way pulling is carried out as well as on the length of the structural elements. The results point to a hierarchy of forces covering a surprisingly large range and to important orientational effects in the response to external stress.

CURVES+ web server for analyzing and visualizing the helical, backbone and groove parameters of nucleic acid structures

Curves+, a revised version of the Curves software for analyzing the conformation of nucleic acid structures, is now available as a web server. This version, which can be freely accessed at http://gbio-pbil.ibcp.fr/cgi/Curves_plus/, allows the user to upload a nucleic acid structure file, choose the nucleotides to be analyzed and after optionally setting a number of input variables, view the numerical and graphic results online or download files containing a set of helical, backbone and groove parameters that fully describe the structure. PDB format files are also provided for offline visualization of the helical axis and groove geometry.

μABC: a systematic microsecond molecular dynamics study of tetranucleotide sequence effects in B-DNA.

We present the results of microsecond molecular dynamics simulations carried out by the ABC group of laboratories on a set of B-DNA oligomers containing the 136 distinct tetranucleotide base sequences. We demonstrate that the resulting trajectories have extensively sampled the conformational space accessible to B-DNA at room temperature. We confirm that base sequence effects depend strongly not only on the specific base pair step, but also on the specific base pairs that flank each step. Beyond sequence effects on average helical parameters and conformational fluctuations, we also identify tetranucleotide sequences that oscillate between several distinct conformational substates. By analyzing the conformation of the phosphodiester backbones, it is possible to understand for which sequences these substates will arise, and what impact they will have on specific helical parameters.

Structure and mechanics of single biomolecules: experiment and simulation

One of the main goals of molecular biology is to understand the structure of biomolecules. With the emergence of single molecule manipulation techniques that structure can now be controlled by the application of stretching and torsional stresses. In this article we review some recent experiments on the stretching and twisting of single biopolymers, testing the elastic properties of DNA and proteins and studying their stress-induced structural transitions. Numerical simulations have emerged as a precious tool to interpret the experimental data and predict the associated structural changes. We shall explain how a combination of these experimental and computational tools open a new vista on the structure of biomolecules.