Refinement of the solution structure of the DNA decamer 5'd(CTGGATCCAG)2: combined use of nuclear magnetic resonance and restrained molecular dynamics

The solution structure of the self-complementary DNA decamer 5'd(CTGGATCCAG)2 comprising the specific target site for the restriction endonuclease BamH1 is investigated by using nuclear magnetic resonance sectroscopy and restrained molecular dynamics. With the exception of the H5'/H5" sugar proton resonances, all the nonexchangeable proton resonances are assigned sequentially by using pure-phase absorption two-dimensional nuclear Overhauser enhancement spectroscopy. From the time dependence of the nuclear Overhauser effects a set of 160 approximate interproton distances is determined and used as the basis of a structure refinement employing restrained molecular dynamics in which the interproton distances are incorporated into the total energy function of the system in the form of an effective potential term. Two restrained dynamics simulations are carried out, starting from classical B- and A-DNA [atomic root mean square (rms) difference 5.7 A]. In both cases convergence is achieved to very similar B-type structures with an atomic rms difference of 0.9 A which is comparable to the rms fluctuations of the atoms about their average positions. In addition, the rms difference between the experimental and calculated values of the interproton distances for both average restrained dynamics structures is approximately 0.3 A. These results suggest that the converged restrained molecular dynamics structures represent reasonable approximations of the solution structure. The average restrained dynamics structures exhibit clear sequence-dependent variations of torsion angles and helical parameters. In addition, the structures exhibit a small bend of around 10-20 degrees at the second (TpG) and eighth (CpA) base pair steps. This can be attributed to the positive base roll angles and large base pair slide values at the two Pyr-Pur steps. The central core of the decamer comprising the six-base recognition site for BamH1 (GGATCC), however, is straight.     

An investigation into the solution structures of two self-complementary DNA oligomers, 5''-d(C-G-T-A-C-G) and 5''-d(A-C-G-C-G-C-G-T), by means of nuclear-Overhauser-enhancement measurements.

A 500 MHz 1H-n.m.r. study on two self-complementary alternating pyrimidine-purine oligodeoxyribonucleotides, 5'-d(C-G-T-A-C-G) and 5'-d(A-C-G-C-G-C-G-T), is presented. By using the proton-proton nuclear Overhauser effect virtually complete assignments are obtained and a large number of interproton distances [113 in the case of 5'-d(C-G-T-A-C-G) and 79 in the case of 5'-d(A-C-G-C-G-C-G-T)], both intra- and inter-nucleotide, are determined. The interproton-distance data are consistent with an overall right-handed B-DNA-type structure for both oligonucleotides, in agreement with their B-type c.d. spectra. However, whereas 5'-d(C-G-T-A-C-G) adopts a conventional B-type structure with a mononucleotide repeating unit, the interproton-distance data provide evidence that 5'-d(A-C-G-C-G-C-G-T) has a dinucleotide repeating unit consisting of alternation in glycosidic bond and sugar pucker conformations.     

Refinement of the solution structure of the ribonucleotide 5'r(GCAUGC)2: combined use of nuclear magnetic resonance and restrained molecular dynamics

The solution structure of the self-complementary hexamer 5'r(GCAUGC)2 is investigated by means of nuclear magnetic resonance spectroscopy and restrained molecular dynamics. The proton resonances are assigned in a sequential manner, and a set of 110 approximate interproton distance restraints are derived from the two-dimensional nuclear Overhauser enhancement spectra. These distances are used as the basis of a structure refinement by restrained molecular dynamics in which the experimental restraints are incorporated into the total energy function of the system in the form of effective potentials. Eight restrained molecular dynamics simulations are carried out, four starting from a structure with regular A-type geometry and four from one with regular B-type geometry. The atomic root mean square (rms) difference between the initial structures is 3.2 A. In the case of all eight simulations, convergence is achieved both globally and locally to a set of very similar A-type structures with an average atomic rms difference between them of 0.8 +/- 0.2 A. Further, the atomic rms differences between the restrained dynamics structures obtained by starting out from the same initial structures but with different random number seeds for the assignment of the initial velocities are the same as those between the restrained dynamics structures starting out from the two different initial structures. These results suggest that the restrained dynamics structures represent good approximations of the solution structure. The converged structures exhibit clear sequence-dependent variation in some of the helical parameters, in particular helix twist, roll, slide, and propellor twist.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structure refinement of oligonucleotides by molecular dynamics with nuclear Overhauser effect interproton distance restraints: application to 5' d(C-G-T-A-C-G)2

The solution structure of the self-complementary DNA hexamer 5' d(C-G-T-A-C-G)2 is refined by restrained molecular dynamics in which 192 interproton distances, determined from pre-steady-state nuclear Overhauser enhancement measurements, are incorporated into the total energy of the system in the form of effective potentials. First the method is tested by applying an idealized set of distance restraints taken from classical B-DNA to a simulation starting off from A-DNA and vice versa. It is shown that in both cases the expected transition between A- and B-DNA occurs. Second, a set of restrained molecular dynamics calculations is carried out starting from both A- and B-DNA with the experimental interproton distances for 5' d(C-G-T-A-C-G)2 as restraints. Convergence to the same B-type structure is achieved with the interproton distances equal to the measured values within experimental error. The root-mean-square atomic difference between the two average restrained dynamics structures (less than 1 A) is approximately the same as the root-mean-square fluctuations of the atoms.     

NMR-based structure of anticancer drug mitoxantrone stacked with terminal base pair of DNA hexamer sequence d-(ATCGAT)2

Mitoxantrone is a promising antitumor drug having considerably reduced cardiotoxicity as compared to anthracyclines. Its binding to deoxyhexanucleotides sequence d-(ATCGAT)₂ has been studied by proton and phosphorous-31 nuclear magnetic resonance spectroscopy. The stoichiometry reveals that 1:1 and 2:1 mitoxantrone-d(ATCGAT)₂ complexes are formed in solution. Significant upfield shifts in 6H/7H, 2H/3H, 11NH, and 12NH protons (～.5?ppm) of mitoxantrone and T6NH imino protons (～.3?ppm) are observed. The phosphorous resonances do not shift significantly indicating that the base pairs do not open at any nucleotide step along the sequence of hexamer. Several inter-molecular Nuclear Overhauser Enhancement connectivities between mitoxantrone and hexanucleotide protons indicate that mitoxantrone chromophore stacks with terminal A1-T6 base pair and side chains involving 12CH₂, 12NH, and 14OH protons are in close proximity of A1, T2, A5, and T6 bases. Absorption and emission spectra show red shift in wavelength maxima, which is characteristic of stacking interaction. At higher mitoxantrone to nucleic acid ratios, electrostatic interactions are dominant. The 2:1 drug/DNA stoichiometric structure obtained by restrained Molecular Dynamics simulations shows considerable distortions in backbone torsional angles and helicoidal parameters although structural fluctuations in 25?ps analysis of trajectory are found to be negligible. Mitoxantrone binds as a monomer at either or both ends of hexamer externally with side chains interacting specifically with DNA. The findings are relevant to the understanding of pharmacological action of drug.     

Three-dimensional structure of a DNA hairpin in solution: two-dimensional NMR studies and distance geometry calculations on d(CGCGTTTTCGCG)

The three-dimensional structure of d(CGCGTTTTCGCG) in solution has been determined from proton NMR data by using distance geometry methods. The rate of dipolar cross-relaxation between protons close together in space is used to calculate distances between proton pairs within 5 A of each other; these distances are used as input to a distance geometry algorithm that embeds this distance matrix in three-dimensional space. The resulting refined structures that best agree with the input distances are all very similar to each other and show that the DNA sequence forms a hairpin in solution; the bases of the loop region are stacked, and the stem region forms a right-handed helix. The advantages and limitations of the technique, as well as the computer requirements of the algorithm, are discussed.     

The structure of the double-stranded RNA pentamer 5'(CACAG) . 5'(CUGUG) determined by nuclear Overhauser enhancement measurements: interproton distance determination and structure refinement on the basis of X-ray coordinates

The structure in solution of the duplex RNA pentamer 5'(CACAG) . 5'(CUGUG), comprising the stem of the T psi C loop of yeast tRNAPhe, has been investigated by means of one- and two-dimensional nuclear Overhauser enhancement measurements. All non-exchangeable base and sugar proton resonances with the exception of the H5'/H5" sugar resonances are assigned in a sequential manner. From the relative intensities of the cross-peaks obtained in the pure-phase absorption two-dimensional nuclear Overhauser enhancement spectra at several mixing times, it is deduced that the RNA pentamer adopts an A-type conformation in solution. Cross-relaxation rates and interproton distances are determined from the time dependence of the nuclear Overhauser effects, principally by one-dimensional measurements. The structure of the RNA pentamer is then refined by restrained least-squares minimization on the basis of both distance and planarity restraints using fibre diffraction data as an initial model. The refined structure of the RNA pentamer is of the A type but exhibits local structural variations in glycosidic bond and backbone torsion angles as well as in propeller twist, base roll and base tilt angles.     

Hybrid DNA i-motif: Aminoethylprolyl-PNA (pC5) enhance the stability of DNA (dC5) i-motif structure

This report describes the synthesis of C-rich sequence, cytosine pentamer, of aep-PNA and its biophysical studies for the formation of hybrid DNA:aep-PNAi-motif structure with DNA cytosine pentamer (dC₅) under acidic pH conditions. Herein, the CD/UV/NMR/ESI-Mass studies strongly support the formation of stable hybrid DNA i-motif structure with aep-PNA even near acidic conditions. Hence aep-PNA C-rich sequence cytosine could be considered as potential DNA i-motif stabilizing agents in vivo conditions.     

Refinement of the solution structure of the RNA-DNA hybrid 5'-[r(GCA)d(TGC)]2. Combined use of nuclear magnetic resonance and restrained molecular dynamics

The solution conformation of the self-complementary RNA-DNA hybrid hexamer 5'-[r(GCA)d(TGC)]2 is investigated by NMR spectroscopy and restrained molecular dynamics. The 1H-NMR spectrum is assigned in a sequential manner using two-dimensional homonuclear Hartmann-Hahn and nuclear Overhauser enhancement spectroscopy. From the latter a set of 178 approximate interproton distance restraints are determined and used as the basis of a structure refinement by restrained molecular dynamics. Eight independent calculations are carried out, four from a classical A-type geometry and four from a classical B-type one. Convergence is achieved to very similar A-type structures with an average atomic root mean square difference between them of 1.0 +/- 0.2 A. The converged structures exhibit variations in helical parameters similar to those found previously for the analogue RNA hexamer 5'-r(GCAUGC)2 [(1988) Biochemistry 27, 1735-1743].     

Crystallization and preliminary X-ray analysis of anti-cancer agent 3-(9-acridinylamino)-5-(hydroxymethyl)aniline complexed with the DNA hexamer d(CGTACG)2

3-(9-Acridinylamino)-5-(hydroxymethyl)aniline (AHMA) is an anti-cancer agent with significant efficacy against murine leukemia and solid tumors. As a DNA topoisomerase inhibitor, AHMA is proposed to form a ternary complex with DNA and topoisomerase and bind to DNA in an intercalative manner. In order to understand the interactions between AHMA and DNA and study the structure-function relationship of amsacrine analogue, the AHMA-d(CGTACG)(2) complex was crystallized using the sitting-drop vapor-diffusion method. The native crystals diffract to 2.9-A resolution and belong to space group P3(1)21 or P3(2)21 with unit-cell parameters a=b=57.52, c=122.17 A when analyzed using Cu Kalpha radiation. Patterson map indicates that in the crystal, the directions of the DNA base stacking are nearly perpendicular to the c-axis of the crystal unit cell.     

Role of B13 Glu in insulin assembly. The hexamer structure of recombinant mutant (B13 Glu-->Gln) insulin.

The assembly of the insulin hexamer brings the six B13 glutamate side-chains at the centre into close proximity. Their mutual repulsion is unfavourable and zinc co-ordination to B10 histidine is necessary to stabilize the well known zinc-containing hexamers. Since B13 is always a carboxylic acid in all known sequences of hexamer forming insulins, it is likely to be important in the hormone's biology. The mutation of B13 Glu-->Gln leads to a stable zinc-free hexamer with somewhat reduced potency. The structures of the zinc-free B13 Gln hexamer and the 2Zn B13 insulin hexamer have been determined by X-ray analysis and refined with 2.5 A and 2.0 A diffraction data, respectively. Comparisons show that in 2Zn B13 Gln insulin, the hexamer structure (T6) is very like that of the native hormone. On the other hand, the zinc-free hexamer assumes a quaternary structure (T3/R3) seen in the native 4Zn insulin hexamer, and normally associated only with high chloride ion concentrations in the medium. The crystal structures show the B13 Gln side-chains only contact water in contrast to the B13 glutamate in 2Zn insulin. The solvation of the B13 Gln may be associated with this residue favouring helix at B1 to B8. The low potency of the B13 Gln insulin also suggests the residue influences the hormone's conformation.     

Conformation of the DNA undecamer 5'd(A-A-G-T-G-T-G-A-T-A-T) bound to the single-stranded DNA binding protein of Escherichia coli. A time-dependent transferred nuclear Overhauser enhancement study

A time-dependent transferred nuclear Overhauser enhancement study of the conformation of the single-stranded DNA 11mer 5'd(A-A-G-T-G-T-G-A-T-A-T) bound to the single-stranded DNA binding protein of Escherichia coli (SSB) is presented. It is shown that the conformation of the bound 11mer is that of a right-handed B-type helix similar to that of the free 11mer. The observation of internucleotide transferred nuclear Overhauser enhancements for every base step excludes the possibility of intercalation by aromatic protein residues. In addition, it is shown that the effective correlation time of the bases (80 ns) corresponds to that of a complex of molecular weight approximately 170,000, containing two SSB tetramers. The sugars, on the other hand, exhibit a shorter effective correlation time (40 ns), indicating the presence of internal motion. This suggests that the bases are anchored to the protein surface, possibly by hydrophobic interactions, whereas the sugar-phosphate groups are directed outwards towards the solvent.     

Determination of internal dynamics of deoxyriboses in the DNA hexamer d(CGTACG)2 by 1H NMR

The conformations and internal dynamics of the deoxyriboses of d(CGTACG)₂ have been determined by NMR measurements at 15°C. The conformations of the sugars were determined using coupling constants and time-dependent NOE measurements. The J-splitting patterns of the H1, H2 and H2 resonances show that the sugars exist as mixtures of conformations near C2 endo (south) and C3 endo (north). The population of the south conformation was larger for the purines than for the pyrimidines. The overall tumbling time of the molecule in ²H₂O was determined from measurements of the cross relaxation rate constant for the H6-H5 vectors of the two cytosine residues. Order parameters were determined for the H1-H2, H2-H2 and H2-H3 vectors from measurements of cross relaxation rate constants, making use of multi-spin analysis of the NOE build up rates. These order parameters are weakly dependent of the base sequence, and except for the terminal Cyt 1 residue, the H2-H2 and H2-H3 vectors are near unity, indicating the absence of rapid pseudorotation on the nanosecond time scale. However, the order parameter for the H1-H2 vector is significantly smaller than expected for rapid pseudorotation indicating the presence of other motions of the sugars. This motion must be about an effective axis parallel to the H2-H vector, and to occur with an angular fluctuation of about 30°.The results show that to obtain highly refined structures for nucleic acids by NMR the effects of spin diffusion and motional averaging cannot be ignored.Some of this work was presented as a poster at the 30^th Experimental NMR Conference at Asilomar, California 1989     

Population-genetic structure of Chuvashia (from data on eight DNA loci in the nuclear genome)

Population-genetic study of indigenous populations representing three ethnic Chuvash group: highland (Cheboksarsk and Morgaush district), lowland (Kanash district) and mid-lowland (Marposad district). Eight polymorphic DNA loci of the nuclear genome (VNTR/PAH, STR/PAH, VNTR/ApoB, VNTR/DAT1, APF, VNTR/eNOS, IVS6aGATT, and KM.19/PstI) were examined in the population of each district. For each of the four population, we estimated the allele and genotype frequency distributions at each polymorphic system, heterozygosities HS and between-population differences FST. In the combined Chuvash sample, HS = 0.464 and FST = 0.006. Loci VNTR(DAT) and VNTR(ApoB) showed highest between-population differentiation (0.009 < or = FST < or = 0.012), and loci IVS6aGATT, APF, VNTR/eNOS, and D7S23 (KM.19), lowest differentiation (0.001 < or = FST < or = 0.003). Analysis of genetic distances revealed somewhat higher genetic similarity between the Cheboksarsk and Morgaush populations belonging to the highland Chuvash group, whereas the highland Chuvash population from the Marposad district, which belong to the mid-lowland group, was more distant from the former populations.     

The crystal structure of the complex between a disaccharide anthracycline and the DNA hexamer d(CGATCG) reveals two different binding sites involving two DNA duplexes

The crystal structure of the complex formed between the anthracycline antibiotic 3'-deamino-3'- hydroxy-4'-(O-L-daunosaminyl)-4-demethoxydoxo rubicin (MEN 10755), an active disaccharide analogue of doxorubicin, and the DNA hexamer d(CGATCG) has been solved to a resolution of 2.1 A. MEN 10755 exhibits a broad spectrum of antitumor activities, comparable with that of the parent compound, but there are differences in the mechanism of action as it is active in doxorubicin-resistant tumors and is more effective in stimulating topoisomerase DNA cleavage. The structure is similar to previously crystallised anthracycline- DNA complexes. However, two different binding sites arise from drug intercalation so that the two halves of the self-complementary duplex are no longer equivalent. In one site both sugar rings lie in the minor groove. In the other site the second sugar protrudes out from the DNA helix and is linked, through hydrogen bonds, to guanine of a symmetry-related DNA molecule. This is the first structure of an anthracycline-DNA complex where an interaction of the drug with a second DNA helix is observed. We discuss the present findings with respect to the relevance of the amino group for DNA binding and to the potential role played by the second sugar in the interactions with topoisomerases or other cellular targets.     

Molecular basis of the dual functions of 2B4 (CD244)

2B4 belongs to the CD2 family of molecules and is expressed on all NK, gammadelta, and memory CD8(+) (alphabeta) T cells. The murine NK receptor 2B4 exhibits both inhibitory and activating functions, whereas human 2B4 has been reported to be an activating molecule. How murine 2B4 can act both as an activating and inhibitory receptor and what distinguishes its function from human 2B4 have remained largely unknown. We use here a model system that allows the study of human and murine 2B4 under identical and controlled conditions. These studies reveal that both human and mouse 2B4 can activate or inhibit NK cells. We show here that the level of 2B4 expression and the degree of 2B4 cross-linking play a significant role in the regulation of signaling lymphocyte activation molecule-associated protein-mediated activation by 2B4. A high level of 2B4 expression, heavy cross-linking, and relative paucity of signaling lymphocyte activation molecule-associated protein promote inhibitory function. Our studies demonstrate how a single receptor can have opposing function depending on the degree of receptor expression, extent of its ligation, and the relative abundance of certain adaptor molecules. Because the levels of 2B4 and CD48 are dynamically regulated, these findings have implications for the regulation of NK cell function.     

Conformational analysis of the trinucleoside diphosphate 3'd(A2'-5'A2'-5'A). An NMR and CD study.

A 500 MHz and 300 MHz NMR study of the trinucleoside diphosphate 3'd(A2'-5'A2'-5'A) is presented. In addition, circular dichroism is used to study base stacking in the title compound. The complete 1H-NMR spectral assignment of the sugar ring proton signals is given. Information about the sugar ring (N- or S-type conformation) and about the backbone geometry along C4'-C5' and C5'-O5' bonds is obtained from the NMR coupling constants. It is shown that the trimer mainly occurs in the N-N-N stacked state at low temperatures; the presence of a minor amount of N-N-S conformational sequence is indicated.     

The Glu(B13) carboxylates of the insulin hexamer form a cage for Cd2+ and Ca2+ ions

Substitution of Cd2+ for Zn2+ yields a hexameric insulin species containing 3 mol of metal ion per hexamer. The Cd2+ binding loci consist of the two His(B10) sites and a new site involving the Glu(B13) residues located at the center of the hexamer [Sudmeier, J. L., Bell, S. J., Storm, M. C., & Dunn, M. F. (1981) Science (Washington, D.C.) 212, 560-562]. Substitution of Co2+ or Co3+ for Zn2+ gives hexamers containing 2 mol of metal per hexamer. Insulin solutions to which both Cd2+ and Co2+ have been added in a ratio of 6:2:1 [In]:[Co2+]:[Cd2+] followed by oxidation to the exchange-inert Co3+ state yield stable hybrid species containing both Co3+ and Cd2+ with a composition of (In)6(Co3+)2Cd2+. The kinetics of the reaction of 2,2',2"-terpyridine (terpy) with the exchange-labile (In)6(Cd2+)2 and (In)6(Co2+)2 derivatives are biphasic and involve the rapid formation of an intermediate with coordination of one terpy molecule to each protein-bound metal ion; then, in a rate-limiting step the terpy-coordinated metal ion dissociates from the protein, and a second molecule of terpy binds to the metal ion to form a bis complex. Reaction of the exchange-inert Co3+ ions of (In)6(Co3+)2 with terpy is a slow apparent first-order process (t1/2 = 13.1 h). In contrast to the kinetic behavior of (In)6(Co2+)2 and (In)6(Cd2+)2, the Cd2+ ions bound to the hybrid (In)6(Co3+)2Cd2+ react quite slowly with terpy (t1/2 = 1 h at pH 8.0).(ABSTRACT TRUNCATED AT 250 WORDS)