Structural and dynamic studies of two antigenic loops from haemagglutinin: A relaxation matrix approach

Summary We have investigated the dynamics and structural behaviour of two antigenic peptides using ¹H NMR. The two cyclic peptides mimic the antigenic site A of influenza haemagglutinin protein; they only differ in the way they were cyclized and in the size of their respective linkers. Homonuclear relaxation parameters extracted from a complete NOE matrix were interpreted in terms of local dynamics. A set of distance constraints was deduced from these parameters which allowed 3D models to be constructed using distance geometry. NOE back-calculation was used to check the validity of the final models. Strong variations of internal motion amplitude have been found in both peptides along their backbone. Motions with high amplitudes have been localized in the Gly-Pro-Gly sequence which forms a -turn in both structures.Abbreviations DSS 3-(trimethylsilyl)-1-propanesulfonic acid - D-loop aspartic acid loop - ELISA enzyme-linked immunoabsorbent assay - f.i.d free induction decay - HOHAHA homonuclear Hartmann-Hahn spectroscopy - HPLC high pressure liquid chromatography - K-loop lysine loop - NMR nuclear magnetic resonance - NOE nuclear Overhauser enhancement - NOESY nuclear Overhauser enhancement spectroscopy - r.m.s.d. root-mean-square deviation of atomic positions     

A carbon-13 nuclear magnetic resonance investigation of the metabolic fluxes associated withy glucose metabolism in human erythrocytes

We have used [2-¹³C]d-glucose and carbon-13 nuclear magnetic resonance (NMR) spectroscopy to investigate metabolic fluxes through the major pathways of glucose metabolism in intact human erythrocytes and to determine the interactions among these pathways under conditions that perturb metabolism. Using the method described, we have been able to measure fluxes through the pentose phosphate pathway, phosphofructokinase, the 2,3-diphosphoglycerate bypass, and phosphoglycerate kinase, as well as glucose uptake, concurrently and in a single experiment. We have measured these fluxes in normal human erythrocytes under the following conditions: (1) fully oxygenated; (2) treated with methylene blue; and (3) deoxygenated. This method makes it possible to monitor various metabolic effects of stresses in normal and pathological states. Not only has ¹³C-NMR spectroscopy proved to be a useful method for measuring in vivo flux through the pentose phosphate pathway, but it has also provided additional information about the cycling of metabolites through the non-oxidative portion of the pentose phosphate pathway. Our evidence from experiments with [1-¹³C]-, [2-¹³C]-, and [3-¹³C]d-glucoses indicates that there is an observable reverse flux of fructose 6-phosphate through the reactions catalyzed by transketolase and transaldolase, even in the presence of a net flux through the pentose phosphate pathway.     

Metropolis Monte Carlo calculations of DNA structure using internal coordinates and NMR distance restraints: An alternative method for generating a high-resolution solution structure

Summary A new method, a restrained Monte Carlo (rMC) calculation, is demonstrated for generating high-resolution structures of DNA oligonucleotides in solution from interproton distance restraints and bounds derived from complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect (NOE) spectral peak intensities. As in the case of restrained molecular dynamics (rMD) refinement of structures, the experimental distance restraints and bounds are incorporated as a pseudo-energy term (or penalty function) into the mathematical expression for the molecular energy. However, the use of generalized helical parameters, rather than Cartesian coordinates, to define DNA conformation increases efficiency by decreasing by an order of magnitude the number of parameters needed to describe a conformation and by simplifying the potential energy profile. The Metropolis Monte Carlo method is employed to simulate an annealing process. The rMC method was applied to experimental 2D NOE data from the octamer duplex d(GTA-TAATG)·d(CATTATAC). Using starting structures from different locations in conformational space (e.g. A-DNA and B-DNA), the rMC calculations readily converged, with a root-mean-square deviation (RMSD) of <0.3 Å between structures generated using different protocols and starting structures. Theoretical 2D NOE peak intensities were calculated for the rMC-generated structures using the complete relaxation matrix program CORMA, enabling a comparison with experimental intensities via residual indices. Simulation of the vicinal proton coupling constants was carried out for the structures generated, enabling a comparison with the experimental deoxyribose ring coupling constants, which were not utilized in the structure determination in the case of the rMC simulations. Agreement with experimental 2D NOE and scalar coupling data was good in all cases. The rMC structures are quite similar to that refined by a traditional restrained MD approach (RMSD<0.5 Å) despite the different force fields used and despite the fact that MD refinement was conducted with additional restraints imposed on the endocyclic torsion angles of deoxyriboses. The computational time required for the rMC and rMD calculations is about the same. A comparison of structural parameters is made and some limitations of both methods are discussed with regard to the average nature of the experimental restraints used in the refinement.Abbreviations MC Monte Carlo - rMC restrained Monte Carlo - MD molecular dynamics - rMD restrained molecular dynamics - DG distance geometry - EM energy minimization - 2D NOE two-dimensional nuclear Overhauser effect - DQF-COSY double-quantum-filtered correlation spectroscopy - RMSD root-mean-square deviation To whom correspondence should be addressed.     

Do NOE distances contain enough information to assess the relative populations of multi-conformer structures?

Summary The feasibility of determining the relative populations of multi-conformer structures from NOE-derived distances alone is assessed. Without cross-validation of the NOE restraints, any population ratio can be refined to a similar quality of the fit. Complete cross-validation provides a less biased measure of fit and allows the estimation of the correct population ratio when used in conjunction with very tight distance restraints. With the qualitative distance restraints most commonly used in NMR structure determination, cross-validation is unsuccessful in providing the correct answer. Other experimental sources are therefore needed to determine relative populations of multi-conformer structures.To whom correspondence should be addressed.     

Dynamic modelling of a helical peptide in solution using NMR data: Multiple conformations and multi-spin effects

Summary Nuclear Overhauser effect (NOE) measurements on molecules in solution provide information about only the ensemble-averaged properties of these molecules. An algorithm is presented that uses a list of NOEs to produce an ensemble of molecules that on average agrees with these NOEs, taking into account the effect of surrounding spins on the buildup of each NOE (spin diffusion). A simplified molecular dynamics simulation on several copies of the molecule in parallel is restrained by forces that are derived directly from differences between calculated and measured NOEs. The algorithm is tested on experimental NOE data of a helical peptide derived from bovine pancreatic trypsin inhibitor.     

Interproton distance bounds from 2D NOE intensities: Effect of experimental noise and peak integration errors

Summary The effect of experimental and integration errors on the calculations in interproton distances from NOE intensities is examined. It is shown that NOE intensity errors can have a large impact on the distances determined. When multiple spin (spin diffusion) effects are significant, the calculated distances are often underestimated, even when using a complete relaxation matrix analysis. In this case, the bias of distances to smaller values is due to the random errors in the NOE intensities. We show here that accurate upper and lower bounds of the distances can be obtained if the intensity errors are properly accounted for in the complete relaxation matrix calculations, specifically the MARDIGRAS algorithm. The basic MARDIGRAS algorithm has been previously described [Borgias, B.A. and James, T.L. (1990) J. Magn. Reson., 87, 475–487]. It has been shown to provide reasonably good interproton distance bounds, but experimental errors can compromise the quality of the resulting restraints, especially for weak cross peaks. In a new approach introduced here, termed RANDMARDI (random error MARDIGRAS), errors due to random noise and integration errors are mimicked by the addition of random numbers from within a specified range to each input intensity. Interproton distances are then calculated for the modified intensity set using MARDIGRAS. The distribution of distances that define the upper and lower distance bounds is obtained by using N randomly modified intensity sets. RANDMARDI has been used in the solution structure determination of the interstrand cross-link (XL) formed between 4-hydroxymethyl-4,5,8-trimethylpsoralen (HMT) and the DNA oligomer d(5-GCGTACGC-3)₂ [Spielmann, H.P. et al. (1995) Biochemistry, 34, 12937–12953]. RANDMARDI generates accurate distance bounds from the experimental NOESY cross-peak intensities for the fixed (known) interproton distances in XL. This provides an independent internal check for the ability of RANDMARDI to accurately fit the experimental data. The XL structure determined using RANDMARDI-generated restrains is in good agreement with other biophysical data that indicate that there is no bend introduced into the DNA by the cross-link. In contrast, isolated spin-pair approximation calculations give distance restraints that, when applied in a restrained molecular dynamics protocol, produce a bent structure.Abbreviations NOE nuclear Overhauser effect - SD standard deviation - HMT 4-hydroxymethyl-4,5,8-trimethylpsoralen - XL psoralen-DNA interstrand cross-link     

Sequential resonance assignments as a basis for determination of spatial protein structures by high resolution proton nuclear magnetic resonance

A general scheme is proposed for the determination of spatial protein structures by proton nuclear magnetic resonance. The scheme relies on experimental observation by two-dimensional nuclear magnetic resonance techniques of complete throughbond and through-space proton-proton connectivity maps. These are used to obtain sequential resonance assignments for the individual residues in the amino acid sequence and to characterize the spatial polypeptide structure by a tight network of semi-quantitative, intramolecular distance constraints.     

Structure and Disorder in an Unfolded State under Nondenaturing Conditions from Ensemble Models Consistent with a Large Number of Experimental Restraints

Obtaining detailed structural models of disordered states of proteins under nondenaturing conditions is important for a better understanding of both functional intrinsically disordered proteins and unfolded states of folded proteins. Extensive experimental characterization of the drk N-terminal SH3 domain unfolded state has shown that, although it appears to be highly disordered, it possesses significant nonrandom secondary and tertiary structure. In our previous attempts to generate structural models of the unfolded state using the program ENSEMBLE, we were limited by insufficient experimental restraints and conformational sampling. In this study, we have vastly expanded our experimental restraint set to include ¹H-¹⁵N residual dipolar couplings, small-angle X-ray scattering measurements, nitroxide paramagnetic relaxation enhancements, O₂-induced ¹³C paramagnetic shifts, hydrogen-exchange protection factors, and ¹⁵N R₂ data, in addition to the previously used nuclear Overhauser effects, amino terminal Cu²⁺-Ni²⁺ binding paramagnetic relaxation enhancements, J-couplings, chemical shifts, hydrodynamic radius, and solvent accessibility restraints. We have also implemented a new ensemble calculation methodology that uses iterative conformational sampling and seeks to calculate the simplest possible ensemble models. As a result, we can now generate ensembles that are consistent with much larger experimental data sets than was previously possible. Although highly heterogeneous and having broad molecular size distributions, the calculated drk N-terminal SH3 domain unfolded-state ensembles have very different properties than expected for random or statistical coils and possess significant nonnative α-helical structure and both native-like and nonnative tertiary structure.