Simulation of interactions between nucleic acid bases by refined atom-atom potential functions

Energy of interaction between nitrogen bases of nucleic acid has been calculated as a function of parameters determining the mutual position of two bases. Refined atom-atom potential functions are suggested. These functions contain terms proportional to the first (electrostatics), sixth (or tenth for the atoms forming a hydrogen bond) and twelfth (repulsion of all atoms) powers of interatomic distance. Calculations have shown that there are two groups of minima of the base interaction energy. The minima of the first group correspond to coplanar arrangement of the base pairs and hydrogen bond formation. The minima of the second group correspond to the position of bases one above the other in almost parallel planes. There are 28 energy minima corresponding to the formation of coplanar pairs with two (three for the G:C pair) almost linear N-H . . . O and (or) N-H . . . N hydrogen bonds. The position of nitrogen bases paired by two such H-bonds in any crystal of nucleic acid component in polynucleotide complexes and in tRNA is close to the position in one of these minima. Besides, for each pair there are energy minima corresponding to the formation of a single N-H . . . O or N-H . . . N and one C-H . . . O or C-H . . . N hydrogen bond. The form of potential surface in the vicinity of minima has been characterized. The results of calculations agree with the experimental data and with more rigorous calculations based on quantum-mechanical approach.     

Simulation of Interactions between Nucleic Acid Bases by Refined Atom-Atom Potential Functions

Abstract

Energy of interaction between nitrogen bases of nucleic acids has been calculated as a function of parameters determining the mutual position of two bases. Refined atom-atom potential functions are suggested. These functions contain terms proportional to the first (electrostatics), sixth (or tenth for the atoms forming a hydrogen bond) and twelfth (repulsion of all atoms) powers of interatomic distance. Calculations have shown that there are two groups of minima of the base interaction energy. The minima of the first group correspond to coplanar arrangement of the base pairs and hydrogen bond formation. The minima of the second group correspond to the position of bases one above the other in almost parallel planes. There are 28 energy minima corresponding to the formation of coplanar pairs with two (three for the G:C pair) almost linear N-H … O and (or) N-H … N hydrogen bonds. The position of nitrogen bases paired by two such H-bonds in any crystal of nucleic acid component, in polynucleotide complexes and in tRNA is close to the position in one of these minima. Besides, for each pair there are energy minima corresponding to the formation of a single N-H … O or N-H … N and one C-H … O or C-H … N hydrogen bond. The form of potential surface in the vicinity of minima has been characterized. The results of calculations agree with the experimental data and with more rigorous calculations based on quantum- mechanical approach.     

The linear dichroism of oriented helical and superhelical polymers

Expressions are derived to relate the linear dichroism of oriented helical and super-helical polymers to the vector equation describing a transition moment in a reference frame based on the position of a single monomer unit in a simple helix and the super-helix tilt angle β. The reduced dichroism of a perfectly oriented superhelical polymer, (Δε/ε)₀, is related to the reduced dichroism of a perfectly oriented helical polymer, (Δε/ε)₁, by the expression (Δε/ε)₀ = ½ (Δε/ε)₁(3 cos² β ? 1). Hence the linear dichroism of a helical polymer can be significantly altered by supercoiling. The application of the general expressions is illustrated by the derivation of specific equations which relate the dichroism of helical and superhelical DNA to the tilt and twist of the DNA bases, the directions of the transition moments within the base planes, and the superhelix tilt angle. Calculations of the reduced dichroism of perfectly oriented, non-super-helical DNA in the A and B configurations suggest that it will be difficult to distinguish between these two forms on the basis of dichroism measurements.     

DNA Modeller: an interactive program for modelling stacks of DNA base pairs on a microcomputer

DNA Modeller is a microcomputer program for interactively manipulatingup to 20 bp in a DNA double helical arrangement. It calculatesthe van der Waals and electrostatic energies of base-base interactionsusing the AMBER potential, minimizes the energy with respectto the pair (buckle, propeller, opening, shear, stretch, stagger)and step (tilt, roll, twist, shift, slide, rise) parameters,calculates lengths of the canonical hydrogen bonds between thecomplementary bases, and calculates interatomic distances betweenthe successive base pairs. Input/output files are simple listsof the step and pair parameters or lists of the atom specifications(N1, C2, etc.) and their Cartesian coordinates (compatible withthe Desktop Molecular Modeller ^*.mol files). The program issupplied with a readbrk utility which transforms PDB/NDB tothe^*.mol format readable by DNA Modeller. The DNA crystal structuresdeposited in the PDB or NDB databases can thus be analyzed,and their bases visualized and interactively manipulated. Inaddition, DNA Modeller can calculate the base pair and stepgeometrical parameters and interaction energies. A plotter utilitycreates wire mono or stereo pictures of the bases. This programis designed for IBM-compatible computers working under DOS orcan run as a DOS application under MS Windows 3.x or Merge (SCOUnix DOS emulator).     

Multiplicity of cadmium binding sites in nucleotides: X-ray evidence for the involvement of O2'' and O3'' as well as phosphate and N7 in inosine 5''-monophosphate.

Single-crystal X-ray methods have been used to characterize a hydrated polymeric cadmium derivative of inosine 5'-monophosphate. In the structure there are two independent cadmium atoms, one of which binds to two ribose oxygen atoms, an N7 position on a base, and to three water molecules. The second metal atom binds to a phosphate oxygen, three water molecules, and to two N7 atoms, which are in cis-positions. For these last; the Cd-N bonds are appreciably out of the planes of the hypoxanthine bases so that the angle between these planes is only 31.4 degrees.     

Interaction of nucleic acid bases with water molecules and formation of mismatched nucleotide pairs

The possibility of the inclusion of water molecules in the formation of mismatched nucleotide pairs was considered in relation to the mechanisms of point errors in template directed biosynthesis. Calculations of the intermolecular interaction energy for systems containing two bases and one water molecule were carried out by the method of atom-atom potential functions. There exist energy minima for each base pair, corresponding to a single N--H...O or N--H...N H-bond between the bases and H-bonding of the water molecule with both bases. The relative positions of glycosyl bonds in some of these minima are closer to those for Watson--Crick pairs, than the positions of minima for these pairs without water. For other minima, the H-bond formation between the water molecule and the two bases additionally stabilizes the relative base position in wobble-pairs with two H-bonds between the bases. The base and water positions in energy minima are compared with the positions in some pairs proposed on the basis of NMR and X-ray data for double helical oligonucleotides.     

Recognition of tRNAs by aminoacyl-tRNA synthetases: Escherichia coli tRNAMet and E. coli methionyl-tRNA synthetase

In previous work we identified several specific sites in Escherichia coli tRNAfMet that are essential for recognition of this tRNA by E. coli methionyl-tRNA synthetase (MetRS) (EC 6.1.1.10). Particularly strong evidence indicated a role for the nucleotide base at the wobble position of the anticodon in the discrimination process. We have now investigated the aminoacylation activity of a series of tRNAfMet derivatives containing single base changes in each position of the anticodon. In addition, derivatives containing permuted sequences and larger and smaller anticodon loops have been prepared. The variant tRNAs have been enzymatically synthesized in vitro by using T4 RNA ligase (EC 6.5.1.3). Base substitutions in the wobble position have been found to reduce aminoacylation rates by at least five orders of magnitude. Derivatives having base substitutions in the other two positions of the anticodon are aminoacylated 55-18,500 times slower than normal. Nucleotides that have specific functional groups in common with the normal anticodon bases are better tolerated at each of these positions than those that do not. A tRNAfMet variant having a six-membered loop containing only the CA sequence of the anticodon is aminoacylated still more slowly, and a derivative containing a five-membered loop is not measurably active. The normal loop size can be increased by one nucleotide with a relatively small effect on the rate of aminoacylation, which indicates that the spatial arrangement of the nucleotides is less critical than their chemical nature. We conclude from these data that recognition of tRNAfMet requires highly specific interactions of MetRS with functional groups on the nucleotide bases of the anticodon sequence. Several other aminoacyl-tRNA synthetases are known to require one or more anticodon bases for efficient aminoacylation of their tRNA substrates, and data from other laboratories suggest that anticodon sequences may be important for accurate discrimination between cognate and noncoagnate tRNAs by these enzymes.     

Solution structure of a five-adenine bulge loop within a DNA duplex.

The three-dimensional solution structure of a DNA molecule of the sequence 5'-d(GCATCGAAAAAGCTACG)-3' paired with 5'-d(CGTAGCCGATGC)-3' containing a five-adenine bulge loop (dA(5)-bulge) between two double helical stems was determined by 2D (1)H and (31)P NMR, infrared, and Raman spectroscopy. The DNA in both stems adopt a classical B-form double helical structure with Watson-Crick base pairing and C2'-endo sugar conformation. In addition, the two dG/dC base pairs framing the dA(5)-bulge loop are formed and are stable at least up to 30 degrees C. The five adenine bases of the bulge loop are localized at intrahelical positions within the double helical stems. Stacking on the double helical stem is continued for the first four 5'-adenines in the bulge loop. The total rise (the height) of these four stacked adenines roughly equals the diameter of the double helical stem. The stacking interactions are broken between the last of these four 5'-adenines and the fifth loop adenine at the 3'-end. This 3'-adenine partially stacks on the other stem. The angle between the base planes of the two nonstacking adenines (A10 and A11) in the bulge loop reflects the kinking angle of the global DNA structure. The neighboring cytosines opposite the dA(5)-bulge (being parts of the bulge flanking base pairs) do not stack on one another. This disruption of stacking is characterized by a partial shearing of these bases, such that certain sequential NOEs for this base step are preserved. In the base step opposite the loop, an extraordinary hydrogen bond is observed between the phosphate backbone of the 5'-dC and the amino proton of the 3'-dC in about two-thirds of the conformers. This hydrogen bond probably contributes to stabilizing the global DNA structure. The dA(5)-bulge induces a local kink into the DNA molecule of about 73 degrees (+/-11 degrees ). This kinking angle and the mutual orientation of the two double helical stems agree well with results from fluorescence resonance energy transfer measurements of single- and double-bulge DNA molecules.     

Structure of a B-DNA dodecamer. II. Influence of base sequence on helix structure

Detailed examination of the structure of the B-DNA dodecamer C-G-C-G-A-A-T-T-C-G-C-G, obtained by single-crystal X-ray analysis (Drew et al., 1981), reveals that the local helix parameters, twist, tilt and roll, are much more strongly influenced by base sequence than by crystal packing or any other external forces. The central EcoRI restriction endonuclease recognition site, G-A-A-T-T-C, is a B helix with an average of 9.8 base-pairs per turn. It is flanked on either side by single-base-pair steps having aspects of an A-like helix character. The dodecamer structure suggests several general principles, whose validity must be tested by other B-DNA analyses. (1) When an external bending moment is applied to a B-DNA double helix, it bends smoothly, without kinks or breaks, and with relatively little effect on local helix parameters. (2) Purine-3′,5′-pyrimidine steps open their base planes towards the major groove, pyrimidine-purine steps open toward the minor groove, and homopolymer (Pur-Pur, Pyr-Pyr) steps resist rolling in either direction. This behavior is related to the preference of pyrimidines for more negative glycosyl torsion angles. (3) CpG steps have smaller helical twist angles than do GpC, as though in compensation for their smaller intrinsic base overlap. Data on A-T steps are insufficient for generalization. (4) G.C base-pairs have smaller propellor twist than A · T, and this arises mainly from interstrand base overlap rather than the presence of the third hydrogen bond. (5) DNAase I cuts preferentially at positions of high helical twist, perhaps because of increased exposure of the backbone to attack. The correlation of the digestion patterns in solution and helical twist in the crystal argues for the essential identity of the helix structure in the two environments. (6) In the two places where the sequence TpCpG occurs, the C slips from under T in order to stack more efficiently over G. At the paired bases of this CpG step, the G and C are tilted so the angle between base planes is splayed out to the outside of the helix. This TpC is the most favored cutting site for DNAase I by a factor of 4.5 (Lomonossoff et al., 1981). (7) The EcoRI restriction endonuclease and methylase both appear to prefer a cutting site of the type purine-purine-A-T-T-pyrimidine, involving two adjacent homopolymer triplets, and this may be a consequence of the relative stiffness of homopolymer base-stacking observed in the dodecamer.     

Theoretical conformational analysis of doulble-stranded polynucleotides]

Calculations of intramolecular interaction energy of two-stranded helical homopolynucleotide in the function of nine conformational variables have been carried out by the method of atom-atom potential functions. Four of these variables determine mutual position of base pairs, other four--deoxiribose ring conformation and other one--orientation of this ring with respect to the base. For this purpose an algorythm connecting dependent variables with independent ones has been developed. The investigation of energy function has shown that in the space of conformational parameters there are two valleys, which correspond to A-and B-families of conformations. Experimentaly determined conformations of two-stranded helical polynucleotides are located along the bottoms of these valleys. Along the bottom of each valley the intramolecular interaction energy changes rather little when conformational parameters change within a wide range. The valleys are separated by an energetical barrier.     

Base interaction and conformation manifestations of repetitive nucleotide sequences

To understand why different nucleotide sequences prefer different double helical conformations and to predict conformational behaviour of definite sequences the base-base interaction energy in regular helices consisting of A:U, A:T, G:C and I:C (hypoxanthine-cytosine) base pairs was calculated. Interaction energy was assumed to be a function of eight conformational parameters: H, the distance between adjacent pairs along helix axes; tau, turn angle of one pair relative to the neighbouring one; angles between base planes in a pair (TW, propeller twist and BL, buckle) and position of pairs with respect to helix axes (D and SL, displacements in the plane normal to helix axes, and TL and RL, inclinations to this plane, tilt and roll, respectively). For H and tau characteristic of A- and B-families of nucleic acid conformations (2.5 A less than H less than or equal to 3.5 A, 30 degrees less than or equal to tau less than or equal to 45 degrees) the ranges of conformational parameters corresponding to energy values close to minimal ones (valleys) and correlations between conformational parameters were revealed. Valleys for different sequences largely coincide but have distinctive characteristics for each sequence. Reasons for base pair planarity distortion in double-stranded helices were considered. The calculations permit to account for A-phility of G:C sequences and B-phility for A:T sequences. The valley for I:C sequence branches. This corresponds to A:T-like behaviour in some cases and G:C-like in the others.     

Self base pairing in a complementary deoxydinucleoside monophosphate duplex: crystal and molecular structure of deoxycytidylyl-(3'-5')-deoxyguanosine

The molecular structure of ammonium deoxycytidylyl-(3'-5')-deoxyguanosine, crystallized from aqueous acetone near pH 4, was determined for X-ray diffraction data. The crystals were tetragonal, space group P43212 with a = b = 11.078 (1) A and c = 45.826 (4) A. The structure was solved by tangent expansion of phases based on a derived phosphorus position and refined to R = 0.060 by full matrix least squares. Molecules related by a 2-fold symmetry axis are connected by hydrogen bonds between the bases and form parallel right-handed duplexes. Pairs of cytosines share a proton at N(3) and are joined by three hydrogen bonds: N(4)-H...O(2)...H-N(4), and N(3)-H...N(3). Guanines are joined by two hydrogen bonds: N(2)-H...N(3) and N(3)...H-N(2). Base-stacking interactions within the duplex are weak with the cytosine and guanine ring planes inclined at 24 degrees to each other in each monomer. Despite the unusual arrangement of the molecules, the sugar phosphate backbone has the g-g- conformation normally associated with right-handed double helical structures. Conformational parameters of the nucleosides are also typical with both sugars C(2')-endo and glycosidic torsion angles 55 degrees for cytidine and 94 degrees for guanosine. The bonding geometry of the bases is influenced by hydrogen bonding and charge-transfer networks in the crystal lattice. The solvent molecules interact with the dimer in three fused circular hydrogen bonding domains with a single disordered ammonium cation per d(CpG) dimer. Parallels with the formation of self base pairs and their implications in molecular biology are discussed.     

Isolated transmembrane helices arranged across a membrane: computational studies.

A computational procedure for predicting the arrangement of an isolated helical fragment across a membrane was developed. The procedure places the transmembrane helical segment into a model triple-phase system 'water-octanol-water'; pulls the segment through the membrane, varying its 'global' position as a rigid body; optimizes the intrahelical and solvation energies in each global position by 'local' coordinates (dihedral angles of side chains); and selects the lowest energy global position for the segment. The procedure was applied to 45 transmembrane helices from the photosynthetic reaction center from Rhodopseudomonas viridis, cytochrome c oxidase from Paracoccus denitrificans and bacteriorhodopsin. In two thirds of the helical fragments considered, the procedure has predicted the vertical shifts of the fragments across the membrane with an accuracy of -0.15 +/- 3.12 residues compared with the experimental data. The accuracy for the remaining 15 fragments was 2.17 +/- 3.07 residues, which is about half of a helix turn. The procedure predicts the actual membrane boundaries of transmembrane helical fragments with greater accuracy than existing statistical methods. At the same time, the procedure overestimates the tilt values for the helical fragments.     

Nucleic Acid Structure Analysis: A Users Guide to a Collection of New Analysis Programs

Abstract

Common nomenclature describing the geometry of nucleic acid structures was established at a 1988 EMBO Workshop on DNA Curvature and Bending (Diekmann, S. (1988) J. Mol. Biol. 208, 787–791; Diekmann, S. (1989) The EMBO Journal 8, 1–4; Sarma, RH. (1988) J. Biomol. Structure & Dynamics 6, 391–395; Dickerson, R.E. (1989) J. Biomol. Structured Dynamics 6, 627–634; Dickerson, RE. et al. (1989)Nuc. Acids Res. 17, 1979–1803). We have subsequently developed and incorporated sophisticated mathematics in a computer program to calculate the parameters described by the guidelines. The program calculates all the local parameters relating complementary bases and neighboring base and base pairs in both Cartesian and helical coordinate frames. In addition, the main mathematical property requested by the EMBO guidelines—that the magnitude of the parameters be independent of strand or direction of measurement—is accomplished without the use of a midway coordinate frame for the rotational parameters. The mathematics preserve the physical intuition used in defining the parameters; in particular, the rotational parameters are true rotations based on a simple physical model (rotation at constant angular velocity for a unit amount of time), not Euler angles or angles between vectors and planes as is the case with other approaches. As a result, the mathematical equations are symmetric with the property that a 5° tilt is the same as a 5° roll or a 5° twist, except that the rotations take place about different axes. In other approaches, a 5° tilt can mean a different amount of net rotation from a 5° roll or a 5° twist. In addition, a great deal of flexibility is built into the program so that the user has control over the analysis, including the input format, the coordinate frame used for the base pairing relationship, the point about which the rotations are performed, and which geometric relationships are analyzed. While there is a great deal of flexibility, the program is easy to use. Interactive queries and user accessible files make the options in the program very convenient to tailor to individual needs. In addition, there is also a program that calculates bond lengths, valence angles, and torsion angles along the nucleic acid backbone, and within the sugar and base rings. Another program ‘learns’ the identities of the bond lengths, valence angles, and torsion angles that the user would like to determine. This last program is especially useful for calculating the hydrogen bonds between atoms in complementary strands as well as the unusual hydrogen bonds found in recently determined nucleic acid NMR structures or within protein/nucleic acid complexes.     

Packing features of the all-AT oligonucleotide d(AAATTT)

We describe the packing features of the oligonucleotide duplex d(AAATTT)2, as determined by X-ray diffraction. There is little information on sequences that only contain A and T bases. The present structure confirms that these sequences tend to pack as a helical arrangement of stacked oligonucleotides in a B conformation with Watson-Crick hydrogen bonding. Our results demonstrate that the virtual TA base step between stacked duplexes has a negative twist that improves base stacking. This observation is consistent with the low stability of TA base steps in B-form DNA.     

On the structure of poly d(A-S4T).

The helical twist of poly d(A-s4T) was determined from the periodicity of the cleavage patterns of the double stranded polydeoxynucleotide adsorbed on calcium phosphate and found to be 14 bp per turn. Both cleavage patterns and 31P NMR spectra indicate a mononucleotide structure rather than an alternating B DNA like poly d(A-T). The failure of nucleosome formation excludes a B type structure. The discrepancy of the mononucleotide structure found in 31P NMR spectra and the dinucleotide structure given by X ray fiber diffraction is explained by an alternating tilt of the planes of the base pairs (base roll) as a consequence of a strong propeller twist. The importance of interstrand stacking interactions of adjacent 4-thiothymidines for the helical stability is discussed.     

Structures of mismatched base pairs in DNA and their recognition by the Escherichia coli mismatch repair system.

The Escherichia coli mismatch repair system does not recognize and/or repair all mismatched base pairs with equal efficiency: whereas transition mismatches (G X T and A X C) are well repaired, the repair of some transversion mismatches (e.g. A X G or C X T) appears to depend on their position in heteroduplex DNA of phage lambda. Undecamers were synthesized and annealed to form heteroduplexes with a single base-pair mismatch in the centre and with the five base pairs flanking each side corresponding to either repaired or unrepaired heteroduplexes of lambda DNA. Nuclear magnetic resonance (n.m.r.) studies show that a G X A mismatch gives rise to an equilibrium between fully helical and a looped-out structure. In the unrepaired G X A mismatch duplex the latter predominates, while the helical structure is predominant in the case of repaired G X A and G X T mismatches. It appears that the E. coli mismatch repair enzymes recognize and repair intrahelical mismatched bases, but not the extrahelical bases in the looped-out structures.     

盐肤木茎的木材解剖学研究

对盐肤木茎进行木材解剖学研究,发现其导管分子穿孔板为单穿孔板。根据穿孔板在导管分子上的位置和数量,导管分子可分为4种类型,即只有一个单穿孔板,位于导管中央,且具有螺纹增厚;两端均为单穿孔板,且具螺纹增厚;两端均为单穿孔板,且具有孔纹增厚;具三个单穿孔板,一端有2个,另一端有1个,且具螺纹增厚。组成细胞中还有分隔木纤维、螺纹管胞。盐肤木为环孔材,射线为异形射线。     

Probing the helical tilt and dynamic properties of membrane-bound phospholamban in magnetically aligned bicelles using electron paramagnetic resonance spectroscopy

Wild-type phospholamban (WT-PLB), a Ca(2+)-ATPase (SERCA) regulator in the sarcoplasmic reticulum membrane, was studied using TOAC nitroxide spin labeling, magnetically aligned bicelles, and electron paramagnetic resonance (EPR) spectroscopy to ascertain structural and dynamic information. Different structural domains of PLB (transmembrane segment: positions 42 and 45, loop region: position 20, and cytoplasmic domain: position 10) were probed with rigid TOAC spin labels to extract the transmembrane helical tilt and structural dynamic information, which is crucial for understanding the regulatory function of PLB in modulating Ca(2+)-ATPase activity. Aligned experiments indicate that the transmembrane domain of wild-type PLB has a helical tilt of 13°±4° in DMPC/DHPC bicelles. TOAC spin labels placed on the WT-PLB transmembrane domain showed highly restricted motion with more than 100ns rotational correlation time (τ(c)); whereas the loop, and the cytoplasmic regions each consists of two distinct motional dynamics: one fast component in the sub-nanosecond scale and the other component is slower dynamics in the nanosecond range.