Many residues in cytochromec populate alternative states under equilibrium conditions

Characterizing the Conformational Dynamics of Human SUMO2: Insights into its Interaction with Metal Ions and SIMs

SUMO (Small Ubiquitin-like Modifiers) proteins are involved in a crucial post-translational modification commonly termed as SUMOylation. In this work, we have investigated the native-state conformational flexibility of human SUMO2 and its interaction with Cu2+ and Zn2+ ions using 15N-1H based 2D NMR spectroscopy. After SUMO1, SUMO2 is the most studied SUMO isoform in humans which shares 45 % and ~80 % similarity with SUMO1 in terms of sequence and structure, respectively. In this manuscript, we demonstrate that compared to SUMO1, several amino acids around the α1-helix region of SUMO2 access energetically similar near-native conformations. These conformations could play a crucial role in SUMO2's non-covalent interactions with SUMO interaction motifs (SIMs) on other proteins. The C-terminal of SUMO2 was found to bind strongly with Cu2+ ions resulting in a trimeric structure as observed by gel electrophoresis. This interaction seems to interfere in its non-covalent interaction with a V/I-x-V/I-V/I based SIM in Daxx protein.

Negative Thermal Expansion and Disorder-to-Order Collapse of an Intrinsically Disordered Protein under Marginally Denaturing Conditions

Recent work with intrinsically disordered proteins (IDPs) has projected a myriad of their survival instincts based mainly on the total charge content, the abundance of polar residues, and the paucity of hydrophobic amino acids. This work uses a plant IDP AtPP16-1 (Arabidopsis thaliana phloem protein class 16-1), whose solution NMR structure was determined by us recently, to show legitimate negative thermal expansion (NTE) of the native state. The thermal expansion continues to be negative even as the tertiary structure is perturbed by ultralow levels of urea up to 0.4 M. The NTE of these subdenatured states is called apparent NTE because they are prone to undergo conformational changes with temperature. Hydrodynamic shrinkage of the NTE IDP is also observed by dynamic light scattering (DLS) and NMR-measured global rotational correlation time (τ_c). The protein with denatured tertiary structure but otherwise preserved native-state secondary structure collapses to a dynamically rigid state. The data are mainly based on thermal coefficients of chemical shift and nuclear relaxation measured by heteronuclear NMR. The hydrodynamic shrinkage and collapse under marginally varying solvent compositions that may arise from unstable tertiary structure and dynamic disorder of chain segments across the backbone could be a generic property of IDPs.

A fine balance of hydrophobic-electrostatic communication pathways in a pH-switching protein

Allostery is the phenomenon of coupling between distal binding sites in a protein. Such coupling is at the crux of protein function and regulation in a myriad of scenarios, yet determining the molecular mechanisms of coupling networks in proteins remains a major challenge. Here, we report mechanisms governing pH-dependent myristoyl switching in monomeric hisactophilin, whereby the myristoyl moves between a sequestered state, i.e., buried within the core of the protein, to an accessible state, in which the myristoyl has increased accessibility for membrane binding. Measurements of the pH and temperature dependence of amide chemical shifts reveal protein local structural stability and conformational heterogeneity that accompany switching. An analysis of these measurements using a thermodynamic cycle framework shows that myristoyl-proton coupling at the single-residue level exists in a fine balance and extends throughout the protein. Strikingly, small changes in the stereochemistry or size of core and surface hydrophobic residues by point mutations readily break, restore, or tune myristoyl switch energetics. Synthesizing the experimental results with those of molecular dynamics simulations illuminates atomistic details of coupling throughout the protein, featuring a large network of hydrophobic interactions that work in concert with key electrostatic interactions. The simulations were critical for discerning which of the many ionizable residues in hisactophilin are important for switching and identifying the contributions of nonnative interactions in switching. The strategy of using temperature-dependent NMR presented here offers a powerful, widely applicable way to elucidate the molecular mechanisms of allostery in proteins at high resolution.

Structure, dynamics, and function of PsDef2 defensin from Pinus sylvestris

Plant defensins demonstrate high structural stability at extreme temperatures and pH values and, in general, are non-toxic to mammalian cells. These properties make them attractive candidates for use in biotechnology and biomedicine. Knowing the structure-function relationship is desirable to guide the design of plant defensin-based applications. Thus far, the broad range of biological activities was described only for one defensin from gymnosperms, the defensin PsDef1 from Scots pine. Here, we report that closely related defensin from the same taxonomy group, PsDef2, differing from PsDef1 by six amino acids, also possesses antimicrobial, antibacterial, and insect α-amylase inhibitory activities. We also report the solution structure and dynamics properties of PsDef2 assessed using a combination of experimental nuclear magnetic resonance (NMR) techniques. Lastly, we perform a comparative analysis of PsDef2 and PsDef1 gaining a molecular-level insight into their structure-dynamics-function relationship.

Destabilization of polar interactions in the prion protein triggers misfolding and oligomerization

The prion protein (PrP) misfolds and oligomerizes at pH 4 in the presence of physiological salt concentrations. Low pH and salt cause structural perturbations in the monomeric prion protein that lead to misfolding and oligomerization. However, the changes in stability within different regions of the PrP prior to oligomerization are poorly understood. In this study, we have characterized the local stability in PrP at high resolution using amide temperature coefficients (T_C ) measured by nuclear magnetic resonance (NMR) spectroscopy. The local stability of PrP was investigated under native as well as oligomerizing conditions. We have also studied the rapidly oligomerizing PrP variant (Q216R) and the protective PrP variant (A6). We report that at low pH, salt destabilizes PrP at several polar residues, and the hydrogen bonds in helices α2 and α3 are weakened. In addition, salt changes the curvature of the α3 helix, which likely disrupts α2-α3 contacts and leads to oligomerization. These results are corroborated by the T_C values of rapidly oligomerizing Q216R-PrP. The poly-alanine substitution in A6-PrP stabilizes α2, which prevents oligomerization. Altogether, these results highlight the importance of native polar interactions in determining the stability of PrP and reveal the structural disruptions in PrP that lead to misfolding and oligomerization.

The N-terminal Tails of Histones H2A and H2B Adopt Two Distinct Conformations in the Nucleosome with Contact and Reduced Contact to DNA

The nucleosome comprises two histone dimers of H2A-H2B and one histone tetramer of (H3-H4)₂, wrapped around by ~145 bp of DNA. Detailed core structures of nucleosomes have been established by X-ray and cryo-EM, however, histone tails have not been visualized. Here, we have examined the dynamic structures of the H2A and H2B tails in 145-bp and 193-bp nucleosomes using NMR, and have compared them with those of the H2A and H2B tail peptides unbound and bound to DNA. Whereas the H2A C-tail adopts a single but different conformation in both nucleosomes, the N-tails of H2A and H2B adopt two distinct conformations in each nucleosome. To clarify these conformations, we conducted molecular dynamics (MD) simulations, which suggest that the H2A N-tail can locate stably in either the major or minor grooves of nucleosomal DNA. While the H2B N-tail, which sticks out between two DNA gyres in the nucleosome, was considered to adopt two different orientations, one toward the entry/exit side and one on the opposite side. Then, the H2A N-tail minor groove conformation was obtained in the H2B opposite side and the H2B N-tail interacts with DNA similarly in both sides, though more varied conformations are obtained in the entry/exit side. Collectively, the NMR findings and MD simulations suggest that the minor groove conformer of the H2A N-tail is likely to contact DNA more strongly than the major groove conformer, and the H2A N-tail reduces contact with DNA in the major groove when the H2B N-tail is located in the entry/exit side.

Structural Dynamic Heterogeneity of Polyubiquitin Subunits Affects Phosphorylation Susceptibility

Polyubiquitin is a multifunctional protein tag formed by the covalent conjugation of ubiquitin molecules. Due to the high rigidity of the ubiquitin fold, the ubiquitin moieties in a polyubiquitin chain appear to be structurally equivalent to each other. It is therefore unclear how a specific ubiquitin moiety in a chain may be preferentially recognized by some proteins, such as the kinase PINK1. Here we show that there is structural dynamic heterogeneity in the two ubiquitin moieties of K48-linked diubiquitin by NMR spectroscopic analyses. Our analyses capture subunit-asymmetric structural fluctuations that are not directly related to the closed-to-open transition of the two ubiquitin moieties in diubiquitin. Strikingly, these newly identified heterogeneous structural fluctuations may be linked to an increase in susceptibility to phosphorylation by PINK1. Coupled with the fact that there are almost no differences in static tertiary structure among ubiquitin moieties in a chain, the observed subunit-specific structural fluctuations may be an important factor that distinguishes individual ubiquitin moieties in a chain, thereby aiding both efficiency and specificity in post-translational modifications.

Molecular insights into the differential structure-dynamics-stability features of interleukin-8 orthologs: Implications to functional specificity

Chemokines are a sub-group of chemotactic cytokines that regulate the leukocyte migration by binding to G-protein coupled receptors (GPCRs) and cell surface glycosaminoglycans (GAGs). Interleukin-8 (CXCL8/IL8) is one of the most essential CXC chemokine that has been reported to be involved in various pathophysiological conditions. Structure-function relationships of human IL8 have been studied extensively. However, no such detailed information is available on IL8 orthologs, although they exhibit significant functional divergence. In order to unravel the differential structure-dynamics-stability-function relationship of IL8 orthologs, comparative molecular analysis was performed on canine (laurasians) and human (primates) IL8 proteins using in-silico molecular evolutionary analysis and solution NMR spectroscopy methods. The residue level NMR studies suggested that, although the overall structural architecture of canine IL8 is similar to that of human IL8, systematic differences were observed in their backbone dynamics and low-energy excited states due to amino acid substitutions. Further, these substitutions also resulted in attenuation of stability and heparin binding affinity in the canine IL8 as compared to its human counterpart. Indeed, structural and sequence analysis evidenced for specificity of molecular interactions with cognate receptor (CXCR1) and glycosaminoglycan (heparin), thus providing evidence for a noticeable functional specificity and divergence between the two IL8 orthologs.

Evidence that the TRPV1 S1-S4 membrane domain contributes to thermosensing

Sensing and responding to temperature is crucial in biology. The TRPV1 ion channel is a well-studied heat-sensing receptor that is also activated by vanilloid compounds, including capsaicin. Despite significant interest, the molecular underpinnings of thermosensing have remained elusive. The TRPV1 S1-S4 membrane domain couples chemical ligand binding to the pore domain during channel gating. Here we show that the S1-S4 domain also significantly contributes to thermosensing and couples to heat-activated gating. Evaluation of the isolated human TRPV1 S1-S4 domain by solution NMR, far-UV CD, and intrinsic fluorescence shows that this domain undergoes a non-denaturing temperature-dependent transition with a high thermosensitivity. Further NMR characterization of the temperature-dependent conformational changes suggests the contribution of the S1-S4 domain to thermosensing shares features with known coupling mechanisms between this domain with ligand and pH activation. Taken together, this study shows that the TRPV1 S1-S4 domain contributes to TRPV1 temperature-dependent activation.

Conformational heterogeneity of Savinase from NMR, HDX-MS and X-ray diffraction analysis

Background

Several examples have emerged of enzymes where slow conformational changes are of key importance for function and where low populated conformations in the resting enzyme resemble the conformations of intermediate states in the catalytic process. Previous work on the subtilisin protease, Savinase, from Bacillus lentus by NMR spectroscopy suggested that this enzyme undergoes slow conformational dynamics around the substrate binding site. However, the functional importance of such dynamics is unknown.

Methods

Here we have probed the conformational heterogeneity in Savinase by following the temperature dependent chemical shift changes. In addition, we have measured changes in the local stability of the enzyme when the inhibitor phenylmethylsulfonyl fluoride is bound using hydrogen-deuterium exchange mass spectrometry (HDX-MS). Finally, we have used X-ray crystallography to compare electron densities collected at cryogenic and ambient temperatures and searched for possible low populated alternative conformations in the crystals.

Results

The NMR temperature titration shows that Savinase is most flexible around the active site, but no distinct alternative states could be identified. The HDX shows that modification of Savinase with inhibitor has very little impact on the stability of hydrogen bonds and solvent accessibility of the backbone. The most pronounced structural heterogeneities detected in the diffraction data are limited to alternative side-chain rotamers and a short peptide segment that has an alternative main-chain conformation in the crystal at cryo conditions. Collectively, our data show that there is very little structural heterogeneity in the resting state of Savinase and hence that Savinase does not rely on conformational selection to drive the catalytic process.

Temperature dependence of NMR chemical shifts: Tracking and statistical analysis

Isotropic chemical shifts measured by solution nuclear magnetic resonance (NMR) spectroscopy offer extensive insights into protein structure and dynamics. Temperature dependences add a valuable dimension; notably, the temperature dependences of amide proton chemical shifts are valuable probes of hydrogen bonding, temperature‐dependent loss of structure, and exchange between distinct protein conformations. Accordingly, their uses include structural analysis of both folded and disordered proteins, and determination of the effects of mutations, binding, or solution conditions on protein energetics. Fundamentally, these temperature dependences result from changes in the local magnetic environments of nuclei, but correlations with global thermodynamic parameters measured via calorimetric methods have been observed. Although the temperature dependences of amide proton and nitrogen chemical shifts are often well approximated by a linear model, deviations from linearity are also observed and may be interpreted as evidence of fast exchange between distinct conformational states. Here, we describe computational methods, accessible via the Shift‐T web server, including an automated tracking algorithm that propagates initial (single temperature) ¹H—¹⁵N cross peak assignments to spectra collected over a range of temperatures. Amide proton and nitrogen temperature coefficients (slopes determined by fitting chemical shift vs. temperature data to a linear model) are subsequently calculated. Also included are methods for the detection of systematic, statistically significant deviation from linearity (curvature) in the temperature dependences of amide proton chemical shifts. The use and utility of these methods are illustrated by example, and the Shift‐T web server is freely available at http://meieringlab.uwaterloo.ca/shiftt.

On identifying low energy conformational excited states with differential ruggedness in human γS-crystallin promoting severe infantile cataracts

Transient excited states in proteins can be accurately probed from temperature dependence of amide proton (¹H^N) chemical shifts displaying significant curvatures. Characterizing these near-native alternative states is of high therapeutic relevance in conformational diseases wherein missense mutations promote structural instability that leads to conformational heterogeneity. Extending the structure-function paradigm from physiology to pathology, we recently reported the solution NMR structure and dynamics of a severe congenital cataract variant, G57W of human γS-crystallin (abbreviated as γS-G57W) which is resistant towards crystallization. In an endeavour to explore the functional consequences of this mutation, here we report for the first time, native state conformational ruggedness in human γS-G57W as compared to its wild-type counterpart from residue resolved nonlinear temperature dependence of backbone ¹H^N chemical shifts using solution NMR spectroscopy. Our calculations suggest that the simulated chemical shift curvatures are indicative of low energy excited states within 2-4 kcal mol^-1 from the native state. Residues accessing alternative conformations populate the N-terminal domain of γS-G57W more than its C-terminal counterpart. Collectively, curvatures retaining native state ensemble on mild denaturation suggest that the free energy landscape in human γS-G57W at the bottom of the folding funnel is sufficiently robust and malleable against such perturbations. Overall, this critical study highlights the functional aspects of such structural malleability promoting infantile cataracts as a global health risk marker.

NMR characterization of conformational fluctuations and noncovalent interactions of SUMO protein from Drosophila melanogaster (dSmt3)

Structural heterogeneity in the native-state ensemble of dSmt3, the only small ubiquitin-like modifier (SUMO) in Drosophila melanogaster, was investigated and compared with its human homologue SUMO1. Temperature dependence of amide proton's chemical shift was studied to identify amino acids possessing alternative structural conformations in the native state. Effect of small concentration of denaturant (1M urea) on this population was also monitored to assess the ruggedness of near-native energy landscape. Owing to presence of many such amino acids, especially in the β₂ -loop-α region, the native state of dSmt3 seems more flexible in comparison to SUMO1. Information about backbone dynamics in ns-ps timescale was quantified from the measurement of ¹⁵ N-relaxation experiments. Furthermore, the noncovalent interaction of dSmt3 and SUMO1 with Daxx12 (Daxx⁷²⁹ DPEEIIVLSDSD⁷⁴⁰ ), a [V/I]-X-[V/I]-[V/I]-based SUMO interaction motif, was characterized using Bio-layer Interferometery and NMR spectroscopy. Daxx12 fits itself in the groove formed by β₂ -loop-α structural region in both dSmt3 and SUMO1, but the binding is stronger with the former. Flexibility of β₂ -loop-α region in dSmt3 is suspected to assist its interaction with Daxx12. Our results highlight the role of native-state flexibility in assisting noncovalent interactions of SUMO proteins especially in organisms where a single SUMO isoform has to tackle multiple substrates single handedly.

Why the Energy Landscape of Barnase Is Hierarchical

We have used NMR and computational methods to characterize the dynamics of the ribonuclease barnase over a wide range of timescales in free and inhibitor-bound states. Using temperature- and denaturant-dependent measurements of chemical shift, we show that barnase undergoes frequent and highly populated hinge bending. Using relaxation dispersion, we characterize a slower and less populated motion with a rate of 750 ± 200 s⁻¹, involving residues around the lip of the active site, which occurs in both free and bound states and therefore suggests conformational selection. Normal mode calculations characterize correlated hinge bending motions on a very rapid timescale. These three measurements are combined with previous measurements and molecular dynamics calculations on barnase to characterize its dynamic landscape on timescales from picoseconds to milliseconds and length scales from 0.1 to 2.5 nm. We show that barnase has two different large-scale fluctuations: one on a timescale of 10⁻⁹−10⁻⁶ s that has no free energy barrier and is a hinge bending that is determined by the architecture of the protein; and one on a timescale of milliseconds (i.e., 750 s⁻¹) that has a significant free energy barrier and starts from a partially hinge-bent conformation. These two motions can be described as hierarchical, in that the more highly populated faster motion provides a platform for the slower (less probable) motion. The implications are discussed. The use of temperature and denaturant is suggested as a simple and general way to characterize motions on the intermediate ns-μs timescale.

Characterization of low-lying excited states of proteins by high-pressure NMR

Hydrostatic pressure alters the free energy of proteins by a few kJ mol^-1, with the amount depending on their partial molar volumes. Because the folded ground state of a protein contains cavities, it is always a state of large partial molar volume. Therefore pressure always destabilises the ground state and increases the population of partially and completely unfolded states. This is a mild and reversible conformational change, which allows the study of excited states under thermodynamic equilibrium conditions. Many of the excited states studied in this way are functionally relevant; they also seem to be very similar to kinetic folding intermediates, thus suggesting that evolution has made use of the 'natural' dynamic energy landscape of the protein fold and sculpted it to optimise function. This includes features such as ligand binding, structural change during the catalytic cycle, and dynamic allostery.

Concurrent Increases and Decreases in Local Stability and Conformational Heterogeneity in Cu, Zn Superoxide Dismutase Variants Revealed by Temperature-Dependence of Amide Chemical Shifts

The chemical shifts of backbone amide protons in proteins are sensitive reporters of local structural stability and conformational heterogeneity, which can be determined from their readily measured linear and nonlinear temperature-dependences, respectively. Here we report analyses of amide proton temperature-dependences for native dimeric Cu, Zn superoxide dismutase (holo pWT SOD1) and structurally diverse mutant SOD1s associated with amyotrophic lateral sclerosis (ALS). Holo pWT SOD1 loses structure with temperature first at its periphery and, while having extremely high global stability, nevertheless exhibits extensive conformational heterogeneity, with ∼1 in 5 residues showing evidence for population of low energy alternative states. The holo G93A and E100G ALS mutants have moderately decreased global stability, whereas V148I is slightly stabilized. Comparison of the holo mutants as well as the marginally stable immature monomeric unmetalated and disulfide-reduced (apo(2SH)) pWT with holo pWT shows that changes in the local structural stability of individual amides vary greatly, with average changes corresponding to differences in global protein stability measured by differential scanning calorimetry. Mutants also exhibit altered conformational heterogeneity compared to pWT. Strikingly, substantial increases as well as decreases in local stability and conformational heterogeneity occur, in particular upon maturation and for G93A. Thus, the temperature-dependence of amide shifts for SOD1 variants is a rich source of information on the location and extent of perturbation of structure upon covalent changes and ligand binding. The implications for potential mechanisms of toxic misfolding of SOD1 in disease and for general aspects of protein energetics, including entropy-enthalpy compensation, are discussed.

Pressure-Dependent Conformation and Fluctuation in Folded Protein Molecules

Hydrostatic pressure leads to nonuniform compression of proteins. The structural change is on average only about 0.1 Å kbar(-1), and is therefore within the range of fluctuations at ambient pressure. The largest changes are around cavities and buried water molecules. Sheets distort much more than helices. Hydrogen bonds compress about 0.012 Å kbar(-1), although there is a wide range, including some hydrogen bonds that lengthen. In the presence of ligands and inhibitors, structural changes are smaller. Pressure has little effect on rapid fluctuations, but with larger scale slower motions, pressure increases the population of excited states (if they have smaller overall volume), and slows the fluctuations. In barnase, pressure is shown to be a useful way to characterise fluctuations on the timescale of microseconds, and helps to show that fluctuations in barnase are hierarchical, with the faster fluctuations providing a platform for the slower ones. The excited states populated at high pressure are probably functionally important.

Characterization of the near native conformational states of the SAM domain of Ste11 protein by NMR spectroscopy

The sterile alpha motif or SAM domain is one of the most frequently present protein interaction modules with diverse functional attributions. SAM domain of the Ste11 protein of budding yeast plays important roles in mitogen-activated protein kinase cascades. In the current study, urea-induced, at subdenaturing concentrations, structural, and dynamical changes in the Ste11 SAM domain have been investigated by nuclear magnetic resonance spectroscopy. Our study revealed that a number of residues from Helix 1 and Helix 5 of the Ste11 SAM domain display plausible alternate conformational states and largest chemical shift perturbations at low urea concentrations. Amide proton (H/D) exchange experiments indicated that Helix 1, loop, and Helix 5 become more susceptible to solvent exchange with increased concentrations of urea. Notably, Helix 1 and Helix 5 are directly involved in binding interactions of the Ste11 SAM domain. Our data further demonstrate that the existence of alternate conformational states around the regions involved in dimeric interactions in native or near native conditions.

NMR characterization of the near native and unfolded states of the PTB domain of Dok1: alternate conformations and residual clusters

BACKGROUND

Phosphotyrosine binding (PTB) domains are critically involved in cellular signaling and diseases. PTB domains are categorized into three distinct structural classes namely IRS-like, Shc-like and Dab-like. All PTB domains consist of a core pleckstrin homology (PH) domain with additional structural elements in Shc and Dab groups. The core PH fold of the PTB domain contains a seven stranded β-sheet and a long C-terminal helix.

PRINCIPAL FINDINGS

In this work, the PTB domain of Dok1 protein has been characterized, by use of NMR spectroscopy, in solutions containing sub-denaturing and denaturing concentrations of urea. We find that the Dok1 PTB domain displays, at sub-denaturing concentrations of urea, alternate conformational states for residues located in the C-terminal helix and in the β5 strand of the β-sheet region. The β5 strand of PTB domain has been found to be experiencing significant chemical shift perturbations in the presence of urea. Notably, many of these residues in the helix and in the β5 strand are also involved in ligand binding. Structural and dynamical analyses at 7 M urea showed that the PTB domain is unfolded with islands of motionally restricted regions in the polypeptide chain. Further, the C-terminal helix appears to be persisted in the unfolded state of the PTB domain. By contrast, residues encompassing β-sheets, loops, and the short N-terminal helix lack any preferred secondary structures. Moreover, these residues demonstrated an intimate contact with the denaturant.

SIGNIFICANCE

This study implicates existence of alternate conformational states around the ligand binding pocket of the PTB domain either in the native or in the near native conditions. Further, the current study demonstrates that the C-terminal helical region of PTB domain may be considered as a potential site for the initiation of folding.

Conformational heterogeneity and dynamics in a βγ-crystallin from Hahella chejuensis

Most of the βγ-crystallins are structural proteins with high intrinsic stability, which gets enhanced by Ca(2+)-binding in microbial members. Functions of most of these proteins are yet to be known. However, a few of them were reported to be involved in Ca(2+)-dependent and stress-related functions. Hahellin, a microbial homolog, is a natively unfolded protein that acquires a well-folded structure upon Ca(2+) binding. Although the structure of βγ-crystallin domains is well understood, the dynamical features are yet to be explored. We have investigated for the first time the equilibrium dynamics, conformational heterogeneity and associated low-lying free-energy states of hahellin in its Ca(2+)-bound form using NMR spectroscopy to understand the dynamics of a βγ-crystallin domain. Hahellin shows large conformational heterogeneity with nearly 40% of the residues, some of which are part of Ca(2+)-binding loops, accessing alternative states. Further, out of the two Greek key motifs, which together constitute the βγ-crystallin domain, the second Greek key motif is floppy as compared to its relatively rigid counterpart. Taken together, we believe that these characteristics might be of importance to understand the stability and functions of βγ-crystallin domains.

Local cooperativity in an amyloidogenic state of human lysozyme observed at atomic resolution

The partial unfolding of human lysozyme underlies its conversion from the soluble state into amyloid fibrils observed in a fatal hereditary form of systemic amyloidosis. To understand the molecular origins of the disease, it is critical to characterize the structural and physicochemical properties of the amyloidogenic states of the protein. Here we provide a high-resolution view of the unfolding process at low pH for three different lysozyme variants, the wild-type protein and the mutants I56T and I59T, which show variable stabilities and propensities to aggregate in vitro. Using a range of biophysical techniques that includes differential scanning calorimetry and nuclear magnetic resonance spectroscopy, we demonstrate that thermal unfolding under amyloidogenic solution conditions involves a cooperative loss of native tertiary structure, followed by progressive unfolding of a compact, molten globule-like denatured state ensemble as the temperature is increased. The width of the temperature window over which the denatured ensemble progressively unfolds correlates with the relative amyloidogenicity and stability of these variants, and the region of lysozyme that unfolds first maps to that which forms the core of the amyloid fibrils formed under similar conditions. Together, these results present a coherent picture at atomic resolution of the initial events underlying amyloid formation by a globular protein.

Characterization of low-energy excited states in the native state ensemble of non-myristoylated and myristoylated neuronal calcium sensor-1

Information on the low-energy excited states of a given protein is important as this controls the structural adaptability and various biological functions of proteins such as co-operativity, response towards various external perturbations. In this article, we characterized individual residues in both non-myristoylated (non-myr) and myristoylated (myr) neuronal calcium sensor-1 (NCS-1) that access alternate states by measuring nonlinear temperature dependence of the backbone amide-proton (¹H(N)) chemical shifts. We found that ~20% of the residues in the protein access alternative conformations in non-myr case, which increases to ~28% for myr NCS-1. These residues are spread over the entire polypeptide stretch and include the edges of α-helices and β-strands, flexible loop regions, and the Ca²(+)-binding loops. Besides, residues responsible for the absence of Ca²(+)-myristoyl switch are also found accessing alternative states. The C-terminal domain is more populated with these residues compared to its N-terminal counterpart. Individual EF-hands in NCS-1 show significantly different number of alternate states. This observation prompts us to conclude that this may lead to differences in their individual conformational flexibility and has implications on the functionality. Theoretical simulations reveal that these low-energy excited states are within an energy band of 2-4 kcal/mol with respect to the native state.

Structural change in a B-DNA helix with hydrostatic pressure

Study of the effects of pressure on macromolecular structure improves our understanding of the forces governing structure, provides details on the relevance of cavities and packing in structure, increases our understanding of hydration and provides a basis to understand the biology of high-pressure organisms. A study of DNA, in particular, helps us to understand how pressure can affect gene activity. Here we present the first high-resolution experimental study of B-DNA structure at high pressure, using NMR data acquired at pressures up to 200 MPa (2 kbar). The structure of DNA compresses very little, but is distorted so as to widen the minor groove, and to compress hydrogen bonds, with AT pairs compressing more than GC pairs. The minor groove changes are suggested to lead to a compression of the hydration water in the minor groove.

Differential native state ruggedness of the two Ca2+-binding domains in a Ca2+ sensor protein

Characterization of near-native excited states of a protein provides insights into various biological functions such as co-operativity, protein-ligand, and protein-protein interactions. In the present study, we investigated the ruggedness of the native state of EhCaBP using nonlinear temperature dependence of backbone amide-proton chemical shifts. EhCaBP is a two-domain EF-hand calcium sensor protein consisting of two EF-hands in each domain and binds four Ca2+ ions. It has been observed that approximately 30% of the residues in the protein access alternative conformations. Theoretical modeling suggested that these low-energy excited states are within 2-3 kcal/mol from the native state. Further, it is interesting to note that the residues accessing alternative conformations are more dominated in the C-terminal domain compared with its N-terminal counterpart suggesting that the former is more rugged in its native state. These distinct characteristics of N- and C-terminal domains of a calcium sensor protein belonging to the super family of calmodulin would have implications for domain dependent Ca2+ signaling pathways.

Structural and thermodynamic effects of ANS binding to human interleukin-1 receptor antagonist

Although 8-anilinonaphthalene-1-sulfonic acid (ANS) is frequently used in protein folding studies, the structural and thermodynamic effects of its binding to proteins are not well understood. Using high-resolution two-dimensional NMR and human interleukin-1 receptor antagonist (IL-1ra) as a model protein, we obtained detailed information on ANS-protein interactions in the absence and presence of urea. The effects of ambient to elevated temperatures on the affinity and specificity of ANS binding were assessed from experiments performed at 25 degrees C and 37 degrees C. Overall, the affinity of ANS was lower at 37 degrees C compared to 25 degrees C, but no significant change in the site specificity of binding was observed from the chemical shift perturbation data. The same site-specific binding was evident in the presence of 5.2 M urea, well within the unfolding transition region, and resulted in selective stabilization of the folded state. Based on the two-state denaturation mechanism, ANS-dependent changes in the protein stability were estimated from relative intensities of two amide resonances specific to the folded and unfolded states of IL-1ra. No evidence was found for any ANS-induced partially denatured or aggregated forms of IL-1ra throughout the experimental conditions, consistent with a cooperative and reversible denaturation process. The NMR results support earlier observations on the tendency of ANS to interact with solvent-exposed positively charged sites on proteins. Under denaturing conditions, ANS binding appears to be selective to structured states rather than unfolded conformations. Interestingly, the binding occurs within a previously identified aggregation-critical region in IL-1ra, thus providing an insight into ligand-dependent protein aggregation.

NMR characterization of the energy landscape of SUMO-1 in the native-state ensemble

Characterizing the low energy excited states in the energy landscape of a protein is one of the exciting and demanding problems in structural biology at the present time. These describe the adaptability of the protein structure to external perturbations. In this context, we used here non-linear dependence of amide proton chemical shifts on temperature to identify residues accessing alternative conformations in SUMO-1 in the native state as well as in the near-native states created by sub-denaturing concentrations of urea. The number of residues accessing alternative conformations increases and the profiles of curved temperature dependence also change with increasing urea concentration. In every case these alternative conformations lie within 2 kcal/mol from the ground state, and are separated from it by low energy barriers. The residues that access alternative conformations span the length of the protein chain but are located at particular regions on the protein structure. These include many of the loops, beta2 and beta5 strands, and some edges of the helices. We observed that some of the regions of the protein structure that exhibit such fluctuations coincide with the protein's binding surfaces with different substrate like GTPase effector domain (GED) of dynamin, SUMO binding motifs (SBM), E1 (activating enzyme, SAE1/SAE2) and E2 (conjugating enzyme, UBC9) enzymes of sumoylation machinery, reported earlier. We speculate that this would have significant implications for the binding of diversity of targets by SUMO-1 for the variety of functions it is involved in.

Conformational stability and dynamics of cytochrome c affect its alkaline isomerization

The alkaline isomerization of horse heart ferricytochrome c (cyt c) has been studied by electronic absorption spectroscopy in the presence of the Hofmeister series of anions: chloride, bromide, rhodanide and perchlorate. The anions significantly affect the apparent pK (a) value of the transition in a concentration-dependent manner according to their position in the Hofmeister series. The Soret region of the absorption spectra is not affected by the presence of the salts and shows no significant structural perturbation of the heme crevice. In the presence of perchlorate and rhodanide anions, the cyanide exchange rate between the bulk solvent and the binding site is increased. These results imply higher flexibility of the protein structure in the presence of chaotropic salts. The thermal and isothermal denaturations monitored by differential scanning calorimetry and circular dichroism, respectively, showed a decrease in the conformational stability of cyt c in the presence of the chaotropic salts. A positive correlation between the stability, DeltaG, of cyt c and the apparent pK (a) values that characterize the alkaline transition indicates the presence of a thermodynamic linkage between these conformational transitions. In addition, the rate constant of the cyanide binding and the partial molar entropies of anions negatively correlate with the pK (a) values. This indicates the important role of anion-induced solvent reorganization on the structural flexibility of cyt c in the alkaline transitions.

Heteronuclear NMR Identifies a Nascent Helix in Intrinsically Disordered Dynein Intermediate Chain: Implications for Folding and Dimerization

The intermediate chain of dynein forms a tight subcomplex with dimeric light chains LC8 and Tctex-1, and together they constitute the cargo attachment complex. There is considerable interest in identifying the role of these light chains in the assembly of the two copies of the intermediate chain. The N-terminal domain of the intermediate chain, IC1-289, contains the binding sites for the light chains, and is a highly disordered monomer but gains helical structure upon binding to light chains LC8 and Tctex-1. To provide insights into the structural and dynamic changes that occur in the intermediate chain upon light chains binding, we have used NMR spectroscopy to compare the properties of two distinct sub-domains of IC1-289: IC84-143 which is the light chains binding domain, and IC198-237, which contains a predicted coiled coil necessary for the increase in ordered structure upon light chain binding. Neither construct has stable secondary structure when probed by circular dichroism and amide chemical shift dispersion. Specific residues of IC84-143 involved in binding to the light chains were identified by their increase in resonance line broadening and the corresponding large intensity reduction in 1H-15N HSQC spectra. Interestingly, IC84-143 shows no sign of structure formation after binding to either LC8 or Tctex-1 or to both. IC198-237, on the other hand, contains a population of a nascent helix at low temperature as identified by heteronuclear NMR relaxation measurements, secondary chemical shifts, and sequential amide-amide connectivities. These data are consistent with a model for light chain binding coupled to intermediate chain dimerization through forming a coiled coil distant from the binding site.

The inverted chevron plot measured by NMR relaxation reveals a native-like unfolding intermediate in acyl-CoA binding protein

The folding kinetics of bovine acyl-CoA binding protein was studied by 15N relaxation dispersion measurements under equilibrium conditions. Relaxation dispersion profiles were measured at several concentrations of guanidine hydrochloride (GuHCl). The unfolding rate constant (k(u)) was determined under conditions favoring folding, for which the folding rate constant (k(f)) dominates the relaxation in stopped-flow kinetic measurements. Conversely, k(f) was determined under conditions favoring unfolding, for which k(u) dominates stopped-flow data. The rates determined by NMR therefore complement those from stopped-flow kinetics and define an "inverted chevron" plot. The combination of NMR relaxation and stopped-flow kinetic measurements allowed determination of k(f) and k(u) in the range from 0.48 M GuHCl to 1.28 M GuHCl. Individually, the stopped-flow and NMR data fit two-state models for folding. However, although the values of k(f) determined by the two methods agree, the values of k(u) do not. As a result, a combined analysis of all data does not comply with a two-state model but indicates that an unfolding intermediate exists on the native side of the dominant energy barrier. The denaturant and temperature dependencies of the chemical shifts and k(u) indicate that the intermediate state is structurally similar to the native state. Equilibrium unfolding monitored by optical spectroscopy corroborate these conclusions. The temperature dependence of the chemical shifts identifies regions of the protein that are selectively destabilized in the intermediate. These results illustrate the power of combining stopped-flow kinetics and NMR spectroscopy to analyze protein folding.

An Experimental Investigation of Conformational Fluctuations in Proteins G and L

The B1 domains of streptococcal proteins G and L are structurally similar, but they have different sequences and they fold differently. We have measured their NMR spectra at variable temperature using a range of concentrations of denaturant. Many residues have curved amide proton temperature dependence, indicating that they significantly populate alternative, locally unfolded conformations. The results, therefore, provide a view of the locations of low-lying, locally unfolded conformations. They indicate approximately 4-6 local minima for each protein, all within ca. 2.5 kcal/mol of the native state, implying a locally rough energy landscape. Comparison with folding data for these proteins shows that folding involves most molecules traversing a similar path, once a transition state containing a beta hairpin has been formed, thereby defining a well-populated pathway down the folding funnel. The hairpin that directs the folding pathway differs for the two proteins and remains the most stable part of the folded protein.

Combining prediction, computation and experiment for the characterization of protein disorder

Several computational and experimental methods exist for identifying disordered residues within proteins. Computational algorithms can now identify these disordered sequences and predict their occurrence within genomes with relatively high accuracy. Recent advances in NMR and mass spectroscopy permit faster and more detailed studies of disordered states at atomic resolutions. Combining prediction, computation and experimentation is proposed to accelerate and enhance the characterization of intrinsically disordered protein.