Interactions of molecules with nucleic acids. XI. Generalization of techniques to generate nucleic acid structures with applications to intercalation sites and kinked structures

The study of possible A and B conformations of alternating DNA using a new program for conformational analysis of duplexes (CONAN)

A new program, CONAN has been designed for CONformational ANalysis of oligonucleotide duplexes with natural and modified bases. It allows to model both regular DNA fragments with different types of symmetry and irregular ones including bends, junctions, mismatched pairs and base lesions. Computations and minimization of the energy are performed in a space of internal structural variables chosen to build start structure easier and conveniently analyze the results obtained. These internal structural variables determine mutual base-base and base-sugar arrangement and sugar puckering. The analytical closure procedure is applied both to sugar rings and to backbone fragments between adjacent sugars. For more effective energy minimization, analytical gradient is calculated. The CONAN was applied to the search for low-energy conformations of poly(dA-dT).poly(dA-dT) and poly(dG-dC).poly(dG-dC) duplexes. Extended regions of low-energy A and B conformations are revealed and characterized. These regions contain structures with different relative values of helical twist, tau, for pur-pyr and pyr-pur steps, namely, conformations with tau (pur-pyr) > tau (pyr-pur) and with tau (pur-pyr) < tau (pyr-pur). Two types of sugar puckering were found for B-form low-energy conformations, the first type with all C2'-endo sugar residues and the second one - with C2'-endo purines and O1'-endo pyrimidines. The calculated conformations are compared with X-ray diffraction data for crystals and fibers and NMR data for solution.

Aflatoxin B1 and DNA adducts. Proposed model for surface noncovalent and covalent complexes with N(7) of guanine. II

An alternate model for surface noncovalent and surface covalent binding of aflatoxin B1 to N(7) of guanine in DNA is proposed. This model considers the out-of-plane motions of C(8) of aflatoxin B1 in those interactions. The covalent intercalated fit of aflatoxin B1 into DNA arises from steric adjustments made by DNA at the covalent intercalation site as well as local strain in the bond angles about N(7) of guanine and C(8) of aflatoxin B1. The bond angle about N(7) deviates modestly from the sp2 value toward the sp3 value. This study suggests that the surface covalent aflatoxin B1-DNA complex serves only a minor role in aflatoxin's precarcinogenic interaction with DNA and is a likely correctable error.

Theoretical studies of cis-Pt(II)-diammine binding to duplex DNA

The binding of cis-Pt(II) diammine (cis-DP) to double-stranded DNA was studied with several kinked conformations that can accommodate the formation of a square planar complex. Molecular mechanics (MM) calculations were performed to optimize the molecular fit. These results were combined with quantum mechanical (QM) calculations to ascertain the relative energetics of ligand binding through water vs direct binding of the phosphate to the ammine and platinum, and to guide the selection of DNA conformations to model complex formation. Based on QM and MM calculations, models are proposed that may be characterized by several general features. A structure involving hydrogen bonding between each ammine and distinct adjacent phosphate groups, referred to as closed conformation (CC), has already been reported. This is also found in the crystal structure of small dimers. We report alternative conformations that may be important in platination of duplex DNA. They are characterized by an intermediate conformation (IC), involving hydrogen bonding between one ammine and phosphate group, and an open conformation (OC), without ammine phosphate hydrogen bonding. The IC and OC can be stabilized by water bridges in the space between the ammine and the phosphate groups. Sugar puckers alternate from the type C(2')-endo or C(1')-exo (S), to the type C(3')-endo or C(2')-exo (N), with intermediate types near O(1')-endo (O). In general, the sugar puckers alternate from S to N to S through the platinated region (3'-TpG*pG*p-5'), with the complexed strand exhibiting, (3')-S*-N*-S-(5') alternation, while the complementary strand shows either (3')-S*-N*-S-(5') or (3')-S*-N*-O-(5') alternation. In both the OC and IC, a hydrogen bond is found between the ammine and O4(T) on thymine (T) at the (3') end, adjacent to the complex site. There is a continuous range of backbone conformations through the platinated region which relate the OC to the IC. The models presented suggest that the dynamics of the binding of the cis-Pt(II)-diammines to adjacent N7(G) in double-stranded DNA may encompass several conformational possibilities, and that water bridges may play a roll in supporting open and intermediate conformations. Proton-proton distances are reported to assist in the experimental determination of conformations.

Distances as degrees of freedom

Tertiary contact distance information of varying resolution for large biological molecules abounds in the literature. The results provided herein develop a framework by which information of this type can be used to reduce the allowable configuration space of a macromolecule. The approach combines graph theory and distance geometry. Large molecules are represented as simple, undirected graphs, with atoms, or groups, as vertices, and distances between them as edges. It is shown that determination of the exact structure of a molecule in three dimensions only requires the specification of all the distances in a single tetrahedron, and four distances to every other atom. This is 4N-10 distances which is a subset of the total N(N-1)/2 unique distances in a molecule consisting of N atoms. This requirement for only 4N-10 distances has serious implications for distance geometry implementations in which all N(N-1)/2 distances are specified by bounded random numbers. Such distance matrices represent overspecified systems which when solved lead to non-obvious distribution of any error caused by inherent contradictions in the input data. It is also shown that numerous valid subsets of 4N-10 distances can be constructed. It is thus possible to tailor a subset of distances using all known distances as degrees of freedom, and thereby reduce the configuration space of the molecule. Simple algebraic relationships are derived that relate sets of distances, and complicated rotations are avoided. These relationships are used to construct minimum, complete sets of distances necessary to specify the exact structure of the entire molecule in three dimensions from incomplete distance information, and to identify sets of inconsistent distances. The method is illustrated for the flexible structural types present in large ribosomal RNAs: 1.) A five-membered ring; 2.) a chemically bonded chain with its ends in contact (i.e., a hairpin loop); 3.) the spatial orientation of two separate molecules, and; 4.) an RNA helix that can have variation in individual base pairs, giving rise to global deviation from standardized helical forms.

The structure of the aflatoxin B1-DNA adduct at N7 of guanine. Theoretical intercalation and covalent adduct models

Two theoretical models are proposed for the conformational structure of both intercalated and covalent adduct complexes of aflatoxin B1, designated AFB1, with N7 of guanine of DNA. The covalent adduct model requires the DNA to kink a minimum of 39 degrees about the covalent site of the C8 and N7 atoms comprising the bond of the covalent complex. The preference of AFB1 for specific G bases within a sequence of GC content followed that of experimental studies with the added feature that for binding to the third G base of a tetramer sequence from the 3'-end, the AFB1 displayed enhanced binding at the 3' site of the targeted guanosine. Binding of AFB1 to the second G base of a tetramer sequence from the 3'-end leads to preference for a 5' site of the targeted guanosine. Inhibition of AFB1's interaction with the targeted DNA in the presence of intercalated ethidium bromide is explained by these proposed models.

Conformational studies of (2'-5') polynucleotides: theoretical computations of energy, base morphology, helical structure, and duplex formation

A detailed theoretical analysis has been carried out to probe the conformational characteristics of (2'-5') polynucleotide chains. Semi-empirical energy calculations are used to estimate the preferred torsional combinations of the monomeric repeating unit. The resulting morphology of adjacent bases and the tendency to form regular single-stranded structures are determined by standard computational procedures. The torsional preferences are in agreement with available nmr measurements on model compounds. The tendencies to adopt base stacked and intercalative geometries are markedly depressed compared to those in (3'-5') chains. Very limited families of regular monomerically repeating single-stranded (2'-5') helices are found. Base stacking, however, can be enhanced (but helix formation is at the same time depressed) in mixed puckered chains. Constrained (2'-5') duplex structures have been constructed from a search of all intervening glycosyl and sugar conformations that form geometrically feasible phosphodiester linkages. Both A- and B-type base stacking are found to generate non-standard backbone torsions and mixed glycosyl/sugar combinations. The 2'- and 5'-residues are locked in totally different arrangements and are thereby prevented from generating long helical structures.

The flexibility of the nucleic acids: (I). "SIR", a novel approach to the variation of polymer geometry in constrained systems

A novel and powerful methodology is developed which allows the alteration of molecular structures subjected to constraints and its application to polynucleotides with mononucleotide repeat symmetry, including the treatment of the flexible sugar rings is described. In contrast to procedures proposed by other authors, the constraints are formulated as differential equations which are linear with respect to the differentials of the geometrical variables. These equations can be solved easily by stepwise numerical integration involving sucessive infinitesimal rotations (SIR). Moreover, these equations define a set of independent curvilinear coordinates which can be used directly as the parameters of the energy functional in an energy minimisation procedure. This methodology allows the scanning of the full configurational space of a complex macromolecule, with direct access to the helicoidal variables in the case of periodic systems. Through this approach many problems involving biomacromolecular conformation, which would otherwise be intractable, may be studied with considerable ease.

A theoretical model for the binding of cis-Pt(NH3)2(+2) to DNA

The binding of cis-Pt(NH3)2B1B2 to the bases B1 and B2, i.e., guanine (G), cytosine (C), adenine (A), and thymine (T), of DNA is studied theoretically. The components of the binding are analyzed and a model structure is proposed for the intrastrand binding to the dB1pdB2 sequence of a kinked double helical DNA. Quantum mechanical calculations of the ligand binding energy indicates that cis-Pt(NH3)2(+2) (cis-PDA) binds to N7(G), N3(C), O2(C), O6(G), N3(A), N7(A), O4(T) and O2(T) in order of decreasing binding energy. Conformational analysis provides structures of kinked DNA in which adjacent bases chelate to cis-PDA. Only bending toward the major groove allows the construction of acceptable square planar complexes. Examples are presented for kinks of -70 degrees and -40 degrees at the receptor site to orient the base pairs for ligand binding to B1 and B2 to form a nearly square planar complex. The energies for complex formation of cis-PDA to the various intra-strand base sites in double stranded DNA are estimated. At least 32 kcal/mole separates the energetically favorable dGpdG.cis-PDA chelate from the dCpdG.cis-PDA chelate. All other possible chelate structures are much higher in energy which correlates with their lack of observation in competition with the preferred dGpdG chelate. The second most favorable ligand energy occurs with N3(C). A novel binding site involving dC(N3)pdG(N7) is examined. Denaturation can result in an anti----syn rotation of C about its glycosidic bond to place N3(C) in the major groove for intrastrand binding in duplex DNA. This novel intrastrand dCpdG complex and the most favored dGpdG structure are illustrated with stereographic projections.