The conformation and stability of ribonucleic acids: modeling base sequence effects in double stranded helices

Static and dynamic conformational properties of AT sequences in B-DNA

A theoretical study of the optimal conformations of nucleic acid oligomers containing tracts of AT base pairs is presented. The oligomers are studied in isolation and complexed with netropsin, a minor groove binding ligand. The flexibility of the oligomers and of their complexes is calculated by adiabatic mapping with respect to the total winding angle. The results of this study show that in uncomplexed oligomers the dinucleotide junctions AA, AT and TA have very different structural parameters and different responses to winding stress. The TA junction is clearly the most flexible and is the principal site for accommodating the imposed overwinding. Complexation by netropsin leads to two important effects: firstly, the three junctions adopt more uniform structures, the largest changes again being observed for TA, secondly, the differences in flexibility as a function of sequence are strongly attenuated.

Modeling of nucleic acid complexes with cationic ligands: a specialized molecular mechanics force field and its application

A potential energy force field designed for modeling nucleic acids and particularly their complexes with cationic ligands is presented. The force field is a modified version of that developed by Weiner, S.J., Kollman, P.A., Nguyen, D.T. and Case, D.A.,J. Comp. Chem. 7,230-252 (1986) and is based upon the use of a distance dependent dielectric constant, epsilon = 4rij, and partially neutralized phosphates to represent solvent and counterion. Changes from the Weiner et al. force field include additional atom types and modifications to van der Waals, electrostatic, hydrogen bonding and torsional parameters. Molecular modeling test cases of the force field are presented for a number of simple small molecules, as well as uracil and benzene dimerization, thymine-adenine and cytosine-guanine base pair formation, and adenosine/deoxyadenosine pseudorotation. Several DNA and RNA oligomers and DNA/RNA intercalation complexes with ethidium are also modeled with the force field. In all cases, the modeling results compare favorably with available experimental results. Additionally, conformational trends observed experimentally for nucleic acids by NMR and X-ray crystallographic techniques are reproduced. The modeling results for ethidium intercalation indicate a complex in which the favorable interactions are primarily van der Waals contacts, and in which electrostatic interactions are a relatively minor component. We feel the force field is particularly useful for molecular mechanics aided drug design, and an analysis of modeling results with respect to design of drugs which bind selectively to RNA is presented.