Antioxidant effect of ethyl pyruvate in respiring neonatal cerebrocortical slices after H(2)O(2) stress

Administration of ethyl pyruvate, which is formed from pyruvate and ethanol, has been found capable of rescuing cells injured by oxidative stress. In one perspective the rescue has been postulated to be metabolic, with the resulting intracellular delivery of pyruvate seen as providing substrate for the TCA Cycle, making it possible to counteract sequela of poly(ADP-ribose)ribosylation, such as depletion of cytosolic NAD(+), glycolytic arrest, and mitochondrial deprivation of pyruvate. The rescue has also been attributed to radical scavenging via the carbonyl groups in ethyl pyruvate and pyruvate. In a previous study we exposed superfused neonatal (P7) brain slices for 60min to 2mM H(2)O(2) and found evidence for both rescue mechanisms. To see if ethyl pyruvate's actions stemmed more from being an antioxidant than from being a nutrient we conducted six new experiments using the same H(2)O(2) protocol, but with two new rescue solutions: [10mM] glucose (glc) plus one of the following: ethyl pyruvate [20mM], or the nonmetabolizable radical scavenger N-tert-butyl-alpha-phenylnitrone (PBN, 1mM). Final ATP values compared to initial, measured in 14.1T (31)P NMR spectra of PCA extracts, were the same: 0.70+/-0.08 for the former (N=3), and 0.64+/-0.08 for the latter (N=3). Quantifications of this study's (1)H NMR metabolites, also measured at 14.1T, exhibited separate clustering when pooled with data from the previous study and compared in a metabolomic multivariate analyses. Because the addition of ethyl pyruvate provided the same ATP protection as the addition of a nonmetabolizable antioxidant, antioxidant protection was its prominent protective mechanism in the chosen, high glucose protocol. Having distinct clusters in the Scores Plot of a Partial Least Squares-Discriminant Analysis suggests the feasibility of constructing statistical models that are predictive.

Photoluminescent characterization of atomic diffusion in core-shell nanoparticles

Eu3+ doped LaF3 nanoparticles with core/shell morphologies were synthesized and selected spectroscopic properties were measured as a function of heat treatment times and temperatures. More specifically, the relative intensity of photoluminescence spectra, both through direct excitation of the lanthanide as well as phonon sideband spectra were evaluated with increasing amounts of time held at specific temperatures. A one dimensional approximation was used to compute an effective diffusion coefficient for the rare earth dopants in LaF3. Despite the simplicity of the model employed, the calculated diffusion coefficients based on the spectroscopic results are accurate within an order of magnitude in comparison to other fluoride crystals yielding a simplified approach to estimating kinetic and diffusion effects in optical materials.

NMR determination of oligonucleotide structure

This unit provides an overview of the use of NMR to determine oligonucleotide structure. It covers basic NMR spectral properties, acquisition of interproton distance restraints and torsion angle restraints, structure refinement, assessment of the quality of the structure obtained. Software programs used in the process are also described.

NTP and PCr responses to hypoxia by hypothermic and normothermic respiring, superfused, neonatal rat cerebrocortical slices: an NMR spectroscopy study at 14.1 Tesla

Although mechanisms of hypothermic neuroprotection during oxygen deprivation have long been investigated, further characterizations of various molecular mechanisms are appropriate. Anticipating future studies of hypothermia and hypoxia/ischemia, we investigated the extent to which our ex vivo, NMR-based, superfused brain slice model might be helpful. (Slices are approximately 350 microm thick, with 18 slices per 8 mm NMR tube.) 31P NMR spectroscopic measurements were made of hypothermia-induced changes in high energy phosphates, while simultaneously monitoring and controlling tissue temperature, using 1H NMR, the high spectroscopic resolution available at 14.1 Tesla (600 MHz for protons), and a recently published protocol. NTP and PCr concentrations in healthy, well-oxygenated slices decreased to (55 +/- 15)% and (66 +/- 30)% of their respective values at 28.0 degrees C when warmed to 38.0 degrees C, in approximate agreement with earlier in vivo studies by others. During 30 min hypoxia NTP and PCr decreased to non-observable values, regardless of temperature. After reoxygenation, NTP and PCr recoveries as percentages of respective prehypoxia values were (63% +/- 16%; 70%) +/- 5%) for hypothermic slices (28.0 degrees C), and (46% +/- 13%; 41% +/- hypothermic neuroprotection during oxygen deprivation in this model, which appears suitable for use in further studies.

Akt phosphorylation and cell survival after hypoxia-induced cytochrome c release in superfused respiring neonatal rat cerebrocortical slices

Phosphorylation of Akt before hypoxia (30 min) and during reoxygenation (4 h) was evaluated in superperfused neonatal rat cerebrocortical slices (350 microm, P7, Sprague-Dawley). Cytosolic cytochrome c intensities in Western blots, which were increased at the end of hypoxia. were decreased during reoxygenation. Western blot intensities of phosphorylated Akt (phospho-Akt), nearly undetectable at the end of hypoxia, recovered quickly during reoxygenation, in a trend opposite that for cytochrome c. At 1.5 h and 4 h after hypoxia they became larger or the same as before hypoxia. Total Akt was unchanged by hypoxia and reoxygenation. Phosphocreatine (PCr) and nucleotide triphosphates (NTP) were measured in parallel 14.1 Tesla ex vivo 31P NMR superfused brain slice studies. PCr and alpha-NTP were nearly undetectable at the end of hypoxia. Although they recovered quickly after hypoxia, they were lower than before hypoxia. Reductions in phospho-Akt during hypoxia were consistent with the general unavailability of basic high energy phosphates. Preferential phosphorylation of Akt after hypoxia suggested that substantial reductions in intracellular energy, as indicated by PCr and NTP, might be tolerated by processes important for generating phospho-Akt. Additionally, the post-hypoxia increase in phospho-Akt might have contributed to concomitant reductions in cytosolic cytochrome c.

The solution structure of the viral binding domain of Tva, the cellular receptor for subgroup A avian leukosis and sarcoma virus

The cellular receptor for subgroup A avian leukosis and sarcoma virus (ALSV-A) is Tva, which contains a motif related to repeats in the low density lipoprotein receptor (LDLR) ligand binding repeat (LBr) and which is necessary for viral entry. As observed with LBr repeats of LDLR, the 47 residue LBr domain of Tva (sTva47) requires calcium during oxidative folding to form the correct disulfide bonds, and calcium enhances the structure of correctly oxidized sTva47, as well as its ability to bind the viral envelope protein (Env). However, solution nuclear magnetic resonance studies indicate that, even in the presence of excess calcium, sTva47 exists in an ensemble of conformations. Nonetheless, as reported here, the structure of the predominant sTva47 solution conformer closely resembles that of other LBr repeats, with identical S-S binding topology and octahedral calcium coordination. The location of W48 and other critical residues on the surface suggests a region of the molecule necessary for Env binding and to mediate post-binding events important for ALSV-A cell entry.

Slow motion in the CAA*TTG sequence of a DNA decamer duplex studied by NMR

In DNA duplexes, pyrimidine-purine steps are believed to be flexible or conformationally unstable. Indeed, several DNA crystal structures exhibit a multitude of conformations for CpA*TpG steps. The question arises of whether this structural flexibility is accompanied by dynamical flexibility, i.e., a question pertaining to the energy barrier between conformations. Except for TpA steps, slow motions on the microsecond-to-millisecond time scale have not been detected in duplexes until now. In the present study, such slow motion was investigated by 1H, 13C, and 15N NMR relaxation measurements on a DNA decamer d(CATTTGCATC)*d(GATGCAAATG). The DNA decamer was enriched with 15% 13C and 98% 15N isotopes for each adenosine and guanosine residue. Three lines of evidence support the notion of slow motion in the CAA*TTG moiety. Analysis of (15)N relaxation showed that the order parameter, S2, of guanosine imino NH groups was about 0.8, similar to that of CH groups for this oligomer. The strong temperature dependence of guanosine NH S2 in the CAA*TTG sequence indicated the presence of a large-amplitude motion. Signals of adenosine H8 protons in the CAA*TTG sequence were broadened in 2D 1H NOESY spectra, which also suggested the existence of slow motion. As well as being smaller than for other adenine residues, the 1H T2 values exhibited a magnetic field strength dependence for all adenosine H8 signals in the ATTTG*CAAAT region, suggesting slow motions more pronounced at the first adenosine in the CAA*TTG sequence but extending over the CAAAT*ATTTG region. This phenomenon was further examined by the pulse field strength dependence of the 1H, 13C, and 15N T1rho values. 1H and 13C T1rho values showed a pulse field strength dependence, but 15N T1rho did not. Assuming a two-site exchange process, an exchange time constant of 20-300 micros was estimated for the first adenosine in the CAA sequence. The exact nature of this motion remains unknown.

Molten-globule structure and membrane binding of the N-terminal protease-resistant domain (63-193) of the steroidogenic acute regulatory protein (StAR)

The first step in steroidogenesis is the movement of cholesterol from the outer to inner mitochondrial membrane; this movement is facilitated by the steroidogenic acute regulatory protein (StAR). StAR has molten-globule properties at low pH and a protease-resistant N-terminal domain at pH 4 and pH 8 comprising residues 63-193. To explore the mechanism of action of StAR we investigated the structural properties of the bacterially expressed N-terminal domain (63-193 StAR) using CD, limited proteolysis and NMR. Far- and near-UV CD showed that the amount of secondary structure was greater at acidic than at neutral pH, but there was little tertiary structure at any pH. Unlike 63-193 StAR liberated from N-62 StAR by proteolysis, biosynthetic 63-193 StAR was no longer resistant to trypsin or proteinase K at pH 7, or to pepsin at pH 4. Addition of trifluoroethanol and SDS increased secondary structure at pH 7, and dodecylphosphocholine and CHAPS increased secondary structure at pH 2, pH 4 and pH 7. However, none of these conditions induced tertiary structure, as monitored by near-UV CD or NMR. Liposomes of phosphatidylcholine, phosphatidylserine and their mixture increased secondary structure of 63-193 StAR at pH 7, as monitored by far-UV CD, and stable protein-liposome complexes were identified by gel-permeation chromatography. These results provide further evidence that the N-terminal domain of StAR is a molten globule, and provide evidence that this conformation facilitates the interaction of the N-terminal domain of StAR with membranes. We suggest that this interaction is the key to understanding the mechanism of StAR's action.

Preservation of the morphological and molecular stability of embryonic tissues

There is presently great interest in using early embryonic tissues, particularly human tissue, for studies of protein and gene expression. Embryonic human tissue is very fragile, and delays often occur before it can be properly prepared for scientific study. Using chick embryos, we have studied the effects of delaying fixation or biochemical isolation on the preservation of cytological characteristics and biochemical molecules. Our study shows that by 60 min post-harvest, tissue morphology and immunofluorescence staining degrades, but the total mRNA profile remains stable. This study suggests that the time between removal of the tissue and fixation is critical to the results and that the critical time is much shorter for embryonic tissues than for more developed tissues. Our results have implications for all research where embryonic tissues are harvested but not processed immediately.

Identification of ligands for RNA targets via structure-based virtual screening: HIV-1 TAR

Binding of the Tat protein to TAR RNA is necessary for viral replication of HIV-1. We screened the Available Chemicals Directory (ACD) to identify ligands to bind to a TAR RNA structure using a four-step docking procedure: rigid docking first, followed by three steps of flexible docking using a pseudobrownian Monte Carlo minimization in torsion angle space with progressively more detailed conformational sampling on a progressively smaller list of top-ranking compounds. To validate the procedure, we successfully docked ligands for five RNA complexes of known structure. For ranking ligands according to binding avidity, an empirical binding free energy function was developed which accounts, in particular, for solvation, isomerization free energy, and changes in conformational entropy. System-specific parameters for the function were derived on a training set of RNA/ligand complexes with known structure and affinity. To validate the free energy function, we screened the entire ACD for ligands for an RNA aptamer which binds L-arginine tightly. The native ligand ranked 17 out of ca. 153,000 compounds screened, i.e., the procedure is able to filter out >99.98% of the database and still retain the native ligand. Screening of the ACD for TAR ligands yielded a high rank for all known TAR ligands contained in the ACD and suggested several other potential TAR ligands. Eight of the highest ranking compounds not previously known to be ligands were assayed for inhibition of the Tat-TAR interaction, and two exhibited a CD50 of ca. 1 microM.

Determination of the populations and structures of multiple conformers in an ensemble from NMR data: multiple-copy refinement of nucleic acid structures using floating weights

A new algorithm is presented for determination of structural conformers and their populations based on NMR data. Restrained Metropolis Monte Carlo simulations or restrained energy minimizations are performed for several copies of a molecule simultaneously. The calculations are restrained with dipolar relaxation rates derived from measured NOE intensities via complete relaxation matrix analysis. The novel feature of the algorithm is that the weights of individual conformers are determined in every refinement step, by the quadratic programming algorithm, in such a way that the restraint energy is minimized. Its design ensures that the calculated populations of the individual conformers are based only on experimental restraints. Presence of internally inconsistent restraints is the driving force for determination of distinct multiple conformers. The method is applied to various simulated test systems. Conformational calculations on nucleic acids are carried out using generalized helical parameters with the program DNAminiCarlo. From different mixtures of A- and B-DNA, minor fractions as low as 10% could be determined with restrained energy minimization. For B-DNA with three local conformers (C2'-endo, O4'-exo, C3'-endo), the minor O4'-exo conformer could not be reliably determined using NOE data typically measured for DNA. The other two conformers, C2'-endo and C3'-endo, could be reproduced by Metropolis Monte Carlo simulated annealing. The behavior of the algorithm in various situations is analyzed, and a number of refinement protocols are discussed. Prior to application of this algorithm to each experimental system, it is suggested that the presence of internal inconsistencies in experimental data be ascertained. In addition, because the performance of the algorithm depends on the type of conformers involved and experimental data available, it is advisable to carry out test calculations with simulated data modeling each experimental system studied.

Structure of the phylogenetically most conserved domain of SRP RNA

The signal recognition particle (SRP) is a phylogenetically conserved ribonucleoprotein required for cotranslational targeting of proteins to the membrane of the endoplasmic reticulum of the bacterial plasma membrane. Domain IV of SRP RNA consists of a short stem-loop structure with two internal loops that contain the most conserved nucleotides of the molecule. All known essential interactions of SRP occur in that moiety containing domain IV. The solution structure of a 43-nt RNA comprising the complete Escherichia coli domain IV was determined by multidimensional NMR and restrained molecular dynamics refinement. Our data confirm the previously determined rigid structure of a smaller subfragment containing the most conserved, symmetric internal loop A (Schmitz et al., Nat Struct Biol, 1999, 6:634-638), where all conserved nucleotides are involved in nucleotide-specific structural interactions. Asymmetric internal loop B provides a hinge in the RNA molecule; it is partially flexible, yet also uniquely structured. The longer strand of internal loop B extends the major groove by creating a ledge-like arrangement; for loop B however, there is no obvious structural role for the conserved nucleotides. The structure of domain IV suggests that loop A is the initial site for the RNA/protein interaction creating specificity, whereas loop B provides a secondary interaction site.

NMR structure of the mature dimer initiation complex of HIV-1 genomic RNA

The two identical genomic RNA strands inside each HIV-1 viral particle are linked through homodimerization of an RNA stem-loop, termed SL1, near their 5' ends. SL1 first dimerizes through a palindromic sequence in its loop, forming a transient kissing-loop complex which then refolds to a mature, linear duplex. We previously reported the NMR structure of a 23-base truncate of SLI in kissing-dimer form, and here report the high-resolution structure of its linear isoform. This structure comprises three short duplex regions--derived from the central palindrome and two stem regions of each strand, respectively--separated by two bulges that each encompass three unpaired adenines flanking the palindromes. The stacking pattern of these adenines differs from that seen in the kissing-loop complex, and leads to greater colinear base stacking overall. Moreover, the mechanical distortion of the palindrome helix is reduced, and base pairs ruptured during formation of the kissing-loop complex are re-established, so that all potential Watson-Crick pairs are intact. These features together likely account for the greater thermodynamic stability of the mature dimer as compared to its kissing-loop precursor.

A Small-Molecule Guanidinium Receptor: The Arginine Cork

A highly preorganized, artificial receptor for guanidinium cations binds arginine with high affinity in polar solvents, even in water! The X-ray crystal structure of its bis(N-ethylguanidinium) complex shows the formation of hydrogen bonds between the receptor and the guest (see picture). This observation explains the strength of the complex formed with the dipeptide diarginine (K(d)=50 µM in H(2)O).

Structure of the most conserved internal loop in SRP RNA

The signal recognition particle (SRP) directs translating ribosomes to the protein translocation apparatus of endoplasmic reticulum (ER) membrane or the bacterial plasma membrane. The SRP is universally conserved, and in prokaryotes consists of two essential subunits, SRP RNA and SRP54, the latter of which binds to signal sequences on the nascent protein chains. Here we describe the solution NMR structure of a 28-mer RNA composing the most conserved part of SRP RNA to which SRP54 binds. Central to this function is a six-nucleotide internal loop that assumes a novel Mg2+-dependent structure with unusual cross-strand interactions; besides a cross-strand A/A stack, two guanines form hydrogen bonds with opposite-strand phosphates. The structure completely explains the phylogenetic conservation of the loop bases, underlining its importance for SRP54 binding and SRP function.

Evaluation of individual and combined neurotoxicity of the immunosuppressants cyclosporine and sirolimus by in vitro multinuclear NMR spectroscopy

Neurotoxicity, a crucial side effect of immunosuppressive therapy with cyclosporine, also has been demonstrated in vitro for sirolimus, a novel macrolide immunosuppressant, which is under clinical investigation in combination with cyclosporine. NMR spectroscopy was used to study the separate and combined effects of cyclosporine and sirolimus on cerebral metabolism, both in brain cells and in perfused rat brain slices. The high-energy phosphate metabolism was already affected significantly at cyclosporine concentrations as low as 100 micrograms/liter: phosphocreatine was reduced by 10 +/- 2% [half-maximal inhibition concentration (IC50) = 1850 +/- 600 micrograms/liter], and nucleoside triphosphate was reduced by 11 +/- 5% (IC50 = 1110 +/- 420 micrograms/liter; n = 4, P <.05). At 500 micrograms/liter cyclosporine, N-acetylaspartate and glutamate were decreased by 13 +/- 7% (IC50 = 1100 +/- 330 micrograms/liter) and 22 +/- 9% (IC50 = 360 +/- 220 micrograms/liter; n = 4, P <.05), respectively. As evaluated using an algorithm based on Loewe isobolograms, combination of cyclosporine and sirolimus resulted in a synergetic reduction of high-energy phosphate metabolites. Addition of sirolimus to the perfusion medium increased brain slice concentrations of cyclosporine. It is concluded that cyclosporine significantly reduced high-energy phosphate metabolism in brain tissue at in vivo relevant concentrations. Combination with sirolimus resulted in synergism, which, in part, is explained by a greater distribution of cyclosporine into the brain tissue in the presence of sirolimus.

NOS inhibitors decrease hypoxia-induced ATP reductions in respiring cerebrocortical slices

BACKGROUND

Excess neuronal nitric oxide (NO) production might cause adenosine triphosphate loss and cellular damage in hypoxic brain parenchyma. 31P nuclear magnetic resonance spectroscopy was used to study hypoxic intracellular responses in perfused respiring cerebrocortical slices, in which NO scavenging by hemoglobin is absent, during NO synthase blockade and NO augmentation.

METHODS

Adenosine triphosphate concentrations were monitored at 4.7 Tesla in respiring slices before, during, and after 60 min of hypoxia (oxygen tension < 5 mmHg). Slices were not treated or were pretreated with 27 microM L-nitroarginine methyl ester (L-NAME), 27 microM 7-nitroindozole (7-NI), or 27 microM L-nitroarginine. Nitrotyrosine:tyrosine ratios of slice extracts were measured using high-performance liquid chromatography. Cresyl violet-stained sections (2 microm) from random slices were examined histologically.

RESULTS

After 60 min of hypoxia, adenosine triphosphate decreased to < or = 3, < or = 3, 65 +/- 6, and 25 +/- 4% of control in slices that were untreated or treated with L-nitroarginine, L-NAME, and 7-NI, respectively. After 120 min of hyperoxic recovery, adenosine triphosphate levels returned to control values in slices pretreated with L-NAME and 7-NI, but to only 30% of control in untreated or L-nitroarginine-treated slices. Nitric oxide donors administered during posthypoxic recovery partially antagonized the adenosine triphosphate recovery found with L-NAME and 7-NI. Nitric oxide synthase activity in slice homogenates, assayed via conversion of L-arginine to citrulline, was < or = 2% of control after all inhibitory treatments. The nitrotyrosine:tyrosine ratio increased by 52% in slices treated with 7-NI and by 200-300% in all other groups. Pretreatment with L-NAME and 7-NI reduced histologic evidence of cell swelling.

CONCLUSION

Neuronal NO is associated with rapid adenosine triphosphate reductions and peroxynitrite formation in acutely hypoxic cerebrocortical slices.

Solution structure of Syrian hamster prion protein rPrP(90-231)

NMR has been used to refine the structure of Syrian hamster (SHa) prion protein rPrP(90-231), which is commensurate with the infectious protease-resistant core of the scrapie prion protein PrPSc. The structure of rPrP(90-231), refolded to resemble the normal cellular isoform PrPC spectroscopically and immunologically, has been studied using multidimensional NMR; initial results were published [James et al. (1997) Proc. Natl. Acad. Sci. U.S.A. 94, 10086-10091]. We now report refinement with better definition revealing important structural and dynamic features which can be related to biological observations pertinent to prion diseases. Structure refinement was based on 2778 unambiguously assigned nuclear Overhauser effect (NOE) connectivities, 297 ambiguous NOE restraints, and 63 scalar coupling constants (3JHNHa). The structure is represented by an ensemble of 25 best-scoring structures from 100 structures calculated using ARIA/X-PLOR and further refined with restrained molecular dynamics using the AMBER 4.1 force field with an explicit shell of water molecules. The rPrP(90-231) structure features a core domain (residues 125-228), with a backbone atomic root-mean-square deviation (RMSD) of 0.67 A, consisting of three alpha-helices (residues 144-154, 172-193, and 200-227) and two short antiparallel beta-strands (residues 129-131 and 161-163). The N-terminus (residues 90-119) is largely unstructured despite some sparse and weak medium-range NOEs implying the existence of bends or turns. The transition region between the core domain and flexible N-terminus, i.e., residues 113-128, consists of hydrophobic residues or glycines and does not adopt any regular secondary structure in aqueous solution. There are about 30 medium- and long-range NOEs within this hydrophobic cluster, so it clearly manifests structure. Multiple discrete conformations are evident, implying the possible existence of one or more metastable states, which may feature in conversion of PrPC to PrPSc. To obtain a more comprehensive picture of rPrP(90-231), dynamics have been studied using amide hydrogen-deuterium exchange and 15N NMR relaxation times (T1 and T2) and 15N{1H} NOE measurements. Comparison of the structure with previous reports suggests sequence-dependent features that may be reflected in a species barrier to prion disease transmission.

Quantitative measurement of transverse and longitudinal cross-correlation between 13C-1H dipolar interaction and 13C chemical shift anisotropy: application to a 13C-labeled DNA duplex

Measurement of both longitudinal and transverse relaxation interference (cross-correlation) between 13C chemical shift anisotropy and 13C-1H dipolar interactions is described. The ratio of the transverse to longitudinal cross-correlation rates readily yields the ratio of spectral densities J(0)/J(omegaC), independent of any structural attributes such as internuclear distance or chemical shift tensor. The spectral density at zero frequency J(0) is also independent of chemical exchange effects. With limited internal motions, the ratio also enables an accurate evaluation of the correlation time for overall molecular tumbling. Applicability of this approach to investigating dynamics has been demonstrated by measurements made at three temperatures using a DNA decamer duplex with purines randomly enriched to 15% in 13C.

Restrained molecular dynamics of solvated duplex DNA using the particle mesh Ewald method

Restrained and unrestrained aqueous solution molecular dynamics simulations applying the particle mesh Ewald (PME) method to DNA duplex structures previously determined via in vacuo restrained molecular dynamics with NMR-derived restraints are reported. Without experimental restraints, the DNA decamer, d(CATTTGCATC).d(GATGCAAATG) and trisdecamer, d(AGCTTGCCTTGAG).d(CTCAAGGCAAGCT), structures are stable on the nanosecond time scale and adopt conformations in the B-DNA family. These free DNA simulations exhibit behavior characteristic of PME simulations previously performed on DNA sequences, including a low helical twist, frequent sugar pucker transitions, BI-BII(epsilon-zeta) transitions and coupled crakshaft (alpha-gamma) motion. Refinement protocols similar to the original in vacuo restrained molecular dynamics (RMD) refinements but in aqueous solution using the Cornell et al. force field [Cornell et al. (1995) J. Am. Chem. Soc., 117, 5179-5197] and a particle mesh Ewald treatment produce structures which fit the restraints very well and are very similar to the original in vacuo NMR structure, except for a significant difference in the average helical twist. Figures of merit for the average structure found in the RMD PME decamer simulations in solution are equivalent to the original in vacuo NMR structure while the figures of merit for the free MD simulations are significantly higher. The free MD simulations with the PME method, however, lead to some sequence-dependent structural features in common with the NMR structures, unlike free MD calculations with earlier force fields and protocols. There is some suggestion that the improved handling of electrostatics by PME improves long-range structural aspects which are not well defined by the short-range nature of NMR restraints.

Structural effect of complete [Rp]-phosphorothioate and phosphorodithioate substitutions in the DNA strand of a model antisense inhibitor-target RNA complex

Chemically modified DNA oligonucleotides have been crucial to the success of antisense therapeutics. Although such modifications are ubiquitous in the clinic, high-resolution structural studies of pharmaceutically relevant derivatives have been limited to only a few molecules. We have completed a high-resolution NMR structural study of three DNA.RNA hybrids with the sequence d(CCTATAATCC). r(GGAUUAUAGG). All hybrids contain an unmodified RNA strand, whereas the DNA strand of each hybrid contains one of three different sugar-phosphate backbone linkages at each nucleotide: (1) phosphate, (2) [Rp]-phosphorothioate, or (3) phosphorodithioate. The UV and NMR melting profiles revealed that the normal hybrid is more stable than the [Rp]-phosphorothioate, which in turn is more stable than the phosphorodithioate. Homonuclear two-dimensional nuclear Overhauser effect spectroscopy and double quantum-filtered correlation spectroscopy afforded nearly complete non-labile proton assignments. The three molecules show nearly equivalent chemical shifts, with the exception of H3' protons, which are shifted downfield in a manner that appears correlated with the degree of sulfur substitution at phosphate. All three hybrids exhibit unusually broad linewidths for deoxyribose protons H2' and H2".Distance restraints were calculated from NOE cross-peak intensities via a complete relaxation matrix approach using the program RANDMARDI. Detailed comparison of interproton distances from each hybrid indicates that the three molecules share a common structure, with neither strand in canonical A or B form. Correlation of R factors, calculated using the program CORMA with DNA H2'-base and H3'-base distances, revealed a relative increase in the population of B-type sugar conformations for deoxyriboses in the A+T-rich center of the hybrid sequence. It is widely known that the activity of enzymes which act upon DNA.RNA hybrid substrates (e.g. ribonuclease H) is impacted when the hybrids contain phosphorothioate or phosphorodithioate substitutions. The structural similarity of the three hybrids examined here suggests that factors other than global structure may mediate the activity of these enzymes.

Effect of cytarabine on the NMR structure of a model okazaki fragment from the SV40 genome

Okazaki fragments occur as intermediates during lagging strand DNA replication. Alterations in Okazaki fragment structure may contribute to the anticancer activities of nucleoside analogues such as cytarabine, a potent anti-leukemic agent that inhibits lagging strand replication. We have determined the solution structures for two model Okazaki fragments, [OKA] and [ARAC]. These sequences are derived from a frequent initiation site for primase during replication of the SV-40 viral genome. The sequence of [ARAC] differs from [OKA] only by substitution of cytarabine for one deoxycytidine. The structure of each model Okazaki fragment was elucidated using NMR spectroscopy and restrained molecular dynamics simulations. The solution structures of [OKA] and [ARAC] each consist of two distinct domains: a DNA duplex region (DDR) and an RNA-DNA hybrid duplex region (HDR). The DDR of [OKA] adopts geometry similar to B-form except for variations in helical parameters, especially twist and roll, which occur in the purine tract, increasing base overlap among the five consecutive purines. The helical axes for the DDR and HDR of [OKA] are bent 22 degrees relative to one another. Although the local structures for the DDR and HDR of [ARAC] are similar to those in [OKA] (root-mean-square deviation (rmsd) approximately 0.8, 1.7 A), the bending at the junction is different (41 degrees for [ARAC] vs 22 degrees for [OKA]). Increased helical bending of cytarabine-substituted Okazaki fragments may contribute to the propensity of cytarabine to inhibit elongation of the lagging strand during DNA replication, and in effecting anticancer activity.

DNA duplex dynamics: NMR relaxation studies of a decamer with uniformly 13C-labeled purine nucleotides

Dynamics in a DNA decamer duplex, d(CATTTGCATC). d(GATGCAAATG), were investigated via a detailed 13C NMR relaxation study. Every 2'-deoxyadenosine and 2'-deoxyguanidine was chemically enriched with 15% 13C and 98% 15N isotopes. Six nuclear relaxation parameters [R(13Cz), R(1Hz), R(2(1)Hz13Cz), R(13Cx), R(2(1)Hz13Cx) and steady-state 13C¿1H¿ NOE] were measured at 600 MHz and three were measured at 500 MHz (1H frequency) for the CH spin systems of sugar 1', 3', and 4' as well as base 8 and 2 positions. A dependence of relaxation parameter values on chemical position was clearly observed; however, no sequence-specific variation was readily evident within our experimental error of approximately 5-10%, except for 3' and 5' termini. It was demonstrated that the random 15% 13C enrichment effectively suppressed both scalar and dipolar contributions of the neighboring carbons and protons on the relaxation parameters. To analyze dynamics via all observed relaxation parameters, full spectral density mapping (1992, J. W. Peng and G. Wagner, J. Magn. Reson. 98, 308) and the "model-free" approach (1982, Lipari and Szabo, J. Am. Chem. Soc. 104, 4546) were applied complementarily. A linear correlation between three spectral density values, J(omegaC), J(omegaH - omegaC), and J(omegaH + omegaC) was observed in plots containing all measured values, but not for the other spectral density terms including J(0). These linear correlations reflect the effect of overall motion and similar internal motions for each CH vector in the decamer. The correlations yielded two correlation times, 3-4 ns and 10-200 ps. One value, 3-4 ns, corresponds to the value of 3.3 ns obtained for the overall isotropic tumbling correlation time determined from analysis of 13C T1/T2 ratios. The possibility of overall anisotropic tumbling was examined, but statistical analysis showed no advantage over the assumption of simple isotropic tumbling. Lack of correlations entailing J(0) implies that a relatively slow chemical exchange contributes to yielding of effective Jeff(0) values. Based on spectral density mapping and the T1/T2 ratio analysis, three basic assumptions were initially employed (and subsequently justified) for the model-free calculation: isotropic overall tumbling, one internal motion, and the presence of chemical exchange terms. Except for terminal residues, the order parameter S2 and the corresponding fast internal motion correlation time were determined to be about 0.8 +/- 0.1 and 20 +/- 20 ps, respectively, for the various CH vectors. Only a few differences were observed between or within sugars and bases. The internal motion is very fast (ps-ns time scale) and its amplitude restricted; e.g., assuming a simple wobble-in-a-cone model, the internal motion is restricted to an angular amplitude of +/-22. 5 degrees for each of the 1', 3', 4', 2, and 8 positions in the purine nucleotides in the entire duplex.

Small structural ensembles for a 17-nucleotide mimic of the tRNA T psi C-loop via fitting dipolar relaxation rates with the quadratic programming algorithm

Solution structures are typically average structures determined with the help of nmr-derived distance and torsion angle information. However, when a biomolecule populates significantly different conformations, the average structure might be prone to artifacts, and other refinement strategies are necessary. For example, when experimental restraints are used in molecular dynamics simulations in a time-averaged fashion (MDtar), the experimental structural information does no longer need to be satisfied at each step of the simulation; instead, the whole trajectory must agree with the restraints. However, the resulting structural ensembles are large and not unique and it is not trivial to extract the essential dynamic features for a system. Here we demonstrate that large MDtar ensembles can be simplified substantially by reducing the number of members to just a few on the basis of adjusting the individual probabilities of the members with the PDQPRO program [N. B. Ulyanov et al. Biophysical Journal (1995), Vol. 68, p. 13]. This algorithm finds the global minimum for a search function that represents the best match of a given ensemble with the experimental dipolar interproton relaxation rates. We have applied this strategy to a 17-residue RNA hairpin, whose loop exhibited considerable flexibility evident from nmr data. This 17mer is a mimic of the T psi C-loop of tRNA, where nucleotide 54 is usually a ribosylthymidine. The methylation of U54, which is completely buried in transfer ribonucleic acid, is administered by tRNA (m5 U54)-methyltransferase (RUMT). Since the 17mer is a good substrate for RUMT, we previously concluded that the flexibility of the 17mer's loop is a key to how RUMT gains access to the methylation site [L. J. Yao et al. Journal of Biomolecular NMR (1996) Vol. 9. p. 229]. Application of the PDQPRO algorithm to the previously acquired MDtar trajectories allowed us to reduce the number of conformations from several hundred to one major and five or six minor conformations with individual populations from approximately 5% to approximately 50% without any deterioration in the match with the experimental data. The major conformation exhibits a continuation of A-form helicity through part of the loop, involving C60 and U59. In this and most other conformations the methylation site in U54 is no longer buried.

A pseudosquare knot structure of DNA in solution

We report a high-resolution NMR structure of a homodimer formed by a synthetic 25 residue DNA oligonucleotide GCTCCCATGGTTTTTGTGCACGAGC. This structure presents a novel structural motif for single-stranded nucleic acids, called a pseudosquare knot (PSQ). The oligonucleotide was originally designed to mimic a slipped-loop structure (SLS), another "unusual" DNA structure postulated as an alternative conformation for short direct repeats in double-stranded DNA. The design of the sequence is compatible with both SLS and PSQ structures, both of which possess identical sets of base-paired and unpaired nucleotides but different tertiary folds. We used deuteration of the H8 positions of purines to ascertain that the PSQ is actually formed under the conditions used. The PSQ structure was solved based on homonuclear proton nuclear Overhauser effect data using complete relaxation matrix methods. The structure essentially consists of two side-by-side helices connected by single-stranded loops. Each of the helices is well-defined; however, the relative orientation of the two remains undetermined by the NMR data. The sequences compatible with the PSQ formation are frequent in single-stranded genomes; this structure may play a role as a dimerization motif.

Nuclear magnetic resonance structure of d(GCATATGATAG). d(CTATCATATGC): a consensus sequence for promoters recognized by sigma K RNA polymerase

The three-dimensional structure of d(GCATATGATAG).d(CTATCATATGC), from the promoter region of a gene regulating sporulation in Bacillus subtilis mother cells, was determined utilizing two-dimensional nuclear Overhauser effect (2D NOE) and double-quantum-filtered COSY (2QF-COSY) spectra. To minimize the effect of methods used to obtain restraints and refine structure, several variables were studied. Interproton distance bounds were calculated very conservatively by running the complete relaxation matrix program MARDIGRAS hundreds of times using 2D NOE spectra for exchangeable and for nonexchangeable protons at different mixing times, assuming different overall correlation times and different starting structures. The 435 distance restraints were used with two different structural refinement methods: restrained molecular dynamics (rMD) and restrained Monte Carlo calculations (rMC). Refinement using different procedures and starting structures resulted in essentially the same structure (<0.8 A rmsd), indicating that the structure is defined by experimental restraints and not the refinement method or variables used. R factors indicate the structures fit the experimental NOE data very well. Some helical parameters, notably large negative X displacement, are characteristic of A-DNA, but others are characteristic of B-DNA. As with TG.CA steps in other duplex DNA sequences studied in our laboratory, the two TG.CA steps have a positive roll, with T6-G7 exhibiting the largest, and consequently a bent helix axis. The converged structure represents a time-averaged structure. However, multiple conformations, especially in deoxyriboses, were evident from vicinal coupling constants obtained from quantitative simulations of 2QF-COSY cross-peaks and from persistent inconsistencies in experimental distances due to nonlinear conformational averaging.

Insights into the dynamic nature of DNA duplex structure via analysis of nuclear Overhauser effect intensities

Sequence-dependent structures of DNA duplexes in solution can be reliably determined using NMR data if care is taken to determine restraint bounds accurately. This entails use of complete relaxation matrix methods to analyze nuclear Overhauser effect (NOE) spectroscopic cross-peak intensities, yielding accurate distance restraints. NMR studies of various DNA duplexes have suggested that there may be some limited internal motions. First, it is typically not possible to reconcile all vicinal proton coupling constants in deoxyribose rings with a single conformer. In addition, with the increased accuracy of interproton distance measurements afforded by the complete relaxation matrix algorithm MARDIGRAS, we find inconsistencies in certain distances which can most readily be ascribed to limited conformational flexibility, since conformational averaging is nonlinear. As base-sugar interproton distances depend on both sugar pucker and glycosidic torsion angle chi, motion involving these structural variables should be reflected by experimental data. Possible motional models have been considered to account for all of the data for three DNA duplexes. Analysis of intraresidue base-sugar interproton NOE bounds patterns suggests a motional model with individual sugars in equilibrium between S (2'-endo) and N (3'-endo) conformations, with S being preferred. As sugar repuckering is correlated with changes in glycosidic torsion angle chi, different sugar conformers imply different values for chi, but this is insufficient to account for all data. A two-state jump between anti and syn glycosidic conformers was considered, but it was incapable of accounting for all data. However, a model with restricted diffusion (rocking) about the glycosidic bond in addition to sugar repuckering was capable of accommodating all experimental data. This motional model is in qualitative agreement with experimental 13C relaxation-derived order parameter values in a DNA duplex.

Structure of the dimer initiation complex of HIV-1 genomic RNA

Retroviral genomes must dimerize to be fully infectious. Dimerization is directed by a unique RNA hairpin structure with a palindrome in its loop: hairpins of two strands first associate transiently through their loops, and then refold to a more stable, linear duplex. The structure of the initial, kissing-loop dimer from HIV-1, solved using 2D NMR, is bent and metastable, its interface being formed not only by standard basepairing between palindromes, but also by a distinctive pattern of interstrand stacking among bases at the stem-loop junctions. This creates mechanical distortions that partially melt both stems, which may facilitate spontaneous refolding of this RNA complex into linear form.

Structure-based design of ligands for protein basic domains: application to the HIV-1 Tat protein

A methodology has been developed for designing ligands to bind a flexible basic protein domain where the structure of the domain is essentially known. It is based on an empirical binding free energy function developed for highly charged complexes and on Monte Carlo simulations in internal coordinates with both the ligand and the receptor being flexible. HIV-1 encodes a transactivating regulatory protein called Tat. Binding of the basic domain of Tat to TAR RNA is required for efficient transcription of the viral genome. The structure of a biologically active peptide containing the Tat basic RNA-binding domain is available from NMR studies. The goal of the current project is to design a ligand which will bind to that basic domain and potentially inhibit the TAR-Tat interaction. The basic domain contains six arginine and two lysine residues. Our strategy was to design a ligand for arginine first and then a superligand for the basic domain by joining arginine ligands with a linker. Several possible arginine ligands were obtained by searching the Available Chemicals Directory with DOCK 3.5 software. Phytic acid, which can potentially bind multiple arginines, was chosen as a building block for the superligand. Calormetric binding studies of several compounds to methylguanidine and Arg-/Lys-containing peptides were performed. The data were used to develop an empirical binding free energy function for prediction of affinity of the ligands for the Tat basic domain. Modeling of the conformations of the complexes with both the superligand and the basic domain being flexible has been carried out via Biased Probability Monte Carlo (BPMC) simulations in internal coordinates (ICM 2.6 suite of programs). The simulations used parameters to ensure correct folding, i.e., consistent with the experimental NMR structure of a 25-residue Tat peptide, from a random starting conformation. Superligands for the basic domain were designed by joining together two molecules of phytic acid with peptidic and peptidomimetic linkers. The linkers were refined by varying the length and side chains of the linking residues, carrying out BPMC simulations, and evaluation of the binding free energy for the best energy conformation. The dissociation constant of the best ligand designed is estimated to be in the low- to mid-nanomolar range.

NMR spectroscopy and MRI investigation of a potential bioartificial liver

NMR feasibility was established for a coaxial hydrophobic-membrane bioreactor containing isolated rat hepatocytes with features designed to mimic the human liver. A novel triple-tuned NMR probe and a perfusion system controlling temperature, gas concentrations, flow-rate, and pH were used. We determined the optimum coaxial interfiber distance (i.e. diffusion distance) for maintaining hepatocyte viability in two bioreactor prototypes. Prototype no. 1 and no. 2 had diffusion distances of 500 microns and 200 microns, respectively. Cell viability was established by 31P NMR and trypan blue exclusion. Only prototype no. 2 maintained cell viability for more than 6 h, indicating the importance of diffusion distance. 31P spectra obtained over this 6 h time period were similar to in vivo spectra of rat liver. The 31P spectra were found to be more sensitive to subacute cell viability than trypan blue exclusion. In the 1H and 31P spectra, 1H2O and inorganic phosphate signals were split in two at all flow-rates, probably due to bulk magnetic susceptibility effects originating from the three bioreactor compartments. MRI was useful for quality control and determining flow dynamics, fiber integrity, and cell inoculate distribution. MRI revealed that the inner fibers were not centered in either prototype. Although an increased flow-rate did not influence spectral resolution or chemical shifts, significant degradation of MRI quality occurred above 50 mL/min. NMR spectroscopy and imaging provide valuable, real-time information on cell biochemistry and flow dynamics which can be used in development and monitoring of bioreactors designed as artificial livers.

Fructose-1,6-bisphosphate preserves adenosine triphosphate but not intracellular pH during hypoxia in respiring neonatal rat brain slices

BACKGROUND

Fructose-1,6-bisphosphate (FBP) sometimes provides substantial cerebral protection during hypoxia or ischemia. 31P/1H nuclear magnetic resonance spectroscopy of cerebrocortical slices was used to study the effects of FBP on hypoxia-induced metabolic changes. In addition, 13C-labeled glucose was administered and 13C nuclear magnetic resonance spectroscopy was used to search for FBP-induced modulations in glycolysis and the pentose-phosphate pathway.

METHODS

In each experiment, 80 slices (350 microm) obtained from ten 7-day-old Sprague-Dawley rat litter mates were placed together in a 20-mm nuclear magnetic resonance tube, perfused, and subjected to 30 min of hypoxia (PO2 < 3 mmHg). Nine experiments were performed, with n = 3 in each of three groups: (1) no treatment with FBP; (2) 60 min of prehypoxia treatment with FBP (2 mM); and (3) 60 min of posthypoxia treatment with FBP (2 mM). 31P/1H Interleaved nuclear magnetic resonance spectra at 4.7 T provided average adenosine triphosphate, intracellular pH, and lactate. Cresyl violet stains of random slices taken at predetermined time points were studied histologically. Some experiments had [2-13C]glucose in the perfusate. Slices from these studies were frozen for perchloric acid extraction of intracellular metabolites and studied with high-resolution 13C nuclear magnetic resonance spectroscopy at 11.75 T.

RESULTS

With no pretreatment with FBP, hypoxia caused an approximately 50% loss of adenosine triphosphate, an approximately 700% increase in lactate, and a decrease in intracellular pH to approximately 6.4. Pretreatment with FBP resulted in no detectable loss of adenosine triphosphate, no increase in lactate, and minimal morphologic changes but did not alter decreases in intracellular pH. 13C Nuclear magnetic resonance spectra of extracted metabolites showed that pretreatment caused accumulation of [1-13C]fructose-6-phosphate, an early pentose-phosphate pathway metabolite. Posthypoxic treatment with FBP had no effects compared with no treatment.

CONCLUSIONS

During severe hypoxia, pretreatment with FBP completely preserves adenosine triphosphate and almost completely preserves cell morphology but does not alter hypoxia-induced decreases in intracellular pH. Pretreatment also substantially augments the flux of glucose into the pentose-phosphate pathway.

Structure of the recombinant full-length hamster prion protein PrP(29-231): the N terminus is highly flexible

The prion diseases seem to be caused by a conformational change of the prion protein (PrP) from the benign cellular form PrPC to the infectious scrapie form PrPSc; thus, detailed information about PrP structure may provide essential insights into the mechanism by which these diseases develop. In this study, the secondary structure of the recombinant Syrian hamster PrP of residues 29-231 [PrP(29-231)] is investigated by multidimensional heteronuclear NMR. Chemical shift index analysis and nuclear Overhauser effect data show that PrP(29-231) contains three helices and possibly one short beta-strand. Most striking is the random-coil nature of chemical shifts for residues 30-124 in the full-length PrP. Although the secondary structure elements are similar to those found in mouse PrP fragment PrP(121-231), the secondary structure boundaries of PrP(29-231) are different from those in mouse PrP(121-231) but similar to those found in the structure of Syrian hamster PrP(90-231). Comparison of resonance assignments of PrP(29-231) and PrP(90-231) indicates that there may be transient interactions between the additional residues and the structured core. Backbone dynamics studies done by using the heteronuclear [1H]-15N nuclear Overhauser effect indicate that almost half of PrP(29-231), residues 29-124, is highly flexible. This plastic region could feature in the conversion of PrPC to PrPSc by template-assisted formation of beta-structure.

A conformational transition at the N terminus of the prion protein features in formation of the scrapie isoform

The scrapie prion protein (PrPSc) is formed from the cellular isoform (PrPC) by a post-translational process that involves a profound conformational change. Linear epitopes for recombinant antibody Fab fragments (Fabs) on PrPC and on the protease-resistant core of PrPSc, designated PrP 27-30, were identified using ELISA and immunoprecipitation. An epitope region at the C terminus was accessible in both PrPC and PrP 27-30; in contrast, epitopes towards the N-terminal region (residues 90 to 120) were accessible in PrPC but largely cryptic in PrP 27-30. Denaturation of PrP 27-30 exposed the epitopes of the N-terminal domain. We argue from our findings that the major conformational change underlying PrPSc formation occurs within the N-terminal segment of PrP 27-30.

Evidence for protein X binding to a discontinuous epitope on the cellular prion protein during scrapie prion propagation

Studies on the transmission of human (Hu) prions to transgenic (Tg) mice suggested that another molecule provisionally designated protein X participates in the formation of nascent scrapie isoform of prion protein (PrPSc). We report the identification of the site at which protein X binds to the cellular isoform of PrP (PrPC) using scrapie-infected mouse (Mo) neuroblastoma cells transfected with chimeric Hu/MoPrP genes even though protein X has not yet been isolated. Substitution of a Hu residue at position 214 or 218 prevented PrPSc formation. The side chains of these residues protrude from the same surface of the C-terminal alpha-helix and form a discontinuous epitope with residues 167 and 171 in an adjacent loop. Substitution of a basic residue at positions 167, 171, or 218 also prevented PrPSc formation: at a mechanistic level, these mutant PrPs appear to act as "dominant negatives" by binding protein X and rendering it unavailable for prion propagation. Our findings seem to explain the protective effects of basic polymorphic residues in PrP of humans and sheep and suggest therapeutic and prophylactic approaches to prion diseases.

Solution structure of a 142-residue recombinant prion protein corresponding to the infectious fragment of the scrapie isoform

The scrapie prion protein (PrPSc) is the major, and possibly the only, component of the infectious prion; it is generated from the cellular isoform (PrPC) by a conformational change. N-terminal truncation of PrPSc by limited proteolysis produces a protein of approximately 142 residues designated PrP 27-30, which retains infectivity. A recombinant protein (rPrP) corresponding to Syrian hamster PrP 27-30 was expressed in Escherichia coli and purified. After refolding rPrP into an alpha-helical form resembling PrPC, the structure was solved by multidimensional heteronuclear NMR, revealing many structural features of rPrP that were not found in two shorter PrP fragments studied previously. Extensive side-chain interactions for residues 113-125 characterize a hydrophobic cluster, which packs against an irregular beta-sheet, whereas residues 90-112 exhibit little defined structure. Although identifiable secondary structure is largely lacking in the N terminus of rPrP, paradoxically this N terminus increases the amount of secondary structure in the remainder of rPrP. The surface of a long helix (residues 200-227) and a structured loop (residues 165-171) form a discontinuous epitope for binding of a protein that facilitates PrPSc formation. Polymorphic residues within this epitope seem to modulate susceptibility of sheep and humans to prion disease. Conformational heterogeneity of rPrP at the N terminus may be key to the transformation of PrPC into PrPSc, whereas the discontinuous epitope near the C terminus controls this transition.

Rotational correlation times of internuclear vectors in a DNA duplex with G-A mismatch determined in aqueous solution by complete relaxation matrix analysis of off-resonance ROESY (O-ROESY) spectra

Complete relaxation matrix analysis of the off-resonance ROESY (O-ROESY) spectra is presented and demonstrated for a synthetic DNA duplex with two G-A mismatches, d(GCTGTC-GAAAGC)2, in solution. The internuclear distance and the rotational correlation time of the internuclear vector could be garnered simultaneously using complete relaxation matrix analysis of O-ROESY, by which spin diffusion effects could be accommodated. Correlation times in the terminal and the mismatched regions were significantly reduced compared to those in other regions, indicating the conformational flexibility of the mismatched pair. The average structure obtained by restrained molecular dynamics simulation with inclusion of variations of the rotational correlation times also indicated a general tendency of the mismatched and contiguous bases to flip to the outside of the double strand. Off-resonance ROESY combined with the complete relaxation matrix analysis method may offer an alternative way to investigate the structures and dynamics of biological macromolecules.

Solution structure and backbone dynamics of the human alpha3-chain type VI collagen C-terminal Kunitz domain,

The solution structure and backbone dynamics of the 58-residue C-terminal Kunitz domain fragment [alpha3(VI)] of human alpha3-chain type VI collagen has been studied by two-dimensional 1H-1H and 1H-15N nuclear magnetic resonance spectroscopy at 303 K. The solution structure is represented by an ensemble of 20 structures calculated with X-PLOR using 612 distance and 47 dihedral angle restraints. The distance restraints were obtained by a complete relaxation matrix analysis using MARDIGRAS. The root mean squared (rms) deviation is 0.91 A for the backbone atoms of the residues Thr2(8)-Gly12(18), Arg15(21)-Tyr35(41), and Gly40(46)-Pro57(63). The central beta-sheet [residues Ile18(24)-Tyr35(41)] and the C-terminal alpha-helix [residues Gln48(54)-Cys55(61)] are better defined with a backbone rms deviation of 0.46 A. The solution structure of alpha3(VI) is virtually identical to the crystal structure of alpha3(VI) and to the solution structure of bovine pancreatic trypsin inhibitor (BPTI). The 15N spin-lattice and spin-spin relaxation rates and the 1H-15N heteronuclear nuclear Overhauser enhancement (NOE) were analyzed using both the "model-free" formalism [Lipari, G., & Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546-4559 and 4559-4570] and the reduced spectral density mapping procedure [Farrow, N. A., Szabo, A., Torchia, D. A., & Kay, L. E. (1995) J. Biomol.NMR 6, 153-162]. The results obtained from the "model-free" analysis include an overall correlation time tauc of 3. 00 ns and backbone order parameters S2 in the range from 0.28 to 0. 93. The necessity of including an exchange term in the analysis of the relaxation data from 14 residues indicated that these residues are involved in motions on the micro- to millisecond time scale. The majority of the 14 residues are located in the vicinity of the Cys14(20)-Cys38(44) disulfide bond, suggesting the presence of a disulfide bond isomerization similar to the one observed in BPTI [Otting, G., Liepinsh, E., & Wüthrich, K. (1993) Biochemistry 32, 3571-3582]. It is suggested that this disulfide bond isomerization is the main reason for the surprisingly small effect on trypsin inhibition observed when Thr13(19) of alpha3(VI) is substituted with Pro.

The dynamic NMR structure of the T psi C-loop: implications for the specificity of tRNA methylation

tRNA (m5U54)-methyltransferase (RUMT) catalyzes the S-adenosylmethionine-dependent methylation of uridine-54 in the T psi C-loop of all transfer RNAs in E. coli to form the 54-ribosylthymine residue. However, in all tRNA structures, residue 54 is completely buried and the question arises as to how RUMT gains access to the methylation site. A 17-mer RNA hairpin consisting of nucleotides 49-65 of the T psi-loop is a substrate for RUMT. Homonuclear NMR methods in conjunction with restrained molecular dynamics (MD) methods were used to determine the solution structure of the 17-mer T-arm fragment. The loop of the hairpin exhibits enhanced flexibility which renders the conventional NMR average structure less useful compared to the more commonly found situation where a molecule exists in predominantly one major conformation. However, when resorting to softer refinement methods such as MD with time-averaged restraints, the conflicting restraints in the loop can be satisfied much better. The dynamic structure of the T-arm is represented as an ensemble of 10 time-clusters. In all of these, U54 is completely exposed. The flexibility of the T psi-loop in solution in conjunction with extensive binding studies of RUMT with the T psi C-loop and tRNA suggest that the specificity of the RUMT/ tRNA recognition is associated with tRNA tertiary structure elements. For the methylation, RUMT would simply have to break the tertiary interactions between the D- and T-loops, leading to a melting of the T-arm structure and making U54 available for methylation.

Physical studies of conformational plasticity in a recombinant prion protein

PrP(Sc) is known to be the major, if not the only, component of the infectious prion. Limited proteolysis of PrP(Sc) produces an N-terminally truncated polypeptide of about 142 residues, designated PrP 27-30. Recently, a recombinant protein (rPrP) of 142 residues corresponding to the Syrian hamster PrP 27-30 was expressed in Escherichia coli and purified (Mehlhorn et al., 1996). rPrP has been refolded into both alpha-helical and beta-sheet structures as well as various intermediates in aqueous buffers. The beta-sheet state and two pH-dependent alpha-helical states were characterized by CD and NMR. The alpha-helical conformation occurred only after the formation of an intramolecular disulfide bond, whereas the beta-sheet form was accessible either with or without the disulfide. Of the different alpha-helical forms studied, only those refolded in the pH range 5-8 were substantially soluble at physiological pH, exhibiting similar conformations and monomeric analytical sedimentation profiles throughout the above pH range. Furthermore, refolded alpha-rPrP showed NMR chemical shift dispersion typical of proteins with native conformations, although 2D NMR indicated large segments of conformational flexibility. It displayed a cooperative thermal denaturation transition; at elevated temperatures, it converted rapidly and irreversibly to the thermodynamically more stable beta-sheet form. Unfolding of alpha-rPrP by GdnHCl revealed a two-phase transition with a relatively stable folding intermediate at 2 M GdnHCl. The deltaG values were estimated to be 1.9 +/- 0.4 kcal/mol for the first phase and 6.5 +/- 1.2 kcal/mol for the second, consistent with a folding core surrounded by significant segments of flexible conformation. By NMR, alpha-rPrP(acid) isolated at pH 2 without refolding exhibited heterogeneous line widths, consistent with an acid-denatured molten globular state. We conclude that to the extent that rPrP constitutes a relevant folding domain of PrP(C), the various conformations exhibited by rPrP suggest that the PrP sequence may be intrinsically plastic in its conformations; indeed, portions of PrP(C) may possess a relatively open conformation which makes it susceptible to conversion into PrP(Sc) under appropriate conditions.

High-resolution NMR of an antisense DNA x RNA hybrid containing alternating chirally pure Rp methylphosphonates in the DNA backbone

A high-resolution proton NMR study has been performed on a hybrid duplex formed by a methylphosphonate (MP) oligodeoxyribonucleotide (MPO) and its target oligoribonucleotide, d(T(MP)CC(MP)-TT(MP)AG(MP)CT(MP)CC(MP)TG) x r(CAGGAGCUAAGGA), where MP corresponds to positions of methylphosphonate linkages in the pure R(p) stereoconfiguration. MP-containing analogs of DNA are reported to be effective antisense agents capable of specifically inhibiting protein synthesis with the R(p) chiral MPOs exhibiting greater affinity for the target mRNA than their S(p) counterparts. Nearly complete proton resonance assignments of the hybrid duplex have been made using two-dimensional nuclear Overhauser effect (2D NOE) spectra, at three different mixing times, and double quantum-filtered COSY (2QF-COSY) spectra. The 2QF-COSY cross-peak patterns which are resolved have been analyzed qualitatively to suggest sugar conformations. Distance restraints have been obtained from the 2D NOE spectra of the duplex in D(2)O. These interproton distance restraints were determined using a complete relaxation matrix method to improve accuracy. Specifically, a new approach termed RANDMARDI has been utilized to calculate these distance restraints, accounting for spectral noise and errors in 2D NOE peak volume integration. The calculated interproton distances and sugar puckers have been analyzed to assess the solution conformation of the hybrid. The hybrid duplex appears to have an overall solution structure which is distinct from standard B- and A-forms, but the RNA strand exhibits features of the A-form. The absence of H1'-H2' cross-peaks in the 2QF-COSY spectrum indicates a C3'-endo type of conformation for ribose sugars in the RNA strand. The deoxyriboses in the antisense DNA strand exhibit a mixed behavior with almost equal scalar coupling constant values for H1'-H2' and H1'-H2" and a strong H3'-H4' 2QF-COSY peak pattern. Variations in calculated values of interproton distances and sixth-root R factor analysis of experimental intensities indicate that the hybrid duplex may have a DNA strand with significant conformational plasticity.

NMR studies of the most conserved RNA domain of the mammalian signal recognition particle (SRP)

Mammalian signal recognition particle (SRP) and its homologues exhibit a phylogenetically conserved RNA domain, whose predicted secondary structure exhibits a hairpin motif with two bulged regions. Two RNA fragments comprising one (24 nt) or two (43 nt) of the conserved bulges were studied. Each fragment binds specifically to the domain of the Escherichia coli homologue of the SRP54 protein, which is involved in signal sequence recognition. The SRP RNA fragments exhibited a pronounced structural stabilization in the presence of Mg2+. Assignments of all base, H1', H2', and most imino proton resonances in the presence of Mg2+ were obtained for the 24mer RNA via NOE spectroscopy and correlated homonuclear NMR methods. 2D NOE patterns permitted a coarse structural description, revealing a relatively compact A-type geometry for the 24mer without any indications of looped-out nucleotides, syn-oriented bases, or base triplets. The GGAA-loop is structurally very similar to that of the GCAA tetraloop [Heus HA, Pardi A, 1991, Science 253:191-194]. Mg2+ seems to stabilize the structure of the conserved bulged region, which involves G:A and C:A mismatch pairs. Deviations from ideal A-type helicity are found for a larger region than the predicted secondary structure implies. Although no explicit assignment effort has been dedicated to the 43mer yet, striking similarity in chemical shift changes upon addition of Mg2+ allowed some structural conclusions. The bulge present in both RNA fragments exhibits a similar, pronounced flexibility in the absence of Mg2+, indicating that the additional bulge in the 43mer does not stabilize the other bulge.