An NMR Study of the Spatial Structure and Intramolecular Dynamics of Modified Analogues of Steroid Hormones

Special features of the use of homo- and heteronuclear correlation methods of NMR in one and two dimensions for studying the spatial structure and intramolecular dynamics of modified analogues of steroid hormones (MASH) are considered. The application of these methods to the assignment of resonances in the high-field ¹H NMR region and to the determination of the most stereospecifically important parameters, such as the vicinal constants of spin–spin coupling (³ J _H–H) and nuclear Overhauser effects (NOE), are discussed using the example of NMR studies of some estrogens and androgens at 300 MHz and on the basis of literature data. The most efficient combination of the methods and the necessary modification of each of them may be chosen considering the spectral and relaxation parameters of MASH in liquid medium, including the anisotropy of the overall diffusive motion. The characteristics of MASH are the wide use of correlations through long-range couplings (COSY-45 and DQF-COSY), the application of the ^4,5 J _H–H constants for the determination of spatial structure, and the advantage of heteronuclear HSQC methods with and without ¹³C decoupling over the corresponding HMQC methods in both resolution and sensitivity. In the conformationally rigid MASH molecules, the anisotropy of the MASH diffusive motion in liquid adversely affects the determination of interproton distances by the calibrating processing method for the NOE difference and NOESY spectra: it results in both overestimated and underestimated distance values depending on the polar angle ratios of the reference and the determined distances. Under certain conditions, conformationally mobile MASH demonstrate the additional contribution of the scalar relaxation mechanism between the indirectly (scalarly) bound protons. This mechanism is responsible for the underestimated values of NOE and the corresponding errors in the distance determination.     

31P NMR relaxation measurements of the phosphate backbone of a double stranded hexadeoxynucleotide in solution: determination of the chemical shift anisotropy

The five phosphates of the deoxynucleotide d(CpGpTpApCpG)₂ have been assigned by two-dimensional heteronuclear NMR spectroscopy. The chemical shift anisotropy and correlation time of each phosphate group has been determined from measurements of the spin-lattice, spin-spin relaxation rate constants and the ³¹P-{¹H} nuclear Overhauser enhancement (NOE) at three magnetic field strengths (4.7 T, 9.4 T, and 11.75 T) and two temperatures (288 K and 298 K). As expected, the relaxation data require two mechanisms to account for the observed rate constants, i.e. dipole-dipole and chemical shift anisotropy. At 9.4 T and 11.75 T, the latter mechanism dominates the relaxation, leading to insignificant NOE intensities. The correlation time, chemical shift anisotropy and effective P-H distance were obtained from least-squares fitting to the data. Comparison of the fitted value for the correlation time with that obtained from ¹H measurements shows that the molecule behaves essentially as rigid rotor on the nanosecond timescale. Large amplitude motions observed in long segments of DNA are due to bending motions that do not contribute significantly to relaxation in short oligonucleotides.Abbreviations CSA chemical shift anisotropy - NOE nuclear Overhauser enhancement Offprint requests to: A. N. Lane     

Distance evaluation via heteronuclear SQC-NOESY experiments

Summary A method for the quantitative determination of interproton distances from¹H NOE relayed heteronuclear correlation is presented. Model compounds are investigated. Accurate distances are obtained if all the factors affecting such distance measurements, such as the local mobility and the presence of strong coupling in heteronuclear systems, are properly taken into account.     

Computation of relaxation matrix elements from incomplete NOESY data sets

Summary The structural determination of biological molecules in solution by NMR relies on the determination of a set of interatomic distances obtained by measurement of intramolecular nuclear Overhauser effects (NOE). It is shown in this paper that it is possible to obtain the accurate relaxation rate (and hence the interatomic distance) from the direct measurement of a single NOE signal. The precise analysis of a NOESY peak evolution with respect to the mixing time allows the evaluation of the relaxation parameters for the pair of spins under consideration. This is done without any assumption on the relaxation of unmeasured spins, or on the movement of the molecule. The theoretical basis of this method is presented. In order to evaluate the proposed method, a simulated case on the protein BPTI is studied, which shows that the method performs very well even in the case of noisy data sets.     

Relation between the spectral NMR parameters and conformation characteristics of the amino acid residue backbone. Analysis of elements of the secondary structure of proteins

Based on the analysis of the proton-proton distance dependences from the conformational characteristics of the L-amino acid residues, the correlation diagram of the NOE cross peak intensity waited values with the regions of the sterically allowed space (phi, psi) was proposed. The method for determining the dihedral angles phi, psi values using the information about NOE cross peak intensities was elaborated. By the model spectral NMR parameters of the bovine pancreatic trypsin inhibitor, it is shown that the accuracy of the angles phi, psi determination exceed the corresponding accuracy provided by other methods of the structural interpretation of the two-dimensional NMR spectroscopy data. The analysis of the waited spectral NMR parameters for the different types of protein regular secondary structures and beta-turns was performed.     

Heteronuclear three-dimensional NMR spectroscopy of the inflammatory protein C5a

The utility of three-dimensional heteronuclear NMR spectroscopy for the assignment of 1H and 15N resonances of the inflammatory protein C5a (MW 8500), uniformly labeled with 15N, is demonstrated at a protein concentration of 0.7 mM. It is shown that dramatic simplification of the 2D nuclear Overhauser effect spectrum (NOESY) is obtained by editing with respect to the frequency of the 15N heteronucleus in a third dimension. The improved resolution in the 3D experiment largely facilitates the assignment of protein NMR spectra and allows for the determination of distance constraints from otherwise overlapping NOE cross peaks for purposes of 3D structure determination. The results show that 15N heteronuclear 3D NMR can facilitate the structure determination of small proteins and promises to be a useful tool for the study of larger systems that cannot be studied by conventional 2D NMR techniques.     

High resolution nuclear magnetic resonance studies of nigeran

Polymer motion in solution can be studied by ¹³CNMR relaxation methods, which provide information about the correlation time for C-H vectors. ¹³C-Relaxation and Nuclear Overhauser Enhancement (NOE) data may frequently be combined to determine the dipole-dipole relaxation contribution. An alternative method is proposed based on a comparison of the proton spin-lattice relaxation rates of the centre proton resonances of an unlabelled molecule with the relaxation rates of the ¹³C satellites (from ¹³C labelled molecules).Selectively labelled nigeran which is an alternating $\text{1 → 3 and 1 → 4 α-}\text{d}\text{-glucan}$ has been investigated. The discussion in terms of the occurrence of different motions for each of the two units of the polymer requires an unambiguous assignment of the two anomeric carbons. For this reason a detailed assignment of the ¹H and ¹³C Nuclear Magnetic Resonance (NMR) spectra of nigeran in dimethylsulphoxide-d₆ is described, based on T₁ and NOE measurements in addition to selective homonuclear and heteronuclear spin decoupling experiments. These values are correlated with a conformation estimated by HSEA hard-spheres calculation. The measurements of the relaxation parameters for labelled and unlabelled compounds which provide an alternative determination of the ¹³C-¹H dipole-dipole relaxation contribution in a macromolecule agree well with ¹³C-{¹H} NOE experiments.     

Relaxation data in NMR structure determination: model calculations for the lysozyme-Gd3+ complex

The effect of including paramagnetic relaxation data as additional restraints in the determination of protein tertiary structures from NMR data has been explored by a systematic series of model calculations. The system used for testing the method was the 2.0 A resolution tetragonal crystal structure of hen egg white lysozyme (129 amino acid residues) and structures were generated using a version of the hybrid "distance geometry-dynamic simulated annealing" procedure. A limited set of 769 NOEs was used as restraints in all the calculations; the strengths of these were categorized into three classes on the basis of distances observed in the crystal structure. The values of 50 phi angles were also restrained on the basis of amide-alpha coupling constants calculated from the X-ray structure. Five sets of 12 structures were determined using differing sets of paramagnetic relaxation data as restraints additional to those involving the NOE and coupling constant data. The paramagnetic relaxation data were modeled on the basis of the distances of defined protons from the crystallographic binding site of Gd3+ in lysozyme. Analysis of the results showed that the relaxation data significantly improved the correspondence between the set of generated structures and the crystal structure, and that the more well defined the relaxation data, the more significant the improvement in the quality of the structures. The results suggest that the inclusion of paramagnetic relaxation restraints could be of significant value for the experimental determination of protein structures from NMR data.     

High-resolution NMR structure of an RNA model system: the 14-mer cUUCGg tetraloop hairpin RNA

We present a high-resolution nuclear magnetic resonance (NMR) solution structure of a 14-mer RNA hairpin capped by cUUCGg tetraloop. This short and very stable RNA presents an important model system for the study of RNA structure and dynamics using NMR spectroscopy, molecular dynamics (MD) simulations and RNA force-field development. The extraordinary high precision of the structure (root mean square deviation of 0.3 Å) could be achieved by measuring and incorporating all currently accessible NMR parameters, including distances derived from nuclear Overhauser effect (NOE) intensities, torsion-angle dependent homonuclear and heteronuclear scalar coupling constants, projection-angle-dependent cross-correlated relaxation rates and residual dipolar couplings. The structure calculations were performed with the program CNS using the ARIA setup and protocols. The structure quality was further improved by a final refinement in explicit water using OPLS force field parameters for non-bonded interactions and charges. In addition, the 2′-hydroxyl groups have been assigned and their conformation has been analyzed based on NOE contacts. The structure currently defines a benchmark for the precision and accuracy amenable to RNA structure determination by NMR spectroscopy. Here, we discuss the impact of various NMR restraints on structure quality and discuss in detail the dynamics of this system as previously determined.     

Gradient enhanced selective experiments in the 1H NMR chemical shift assignment of the skeleton and side-chain resonances of stigmasterol, a phytosterol derivative

The applicability of homonuclear gradient enhanced NMR experiments is demonstrated in the structure determination of steroid derivatives using stigmasterol as a model compound. High resolution 1H NMR spectra of steroids very often display well resolved multiplets usually in the low-field region, and these signals can be used as starting points in several selective NMR experiments to study scalar (J coupling), and dipolar (NOE) interactions. Selective excitation was carried out using a double pulsed-field gradient spin-echo sequence (DPFGSE) in which 180 degrees Gaussian pulses are sandwiched between sine shaped z-gradients. Scalar interactions were studied by homonuclear DPFGSE-COSY, DPFGSE-relay-COSY and DPFGSE-TOCSY experiments, while DPFGSE-NOESY was used to monitor spatial environment of the selectively excited proton. These methods provided unambiguous assignments for signals of the main skeleton and the side-chain of the steroid molecule. In addition, they allowed determination of the conformationally important homonuclear proton-proton coupling constants (J).     

Assessing potential bias in the determination of rotational correlation times of proteins by NMR relaxation

The various factors that influence the reliable and efficient determination of the correlation time describing molecular reorientation of proteins by NMR relaxation methods are examined. Nuclear Overhauser effects, spin-lattice, and spin-spin relaxation parameters of 15N NMR relaxation in ubiquitin have been determined at 17.6, 14.1, 11.7 and 9.4 Tesla. This unusually broad set of relaxation parameters has allowed the examination of the influence of chemical shift anisotropy, the functional form of the model-free spectral density, and the reliability of determined spin- spin relaxation parameters on the characterization of global tumbling of the protein. Treating the 15N chemical shift anisotropy (CSA) as an adjustable parameter, a consensus value of –170 ± 15ppm for the breadth of the chemical shift tensor and a global isotropic correlation time of 4.1ns are found when using the model-free spectral density to fit T1 and NOE data from all fields. The inclusion of T2 relaxation parameters in the determination of the global correlation time results in its increase to 4.6ns. This apparent inconsistency may explain a large portion of the discrepancy often found between NMR- and fluorescence-derived m values for proteins. The near identity of observed T2 and T1 values suggests that contributions from slow motions are not the origin of the apparent inconsistency with obtained T1 and NOE data. Various considerations suggest that the origin of this apparent discrepancy may reside in a contribution to the spectral density at zero frequency that is not represented by the simple model-free formalism in addition to the usual experimental difficulties associated with the measurement of these relaxation parameters. Finally, an axially symmetric diffusion tensor for ubiquitin is obtained using exclusively T1 and NOE data. A recommendation is reached on the types and combinations of relaxation data that can be used to reliably determine m values. It is also noted that the reliable determination of m values from 15N T1 and NOE relaxation parameters will become increasingly difficult as m increases.     

Protein dynamics and distance determination by NOE measurements

Present analysis procedures for NMR structure determination of macromolecules presuppose fixed internuclear distances. Improvement of the precision of the requisite NOE information has stimulated the use of more quantitative distance constraints thus necessitating examination as to whether the assumption of a rigid model systematically biases the distance estimates. Analysis using the simple (r-6) dependence of NOE buildup rates seriously underestimates the correct distance for spatially proximal proton pairs having fluctuations comparable to those observed in X-ray temperature factor analysis. However, by calculating the proper generalized order parameter it is shown that for nuclei undergoing rapid isotropic uncorrelated fluctuations the effective distance is identical to the distance between the mean positions of the nuclei. Similar analysis of molecular dynamics simulation data from bovine pancreatic trypsin inhibitor indicates that the distance obtained from the generalized order parameter predicts the distance between the mean positions to within a few percent regardless of the degree of correlation of the pairwise motion for virtually all main chain and dynamically constrained side chain protons.     

Unique backbone-water interaction detected in sphingomyelin bilayers with 1H/31P and 1H/13C HETCOR MAS NMR spectroscopy

Two-dimensional ¹H/³¹P dipolar heteronuclear correlation (HETCOR) magic-angle spinning nuclear magnetic resonance (NMR) is used to investigate the correlation of the lipid headgroup with various intra- and intermolecular proton environments. Cross-polarization NMR techniques involving ³¹P have not been previously pursued to a great extent in lipid bilayers due to the long ¹H-³¹P distances and high degree of headgroup mobility that averages the dipolar coupling in the liquid crystalline phase. The results presented herein show that this approach is very promising and yields information not readily available with other experimental methods. Of particular interest is the detection of a unique lipid backbone-water intermolecular interaction in egg sphingomyelin (SM) that is not observed in lipids with glycerol backbones like phosphatidylcholines. This backbone-water interaction in SM is probed when a mixing period allowing magnetization exchange between different ¹H environments via the nuclear Overhauser effect (NOE) is included in the NMR pulse sequence. The molecular information provided by these ¹H/³¹P dipolar HETCOR experiments with NOE mixing differ from those previously obtained by conventional NOE spectroscopy and heteronuclear NOE spectroscopy NMR experiments. In addition, two-dimensional ¹H/¹³C INEPT HETCOR experiments with NOE mixing support the ¹H/³¹P dipolar HETCOR results and confirm the presence of a H₂O environment that has nonvanishing dipolar interactions with the SM backbone.     

Intermolecular relaxation has little effect on intra-peptide exchange-transferred NOE intensities

Exchange-transferred nuclear Overhauser enhancement (etNOE) provides a useful method for determining the 3-dimensional structure of a ligand bound to a high-molecular-weight complex. Some concern about the accuracy of such structures has arisen because indirect relaxation can occur between the ligand and macromolecule. Such indirect relaxation, or spin diffusion, would lead to errors in interproton distances used as restraints in structure determination. We address this concern by assessing the extent of intermolecular spin diffusion in nineteen peptide-protein complexes of known structure. Transferred NOE intensities were simulated with the program CORONA (Calculated OR Observed NOESY Analysis) using the rate-matrix approach to include contributions from indirect relaxation between protein-ligand and intraligand proton pairs. Intermolecular spin diffusion contributions were determined by comparing intensities calculated with protonated protein to those calculated with fully deuterated protein. The differences were found to be insignificant overall, and to diminish at short mixing times and high mole ratios of peptide to protein. Spin diffusion between the peptide ligand and the protein contributes less to the etNOE intensities and alters fewer cross peaks than the well-studied intramolecular spin diffusion effects. Errors in intraligand interproton distances due to intermolecular relaxation effects were small on average and can be accounted for with the restraint functions commonly used in NMR structure determination methods. In addition, a rate-matrix approach to calculate distances from etNOESY intensities using a volume matrix comprising only intraligand intensities was found to give accurate values. Based on these results, we conclude that structures determined from etNOESY data are no less accurate due to spin diffusion than structures determined from conventional NOE intensities.     

Structure determination of a DNA octamer in solution by NMR spectroscopy. Effect of fast local motions.

NMR structures of biomolecules are primarily based on nuclear Overhauser effects (NOEs) between protons. For the interpretation of NOEs in terms of distances, usually the assumption of a single rotational correlation time corresponding to a rigid molecule approximation is made. Here we investigate the effect of fast internal motions of the interproton vectors in the context of the relaxation matrix approach for structure determination of biomolecules. From molecular dynamics simulations generalized order parameters were calculated for the DNA octamer d(GCGTTCGC).d(CGCAACGC), and these were used in the calculation of NOE intensities. The magnitudes of the order parameters showed some variation for the different types of interproton vectors. The lowest values were observed for the interresidue base H6/H8-H2" proton vectors (S2 = 0.60), while the cytosine H5-H6 interproton vectors were among the most motionally restricted (S2 = 0.92). Inclusion of the motion of the interproton vectors resulted in a much better agreement between theoretically calculated NOE spectra and the experimental spectra measured by 2D NOE spectroscopy. The interproton distances changed only slightly, with a maximum of 10%; nevertheless, the changes were significant and resulted in constraints that were better satisfied. The structure of the DNA octamer was determined by using restrained molecular dynamics simulations with H2O as a solvent, with and without the inclusion of local internal motions. Starting from A- or B-DNA, the structures showed a high local convergence (0.86 A), while the global convergence for the octamer was ca. 2.6 A.     

Refinement of the NMR solution structure of a protein to remove distortions arising from neglect of internal motion.

The effect of internal motion on the quality of a protein structure derived from nuclear magnetic resonance (NMR) cross relaxation has been investigated experimentally. Internal rotation of the tyrosine-31 ring of turkey ovomucoid third domain was found to mediate magnetization transfer; the effect led to underestimation of proton-proton distances in its immediate neighborhood. Experimental methods that distinguish pure cross relaxation from chemical exchange mediated cross relaxation were used to separate true distances from distorted ones. Uncorrected and corrected sets of distances, where the corrections took internal motion into account, each were used as input to a distance geometry program for structural modeling. Each set of distances yielded a family of similar (converged) structures. The two families of structures differed considerably (2 A) in the region of tyrosine-31. In addition, differences as large as 1 A were observed at other positions throughout the structure. These results emphasize the importance of analyzing the effects of internal motions in order to obtain more accurate NMR solution structures.     

Heteronuclear 1H-15N nuclear magnetic resonance studies of the c subunit of the Escherichia coli F1F0 ATP synthase: assignment and secondary structure.

Nuclear magnetic resonance (NMR) studies of the c subunit of F1F0 ATP synthase from Escherichia coli are presented. A combination of homonuclear (1H-1H) and heteronuclear (1H-15N) 2D and 3D methods was applied to the 79-residue protein, dissolved in trifluoroethanol. Resonance assignment for all the backbone amide groups and many C alpha H side-chain protons was achieved. Analysis of inter- and intraresidue 1H-1H nuclear Overhauser effect (NOE) data and scalar coupling constant information indicates that this protein contains two extended regions of predominant alpha-helical character (residues 10-40 and 48-77) separated by an eight-residue segment which displays little evidence of ordered secondary structure. This model is consistent with information about the molecular motion of the protein deduced from 15N-1H heteronuclear NOE data and observed pKa values of carboxylic acid groups.     

Complete relaxation matrix refinement of NMR structures of proteins using analytically calculated dihedral angle derivatives of NOE intensities

Summary A new method for refining three-dimensional (3D) NMR structures of proteins is described, which takes account of the complete relaxation pathways. Derivatives of the NOE intensities with respect to the dihedral angles are analytically calculated, and efficiently evaluated with the use of a filter technique for identifying the dominant terms of these derivatives. This new method was implemented in the distance geometry program DIANA. As an initial test, we refined 30 rigid distorted helical structures, using a simulated data set of NOE distance constraints for a rigid standard -helix. The final root-mean-square deviations of the refined structures relative to the standard helix were less than 0.1 Å, and the R-factors dropped from values between 7% and 32% to values of less than 0.5% in all cases, which compares favorably with the results from distance geometry calculations. In particular, because spin diffusion was not explicitly considered in the evaluation of exact¹H–¹H distances corresponding to the simulated NOE intensities, a group of nearly identical distance geometry structures was obtained which had about 0.5 Å root-mean-square deviation from the standard -helix. Further test calculations using an experimental NOE data set recorded for the protein trypsin inhibitor K showed that the complete relaxation matrix refinement procedure in the DIANA program is functional also with systems of practical interest.Abbreviations RMSD root-mean-square deviation - NOE nuclear Overhauser enhancement - NOESY 2-dimensional nuclear Overhauser enhancement spectroscopy - CPU central processing unit     

The NMR side-chain assignments and solution structure of enzyme IIBcellobiose of the phosphoenolpyruvate-dependent phosphotransferase system of Escherichia coli.

The assignment of the side-chain NMR resonances and the determination of the three-dimensional solution structure of the C10S mutant of enzyme IIBcellobiose (IIBcel) of the phosphoenolpyruvate-dependent phosphotransferase system of Escherichia coli are presented. The side-chain resonances were assigned nearly completely using a variety of mostly heteronuclear NMR experiments, including HCCH-TOCSY, HCCH-COSY, and COCCH-TOCSY experiments as well as CBCACOHA, CBCA(CO)NH, and HBHA(CBCA)(CO)NH experiments. In order to obtain the three-dimensional structure, NOE data were collected from 15N-NOESY-HSQC, 13C-HSQC-NOESY, and 2D NOE experiments. The distance restraints derived from these NOE data were used in distance geometry calculations followed by molecular dynamics and simulated annealing protocols. In an iterative procedure, additional NOE assignments were derived from the calculated structures and new structures were calculated. The final set of structures, calculated with approximately 2000 unambiguous and ambiguous distance restraints, has an rms deviation of 1.1 A on C alpha atoms. IIBcel consists of a four stranded parallel beta-sheet, in the order 2134. The sheet is flanked with two and three alpha-helices on either side. Residue 10, a cysteine in the wild-type enzyme, which is phosphorylated during the catalytic cycle, is located at the end of the first beta-strand. A loop that is proposed to be involved in the binding of the phosphoryl-group follows the cysteine. The loop appears to be disordered in the unphosphorylated state.     

Practical applications of time-averaged restrained molecular dynamics to ligand-receptor systems: FK506 bound to the Q50R,A95H,K98I triple mutant of FKBP-13

Summary The ability of time-averaged restrained molecular dynamics (TARMD) to escape local low-energy conformations and explore conformational space is compared with conventional simulated-annealing methods. Practical suggestions are offered for performing TARMD calculations with ligand-receptor systems, and are illustrated for the complex of the immunosuppressant FK506 bound to Q50R,A95H,K98I triple mutant FKBP-13. The structure of ¹³C-labeled FK506 bound to triple-mutant FKBP-13 was determined using a set of 87 NOE distance restraints derived from HSQC-NOESY experiments. TARMD was found to be superior to conventional simulated-annealing methods, and produced structures that were conformationally similar to FK506 bound to wild-type FKBP-12. The individual and combined effects of varying the NOE restraint force constant, using an explicit model for the protein binding pocket, and starting the calculations from different ligand conformations were explored in detail.Abbreviations DG distance geometry - dmFKBP-12 double-mutant (R42K,H87V) FKBP-12 - FKBP-12 FK506-binding protein (12 kDa) - FKBP-13 FK506-binding protein (13 kDa) - HSQC heteronuclear single-quantum coherence - K_NOE force constant (penalty) for NOE-derived distance restraints - MD molecular dynamics - NOE nuclear Overhauser effect - SA simulated annealing - TARMD molecular dynamics with time-averaged restraints - tmFKBP-13 triple-mutant (Q50R,A95H,K98I) FKBP-13 - wtFKBP-12 wild-type FKBP-12