Carbohydrates and drug discovery--the role of computer simulation.

Recent advances in the molecular modelling of carbohydrates have brought this technique to a level comparable with that of protein and nucleic acid simulations. After a brief introduction to the techniques used in the computer simulation of carbohydrates and carbohydrate interactions, an overview of applications in the field of carbohydrate-related drug discovery is presented.

Effects of base sequence on the loop folding in DNA hairpins.

High-resolution NMR and UV-melting experiments have been used to study the hairpin formation of partly self-complementary DNA fragments in an attempt to derive rules that describe the folding in these molecules. Earlier experiments on the hexadecanucleotide d(ATCCTA-TTTT-TAGGAT) had indicated that within the loop of four thymidines a wobble T-T pair is formed (Blommers et al., 1987). In the present paper it is shown that if the first and the last thymines of the intervening sequence are replaced by complementary bases, sometimes base pairs can be formed. Thus for the intervening sequences -CTTG- and -TTTA- with the pyrimidine in the 5'-position and the purine in the 3'-position, a base pair is formed leading to a loop consisting of two residues. For the intervening sequences -GTTC- and -ATTT- with the purine in the 5'-position and the pyrimidine in the 3'-position, this turns out not to be the case. It was found that it made no difference when the four-membered sequence was closed by a G-C base pair or an A-T base pair. Replacement of the two central thymidine residues by the more bulky adenine residues limits the hairpin to a four-membered loop scheme. Very surprisingly, it was found from 2D NOE experiments that the T-A base pair, formed in the loop consisting of the -TTTA- sequence, is a Hoogsteen pair. It is argued that the pairing of the bases in this scheme may facilitate the formation of a loop of two residues, since the distance of the C1' atoms in this base pair is 8.6 A instead of 10.4 A found in the canonical Watson-Crick base pair. Combination of the data obtained for the series of DNA fragments studied shows that the results can be explained by a simple, earlier proposed, loop folding principle which assumes that the folding of the four-membered loop is dictated by the stacking of the double-helical stem of the hairpin.

Solution structure of the 3'-5' cyclic dinucleotide d(pApA). A combined NMR, UV melting, and molecular mechanics study.

The 3'-5' cyclic dinucleotide d(pApA) was studied by means of 1H and 13C NMR experiments, UV-melting experiments, and molecular mechanics calculations. The 1H and 13C NMR spectra were analyzed by means of 2-dimensional NMR experiments. J-Coupling analysis of the 1D and 2D 1H and 13C spectra was used to determine the conformation of the ring systems in the molecule. It appeared that at low temperature (283 K) the deoxyribose sugars adopt a N-type conformation. The geometry is best described by an intermediate between the 3(2)T and 3E forms. In addition, we were able to derive all other torsion angles in the phosphate backbone ring system, i.e., alpha +, beta t, gamma +, delta (= 89 degrees), epsilon t, and zeta +. When the molecule is subjected to an energy minimization procedure (using the program AMBER), the sugar ring system retains, practically speaking, the torsion angles found from the NMR experiments, while the torsion angles around the glycosidic bond adopt a value of 175 degrees in the minimum energy conformation. UV-melting experiments indicate that two molecules can form a dimer in which the adenine bases are intercalated. The feasibility of this structure is indicated by molecular mechanics calculations. At higher temperatures the dimer is converted into separate monomers. In the monomer form the sugars exhibit S-pucker 20% of the time. Concomitantly with the conversion of the N- to the S-conformation, the torsion angles alpha and gamma change.

Studies of the solution structure of the bleomycin-A2-zinc complex by means of two-dimensional NMR spectroscopy and distance geometry calculations.

Application of various two-dimensional NMR techniques (SECSY, COSY and NOESY) enabled the complete assignment of the 1H-NMR spectrum of the bleomycin-A2-zinc complex in H2O and D2O at pH 6.7. The spectra were interpreted at 277 K as well as at 300 K. Identification of the resonances permitted a vicinal coupling constant analysis which revealed that the conformation around the C alpha and C beta bond of the beta-aminoalanine and beta-hydroxyhistidine residues is fixed. From this finding it was concluded that both amino functions of the beta-aminoalanine fragment and the amide and imidazole groups of the beta-hydroxyhistidine moiety are involved in zinc coordination. Also, for the mannose carbamoyl group and the pyrimidine ring active participation in zinc coordination could be established. NOE data together with the six coordination sites proposed above were used as interpoint distance constraints in distance geometry calculations for the bleomycin-A2-zinc complex in H2O. Sets of ten structures, randomly chosen within the distance constraints, were calculated (with and without the zinc ion). The calculated structures were very similar but in case that the zinc ion was omitted some flexibility was observed, within the distance constraints, in the pyrimidine-aminoalanine region. Because of the great overall similarity between the structures, a reliable representation of the solution conformation of the bleomycin-zinc complex was reached. Surprisingly, no regular symmetry around the zinc ion was found to be present in the generated structures.

The dynamical structure of the RNA in alfalfa mosaic virus studied by 31P-nuclear magnetic resonance.

The structure of the viral RNA in alfalfa mosaic virus (AlMV) was investigated by means of 31P-nuclear magnetic resonance (NMR). It was found that the 31P-NMR line width of AlMV Top a particles is significantly smaller than that of the larger Bottom particles. At low temperatures, the totational correlation time of the 31P nuclei essentially equals the tumbling rate of the virus particle, indicating that the RNA is contained rigidly inside the virion. At more elevated temperatures, the NMR line width sharpens more than expected on the basis of viscosity changes and the RNA exhibits internal mobility. The occurrence of internal mobility is paralleled by an increased internal mobility of the N-terminal part of the coat protein, as could be observed by 1H-NMR spectroscopy. The influence of EDTA on the 31P-NMR line width appeared to be negligible, which is in agreement with the idea that AlMV does not 'swell' like several other RNA-containing plant viruses.

Circular dichroism and 500-MHz proton magnetic resonance studies of the interaction of Escherichia coli translational initiation factor 3 protein with the 16S ribosomal RNA 3' cloacin fragment.

The RNA helix destabilizing properties of Escherichia coli initiation factor 3 protein (IF3), and its affinity for an evolutionarily conserved sequence at the 3' end of 16S rRNA, led us to examine the details of the protein-nucleic acid interactions upon IF3 binding to the 49-nucleotide 3'-terminal cloacin DF13 fragment of 16S rRNA by studying the circular dichroism (CD) and proton magnetic resonance spectra of the RNA, the protein, and their complex. In a physiological tris(hydroxymethyl)aminomethane buffer, where the interaction is primarily nonionic and sequence specific, addition of IF3 decreases the RNA 268-nm CD peak hyperbolically by 19% to an end point of about one IF3 per RNA strand. The titration curve is best fit by an association constant of (1.80 +/- 0.05) X 10(7) M-1, within the range estimated by a nuclease mapping study of the same system [Wickstrom, E. (1983) Nucleic Acids Res. 11, 2035-2052]. In a low-salt phosphate buffer without Mg2+, where the interaction is primarily ionic and nonspecific, titration with IF3 decreases the peak CD sigmoidally by 35% to an end point of two IF3 per strand. The titration curve is best fit by an intrinsic association constant of (1.7 +/- 0.7) X 10(6) M-1 for each IF3 and a cooperativity constant of 33 +/- 6. In a physiological phosphate buffer lacking Mg2+, the dispersion of aromatic proton magnetic resonance peaks and upfield-shifted methyl proton resonances indicates a high degree of secondary and tertiary structure in the protein. In an equimolar mixture of IF3 and RNA cloacin fragment, several changes in identifiable IF3 and RNA resonances are observed.(ABSTRACT TRUNCATED AT 250 WORDS)

On loop folding in nucleic acid hairpin-type structures.

In a series of studies, combining NMR, optical melting and T-jump experiments, it was found that DNA hairpins display a maximum stability when the loop part of the molecule comprises four or five nucleotide residues. This is in contrast with the current notion based on RNA hairpin studies, from which it had been established that a maximum hairpin stability is obtained for six or seven residues in the loop. Here we present a structural model to rationalize these observations. This model is based on the notion that to a major extent base stacking interactions determine the stability of nucleic acid conformations. The model predicts that loop folding in RNA is characterized by an extension of the base stacking at the 5'-side of the double helix by five or six bases; the remaining gap can then easily be closed by two nucleotides. Conversely, loop folding in DNA is characterized by extending base stacking at the 3'-side of the double helical stem by two or three residues; again bridging of the remaining gap can then be achieved by one or two nucleotides. As an example of loop folding in RNA the anticodon loop of yeast tRNAPhe is discussed. For the DNA hairpin formed by d(ATCCTAT4TAGGAT) it is shown that the loop structure obtained from molecular mechanics calculations obeys the above worded loop folding principles.

Structure elucidation of two isomeric steroids: photolytical and thermal reaction products from norethisterone studied by two-dimensional nuclear magnetic resonance.

Two-dimensional nuclear magnetic resonance was used for the structure elucidation of two isomeric photoproducts of norethisterone, a commonly used progestogen in oral contraceptives. The predominant one of the two isolated products derived from photochemical decomposition of norethisterone upon irradiation with UV-B light (280-320 nm) was 5 alpha, 17 beta-dihydroxy-19-nor-17 alpha-pregn-20-yn-3-one. The minor photoproduct appeared to be the analogous 5 beta-isomer, i.e. 5 beta, 17 beta-dihydroxy-19-nor-17 alpha-pregn-20-yn-3-one. The latter compound was also obtained from the thermal reduction of norethisterone-4 beta, 5 beta-epoxide using aluminium amalgam in isopropanol. Two-dimensional NMR appeared to be superior to mass and IR spectrometry in identifying the isomers.

The photochemical decomposition of the progestogenic 19-norsteroid, norethisterone, in aqueous medium.

Norethisterone, a contraceptive 19-norsteroid, was decomposed in aqueous medium (pH 7.4) by UV-B radiation (280-320 nm). This 4-en-3-oxo-19-norsteroid was not prone to the skeletal rearrangement reactions usually observed in steroids possessing a C10-methyl group. Under the reaction conditions applied, products were formed by addition of molecules, such as solvent molecules or a second steroid molecule, and by reduction of the double bond. The prevalence of addition type reactions may have consequences for the application of norethisterone-like steroids in subdermal contraceptive devices.

Hairpin formation in synthetic oligonucleotides.

The structure and dynamics of the homologous series of the (partly) self-complementary DNA fragments, d(ATCCTATnTAGGAT) n = 0-7, were investigated in a combined NMR, T-jump, and optical melting study. It is shown that all compounds in the series may adopt hairpin like conformations, even for n less than 3, although for these smaller n values this only occurs in significant amounts at relatively low concentrations (approximately 10 microM). The enthalpy change accompanying the hairpin-coil melting transition turns out to depend on the number of intervening thymidines, n. It is shown that this does not mean that the enthalpy of loop closure is significantly different from zero, but that loop formation stabilizes the base pair closing the loop. The results indicate that for DNA the optimal loop consists of four or five residues. The observation that hairpins are formed for n less than 3 and that the stability of DNA hairpins is at its maximum for loop lengths of four to five residues is at variance with earlier findings for RNA. In the latter case the optimal loop size consists of six to seven residues, whereas for less than three intervening residues only, dimer, and no hairpin formation, was observed [17, 20]. A direct comparison with RNA behaviour was made by studying r(AUCCUAUT4UAGGAU), T = ribothymidine. In contrast to its DNA analogue, d(ATCCTAT4TAGGAT), the ribo-fragment forms a dimer as well as a hairpin at low (10 microM) concentrations. With the thermodynamic melting parameters deduced from the present experiments the differences between DNA and RNA melting behaviour can be explained.

cis-diamminedichloroplatinum(II) induced distortion of a single and double stranded deoxydecanucleosidenonaphosphate studied by nuclear magnetic resonance.

The structural distortion of a single- and a double-stranded decadeoxynucleotide upon binding of cis-PtCl2(NH3)2 was studied by 1H-NMR. After selective platination of d(T-C-T-C-G-G-T-C-T-C) (I) at the central d(-GpG-) site (resulting in I-Pt), several non-exchangeable base protons as well as H1', H2', H2" and H3' protons could be assigned by means of conventional NMR double-resonance techniques. Addition of the complementary decamer strand to I and I-Pt yielded the double-stranded III and III-Pt, respectively. All non-exchangeable base, H1', and most of the H2' and H2" protons in the two double stranded compounds could be assigned using 2D-chemical shift correlation (COSY) and nuclear Overhauser enhancement (NOESY) techniques. The double stranded compound III appears to adopt a B-DNA like structure. Comparison of NOEs and proton-proton coupling constants in the d(-GpG-).cisPt part in I-Pt and III-Pt reveals that their structure displays large similarity. Significant chemical shift changes (i.e. larger than 0.1 ppm) between III and III-Pt are restricted to the central four base pairs. It follows that the outer three base pairs, located on either side of the central four base pairs in III-Pt are likely to adopt a regular B-DNA type helix. The observed large upfield and downfield chemical shifts in the d(-CpGpG-) part of III with respect to III-Pt can be rationalized by describing the distortion of the double helix as a kink. A discussion of the observed physical effects upon platination of a double-stranded oligonucleotide is presented.

Imino-proton resonances of yeast tRNAPhe studied by two-dimensional nuclear Overhauser enhancement spectroscopy.

Application of two-dimensional nuclear Overhauser enhancement (NOE) spectroscopy to yeast tRNAPhe in H2O solution demonstrates that all imino-proton resonances, related to the secondary structure, and nearly all imino proton resonances, originating from the tertiary structure, can be assigned efficiently by this method. The results corroborate the assignments of the imino-proton resonances of this tRNA as established previously by one-dimensional NOE experiments (only the assignment of base pairs G1 X C72 and C2 X G71 should be reversed). The advantages of two-dimensional NOE spectroscopy over one-dimensional NOE spectroscopy for the assignments of imino-proton resonances and the structure elucidation of tRNA are illustrated and discussed. Furthermore, the use of non-exchangeable proton resonances as probes of the molecular structure is explored.

Carbon-13 NMR in conformational analysis of nucleic acid fragments. 4. The torsion angle distribution about the C3'-O3' bond in DNA constituents.

Carbon-13 and proton NMR spectra of a series of oligodeoxynucleotides (d(CT), d(CC), d(TA), d(AT), d(CG), d(GC), d(AG), d(AAA), d(TATA) and d(GGTAAT] were measured at various temperatures. The three coupling constants that are related to the magnitude of backbone angle epsilon (J(C4'-P), J(C2'-P) and J(H3'-P] are analyzed in terms of a three-state equilibrium about this bond. Two epsilon (trans) angles occur, which differ in magnitude depending on the conformation (N or S) of the adjoining deoxyribose ring. The S-type deoxyribose ring is associated with a smaller epsilon (trans) angle: epsilon (t,S) = 192 degrees. The N-type deoxyribose ring is associated with a larger epsilon (trans) angle epsilon (t,N) = 212 degrees. The third rotamer participating in the conformational equilibrium, is a gauche(-) (epsilon (-] conformer and occurs exclusively in combination with the S-type sugar ring (epsilon (-,S) = 266 degrees). Within the limits of experimental error, the magnitude of these three angles appears to be independent of the particular base sequence, except in the case of d(CG) where a slightly larger epsilon (t,S) angle (197 degrees) is indicated. A simple equation is proposed which may be used to calculate the population of epsilon (t,S) conformer in cases where only J(H3'-P) is known.

Complete assignment of the 500 MHz 1H-NMR spectra of bleomycin A2 in H2O and D2O solution by means of two-dimensional NMR spectroscopy.

1H-NMR spectra of bleomycin A2 recorded at 500 MHz in D2O and H2O at 24 degrees C and 3 degrees C were investigated. Resonances of the individual spin systems were identified by using two-dimensional correlation spectroscopy (COSY), two-dimensional spin echo correlated spectroscopy (SECSY) and by the application of two-dimensional Nuclear Overhauser Effect spectroscopy (NOESY). Employment of these techniques allowed the assignment of 113 exchangeable and 59 non-exchangeable protons in the 1H NMR spectrum of bleomycin A2. By means of 2D NOE spectroscopy also interresidual connectivities could be observed. Comparison of the NOESY spectra at 3 degrees C and 24 degrees C suggest that at low temperatures the central party of the bleomycin A2 molecule tends to adopt an extended conformation.

Discrimination between A-type and B-type conformations of double helical nucleic acid fragments in solution by means of two-dimensional nuclear Overhauser experiments.

The solution structure of two double helical nucleic acid fragments, viz, r(CGCGCG) and d(CGCGCG), was probed by means of two-dimensional nuclear Overhauser effect spectroscopy. The two compounds were selected as models for the A-type and B-type double helical conformations, respectively, and it is shown that for each of the two model compounds the intensities of the NOE cross peaks between base- and H2' (deoxy)ribose proteins are qualitatively in correspondence with the relative NOE intensities expected on basis of the supposed duplex conformations. Thus our results indicate that NOE-data can be used to differentiate between A-and B-type double helical conformations in solution. Coupling constant data show that, except for G(6), all ribose rings in r(CGCGCG) adopt pure N (C3'-endo) conformations thereby manifesting that this molecule takes up a regular A-type double helical conformation in solution. In contrast, the deoxyribose rings in d(CGCGCG) retain conformational freedom in the duplex state, albeit that the N/S-equilibrium is biased towards the S (C2'-endo) sugar conformation. This finding indicates that in solution the B-DNA backbone is highly dynamic.

Sequence-dependent structural variation in single-helical DNA. Proton NMR studies of d(T-A-T-A) and d(A-T-A-T) in aqueous solution.

The two deoxyribotetranucleoside triphosphates d(T-A-T-A) and d(A-T-A-T) were investigated in aqueous solution by one- and two-dimensional proton NMR at 300 and 500 MHz. It is demonstrated that both compounds occur predominantly in the single-helical form. Accurate coupling constants are obtained by computer simulation of several 500-MHz spectra. The data are interpreted in terms of N and S pseudorotational ranges. The geometry of the major S-type conformers displays a clear sequence dependence, as expressed by variation of the endocyclic backbone angle delta (C5'-C4'-C3'-O3'). A simple sum rule is proposed to predict delta variation in single-helical DNA fragments. Comparisons are made with other sequence-dependent geometries as observed in a double-helical B-DNA fragment in the crystalline state. Furthermore, one- and two-dimensional nuclear Overhauser effect (NOE) spectroscopy was carried out on d(T-A-T-A). An inventory is made of the observed intra- and inter-residue NOEs. The NOE data confirm the presence of a highly stacked single-helical conformation of d(T-A-T-A) in solution. No indications are found for the formation of a bulge-out structure as observed for analogous alternating purine-pyrimidine oligoribonucleotides.

Carbon-13 NMR in conformational analysis of nucleic acid fragments. 3. The magnitude of torsional angle epsilon in d(TpA) from CCOP and HCOP NMR coupling constants.

Carbon-13 NMR spectra of the deoxyribonucleotide d(TpA), 3',5'-cyclic AMP and 3',5'-cyclic dAMP were measured. It is shown that the different substitution of C2' in deoxyribonucleotides versus ribonucleotides does not affect the vicinal C2'-C3'-O3'-P coupling to a measurable extent. Therefore, the same set of Karplus parameters may be used for the C2'-C3-O3'-P couplings in ribonucleotides and in deoxyribonucleotides. Vicinal carbon-phosphorus and proton-phosphorus coupling constants are used to calculate the magnitude of the torsion angle epsilon (C4'-C3'-O3'-P), which amounts to 195(0) in the trans conformer and to 261(0) in the gauche(-) conformer.

Carbon-13 NMR in conformational analysis of nucleic acid fragments. 2. A reparametrization of the Karplus equation for vicinal NMR coupling constants in CCOP and HCOP fragments.

13C-31P coupling constants of 10 oligoribonucleoside phosphates, measured at a number of temperatures, are presented. The combination of these data with 1H-31P couplings of the same compounds leads to the derivation of two new and mutually consistent sets of Karplus parameters: J(CCOP) = 6.9cos2 phi--3.4cos phi + 0.7 J(HCOP) = 15.3cos2 phi--6.1cos phi + 1.6 At the same time new values for the base sequence dependent magnitude of the trans conformer of the backbone angle epsilon (C4'-C3'-O3'-P) are calculated. The present results show that the magnitude of epsilon(t) in right-handed ribo helices is confined to the range 214 degrees-226 degrees (average 219 degrees), which is in much better agreement with single crystal X-ray studies (average 218 degrees) than were previous deductions from NMR spectroscopic results (average 208 degrees).

Conformational analysis of the single-helical DNA fragment d(T-A-A-T) in aqueous solution. The combined use of NMR proton chemical shifts and coupling constants obtained at 300 MHz and 500 MHz.

Proton NMR studies at 300 MHz and 500 MHz were carried out on the tetranucleoside trisphosphate d(T-A-A-T). The thermodynamics of the three stacking interactions, derived from chemical shift versus temperature profiles, were used to extrapolate the observed coupling constants, measured at a range of temperatures, to values appropriate to the fully stacked forms of the molecule. The data were interpreted in terms of N and S pseudorotational ranges [ Altona , C. and Sundaralingham , M. (1972) J. Am. Chem. Soc. 94, 8205-8212]. It is shown that the stacked state of the molecule cannot be described by one conformer, but consists of one major structure (60%) in which all sugar rings have S-type geometry and another structure (30%) in which residue dT(4) has an N-type sugar. The remainder of the stacked states consists of one or more conformers with two or three sugar residues in the N-type pseudorotational range. Detailed geometrical models are proposed for the major stacked conformers encountered in aqueous solution.