On the use of ultraviolet resonance Raman intensities to elaborate molecular force fields: application to nucleic acid bases and aromatic amino acid residues models.

Normal modes analyses for different molecules with biological interest have been performed and checked via the calculation of resonance Raman intensities. For this purpose, molecular orbital calculations were used to determine bond order changes in the lowest-lying electronic transitions. These bond order changes were used to calculate resonance Raman intensities in order to obtain correct vibrational assignments and reliable force fields.

Fractal-like patterns in DNA films, B form at 0% relative humidity, and antiheteronomous DNA: an IR study.

This report details the observation of (a) the intact double helix of DNA at 0% relative humidity in poly(dA).poly(dT) and poly(dA-dT).poly(dA-dT) films, (b) the fractal-like growth of DNA crystals in films of poly(dA).poly(dT), and (c) poly(dA).poly(dT) with adenines in the B form and thymines in the A form. Observation (a) is based on the behavior of the middle ir signature of double-helical base stacking, the 1714 cm-1 peak, at low water activity. This observation is modeled as a trapping of water in the polycrystalline part of these films. We interpret the glycosidic region of the middle ir spectra of the polycrystalline films of poly(dA).poly(dT) at 0% relative humidity to indicate that at least part of the adenine strand in is the B conformation and the thymine strand is in the A form, whereas previous assignments of poly(dA).poly(dT) films at 75% relative humidity indicated the opposite strand conformation assignment [E. Taillandier et al. (1987) Biochemistry, Vol. 26 p. 3361].

Details of the acyl-enzyme intermediate and the oxyanion hole in serine protease catalysis.

Raman, absorbance, and kinetic measurements were used to determine how the serine protease active site feature known as the oxyanion hole interacts with an acyl-enzyme intermediate. The substrate, p-(dimethylamino)benzoylimidazolide (DAB-Im), was synthesized and used to prepare DAB-acyl-enzymes of wild-type (WT) and N155G subtilisin-BPN' (the N155G mutant lacks a fully functioning oxyanion hole), alpha-chymotrypsin (CHT), and bovine trypsin (TRY). DAB-acyl-enzyme deacylation rate constants, k3, were found to span a 720-fold range at pH 7.8 (DAB-WT > DAB-TRY > DAB-N155G > DAB-CHT). DAB-N155G was found to deacylate 80-fold slower than DAB-WT, indicating a 2.6 kcal/mol loss of transition-state binding energy due to this mutation. Absorbance spectra revealed strongly red-shifted absorbance lambda max values for all of the DAB-acyl-enzymes. The red shift was found to be 2.0 nm less in DAB-N155G, indicating that the oxyanion hole is partially responsible for this electronic perturbation of the DAB chromophore at the active site. Raman difference spectra of the DAB-acyl-enzymes measured at pH 5.0 and 8.6, with 18O-labeling of the carbonyl, show that the molecular motions most perturbed by the active site are three associated with the scissile acyl bond. Most interesting is the carbonyl stretching vibration, v(C = O), whose motion extends into the hydrolytic reaction coordinate. Comparison of the v(C = O) of DAB-WT and DAB-N155G reveals that the oxyanion hole does indeed form a hydrogen-bonding interaction with the carbonyl oxygen, the strength of which increases at pH 8.6. Interestingly, the DAB-TRY carbonyl forms very strong hydrogen bonds, even at pH 5.0, but DAB-CHT does not, even at pH 8.6. The low-frequency (1661 cm-1) v(C = O)'s of pH 5.0 DAB-TRY and pH 8.6 DAB-WT are proposed to correspond to a tetrahedrally distorted carbonyl center like that observed in the crystal structure of guanidinobenzoyl-TRY (Mangel et al., 1990). The strength of hydrogen bonding between the DAB-acyl-enzyme's carbonyl and the oxyanion hole, as gauged by the v(C = O) frequency, was found to correlate positively with an increased deacylation rate. This correlation, as well as calculated acyl-enzyme carbonyl bond lengths, which indicate a 0.015-A lengthening due to the oxyanion hole interaction, was found to be in good agreement with previously published resonance Raman data of alpha, beta-unsaturated arylacryloyl-acyl-enzymes (Tonge & Carey, 1990b, 1992).

Metabolite-modulated complex formation between alpha-glycerophosphate dehydrogenase and lactate dehydrogenase.

A modified Hummel-Dreyer equilibrium chromatography technique was used to test the hypothesis that NADH induces the molecular association of lactate dehydrogenase (LDH) and alpha-glycerol-3-phosphate dehydrogenase (alpha-GDH). In the presence of a very limited NADH concentration, a unique elution profile with a new peak running immediately ahead of a trough at the free alpha-GDH elution position is obtained. The appearance of this peak-trough profile is physical evidence that reversible association between LDH and alpha-GDH occurs over a very limited range of free NADH concentrations. The association constant for this complex formation between LDH and alpha-GDH is estimated to be 2.0 microM-1. With the NADH concentration increased to saturation level, no evidence of binding is observed. Such concentration-dependent behavior suggests that the strong competition between LDH and alpha-GDH for the limited amount of NADH tends to promote the enzyme-enzyme contact in order to make the most efficient use of the shared metabolite. The experimental results described in this article make a convincing argument for a metabolite-modulated enzyme-enzyme interaction along this metabolic pathway.

Base sequence criteria and Cartesian coordinates for stable B/Z and B/Z/B junctions in relaxed DNA.

It seems increasingly evident that if the Z form of DNA exists in the genome it must exist as short sections of alternating pyrimidine-purine sequences in the midst of very long sections of B-form DNA. We have determined the minimum length of a string of alternating CG base pairs that can go into the Z form in the middle of a long B form. Self-complimentary oligomers of the form T(M)(CG)(N)A(M) were synthesized. The conformation of the resulting duplex was determined in 6M aqueous NaCl solution by Raman scattering. We have found that 12 alternating CG base pairs is the minimum length required to form a stable Z form of DNA inside of a long B form section. Only the 4 center CG base pairs go into the Z form. These 4 CG base pairs in the Z form are flanked on each side by 4 CG base pairs in a non-Z (probably B) form as well as the ..TT.. ..AA.. sequences in the B form. We propose a model of the B/Z junction in which the double helix flips directly from the B form to the Z form so that there are no base pairs in the junction. In this model the B form is nucleated in the AT base pairs on each end and is propagated into the CG base pairs in the center. This model is supported by isotopic H/D exchange experiments that shows that the H/D exchange of the non-Z form CG base pairs is highly retarded and indicates that they remain in the B form. A Thermodynamic analysis of the concentration dependence of the melting point of the duplexes in both low and high salt, supports our model and rules out the possibility of hairpin formation. The enthalpy for the formation of a B/Z junction is determined to be about +16 kcal/junction. A comparison of these results with recent results on B/Z junctions in super-coiled DNA is given. Molecular modeling calculations permit us to obtain values for the coordinates and torsional angles of the oligomers showing both B/Z and B/Z/B junctions. The Cartesian coordinates for these oligomers as well as stereo figures of these models in color are available from the authors.

The time dependent UV resonance Raman spectra, conformation, and biological activity of acetylcholine analogues upon binding to acetylcholine binding proteins.

In order to obtain quantitative data on the relation between the conformation of acetylcholine and its interaction with biologically significant proteins, a series of acetylcholine analogues with absorption bands in the region 200-300 nm have been synthesized or obtained commercially. Each of these compounds were assayed to measure its activity as an ion channel activator of the nicotinic acetylcholine receptor protein (AChR). In addition, the suitability of some of these compounds as substrates for hydrolysis by acetylcholine esterase (AChE) was determined. One of these analogues, dimethylthionocarbamylcholine (DMTC-Ch), has the ester carbonyl oxygen replaced by a thionyl sulfur. DMTC-Ch has been found to be quite active as an ion channel activator when bound to AChR and was found to react with the enzyme AChE as a suicide substrate. It forms a thionoester of the serine at the AChE active site by an ester exchange reaction that releases the choline as the first product. However, the second or acid product is not released even at pH 7.5 over a period of days. This acetylcholine analog has an absorption band at about 240 nm and exhibits very strong ultraviolet resonance Raman (UVRR) spectra using 239 nm excitation from a frequency modified Nd:YAG laser. This technique allows observation of both conformational changes of the ligand molecule that result in frequency changes as well as changes in the excited state electronic structure that results in changes in the relative intensity of the Raman bands. The time dependence of the UVRR spectrum of the ligand upon binding to both AChE and AChR has been studied from 0.1 msec to minutes. Some time dependence in the conformation of DMTC-Ch upon binding to AChE has been found for very short (0.1-0.5 msec) times. However, no change in the conformation of this neurotransmitter analog is found in the available time range upon binding to AChR. From these data it is concluded that a previous suggestion that acetylcholine has a conformational change upon binding to AChR may be incorrect since the solution behavior of the carbamyl cholines and acetylcholine are similar. Even if acetylcholine does change conformation upont binding to AChR, it is unlikely that such a conformational change plays a significant role in channel activation. We present strong evidence that acetylcholine and its analogues can be active in a variety of conformations.(ABSTRACT TRUNCATED AT 400 WORDS)

Deuterium quadrupole coupling in N-acetylglycine and librational dynamics in solid poly(gamma-benzyl-L-glutamate).

To study the dynamics of peptide groups in solid proteins, we have accurately determined the principal components and molecular orientation of the electric field gradient tensor for the exchangeable deuterons in monoclinic N-acetylglycine by single-crystal deuterium nuclear magnetic resonance. These results are compared with the principal components of the amide deuterons in solid poly(gamma-benzyl-L-glutamate) measured in powder samples over a wide temperature range (140-400 K). The comparison indicates that in the solid polypeptide the N-D bonds undergo a small-amplitude torsional reorientation (libration) perpendicular to the peptide plane. To estimate dynamic rates, longitudinal relaxation times (T1 values) are reported for N-acetylglycine and poly(gamma-benzyl-L-glutamate). T1 values for the carboxyl and amide deuterons in N-acetylglycine are approximately 100 s, whereas for the amide deuterons in the polypeptide T1 approximately 1 s, also indicating that the N-D bonds are not stationary in the polypeptide. We determine from the reduced quadrupole coupling tensor the mean-square amplitude for the libration and show that it increases linearly with temperature. A simple qualitative theory for the relaxation times is presented on the basis of the assumption that the N-D reorientation is described either as a diffusion process in a square well or as a damped Langevin oscillator with a harmonic restoring force. The conclusion is that the short relaxation times of the polypeptide amide deuterons result from substantial frictional effects on reorientation that increase with temperature.

The hydrogen-bonding structure in parallel-stranded duplex DNA is reverse Watson-Crick.

Raman spectra of the parallel-stranded duplex formed from the deoxyoligonucleotides 5'-d-[(A)10TAATTTTAAATATTT]-3' (D1) and 5'-d[(T)10ATTAAAATTTATAAA]-3' (D2) in H2O and D2O have been acquired. The spectra of the parallel-stranded DNA are then compared to the spectra of the antiparallel double helix formed from the deoxyoligonucleotides D1 and 5'-d(AAATATTTAAAATTA-(T)10]-3' (D3). The Raman spectra of the antiparallel-stranded (aps) duplex are reminiscent of the spectra of poly[d(A)].poly[d(T)] and a B-form structure similar to that adopted by the homopolymer duplex is assigned to the antiparallel double helix. The spectra of the parallel-stranded (ps) and antiparallel-stranded duplexes differ significantly due to changes in helical organization, i.e., base pairing, base stacking, and backbone conformation. Large changes observed in the carbonyl stretching region (1600-1700 cm-1) implicate the involvement of the C(2) carbonyl of thymine in base pairing. The interaction of adenine with the C(2) carbonyl of thymine is consistent wtih formation of reverse Watson-Crick base pairing in parallel-stranded DNA. Phosphate-furanose vibrations similar to those observed for B-form DNA of heterogenous sequence and high A,T content are observed at 843 and 1092 cm-1 in the spectra of the parallel-stranded duplex. The 843-cm-1 band is due to the presence of a sizable population of furanose rings in the C2'-endo conformation. Significant changes observed in the regions from 1150 to 1250 cm-1 and from 1340 to 1400 cm-1 in the spectra of the parallel-stranded duplex are attributed to variations in backbone torsional and glycosidic angles and base stacking.

Raman spectroscopic studies of the DNA cro binding site conformation, free and bound to cro protein.

Raman spectra of the DNA binding site for cro repressor protein were obtained in the presence and absence of bound cro protein. The 17 base pair fragment is a consensus sequence of the six cro binding sites in phage lambda, except that the second base to the right of the center of pseudosymmetry is altered. Analysis of the spectrum of the free DNA indicates that the molecule exists in a B-like conformation with deviations from the usual B form occurring mainly in the bands assigned to A-T vibrations. The spectrum of the bound DNA was obtained by subtracting the spectrum of free cro from the spectrum of the complex which was estimated to be 90% bound. The DNA undergoes significant structural changes upon binding to the protein; most notable of these changes is a destacking of the G-C bases reflected by increases in the 1240, 1262, and 1320 cm-1 bands. A decrease in the 1361 cm-1 band that occurs has also been assigned to a destacking in guanine bases. The appearance of a 705 cm-1 band and the decrease and downshift of the 670 cm-1 band are consistent with the appearance of A-like character in the A-T region of the binding site when the protein binds; however, the spectra indicate that the entire binding site remains in a distorted B-like conformation. We use the 705 cm-1 band to estimate A-like character because the 800-850 cm-1 region is obscured by interference from strong protein bands. Other shifts in both intensity and position cannot be assigned to characteristic changes in conformation and therefore must be attributed to the protein influencing the structure in a novel way.

The infrared and Raman spectra of the duplex of d(GGTATACC) in the crystal show bands due to both the A-form and the B-form of DNA.

The deoxyoligonucleotide, d(GGTATACC), forms a duplex structure that crystallizes in the DNA A form. This has been shown by both X-ray diffraction studies and Raman spectroscopy (1,2). The presence of the DNA B form has been reported using diffuse X-ray scattering from a crystal of the closely related sequence d(GGBrUABrUACC)(3). In this paper the infrared spectrum of the d(GGTATACC) crystal is presented and curve resolution of both the Raman and IR spectra have been carried out. The percentage of A and B forms have been estimated. The %B form in the crystal has been estimated from the IR spectra to be about 15% and from Raman to be about 20%. Moreover the IR spectrum of the A conformation in the crystal is slightly different from the IR spectrum of the A conformation in polynucleotide fibers in particular in the region of the phosphate stretching vibrations and of the in-plane double bond vibrations of the bases. We show that it is feasible to obtain IR as well as Raman spectra of small crystals of oligonucleotides and that this is a good method of identifying all of the different conformations that may be in the crystal.

A Raman study of low frequency intrahelical modes in A-, B-, and C-DNA.

We have obtained low frequency (less than 200 cm-1) Raman spectra of calf-thymus DNA and poly(rI).poly(rC) as a function of water content and counterion species and of d(GGTATACC)2 and d(CGCGAATTCGCG)2 crystals. We have found that the Raman scattering from water in the first and second hydration shells does not contribute directly to the Raman spectra of DNA. We have determined the number of strong Raman active modes by comparing spectra for different sample orientations and polarizations and by obtaining fits to the spectra. We have found at least five Raman active modes in the spectra of A- and B-DNA. The frequencies of the modes above 40 cm-1 do not vary with counterion species, and there are only relatively small changes upon hydration. These modes are, therefore, almost completely internal. The mode near 34 cm-1 in A-DNA is mostly internal, whereas the mode near 25 cm-1 is dominated by interhelical interactions. The observed intensity changes upon dehydration were found to be due to the decrease in interhelical distance. Polymer length appears to play a role in the lowest frequency modes.

Low frequency Raman spectra of Z-DNA.

We have performed a Raman study of the low frequency modes in three oligo- and polynucleotides in Z-conformation, and we compare the spectra of these samples to those of two polynucleotides in B-conformation. In Z-DNA we find 5 intrahelical modes below 200 cm-1, in addition to the interhelical mode near 30 cm-1 which is only observed in crystalline samples. The most prominent intrahelical mode has a frequency of about 105 cm-1, close to the frequency of the strongest intrahelical mode in A-and B-DNA. The sequence dependence of the frequency of this mode is considerably larger than for the same mode in B-DNA. The other modes are less pronounced, and their frequency variations with base sequence are within the experimental accuracy.

A temperature-dependent Z to B to single-strand transition in d(CGCG).

A thermally induced left- to right-handed helical interconversion for the self-complementary tetradeoxynucleotide d(CGCG) has been observed in the presence of NaCl at concentrations from 1.5 to 2.25M. Analysis of the Raman spectrum of d(CGCG) in aqueous solution as a function of temperature indicates that at low temperature a left-handed Z-helical form is the predominant species. An increase in temperature results in a decrease in the population of the left-handed form and an increase in the population of the right-handed form. Further elevation in temperature results in extensive disruption of base stacking and a loss of secondary structure. This unstacking presumably represents dissociation to a single-stranded structure. The data suggest that a temperature-dependent Z to B to single-strand transition occurs under the conditions employed. A direct Z to single-strand conversion was not identified. In contrast to previous examples, the thermally induced left- to right-handed helical conversion for d(CGCG) can occur in the absence of chemically modified nucleic acid bases, alcohol solutions, or divalent ions.

The length of a junction between the B and Z conformations in DNA is three base pairs or less.

Recently it has been suggested that double-helical complexes formed between the DNA sequences (CG)n(A)m and their conjugates, (T)m(CG)n, would be candidates for the formation of a B-Z junction in aqueous solution at high salt concentrations [Peticolas et al. (1988) Proc. Natl. Acad. Sci. U.S.A. 85, 2579-2583]. The junction was predicted to occur between a B-type helix in the d(A)m.d(T)m section and a Z-type helix in the self-complementary (CG)n.(CG)n sequence. In this paper we report Raman experiments on the deoxyoligonucleotides d(CGCGCGCGCGCGAAAAA) and d(CGCGCGAAAAA) and their complements. It is found the latter compound cannot be induced into the Z form in saturated salt solution but that the former sequence goes into a B-Z junction at 5.5 M salt. From a comparison of the relative intensity of the Raman conformational marker bands for B and Z DNA for both the A-T and C-G base pairs, it is shown that in 5.5 M NaCl solution none of the A-T base pairs are in the Z form, but nine of the C-G base pairs are in the Z form. The remaining three C-G base pairs are either in the junction or in the B form. Thus, the junction is formed from three or less C-G base pairs. If the solution is made 95 microM with NiCl2, then the entire duplex goes into the Z form and the Raman bands of the adenine are completely changed into those of the Z form.(ABSTRACT TRUNCATED AT 250 WORDS)

A duplex of the oligonucleotides d(GGGGGTTTTT) and d(AAAAACCCCC) forms an A to B conformational junction in concentrated salt solutions.

The coexistence of both A form and B form tracts and formation of an A-B junction in the oligomer d(GGGGGTTTTT).d(AAAAACCCCC) in saturated sodium chloride solution have been detected by Raman spectroscopy. The entire duplex adopts the familiar B-form conformation in aqueous solution at low salt concentrations (0.1M NaCl). In 6M NaCl the adoption of an A form is observed within the G,C tract while a B-form is maintained in the A.T tract. The experimental results indicate that two different helical forms can co-exist in a rather short oligonucleotide and that formation of an A-B junction can occur over a fairly small span of bases. This is in agreement with recent rules governing the relation between base sequence and secondary structure of DNA published from this laboratory. The conformational preferences of each of the individual oligomers d(AAAAACCCCC) and d(GGGGGTTTTT) have also been investigated. The oligomer d(AAAAACCCCC) is single stranded but some evidence for base stacking is observed at 2 degrees C. In contrast, a double stranded B-form structure characterized by wobble G-T base pairing is observed for d(GGGGGTTTTT) in 0.1M and 6M NaCl.

Environmentally induced conformational changes in B-type DNA: comparison of the conformation of the oligonucleotide d(TCGCGAATTCGCG) in solution and in its crystalline complex with the restriction nuclease EcoRI.

Raman spectroscopic analysis of the secondary structure of the crystalline restriction endonuclease EcoRI, the oligonucleotide d(TCGCGAATTCGCG) in solution, and the corresponding crystalline EcoRI-oligonucleotide complex reveals structural differences between the complexed and uncomplexed protein and oligonucleotide components that appear to be linked to complex formation. Structural differences that are spectroscopically identified include (1) an increase in the population of furanose rings adopting the C3'-endo conformation and (2) spectroscopically observed changes in base stacking which are probably associated with the crystallographically observed distortion of the phosphate backbone about positions C(3)-G(4) and C(9)-G(10) and unwinding between the symmetry-related segments GAA-TTC which make up the central recognition core (McClarin et al., 1986). Changes in base stacking due to distortions and unwinding along the oligonucleotide result in differences in the base vibrational region between the spectra of the complex and the oligonucleotide in solution. The spectroscopic analysis indicates that the C2'-endo population is similar for the oligonucleotide in solution and in the complex. The additional C3'-endo population in the complex appears to arise from the conversion of rings adopting alternative conformations such as C1'-exo and O1'-endo. Analysis of the vibrational bands derived from guanine indicates that the population of guanine residues associated with furanose rings in a C2'-endo conformation is similar for the oligonucleotide in solution and in the crystalline complex. This implies that the increase in C3'-endo population is not associated with guanine residues. Large conformational distortions such as those observed in the crystal distortions are not observed in either the crystal or the solution of the oligomer d(CGCGAATTCGCG).(ABSTRACT TRUNCATED AT 250 WORDS)

Some rules for predicting the base-sequence dependence of DNA conformation.

Two tables have been constructed showing the crystal and solution conformations of short sequences of DNA. Each of these DNAs has been found to be in one of three different conformations--the A, B, or Z form--depending upon the base sequence and the environmental conditions. A set of rules is presented showing the tendency of certain base pairs to direct the DNA conformation into the A, B, or Z genus in saturated salt solutions and in crystals. These rules are based on a consideration of nearest-neighbor interactions that are interpreted in terms of 10 different two-letter code words made from the letters denoting the bases guanine (G), cytosine (C), adenine (A), and thymine (T). One table discusses the effect on DNA conformation of 3 strong words that tend to direct a DNA oligomer into either the A, B, or Z genus in crystals or in aqueous solutions containing a high salt concentration (6 M). The second table discusses the remaining 7 code words that appear to have a much weaker effect on conformation. The sequences that are most likely to lead to A-Z, B-Z, and A-B junctions are discussed, as is the possible biological significance of these rules.