Analysis of phi and chi 1 torsion angles for hen lysozyme in solution from 1H NMR spin-spin coupling constants

Three-bond 3JHN alpha coupling constants have been determined for 106 residues and 3J alpha beta coupling constants have been measured for 57 residues of the 129-residue protein hen egg white lysozyme. These NMR data have been compared with torsion angles defined in the tetragonal and the triclinic crystal forms of the protein. For most residues the measured 3JHN alpha values were consistent with the phi torsion angles found in both crystal forms; the RMS difference between the coupling constants calculated by using the tetragonal crystal structure phi angles and the experimental 3JHN alpha values is 0.88 Hz. Thus there appears to be no significant averaging of the phi torsion angle either in the interior or at the surface of the protein. For 41 of the residues where 3J alpha beta coupling constants have been determined, the values are consistent with a single staggered conformation about the chi 1 torsion angle and there is complete agreement between the NMR data in solution and the torsion angles defined in the crystalline state. In contrast, for the other 16 residues where 3J alpha beta coupling constant values have been measured, the data indicate extensive motional averaging about the chi 1 torsion angle. These residues occur largely on the surface of the protein and examination of the crystal structures shows that many of these residues adopt a different conformation in the triclinic and tetragonal crystal forms and have high crystallographic temperature factors. It appears, however, that in solution conformational flexibility of the side chains of surface residues is significantly more pronounced than in individual crystal structures.     

alpha and beta-glycopyranosyl phosphates and 1.2-phosphates. Assignments of conformations in solution by 13C and 1H NMR.

The 1H and 13C NMR parameters of the anomeric pairs of aldopyranosyl phosphates and their rigid 1,2-phosphate derivatives are reported.The derivatives of D-glucose, D-galactose, and D-mannose exist in the 4C1 conformation while the L-fuco derivatives are in the C4 conformation. As judged by 31P--1H and 31P--13C coupling constants, all of the alpha anomers of the aldopyranosyl phosphates have the phosphate moiety predominantly trans to C(2) while in the beta anomers other rotamers make significant contributions. This relationship remains the same for the biologically important nucleoside diphosphate sugars (UDPGlc, UDPGal, GDPMan, and GDPFuc). From the pH dependence of 13C chemical shifts, observed in 0.5 M solutions, the pK'a2 of the alpha anomers is 6.1 while the pK'a2 of the beta anomers is 0.6--0.8 pH unit lower. In the 1.2-phosphates, the chair conformation of the parent aldose is retained while an envelope conformation is formed by the cyclic phosphate. In the alpha anomers, the plane is formed between C(2), C(1), O(1), and P while O(2) is above the plane. In the beta anomers, O(1) is out of the plane formed by the other atoms. The beta anomers have phosphorus coupled to C(3) with coupling constants of 10.8--11.7 Hz, approximately 2 Hz greater than the maximum reported for trans coupling (Lapper, R. D., & Smith, I. C. P. (1973) J. Am. Chem. Soc. 95, 2880).     

1H NMR assignments and conformational analysis of the oligoribonucleotides CA, CAU, CAUG, ACAUG, and UCAUG: observation of pyrimidine H5-H1' long-range scalar couplings

Conformational analysis and 1H NMR spectral assignments have been carried out using COSY and RELAY methods for a series of related oligoribonucleotides including two pentamers with 5'-dangling bases. Intraresidue long-range five bond scalar coupling was observed between pyrimidine H5 and H1' protons in the COSY-45 spectra and this feature was useful for both assignment purposes and conformational analysis. The ribose ring conformations were predominantly C3'-endo with the C2'-endo population increasing at the 3'-terminus. The 5'-dangling bases were not stacked efficiently, exhibiting lower % C3'-endo values than their 3'-nearest neighbors. Backbone torsion angle population. beta t, gamma +, epsilon t, were determined using 1H-1H, 1H-31P, and 13C-31P coupling constants. From beta t and gamma + populations the U3-G4 step in CAUG was found to be less efficiently stacked than the C1-A2 and A2-U3 steps. This observation in solution is consistent with the fiber diffraction A-RNA model (S. Arnott, D.W.L. Hukins, S.D. Dover, W. Fuller and A.R. Hodgson, J. Mol. Biol. 81, 107-122, 1973) which also predicts poor stacking in a U-G dinucleotide. The epsilon t populations were greater than 65% for all C3'-O3' bonds and consistent with a right-handed A-RNA helix.     

Nuclear magnetic resonance analyses of side chain conformations of histidine and aromatic amino acid derivatives. Solvent and pH dependence

Stereoselectively beta-deuterated species were synthesized of Ac-His-NHMe, Ac-His-OEt, Ac-His-OH and H-His-NHMe, which are useful as models of histidine residues in peptides. From the spectral comparison of 1H n.m.r., the beta-proton resonances of the normal species were unambiguously assigned. In (C2H3)2SO, C2(2)H5O2H, C2H3O2H, and C5(2)H5N solution and in aqueous solution, the lower-field and higher-field components of beta-proton resonances of the four histidine derivatives are assigned to the pro-R and pro-S protons, respectively. The alternative assignments apply for Ac-His-NHMe, Ac-His-OEt and Ac-His-OH in non-polar solvents such as C2HCl3. Vicinal coupling constants 3J alpha beta S and 3J alpha beta R were obtained for calculating the fractional populations of rotamers about the C alpha-C beta bond. The rotamer populations depend little on the ionization states of the alpha-amino and carboxyl groups or the imidazole ring. The rotamer populations depend significantly on the solvent polarity, similar to those of Phe, Tyr and Trp derivatives. For the two beta-proton resonances of His, Phe, Tyr, and Trp derivatives in a variety of solvents, linear relationships are found between the differences in chemical shifts and the differences in vicinal coupling constants.     

The solution structure of a cyclic endothelin antagonist, BQ-123, based on 1H-1H and 13C-1H three bond coupling constants.

A cyclic pentapeptide endothelin antagonist, cyclo(dTrp-dAsp-Pro-dVal-Leu), recently reported (K. Ishikawa et al., 13th Am. Pept. Symp., Cambridge MA, 1991) has been studied by NMR spectroscopy and molecular modeling. A stable structure has been determined without the use of nuclear Overhauser effects and is based primarily on homonuclear and heteronuclear three bond coupling constants. The 13C-edited TOCSY experiment is demonstrated at natural abundance and approximately 30 mM peptide concentrations. Three bond 13C-1H coupling constants obtained by this method are shown to reduce the ambiguity in phi angle determination which exists when only interproton coupling constants are used. Three out of four phi angles were determined uniquely by this method and the fourth was reduced to two possible values. The proline phi angle was determined to be -78 degrees based on the 3JH alpha, H beta and 3JH alpha, H beta coupling constants. Comparison of amide proton temperature dependence, chemical shifts and vicinal proton coupling constants in a 20% acetonitrile/80% water solvent mixture and in (CD3)2SO indicates that the structure is similar in both solvents.     

The binding of human complement proteins C5, factor B, beta 1H and properdin to complement fragment C3b on zymosan.

The covalent binding of complement fragment C3b to zymosan by the action of the alternative-pathway C3 convertase and the reversible binding of several complement proteins (component C5, factor B, beta 1H and properdin) to C3b on zymosan have been investigated. When C3b is deposited on zymosan after activation by a surface-bound C3 convertase, the C3b molecules are deposited in foci around the C3 convertase site, with an average of 30 C3b molecules per site. The association constants of C5, factor B, beta 1H, and properdin for C3b bound to zymosan have been determined. The association constants ranged from 6.5 x 10(-5) M-1 for factor B to 2.9 x 10(7) M-1 for properdin. An approximate stoichiometry of 1 : 1 for C5, factor B, and properdin binding to C3b has been observed. Curvilinear Scatchard plots were observed for beta 1H binding to C3b, with the maximal extrapolated ratio of beta 1H to C3b of 0.32. Physiological amounts of properdin increase by 7-fold the affinity constant for factor B binding to C3b with no alteration in the stoichiometry. Similarly, physiological amounts of factor B increase the affinity constant of properdin to C3b about 4-fold with only a small measured difference in stoichiometry. Competition binding studies and protein modification suggest that C5, factor B, beta 1H, and properdin each bind to a distinct region on C3b.     

Determination of the local structure of the protein insectotoxin I5A from the scorpion Buthus eupeus from 1H-NMR spectroscopy data

The local structure (torsion angles phi, psi and chi 1 of amino acid residues) of insectotoxin I5A (35 residues) of scorpion Buthus eupeus has been determined from cross-peak integral intensities in two-dimensional nuclear Overhauser enhancement (NOESY) spectra and spin coupling constants of vicinal H--NC alpha--H and H--C alpha C beta--H protons. The local structure determination was carried out by fitting complete relaxation matrix of peptide unit protons (protons of a given residue and NH proton of the next residue in the amino acid sequence) with experimental NOESY cross-peak intensities. The obtained intervals of backbone torsional angles phi and psi consistent with NMR data were determined for all but Gly residues. The predominant C alpha--C beta rotamer of the side chain has been unambiguously determined for 42% of the insectotoxin amino acid residues whereas for another 46% residues experimental data are fitted equally well with two rotamers. Stereospecific assignments were obtained for 38% of beta-methylene groups. The determined torsional angles phi, psi and chi 1 correspond to the sterically allowed conformations of the amino acid residues and agree with the insectotoxin secondary structure established earlier by 1H NMR spectroscopy.     

Two-dimensional 1H and 31P NMR spectra and restrained molecular dynamics structure of a mismatched GA decamer oligodeoxyribonucleotide duplex

Assignment of the 1H and 31P NMR spectra of a tandem G.A mismatched base pair decamer oligodeoxyribonucleotide duplex, d(CCAAGATTGG)2, has been made by two-dimensional 1H-1H and heteronuclear 31P-1H correlated spectroscopy. Unusual downfield 31P resonances have been assigned by a pure absorption phase constant-time heteronuclear 31P-1H correlated spectrum to be associated with the phosphates on the 5'- and 3'-sides of the mismatched guanosine residue. JH3'-P coupling constants for each of the phosphates of the decamer were obtained from the 1H-31P J-resolved selective proton-flip 2D spectrum. The two most downfield-shifted 31P resonances each appear to consist of two overlapping signals that can be resolved into two distinct doublets with different coupling constants in the J-resolved spectrum. This as well as the temperature dependence of the 31P spectra demonstrates that two distinct conformations exist at lower temperatures. By use of a modified Karplus relationship, the C4'-C3'-O3'-P torsional angles (epsilon) were obtained. A linear correlation between 31P chemical shifts and the measured coupling constants is quite good (only when the larger set of coupling constants of the two most downfield 31P signals is included). The 31P chemical shifts as well as the measured coupling constants tend to follow the positional variation seen in other duplexes of interior phosphates resonating more upfield than terminal residues and of interior phosphates exhibiting smaller coupling constants; however, this pattern is disrupted at the site of the mismatch. Modeling and initial NOESY distance restrained molecular mechanics energy minimization and restrained molecular dynamics support previous observations that the mismatched guanine and adenine bases are both in anti conformations. Most significantly, the epsilon backbone torsional angle variaions calculated from the NOESY distance restrained structures are in agreement with both the crystal structure values and the measured JH3'-P coupling constants.     

1H NMR assignments and secondary structure of human beta 2-microglobulin in solution.

Sequence-specific resonance assignments of human beta 2-microglobulin (M(r) 12,000) and its secondary structure are determined by 2D NMR techniques. The protein is found to contain two antiparallel beta-sheets each of four beta-strands with the beta-sheets being connected by a single disulfide linkage. No evidence for any regular helical structure is found. Amide proton-solvent-exchange rate constants and 3JHN alpha coupling constants are evaluated.     

1H, 13C,and 15N assignments and secondary structure of the FK506 binding protein when bound to ascomycin

The ¹H, ¹³C, and ¹⁵N resonances of FKBP when bound to the immunosuppressant, ascomycin, were assigned using a computer-aided analysis of heteronuclear double and triple resonance three-dimensional nmr spectra of [U-¹⁵N] FKBP/ascomycin and [U-¹⁵N, ¹³C] FKBP/ascomycin. In addition, from a preliminary analysis of two heteronuclear four-dimensional data sets, ³J coupling constants, amide exchange data, and the differences between the C^α and C^β chemical shifts of FKBP to random coil values, the secondary structure of FKBP when bound to ascomycin was determined. The secondary structure of FKBP when bound to ascomycin in solution closely resembled the x-ray structure of the FKBP/FK506 complex but differed in some aspects from the structure of uncomplexed FKBP in solution. © 1993 John Wiley & Sons, Inc.     

1H NMR conformational study of a variety of alpha-anomers of C5-substituted 2'-deoxyuridines: comparison to their antiherpetic beta counterparts

Although alpha-nucleosides are not found in nucleic acid, they do occur as constituents of smaller molecules in living cells, e.g., in vitamin B(12). There are now several examples of alpha-nucleosides exerting a biological activity in some instances equal to, or even exceeding, that of the corresponding beta-anomer. Examples include growth inhibitory properties against mouse leukemia cells and antitumor activity. From stereochemical point of view, alpha-anomers serve as references for studying of interaction of the base with the sugar moiety in beta-anomers and may help in better understanding of structure-activity relationships. One important problem preventing conformational analysis of alpha nucleosides is uncertainty in the determination of vicinal coupling constants from simulation of overlapping sugar proton resonances of strongly coupled spin systems. A successful resolution of near-isochronous H3' and H4' resonances made possible a full conformational analysis for a series of alpha-anomers C5-substituted 2'-deoxyuridines, including methyl, ethyl, isopropyl, fluor, vinyl, and bromovinyl, in comparison to their beta counterparts. Conformation of the sugar ring is determined from proton-proton coupling constants and described in terms of pseudorotation between two main puckering domains C2'endo (S) and C3'endo (N). A thorough analysis of chemical shifts as well as conformation of the sugar ring and C4'-C5' rotamers made possible determination of conformational preferences in equilibrium about the glycosidic bond between two regions, anti and syn. This work provides insights into the role of anomeric configuration of the base in conformational behavior of the sugar moiety, a link in the backbone of nucleic acids.     

Solution conformation of sialosylcerebroside (GM4) and its NeuAc(alpha 2----3)Gal beta sugar component

The solution conformations of GM4 ganglioside [NeuAc(alpha 2----3)Gal (beta 1----1)Cer] in (2H3C)2SO and its component disaccharide in 2H2O were investigated with the aid of 1H-nuclear magnetic resonance spectroscopy (nuclear Overhauser effect, analysis of coupling constants) and by energy-minimum calculations. The existence of three low-energy conformers obtained by theoretical calculations was supported by experimental findings in the case of GM4, whereas the disaccharide appears to exist as a mixture of two conformers.     

Correlation of tryptophan fluorescence intensity decay parameters with 1H NMR-determined rotamer conformations: [tryptophan2]oxytocin.

While the fluorescence decay kinetics of tyrosine model compounds [Laws, W. R., Ross, J. B. A., Wyssbrod, H. R., Beechem, J. M., Brand, L., & Sutherland, J. C. (1986) Biochemistry 25, 599-607] and the tyrosine residue in oxytocin [Ross, J. B. A., Laws, W. R., Buku, A., Sutherland, J. C., & Wyssbrod, H. R. (1986) Biochemistry 25, 607-612] can be explained in terms of heterogeneity derived from the three ground-state chi 1 rotamers, a similar correlation has yet to be directly observed for a tryptophan residue. In addition, the asymmetric indole ring might also lead to heterogeneity from chi 2 rotations. In this paper, the time-resolved and steady-state fluorescence properties of [tryptophan2]oxytocin at pH 3 are presented and compared with 1H NMR results. According to the unrestricted analyses of individual fluorescence decay curves taken as a function of emission wavelength and a global analysis of these decay curves for common emission wavelength-independent decay constants, only three exponential terms are required. In addition, the preexponential weighting factors (amplitudes) have the same relative relationship (weights) as the 1H NMR-determined chi 1 rotamer populations of the indole side chain. 15N was used in heteronuclear coupling experiments to confirm the rotamer assignments. Inclusion of a linked function restricting the decay amplitudes to the chi 1 rotamer populations in the individual decay curve analyses and in the global analysis confirms this correlation. According to qualitative nuclear Overhauser data, there are two chi 2 populations. Depending upon the degree of correlation between chi 2 and chi 1, there may be from three to six side-chain conformations for the tryptophan residue. The combined fluorescence and NMR results are consistent with a rotamer model in which either (i) the chi 2 rotations are fast compared to the fluorescence intensity decay of the tryptophan residue, (ii) environmental factors affecting fluorescence intensity decay properties are dominated by chi 1 interactions, or (iii) the chi 2 and chi 1 rotations are highly correlated.     

Long-range coupling between separate docking sites in interleukin-1beta

The human cytokine interleukin-1beta (IL-1beta) interacts with the interleukin type I receptor using two large docking surfaces designated A and B. Crystallographic studies reveal that a single histidine residue (His30) in IL-1beta makes critical electrostatic interactions at the receptor/ligand interface. To study the function of this residue at site A, four mutant forms of IL-1beta (H30A, H30D, H30F and H30R) were investigated. The mutation that introduces charge repulsion at His30 destabilizes the protein, but paradoxically causes the least effect on receptor binding (H30D). Mutations that enhance hydrophobic or electrostatic interactions have little effect on protein stability yet markedly impair receptor binding (H30F, H30R). All mutations can transmit effects from site A to site B, as evidenced by changes in the binding of a single-chain antibody highly specific for site B. Dihedral scalar coupling constants for the wild-type IL-1beta and the four His mutant proteins showed changes in backbone angles in residues located around site B, some approximately 30 angstroms away from His30 in site A. A comparison of native solvent exchange in wild-type and mutated IL-1beta shows transmission of local destabilization along the hydrogen bond network of the beta-sheet. Taken together, the data indicate that a single residue in site A of IL-1beta can impact stability and function through perturbations in both local and long-range contacts.     

Determination of the secondary structure and molecular topology of interleukin-1 beta by use of two- and three-dimensional heteronuclear 15N-1H NMR spectroscopy

A study of the regular secondary structure elements of recombinant human interleukin-1 beta has been carried out using NMR spectroscopy. Using a randomly 15N labeled sample, a number of heteronuclear three- and two-dimensional NMR experiments have been performed, which have enabled a complete analysis of short-, medium-, and long-range NOEs between protons of the polypeptide backbone, based on the sequence-specific resonance assignments that have been reported previously [Driscoll, P. C., Clore, G. M., Marion, D., Wingfield, P. T., & Gronenborn, A. M. (1990) Biochemistry 29, 3542-3556]. In addition, accurate measurements of a large number of 3JHN alpha coupling constants have been carried out by two-dimensional heteronuclear multiple-quantum-coherence-J spectroscopy. Amide NH solvent exchange rates have been measured by following the time dependence of the 15N-1H correlation spectrum of interleukin-1 beta on dissolving the protein in D2O solution. Analysis of these data indicate that the structure of interleukin-1 beta consists of 12 extended beta-strands aligned in a single extended network of antiparallel beta-sheet structure that in part folds into a skewed six-stranded beta-barrel. In the overall structure the beta-strands are connected by tight turns, short loops, and long loops in a manner that displays approximate pseudo-three-fold symmetry. The secondary structure analysis is discussed in the light of the unrefined X-ray structure of interleukin-1 beta at 3-A resolution [Priestle, J. P., Sch?r, H.-P., & Grütter, M. G. (1988) EMBO J. 7, 339-343], as well as biological activity data. Discernible differences between the two studies are highlighted. Finally, we have discovered conformational heterogeneity in the structure of interleukin-1 beta, which is characterized by an exchange rate that is slow on the NMR chemical shift time scale.     

New Karplus equations for 2JHH, 3JHH, 2JCH, 3JCH, 3JCOCH, 3JCSCH, and 3JCCCH in some aldohexopyranoside derivatives as determined using NMR spectroscopy and density functional theory calculations

The (1)H-(13)C coupling constants of methyl alpha- and beta-pyranosides of D-glucose and D-galactose have been measured by one-dimensional and two-dimensional (1)H-(13)C heteronuclear zero and double quantum, phase sensitive J-HMBC spectra to determine a complete set of coupling constants ((1)J(CH), (2)J(CH), (3)J(CH), (2)J(HH), and (3)J(HH)) within the exocyclic hydroxymethyl group (CH(2)OH) for each compound. In parallel with these experimental studies, structure, energy, and potential energy surfaces of the hydroxymethyl group for these compounds were determined employing quantum mechanical calculations at the B3LYP level using the 6-311++G( * *) basis set. Values of the vicinal coupling constants involving (1)H and (13)C in the C5-C6 (omega) and C6-O6 (theta) torsion angles in the aldohexopyranoside model compounds were calculated with water as the solvent using the PCM method. To test the relationship between (3)J(CXCH) (X=C, O, S) and torsion angle C1-X (phi) around the anomeric center, the conformations of 24 derivatives of glucose and galactose, which represent sequences of atoms at the anomeric center of C-glycosides (C-C bond), O-glycosides (C-O bond), thioglycosides (C-S bond), glycosylamines (C-N bond), and glycosyl halides (C-halogen (F/Cl) bond) have been calculated. Nonlinear regression analysis of the coupling constants (1)J(C1,H1), (2)J(C5,H6R), (2)J(C5,H6S), (2)J(C6,H5), (3)J(C4,H6R), (3)J(C4,H6S), (2)J(H6R,H5), and (3)J(H5,H6R) as well as (3)J(CXCH) (X=C, O, S) on the dihedral angles omega, theta, and phi have yielded new Karplus equations. Good agreement between calculated and experimentally measured coupling constants revealed that the DFT method was able to accurately predict J-couplings in aqueous solutions.     

Observation of H-bond mediated 3hJH2H3 coupling constants across Watson-Crick AU base pairs in RNA

^3hJ_H2H3trans-hydrogen bond scalar coupling constants have been observed for the first time in Watson-Crick AU base pairs in uniformly ¹⁵N-labeled RNA oligonucleotides using a new ^2hJ_NN-HNN-E. COSY experiment. The experiment utilizes adenosine H2 (AH2) for original polarization and detection, while employing ^2hJ_NNcouplings for coherence transfer across the hydrogen bonds (H-bonds). The H3 protons of uracil bases are unperturbed throughout the experiment so that these protons appear as passive spins in E. COSY patterns. ^3hJ_H2H3coupling constants can therefore be accurately measured in the acquisition dimension from the displacement of the E. COSY multiplet components, which are separated by the relatively large ¹J_H3N3coupling constants in the indirect dimension of the two-dimensional experiment. The ^3hJ_H2H3scalar coupling constants determined for AU base pairs in the two RNA hairpins examined here have been found to be positive and range in magnitude up to 1.8 Hz. Using a molecular fragment representation of an AU base pair, density functional theory/finite field perturbation theory (DFT/FPT) methods have been applied to attempt to predict the relative contributions of H-bond length and angular geometry to the magnitude of ^3hJ_H2H3coupling constants. Although the DFT/FPT calculations did not reproduce the full range of magnitude observed experimentally for the ^3hJ_H2H3coupling constants, the calculations do predict the correct sign and general trends in variation in size of these coupling constants. The calculations suggest that the magnitude of the coupling constants depends largely on H-bond length, but can also vary with differences in base pair geometry. The dependency of the ^3hJ_H2H3coupling constant on H-bond strength and geometry makes it a new probe for defining base pairs in NMR studies of nucleic acids.     

Nuclear magnetic resonance conformational studies on the chemotactic tripeptide formyl-L-methionyl-L-leucyl-L-phenylalanine. A small beta sheet.

Previous work by several groups has shown that the combination of spin--spin coupling constants and spectral density components (derived from spin--lattice relaxation and/or nuclear Overhauser measurements) may aid in the task of conformational determination of peptides in solution. Using the peptide formyl-L-methionyl-L-leucyl-L-phenylalanine, which is a potent specific chemotactic agent for leucocytes, we show the following: (a) that 3JNHCH coupling constants are consistent with a high degree of rigidity in the peptide backbone in solution, (b) that 3H isotopic substitution in combination with relaxation data taken at different Larmor frequencies enables spectral density, and thence conformational, information to be obtained, (c) that side-chain conformations for this molecule mirror, in some aspects, those found in the solid state for other peptides containing the same residues, and (d) that temperature dependence of amide chemical shifts does not have direct implication concerning the existence of intramolecular hydrogen bonds in peptides. We are able to propose a family of conformations which appear to interchange rapidly on the NMR time scale and are characterized by a distribution of side-chain rotamers. The basic backbone conformation is, or closely approximates, a small beta antiparallel pleated sheet and as such suggests a possible mode of receptor--chemotactic peptide interaction.