15N R2 Research Articles

Cross-Polarization of Insensitive Nuclei from Water Protons for Detection of Protein-Ligand Binding.

Hyperpolarization derived from water protons enhances the NMR signal of 15N nuclei in a small molecule, enabling the sensitive detection of a protein-ligand interaction. The water hyperpolarized by dissolution dynamic nuclear polarization (D-DNP) acts as a universal signal enhancement agent. The 15N signal of benzamidine was increased by 1480-fold through continuous polarization transfer by J-coupling-mediated cross-polarization (J-CP) via the exchangeable protons. The signal enhancement factor favorably compares to factors of 110- or 17-fold using non-CP-based polarization transfer mechanisms. The hyperpolarization enabled detection of the binding of benzamidine to the target protein trypsin with a single-scan measurement of 15N R2 relaxation. J-CP provides an efficient polarization mechanism for 15N or other low-frequency nuclei near an exchangeable proton. The hyperpolarization transfer sustained within the relaxation time limit of water protons additionally can be applied for the study of macromolecular structure and biological processes.

Dilated Cardiomyopathy Mutations and Phosphorylation disrupt the Active Orientation of Cardiac Troponin C

Cardiac troponin (cTn) is made up of three subunits, cTnC, cTnI, and cTnT. The regulatory N-terminal domain of cTnC (cNTnC) controls cardiac muscle contraction in a calcium-dependent manner. We show that calcium-saturated cNTnC can adopt two different orientations, with the “active” orientation consistent with the 2020 cryo-EM structure of the activated cardiac thin filament by Yamada et al. Using solution NMR 15N R2 relaxation analysis, we demonstrate that the two domains of cTnC tumble independently (average R2 10 s−1), being connected by a flexible linker. However, upon addition of cTnI1-77, the complex tumbles as a rigid unit (R2 30 s−1). cTnI phosphomimetic mutants S22D/S23D, S41D/S43D and dilated cardiomyopathy- (DCM-)associated mutations cTnI K35Q, cTnC D75Y, and cTnC G159D destabilize the active orientation of cNTnC, with intermediate 15N R2 rates (R2 17–23 s−1). The active orientation of cNTnC is stabilized by the flexible tails of cTnI, cTnI1-37 and cTnI135-209. Surprisingly, when cTnC is incorporated into complexes lacking these tails (cTnC-cTnI38-134, cTnC-cTnT223-288, or cTnC-cTnI38-134-cTnT223-288), the cNTnC domain is still immobilized, revealing a new interaction between cNTnC and the IT-arm that stabilizes a “dormant” orientation. We propose that the calcium sensitivity of the cardiac troponin complex is regulated by an equilibrium between active and dormant orientations, which can be shifted through post-translational modifications or DCM-associated mutations.

CALX-CBD1 Ca2+-Binding Cooperativity Studied by NMR Spectroscopy and ITC with Bayesian Statistics

The Na+/Ca2+ exchanger of Drosophila melanogaster, CALX, is the main Ca2+-extrusion mechanism in olfactory sensory neurons and photoreceptor cells. Na+/Ca2+ exchangers have two Ca2+ sensor domains, CBD1 and CBD2. In contrast to the mammalian homologs, CALX is inhibited by Ca2+ binding to CALX-CBD1, whereas CALX-CBD2 does not bind Ca2+ at physiological concentrations. CALX-CBD1 consists of a β-sandwich and displays four Ca2+-binding sites at the tip of the domain. In this study, we used NMR spectroscopy and isothermal titration calorimetry (ITC) to investigate the cooperativity of Ca2+ binding to CALX-CBD1. We observed that this domain binds Ca2+ in the slow exchange regime at the NMR chemical shift timescale. Ca2+ binding restricts the dynamics in the Ca2+-binding region. Experiments of 15N chemical exchange saturation transfer and 15N R2 dispersion allowed the determination of Ca2+ dissociation rates (∼30 s−1). NMR titration curves of residues in the Ca2+-binding region were sigmoidal because of the contribution of chemical exchange to transverse magnetization relaxation rates, R2. Hence, a novel, to our knowledge, approach to analyze NMR titration curves was proposed. Ca2+-binding cooperativity was examined assuming two different stoichiometric binding models and using a Bayesian approach for data analysis. Fittings of NMR and ITC binding curves to the Hill model yielded nHill ∼2.9, near maximal cooperativity (nHill = 4). By assuming a stepwise model to interpret the ITC data, we found that the probability of binding from 2 up to 4 Ca2+ is approximately three orders of magnitude higher than that of binding a single Ca2+. Hence, four Ca2+ ions bind almost simultaneously to CALX-CBD1. Cooperative Ca2+ binding is key to enable this exchanger to efficiently respond to changes in the intracellular Ca2+ concentration in sensory neuronal cells.

The Solution Structure and Dynamics of Cd-Metallothionein from Helix pomatia Reveal Optimization for Binding Cd over Zn.

Metallothioneins (MTs) are cysteine-rich polypeptides that are naturally found coordinated to monovalent and/or divalent transition metal ions. Three metallothionein isoforms from the Roman snail Helix pomatia are known. They differ in their physiological metal load and in their specificity for transition metal ions such as Cd2+ (HpCdMT isoform) and Cu+ (HpCuMT isoform) or in the absence of a defined metal specificity (HpCd/CuMT isoform). We have determined the solution structure of the Cd-specific isoform (HpCdMT) by nuclear magnetic resonance spectroscopy using recombinant isotopically labeled protein loaded with Zn2+ or Cd2+. Both structures display two-domain architectures, where each domain comprises a characteristic three-metal cluster similar to that observed in the β-domains of vertebrate MTs. The polypeptide backbone is well-structured over the entire sequence, including the interdomain linker. Interestingly, the two domains display mutual contacts, as observed before for the metallothionein of the snail Littorina littorea, to which both N- and C-terminal domains are highly similar. Increasing the length of the linker motionally decouples both domains and removes mutual contacts between them without having a strong effect on the stability of the individual domains. The structures of Cd6- and Zn6-HpCdMT are nearly identical. However, 15N relaxation, in particular 15N R2 rates, is accelerated for many residues of Zn6-HpCdMT but not for Cd6-HpCdMT, revealing the presence of conformational exchange effects. We suggest that this snail MT isoform is evolutionarily optimized for binding Cd rather than Zn.

Millisecond Timescale Motions Connect Amino Acid Interaction Networks in Alpha Tryptophan Synthase.

Tryptophan synthase is a model system for understanding allosteric regulation within enzyme complexes. Amino acid interaction networks were previously delineated in the isolated alpha subunit (αTS) in the absence of the beta subunit (βTS). The amino acid interaction networks were different between the ligand-free enzyme and the enzyme actively catalyzing turnover. Previous X-ray crystallography studies indicated only minor localized changes when ligands bind αTS, and so, structural changes alone could not explain the changes to the amino acid interaction networks. We hypothesized that the network changes could instead be related to changes in conformational dynamics. As such, we conducted nuclear magnetic resonance relaxation studies on different substrate- and products-bound complexes of αTS. Specifically, we collected 15N R2 relaxation dispersion data that reports on microsecond-to-millisecond timescale motion of backbone amide groups. These experiments indicated that there are conformational exchange events throughout αTS. Substrate and product binding change specific motional pathways throughout the enzyme, and these pathways connect the previously identified network residues. These pathways reach the αTS/βTS binding interface, suggesting that the identified dynamic networks may also be important for communication with the βTS subunit.

15N transverse relaxation measurements for the characterization of µs–ms dynamics are deteriorated by the deuterium isotope effect on 15N resulting from solvent exchange

15N R2 relaxation measurements are key for the elucidation of the dynamics of both folded and intrinsically disordered proteins (IDPs). Here we show, on the example of the intrinsically disordered protein α-synuclein and the folded domain PDZ2, that at physiological pH and near physiological temperatures amide-water exchange can severely skew Hahn-echo based 15N R2 relaxation measurements as well as low frequency data points in CPMG relaxation dispersion experiments. The nature thereof is the solvent exchange with deuterium in the sample buffer, which modulates the 15N chemical shift tensor via the deuterium isotope effect, adding to the apparent relaxation decay which leads to systematic errors in the relaxation data. This results in an artificial increase of the measured apparent 15N R2 rate constants-which should not be mistaken with protein inherent chemical exchange contributions, Rex, to 15N R2. For measurements of 15N R2 rate constants of IDPs and folded proteins at physiological temperatures and pH, we recommend therefore the use of a very low D2O molar fraction in the sample buffer, as low as 1%, or the use of an external D2O reference along with a modified 15N R2 Hahn-echo based experiment. This combination allows for the measurement of Rex contributions to 15N R2 originating from conformational exchange in a time window from µs to ms.

P15PAF Is an Intrinsically Disordered Protein with Nonrandom Structural Preferences at Sites of Interaction with Other Proteins

We present to our knowledge the first structural characterization of the proliferating-cell-nuclear-antigen-associated factor p15PAF, showing that it is monomeric and intrinsically disordered in solution but has nonrandom conformational preferences at sites of protein-protein interactions. p15PAF is a 12 kDa nuclear protein that acts as a regulator of DNA repair during DNA replication. The p15PAF gene is overexpressed in several types of human cancer. The nearly complete NMR backbone assignment of p15PAF allowed us to measure 86 N-HN residual dipolar couplings. Our residual dipolar coupling analysis reveals nonrandom conformational preferences in distinct regions, including the proliferating-cell-nuclear-antigen-interacting protein motif (PIP-box) and the KEN-box (recognized by the ubiquitin ligase that targets p15PAF for degradation). In accordance with these findings, analysis of the 15N R2 relaxation rates shows a relatively reduced mobility for the residues in these regions. The agreement between the experimental small angle x-ray scattering curve of p15PAF and that computed from a statistical coil ensemble corrected for the presence of local secondary structural elements further validates our structural model for p15PAF. The coincidence of these transiently structured regions with protein-protein interaction and posttranslational modification sites suggests a possible role for these structures as molecular recognition elements for p15PAF.

Proton-decoupled CPMG: A better experiment for measuring 15N R2 relaxation in disordered proteins

15N R2 relaxation is one of the most informative experiments for characterization of intrinsically disordered proteins (IDPs). Small changes in nitrogen R2 rates are often used to determine how IDPs respond to various biologically relevant perturbations such as point mutations, posttranslational modifications and weak ligand interactions. However collecting high-quality 15N relaxation data can be difficult. Of necessity, the samples of IDPs are often prepared with low protein concentration and the measurement time can be limited because of rapid sample degradation. Furthermore, due to hardware limitations standard experiments such as 15N spin-lock and CPMG can sample the relaxation decay only to ca. 150ms. This is much shorter than 15N T2 times in disordered proteins at or near physiological temperature. As a result, the sampling of relaxation decay profiles in these experiments is suboptimal, which further lowers the precision of the measurements.Here we report a new implementation of the proton-decoupled (PD) CPMG experiment which allows one to sample 15N R2 relaxation decay up to ca. 0.5–1s. The new experiment has been validated through comparison with the well-established spin-lock measurement. Using dilute samples of denatured ubiquitin, we have demonstrated that PD-CPMG produces up to 3-fold improvement in the precision of the data. It is expected that for intrinsically disordered proteins the gains may be even more substantial. We have also shown that this sequence has a number of favorable properties: (i) the spectra are recorded with narrow linewidth in nitrogen dimension; (ii) 15N offset correction is small and easy to calculate; (iii) the experiment is immune to various spurious effects arising from solvent exchange; (iv) the results are stable with respect to pulse miscalibration and rf field inhomogeneity; (v) with minimal change, the pulse sequence can also be used to measure R2 relaxation of 15Nε spins in arginine side chains. We anticipate that the new experiment will be a valuable addition to the NMR toolbox for studies of IDPs.

Crystallographic and NMR Evidence for Flexibility in Oligosaccharyltransferases and Its Catalytic Significance

Oligosaccharyltransferase (OST) is a membrane-bound enzyme that catalyzes the transfer of an oligosaccharide to an asparagine residue in glycoproteins. It possesses a binding pocket that recognizes Ser and Thr residues at the+2 position in theN-glycosylation consensus, Asn-X-Ser/Thr. We determined the crystal structures of the C-terminal globular domains of the catalytic subunits of two archaeal OSTs. A comparison with previously determined structures identified a segment with remarkable conformational plasticity, induced by crystal contact effects. We characterized its dynamic properties in solution by (15)N NMR relaxation analyses. Intriguingly, the mobile region contains the+2 Ser/Thr-binding pocket. In agreement, the flexibility restriction forced by an engineered disulfide crosslink abolished the enzymatic activity, and its cleavage fully restored activity. These results suggest the necessity of multiple conformational states in the reaction. The dynamic nature of the Ser/Thr pocket could facilitate the efficient scanning of N-glycosylation sequons along nascent polypeptide chains.

Structure and Dynamics of the First Archaeal Parvulin Reveal a New Functionally Important Loop in Parvulin-type Prolyl Isomerases

Parvulins are a group of peptidyl-prolyl isomerases (PPIases) responsible for important biological processes in all kingdoms of life. The PinA protein from the psychrophilic archaeon Cenarchaeum symbiosum is a parvulin-like PPIase. Due to its striking similarity to the human parvulins Pin1 and Par14, PinA constitutes an interesting subject for structural and functional studies. Here, we present the first high resolution NMR structure of an archaeal parvulin, PinA, based on 1798 conformational restraints. Structure calculation yields an ensemble of 20 convergent low energy structures with a backbone r.m.s.d. value of 0.6 Å within the secondary structure elements. The overall fold of PinA comprises the β-α3-β-α-β2 fold typical for all parvulin structures known so far, but with helix III being a short 310-helix. A detailed comparison of this high resolution structure of the first archaeal PinA protein with bacterial and eukaryotic parvulin PPIase structures reveals an atypically large catalytic binding site. This feature provides an explanation for cold-adapted protein function. Moreover, the residues in and around 310-helix III exhibit strong intramolecular dynamics on a microsecond to millisecond timescale and display structural heterogeneity within the NMR ensemble. A putative peptide ligand was found for PinA by phage display and was used for 1H-15N-HSQC titrations. Again, the flexible region around 310-helix III as well as residues of the peptide binding pocket showed the strongest chemical shift perturbations upon peptide binding. The local flexibility of this region also was modulated by ligand binding. A glycine and two positively charged residues are conserved in most parvulin proteins in this flexible loop region, which may be of general functional importance for parvulin-type PPIases.

Native like structure in the unfolded state of the villin headpiece helical subdomain, an ultrafast folding protein

The villin headpiece subdomain, HP36, is the smallest naturally occurring protein that folds cooperatively. Its small size, rapid folding, and simple three-helix topology have made it an extremely popular system for computational studies of protein folding. The role of unfolded state structure in rapid folding is an area of active investigation, but relatively little is known about the properties of unfolded states under native conditions. A peptide fragment, HP21, which contains the first and second helices of HP36 has been shown to be a good model for structure in the unfolded state of the intact domain but a detailed description of the conformational propensities of HP21 is lacking and the balance between native and nonnative interactions is not known. A series of three-dimensional NMR experiments were performed on (13)C, (15)N-labeled HP21 to investigate in detail its conformational propensities. Analysis of (13)C(alpha), (13)C(beta), (13)CO chemical shifts, Deltadelta(13)C(alpha) - Deltadelta(13)C(beta) secondary shifts, the secondary structure propensity scores, NOEs, (15)N R(2) values and comparison of experimental chemical shifts with those of HP36 and with chemical shifts calculated using the SHIFTS and SHIFTX programs all indicate that there is significant native like structure in the HP21 ensemble, and thus by implication in the unfolded state of HP36.

The Low-pH Unfolded State of the C-Terminal Domain of the Ribosomal Protein L9 Contains Significant Secondary Structure in the Absence of Denaturant but Is No More Compact Than the Low-pH Urea Unfolded State

There is considerable interest in the properties of the unfolded states of proteins, particularly unfolded states which can be populated in the absence of high concentrations of denaturants. Interest in the unfolded state ensemble reflects the fact that it is the starting point for protein folding as well as the reference state for protein stability studies and can be the starting state for pathological aggregation. The unfolded state of the C-terminal domain (residues 58-149) of the ribosomal protein L9 (CTL9) can be populated in the absence of denaturant at low pH. CTL9 is a 92-residue globular alpha, beta protein. The low-pH unfolded state contains more secondary structure than the low-pH urea unfolded state, but it is not a molten globule. Backbone ( (1)H, (13)C, and (15)N) NMR assignments as well as side chain (13)C beta and (1)H beta assignments and (15)N R 2 values were obtained for the pH 2.0 unfolded form of CTL9 and for the urea unfolded state at pH 2.5. Analysis of the deviations of the chemical shifts from random coil values indicates that residues that comprise the two helices in the native state show a clear preference for adopting helical phi and psi angles in the pH 2.0 unfolded state. There is a less pronounced but nevertheless clear tendency for residues 107-124 to preferentially populate helical phi and psi values in the unfolded state. The urea unfolded state has no detectable tendency to populate any type of secondary structure even though it is as compact as the pH 2.0 unfolded state. Comparison of the two unfolded forms of CTL9 provides direct experimental evidence that states which differ significantly in their secondary structure can have identical hydrodynamic properties. This in turn demonstrates that global parameters such as R h or R g are very poor indicators of "random coil" behavior.

Entropy-driven cAMP-dependent Allosteric Control of Inhibitory Interactions in Exchange Proteins Directly Activated by cAMP

Exchange proteins directly activated by cAMP (EPACs) are guanine nucleotide-exchange factors for the small GTPases Rap1 and Rap2 and represent a key receptor for the ubiquitous cAMP second messenger in eukaryotes. The cAMP-dependent activation of apoEPAC is typically rationalized in terms of a preexisting equilibrium between inactive and active states. Structural and mutagenesis analyses have shown that one of the critical determinants of the EPAC activation equilibrium is a cluster of salt bridges formed between the catalytic core and helices alpha1 and alpha2 at the N terminus of the cAMP binding domain and commonly referred to as ionic latch (IL). The IL stabilizes the inactive states in a closed topology in which access to the catalytic domain is sterically occluded by the regulatory moiety. However, it is currently not fully understood how the IL is allosterically controlled by cAMP. Chemical shift mapping studies consistently indicate that cAMP does not significantly perturb the structure of the IL spanning sites within the regulatory region, pointing to cAMP-dependent dynamic modulations as a key allosteric carrier of the cAMP-signal to the IL sites. Here, we have therefore investigated the dynamic profiles of the EPAC1 cAMP binding domain in its apo, cAMP-bound, and Rp-cAMPS phosphorothioate antagonist-bound forms using several 15N relaxation experiments. Based on the comparative analysis of dynamics in these three states, we have proposed a model of EPAC activation that incorporates the dynamic features allosterically modulated by cAMP and shows that cAMP binding weakens the IL by increasing its entropic penalty due to dynamic enhancements.

Tailoring Relaxation Dispersion Experiments for Fast-Associating Protein Complexes

NMR relaxation dispersion spectroscopy is a powerful technique to elucidate the mechanism of protein−protein binding reactions. However, it is difficult to optimize the concentration ratios that give relaxation dispersions of appropriate amplitude to determine accurate kinetic and thermodynamic parameters, especially in cases of very tight binding. In this study, we have obtained 15N R2 dispersions of Asn803-hydroxylated hypoxia-inducible factor-1α (HIF−OH) in the presence of a substoichiometric amount of its target protein, the transcriptional adapter zinc-binding (TAZ1) domain of CREB binding protein, whereas no R2 dispersion was observed for the bound state of HIF−OH at 1:1 concentration ratio because the binding is too tight. Although the R2 dispersions were measured for the free peaks of HIF−OH, they enabled us to quantitate the kinetic and thermodynamic parameters of HIF−OH/TAZ1 binding process. Simulations of effective R2 rates revealed that the association rate is the key factor to determine the amplitude of R2 dispersions. By careful optimization of the concentration ratio, the R2 dispersion method should be generally applicable for studying a wide range of protein−protein, protein−nucleic acid, and protein−small molecule interactions.

Error estimation and global fitting in transverse-relaxation dispersion experiments to determine chemical-exchange parameters

Off-resonance effects can introduce significant systematic errors in R2 measurements in constant-time Carr-Purcell-Meiboom-Gill (CPMG) transverse relaxation dispersion experiments. For an off-resonance chemical shift of 500 Hz, 15N relaxation dispersion profiles obtained from experiment and computer simulation indicated a systematic error of ca. 3%. This error is three- to five-fold larger than the random error in R2 caused by noise. Good estimates of total R2 uncertainty are critical in order to obtain accurate estimates in optimized chemical exchange parameters and their uncertainties derived from chi2 minimization of a target function. Here, we present a simple empirical approach that provides a good estimate of the total error (systematic + random) in 15N R2 values measured for the HIV protease. The advantage of this empirical error estimate is that it is applicable even when some of the factors that contribute to the off-resonance error are not known. These errors are incorporated into a chi2 minimization protocol, in which the Carver-Richards equation is used fit the observed R2 dispersion profiles, that yields optimized chemical exchange parameters and their confidence limits. Optimized parameters are also derived, using the same protein sample and data-fitting protocol, from 1H R2 measurements in which systematic errors are negligible. Although 1H and 15N relaxation profiles of individual residues were well fit, the optimized exchange parameters had large uncertainties (confidence limits). In contrast, when a single pair of exchange parameters (the exchange lifetime, tau(ex), and the fractional population, p(a)), were constrained to globally fit all R2 profiles for residues in the dimer interface of the protein, confidence limits were less than 8% for all optimized exchange parameters. In addition, F-tests showed that quality of the fits obtained using tau(ex), p(a) as global parameters were not improved when these parameters were free to fit the R2 profiles of individual residues. Finally, nearly the same optimized global tau(ex), p(a) values were obtained, when the 1H and 15N data sets for residues in the dimer interface, were fit independently; the difference in optimized global parameters, ca. 10%, was of marginal significance according to the F-test.

Microsecond time scale dynamics in the RXR DNA-binding domain from a combination of spin-echo and off-resonance rotating frame relaxation measurements.

Slow protein dynamics can be studied by 15N spin-echo (CPMG) and off-resonance rotating frame relaxation through the effective field dependence of the exchange-mediated relaxation contribution. It is shown that, by a combination of these complementary techniques, a more extended sampling of the microsecond time scale processes is achieved than by either method alone. 15N R2 and improved off-resonance R1 rho experiments [Mulder et al. (1998) J. Magn. Reson., 131, 351-357] were applied to the 9-cis-retinoic acid receptor DNA-binding domain and allowed the identification of 14 residues exhibiting microsecond time scale dynamics. Assuming exchange between two conformational substates, average lifetimes ranging from 37 to 416 microseconds, and chemical shift differences of up to 3 ppm were obtained. The largest perturbation of tertiary structure was observed for the second zinc finger region, which was found to be disordered in the solution structure [Holmbeck et al. (1998) J. Mol. Biol., 281, 271-284]. Since this zinc-coordinating domain comprises the principal dimerization interface for RXR in a wide repertoire of complexes with different hormone receptors to their cognate response elements, this finding has important implications for our understanding of nuclear receptor assembly on DNA direct repeats. The flexibility observed for the dimerization domain may explain how RXR, through the ability to adaptively interact with a wide variety of highly homologous partner molecules, demonstrates such a versatile DNA-binding repertoire.