Calculations by molecular mechanics methods showed that diglycerides containing stearic, oleic, linoleic, or linolenic acids can bind to DNA. The binding energy of diglycerides to the DNA minor groove is higher than that to the major groove. The bond energies of diglycerides to DNA were calculated as the sums of the contributions of the binding energies between DNA and two fatty acid residues. These energies depend on the structure of the fatty acid and the nucleotide composition of DNA and vary from 20 to 120 kcal mol -1 .In our previous studies, 1-3 we have used molecular mechanics methods to determine the structures and estimate the bond energies of DNA to saturated and unsaturated fatty acids containing 18 carbon atoms, viz., stearic (18:0), cis and trans isomers of oleic (18:1), linoleic (18:2), and linolenic (18:3) acids, as well as to cholesterol. 4 Fatty acids were shown to be stronger bound to the DNA minor groove than to the major groove. The energy of interactions between fatty acids and DNA was found to depend on the number of double bonds and the isomerism of fatty acids, as well as on the nucleotide composition of DNA. In another series of studies, we calculated the bond energies of phos phatidylcholine or sphingomyeline to 64 oligodeoxyribo nucleotides corresponding to triplets of the genetic code. 5,6 It was demonstrated that GC rich triplets form more stable complexes with phospholipids and that the complexes with sphingomyeline are characterized by higher energies.The present study continues our research on the com plexation of DNA with lipids by molecular mechanics methods. The aim was to investigate the character of in teractions and the bond energies of DNA to diglycerides containing stearic, trans oleic, cis oleic, trans,trans li noleic, cis,cis linoleic, trans,trans,trans linolenic, and cis,cis,cis linolenic acids.
Calculation procedureAs in our earlier studies, 1-6 all calculations were per formed by the molecular mechanics method using the MM + force field and the Hyperchem 5.01 program package. In this method, the total energy E total of systems (complexes) is calculated as the sum of the contributions from bond length deformations, bond angle deformations, the torsional strain, the electrostatic interaction energy, the van der Waals (nonbonded) interaction energy, and the hydrogen bond energy. The full geometry optimization was carried out for fragments of the DNA double strand, diglycerides (Fig. 1), and the oligodeoxyribonucleotide (ON)-ligand complexes. The bond energies E bond of the ligands with ON were evaluated. The B form of ON was considered, and sodium ions were added to fulfill the electroneutrality condition.The bond energies of ON with the ligands were evalu ated as the difference between the total energy of the ON-diglyceride complex and the sum of the energies of the isolated components of the complex. Taking into account the geometric features of the DNA B form in terms of the double stranded model, the ligand (diglyceride) molecules were located in the minor or major gr...