Molecular Mechanics Studies of Homopolymeric and Mixed Sequence Py.Pu*Pu Triple Helices,

DNA triple helices containing two purine strands and one pyrimidine strand (C.G*G and T.A*A) have been studied, using model building followed by energy minimisation, for different orientations of the third strand resulting from variation in the hydrogen bonding between the Watson-Crick duplex and the third strand and the glycosidic torsion angle in the third strand. Our results show that in the C.G*G case the structure with a parallel orientation of the third strand, resulting from Hoogsteen hydrogen bonds between the third strand and the Watson-Crick duplex, is energetically the most favourable while in the T.A*A case the antiparallel orientation of the third strand, resulting from reverse Hoogsteen hydrogen bonds, is energetically the most favourable. These studies when extended to the mixed sequence triplexes, in which the second strand is a mixture of G and A, correspondingly the third strand is a mixture of G and APT, show that though the parallel orientation is still energetically more favourable, the antiparallel orientation becomes energetically comparable with an increasing number of thymines in the third strand. Structurally, for the mixed triplexes containing G and T in the third strand, it is seen that the basepair non-isomorphism between the C.G*G and the T.A*T triplets can be overcome with some changes in the base pair parameters without much distortion of either the backbone or the hydrogen bonds.

Structural Polymorphism in d(T).d(A)*d(T) Triple Helices

DNA triple helices containing two thymine strands and one adenine strand have been studied, using model building followed by energy minimisation, for different orientations of the third strand resulting from variation in the hydrogen bonding between the Watson-Crick duplex and the third strand and the glycosidic torsion angle in the third strand. Our results show that the structure with a parallel orientation of the third strand, in which the third strand base forms Hoogsteen hydrogen bonds with the adenine base in the Watson-Crick duplex, is energetically the most favourable. An antiparallel orientation of the third strand is also possible, in which the third strand base hydrogen bonds to both the bases in the Watson-Crick duplex. This structure is energetically comparable to the parallel structure. For the parallel triplex a 200ps molecular dynamics simulation starting from two different starting structures indicates that at 300K significant structural heterogeneity exists in this tripler structure. The results are compared with existing structural data on this class of triplexes derived from theoretical and NMR techniques.

Analysis of sequence dependent variations in secondary and tertiary structure of tRNA molecules

The double helical regions of the five tRNA(Phe) and two tRNA(Asp) crystal structures have been analyzed using the local basepair step parameters. The sequence dependent effects in the mini double helices of tRNA are very similar to those observed in the crystal structures of oligonucleotides in the A-form, the purine-pyrimidine and purine-purine steps have small roll angles when compared to the fiber models of A-DNA as well as A-RNA, while the pyrimidine-purine doublet steps have large roll angles. The orientation of the basepairs in the D-stem is unusual and invariant i.e. they are different from the other three stems but are very similar in all the five tRNA(Phe) crystal structures, presumably due to tertiary interaction of the Watson-Crick basepairs with other bases, with all bases being highly conserved. The origin of the differences between the tertiary structures of tRNA(Phe) and tRNA(Asp) from yeast has also been investigated. It is found that even though the angle between the acceptor arm and the D-stem is very similar in the two structures, the angle subtended by the acceptor arm and the anticodon arm is smaller in the tRNA(Phe) structure (by more than 10 degrees). This is due to differences in the orientation of the two mini helices constituting the anticodon arm, which are inclined to each other by approximately 25 degrees in tRNA(Phe) and 16 degrees in tRNA(Asp). In addition, the acceptor arm, the D-stem and the anticodon stem are nearly coplanar in tRNA(Phe), while in tRNA(Asp) the anticodon stem projects out of the plane defined by the acceptor arm and the anticodon stem. These two features together lead to a larger separation between the acceptor and anticodon ends in tRNA(Asp) and indicate that the junction between the D-stem and the anticodon stem is quite variable, with features characteristic of a ball-and-socket type joint and determined for each tRNA molecule by the base sequence at the junction.

Chain folding and A:T pairing in human telomeric DNA: a model-building and molecular dynamics study

The various types of chain folding and possible intraloop as well as interloop base pairing in human telomeric DNA containing d(TTAG(3)) repeats have been investigated by model-building, molecular mechanics, and molecular dynamics techniques. Model-building and molecular mechanics studies indicate that it is possible to build a variety of energetically favorable folded-back structures with the two TTA loops on same side and the 5' end thymines in the two loops forming TATA tetrads involving a number of different intraloop as well as interloop A:T pairing schemes. In these folded-back structures, although both intraloop and interloop Watson-Crick pairing is feasible, no structure is possible with interloop Hoogsteen pairing. MD studies of representative structures indicate that the guanine-tetraplex stem is very rigid and, while the loop regions are relatively much more flexible, most of the hydrogen bonds remain intact throughout the 350-ps in vacuo simulation. The various possible TTA loop structures, although they are energetically similar, have characteristic inter proton distances, which could give rise to unique cross-peaks in two-dimensional nuclear Overhauser effect spectroscopy (NOESY) experiments. These folded-back structures with A:T pairings in the loop region help in rationalizing the data from chemical probing and other biochemical studies on human telomeric DNA.

Implications of the recent high statistics determination of the pion electromagnetic form factor in the timelike region

The recently evaluated two-pion contribution to the muon g - 2 and the phase of the pion electromagnetic form factor in the elastic region, known from pi pi scattering by Fermi-Watson theorem, are exploited by analytic techniques for finding correlations between the coefficients of the Taylor expansion at t = 0 and the values of the form factor at several points in the spacelike region. We do not use specific parametrizations, and the results are fully independent of the unknown phase in the inelastic region. Using for instance, from recent determinations, < r(pi)(2)> = (0.435 +/- 0.005) fm(2) and F(-1.6 GeV2) = 0.243(-0.014)(+0.022), we obtain the allowed ranges 3.75 GeV-4 less than or similar to c less than or similar to 3.98 GeV-4 and 9.91 GeV-6 less than or similar to d less than or similar to 10.46 GeV-6 for the curvature and the next Taylor coefficient, with a strong correlation between them. We also predict a large region in the complex plane where the form factor cannot have zeros.

Mutagenic significance of proton acidities in methylated guanine and thymine bases and deoxynucleosides: A theoretical study

Proton changes have been advanced as being the key molecular basis for the mutagenecity of alkylated DNA bases and nucleosides, leading to questions as to which protons are involved and whether the protic changes are tautomeric shifts or abstractions. This semiempirical molecular orbital study seeks to clarify the issue by examining the various possibilities open for these protic changes in a number of methylated guanines and thymines and their deoxynucleosides. Proton shifts leading to tautomer formation are not predicted as being thermodynamically favourable in most cases. The most feasible proton abstractions are predicted to involve the Watson-Crick protons in all cases, which corroborates Watson-Crick proton loss as providing the key molecular basis for the induction of point mutations. The calculated proton acidities correlate well with experimental data. The gas-phase deprotonation enthalpies for a number of alkylated nucleosides are found to correlate linearly with the solvent-phase pK(a) values. The theoretically calculated enthalpies in a simulated aqueous solvent phase of the deprotonation reactions of various nucleic acid bases are also found to have good linear correlations with experimental pK(a) values. The consensus of these calculations is that O-6-alkyldeoxyguanosines, and O-2- and O-4-alkyldeoxythymidines would be mutagenic while N-7-alkyldeoxyguanosines would not be mutagenic (as experiment indicates). The untested N-3-methyldeoxyguanosine is predicted to be mutagenic. (C) 1997 Elsevier Science B.V.

Conjugacy invariant pseudo-norms, representability and RNA secondary structures

To a reasonable approximation, a secondary structures of RNA is determined by Watson-Crick pairing without pseudo-knots in such a way as to minimise the number of unpaired bases: We show that this minimal number is determined by the maximal conjugacy-invariant pseudo-norm on the free group on two generators subject to bounds on the generators. This allows us to construct lower bounds on the minimal number of unpaired bases by constructing conjugacy invariant pseudo-norms. We show that one such construction, based on isometric actions on metric spaces, gives a sharp lower bound. A major goal here is to formulate a purely mathematical question, based on considering orthogonal representations, which we believe is of some interest independent of its biological roots.