Conformational flexibility of DNA: an extension of the stereochemical guidelines

A consistent potential energy parameter set for lipids: dipalmitoylphosphatidylcholine as a benchmark of the GROMOS96 45A3 force field

The performance of the GROMOS96 parameter set 45A3 developed for aliphatic alkanes is tested on a bilayer of dipalmitoylphosphatidylcholine (DPPC) in water in the liquid-crystalline L(alpha) phase. Variants of the force-field parameter set as well as different sets of simulation conditions or simulation parameter sets are evaluated. In the case of the force-field parameters, the van der Waals constants for the non-bonded interaction of the ester carbonyl carbon and the partial charges and charge group definition of the phosphatidylcholine head group are examined. On the methodological side, different cut-off distances for the non-bonded interactions, use of a reaction-field force due to long-range electrostatic interactions, the frequency of removal of the centre of mass motion and the strength of the coupling of the pressure of the system to the pressure bath are tested. The area per lipid, as a measure of structure, the order parameters of the chain carbons, as a measure of membrane fluidity, and the translational diffusion of the lipids in the plane of the bilayer are calculated and compared with experimental values. An optimal set of simulation parameters for which the GROMOS96 parameter set 45A3 yields a head group area, chain order parameters and a lateral diffusion coefficient in accordance with the experimental data is listed.

A general procedure for generation of curved DNA molecules

A general method for generation of base-pairs in a curved DNA structure, for any prescribed values of helical parameters--unit rise (h), unit twist (theta), wedge roll (theta R) and wedge tilt (theta T), propeller twist (theta p) and displacement (D) is described. Its application for generation of uniform as well curved structures is also illustrated with some representative examples. An interesting relationship is observed between helical twist (theta), base-pair parameters theta x, theta y and the wedge parameters theta R, theta T, which has important consequences for the description and estimation of DNA curvature.

Energetics of left and right handed models of DNA

It has been shown by model building studies that various right handed and left handed models are compatible with X-ray data of B-DNA and C-DNA. These models are also found to be in good agreement with infrared dichroism data. Detailed potential energy calculations have now been carried out for these models, viz., right and left handed B-DNA and right and left handed C-DNA. It is found that base sugar stacking and interactions involving the phosphate groups are the dominant forces for stabilizing a particular structure. For some sequences, viz., A-A, T-A and C-A, left handed stacking is quite favourable in both B and C structures. But intranucleotide interactions make the left B-DNA unfavourable while the left C-DNA structure is more stable, for all the sequences, than the right C-DNA structure, proposed from fibre data. For the hexanucleoside pentaphosphate fragments the same trend is observed, with the right handed B-DNA being the most stable of the four models studied. However, the left C-DNA structure is only marginally higher in energy, particularly if the shielding effect of the counter ions, on the phosphate group is taken into consideration.