Thermodynamics of stacking and of self-association of the dinucleoside monophosphate m2(6)A-U from proton NMR chemical shifts: differential concentration temperature profile method

Structure and dynamics of the proton-selective histidine and the gating tryptophan in an inward rectifying hybrid influenza B and A virus M2 proton channel

The M2 proteins of influenza A and B viruses form acid-activated proton channels that are essential for the virus lifecycle. Proton selectivity is achieved by a transmembrane (TM) histidine whereas gating is achieved by a tryptophan residue. Although this functional apparatus is conserved between AM2 and BM2 channels, AM2 conducts protons exclusively inward whereas BM2 conducts protons in either direction depending on the pH gradient. Previous studies showed that in AM2, mutations of D44 abolished inward rectification of AM2, suggesting that the tryptophan gate is destabilized. To elucidate how charged residues C-terminal to the tryptophan regulates channel gating, here we investigate the structure and dynamics of H19 and W23 in a BM2 mutant, GDR-BM2, in which three BM2 residues are mutated to the corresponding AM2 residues, S16G, G26D and H27R. Whole-cell electrophysiological data show that GDR-BM2 conducts protons with inward rectification, identical to wild-type (WT) AM2 but different from WT-BM2. Solid-state NMR 15N and 13C spectra of H19 indicate that the mutant BM2 channel contains higher populations of cationic histidine and neutral τ tautomers compared to WT-BM2 at acidic pH. Moreover, 19F NMR spectra of 5-19F-labeled W23 resolve three peaks at acidic pH, suggesting three tryptophan sidechain conformations. Comparison of these spectra with the tryptophan spectra of other M2 peptides suggests that these indole sidechain conformations arise from interactions with the C-terminal charged residues and with the N-terminal cationic histidine. Taken together, these solid-state NMR data show that inward rectification in M2 proton channels is accomplished by tryptophan interactions with charged residues on both its C-terminal and N-terminal sides. Gating of these M2 proton channels is thus accomplished by a multi-residue complex with finely tuned electrostatic and aromatic interactions.

Disulfide engineering on temporin-SHf: Stabilizing the bioactive conformation of an ultra-short antimicrobial peptide

In Silico searching for short antimicrobial peptides has revealed temporin-SHf as the short (8AA), hydrophobic, broad spectrum, and natural antimicrobial peptide. Important drawback associated with temporin-SHf is the susceptibility of its bioactive conformation for denaturation and proteolytic degradation. In the current report, disulfide engineering strategy has been adopted to improve the stability of bioactive conformation of temporin-SHf. The functionally non-critical Leu4 and Ile7 residues at i and i + 3 position of helical conformation of temporin-SHf were mutated with cysteine disulfide. Designed [L4C, I7C]temporin-SHf was synthesized, characterized using NMR spectroscopy, and accessed for antimicrobial activity. [L4C, I7C]Temporin-SHf adopts helical conformation from Phe3 to Phe8 in the absence of membrane-mimetic environment and retains broad spectrum antimicrobial activity. The reduction potential of cysteine disulfide of [L4C, I7C]temporin-SHf is -289 mV. Trypsin-induced digestion and serum-induced digestion have confirmed the advantage of cysteine disulfide in imparting proteolytic stability to temporin-SHf. Disulfide-stabilized temporin-SHf may serve as a good model for the rational design of temporin-SHf based antibiotics for treatment of infectious diseases.

Proton detection for signal enhancement in solid-state NMR experiments on mobile species in membrane proteins

Direct proton detection is becoming an increasingly popular method for enhancing sensitivity in solid-state nuclear magnetic resonance spectroscopy. Generally, these experiments require extensive deuteration of the protein, fast magic angle spinning (MAS), or a combination of both. Here, we implement direct proton detection to selectively observe the mobile entities in fully-protonated membrane proteins at moderate MAS frequencies. We demonstrate this method on two proteins that exhibit different motional regimes. Myelin basic protein is an intrinsically-disordered, peripherally membrane-associated protein that is highly flexible, whereas Anabaena sensory rhodopsin is composed of seven rigid transmembrane α-helices connected by mobile loop regions. In both cases, we observe narrow proton linewidths and, on average, a 10× increase in sensitivity in 2D insensitive nuclear enhancement of polarization transfer-based HSQC experiments when proton detection is compared to carbon detection. We further show that our proton-detected experiments can be easily extended to three dimensions and used to build complete amino acid systems, including sidechain protons, and obtain inter-residue correlations. Additionally, we detect signals which do not correspond to amino acids, but rather to lipids and/or carbohydrates which interact strongly with membrane proteins.

In vitro biosynthesis and chemical identification of UDP-N-acetyl-d-quinovosamine (UDP-d-QuiNAc)

N-acetyl-d-quinovosamine (2-acetamido-2,6-dideoxy-d-glucose, QuiNAc) occurs in the polysaccharide structures of many Gram-negative bacteria. In the biosynthesis of QuiNAc-containing polysaccharides, UDP-QuiNAc is the hypothetical donor of the QuiNAc residue. Biosynthesis of UDP-QuiNAc has been proposed to occur by 4,6-dehydration of UDP-N-acetyl-d-glucosamine (UDP-GlcNAc) to UDP-2-acetamido-2,6-dideoxy-d-xylo-4-hexulose followed by reduction of this 4-keto intermediate to UDP-QuiNAc. Several specific dehydratases are known to catalyze the first proposed step. A specific reductase for the last step has not been demonstrated in vitro, but previous mutant analysis suggested that Rhizobium etli gene wreQ might encode this reductase. Therefore, this gene was cloned and expressed in Escherichia coli, and the resulting His6-tagged WreQ protein was purified. It was tested for 4-reductase activity by adding it and NAD(P)H to reaction mixtures in which 4,6-dehydratase WbpM had acted on the precursor substrate UDP-GlcNAc. Thin layer chromatography of the nucleotide sugars in the mixture at various stages of the reaction showed that WbpM converted UDP-GlcNAc completely to what was shown to be its 4-keto-6-deoxy derivative by NMR and that addition of WreQ and NADH led to formation of a third compound. Combined gas chromatography-mass spectrometry analysis of acid hydrolysates of the final reaction mixture showed that a quinovosamine moiety had been synthesized after WreQ addition. The two-step reaction progress also was monitored in real time by NMR. The final UDP-sugar product after WreQ addition was purified and determined to be UDP-d-QuiNAc by one-dimensional and two-dimensional NMR experiments. These results confirmed that WreQ has UDP-2-acetamido-2,6-dideoxy-d-xylo-4-hexulose 4-reductase activity, completing a pathway for UDP-d-QuiNAc synthesis in vitro.

Molecular promoting of aluminum metal-organic framework topology MIL-101 by N,N-dimethylformamide

In situ NMR and DFT modeling demonstrate that N,N-dimethylformamide (DMF) promotes the formation of metal-organic framework NH2-MIL-101(Al). In situ NMR studies show that upon dissociation of an aluminum-coordinated aqua ligand in NH2-MOF-235(Al), DMF forms a H-Cl-DMF complex during synthesis. This reaction induces a transformation from the MOF-235 topology into the MIL-101 topology. Electronic structure density functional theory (DFT) calculations show that the use of DMF instead of water as the synthesis solvent decreases the energy gap between the kinetically favored MIL-101 and thermodynamically favored MIL-53 products. DMF therefore promotes MIL-101 topology both kinetically and thermodynamically.

Measurement of pH by NMR Spectroscopy in Concentrated Aqueous Fluoride Buffers

An NMR spectroscopic technique has been developed to give rapid, accurate pH measurements on tenth-milliliter samples of concentrated acidic aqueous solutions buffered by fluoride ion in the pH 1.5 - 4.5 range. The fluoride (19)F chemical shift has been calibrated as a function of pH at 0.1 and 1.0 M concentration by reference to an internal 3-fluoropyridine standard. Subsequent measurements of fluoride buffer pH required no additives and only two NMR spectra in the presence of an external reference standard.

Probing microsecond time scale dynamics in proteins by methyl (1)H Carr-Purcell-Meiboom-Gill relaxation dispersion NMR measurements. Application to activation of the signaling protein NtrC(r)

To study microsecond processes by relaxation dispersion NMR spectroscopy, low power deposition and short pulses are crucial and encourage the development of experiments that employ (1)H Carr-Purcell-Meiboom-Gill (CPMG) pulse trains. Herein, a method is described for the comprehensive study of microsecond to millisecond time scale dynamics of methyl groups in proteins, exploiting their high abundance and favorable relaxation properties. In our approach, protein samples are produced using [(1)H, (13)C]-d-glucose in ∼100% D(2)O, which yields CHD(2) methyl groups for alanine, valine, threonine, isoleucine, leucine, and methionine residues with high abundance, in an otherwise largely deuterated background. Methyl groups in such samples can be sequence-specifically assigned to near completion, using (13)C TOCSY NMR spectroscopy, as was recently demonstrated (Otten, R.; et al. J. Am. Chem. Soc. 2010, 132, 2952-2960). In this Article, NMR pulse schemes are presented to measure (1)H CPMG relaxation dispersion profiles for CHD(2) methyl groups, in a vein similar to that of backbone relaxation experiments. Because of the high deuteration level of methyl-bearing side chains, artifacts arising from proton scalar coupling during the CPMG pulse train are negligible, with the exception of Ile-δ1 and Thr-γ2 methyl groups, and a pulse scheme is described to remove the artifacts for those residues. Strong (13)C scalar coupling effects, observed for several leucine residues, are removed by alternative biochemical and NMR approaches. The methodology is applied to the transcriptional activator NtrC(r), for which an inactive/active state transition was previously measured and the motions in the microsecond time range were estimated through a combination of backbone (15)N CPMG dispersion NMR spectroscopy and a collection of experiments to determine the exchange-free component to the transverse relaxation rate. Exchange contributions to the (1)H line width were detected for 21 methyl groups, and these probes were found to collectively report on a local structural rearrangement around the phosphorylation site, with a rate constant of (15.5 ± 0.5) × 10(3) per second (i.e., τ(ex) = 64.7 ± 1.9 μs). The affected methyl groups indicate that, already before phosphorylation, a substantial, transient rearrangement takes place between helices 3 and 4 and strands 4 and 5. This conformational equilibrium allows the protein to gain access to the active, signaling state in the absence of covalent modification through a shift in a pre-existing dynamic equilibrium. Moreover, the conformational switching maps exactly to the regions that differ between the solution NMR structures of the fully inactive and active states. These results demonstrate that a cost-effective and quantitative study of protein methyl group dynamics by (1)H CPMG relaxation dispersion NMR spectroscopy is possible and can be applied to study functional motions on the microsecond time scale that cannot be accessed by backbone (15)N relaxation dispersion NMR. The use of methyl groups as dynamics probes extends such applications also to larger proteins.

Chemical-shift referencing and resolution stability in gradient LC-NMR (acetonitrile:water)

Wide ranges of solvent conditions are generated during solvent-gradient LC-NMR. This complicates the referencing of the chemical-shift scale of the resulting NMR data. The problems that arise when performing LC-NMR in acetonitrile:water - particularly when the mobile-phase composition can range anywhere from 0% to 100% - are examined here, and the reliability of the secondary reference signals are evaluated. It is shown that under these conditions the use of the acetonitrile signal is superior to the use of the water signal in any form (either the (1)H or the (2)H signal) as a secondary reference, a lock signal, and a signal for shimming. The limitations of the referencing methods and other experimental parameters, and the limitations on the solvent-gradient ramp parameters, primarily as they affect lineshapes, are all shown. These results are compared to the way some other publications have referenced the (1)H chemical-shift axis (when using acetonitrile:water mixtures to perform reversed-phase chromatography LC-NMR).

The T-lock: automated compensation of radio-frequency induced sample heating

Modern high-field NMR spectrometers can stabilize the nominal sample temperature at a precision of less than 0.1 K. However, the actual sample temperature may differ from the nominal value by several degrees because the sample heating caused by high-power radio frequency pulses is not readily detected by the temperature sensors. Without correction, transfer of chemical shifts between different experiments causes problems in the data analysis. In principle, the temperature differences can be corrected by manual procedures but this is cumbersome and not fully reliable. Here, we introduce the concept of a "T-lock", which automatically maintains the sample at the same reference temperature over the course of different NMR experiments. The T-lock works by continuously measuring the resonance frequency of a suitable spin and simultaneously adjusting the temperature control, thus locking the sample temperature at the reference value. For three different nuclei, 13C, 17O and 31P in the compounds alanine, water, and phosphate, respectively, the T-lock accuracy was found to be <0.1 K. The use of dummy scan periods with variable lengths allows a reliable establishment of the thermal equilibrium before the acquisition of an experiment starts.

Dinucleotide TpT and its 2'-O-Me analogue possess different backbone conformations and flexibilities but similar stacked geometries

UV irradiation at 254 nm of 2'-O,5-dimethyluridylyl(3'-5')-2'-O,5-dimethyluridine (1a) and of natural thymidylyl(3'-5')thymidine (1b) generates the same photoproducts (CPD and (6-4)PP; responsible for cell death and skin cancer). The ratios of quantum yields of photoproducts obtained from 1a (determined herein) to that from 1b are in a proportion close to the approximately threefold increase of stacked dinucleotides for 1a compared with those of 1b (from previous circular dichroism results). 1a and 1b however are endowed with different predominant sugar conformations, C3'-endo (1a) and C2'-endo (1b). The present investigation of the stacked conformation of these molecules, by unrestrained state-of-the-art molecular simulation in explicit solvent and salt, resolves this apparent paradox and suggests the following main conclusions. Stacked dinucleotides 1a and 1b adopt the main characteristic features of a single-stranded A and B form, respectively, where the relative positions of the backbone and the bases are very different. Unexpectedly, the geometry of the stacking of two thymine bases, within each dinucleotide, is very similar and is in excellent agreement with photochemical and circular dichroism results. Analyses of molecular dynamics trajectories with conformational adiabatic mapping show that 1a and 1b explore two different regions of conformational space and possess very different flexibilities. Therefore, even though their base stacking is very similar, these molecules possess different geometrical, mechanical, and dynamical properties that may account for the discrepancy observed between increased stacking and increased photoproduct formations. The computed average stacked conformations of 1a and 1b are well-defined and could serve as starting models to investigate photochemical reactions with quantum dynamics simulations.

Structural Evaluation of a Novel Pro-apoptotic Peptide Coupled to CNGRC Tumor Homing Sequence by NMR

Hunter-killer peptides (HKPs) are synthetic peptides that target specific cell types for apoptosis. These studies report functional and structural characteristics of HKP9, an hunter-killer peptide that specifically targets tumor vasculature with a new apoptotic sequence. Vesicle leakage experiments were performed as a model for membrane perturbing activity. Placement of the homing sequence reduces both cell toxicity and vesicle leakage activity. NMR studies elucidate the conformation and orientation of HKP9 in micelles. The positively charged end of the HKP9 killing sequence is solvent exposed; however, the central portion of the peptide is helical and buried in dodecylphosphorylcholine micelles. The homing sequence is less solvent exposed than in a previously reported tumor-homing peptide. The results suggest that solvent accessibility of the homing sequence should be considered in design of future peptides.

Calcium-modulated S100 protein-phospholipid interactions. An NMR study of calbindin D9k and DPC

The cellular functions of several S100 proteins involve specific interactions with phospholipids and the cell membrane. The interactions between calbindin D(9k) (S100D) and the detergent dodecyl phosphocholine (DPC) were studied using NMR spectroscopy. In the absence of Ca(2+), the protein associates with DPC micelles. The micelle-associated state has intact helical secondary structures but no apparent tertiary fold. At neutral pH, Ca(2+)-loaded calbindin D(9k) does not associate with DPC micelles. However, a specific interaction is observed with individual DPC molecules at a site close to the linker between the two EF-hands. Binding to this site occurs only when Ca(2+) is bound to the protein. A reduction in pH in the absence of Ca(2+) increases the stability of the micelle-associated state. This along with the corresponding reduction in Ca(2+) affinity causes a transition to the micelle-associated state also in the presence of Ca(2+) when the pH is lowered. Site-specific analysis of the data indicates that calbindin D(9k) has a core of three tightly packed helices (A, B, and D), with a dynamic fourth helix (C) more loosely associated. Evidence is presented that the Ca(2+)-binding characteristics of the two EF-hands are distinctly different in a micelle environment. The role of calbindin D(9k) in the cell is discussed, along with the broader implications for the function of the S100 protein family.

Orientation and helical conformation of a tissue-specific hunter-killer peptide in micelles

Hunter-killer peptides are chimeric synthetic peptides that selectively target specific cell types for an apoptotic death. These peptides, which are models for potential therapeutics, contain a homing sequence for receptor-mediated interactions and a pro-apoptotic sequence. Homing domains have been designed to target angiogenic tumor cells, prostate cells, arthritic tissue and, most recently, adipose tissue. After a receptor-mediated internalization, the apoptotic sequence, which contains D-enantiomer amino acids, initiates apoptosis through mitochondrial membrane disruption. We have begun structure and functional studies on a peptide (HKP1) that specifically targets angiogenic tumor cells for apoptosis. As a model for mitochondrial membrane disruption, we have examined peptide-induced leakage of a calcein fluorophore from large unilamellar vesicles. These experiments demonstrate more potent leakage activity by HKP1 than the peptide lacking the homing domain. Circular dichroism and 2D homonuclear NMR experiments demonstrate that this tumor-specific HKP adopts a left-handed amphipathic helix in association with dodecylphosphorylcholine micelles in a parallel orientation to the lipid-water interface with the homing domain remaining exposed to solvent. The amphipathic helix of the apoptotic domain orients with nonpolar leucine and alanine residues inserting most deeply into the lipid environment.

Correlation of three-dimensional structures with the antibacterial activity of a group of peptides designed based on a nontoxic bacterial membrane anchor

To understand the functional differences between a nontoxic membrane anchor corresponding to the N-terminal sequence of the Escherichia coli enzyme IIA(Glc) and a toxic antimicrobial peptide aurein 1.2 of similar sequence, a series of peptides was designed to bridge the gap between them. An alteration of a single residue of the membrane anchor converted it into an antibacterial peptide. Circular dichroism spectra indicate that all peptides are disordered in water but helical in micelles. Structures of the peptides were determined in membrane-mimetic micelles by solution NMR spectroscopy. The quality of the distance-based structures was improved by including backbone angle restraints derived from a set of chemical shifts ((1)H(alpha), (15)N, (13)C(alpha), and (13)C(beta)) from natural abundance two-dimensional heteronuclear correlated spectroscopy. Different from the membrane anchor, antibacterial peptides possess a broader and longer hydrophobic surface, allowing a deeper penetration into the membrane, as supported by intermolecular nuclear Overhauser effect cross-peaks between the peptide and short chain dioctanoyl phosphatidylglycerol. An attempt was made to correlate the NMR structures of these peptides with their antibacterial activity. The activity of this group of peptides does not correlate exactly with helicity, amphipathicity, charge, the number of charges, the size of the hydrophobic surface, or hydrophobic transfer free energy. However, a correlation is established between the peptide activity and membrane perturbation potential, which is defined by interfacial hydrophobic patches and basic residues in the case of cationic peptides. Indeed, (31)P solid state NMR spectroscopy of lipid bilayers showed that the extent of lipid vesicle disruption by these peptides is proportional to their membrane perturbation potential.

NMR Structure of the (+)-CPI-indole/d(GACTAATTGAC)-d(GTCAATTAGTC) Covalent Complex

We report the NMR solution structure of (+)-CPI-indole (CPI, 1,2,8,8a-tetrahydrocyclopropa[c]pyrrolo[3,2-e]indol-4(5H)-one), an agent belonging to the CC-1065/duocarmycin family of antitumor compounds. This (+)-CPI-indole structure is covalently bound to d(G(1)ACTAATTGTC(11))-d(G(12)TCAATTAGTC(22)), a synthetic DNA duplex containing a high-affinity binding site. The three-dimensional structure has been determined by several cycles of restrained molecular dynamics calculations with a total of 563 NMR-derived constraints, both in vacuo and by using the generalized Born solvent continuum model. In-depth analysis of the structure of this ligand-DNA complex led to a detailed knowledge of the bound state conformation of the CPI-indole, the most simplified agent related to CC-1065 and duocarmycins, the parent members of a family of extremely potent antitumor compounds. Comparison of the CPI-indole bound conformation with those previously found for (+)-duocarmycin SA (DSA), its unnatural enantiomer (-)-DSA, and the demethoxylated analogue (+)-DSI in their DNA complexes provided additional evidence of the tight correlation between the catalytic effect exerted by DNA on the alkylation reaction and the extent of angular twist between the two planar heteroaromatic subunits of these agents. Additionally, comparison of the structural features of the DNA-bound state of a "naked" ligand, such as CPI-indole, with those of various other duocarmycin agents provided useful information for the interpretation of the observed effects on chemical reactivity of the different substitution patterns at the hemispheres of these types of complex.

Solution structure of the N-terminal amphitropic domain of Escherichia coli glucose-specific enzyme IIA in membrane-mimetic micelles

The N-terminal domain of enzyme IIA(Glc) of the Escherichia coli phosphoenolpyruvate:sugar phosphotransferase system confers amphitropism to the protein, allowing IIA(Glc) to shuttle between the cytoplasm and the membrane. To further understand this amphitropic protein, we have elucidated, by NMR spectroscopy, the solution structure of a synthetic peptide corresponding to the N-terminal domain of IIA(Glc). In water, this peptide is predominantly disordered, consistent with previous data obtained in the absence of membranes. In detergent micelles of dihexanoylphosphatidylglycerol (DHPG) or sodium dodecylsulfate (SDS), however, residues Phe 3-Val 10 of the peptide adopt a helical conformation in the ensemble of structures calculated on the basis of NOE-derived distance restraints. The root mean square deviations for superimposing the backbone atoms of the helical region are 0.18 A in DHPG and 0.22 A in SDS. The structure, chemical shifts, and spin-spin coupling constants all indicate that, of the four lysines in the N-terminal domain of IIA(Glc), only Lys 5 and Lys 7 in the amphipathic helical region interact with DHPG. In addition, the peptide-detergent interactions were investigated using intermolecular NOESY experiments. The aliphatic chains of anionic detergents DHPG, SDS, and 2,2-dimethyl-2-silapentane-5-sulfonate sodium salt (DSS) all showed intermolecular NOE cross-peaks to the peptide, providing direct evidence for the putative membrane anchor of IIA(Glc) in binding to the membrane-mimicking micelles.

Mechanism of flavin mononucleotide cofactor binding to the Desulfovibrio vulgaris flavodoxin. 2. Evidence for cooperative conformational changes involving tryptophan 60 in the interaction between the phosphate- and ring-binding subsites

A mechanism has been proposed for the binding of flavin mononucleotide (FMN) and riboflavin to the apoflavodoxin from Desulfovibrio vulgaris [Murray, T. A., and Swenson, R. P. (2003) Biochemistry 42, 2307-2316]. In this model, the binding of the flavin isoalloxazine ring is dependent on the presence of a phosphate moiety in the phosphate-binding subsite, suggesting a cooperative interaction between that region and the ring-binding subsite. In the absence of inorganic phosphate, FMN can bind through the initial association of its 5'-phosphate group in the phosphate-binding subsite followed by insertion of the flavin ring. Because riboflavin lacks the 5'-phosphate group, it is unable to bind to this apoprotein in the absence of inorganic phosphate in solution. However, inorganic phosphate can enhance the rate of ring binding by occupying the phosphate-binding subsite. In this paper, NMR, near-UV circular dichroism (CD), and fluorescence spectroscopy provide evidence for a phosphate-induced conformational change within the isoalloxazine ring-binding subsite. Phosphate-dependent changes in the chemical shifts of 22 amide groups were observed in (1)H-(15)N HSQC NMR spectra. The majority of these groups are proximal to the phosphate-binding subsite or the loops that constitute the isoalloxazine ring-binding site. Also, a phosphate-dependent change in the environment or position of the Trp60 side chain was apparent in the NMR data and was confirmed by associated changes in the near-UV CD and tryptophan fluorescence spectra when compared to the spectra of the W60A mutant. These data suggest that phosphate, either the 5'-phosphate of the FMN or inorganic phosphate from solution, facilitates the movement of the side chain of Trp60 out of the isoalloxazine ring-binding site and other associated conformational changes, creating a population of apoflavodoxin that is capable of binding the isoalloxazine ring. This conformational switch may explain why some apoflavodoxins cannot bind riboflavin and also supports the "aromatic gate" model proposed from the crystal structure of the Anabaena apoflavodoxin [Genzor, C. G., Perales-Alcon, A., Sancho, J., and Romero, A. (1996) Nat. Struct. Biol. 3, 329-332].

The depletion of protein signals in metabonomics analysis with the WET-CPMG pulse sequence

Nuclear magnetic resonance (NMR) spectroscopy is a powerful analytical tool capable of providing a comprehensive metabolic profile of biofluids such as urine, plasma, and serum. Unfortunately, when measuring serum and plasma, the high protein concentration can obscure the signals originating from low molecular weight metabolites. We evaluated the use of different parameters within the Carr-Purcell-Meiboom-Gill (CPMG) pulse train of fast spin-echoes to remove the macromolecular signal contribution in one-dimensional proton (1H) NMR spectra. Experimental parameters such as the refocusing delay in the CPMG pulse train, pulse miscalibration, and recycle time were examined to assess the ability to remove the protein signals from the spectrum without causing a deleterious effect on the signals originating from free, low molecular weight metabolites. The 1H-NMR spectra of a variety of serum samples spiked with 2'-deoxyadenosine were acquired using various acquisition parameters. Our results show that the delay used in the CPMG spin-echo and the combination of the acquisition pulse flip angle and recycle time are the two major factors affecting the observed metabolite signal amplitudes in the resulting 1H-NMR spectrum.

Ca2+ binding and conformational changes in a calmodulin domain

Calcium activation of the C-terminal domain of calmodulin was studied using 1H and 15N NMR spectroscopy. The important role played by the conserved bidentate glutamate Ca2+ ligand in the binding loops is emphasized by the striking effects resulting from a mutation of this glutamic acid to a glutamine, i.e. E104Q in loop III and E140Q in loop IV. The study involves determination of Ca2+ binding constants, assignments, and structural characterizations of the apo, (Ca2+)1, and (Ca2+)2 states of the E104Q mutant and comparisons to the wild-type protein and the E140Q mutant [Evenäs et al. (1997) Biochemistry 36, 3448-3457]. NMR titration data show sequential Ca2+ binding in the E104Q mutant. The first Ca2+ binds to loop IV and the second to loop III, which is the order reverse to that observed for the E140Q mutant. In both mutants, the major structural changes occur upon Ca2+ binding to loop IV, which implies a different response to Ca2+ binding in the N- and C-terminal EF-hands. Spectral characteristics show that the (Ca2+)1 and (Ca2+)2 states of the E104Q mutant undergo global exchange on a 10-100 micros time scale between conformations seemingly similar to the closed and open structures of this domain in wild-type calmodulin, paralleling earlier observations for the (Ca2+)2 state of the E140Q mutant, indicating that both glutamic acid residues, E104 and E140, are required for stabilization of the open conformation in the (Ca2+)2 state. To verify that the NOE constraints cannot be fulfilled in a single structure, solution structures of the (Ca2+)2 state of the E104Q mutant are calculated. Within the ensemble of structures the precision is good. However, the clearly dynamic nature of the state, a large number of violated distance restraints, ill-defined secondary structural elements, and comparisons to the structures of calmodulin indicate that the ensemble does not provide a good picture of the (Ca2+)2 state of the E104Q mutant but rather represents the distance-averaged structure of at least two distinct different conformations.

High resolution solution structure of a DNA duplex alkylated by the antitumor agent duocarmycin SA

The three-dimensional solution structure of duocarmycin SA in complex with d-(G1ACTAATTGAC11).d-(G12TCATTAGTC22) has been determined by restrained molecular dynamics and relaxation matrix calculations using experimental NOE distance and torsion angle constraints derived from 1H NMR spectroscopy. The final input data consisted of a total of 858 distance and 189 dihedral angle constraints, an average of 46 constraints per residue. In the ensemble of 20 final structures, there were no distance constraint violations >0.06 A or torsion angle violations >0.8 degrees. The average pairwise root mean square deviation (RMSD) over all 20 structures for the binding site region is 0.57 A (average RMSD from the mean: 0.39 A). Although the DNA is very B-like, the sugar-phosphate backbone torsion angles beta, epsilon, and zeta are distorted from standard values in the binding site region. The structure reveals site-specific bonding of duocarmycin SA at the N3 position of adenine 19 in the AT-rich minor groove of the duplex and binding stabilization via hydrophobic interactions. Comparisons have been made to the structure of a closely related complex of duocarmycin A bound to an AT-rich DNA duplex. These results provide insights into critical aspects of the alkylation site selectivity and source of catalysis of the DNA alkylating agents, and the unusual stability of the resulting adducts.

Electron-transfer properties of Pseudomonas aeruginosa [Lys44, Glu64]azurin

In the hydrophobic patch of azurin from Pseudomonas aeruginosa, an electric dipole was created by changing Met44 into Lys and Met64 into Glu. The effect of this dipole on the electron-transfer properties of azurin was investigated. From a spectroscopic characterization (NMR, EPR and ultraviolet-visible) it was found that both the copper site and the overall structure of the [Lys44, Glu64]azurin were not disturbed by the two mutations. A small perturbation of the active site at high pH, similar to that observed for [Lys44]azurin, occurs in the double mutant. At neutral pH the electron-self-exchange rate constant of the double mutant shows a decrease of three orders of magnitude compared with the wild-type value. The possible reasons for this decrease are discussed. Electron transfer with the proposed physiological redox partners cytochrome c551 and nitrite reductase have been investigated and the data analyzed in the Marcus framework. From this analysis it is confirmed that the hydrophobic patch of azurin is the interaction site with both partners, and that cytochrome c551 uses its hydrophobic patch and nitrite reductase a negatively charged surface area for the electron transfer.

Kinetics of DNA hydration

The hydration of the d(CGCGAATTCGCG) B-DNA duplex in solution was studied by nuclear magnetic relaxation dispersion (NMRD) of the water nuclei 1H, 2H, and 17O, and by nuclear Overhauser effects (NOEs) in high-resolution two-dimensional 1H NMR spectra. By comparing results from the free duplex with those from its complex with netropsin, water molecules in the "spine of hydration" in the AATT region of the minor groove could be distinguished from hydration water elsewhere in the duplex. The 2H and 17O relaxation dispersions yield a model-independent residence time of 0.9(+/-0.1) ns at 4 degrees C for five highly ordered water molecules in the spine. When corrected for frequency offset effects, the NOE data yield the same residence time as the NMRD data, giving credence to both methods. At 27 degrees C, the residence time is estimated to 0.2 ns, a factor of 40 shorter than the tumbling time of the duplex. The NMRD data show that all water molecules associated with the duplex, except the five molecules in the spine, have residence times significantly shorter than 1 ns at 4 degrees C. There is thus no long-lived hydration structure associated with the phosphate backbone. In contrast to 2H and 17O, the 1H relaxation dispersion is dominated by labile DNA protons and therefore provides little information about DNA hydration.

NMR studies of the E140Q mutant of the carboxy-terminal domain of calmodulin reveal global conformational exchange in the Ca2+-saturated state

In the present investigation, the Ca2+ activation of the C-terminal domain of bovine calmodulin and the effects of replacing the bidentate Ca2+-coordinating glutamic acid residue in the 12th and last position of loop IV with a glutamine are studied by NMR spectroscopy. The mutation E140Q results in sequential Ca2+ binding in this domain and has far-reaching effects on the structure of (Ca2+)2 TR2C, thereby providing further evidence for the critical role of this glutamic acid residue for the Ca2+-induced conformational change of regulatory EF-hand proteins. Analyses of the NOESY spectra of the mutant under Ca2+-saturated conditions, such that 97% of the protein is in the (Ca2+)2 form, revealed two sets of mutually exclusive NOEs. One set of NOEs is found to be consistent with the closed structure observed in the apo state of the C-terminal domain of the wild-type protein, while the other set supports the open structure observed in the Ca2+-saturated state. In addition, several residues in the hydrophobic core exhibit broadened resonances. We conclude that the (Ca2+)2 form of the mutant experiences a global conformational exchange between states similar to the closed and open conformations of the C-terminal domain of wild-type calmodulin. A population of 65 +/- 15% of the open conformation and an exchange rate of (1-7) x 10(4) s(-1) were estimated from the NMR data and the chemical shifts of the wild-type protein. From a Ca2+ titration of the 15N-labeled mutant, the macroscopic binding constants [log(K1) = 4.9 +/- 0.3 and log(K2) = 3.15 +/- 0.10] and the inherent chemical shifts of the intermediate (Ca2+)1 form of the mutant were determined using NMR. Valuable information was also provided on the mechanism of the Ca2+ activation and the roles of the structural elements in the two Ca2+-binding events. Comparison with the wild-type protein indicates that the (Ca2+)1 conformation of the mutant is essentially closed but that some rearrangement of the empty loop IV toward the Ca2+-bound form has occurred.

Sequence dependence and direct measurement of crossover isomer distribution in model Holliday junctions using NMR spectroscopy

A 32-base-pair model of the Holliday junction (HJ) intermediate in genetic recombination has been prepared and analyzed in-depth by 2D and 3D (1)H NMR spectroscopy. This HJ (J2P1) corresponds to a cyclic permutation of the base pairs at the junction relative to a previously studied HJ [J2; Chen, S.-M., & Chazin, W.J. (1994) Biochemistry 33, 11453-11459], designed to probe the effect of the sequence at the n - 1 position (where n is the residue directly at the branch point) on the stacking geometry. Observation of several interbase nuclear Overhauser effects (NOEs) clearly indicates a strong preference for the isomer opposite that observed for J2, confirming the dependence of stacking isomer preference on the sequence at the junction. As for other model HJs studied, a small equilibrium distribution of the alternate isomer could be identified. A sample of J2P1 was prepared with a single (15)N-labeled thymine residue at the branch point. 1D (15)n-filtered (1)H-detected experiments on this sample at low temperature give strong support for the co-existence of the two stacking isomers and provide a much more direct and accurate measure of the crossover isomer distribution. The comparative analysis of our immobile HJs and a model cruciform structure [Pikkemaat, J.A., van den Elst, H., van Boom, J.H., & Altona, C. (1994) Biochemistry 33, 14896-14907] sheds new light on the issue of the relevance of crossover isomer preference in vivo.

Slow conformational exchange in DNA minihairpin loops: a conformational study of the circular dumbbell d<pCGC-TT-GCG-TT>

In recent years various examples of highly stable two-residue hairpin loops (miniloops) in DNA have been encountered. As the detailed structure and stability of miniloops appear to be determined not only by the nature and sequence of the two bases in the loop, but also by the closing base pair, it is desirable to carry out in-depth studies of especially designed small model DNA compounds. Therefore, a circular DNA dumbbell-like molecule is tailored to consist of a stem of three Watson-Crick base pairs, flanked on each side by a minihairpin loop. The resulting circular DNA decamer 5'-d<pCGC-TT-GCG-TT>-3' (I) is studied in solution by means of nmr spectroscopy. At a temperature of 269 K the molecule occurs in a 50/50 mixture of two dumbbell structures (denoted L2L2 and L2L4). L2L2 contains three Watson-Crick C-G base pairs and two two-residue loops (H2-family type) in opposite parts of the molecule. On raising the temperature from 269 to 314 K, the L2L4 conformer becomes increasingly dominant (95% at 314 K). This conformer has a partially disrupted closing G-C base pair in the 5'-GTTC-3' loop with only one remaining solvent-accessible hydrogen bond between NH alpha of the cytosine C(1) and O6 of the guanine G(8), whereas the opposite 5'-CTTG-3' loop remains stable. The disruption of the C(1)-G(8) base pair in the L2L4 form is correlated with the presence of a syn orientation for the C(1) base at the 5'-3' loop-stem junction in the 5'-GTTC-3' loop. The two conformers, L2L2 and L2L4, occur in slow equilibrium (2-20 s-1). Moderate line broadening of specific 1H, 13C, and 31P resonances of residues C(1), G(8), T(9), and T(10) at low temperatures, due to chemical exchange between L2L2 and L2L4, show that the interconversion from an anti to syn conformer in residue C(1) has a small local effect on the structure of the dumbbell. T1 relaxation measurements, chemical-shift considerations, and complete band-shape calculations of the exchange process of the G(8) imino proton reveal a possibility for the existence of multiconformational states in the anti-syn equilibrium.

Thermodynamics of melting of the circular dumbbell d<pCGC-TT-GCG-TT>

The conformational behavior of DNA minihairpin loops is sensitive to the directionality of the base pair that closes the loop. Especially tailored circular dumbbells, consisting of a stem of three Watson-Crick base pairs capped on each side with a minihairpin loop, serve as excellent model compounds by means of which deeper insight is gained into the relative stability and melting properties of hairpin loops that differ only in directionality of the closing pair: C-G vs G-C. For this reason the thermodynamic properties of the circular DNA decamers 5'-d<pCGC-TT-GCG-TT>-3' (I) and reference compounds 5'-d<pGGC-TT-GCC-TT>-3' (II) and 5'-d(GCG-TC-CGC)-3' (III) are studied by means of nmr spectroscopy. Molecules I and II adopt dumbbell structures closed on both sides by a two-membered hairpin loop. At low temperature I consists of a mixture of two slowly exchanging forms, denoted L2L2 and L2L4. The low-temperature L2L2 form is the fully intact minihairpin structure with three Watson-Crick C-G base pairs. The high-temperature form, L2L4, contains a partially disrupted closing G-C base pair in the 5'-GTTC-3' loop, with the cytosine base placed in a syn orientation. The opposite 5'-CTTG-3' loop remains stable. A study of the noncircular hairpin structure III shows similar conformational behavior for the 5'-GTTC-3' loop as found in I; a syn orientation for C(6) and two slowly exchanging imino proton signals for G(3). The melting point Tm of II was estimated to lie above 365 K. The Tm value of the duplex stem and the 5'-CTTG-3' loop of the L2L4 form of I is 352 +/- 2 K. The delta H0 is calculated as -89 +/- 10 kJ/mol. The Tm value determined for the individual residues of the 5'-GTTC-3' loop lies 4 degrees-11 degrees lower. The enthalpy delta H0 of melting the thymine residues in the 5'-GTTC-3' loop is calculated to be -61 +/- 7 kJ/mol. Thermodynamic data of the equilibrium between the slowly exchanging two- and four-membered loop conformers of I reveal an upper limit for delta H0 of +30 kJ/mol in going from a two-membered to a four-membered loop, in agreement with the enthalpy difference of +28 kJ/mol between the two loops at the Tm midpoint. For hairpin III the upper limit for delta H0 in going from a two-membered to a four-membered loop amounts to +/- 21 kJ/mol.(ABSTRACT TRUNCATED AT 400 WORDS)

1H, 15N, and 13C backbone chemical shift assignments, secondary structure, and magnesium-binding characteristics of the Bacillus subtilis response regulator, Spo0F, determined by heteronuclear high-resolution NMR

Spo0F, sporulation stage 0 F protein, a 124-residue protein responsible, in part, for regulating the transition of Bacillus subtilis from a vegetative state to a dormant endospore, has been studied by high-resolution NMR. The 1H, 15N, and 13C chemical shift assignments for the backbone residues have been determined from analyses of 3D spectra, 15N TOCSY-HSQC, 15N NOESY-HSQC, HNCA, and HN(CO)CA. Assignments for many sidechain proton resonances are also reported. The secondary structure, inferred from short- and medium-range NOEs, 3JHN alpha coupling constants, and hydrogen exchange patterns, define a topology consistent with a doubly wound (alpha/beta)5 fold. Interestingly, comparison of the secondary structure of Spo0F to the structure of the Escherichia coli response regulator, chemotaxis Y protein (CheY) (Volz K, Matsumura P, 1991, J Biol Chem 266:15511-15519; Bruix M et al., 1993, Eur J Biochem 215:573-585), show differences in the relative length of secondary structure elements that map onto a single face of the tertiary structure of CheY. This surface may define a region of binding specificity for response regulators. Magnesium titration of Spo0F, followed by amide chemical shift changes, gives an equilibrium dissociation constant of 20 +/- 5 mM. Amide resonances most perturbed by magnesium binding are near the putative site of phosphorylation, Asp 54.

Thermodynamic and structural properties of r(ACC) as revealed by ultraviolet electronic absorption, circular dichroism, 1H-NMR spectroscopy and Monte Carlo simulations

UV absorption, circular dichroism (CD) and 1H NMR, associated with Monte Carlo (MC) molecular structure simulations have been applied to the study of the trinucleoside diphosphate: r(ACC). The MC study which has been conducted as a function of temperature, is based on random variations of the nucleotide conformational angles, i.e. phosphodiester chain torsional angles and sugar pucker pseudorotational angles. All of the chemical bond lengths and valence angles remained fixed during the structural simulation, except those of the sugar pucker. Six different initial structures have been selected in order to explore the molecular conformational space as completely as possible. This simulation procedure led to distinct families of equilibrium conformations at 283, 298 and 318 K. The thermodynamical parameters such as variations in entropy, enthalpy and also melting temperature (delta SX0, delta HX0 and Tm) of the stacking (X) equilibrium were obtained from UV absorption and circular dichroism (CD) spectra recorded over a 80K temperature range. Chemical shifts (delta), vicinal coupling constants (3Jk,l), and cross-relaxation rate (sigma k,l) of trimers were measured at 400.13 MHz over a range of concentrations (2-13 mM) and temperatures (283-333K). Least-squares fitting of the experimental chemical shifts to simple models of association (A) and stacking equilibria allowed separation of the variations in the delta values (delta delta X and delta delta A) due to either phenomenon. The three NMR data sets (delta delta X, 3Jk,l, and sigma k,l) were then evaluated for the minima conformers obtained with the MC stimulations. Theoretical values of delta delta X were estimated using the results of an ab initio study while the coupling constant data were simulated with Karplus-type equations. Finally, the relaxation data were simulated from the distance matrices using treatment for cases of both slow conformational exchange accompanied by rapid small-amplitude fluctuations about the minima structures. A consistent picture of the large amplitude deformations (torsional angle variation) of these trimers has emerged from the present study. Optimized conformational blends at 283,296 and 318K were obtained by least-squares fitting of the experimental data to the theoretical ones, while considering the populations as adjustable parameters. As it would be expected, the right-handed helical conformation (A-RNA type) is found to be the major stacked species, in the temperature range of 283 to 318K. Limited evidence for bulged structures has been obtained, whereas novel reverse-stacked and half-stacked conformers also presented theoretical data compatible with the NMR observables of aqueous r(ACC).

The role of lysine 99 of Thiobacillus versutus cytochrome c-550 in the alkaline transition

The methionine ligand of the heme iron in ferricytochrome c-550 from Thiobacillus versutus is replaced by another residue at high pH. This transition is similar to the alkaline transition in mitochondrial cytochrome c. To investigate the possible role of lysine 99 in this process, this residue has been mutated to a glutamate. The mutation causes the apparent pKa of the transition to decrease from 11.2 in wild type to 10.8 in Lys99Glu cytochrome c-550. This destabilization of the native form is ascribed to the absence of the hydrogen bond between the epsilon-amine group of Lys99 and the carbonyl of Lys54 in the mutant protein. The 1H-NMR spectrum of Lys99Glu ferricytochrome c-550 at alkaline pH still shows resonance positions of the heme methyl peaks that are characteristic of the alkaline form. These results strongly suggest that Lys99 does not act as a ligand in the high-pH form, contrary to the case of yeast iso-1-cytochrome c. Evidence has been presented that in the latter protein the homologous Lys79 can act as a ligand in the alkaline form [1993, J. Am. Chem. Soc. 115, 7507-7508]. In the EPR spectrum of Lys99Glu cytochrome c-550 the species with Met-His coordination (gz = 3.27) is replaced by two forms with gz = 3.45 and 3.20 in the alkaline form (pH > or = 10.6). At pH > 11 yet another form is observed with g-values 2.87, 2.18 and 1.60, tentatively identified as a species with a lysine-histidinate coordination of the heme iron.

Characterization of mutant Met100Lys of cytochrome c-550 from Thiobacillus versutus with lysine-histidine heme ligation

The heme iron in cytochrome c-550 from Thiobacillus versutus has a methionine and a histidine as axial ligands. In order to study the characteristics of a possible lysine-histidine ligation in a heme protein, the methionine has been replaced by a lysine. This residue acts as a ligand between pH 3 and 12. The midpoint potential of the mutant has shifted -329 mV compared to wild type, but apart from this shift the pH dependence of the midpoint potential is unchanged, suggesting that the large drop is caused by specific ligand effects and not by protein refolding. While the EPR spectrum of wild-type cytochrome c-550 shows one species with gz = 3.35, in the spectrum of the mutant two species occur with gz values of 3.53 and 3.30. The intensity ratio of both species depends on the presence of organic cosolvents. In the low frequency region (-4 to -1 ppm) of the 1H NMR spectrum of mutant ferrocytochrome c-550, four one-proton peaks replace the resonances of the ligand methionine side chain protons. Using two-dimensional NMR spectroscopy (COSY and NOESY), these protons and five others have been assigned to the lysine ligand. The spectroscopic results obtained for this mutant show similarities with those observed for the alkaline form of cytochrome c, supporting the Lys-His ligation proposed for this protein. The data are consistent with the evidence for amine ligation in cytochrome f: the EPR spectrum of M100K cytc-550 is similar to that of cytochrome f. However, the NMR spectra show significant differences.(ABSTRACT TRUNCATED AT 250 WORDS)

Solution structure of lariat RNA by 500 MHz NMR spectroscopy and molecular dynamics studies in water

A 500 MHz NMR study of the lariat RNA tetramer 1 and pentamer 2 mimicking the naturally occurring lariat RNA is reported. The conformational properties of 1 and 2 were compared with those of a linear branched RNA tetramer 3, which show that the conformational features of the two lariat RNAs, 1 and 2, are quite constrained and significantly different from those observed for the linear branched RNA tetramer 3. The conformation of all sugar residues forming the lariat ring in 1 and 2 are locked in a rigid South-type conformation. All residues in both lariat RNAs have a high population of gamma+ (67-85%) and beta t (95-100%) rotamers except guanosine where the gamma+ population is low. The conformation around the glycosidic bond is anti for all residues except for guanosine where NOE data indicates an equilibrium of syn<-->anti. In both lariat RNAs, 1 and 2, the temperature dependent 1H and 31P chemical shifts as well as the oligomerization shifts, with respect to adenosine 2',3',5'-triethyl-phosphate (Sund et al., 1992, Tetrahedron 48, 695) suggests that the 3'-->5' linked U4 or C4 residue is stacked on guanosine. Subsequently, 1H-1H, 1H-31P and 13C-31P coupling constants derived torsional constraints were used for molecular dynamics study in water with counter sodium ions for a total of 226 ps. The MD simulations were first carried out with harmonic torsional constraints which were derived from J couplings (0-86 ps) and then completely without constraints (96-226 ps). The lack of any major changes in the conformation of the two lariat-RNA structures upon releasing the NMR constraints indicate that the conformers generated in the MD simulation in water agree well with the structural features suggested by experimental observables. This means that the ensemble of conformers generated during the MD trajectory of 226 ps are not artificially held in these conformations due to the NMR constraints, suggesting that these conformers can be considered to be good representatives of the actual NMR observed solution structures.

An NMR study of the conformation and thermodynamics of the circular dumbbell d [formula: see text] Slow exchange between two- and four-membered hairpin loops

The circular DNA decamer 5'-d [formula: see text] 3' is studied in solution by means of NMR spectroscopy. At low temperature the molecule adopts a dumbbell structure with three Watson-Crick C-G base pairs and two two-residue loops in opposite parts of the molecule. On raising the temperature another conformer appears, in which the closing C-G base pair in the 5'-GTTC-3' loop is disrupted, whereas the opposite 5'-CTTG-3' loop remains stable. The two conformers are in slow equilibrium over a limited temperature range.

Assignment of the 600-MHz 1H-NMR spectrum of amicyanin from Thiobacillus versutus by two-dimensional NMR methods provides information on secondary structure

The nearly complete assignment (pH 6.8; T 310 K) of the 1H-NMR spectrum of reduced amicyanin from Thiobacillus versutus is reported. Experimental evidence is presented, that the structure of the amicyanin contains two beta-sheets, a feature common to plastocyanins and azurins. The loops joining the beta-strands have also been identified. The loop F-G (Thr94-Phe98), together with the flanking residues Cys93 and Met99, comprises three of the four copper ligands and is short compared to similar loops in plastocyanin and azurin. His96 turns out to be the copper ligand that can be protonated. Amicyanin resembles plastocyanin in overall structure but differs from it on account of an N-terminal strand of 22 amino acids in front of strand A, shorter loops A-B, D-E and F-G and the absence of any alpha-helical segments.

Isolation and characterization of cytochrome c550 from the methylamine-oxidizing electron-transport chain of Thiobacillus versutus

The isolation and purification of cytochrome c550 from the methylamine-oxidizing electron-transport chain in Thiobacillus versutus is reported. The cytochrome is a single-heme-containing type I cytochrome c with a relative molecular mass of 16 +/- 1 kDa, an isoelectric point of 4.6 +/- 0.1, a midpoint potential of 272 +/- 3 mV at pH less than 4 and 255 +/- 5 mV at pH = 7.0, and an axial coordination of the Fe by a methionine and a histidine. The midpoint potential decreases with increasing pH due to the deprotonation of a group tentatively identified as a propionate (pKa = 6.5 +/- 0.1 and 6.7 +/- 0.1 in the oxidized and reduced protein, respectively) and a change in the Fe coordination at pH greater than 10. The electron-self-exchange rate appears to depend strongly on the ionic strength of the solution and is relatively insensitive to changes in pH. At 313 K and pH 5.2 the electron-exchange rate amounts to 0.7 x 10(2) M-1 s-1 and 5.3 x 10(2) M-1 s-1 at I = 40 mM and I = 200 mM, respectively. Amino acid composition and molar absorption coefficients at various wavelengths are reported. Resonances of heme protons and the epsilon H3 group of the ligand methionine of the Fe have been identified in the 1H-NMR spectrum of the reduced as well as the oxidized cytochrome.

Conformational consequences of the incorporation of arabinofuranosylcytidine in DNA. An NMR study of the DNA fragments d(CGCTAGCG) and d(CGaCTAGCG) in solution

The self-complementary octamers d(CGCTAGCG) and d(CGaCTAGCG) (aC, arabinofuranosylcytidine) were studied by means of NMR spectroscopy. It is shown that d(CGaCTAGCG), under suitable conditions of oligonucleotide concentration, ionic strength and temperature, exclusively adopts a hairpin structure. However, under the same experimental conditions (5 mM DNA, no added salt, 295 K) d(CGCTAGCG) mainly adopts a B-DNA-type duplex. At lower temperatures (less than or equal to 290 K) the hairpin form of d(CGaCTAGCG) occurs in slow exchange with an intact B-DNA-type duplex. When the DNA concentration of d(CGCTAGCG) is dramatically reduced (less than or equal to 0.5 mM) the hairpin form becomes highly favoured at the expense of the dimer. Moreover, proton-chemical-shift considerations indicate that the structural features of the hairpin structure of d(CGCTAGCG) mimic, in part, those of the modified octamer d(CGaCTAGCG), i.e. a loop comprising only the two central residues with the thymine located into the minor groove (Pieters, J. M. L., de Vroom, E., van der Marel, G. A., van Boom, J. H., Koning, T. M. G., Kaptein, R. and Altona, C. unpublished results). Thermodynamic analysis of d(CGCTAGCG) yields an average Tmd value of 342 K (1 M DNA) and a delta Hod value of -266 kJ/mol for the dimer/coil transition and an average Tmh value of 321 K and delta Hoh - 102 kJ/mol for the hairpin/coil equilibrium. For the duplex/coil equilibrium of d(CGaCTAGCG) an average Tmd value of 336 K (1 M DNA) and delta Hod value of -253 kJ/mol are deduced. The hairpin/coil transition of d(CGaCTAGCG) is characterized by a delta Hoh value of -104 kJ/mol and an average Tmh value of 331 K. It is concluded that incorporation of an arabinofuranosylcytidine in the octamer d(CGaCTAGCG) results in stabilization of the hairpin form, whereas the dimer is destablized by two aC.dG base pairs.

Conformational and thermodynamic consequences of the introduction of a nick in duplexed DNA fragments: an NMR study augmented by biochemical experiments

NMR studies were carried out on various equimolar mixtures consisting of a combination of oligomers: d(ACGGCT) (I). d(pACGGCT) (Ia), d(TGCAGT) (II), d(AGCCGTACTGCA) (III), d(TGCAGTACGGCT) (IV). It is shown that I + II + III (MI) and Ia + II + III (M2) form stable duplexes with nicks in the centre of the respective double helices. A close analysis of the NOESY experiments of M1 and M2 revealed that these fragments form B-DNA type duplex structures. A comparison of the chemical-shift data of the nicked duplexes with those of the intact duplex of III + IV (M3) demonstrated that only small local distortions occur when a nick is introduced. The chemical-shift profiles of M1 and M3 were used to obtain the thermodynamic data for the duplex/coil transitions. The profiles of M1 were analysed by means of a new thermodynamic model (TRIDUP). From the calculated thermodynamic data of M1 and M3 it is concluded that the melting behaviour of M1 occurs cooperatively. A ligation experiment demonstrated that the relatively small substrate (M2) was almost completely joined after an overnight incubation at 14 degrees C.

Thermodynamics of the various forms of the dodecamer d(ATTACCGGTAAT) and of its constituent hexamers from proton nmr chemical shifts and UV melting curves: three-state and four-state thermodynamic models

Chemical shifts of base and H1' protons of the single-stranded hexamers d(ATTACC) and d(GGTAAT), of the 1:1 mixtures of these complementary hexamers, and of the self-complementary dodecamer d(ATTACCGGTAAT) were measured at various temperatures in aqueous solution. Four different sample concentrations were used in the case of the dodecamer and of the mixture of the complementary hexamers; the individual hexamers were measured at two different DNA concentrations. Absorbance temperature profiles at five different NaCl concentrations were measured for the dodecamer in order to quantify the effect of the ionic strength on the duplex formation. Under suitable conditions of nucleotide concentration, temperature, and ionic strength, the dodecamer adopts either a B-DNA duplex or a hairpin-loop structure. Chemical shift vs temperature profiles, constructed for all samples, were used to obtain thermodynamic parameters either for the various stacking interactions in the single strands or for the duplex or the hairpin-loop formation. In the analysis of the duplex formation of the hexamers, a two-state approach appeared too simple, because systematic deviations were revealed. Therefore, a new three-state model (DUPSTAK) was developed. In order to investigate the magnitude of error arising from the use of the two-state approach in cases where the DUPSTAK model appears more appropriate, a series of test calculations was made. The magnitude of error in the enthalpy and in the entropy of duplex melting is found to depend linearly upon the actual melting temperature and not upon the individual delta Hd degrees and delta Sd degrees values. Thermodynamic analysis of the chemical shift vs temperature profiles in D2O solution (no added salt) yields an average Tmd value of 341 K (1M DNA) and delta Hd degrees of - 121 kJ.mol-1 for the dimer/random-coil transition of the hexamer duplex d(ATTACC).d(GGTAAT). For the duplex in equilibrium random-coil transition of the 12-mer d(ATTACCGGTAAT) an average Tmd value of 336 K (1M DNA) and delta Hd degrees of -372 kJ.mol-1 are found. The hairpin/random-coil transition of d(ATTACCGGTAAT) is characterized by a rather large delta Hh degrees value, -130 kJ.mol-1, and an average Tmh value of 304 K.

Alterations in the d(CpGpT) structure in solution as a result of [PtCl(diethylenetriamine)]+ binding

The trinucleotide d(CpGpT) reacts with [PtCl(dien)]Cl (dien = diethylenetriamine) to yield as a single adduct Pt(dien)[d(CpGpT)-N7(2)]. The structure of this adduct in solution has been analysed with the aid of NMR spectroscopy and compared with that of the unmodified trinucleotide. A change in the population of the S conformer of the guanosine deoxyribose ring and a syn preference of the guanine residue are the most important changes occurring upon platination. As a result the dC-dG stack disappears, whereas the dG-dT stack is hardly affected. The CD spectra of both platinated and free d(CpGpT) confirm the different nature of the two molecules.

Influence of uracil on the conformational behaviour of RNA oligonucleotides in solution

NMR studies were carried out on some alternating pyrimidine-purine sequences: the single-stranded tetramers CACA and UGUG and the self-complementary octamer CACAUGUG. Assignments, based upon COSY, homonuclear Hartmann-Hahn, and NOESY experiments, are given for the resonances of all base protons and of several sugar protons. Chemical shift vs temperature profiles were used to obtain thermodynamic parameters for the single-stranded stack in equilibrium with random coil and the duplex in equilibrium with random coil equilibria. The populations of N-type conformer of the ribose rings were estimated from the observed J1'2'. Comparisons with another alternating pyrimidine-purine sequence Um2(6)AUm2(6)A and with the deoxyribose counterparts d(CACA), d(TGTG) and d(CACATGTG) are given. Previous 1H-NMR investigations of Um2(6)AUm2(6)A revealed that the population of bulge-out structure diminishes compared to m2(6)AUm2(6)A due to the U(1)-m2(6)A(2) stacking interaction. In CACA a strong stacking proclivity (Tm = 310 K) together with a clear preference for N-type ribose is observed. However, the stacking interactions in UGUG are relatively less stable (Tm = 288 K) and a bias towards S-type sugar is present. Besides a small amount of stack, a significant contribution of bulge out structure is proposed for UGUG. We conclude that the nature of the pyrimidine base mainly determines the formation of bulge-out structures. The poor stacking properties of uracil now appear to be mainly responsible for this phenomenon. Comparison with the deoxyribose counterparts shows a reasonable agreement between the Tm values of CACA and d(CACA), whereas the Tm of UGUG (288 K) is much lower than the Tm of d(TGTG) (315 K). It is suggested that the absence of bulge-out structures in DNA purine-pyrimidine-purine sequences is related to the relatively strong stacking proclivity of dT residues compared to that of U residues. The Tm values (average 341 K) for the duplex in equilibrium with random coil transition obtained for each residue of CACAUGUG appear very similar. All ribose rings, except the G(8), adopt a pure N conformer in the duplex. This is taken to mean that the differences in conformational behaviour of the constituent tetramers disappear upon duplex formation.

Extra thymidine stacks into the d(CTGGTGCGG).d(CCGCCCAG) duplex. An NMR and model-building study

NMR and model-building studies were carried out on the duplex d(CTGGTGCGG).d(CCGCCCAG), referred to as (9+8)-mer, which contains an unpaired thymidine residue. Resonances of the base and of several sugar protons of the (9+8)-mer were assigned by means of a NOESY experiment. Interresidue NOEs between dG(4) and dT(5) as well as between dT(5) and dG(6) provided evidence that the extra dT is stacked into the duplex. Thermodynamic analysis of the chemical shift vs temperature profiles yielded an average TmD value of 334 K and delta HD of -289 kJmol-1 for the duplex in equilibrium random-coil transition. The shapes of the shift profiles as well as the thermodynamic parameters obtained for the extra dT residue and its neighbours again indicate that the unpaired dT base is incorporated inside an otherwise intact duplex. This conclusion is further supported by (a) the observation of an imino-proton resonance of the unpaired dT; (b) the relatively small dispersion in 31P chemical shifts (approximately 0.5 ppm) for the (9+8)-mer, which indicates the absence of t/g or g/t combinations for the phosphate diester torsion angles alpha/zeta. An energy-minimized model of the (9+8)-mer, which fits the present collection of experimental data, is presented.

Conformational analysis of the tetranucleotides m6(2)A-m6(2)A-U-m6(2)A(m6(2)A = N6-dimethyladenosine) and U-m6(2)A-U-m6(2)A and of the hybrid dA-r(U-A). A one- and two-dimensional NMR study

A 1H-NMR investigation was carried out on the tetranucleotides U-m6(2)A-U-m6(2)A and m6(2)A-m6(2)A-U-m6(2)A (m6(2) = N6-dimethyladenosine) as well as on the hybrid trinucleotide dA-r(U-A). An extensive comparison with m6(2)A-U-m6(2)A and other relevant compounds is made. Previous proton NMR studies on trinucleotides have shown that purine-pyrimidine-purine sequences prefer to adopt a mixture of states which have as a common feature that the interior pyrimidine residue bulges out, whereas the flanking purine residues stack upon each other. A stacking interaction on the 3' side of the bulge is known to have no measurable effect on the bulge population. Chemical-shift data, ribose ring conformational analysis and information from NOE experiments now show unambiguously that the moderate U(1)-m6(2)A(2) stack in U-m6(2)A-U-m6(2)A diminishes the population of bulged-out structures in favour of a regular stack. This tendency towards conformational transmission in the downstream 5'----3' direction is fully confirmed by the fact that the strong m6(2)A(1)-m6(2)A(2) stack in the tetranucleotide m6(2)A-m6(2)A-U-m6(2)A virtually precludes the formation of bulged-out structures. The conformational characteristics of dA-r(U-A) appear comparable with those of m6(2)A-U-m6(2)A, which indicates that the presence of a 2'-hydroxyl group in the first purine residue is not a necessary prerequisite for the formation of a bulge.

Conformational analysis of the dinucleotides 5'-methylphospho-N6-dimethyladenylyl-uridine (mpm62A-U) and 5'methylphospho-uridylyl-N6-dimethyladenosine (mpU-m62A) and of the trinucleotide U-m62A-U. A nuclear magnetic resonance and circular dichroic study

NMR and CD studies were carried out on the dinucleotides 5'-methylphospho-N6-dimethyladenylyl-uridine (mpm62-U) and 5'-methylphospho-uridylyl-N6-dimethyladenosine (mpU-m62A) and on the trinucleotide U-m62A-U. A detailed comparison is given of the conformational features of mpm62A-U and mpU-m62A with the corresponding 5'-nonphosphorylated dinucleotides m62A-U and U-m62A, respectively. The behaviour of the trinucleotide U-m62A-U is compared with the properties of the constituent dinucleotides U-m62A and mpm62A-U. Chemical-shift and CD data were used to determine the amount of stacking interactions. For each compound NMR spectra were recorded at two or three sample concentrations in order to separate intermolecular and intramolecular base-base interactions. The coupling constants of the ribose ring are interpreted in terms of the N/S equilibrium, and population distributions along the backbone angles beta, gamma and epsilon are presented. The combined data indicate a strong similarity between mpm62A-U and m62A-U both in degree and in mode of stacking. In contrast, the existence of different types of stacking interactions in mpU-m62A and U-m62A is suggested in order to explain the NMR and CD data. It is concluded that dinucleoside bisphosphates serve better as a model for the behaviour of trinucleotides than dinucleoside monophosphates. The trinucleotide U-m62A-U adopts a regular single-stranded stacked RNA structure with preference for N-type ribose and gamma+ and beta t backbone torsion angles. The difference in behaviour between the U-m62A- part of U-m62A-U and the dimer U-m62A is seen as a typical example of conformational transmission.

An NMR study of polymorphous behaviour of the mismatched DNA octamer d(m5C-G-m5C-G-A-G-m5C-G) in solution. The B-duplex and hairpin forms

By means of one- and two-dimensional NMR spectroscopy the solution structures of the partly self-complementary octamer d(m5C-G-m5C-G-A-G-m5C-G) were investigated. It is shown that this DNA fragment, under conditions of high DNA concentration (8 mM DNA) and/or high ionic strength prefers to adopt a duplex structure. At low DNA concentration (0.4 mM DNA), the duplex exists in a 1:1 slow equilibrium with a monomeric hairpin form. Addition of salt destabilizes the hairpin structure in favour of the dimer. At high temperatures the hairpin form, as well as the dimer structure, exist in a fast equilibrium with the random-coil form. For the hairpin/random-coil equilibrium a Tm of 329 K and a delta H degree of -121 kJ.mol-1 were deduced. These thermodynamic parameters are independent of the DNA concentration, as is expected for a monomeric structure. For the dimer to coil transition a Tm of 359 K (1 M DNA) and a delta H degree of -285 kJ.mol duplex-1 were derived. The thermodynamic data of the hairpin-coil transition mutually agree with those recently reported for the hairpin to random coil equilibrium of the DNA octamer d(m5C-G-m5C-G-T-G-m5C-G) [Orbons, L. P. M., van der Marel, G. A., van Boom, J. H. & Altona, C. (1987) J. Biomol. Struct. Dyns. 4, 939-963]. It is demonstrated that the dimer structure exhibits B-DNA characteristics, as is witnessed by the NOESY experiments and the analysis of the proton-proton coupling data. It is shown that the base-pair formation of the G x A mismatches is anti-anti. A comparison of 1H and 31P chemical-shift data of the title compound with those of a well-characterized B-DNA structure reveals large differences in the dm5C(3)-dG(4)-dA(5) part of the mismatched dimer structure. These differences apparently indicate some major local structural changes due to the incorporation of the G x A mismatches. Under the most extreme conditions used (i.e. up to 3 M NaCl or 75% CH3OH in the presence of 10 mM MgCl2) no Z-DNA structure was observed. It is shown that the structural features of the hairpin form of the title compound mimic those of the hairpin structure of d(m5C-G-m5C-G-T-G-m5C-G). An energy-minimized model of the hairpin form is given.

Influence of the base sequence on the conformational behaviour of DNA polynucleotides in solution

NMR studies were carried out on samples of the non-self-complementary tetramers d(C-A-C-A), d(T-G-T-G), d(G-A-G-A) and d(T-C-T-C) and of 1:1 mixtures of the complementary tetramers d(C-A-C-A) X d(T-G-T-G) and d(G-A-G-A) X d(T-C-T-C) at two DNA concentrations and of the self-complementary octamers d(C-A-C-A-T-G-T-G) and d(G-A-G-A-T-C-T-C). Assignments, based upon one-dimensional NOE and homonuclear-decoupling and two-dimensional correlated and NOE spectroscopies are given of the resonances of most of the base and sugar protons. Chemical shift vs temperature profiles, constructed for all samples, yielded insight into the temperature- and concentration-dependent conformational behaviour of the compounds and were used to obtain thermodynamic parameters pertaining to the stacked-single-strand----random-coil and duplex----random-coil equilibria. Vicinal proton-proton couplings were analyzed in terms of the conformation of the deoxyribose rings in the single-stranded tetramers and duplexed octamers. The NOE patterns, chemical shift profiles, imino-proton resonances and coupling data revealed that the compounds adopt B-DNA-like structures. The ratio duplexed/stacked-single-strand/random coil depends upon external conditions as well as upon base sequence. The thermodynamic data indicate that: in terms of single-helical stacking, the R-R steps (Tm 321-328 K) appear more stable than the Y-R or R-Y steps (Tm 308-316 K) and the Y-Y steps score least (Tm 290-300 K), and the duplexes consisting of alternating, d(Y-R)n, strands are more stable, in terms of delta H degrees, compared to the d(R-R)n X d(Y-Y)n duplexes. The analyses of the couplings demonstrated that the sugars of the single-stranded tetramers and duplexed octamers occur as a blend of N- and S-type conformers, with a preference for the S-type (C2'-endo) sugar conformation: upon duplex formation, no significant shift in the N-type/S-type ratio was observed. The fraction S-type sugar conformation of a given residue, %S, in the stacked-single strands was found to depend upon the nature of its own base and that of the adjacent residues: sugars in an R-R stretch display high values of %S (90-100), whereas those in Y-Y stretches show relatively low values (approximately equal to 65).(ABSTRACT TRUNCATED AT 400 WORDS)

An NMR study of the polymorphous behavior of the mismatched DNA octamer d(m5C-G-m5C-G-T-G-m5C-G) in solution. The B, Z, and hairpin forms

The polymorphism exhibited by the mismatched octamer d(m5C-G-m5C-G-T-G-m5C-G), as a function of the temperature, DNA concentration and ionic strength, was investigated by means of NMR spectroscopy. It is shown that this partly self-complementary DNA fragment, under conditions of low DNA concentration (0.4 mM) and low ionic strength, exclusively prefers to adopt a monomeric hairpin form, which consists of a stem of three Watson-Crick-type base pairs and a loop of only two residues. This in striking contrast with earlier intimations in literature, which postulated that in oligonucleotides loop formations containing only two residues are sterically impossible. Moreover, the hairpin form displays an unusual stability in comparison with previously reported hairpins. A Tm of 332 K and a delta H degree of -130 kJ.mol-1 were calculated for the hairpin to random coil transition. At high DNA concentration (8 mM) and/or upon the addition of sodium chloride the hairpin form occurs in slow exchange with a B-DNA dimer structure (approximately 20% at 270 K, no added salt), which comprises two central GxT-mismatched base pairs with the bases as major tautomers. At higher ionic strength (greater than 100 mM NaCl), or upon the addition of methanol, a third species appears, which is in slow exchange with both the B dimer and the hairpin form. This third species could be identified with a Z DNA form, comprising two GxT mismatches with the bases as major tautomers, with the guanine bases syn and the cytosine and thymine bases anti.