A high-resolution solid-state 13C-NMR study on crystalline bovine heart cytochrome-c oxidase and lysozyme. Dynamic behavior of protein and detergent in the complex.

We have recorded 100.6-MHz high-resolution solid-state 13C-NMR spectra of crystalline cytochrome-c oxidase from bovine heart muscle and hen egg-white lysozyme, to compare conformation and dynamics of a typical membrane-protein complex with those of lysozyme. The absence of severe interference with the solid-state 13C-NMR spectra, from both the line broadenings from paramagnetic centers and overlapping of intense detergent signals, provided spectral resolution of 13C-NMR feature of cytochrome-c oxidase crystals comparable to that of lysozyme crystal and better than that of dissolved or lyophilized samples. In fact, the observed peak intensities of the polar heads of the detergents BL8SY and Brij 35 were only about 10% and 3% of the anticipated values, respectively. The dynamic behavior of the backbone and side chains of cytochrome-c oxidase was compared with that of lysozyme on the basis of the 13C spin-lattice relaxation times (T1): the backbone of the cytochrome-c oxidase turned out to be more flexible than that of lysozyme. Molecular motions of the detergent molecules attached to the proteins are found to be highly heterogeneous. Detergent molecules undergo rapid tumbling motions in the crystals in about 10 ns as detected by T1. In addition to rapid motions, slow motions were detected by 1H spin-lattice relaxation time in the rotating frame (TH1 rho) and cross-polarization time (TCH), together with data from static spectra, indicating that the aliphatic portion of the detergent interacts more strongly with hydrophobic protein surfaces than do the polar heads.     

Active site dynamics in the lead-dependent ribozyme

Conformational dynamics are an important property of ribozymes and other RNA molecules but there is currently only limited information on the relationship between dynamics and RNA function. A recent structural study of the lead-dependent ribozyme, known as the leadzyme, showed significant dynamics at the active site and indicated that a structural rearrangement is required for the reaction to proceed from the ground to the transition state. In this work, microsecond-to-millisecond dynamics of the leadzyme are probed by analysis of the power dependence of (13)C NMR relaxation times in the rotating frame (T(1)(rho)). These results revealed a wide range of conformational dynamics for various residues in the leadzyme. For residue A25 in the active site, the power dependence of T(1)(rho) yielded an exchange lifetime similar to that previously measured by line-shape analysis, and provides an important calibration of this T(1)(rho) methodology for probing the dynamics of macromolecules. Strong evidence was also found for a previously suggested dynamic network of hydrogen bonds stabilizing the GAAA tetraloop motif. Within the active site of the leadzyme, internal motions are observed on a wide variety of time scales, suggesting a complex landscape of accessible states, and potential correlations between observed motions and catalytic function are discussed. These results demonstrate that the power dependence of (13)C T(1)(rho) relaxation times provides a valuable method for probing dynamics in nucleic acids.     

Dynamic Structure of Proteins in Solid State. 1H and 13C NMR Relaxation Study

Abstract

Temperature dependencies of ¹H non-selective NMR T₁ and T₂ relaxation times measured at two resonance frequencies and natural abundance ^l3C NMR relaxation times T_l and T_lr measured at room temperature have been studied in a set of dry and wet solid proteins—;Bacterial RNase, lysozyme and Bovine serum albumin (BSA). The proton and carbon data were interpreted in terms of a model supposing three kinds of internal motions in a protein. These are rotation of the methyl protons around the axis of symmetry of the methyl group, and fast and slow oscillations of all atoms. The correlation times of these motions in solid state are found around 10^?11, 10^?9 and 10^?6 s, respectively. All kinds of motion are characterized by the inhomogeneous distribution of the correlation times. The protein dehydration affects only the slow internal motion. The amplitude of the slow motion obtained from the carbon data is substantially less than that obtained from the proton data. This difference can be explained by taking into account different relative inter- and intra- chemical group contributions to the proton and carbon second moments. The comparison of the solid state and solution proton relaxation data showed that the internal protein dynamics in these states is different: the slow motion seems to be few orders of magnitude faster in solution.     

Sampling of protein dynamics in nanosecond time scale by 15N NMR relaxation and self-diffusion measurements.

This paper presents a procedure for detection of intermediate nanosecond internal dynamics in globular proteins. The procedure uses 1H-15N relaxation measurements at several spectrometer frequencies and hydrodynamic calculations based on experimental self-diffusion coefficients. New heteronuclear experiments, using pulse field gradients, are introduced for the measurement of translation diffusion coefficients of 15N labeled proteins. An advanced interpretation of recently published (Luginbühl et al., Biochemistry, 36, 7305-7312 (1997)) backbone amide 15N relaxation data, measured at two spectrometers (400 and 750 MHz for 1H) for N-terminal DNA-binding domain (1-63) of 434 repressor, is presented. Non-applicability of commonly used fast (picosecond) dynamics model (FD) was justified by (i) poor fit of relaxation data by the FD model-free spectral density function both for isotropic and anisotropic models of the overall molecular tumbling; (ii) specific dependence of the overall rotation correlation times calculated from T1/T2 ratio on the spectrometer frequency; (iii) mismatch of the ratio of longitudinal 15N relaxation times T1, measured at different spectrometer frequencies, in comparison with that anticipated for the FD model; (iv) significantly underestimated overall rotation correlation time provided by the FD model (5.50+/-0.15 and 5.80+/-0.15 ns for 750 and 400 MHz spectrometer frequency respectively) in comparison with correlation time obtained from hydrodynamics. On the other hand, all relaxation and hydrodynamics data are in good correspondence with the model of intermediate (nanoseconds) dynamics. Overall rotation correlation time of 7.5+/-0.7 ns was calculated from experimental translation self-diffusion rate using hydrodynamics formalism (Garcia de la Torre, J. and Bloomfield, V.A. Quart. Rev. Biophys., 14, 81-139 (1981)). The statistical analysis of 15N relaxation data along with the hydrodynamic consideration clearly revealed that most of the residues in 434(1-63) repressor are involved in the nanosecond internal dynamics characterized by the the mean order parameters of 0.59+/-0.06 and the correlation times of ca. 5 ns.     

Analysis of the relationship between enzyme activity and its internal motion using nuclear magnetic resonance: 15N relaxation studies of wild-type and mutant lysozyme

A mutant lysozyme where R14 and H15 are deleted together has higher activity and a similar binding ability to an inhibitor, trimer of N-acetylglucosamine ((NAG)3), compared with wild-type lysozyme. Since this has been attributed to intrinsic protein dynamic properties, we investigated the relationship between the activity and the internal motions of proteins. Backbone dynamics of the free and the complex forms with the (NAG)3 have been studied by measurement of the 15N T1 and T2 relaxation rates and NOE determinations at 600 MHz. Analysis of the data using the model-free formalism showed that the generalized order parameters (S2) were almost the same in wild-type and mutant lysozyme in unbound state, indicating that the mutation had little effect on the global internal motions. On the other hand, in the presence of (NAG)3, although some signals located around the active site were broadened or decreased in intensity because of strong perturbation by (NAG)3, there were several residues that showed increased or decreased backbone S2 in the complexed lysozymes. A comparison of the internal motions of the wild-type and mutant complexes showed a number of distinct dynamic differences between them. In particular, many residues located at or near active-site regions (turn 1, strand 2, turn 2 and long loop), displayed greater backbone dynamics reflecting the order parameter in mutant complex relative to mutant free. Furthermore, the Rex values at the loop C-D region, which was considered to be important for enzymatic activity, significantly increased. From these results, it was suggested that variations in the dynamics of these regions may play an important role in the enzyme activity.     

Heat capacity and conformation of proteins in the denatured state

Heat capacity, intrinsic viscosity and ellipticity of a number of globular proteins (pancreatic ribonuclease A, staphylococcal nuclease, hen egg-white lysozyme, myoglobin and cytochrome c) and a fibrillar protein (collagen) in various states (native, denatured, with and without disulfide crosslinks or a heme) have been studied experimentally over a broad range of temperatures. It is shown that the partial heat capacity of denatured protein significantly exceeds the heat capacity of native protein, especially in the case of globular proteins, and is close to the value calculated for an extended polypeptide chain from the known heat capacities of individual amino acid residues. The significant residual structure that appears at room temperature in the denatured states of some globular proteins (e.g. myoglobin and lysozyme) at neutral pH results in a slight decrease of the heat capacity, probably due to partial screening of the protein non-polar groups from water. The heat capacity of the unfolded state increases asymptotically, approaching a constant value at about 100 degrees C. The temperature dependence of the heat capacity of the native state, which can be determined over a much shorter range of temperature than that of the denatured state and, correspondingly, is less certain, appears to be linear up to 80 degrees C. Therefore, the denaturational heat capacity increment seems to be temperature-dependent and is likely to decrease to zero at about 140 degrees C.     

Nuclear magnetic resonance studies of amino acids and proteins. Rotational correlation times of proteins by deuterium nuclear magnetic resonance spectroscopy

We show that measurement of the spin-lattice (T1) and spin-spin (T2) relaxation times (or line widths) of irrotationally bound 2H nuclei in macromolecules undergoing isotropic rotational motion outside of the extreme narrowing limit (i.e., for the case omega 02 tau R2 much greater than 1) permits determination of both the rotational correlation time (tau R) of the macromolecule and the electric quadrupole coupling constant (e2qQ/h) of the 2H label. The technique has the advantage over 13C nuclear magnetic resonance (NMR) that no assumptions about bond lengths (which appear to the sixth power in 13C relaxation studies) or relaxation mechanisms need to be made, since relaxation will always be quadrupolar, even for aromatic residues at high field. Asymmetry parameter (eta) uncertainties are shown to cause negligible effects on tau R determinations, and in any case it is shown that both e2qQ/h and eta may readily be determined in separate solid-state experiments. By way of example, we report 2H NMR results on aqueous lysozyme (EC 3.2.1.17) at 5.2 and 8.5 T (corresponding to 2H-resonance frequencies of 34 and 55 MHz). Interpretation of the results in terms of the isotropic rigid-rotor model yields e2qQ/h values of approximately equal to 170 or approximately equal to 190 kHz, respectively, for the imidazolium and free-base forms of [epsilon 1-2H] His-15 lysozyme in solution, in excellent agreement with e2qQ/h values of approximately 167 and approximately 190 kHz obtained for the free amino acids in the solid state. In principle, the method may in suitable cases permit comparison between the dynamic structures of proteins in solution and in the crystalline solid state.     

Low-frequency motion in membranes. The effect of cholesterol and proteins

Nuclear magnetic resonance (NMR) relaxation techniques have been used to study the effect of lipid-protein interactions on the dynamics of membrane lipids. Proton enhanced (PE) 13C-NMR measurements are reported for the methylene chain resonances in red blood cell membranes and their lipid extracts. For comparison similar measurements have been made of phospholipid dispersions containing cholesterol and the polypeptide gramicidin A+. It is found that the spin-lattice relaxation time in the rotating reference frame (T1 rho) is far more sensitive to protein, gramicidin A+ or cholesterol content than is the laboratory frame relaxation time (T1). Based on this data it is concluded that the addition of the second component to a lipid bilayer produces a low-frequency motion in the region of 10(5) to 10(7) Hz within the membrane lipid. The T1 rho for the superimposed resonance peaks derived from all parts of the phospholipid chain are all influenced in the same manner suggesting that the low frequency motion involves collective movements of large segments of the hydrocarbon chain. Because of the molecular co-operativity implied in this type of motion and the greater sensitivity of T1 rho to the effects of lipid-protein interactions generally, it is proposed that these low-frequency perturbations are felt at a greater distance from the protein than those at higher frequencies which dominate T1.     

Studies of individual carbon sites of proteins in solution by natural abundance carbon 13 nuclear magnetic resonance spectroscopy. Relaxation behavior.

The aromatic regions in proton-decoupled natural abundance 13C Fourier transform nuclear magnetic resonance spectra (at 14.2 kG) of small native proteins contain broad methine carbon bands and narrow nonprotonated carbon resonances. Some factors that affect the use of natural abundance 13C Fourier transform NMR spectroscopy for monitoring individual nonprotonated aromatic carbon sites of native proteins in solution are discussed. The effect of protein size is evaluated by comparing the 13C NMR spectra of horse heart ferrocytochrome c, hen egg white lysozyme, horse carbon monoxide myoglobin, and human adult carbon monoxide hemoglobin. Numerous single carbon resonances are observed in the aromatic regions of 13C NMR spectra of cytochrome c, lysozyme, and myoglobin. The much larger hemoglobin yields few resolved individual carbon resonances. Theoretical and some experimental values are presented for the natural linewidths (W), spin-lattice relaxation times (T1), and nuclear Overhauser enhancements (NOE) of nonprotonated aromatic carbons and Czeta of arginine residues. In general, the 13C-1H dipolar mechanism dominates the relaxation of these carbons. 13C-14N dipolar relaxation contributes significantly to 1/T1 of C epsilon2 of tryptophan residues and Czeta of arginine residues of proteins in D2O. The NOE of each nonprotonated aromatic carbon is within experimental error of the calculated value of about 1.2. As a result, integrated intensities can be used for making a carbon count. Theoretical results are presented for the effect of internal rotation on W, T1, and the NOE. A comparison with the experimental T1 and NOE values indicates that if there is internal rotation of aromatic amino acid side chains, it is not fast relative to the over-all rotational motion of the protein.     

Mechanisms for Benzene Dissociation through the Excited State of T4 Lysozyme L99A Mutant

The atomic-level mechanisms that coordinate ligand release from protein pockets are only known for a handful of proteins. Here, we report results from accelerated molecular dynamics simulations for benzene dissociation from the buried cavity of the T4 lysozyme Leu99Ala mutant (L99A). In these simulations, benzene is released through a previously characterized, sparsely populated room-temperature excited state of the mutant, explaining the coincidence for experimentally measured benzene off rate and apo protein slow-timescale NMR relaxation rates between ground and excited states. The path observed for benzene egress is a multistep ligand migration from the buried cavity to ultimate release through an opening between the F/G-, H-, and I-helices and requires a number of cooperative multiresidue and secondary-structure rearrangements within the C-terminal domain of L99A. These rearrangements are identical to those observed along the ground state to excited state transitions characterized by molecular dynamic simulations run on the Anton supercomputer. Analyses of the molecular properties of the residues lining the egress path suggest that protein surface electrostatic potential may play a role in the release mechanism. Simulations of wild-type T4 lysozyme also reveal that benzene-egress-associated dynamics in the L99A mutant are potentially exaggerations of the substrate-processivity-related dynamics of the wild type.     

Influence of hydration on the dynamics of lysozyme

Quasielastic neutron and light-scattering techniques along with molecular dynamics simulations were employed to study the influence of hydration on the internal dynamics of lysozyme. We identified three major relaxation processes that contribute to the observed dynamics in the picosecond to nanosecond time range: 1), fluctuations of methyl groups; 2), fast picosecond relaxation; and 3), a slow relaxation process. A low-temperature onset of anharmonicity at T approximately 100 K is ascribed to methyl-group dynamics that is not sensitive to hydration level. The increase of hydration level seems to first increase the fast relaxation process and then activate the slow relaxation process at h approximately 0.2. The quasielastic scattering intensity associated with the slow process increases sharply with an increase of hydration to above h approximately 0.2. Activation of the slow process is responsible for the dynamical transition at T approximately 200 K. The dependence of the slow process on hydration correlates with the hydration dependence of the enzymatic activity of lysozyme, whereas the dependence of the fast process seems to correlate with the hydration dependence of hydrogen exchange of lysozyme.     

Enzyme-inhibitor complexes of lysozyme with glucosamine inhibitors. A molecular dynamics study through 2H-NMR

2H relaxation measurements coupled with multiple specific 2H labeling have provided insight into the molecular dynamics of N-acetyl-D-glucosamine (GlcNAc) inhibitors bound to lysozyme. Deuteron T1 and T2 data for the bound state of methyl alpha- and -beta-GlcNAc 2H-labeled in the glycosidic methyl and C2 positions have been derived from measurements at different enzyme/inhibitor ratios. Rotational correlation times calculated therefrom for the labeled sites indicate, in both cases, tight binding for the sugar ring (tau(b) = 3.0 x 10(-9) s) accompanied by fast internal rotation, about one axis, of the glycosidic methyl groups (tau(r) = 5.5-7.6 x 10(-11) s). The small but consistent difference in the rates of internal rotation for the alpha- and beta-anomeric inhibitors may be indicative of different solution structures of the enzyme-inhibitor complexes.     

Imazalil-cyclomaltoheptaose (beta-cyclodextrin) inclusion complex: preparation by supercritical carbon dioxide and 13C CPMAS and 1H NMR characterization

An inclusion complex between imazalil (IMZ), a selected fungicide, and cyclomaltoheptaose (beta-cyclodextrin, betaCD) was obtained using supercritical fluid carbon dioxide. The best preparation conditions were determined, and the inclusion complex was investigated by means of 1H NMR spectroscopy in aqueous solution and 13C CPMAS NMR spectroscopy in the solid state. Information on the geometry of the betaCD/IMZ complex was obtained from ROESY spectroscopy, while the dynamics of the inclusion complex in the kilohertz range was obtained from the proton spin-lattice relaxation times in the rotating frame, T(1rho) (1H).     

Characterization of the wound-induced material in Citrus paradisi fruit peel by carbon-13 CP-MAS solid state NMR spectroscopy

Grapefruit, Citrus paradisi, were injured, inoculated with Penicillium digitatum and incubated under conditions favourable for the accumulation of defence related material. Histochemical examination revealed that tissues adjacent to inoculated injuries contained phloroglucinol-HCl (PG-HCl) reactive material. Solvent washed cell wall preparations of intact and injured-inoculated peel were further purified using a mixture of cell wall degrading enzymes. Samples from injured inoculated tissue contained PG-HCl reactive globular material in addition to the fragments of xylem and cuticle found in controls. The principal chemical moieties of the material that accumulates in grapefruit injuries during wound-healing were studied by solid state 13C cross-polarization magic angle spinning NMR. A complete assignment of the NMR signals was made. From the analysis evidence was found that cellulose and hemicellulose are the biopolymers present in the intact peel samples, in addition, relevant quantities of cutin were found in the residues of enzyme digest. The NMR difference spectrum intact- wounded peels showed resonances which were attributed to all major functional groups of the aromatic-aliphatic suberin polyester of new material produced by the wounds. Information on the latter polyester was obtained by analyzing the T(1)rho (1H) relaxation.     

Hydration dependence of backbone and side chain polylysine dynamics: a 13C solid-state NMR and IR spectroscopy study

The molecular dynamics of solid poly-L-lysine has been studied by the following natural abundance (13)C-NMR relaxation methods: measurements of the relaxation times T(1) at two resonance frequencies, off-resonance T(1rho) at two spin-lock frequencies, and proton-decoupled T(1rho). Experiments were performed at different temperatures and hydration levels (up to 17% H(2)O by weight). The natural abundance (13)C-CPMAS spectrum of polylysine provides spectral resolution of all types of backbone and side chain carbons and thus, dynamic parameters could be determined separately for each of them. At the same time, the conformational properties of polylysine were investigated by Fourier transform infrared spectroscopy. The data obtained from the different NMR experiments were simultaneously analyzed using the correlation function formalism and model-free approach. The results indicate that in dry polylysine both backbone and side chains take part in two low amplitude motions with correlation times of the order of 10(-4) s and 10(-9) s. Upon hydration, the dynamic parameters of the backbone remain almost constant except for the amplitude of the slower process that increases moderately. The side chain dynamics reveals a much stronger hydration response: the amplitudes of both slow and fast motions increase significantly and the correlation time of the slow motion shortens by about five orders of magnitude, and at hydration levels of more than 10% H(2)O fast and slow side chain motions are experimentally indistinguishable. These changes in the molecular dynamics cannot be ascribed to any hydration-dependent conformational transitions of polylysine because IR spectra reveal almost no hydration dependence in either backbone or side chain absorption domains. The physical nature of the fast and slow motions, their correlation time distributions, and hydration dependence of microdynamic parameters are discussed.     

Activities of daily living influence tibial cartilage T1rho relaxation times

Quantitative T1rho magnetic resonance imaging (MRI) can potentially help identify early-stage osteoarthritis (OA) by non-invasively assessing proteoglycan concentration in articular cartilage. T1rho relaxation times are negatively correlated with proteoglycan concentration. Cartilage compresses in response to load, resulting in water exudation, a relative increase in proteoglycan concentration, and a decrease in the corresponding T1rho relaxation times. To date, there is limited information on changes in cartilage composition resulting from daily activity. Therefore, the objective of this study was to quantify changes in tibial cartilage T1rho relaxation times in healthy human subjects following activities of daily living. It was hypothesized that water exudation throughout the day would lead to decreased T1rho relaxation times. Subjects underwent MR imaging in the morning and afternoon on the same day and were free to go about their normal activities between scans. Our findings confirmed the hypothesis that tibial cartilage T1rho relaxation times significantly decreased (by 7%) over the course of the day with loading, which is indicative of a relative increase in proteoglycan concentration. Additionally, baseline T1rho values varied with position within the cartilage, supporting a need for site-specific measurements of T1rho relaxation times. Understanding how loading alters the proteoglycan concentration in healthy cartilage may hold clinical significance pertaining to cartilage homeostasis and potentially help to elucidate a mechanism for OA development. These results also indicate that future studies using T1rho relaxation times as an indicator of cartilage health should control the loading history prior to image acquisition to ensure the appropriate interpretation of the data.     

The structural dynamics of myoglobin

Conformational fluctuations in proteins were initially invoked to explain the observation that diffusion of small ligands through the matrix is a global phenomenon. Small globular proteins contain internal cavities that play a role not only in matrix dynamics but also in controlling function, tracing a pathway for the diffusion of the ligand to and from the active site. This is the main point addressed in this Review, which presents pertinent information obtained on myoglobin (Mb). Mb, a simple globular heme protein which binds reversibly oxygen and other ligands. The bond between the heme Fe(II) and gaseous ligands can be photodissociated by a laser pulse, generating a non-equilibrium population of protein structures that relaxes on a picosecond to millisecond time range. This process is associated with migration of the ligand to internal cavities of the protein, which are known to bind xenon. Some of the results obtained by laser photolysis, molecular dynamics simulations, and X-ray diffraction of intermediate states of wild-type and mutant myoglobins are summarized. The extended relaxation of the globin moiety directly observed by Laue crystallography reflects re-equilibration among conformational substates known to play an essential role in controlling protein function.     

Solid state NMR and X-ray diffraction studies of alpha-D-galacturonic acid monohydrate

Crystalline alpha-D-galacturonic acid monohydrate has been studied by 13C CPMAS NMR and X-ray crystallography. The molecular dynamics were investigated by evaluating 13C spin-lattice relaxation in the rotating frame (T1rho) and chemical-shift-anisotropy properties of each carbon. Only limited molecular motions can be detected in the low frequency (< 10(4) Hz) range by 13C relaxation time measurements (T1rho) and changes of chemical shift anisotropy properties as a function of temperature. X-ray analysis (at both ambient temperature and 150 K) shows that the acid has the usual chair-shaped, pyranose ring conformation, and that the acid and water molecules are linked, through all their O-H groups, in an extensively hydrogen-bonded lattice.     

Conformational dynamics of yeast calmodulin in the Ca(2+)-bound state probed using NMR relaxation dispersion

Most calmodulin (CaM) in apo and Ca(2+)-bound states show a dumb-bell-like structure, involving the N- and C-terminal domains, connected with a flexible linker. However, Ca(2+)-bound yeast calmodulin (yCaM) takes on a unique globular structure; the target-binding site of this protein is autoinhibited. We applied NMR relaxation dispersion experiments to yCaM in the Ca(2+)-bound state. The amide (15)N and (1)H(N) relaxation dispersion profiles indicated the presence of conformational dynamics for specific residues at the interface between the N- and C-terminal domains. We conclude that these conformational dynamics were derived from the mobility of the C-terminal domain.