Cross-correlated Rates Research Articles

DNA base order parameter determination without influence of chemical exchange.

Nuclear magnetic resonance (NMR) spectroscopy is a versatile tool used to investigate the dynamic properties of biological macromolecules and their complexes. NMR relaxation data can provide order parameters S2, which represent the mobility of bond vectors reorienting within a molecular frame. Determination of S2 parameters typically involves the use of transverse NMR relaxation rates. However, the accuracy in S2 determination can be diminished by elevation of the transverse relaxation rates through conformational or chemical exchange involving protonation/deprotonation or non-Watson-Crick base-pair states of nucleic acids. Here, we propose an approach for determination of S2 parameters without the influence of exchange processes. This approach utilizes transverse and longitudinal 13C chemical shift anisotropy (CSA) - dipole-dipole (DD) cross-correlation rates instead of 13C transverse relaxation rates. Anisotropy in rotational diffusion is taken into consideration. An application of this approach to nucleotide base CH groups of a uniformly 13C/15N-labeled DNA duplex is demonstrated.

Mobility of Histidine Side Chains Analyzed with 15N NMR Relaxation and Cross-Correlation Data: Insight into Zinc-Finger-DNA Interactions.

Due to chemical exchange, the mobility of histidine (His) side chains of proteins is typically difficult to analyze by NMR spectroscopy. Using an NMR approach that is uninfluenced by chemical exchange, we investigated internal motions of the His imidazole NH groups that directly interact with DNA phosphates in the Egr-1 zinc-finger-DNA complex. In this approach, the transverse and longitudinal cross-correlation rates for 15N chemical shift anisotropy and 15N-1H dipole-dipole relaxation interference were analyzed together with 15N longitudinal relaxation rates and heteronuclear Overhauser effect data at two magnetic field strengths. We found that the zinc-coordinating His side chains directly interacting with DNA phosphates are strongly restricted in mobility. This makes a contrast to the arginine and lysine side chains that retain high mobility despite their interactions with DNA phosphates in the same complex. The entropic effects of side-chain mobility on the molecular association are discussed.

Determining the Dopaminergic Feedback Pathway in Pancreatic β-Cells with Fluorescence Fluctuation Spectroscopy

Tight regulation of insulin and glucagon, two hormones secreted by the pancreatic islet, allows the body to maintain glucose homeostasis. Insulin resistance coupled with insufficient insulin secretion leads to type 2 diabetes. Insulin secretion from pancreatic β-cells is regulated by multiple signaling inputs including the neurotransmitter dopamine. Upon treatment with dopamine, the amplitude and frequency of intracellular free calcium oscillations in β-cells decrease, leading to a concomitant decrease in insulin section. We hypothesize that activation of the dopamine receptor (DRD3) by dopamine releases the Gβγ complex to stimulate a G protein-coupled inwardly-rectifying potassium channel (GIRK3). This inward current keeps the cell polarized, and reduces currents through voltage-gated calcium channels (CaV1.2). This model predicts numerous protein interactions that may only be transient in nature, in particular, we expect close interactions between the DRD3 and GIRK3 as well as each of these with the Gβγ subunit. To test this hypothesis, we are utilizing fluorescence fluctuation spectroscopy to determine protein interactions in stable β-cell lines. Fluorescence fluctuation spectroscopy provides a method to determine interaction between two proteins in a live cell that surpasses many of the limitations of FRET. Single and two color studies allow us to determine diffusion constants of each protein before and after dopamine stimulation and the cross-correlation, or amount of interaction, between two proteins. Additionally, brightness analysis can be used to determine homo- and heterodimerization of the tagged proteins. Diffusion constants, cross-correlation rates, and protein clustering among the DRD3, GIRK3, and Gβγ subunit will be presented. Through this analysis, the dopamine signaling pathway in pancreatic β-cells can be determined, which will reveal novel potential therapeutic targets that can increase insulin secretion during diabetes mellitus.

Cytochrome-P450–Cytochrome-b5 Interaction in a Membrane Environment Changes 15N Chemical Shift Anisotropy Tensors

It has been well realized that the dependence of chemical shift anisotropy (CSA) tensors on the amino acid sequence, secondary structure, dynamics, and electrostatic interactions can be utilized in the structural and dynamic studies of proteins by NMR spectroscopy. In addition, CSA tensors could also be utilized to measure the structural interactions between proteins in a protein-protein complex. To this end, we report the experimentally measured backbone amide-(15)N CSA tensors for a membrane-bound 16.7 kDa full-length rabbit cytochrome-b5 (cytb5), in complexation with a 55.8 kDa microsomal rabbit cytochrome P450 2B4 (cytP4502B4). The (15)N-CSAs, determined using the (15)N CSA/(15)N-(1)H dipolar coupling transverse cross-correlated rates, for free cytb5 are compared with those for the cytb5 bound to cytP4502B4. An overall increase in backbone amide-(15)N transverse cross-correlated rates for the cytb5 residues in the cytb5-cytP450 complex is observed as compared to the free cytb5 residues. Due to fast spin-spin relaxation (T2) and subsequent broadening of the signals in the complex, we could measure amide-(15)N CSAs only for 48 residues of cytb5 as compared to 84 residues of free cytb5. We observed a change in (15)N CSA for most residues of cytb5 in the complex, as compared to free cytb5, suggesting a dynamic interaction between the oppositely charged surfaces of anionic cytb5 and cationic cytP450. The mean values of (15)N CSA determined for residues in helical, sheet, and turn regions of cytb5 in the complex are -184.5, -146.8, and -146.2 ppm, respectively, with an overall average value of -165.5 ppm (excluding the values from residues in more flexible termini). The measured CSA value for residues in helical conformation is slightly larger as compared to previously reported values. This may be attributed to the paramagnetic effect from Fe(III) of the heme in cytb5, which is similar to our previously reported values for the free cytb5.

Probing Local Backbone Geometries in Intrinsically Disordered Proteins by Cross-Correlated NMR Relaxation.

Probing Local Backbone Geometries in Intrinsically Disordered Proteins by Cross-Correlated NMR Relaxation.

Determination of 15N chemical shift anisotropy from a membrane-bound protein by NMR spectroscopy.

Chemical shift anisotropy (CSA) tensors are essential in the structural and dynamic studies of proteins using NMR spectroscopy. Results from relaxation studies in biomolecular solution and solid-state NMR experiments on aligned samples are routinely interpreted using well-characterized CSA tensors determined from model compounds. Since CSA tensors, particularly the (15)N CSA, highly depend on a number of parameters including secondary structure, electrostatic interaction, and the amino acid sequence, there is a need for accurately determined CSA tensors from proteins. In this study, we report the backbone amide-(15)N CSA tensors for a 16.7-kDa membrane-bound and paramagnetic-heme containing protein, rabbit Cytochrome b(5) (cytb(5)), determined using the (15)N CSA/(15)N-(1)H dipolar transverse cross-correlation rates. The mean values of (15)N CSA determined for residues in helical, sheet, and turn regions are -187.9, -166.0, and -161.1 ppm, respectively, with an overall average value of -171.7 ppm. While the average CSA value determined from this study is in good agreement with previous solution NMR experiments on small globular proteins, the CSA value determined for residues in helical conformation is slightly larger, which may be attributed to the paramagnetic effect from Fe(III) of the heme unit in cytb(5). However, like in previous solution NMR studies, the CSA values reported in this study are larger than the values measured from solid-state NMR experiments. We believe that the CSA parameters reported in this study will be useful in determining the structure, dynamics, and orientation of proteins, including membrane proteins, using NMR spectroscopy.

Temperature dependence of fast carbonyl backbone dynamics in chicken villin headpiece subdomain

Temperature-dependence of protein dynamics can provide information on details of the free energy landscape by probing the characteristics of the potential responsible for the fluctuations. We have investigated the temperature-dependence of picosecond to nanosecond backbone dynamics at carbonyl carbon sites in chicken villin headpiece subdomain protein using a combination of three NMR relaxation rates: (13)C' longitudinal rate, and two cross-correlated rates involving dipolar and chemical shift anisotropy (CSA) relaxation mechanisms, (13)C'/(13)C'-(13)C(α) CSA/dipolar and (13)C'/(13)C'-(15)N CSA/dipolar. Order parameters have been extracted using the Lipari-Szabo model-free approach assuming a separation of the time scales of internal and molecular motions in the 2-16°C temperature range. There is a gradual deviation from this assumption from lower to higher temperatures, such that above 16°C the separation of the time scales is inconsistent with the experimental data and, thus, the Lipari-Szabo formalism can not be applied. While there are variations among the residues, on the average the order parameters indicate a markedly steeper temperature dependence at backbone carbonyl carbons compared to that probed at amide nitrogens in an earlier study. This strongly advocates for probing sites other than amide nitrogen for accurate characterization of the potential and other thermodynamics characteristics of protein backbone.

Creation and evolution of net proton hyperpolarization arising from para-hydrogenation

When a hydrogenation reaction is carried out with gaseous hydrogen enriched in its para- isomer in the earth magnetic field (prior to adiabatic insertion of the sample in the NMR magnet), enhanced proton longitudinal order (represented by 2 I z A I z B ) is created but also difference of enhanced polarizations ( I z A - I z B ). In a first part, it is shown theoretically and experimentally that the longitudinal relaxation time of this polarization difference is roughly twice the ones of individual polarizations. The second part is devoted to a pulse sequence designed for transforming this difference into net hyperpolarization. The evolution of this global hyperpolarization is studied experimentally in a third part and it is observed that a fraction of hyperpolarization possesses an effective longitudinal relaxation time similar to the one of the initial polarization difference. Those experimental results are interpreted by numerical calculations based on Solomon-type equations including the longitudinal order and possibly dipolar-csa cross correlation rates.

Comparison of fast backbone dynamics at amide nitrogen and carbonyl sites in dematin headpiece C-terminal domain and its S74E mutant

We perform a detailed comparison of fast backbone dynamics probed at amide nitrogen versus carbonyl carbon sites for dematin headpiece C-terminal domain (DHP) and its S74E mutant (DHPS74E). Carbonyl dynamics is probed via auto-correlated longitudinal rates and transverse C'/C'-C(alpha) CSA/dipolar and C'/C'-N CSA/dipolar cross-correlated rates, while (15)N data are taken from a previous study. Resulting values of effective order parameters and internal correlation times support the conclusion that C' relaxation reports on a different subset of fast motions compared to those probed at N-H bond vectors in the same peptide planes. (13)C' order parameters are on the average 0.08 lower than (15)N order parameters with the exception of the flexible loop region in DHP. The reduction of mobility in the loop region upon the S74E mutation can be seen from the (15)N order parameters but not from the (13)C order parameters. Internal correlation times at (13)C' sites are on the average an order of magnitude longer than those at (15)N sites for the well-structured C-terminal subdomains, while the more flexible N-terminal subdomains have more comparable average internal correlation times.

Probing Side-Chain Dynamics in the Proteasome by Relaxation Violated Coherence Transfer NMR Spectroscopy

A pair of experiments is presented for measuring intra-methyl 1H-1H dipolar cross-correlated spin relaxation rates in highly deuterated, methyl protonated proteins with significantly improved sensitivity relative to previously developed experiments that measure dynamics via 1H spin relaxation. In applications to proteins with correlation times in the macromolecular limit, these cross-correlation rates are related directly to order parameters, characterizing the amplitude of motion of methyl-containing side-chains. The experimental approach is validated by comparing extracted order parameters with those obtained via 2H and 13C spin relaxation methods for both protein L (7.5 kDa) and malate synthase G (82 kDa), with excellent correlations obtained. The methodology is applied to study Ile, Leu, and Val side-chain dynamics in a 360 kDa "half-proteasome" complex. In particular, order parameters obtained from the WT complex and from a second complex where the proteasome gating residues are deleted establish that the relative levels of dynamics in each of the two molecules are very similar. It thus becomes clear that there is no communication between gating residues and other regions of the molecule involving pico- to nanosecond time-scale dynamics of these methyl-containing side-chains.

NMR structures of paramagnetic metalloproteins

Metalloproteins represent a large share of the proteome and many of them contain paramagnetic metal ions. The knowledge, at atomic resolution, of their structure in solution is important to understand processes in which they are involved, such as electron transfer mechanisms, enzymatic reactions, metal homeostasis and metal trafficking, as well as interactions with their partners. Formerly considered as unfeasible, the first structure in solution by nuclear magnetic resonance (NMR) of a paramagnetic protein was obtained in 1994. Methodological and instrumental advancements pursued over the last decade are such that NMR structure of paramagnetic proteins may be now routinely obtained. We focus here on approaches and problems related to the structure determination of paramagnetic proteins in solution through NMR spectroscopy. After a survey of the background theory, we show how the effects produced by the presence of a paramagnetic metal ion on the NMR parameters, which are in many cases deleterious for the detection of NMR spectra, can be overcome and turned into an additional source of structural restraints. We also briefly address features and perspectives given by the use of 13C-detected protonless NMR spectroscopy for proteins in solution. The structural information obtained through the exploitation of a paramagnetic center are discussed for some Cu2+ -binding proteins and for Ca2+ -binding proteins, where the replacement of a diamagnetic metal ion with suitable paramagnetic metal ions suggests novel approaches to the structural characterization of proteins containing diamagnetic and NMR-silent metal ions.

Anisotropic Local Motions and Location of Amide Protons in Proteins

A new method has been developed to obtain dynamic and structural information about peptide planes in proteins by a combination of measurements of weak short-range cross-correlation rates R(H(N)N/NC') that are due to concerted fluctuations of the H(N)-N and N-C' dipole-dipole interactions and stronger long-range cross-correlation rates R(C'H(N)/H(N)N) and R(NH(N)/H(N)C(alpha)). The rates were interpreted using the axially symmetric Gaussian axial fluctuation model (GAF). The oscillation amplitudes as well as the positions of H(N) atoms with respect to peptide planes in ubiquitin were determined. Most N-H(N) bonds were found not to lie exactly along the bisector of the N-C' and N-C(alpha) bonds but to be slightly tilted toward the carbon-terminal side of the peptide.

Temperature-dependent protein backbone dynamics from auto- and cross-correlated NMR relaxation rates

The temperature dependence of nuclear magnetic resonance relaxation rates was investigated for the backbone of15N/13C labeled human ubiquitin in the temperature range of 20–50 °C. The15N autorelaxation rates give evidence that the potential energy functions for15N−HN bonds are not quadratic, in agreement with results for other proteins. Cross-correlation rates arising from correlated fluctuations of two15N−HN dipole-dipole interactions involving successive residues were obtained by the method of Pelupessy et al. (P. Pelupessy, S. Ravindranathan, G. Bodenhausen: J. Biomol. NMR 25, 265–280, 2003). The results suggest the presence of slow internal motions at 50 °C.

Spin-state selection for increased confidence in cross-correlation rates measurements

A new approach is described for measuring chemical shift anisotropy (CSA)/dipolar cross-correlated relaxation (CCR) rates based on the selection of the individual 15N doublet components prior to the relaxation period. The method uses the spin-state-selective element (S3E) of Sørensen and co-authors [Meissner et al. (1997) J. Mag. Reson., 128, 92-97]. The main advantage of the new method compared to other J-resolved experiments is that it does not create problems of additional signal overlap encountered in coupled spectra. At the same time, this approach allows a simpler control of magnetization pathways than the indirect methods. The method is demonstrated for the B3 domain of protein G.

Limited Variations in 15N CSA Magnitudes and Orientations in Ubiquitin Are Revealed by Joint Analysis of Longitudinal and Transverse NMR Relaxation

The site-specific magnitudes and orientations of the chemical shift tensors have been estimated for 70 backbone (15)N-nuclei in human ubiquitin from the field dependence of dynamic independent ratios between relaxation rates, both longitudinal and transverse, measured at 9.4, 11.7, 14.1, and 18.8 T. The results were jointly analyzed with previously published relaxation data [Fushman; Tjandra; Cowburn. J.Am. Chem. Soc. 1998, 120, 10947-10952] [Kover; Batta. J. Mag. Reson. 2001, 150, 137-146]. The effective magnitudes of the anisotropies distribute around 169 ppm with a variability of 5 ppm. The orientation factors, reflecting the orientation of the CSA relative to the NH bond, distribute around -0.80 with a variability of 0.04, which corresponds to an angle between the symmetry axis of an assumed axially symmetric shielding tensor and the NH bond of 21.4 degrees, and a variability of 2.3 degrees. Correlations with the isotropic (15)N-chemical shifts are observed. Variations in the shielding anisotropies add uncertainty to the obtained order parameters proportional to the square of the magnetic field, when data are analyzed using an assumed invariant CSA tensor for all sites. Around 3% additional uncertainty in the order parameters for 800 MHz data is expected. The optimal TROSY field for amide nitrogen TROSY is estimated, with only marginal variations due to site-to-site variations. Variations in the shielding tensors add uncertainty to the exchange terms calculated from cross-correlation rates. An approach for estimating the exchange terms is suggested, where the uncertainty due to CSA-variations is minimized.

HMBC-like experiment based on longitudinal csa/dipolar cross-correlation.

A gradient-based sequence is proposed for efficiently filtering out all quantities except the longitudinal two-spin order, created by csa/dipolar cross-correlation rates. The dipolar interaction is between a proton and a heteronucleus, the csa being generally the one of the heteronucleus. Proton detection is carried out to benefit from maximum sensitivity. The resulting two-dimensional spectrum has the same aspect as an HMBC spectrum, implying the existence of a J coupling between the two considered nuclei, but here cross-peaks indicate the strength of the relevant csa/dipolar cross-correlation rate. The method is especially interesting in the case of medium-sized molecules where a given heteronucleus is subjected to several csa/dipolar cross-correlation effects, along with the corresponding J couplings.

Measurement of15N csa/dipolar cross-correlation rates by means of Spin State Selective experiments

We propose a method for the determination of (15)N csa/dipolar cross-correlation rates based on the measurement of the two apparent transverse (or longitudinal) relaxation rates associated with each component of the nitrogen doublet (N(alpha) and N(beta)). This is achieved by inserting a spin state selective scheme in conventional inverse Carr-Purcell-Meiboom-Gill (or inversion-recovery) pulse sequence which allows for the edition of a HSQC-type spectrum for each of the spin states. Transverse cross-correlation rates necessitate two independent sets of measurements (for N(alpha) and N(beta), respectively), whereas for longitudinal cross correlation rates, besides N(alpha) and N(beta) measurements, the method requires the knowledge of both the (15)N longitudinal auto-relaxation rate and the longitudinal two-spin order (2NzHz) auto-relaxation rate. These additional parameters are mandatory because of the non-exponential behavior of the N(alpha) and N(beta) longitudinal decays. Conversely, the present method does not require any complex manipulation of 2D spectra, the cross-correlation rates being obtained from the difference of the two (N(alpha) and N(beta)) apparent relaxation rates. This approach is applied to (15)N-labelled ubiquitin at two different magnetic fields (9.4 T and 14.1 T).

Measurement of dipolar cross-correlation in methylene groups in uniformly 13C-, 15N-labeled proteins.

A CC(CO)NH TOCSY-based 3D pulse scheme is presented for measuring (1)H-(13)C dipole-dipole cross-correlated relaxation at CH(2) positions in uniformly (13)C-, (15)N-labeled proteins. Simulations based on magnetization evolution under relaxation and scalar coupling interactions show that cross-correlation rates between (1)H-(13)C dipoles in CH(2) groups can be simply obtained from the intensities of (13)C triplets. The normalized cross-correlation relaxation rates are related to cross-correlation order parameters for macromolecules undergoing isotropic motion, which reflect the degrees of spatial restriction of CH(2) groups. The study on human intestinal fatty acid binding protein (131 residues) in the presence of oleic acid demonstrates that side chain dynamics at most CH(2) positions can be characterized for proteins less than 15 kDa in size, with the proposed TOCSY-based approach.

New routes to the detection of relaxation allowed coherence transfer in paramagnetic molecules.

A new pulse sequence is proposed to measure cross correlation rates between 1H Curie spin relaxation and 1H-15N or 1H-13C dipole-dipole coupling (%@mt;sys@%Gamma%@sx@%H,HX%@be@%CS,DD%@sxx@%%@mx@% ) in paramagnetic systems. The new sequence has been conceived to obtain quantitative measurements of cross correlation rates in the close proximity to the paramagnetic center, preventing the loss of information due to fast transverse relaxation. The approach was tested on the dicalcium protein calbindin D9k in which Ca2+ has been substituted at site II with Ce3+. At variance with previously reported experiments, all peaks observed in HSQC experiments tailored to paramagnetic signals give quantitative estimates of %@mt;sys@%Gamma%@sx@%H%@ital@%,%@rsf@%HX%@be@%CS%@ital@%,%@rsf@%DD%@sxx@%%@mx@% . This is crucial to refine the immediate proximity of the metal ion.

Monitoring the early steps of unfolding of dicalcium and mono-Ce3+-substituted forms of P43M calbindin D9k.

Early steps of unfolding of P43M Calbindin D(9k) have been evaluated by NMR spectroscopy on the native dicalcium and on the paramagnetic monocerium-substituted derivative. Although at 2 M GdmHCl the protein core maintains its overall folding and structure, amide (15)N R(2) measurements and cross correlation rates between N-H dipole-dipole relaxation and (15)N CSA relaxation reveal a closer and stronger packing of the hydrophobic interactions in the protein as a response to the presence of denaturing agents in solution. A complete reorientation of the Met43 side chain toward the hydrophobic core is accomplished by the disappearance of the millisecond dynamics observed on the native form of Calbindin D(9k), while cross correlation rates provide evidence that the two-way hydrogen bond between Leu23 and Val61 is broken or substantially weakened. The substitution of the calcium ion in site II with the paramagnetic Ce(3+) ion allowed us to obtain a number of long-range nonconventional constraints, namely, pseudocontact shifts, which were used, together with the NOEs collected on the native state, to monitor subtle structural variations occurring in the non-native state of the protein. Although the average rmsd between the structures of native and non-native states is small (0.48 A), structural rearrangements could be reliably identified. Our results provide unprecedented information about the behavior of Calbindin D(9k) during the early steps of unfolding. Furthermore, they constitute strong evidence of the efficiency of paramagnetism-based constraints in monitoring subtle structural changes that are beyond the sensitivity of an approach based only on NOE.