Effects of Calcium Binding on the Side-chain Methyl Dynamics of Calbindin D9k: A 2H NMR Relaxation Study

Deuterium spin relaxation of fractionally deuterated ribonuclease H using paired 475 and 950 MHz NMR spectrometers

Deuterium (2H) spin relaxation of 13CH2D methyl groups has been widely applied to investigate picosecond-to-nanosecond conformational dynamics in proteins by solution-state NMR spectroscopy. The B0 dependence of the 2H spin relaxation rates is represented by a linear relationship between the spectral density function at three discrete frequencies J(0), J(ωD) and J(2ωD). In this study, the linear relation between 2H relaxation rates at B0 fields separated by a factor of two and the interpolation of rates at intermediate frequencies are combined for a more robust approach for spectral density mapping. The general usefulness of the approach is demonstrated on a fractionally deuterated (55%) and alternate 13C-12C labeled sample of E. coli RNase H. Deuterium relaxation rate constants (R1, R1ρ, RQ, RAP) were measured for 57 well-resolved 13CH2D moieties in RNase H at 1H frequencies of 475 MHz, 500 MHz, 900 MHz, and 950 MHz. The spectral density mapping of the 475/950 MHz data combination was performed independently and jointly to validate the expected relationship between data recorded at B0 fields separated by a factor of two. The final analysis was performed by jointly analyzing 475/950 MHz rates with 700 MHz rates interpolated from 500/900 MHz data to yield six J(ωD) values for each methyl peak. The J(ω) profile for each peak was fit to the original (τM, Sf2, τf) or extended model-free function (τM, Sf2, Ss2, τf, τs) to obtain optimized dynamic parameters.

Effect of monovalent ion binding on molecular dynamics of the S100-family calcium-binding protein calbindin D9k

Calbindin D_9k is a member of the S100 subfamily of EF-hand calcium binding proteins, and has served as an important model system for biophysical studies. The fast timescale dynamics of the calcium-free (apo) state is characterized using molecular dynamics simulations. Order parameters for the backbone NH bond vectors are determined from the simulations and compared with experimentally derived values, with a focus on the dynamics of calcium-binding site I. There is a significant discrepancy between simulated and experimental order parameters for site I residues in the case of no ion bound in site I. However, it was found in the simulations that a Na⁺ ion can bind in site I, and the resulting order parameters determined from the simulations are in excellent agreement with experiment. Comparisons are made to X-ray structures of other S100 family members in which Na⁺ ions were observed or suggested to be bound in site I. © 2019 Wiley Periodicals, Inc.

Enhanced spectral density mapping through combined multiple-field deuterium 13CH2D methyl spin relaxation NMR spectroscopy

Quadrupolar relaxation of ²H (D) nuclear spins is a powerful probe of conformational dynamics in biological macromolecules. Deuterium relaxation rate constants are determined by the spectral density function for reorientation of the C-D bond vector at zero, single-quantum, and double-quantum ²H frequencies. In the present work, ²H relaxation rate constants were measured for an E. coli ribonuclease H [U-²H, ¹⁵N] ILV-[¹³CH₂D] sample using 400, 500, 800, and 900 MHz NMR spectrometers and analyzed by three approaches to determine spectral density values. First, data recorded at each static magnetic field were analyzed independently. Second, data recorded at 400 and 800 MHz were analyzed jointly and data recorded at other fields were analyzed independently. Third, data recorded at 400 and 500 MHz were interpolated to 450 MHz, and the resulting two pairs of data, corresponding to 400 MHz/800 MHz and 450 MHz/900 MHz, were analyzed jointly. The second and third approaches rely on the identity between the double quantum frequency at the lower field and the single quantum frequency at the higher field. Spectral density values for 32 of the 48 resolvable ILV methyl resonances were fit by the Lipari-Szabo model-free formalism and used to validate the three methods. The three spectral density mapping methods performed equally well in cross validation with data recorded at 700 MHz. However, the third method yielded approximately 10-15% more precise estimates of model-free parameters and consequently provides a general strategy for analysis of ²H spin relaxation data in biological macromolecules.

Post-translational S-nitrosylation is an endogenous factor fine tuning the properties of human S100A1 protein

BACKGROUND

S100A1 protein is a proposed target of molecule-guided therapy for heart failure.

RESULTS

S-Nitrosylation of S100A1 is present in cells, increases Ca(2+) binding, and tunes the overall protein conformation.

CONCLUSION

Thiol-aromatic molecular switch is responsible for NO-related modification of S100A1 properties.

SIGNIFICANCE

Post-translational S-nitrosylation may provide functional diversity and specificity to S100A1 and other S100 protein family members. S100A1 is a member of the Ca(2+)-binding S100 protein family. It is expressed in brain and heart tissue, where it plays a crucial role as a modulator of Ca(2+) homeostasis, energy metabolism, neurotransmitter release, and contractile performance. Biological effects of S100A1 have been attributed to its direct interaction with a variety of target proteins. The (patho)physiological relevance of S100A1 makes it an important molecular target for future therapeutic intervention. S-Nitrosylation is a post-translational modification of proteins, which plays a role in cellular signal transduction under physiological and pathological conditions. In this study, we confirmed that S100A1 protein is endogenously modified by Cys(85) S-nitrosylation in PC12 cells, which are a well established model system for studying S100A1 function. We used isothermal calorimetry to show that S-nitrosylation facilitates the formation of Ca(2+)-loaded S100A1 at physiological ionic strength conditions. To establish the unique influence of the S-nitroso group, our study describes high resolution three-dimensional structures of human apo-S100A1 protein with the Cys(85) thiol group in reduced and S-nitrosylated states. Solution structures of the proteins are based on NMR data obtained at physiological ionic strength. Comparative analysis shows that S-nitrosylation fine tunes the overall architecture of S100A1 protein. Although the typical S100 protein intersubunit four-helix bundle is conserved upon S-nitrosylation, the conformation of S100A1 protein is reorganized at the sites most important for target recognition (i.e. the C-terminal helix and the linker connecting two EF-hand domains). In summary, this study discloses cysteine S-nitrosylation as a new factor responsible for increasing functional diversity of S100A1 and helps explain the role of S100A1 as a Ca(2+) signal transmitter sensitive to NO/redox equilibrium within cells.

Probing local structural fluctuations in myoglobin by size-dependent thiol-disulfide exchange

All proteins undergo local structural fluctuations (LSFs) or breathing motions. These motions are likely to be important for function but are poorly understood. LSFs were initially defined by amide hydrogen exchange (HX) experiments as opening events, which expose a small number of backbone amides to (1)H/(2)H exchange, but whose exchange rates are independent of denaturant concentration. Here, we use size-dependent thiol-disulfide exchange (SX) to characterize LSFs in single cysteine-containing variants of myoglobin (Mb). SX complements HX by providing information on motions that disrupt side chain packing interactions. Most importantly, probe reagents of different sizes and chemical properties can be used to characterize the size of structural opening events and the properties of the open state. We use thiosulfonate reagents (126-274 Da) to survey access to Cys residues, which are buried at specific helical packing interfaces in Mb. In each case, the free energy of opening increases linearly with the radius of gyration of the probe reagent. The slope and the intercept are interpreted to yield information on the size of the opening events that expose the buried thiol groups. The slope parameter varies by over 10-fold among Cys positions tested, suggesting that the sizes of breathing motions vary substantially throughout the protein. Our results provide insight to the longstanding question: how rigid or flexible are proteins in their native states?

Response of small-scale, methyl rotors to protein-ligand association: a simulation analysis of calmodulin-peptide binding

Changes in the free energy barrier (DeltaE), entropy, and motional parameters associated with the rotation of methyl groups in a protein (calmodulin (CaM)) on binding a ligand (the calmodulin-binding domain of smooth-muscle myosin (smMLCKp)) are investigated using molecular dynamics simulation. In both the bound and uncomplexed forms of CaM, the methyl rotational free energy barriers follow skewed-Gaussian distributions that are not altered significantly upon ligand binding. However, site-specific perturbations are found. Around 11% of the methyl groups in CaM exhibit changes in DeltaE greater than 0.7 kcal/mol on binding. The rotational entropies of the methyl groups exhibit a nonlinear dependence on DeltaE. The relations are examined between motional parameters (the methyl rotational NMR order parameter and the relaxation time) and DeltaE. Low-barrier methyl group rotational order parameters deviate from ideal tetrahedrality by up to approximately 20%. There is a correlation between rotational barrier changes and proximity to the protein-peptide binding interface. Methyl groups that exhibit large changes in DeltaE are found to report on elements in the protein undergoing structural change on binding.

Detection of nanosecond time scale side-chain jumps in a protein dissolved in water/glycerol solvent

In solution, the correlation time of the overall protein tumbling, tau(R), plays a role of a natural dynamics cutoff-internal motions with correlation times on the order of tau ( R ) or longer cannot be reliably identified on the basis of spin relaxation data. It has been proposed some time ago that the 'observation window' of solution experiments can be expanded by changing the viscosity of solvent to raise the value of tau(R). To further explore this concept, we prepared a series of samples of alpha-spectrin SH3 domain in solvent with increasing concentration of glycerol. In addition to the conventional (15)N labeling, the protein was labeled in the Val, Leu methyl positions ((13)CHD(2) on a deuterated background). The collected relaxation data were used in asymmetric fashion: backbone (15)N relaxation rates were used to determine tau(R) across the series of samples, while methyl (13)C data were used to probe local dynamics (side-chain motions). In interpreting the results, it has been initially suggested that addition of glycerol leads only to increases in tau(R), whereas local motional parameters remain unchanged. Thus the data from multiple samples can be analyzed jointly, with tau(R) playing the role of experimentally controlled variable. Based on this concept, the extended model-free model was constructed with the intent to capture the effect of ns time-scale rotameric jumps in valine and leucine side chains. Using this model, we made a positive identification of nanosecond dynamics in Val-23 where ns motions were already observed earlier. In several other cases, however, only tentative identification was possible. The lack of definitive results was due to the approximate character of the model-contrary to what has been assumed, addition of glycerol led to a gradual 'stiffening' of the protein. This and other observations also shed light on the interaction of the protein with glycerol, which is one of the naturally occurring osmoprotectants. In particular, it has been found that the overall protein tumbling is controlled by the bulk solvent, and not by a thin solvation layer which contains a higher proportion of water.

Frameworks for understanding long-range intra-protein communication

The phenomenon of intra-protein communication is fundamental to such processes as allostery and signaling, yet comparatively little is understood about its physical origins despite notable progress in recent years. This review introduces contemporary but distinct frameworks for understanding intra-protein communication by presenting both the ideas behind them and a discussion of their successes and shortcomings. The first framework holds that intra-protein communication is accomplished by the sequential mechanical linkage of residues spanning a gap between distal sites. According to the second framework, proteins are best viewed as ensembles of distinct structural microstates, the dynamical and thermodynamic properties of which contribute to the experimentally observable macroscale properties. Nuclear magnetic resonance (NMR) spectroscopy is a powerful method for studying intra-protein communication, and the insights into both frameworks it provides are presented through a discussion of numerous examples from the literature. Distinct from mechanical and thermodynamic considerations of intra-protein communication are recently applied graph and network theoretic analyses. These computational methods reduce complex three dimensional protein architectures to simple maps comprised of nodes (residues) connected by edges (inter-residue "interactions"). Analysis of these graphs yields a characterization of the protein's topology and network characteristics. These methods have shown proteins to be "small world" networks with moderately high local residue connectivities existing concurrently with a small but significant number of long range connectivities. However, experimental studies of the tantalizing idea that these putative long range interaction pathways facilitate one or several macroscopic protein characteristics are unfortunately lacking at present. This review concludes by comparing and contrasting the presented frameworks and methodologies for studying intra-protein communication and suggests a manner in which they can be brought to bear simultaneously to further enhance our understanding of this important fundamental phenomenon.

Hydrophobic core mutations in CI2 globally perturb fast side-chain dynamics similarly without regard to position

Protein dynamics is currently an area of intense research because of its importance as complementary information to the huge quantity of available data relating protein structure and function. Because it is usually the influence of dynamics on function that is studied, the physical determinants of the distribution of flexibility in proteins have not been explored as thoroughly. In the present NMR study, an expanded suite of five (2)H relaxation experiments was used to characterize the picosecond-to-nanosecond side-chain dynamics of chymotrypsin inhibitor 2 (CI2) and five hydrophobic core mutants, some of which are members of the folding nucleus. Because CI2 is a homologue of the serine protease inhibitor eglin c, which has already been extensively characterized in terms of its dynamics, it was possible to compare not only side-chain dynamics but also the responses of these dynamics to analogous mutations. Remarkably, each of the five core mutations in CI2 led to similar and reproducible increases in side-chain flexibility throughout the entire structure. Although the expanded suite of (2)H relaxation experiments does not affect model selection for the vast majority of residues, it did enable the detection of increasing levels of nanosecond-scale motions in CI2's reactive site binding loop as the L68 side chain was progressively shortened by mutation. Collectively, we observed that the CI2 mutants are more dynamically similar to each other than to the more rigid wild-type CI2, from which we propose that wild-type CI2 has been optimized to a specific level of rigidity which may aid in its function as a serine protease inhibitor. We also observed that the pattern of side-chain dynamics of CI2 is quantitatively similar to eglin c, but that this similarity is lost upon mutating both proteins at an equivalent position. Finally, (15)N relaxation was used to characterize the backbone dynamics of wild-type and mutant CI2. Interestingly, mutation at folding nucleus positions led to widespread increases in backbone flexibility, whereas non-folding-nucleus positions led to increases in flexibility in the C-terminal half of the protein only.

Determination of NMR lineshape anisotropy of guest molecules within inclusion complexes from molecular dynamics simulations

Nonspherical cages in inclusion compounds can result in non-uniform motion of guest species in these cages and anisotropic lineshapes in NMR spectra of the guest. Herein, we develop a methodology to calculate lineshape anisotropy of guest species in cages based on molecular dynamics simulations of the inclusion compound. The methodology is valid for guest atoms with spin 1/2 nuclei and does not depend on the temperature and type of inclusion compound or guest species studied. As an example, the nonspherical shape of the structure I (sI) clathrate hydrate large cages leads to preferential alignment of linear CO(2) molecules in directions parallel to the two hexagonal faces of the cages. The angular distribution of the CO(2) guests in terms of a polar angle theta and azimuth angle phi and small amplitude vibrational motions in the large cage are characterized by molecular dynamics simulations at different temperatures in the stability range of the CO(2) sI clathrate. The experimental (13)C NMR lineshapes of CO(2) guests in the large cages show a reversal of the skew between the low temperature (77 K) and the high temperature (238 K) limits of the stability of the clathrate. We determine the angular distributions of the guests in the cages by classical MD simulations of the sI clathrate and calculate the (13)C NMR lineshapes over a range of temperatures. Good agreement between experimental lineshapes and calculated lineshapes is obtained. No assumptions regarding the nature of the guest motions in the cages are required.

An improved picture of methyl dynamics in proteins from slowly relaxing local structure analysis of 2H spin relaxation

Protein dynamics is intimately related to biological function. Core dynamics is usually studied with 2H spin relaxation of the 13CDH2 group, analyzed traditionally with the model-free (MF) approach. We showed recently that MF is oversimplified in several respects. This includes the assumption that the local motion of the dynamic probe and the global motion of the protein are decoupled, the local geometry is simple, and the local ordering is axially symmetric. Because of these simplifications MF has yielded a puzzling picture where the methyl rotation axis is moving rapidly with amplitudes ranging from nearly complete disorder to nearly complete order in tightly packed protein cores. Our conclusions emerged from applying to methyl dynamics in proteins the slowly relaxing local structure (SRLS) approach of Polimeno and Freed (Polimeno, A.; Freed, J. H. J. Phys. Chem. 1995, 99, 10995-11006.), which can be considered the generalization of MF, with all the simplifications mentioned above removed. The SRLS picture derived here for the B1 immunoglobulin binding domain of peptostreptococcal protein L, studied over the temperature range of 15-45 degrees C, is fundamentally different from the MF picture. Thus, methyl dynamics is characterized structurally by rhombic local potentials with varying symmetries and dynamically by tenfold slower rates of local motion. On average, potential rhombicity decreases, mode-coupling increases, and the rate of local motion increases with increasing temperature. The average activation energy for local motion is 2.0 +/- 0.2 kcal/mol. Mode-coupling affects the analysis even at 15 degrees C. The accuracy of the results is improved by including in the experimental data set relaxation rates associated with rank 2 coherences.

Substrate recognition reduces side-chain flexibility for conserved hydrophobic residues in human Pin1

Pin1 is a peptidyl-prolyl isomerase consisting of a WW domain and a catalytic isomerase (PPIase) domain connected by a flexible linker. Pin1 recognizes phospho-Ser/Thr-Pro motifs in cell-signaling proteins, and is both a cancer and an Alzheimer's disease target. Here, we provide novel insight into the functional motions underlying Pin1 substrate interaction using nuclear magnetic resonance deuterium ((2)D) and carbon ((13)C) spin relaxation. Specifically, we compare Pin1 side-chain motions in the presence and absence of a known phosphopeptide substrate derived from the mitotic phosphatase Cdc25. Substrate interaction alters Pin1 side-chain motions on both the microsecond-millisecond (mus-ms) and picosecond-nanosecond (ps-ns) timescales. Alterations include loss of ps-ns flexibility along an internal conduit of hydrophobic residues connecting the catalytic site with the interdomain interface. These residues are conserved among Pin1 homologs; hence, their dynamics are likely important for the Pin1 mechanism.