Protein hydrogen exchange mechanism: local fluctuations

Deep Eutectic Solvent-Based Highly Sensitive Turn-On Fluorescent Probe for D2O

Owing to the importance of heavy water in spectroscopy, nuclear energy generation, chemical characterization, and biological industry, a design of a robust, cheap, nontoxic, and sensitive D2O sensor is very important. In this work, taking advantage of the singular emission fluorescence of the deep eutectic solvent prepared in our laboratory, we propose a first of its kind highly sensitive turn-on fluorescent sensor to effectively sense D2O at an ultratrace level based on rapid exchange of the labile DES proton with deuterium. This method can be used as a full-range heavy water detection strategy with a limit of detection of 0.079% (v/v) or 870 ppm. The isotopic purity (IP) obtained from DES fluorescence measurements is also in close agreement with that of the conventional FT-IR method. The current DES-based sensor thus allows both sensing and isotopic purity of D2O and can serve as one of the most sensitive monitoring strategies for heavy water analysis.

Molecular insights on the conformational dynamics of a P76C mutant of human cytochrome c and the enhancement on its peroxidase activity

In apoptotic pathway, the interaction of Cytochrome c (Cytc) with cardiolipin in vivo is a key process to induce peroxidase activity of Cytc and trigger the release of Cytc in the inner mitochondria into cytosol. The peroxidase active form of Cytc occurs due to local conformational changes that support the opening of the heme crevice and the loss of an axial ligand between Met80 and heme Fe. Structural adjustments at the Ω-loop segments of Cytc are required for such process. To study the role of the distal Ω-loop segments comprising residues 71-85 in human Cytc (hCytc), we investigated a cysteine mutation at Pro76, one of the highly conserved residues in this loop. The effect of P76C mutant was explored by the combination of experimental characterizations and molecular dynamics (MD) simulations. The peroxidase activity of the P76C mutant was found to be significantly increased by ∼13 folds relative to the wild type. Experimental data on global denaturation, alkaline transition, heme bleaching, and spin-labeling Electron Spin Resonance were in good agreement with the enhancement of peroxidase activity. The MD results of hCytc in the hexacoordinate form suggest the important changes in P76C mutant occurred due to the unfolding at the central Ω-loop (residues 40-57), and the weakening of H-bond between Tyr67 and Met80. Whereas the experimental data implied that the P76C mutant tend to be in equilibrium between the pentacoordinate and hexacoordinate forms, the MD and experimental information are complementary and were used to support the mechanisms of peroxidase active form of hCytc.

The contribution of electrostatics to hydrogen exchange in the unfolded protein state

Although electrostatics have long been recognized to play an important role in hydrogen exchange (HX) with solvent, the quantitative assessment of its magnitude in the unfolded state has hitherto been lacking. This limits the utility of HX as a quantitative method to study protein stability, folding, and dynamics. Using the intrinsically disordered human protein α-synuclein as a proxy for the unfolded state, we show that a hybrid mean-field approach can effectively compute the electrostatic potential at all backbone amide positions along the chain. From the electrochemical potential, a fourfold reduction in hydroxide concentration near the protein backbone is predicted for the C-terminal domain, a prognosis that is in direct agreement with experimentally derived protection factors from NMR spectroscopy. Thus, impeded HX for the C-terminal region of α-synuclein is not the result of intramolecular hydrogen bonding and/or structure formation.

Activation loop dynamics are controlled by conformation-selective inhibitors of ERK2

Conformational selection by small molecules expands inhibitory possibilities for protein kinases. Nuclear magnetic resonance (NMR) measurements of the mitogen-activated protein (MAP) kinase ERK2 have shown that activation by dual phosphorylation induces global motions involving exchange between two states, L and R. We show that ERK inhibitors Vertex-11e and SCH772984 exploit the small energetic difference between L and R to shift the equilibrium in opposing directions. An X-ray structure of active 2P-ERK2 complexed with AMP-PNP reveals a shift in the Gly-rich loop along with domain closure to position the nucleotide in a more catalytically productive conformation relative to inactive 0P-ERK2:ATP. X-ray structures of 2P-ERK2 complexed with Vertex-11e or GDC-0994 recapitulate this closure, which is blocked in a complex with a SCH772984 analog. Thus, the L→R shift in 2P-ERK2 is associated with movements needed to form a competent active site. Solution measurements by hydrogen-exchange mass spectrometry (HX-MS) reveal distinct binding interactions for Vertex-11e, GDC-0994, and AMP-PNP with active vs. inactive ERK2, where the extent of HX protection correlates with R state formation. Furthermore, Vertex-11e and SCH772984 show opposite effects on HX near the activation loop. Consequently, these inhibitors differentially affect MAP kinase phosphatase activity toward 2P-ERK2. We conclude that global motions in ERK2 reflect conformational changes at the active site that promote productive nucleotide binding and couple with changes at the activation loop to allow control of dephosphorylation by conformationally selective inhibitors.

Characterization of H/D exchange in type 1 pili by proton-detected solid-state NMR and molecular dynamics simulations

Uropathogenic Escherichia coli invades and colonizes hosts by attaching to cells using adhesive pili on the bacterial surface. Although many biophysical techniques have been used to study the structure and mechanical properties of pili, many important details are still unknown. Here we use proton-detected solid-state NMR experiments to investigate solvent accessibility and structural dynamics. Deuterium back-exchange at labile sites of the perdeuterated, fully proton back-exchanged pili was conducted to investigate hydrogen/deuterium (H/D) exchange patterns of backbone amide protons in pre-assembled pili. We found distinct H/D exchange patterns in lateral and axial intermolecular interfaces in pili. Amide protons protected from H/D exchange in pili are mainly located in the core region of the monomeric subunit and in the lateral intermolecular interface, whereas the axial intermolecular interface and the exterior region of pili are highly exposed to H/D exchange. Additionally, we performed molecular dynamics simulations of the type 1 pilus rod and estimated the probability of H/D exchange based on hydrogen bond dynamics. The comparison of the experimental observables and simulation data provides insights into stability and mechanical properties of pili.

Networks of electrostatic and hydrophobic interactions modulate the complex folding free energy surface of a designed βα protein

The successful de novo design of proteins can provide insights into the physical chemical basis of stability, the role of evolution in constraining amino acid sequences, and the production of customizable platforms for engineering applications. Previous guanidine hydrochloride (GdnHCl; an ionic denaturant) experiments of a designed, naturally occurring βα fold, Di-III_14, revealed a cooperative, two-state unfolding transition and a modest stability. Continuous-flow mixing experiments in our laboratory revealed a simple two-state reaction in the microsecond to millisecond time range and consistent with the thermodynamic results. In striking contrast, the protein remains folded up to 9.25 M in urea, a neutral denaturant, and hydrogen exchange (HDX) NMR analysis in water revealed the presence of numerous high-energy states that interconvert on a time scale greater than seconds. The complex protection pattern for HDX corresponds closely with a pair of electrostatic networks on the surface and an extensive network of hydrophobic side chains in the interior of the protein. Mutational analysis showed that electrostatic and hydrophobic networks contribute to the resistance to urea denaturation for the WT protein; remarkably, single charge reversals on the protein surface restore the expected urea sensitivity. The roughness of the energy surface reflects the densely packed hydrophobic core; the removal of only two methyl groups eliminates the high-energy states and creates a smooth surface. The design of a very stable βα fold containing electrostatic and hydrophobic networks has created a complex energy surface rarely observed in natural proteins.

CcdB at pH 4 Forms a Partially Unfolded State with a Dry Core

pH is an important factor that affects the protein structure, stability, and activity. Here, we probe the nature of the low-pH structural form of the homodimeric CcdB (controller of cell death B) protein. Characterization of CcdB protein at pH 4 and 300 K using circular dichroism spectroscopy, 8-anilino-1-naphthalene-sulphonate binding, and Trp solvation studies suggests that it forms a partially unfolded state with a dry core at equilibrium under these conditions. CcdB remains dimeric at pH 4 as shown by multiple techniques, such as size-exclusion chromatography coupled to multiangle light scattering, analytical ultracentrifugation, and electron paramagnetic resonance. Comparative analysis using two-dimensional ¹⁵N-¹H heteronuclear single-quantum coherence NMR spectra of CcdB at pH 4 and 7 suggests that the pH 4 and native state have similar but nonidentical structures. Hydrogen-exchange-mass-spectrometry studies demonstrate that the pH 4 state has substantial but anisotropic changes in local stability with core regions close to the dimer interface showing lower protection but some other regions showing higher protection relative to pH 7.

Interpreting Hydrogen-Deuterium Exchange Events in Proteins Using Atomistic Simulations: Case Studies on Regulators of G-Protein Signaling Proteins

Hydrogen-deuterium exchange (HDX) experiments are widely used in studies of protein dynamics. To predict the propensity of amide hydrogens for exchange with deuterium, several models have been reported in which computations of amide-hydrogen protection factors are carried out using molecular dynamics (MD) simulations. Given significant variation in the criteria used in different models, the robustness and broader applicability of these models to other proteins, especially homologous proteins showing distinct amide-exchange patterns, remains unknown. The sensitivity of the predictions when MD simulations are conducted with different force-fields is yet to tested and quantified. Using MD simulations and experimental HDX data on three homologous signaling proteins, we report detailed studies quantifying the performance of seven previously reported models (M1-M7) of two general types: empirical and fractional-population models. We find that the empirical models show inconsistent predictions but predictions of the fractional population models are robust. Contrary to previously reported work, we find that the solvent-accessible surface area of amide hydrogens is a useful metric when combined with a new metric defining the distances of amide hydrogens from the first polar atoms in proteins. On the basis of this, we report two new models, one empirical (M8) and one population-based (M9). We find strong protection of amide hydrogens from solvent exchange both within the stable helical motifs and also in the interhelical loops. We further observe that the exchange-competent states of amide hydrogens occur on the sub 100 ps time-scale via localized fluctuations, and such states among amides of a given protein do not appear to show any cooperativity or allosteric coupling.

Biochemical and Biophysical Characterisation of Higher Oligomeric Structure of Rat Nucleosome Assembly Protein 1

Nucleosome assembly protein 1 (NAP1) is a histone chaperone that exchanges histone H2A-H2B dimer from chromatin templates. Studies with yeast NAP1 (yNAP1) have revealed its existence as multiple oligomeric species in solution. Here, rat NAP1 (rNAP1), which is 98% identical to the human NAP1 (hNAP1) was used as a model to characterize the oligomeric structures of this protein in higher eukaryotes. Gel filtration chromatography and Dynamic light scattering of recombinant rNAP1 indicated that the protein exists as a complex mixture of multimeric species even at 500 mM ionic strength. The solution-state complexity remains unchanged even at higher ionic strengths. Equilibrium unfolding (ΔG 14.6 kcal mol^- 1) shows that rNAP1, both dimeric and oligomeric, follow the two-state model of unfolding with no detectable intermediates. Homology modelling suggests that rat and yeast NAP1 share an overall similar structure with conserved domains. However, dissimilar substitutions like threonine and lysine with glycine in the β-hairpin involved in oligomerization, possibly leads to the observed differences in the oligomerization propensity of the two proteins. Molecular dynamic simulation (MDS) of the two structures also revealed that rNAP1 dimer is more stable owing to the extensive hydrogen bonding in comparison to yNAP1. Further, in vitro kinase assay showed that the phosphorylation of rNAP1 favors oligomerization with no effect on its histone binding capacity. Our results clearly suggest that there are differences in the in-solution behavior of rNAP1 compared to yNAP1 which may have in vivo functional implications for the regulation of these complexes during chromatin assembly and rearrangement.

Computational method allowing Hydrogen-Deuterium Exchange Mass Spectrometry at single amide Resolution

Hydrogen-deuterium exchange (HDX) coupled with mass spectrometry (HDXMS) is a rapid and effective method for localizing and determining protein stability and dynamics. Localization is routinely limited to a peptide resolution of 5 to 20 amino acid residues. HDXMS data can contain information beyond that needed for defining protein stability at single amide resolution. Here we present a method for extracting this information from an HDX dataset to generate a HDXMS protein stability fingerprint. High resolution (HR)-HDXMS was applied to the analysis of a model protein of a spectrin tandem repeat that exemplified an intuitive stability profile based on the linkage of two triple helical repeats connected by a helical linker. The fingerprint recapitulated expected stability maximums and minimums with interesting structural features that corroborate proposed mechanisms of spectrin flexibility and elasticity. HR-HDXMS provides the unprecedented ability to accurately assess protein stability at the resolution of a single amino acid. The determination of HDX stability fingerprints may be broadly applicable in many applications for understanding protein structure and function as well as protein ligand interactions.

Non-equilibrium hydrogen exchange for determination of H-bond strength and water accessibility in solid proteins

We demonstrate measurement of non-equilibrium backbone amide hydrogen-deuterium exchange rates (HDX) for solid proteins. The target of this study are the slowly exchanging residues in solid samples, which are associated with stable secondary-structural elements of proteins. These hydrogen exchange processes escape methods measuring equilibrium exchange rates of faster processes. The method was applied to a micro-crystalline preparation of the SH3 domain of chicken α-spectrin. Therefore, from a 100% back-exchanged micro-crystalline protein preparation, the supernatant buffer was exchanged by a partially deuterated buffer to reach a final protonation level of approximately 20% before packing the sample in a 1.3 mm rotor. Tracking of the HN peak intensities for 2 weeks reports on site-specific hydrogen bond strength and also likely reflects water accessibility in a qualitative manner. H/D exchange can be directly determined for hydrogen-bonded amides using ¹H detection under fast magic angle spinning. This approach complements existing methods and provides the means to elucidate interesting site-specific characteristics for protein functionality in the solid state.

Targeted mutagenesis of the P22 portal protein reveals the mechanism of signal transmission during DNA packaging

The portal vertex in dsDNA bacteriophage serves as the site for genome encapsidation and release. In several of these viruses, efficient termination of DNA packaging has been shown to be dependent on the density of packaged DNA. The portal protein has been implicated as being part of the sensor that regulates packaging termination through DNA-dependent conformational changes during packaging. The mechanism by which DNA induces these conformational changes remains unknown. In this study, we explore how point mutants in the portal core can result in changes in genome packaging density in P22. Mutations in the portal core that subtly alter the structure or dynamics of the protein result in an increase in the amount of DNA packaged. The magnitude of the change is amino acid and location specific. Our findings suggest a mechanism wherein compression of the portal core is an essential aspect of signal transmission during packaging.

A J-modulated protonless NMR experiment characterizes the conformational ensemble of the intrinsically disordered protein WIP

Intrinsically disordered proteins (IDPs) are multi-conformational polypeptides that lack a single stable three-dimensional structure. It has become increasingly clear that the versatile IDPs play key roles in a multitude of biological processes, and, given their flexible nature, NMR is a leading method to investigate IDP behavior on the molecular level. Here we present an IDP-tailored J-modulated experiment designed to monitor changes in the conformational ensemble characteristic of IDPs by accurately measuring backbone one- and two-bond J(¹⁵N,¹³Cα) couplings. This concept was realized using a unidirectional (H)NCO ¹³C-detected experiment suitable for poor spectral dispersion and optimized for maximum coverage of amino acid types. To demonstrate the utility of this approach we applied it to the disordered actin-binding N-terminal domain of WASp interacting protein (WIP), a ubiquitous key modulator of cytoskeletal changes in a range of biological systems. One- and two-bond J(¹⁵N,¹³Cα) couplings were acquired for WIP residues 2-65 at various temperatures, and in denaturing and crowding environments. Under native conditions fitted J-couplings identified in the WIP conformational ensemble a propensity for extended conformation at residues 16-23 and 45-60, and a helical tendency at residues 28-42. These findings are consistent with a previous study of the based upon chemical shift and RDC data and confirm that the WIP^2-65 conformational ensemble is biased towards the structure assumed by this fragment in its actin-bound form. The effects of environmental changes upon this ensemble were readily apparent in the J-coupling data, which reflected a significant decrease in structural propensity at higher temperatures, in the presence of 8 M urea, and under the influence of a bacterial cell lysate. The latter suggests that crowding can cause protein unfolding through protein-protein interactions that stabilize the unfolded state. We conclude that J-couplings are a useful measureable in characterizing structural ensembles in IDPs, and that the proposed experiment provides a practical method for accurately performing such measurements, once again emphasizing the power of NMR in studying IDP behavior.

Protein Folding-How and Why: By Hydrogen Exchange, Fragment Separation, and Mass Spectrometry

Advanced hydrogen exchange (HX) methodology can now determine the structure of protein folding intermediates and their progression in folding pathways. Key developments over time include the HX pulse labeling method with nuclear magnetic resonance analysis, the fragment separation method, the addition to it of mass spectrometric (MS) analysis, and recent improvements in the HX MS technique and data analysis. Also, the discovery of protein foldons and their role supplies an essential interpretive link. Recent work using HX pulse labeling with MS analysis finds that a number of proteins fold by stepping through a reproducible sequence of native-like intermediates in an ordered pathway. The stepwise nature of the pathway is dictated by the cooperative foldon unit construction of the protein. The pathway order is determined by a sequential stabilization principle; prior native-like structure guides the formation of adjacent native-like structure. This view does not match the funneled energy landscape paradigm of a very large number of folding tracks, which was framed before foldons were known and is more appropriate for the unguided residue-level search to surmount an initial kinetic barrier rather than for the overall unfolded-state to native-state folding pathway.

Estimation of Hydrogen-Exchange Protection Factors from MD Simulation Based on Amide Hydrogen Bonding Analysis

Hydrogen exchange (HX) studies have provided critical insight into our understanding of protein folding, structure, and dynamics. More recently, hydrogen exchange mass spectrometry (HX-MS) has become a widely applicable tool for HX studies. The interpretation of the wealth of data generated by HX-MS experiments as well as other HX methods would greatly benefit from the availability of exchange predictions derived from structures or models for comparison with experiment. Most reported computational HX modeling studies have employed solvent-accessible-surface-area based metrics in attempts to interpret HX data on the basis of structures or models. In this study, a computational HX-MS prediction method based on classification of the amide hydrogen bonding modes mimicking the local unfolding model is demonstrated. Analysis of the NH bonding configurations from molecular dynamics (MD) simulation snapshots is used to determine partitioning over bonded and nonbonded NH states and is directly mapped into a protection factor (PF) using a logistics growth function. Predicted PFs are then used for calculating deuteration values of peptides and compared with experimental data. Hydrogen exchange MS data for fatty acid synthase thioesterase (FAS-TE) collected for a range of pHs and temperatures was used for detailed evaluation of the approach. High correlation between prediction and experiment for observable fragment peptides is observed in the FAS-TE and additional benchmarking systems that included various apo/holo proteins for which literature data were available. In addition, it is shown that HX modeling can improve experimental resolution through decomposition of in-exchange curves into rate classes, which correlate with prediction from MD. Successful rate class decompositions provide further evidence that the presented approach captures the underlying physical processes correctly at the single residue level. This assessment is further strengthened in a comparison of residue resolved protection factor predictions for staphylococcal nuclease with NMR data, which was also used to compare prediction performance with other algorithms described in the literature. The demonstrated transferable and scalable MD based HX prediction approach adds significantly to the available tools for HX-MS data interpretation based on available structures and models.

High-resolution HDX-MS reveals distinct mechanisms of RNA recognition and activation by RIG-I and MDA5

RIG-I and MDA5 are the major intracellular immune receptors that recognize viral RNA species and undergo a series of conformational transitions leading to the activation of the interferon-mediated antiviral response. However, to date, full-length RLRs have resisted crystallographic efforts and a molecular description of their activation pathways remains hypothetical. Here we employ hydrogen/deuterium exchange coupled with mass spectrometry (HDX-MS) to probe the apo states of RIG-I and MDA5 and to dissect the molecular details with respect to distinct RNA species recognition, ATP binding and hydrolysis and CARDs activation. We show that human RIG-I maintains an auto-inhibited resting state owing to the intra-molecular HEL2i-CARD2 interactions while apo MDA5 lacks the analogous intra-molecular interactions and therefore adopts an extended conformation. Our work demonstrates that RIG-I binds and responds differently to short triphosphorylated RNA and long duplex RNA and that sequential addition of RNA and ATP triggers specific allosteric effects leading to RIG-I CARDs activation. We also present a high-resolution protein surface mapping technique that refines the cooperative oligomerization model of neighboring MDA5 molecules on long duplex RNA. Taken together, our data provide a high-resolution view of RLR activation in solution and offer new evidence for the molecular mechanism of RLR activation.

Probing backbone dynamics with hydrogen/deuterium exchange mass spectrometry

Protein dynamics can be probed by the solution technique amide hydrogen/deuterium exchange. The exchange rate of hydrogen for deuterium along a peptide backbone is dependent on the extent of hydrogen bonding from secondary structure, accessibility by D2O, and protein motions. Both global and local conformational changes that alter bonding or structure will lead to changes in the amount of deuterium incorporated. The deuterium can be localized via pepsin digestion of the protein and quantified by electrospray ionization mass spectrometry through the mass shifts of the resulting peptides. The technique is emerging as an essential tool to study protein structure in solution due to the exceptional capability of examining both dynamic and structural changes related to protein function.

Affinity capture of biotinylated proteins at acidic conditions to facilitate hydrogen/deuterium exchange mass spectrometry analysis of multimeric protein complexes

Characterization of conformational and dynamic changes associated with protein interactions can be done by hydrogen/deuterium exchange mass spectrometry (HDX-MS) by comparing the deuterium uptake in the bound and unbound state of the proteins. Investigation of local hydrogen/deuterium exchange in heteromultimeric protein complexes poses a challenge for the method due to the increased complexity of the mixture of peptides originating from all interaction partners in the complex. Previously, interference of peptides from one interaction partner has been removed by immobilizing the intact protein on beads prior to the HDX-MS experiment. However, when studying protein complexes of more than two proteins, immobilization can possibly introduce steric limitations to the interactions. Here, we present a method based on the high affinity biotin-streptavidin interaction that allows selective capture of biotinylated proteins even under the extreme conditions for hydrogen/deuterium exchange quenching i.e. pH 2.5 and 0 °C. This biotin-streptavidin capture strategy allows hydrogen/deuterium exchange to occur in proteins in solution and enables characterization of specific proteins in heteromultimeric protein complexes without interference of peptides originating from other interaction partners in the complex. The biotin-streptavidin strategy has been successfully implemented in a model system with two recombinant monoclonal antibodies that target nonoverlapping epitopes on the human epidermal growth factor receptor (EGFR). We present a workflow for biotinylation and characterization of recombinant antibodies and demonstrate affinity capture of biotinylated antibodies under hydrogen/deuterium exchange quench conditions by the biotin-streptavidin strategy.

Mass spectrometry-based methods to study protein architecture and dynamics

Mass spectrometry is now an indispensable tool in the armamentarium of molecular biophysics, where it is used for tasks ranging from protein sequencing and mapping of post-translational modifications to studies of higher order structure, conformational dynamics, and interactions of proteins with small molecule ligands and other biopolymers. This mini-review highlights several popular mass spectrometry-based tools that are now commonly used for structural studies of proteins beyond their covalent structure with a particular emphasis on hydrogen exchange and direct electrospray ionization mass spectrometry.

Ligand-induced protein mobility in complexes of carbonic anhydrase II and benzenesulfonamides with oligoglycine chains

This paper describes a biophysical investigation of residual mobility in complexes of bovine carbonic anhydrase II (BCA) and para-substituted benzenesulfonamide ligands with chains of 1-5 glycine subunits, and explains the previously observed increase in entropy of binding with chain length. The reported results represent the first experimental demonstration that BCA is not the rigid, static globulin that has been typically assumed, but experiences structural fluctuations upon binding ligands. NMR studies with (15)N-labeled ligands demonstrated that the first glycine subunit of the chain binds without stabilization or destabilization by the more distal subunits, and suggested that the other glycine subunits of the chain behave similarly. These data suggest that a model based on ligand mobility in the complex cannot explain the thermodynamic data. Hydrogen/deuterium exchange studies provided a global estimate of protein mobility and revealed that the number of exchanged hydrogens of BCA was higher when the protein was bound to a ligand with five glycine subunits than when bound to a ligand with only one subunit, and suggested a trend of increasing number of exchanged hydrogens with increasing chain length of the BCA-bound ligand, across the series. These data support the idea that the glycine chain destabilizes the structure of BCA in a length-dependent manner, causing an increase in BCA mobility. This study highlights the need to consider ligand-induced mobility of even "static" proteins in studies of protein-ligand binding, including rational ligand design approaches.

Clusters of branched aliphatic side chains serve as cores of stability in the native state of the HisF TIM barrel protein

Imidazole-3-glycerol phosphate synthase is a heterodimeric allosteric enzyme that catalyzes consecutive reactions in imidazole biosynthesis through its HisF and HisH subunits. The unusually slow unfolding reaction of the isolated HisF TIM barrel domain from the thermophilic bacteria, Thermotoga maritima, enabled an NMR-based site-specific analysis of the main-chain hydrogen bonds that stabilize its native conformation. Very strong protection against exchange with solvent deuterium in the native state was found in a subset of buried positions in α-helices and pervasively in the underlying β-strands associated with a pair of large clusters of isoleucine, leucine and valine (ILV) side chains located in the α7(βα)8(βα)1-2 and α2(βα)3-6β7 segments of the (βα)8 barrel. The most densely packed region of the large cluster, α3(βα)4-6β7, correlates closely with the core of stability previously observed in computational, protein engineering and NMR dynamics studies, demonstrating a key role for this cluster in determining the thermodynamic and structural properties of the native state of HisF. When considered with the results of previous studies where ILV clusters were found to stabilize the hydrogen-bonded networks in folding intermediates for other TIM barrel proteins, it appears that clusters of branched aliphatic side chains can serve as cores of stability across the entire folding reaction coordinate of one of the most common motifs in biology.

Conformational characterization of the charge variants of a human IgG1 monoclonal antibody using H/D exchange mass spectrometry

MAb1, a human IgG1 monoclonal antibody produced in a NS0 cell line, exhibits charge heterogeneity because of the presence of variants formed by processes such as N-terminal glutamate cyclization, C-terminal lysine truncation, deamidation, aspartate isomerization and sialylation in the carbohydrate moiety. Four major charge variants of MAb1 were isolated and the conformations of these charge variants were studied using hydrogen/deuterium exchange mass spectrometry, including the H/D exchange time course (HX-MS) and the stability of unpurified proteins from rates of H/D exchange (SUPREX) techniques. HX-MS was used to evaluate the conformation and solution dynamics of MAb1 charge variants by measuring their deuterium buildup over time at the peptide level. The SUPREX technique evaluated the unfolding profile and relative stability of the charge variants by measuring the exchange properties of globally protected amide protons in the presence of a chemical denaturant. The H/D exchange profiles from both techniques were compared among the four charge variants of MAb1. The two techniques together offered extensive understanding about the local and subglobal/global unfolding of the charge variants of MAb1. Our results demonstrated that all four charge variants of MAb1 were not significantly different in conformation, solution dynamics and chemical denaturant-induced unfolding profile and stability, which aids in understanding the biofunctions of the molecules. The analytical strategy used for conformational characterization may also be applicable to comparability studies done for antibody therapeutics.

Remodeling of the folding free energy landscape of staphylococcal nuclease by cavity-creating mutations

The folding of staphylococcal nuclease (SNase) is known to proceed via a major intermediate in which the central OB subdomain is folded and the C-terminal helical subdomain is disordered. To identify the structural and energetic determinants of this folding free energy landscape, we have examined in detail, using high-pressure NMR, the consequences of cavity creating mutations in each of the two subdomains of an ultrastable SNase, Δ+PHS. The stabilizing mutations of Δ+PHS enhanced the population of the major folding intermediate. Cavity creation in two different regions of the Δ+PHS reference protein, despite equivalent effects on global stability, had very distinct consequences on the complexity of the folding free energy landscape. The L125A substitution in the C-terminal helix of Δ+PHS slightly suppressed the major intermediate and promoted an additional excited state involving disorder in the N-terminus, but otherwise decreased landscape heterogeneity with respect to the Δ+PHS background protein. The I92A substitution, located in the hydrophobic OB-fold core, had a much more profound effect, resulting in a significant increase in the number of intermediate states and implicating the entire protein structure. Denaturant (GuHCl) had very subtle and specific effects on the landscape, suppressing some states and favoring others, depending upon the mutational context. These results demonstrate that disrupting interactions in a region of the protein with highly cooperative, unfrustrated folding has very profound effects on the roughness of the folding landscape, whereas the effects are less pronounced for an energetically equivalent substitution in an already frustrated region.

Structural reorganization triggered by charging of Lys residues in the hydrophobic interior of a protein

Structural consequences of ionization of residues buried in the hydrophobic interior of proteins were examined systematically in 25 proteins with internal Lys residues. Crystal structures showed that the ionizable groups are buried. NMR spectroscopy showed that in 2 of 25 cases studied, the ionization of an internal Lys unfolded the protein globally. In five cases, the internal charge triggered localized changes in structure and dynamics, and in three cases, it promoted partial or local unfolding. Remarkably, in 15 proteins, the ionization of the internal Lys had no detectable structural consequences. Highly stable proteins appear to be inherently capable of withstanding the presence of charge in their hydrophobic interior, without the need for specialized structural adaptations. The extent of structural reorganization paralleled loosely with global thermodynamic stability, suggesting that structure-based pK(a) calculations for buried residues could be improved by calculation of thermodynamic stability and by enhanced conformational sampling.

Sulfamate proton solvent exchange in heparin oligosaccharides: evidence for a persistent hydrogen bond in the antithrombin-binding pentasaccharide Arixtra

Sulfamate groups (NHSO(3)(-)) are important structural elements in the glycosaminoglycans (GAGs) heparin and heparan sulfate (HS). In this work, proton nuclear magnetic resonance (NMR) line-shape analysis is used to explore the solvent exchange properties of the sulfamate NH groups within heparin-related mono-, di-, tetra- and pentasaccharides as a function of pH and temperature. The results of these experiments identified a persistent hydrogen bond within the Arixtra (fondaparinux sodium) pentasaccharide between the internal glucosamine sulfamate NH and the adjacent 3-O-sulfo group. This discovery provides new insights into the solution structure of the Arixtra pentasaccharide and suggests that 3-O-sulfation of the heparin N-sulfoglucosamine (GlcNS) residues pre-organize the secondary structure in a way that facilitates binding to antithrombin-III. NMR studies of the GlcNS NH groups can provide important information about heparin structure complementary to that available from NMR spectral analysis of the carbon-bound protons.

Localized hydration in lyophilized myoglobin by hydrogen-deuterium exchange mass spectrometry. 1. Exchange mapping

The local effects of hydration on myoglobin (Mb) in solid matrices containing mannitol or sucrose (1:1 w/w, protein:additive) were mapped using hydrogen-deuterium exchange with mass spectrometric analysis (HDX-MS) at 5 °C and compared to solution controls. Solid powders were exposed to D₂O(g) at controlled activity (a(w)) followed by reconstitution and analysis of the intact protein and peptides produced by pepsin digestion. HDX varied with matrix type, a(w), and position along the protein backbone. HDX was less in sucrose matrices than in mannitol matrices at all a(w) while the difference in solution was negligible. Differences in HDX in the two matrices were detectable despite similarities in their bulk water content. The extent of exchange in solids is proposed as a measure of the hydration of exchangeable amide groups, as well as protein conformation and dynamics; pepsin digestion allows these effects to be mapped with peptide-level resolution.

A delicate interplay of structure, dynamics, and thermodynamics for function: a high pressure NMR study of outer surface protein A

Outer surface protein A (OspA) is a crucial protein in the infection of Borrelia burgdorferi causing Lyme disease. We studied conformational fluctuations of OspA with high-pressure (15)N/(1)H two-dimensional NMR along with high-pressure fluorescence spectroscopy. We found evidence within folded, native OspA for rapid local fluctuations of the polypeptide backbone in the nonglobular single layer β-sheet connecting the N- and C-terminal domains with τ << ms, which may give the two domains certain independence in mobility and thermodynamic stability. Furthermore, we found that folded, native OspA is in equilibrium (τ >> ms) with a minor conformer I, which is almost fully disordered and hydrated for the entire C-terminal part of the polypeptide chain from β8 to the C-terminus. Conformer I is characterized with ΔG(0) = 32 ± 9 kJ/mol and ΔV(0) = -140 ± 40 mL/mol, populating only ∼0.001% at 40°C at 0.1 MPa, pH 5.9. Because in the folded conformer the receptor binding epitope of OspA is buried in the C-terminal domain, its transition into conformer I under in vivo conditions may be critical for the infection of B. burgdorferi. The formation and stability of the peculiar conformer I are apparently supported by a large packing defect or cavity located in the C-terminal domain.

Residue-specific side-chain packing determines the backbone dynamics of transmembrane model helices

The transmembrane domains (TMDs) of membrane-fusogenic proteins contain an overabundance of β-branched residues. In a previous effort to systematically study the relation among valine content, fusogenicity, and helix dynamics, we developed model TMDs that we termed LV-peptides. The content and position of valine in LV-peptides determine their fusogenicity and backbone dynamics, as shown experimentally. Here, we analyze their conformational dynamics and the underlying molecular forces using molecular-dynamics simulations. Our study reveals that backbone dynamics is correlated with the efficiency of side-chain to side-chain van der Waals packing between consecutive turns of the helix. Leu side chains rapidly interconvert between two rotameric states, thus favoring contacts to its i±3 and i±4 neighbors. Stereochemical restraints acting on valine side chains in the α-helix force both β-substituents into an orientation where i,i±3 interactions are less favorable than i,i±4 interactions, thus inducing a local packing deficiency at VV3 motifs. We provide a quantitative molecular model to explain the relationship among chain connectivity, side-chain mobility, and backbone flexibility. We expect that this mechanism also defines the backbone flexibility of natural TMDs.

Strategies for the thermodynamic characterization of linked binding/local folding reactions within the native state application to the LID domain of adenylate kinase from Escherichia coli

Conformational fluctuations in proteins have emerged as an important aspect of biological function, having been linked to processes ranging from molecular recognition and catalysis to allostery and signal transduction. In spite of the realization of their importance, however, the connections between fluctuations and function have largely been empirical, even when they have been quantitative. Part of the problem in understanding the role of fluctuations in function is the fact that the mere existence of fluctuations complicates the interpretation of classic mutagenesis approaches. Namely, mutagenesis, which is typically targeted to an internal position (to elicit an effect), will change the fluctuations as well as the structure of the native state. Decoupling these effects is essential to an unambiguous understanding of the role of fluctuations in function. Here, we use a mutation strategy that targets surface-exposed sites in flexible parts of the molecule for mutation to glycine. Such mutations leave the ground-state structure unaffected. As a result, we can assess the nature of the fluctuations, develop a quantitative model relating fluctuations to function (in this case, molecular recognition), and unambiguously resolve the probabilities of the fluctuating states. We show that when this approach is applied to Escherichia coli adenylate kinase (AK), unique thermodynamic and structural insights are obtained, even when classic mutagenesis approaches targeted to the same region yield ambiguous results.

Quantification of protein backbone hydrogen-deuterium exchange rates by solid state NMR spectroscopy

We present the quantification of backbone amide hydrogen-deuterium exchange rates (HDX) for immobilized proteins. The experiments make use of the deuterium isotope effect on the amide nitrogen chemical shift, as well as on proton dilution by deuteration. We find that backbone amides in the microcrystalline α-spectrin SH3 domain exchange rather slowly with the solvent (with exchange rates negligible within the individual (15)N-T (1) timescales). We observed chemical exchange for 6 residues with HDX exchange rates in the range from 0.2 to 5 s(-1). Backbone amide (15)N longitudinal relaxation times that we determined previously are not significantly affected for most residues, yielding no systematic artifacts upon quantification of backbone dynamics (Chevelkov et al. 2008b). Significant exchange was observed for the backbone amides of R21, S36 and K60, as well as for the sidechain amides of N38, N35 and for W41ε. These residues could not be fit in our previous motional analysis, demonstrating that amide proton chemical exchange needs to be considered in the analysis of protein dynamics in the solid-state, in case D(2)O is employed as a solvent for sample preparation. Due to the intrinsically long (15)N relaxation times in the solid-state, the approach proposed here can expand the range of accessible HDX rates in the intermediate regime that is not accessible so far with exchange quench and MEXICO type experiments.

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?

Neutralizing positive charges at the surface of a protein lowers its rate of amide hydrogen exchange without altering its structure or increasing its thermostability

This paper combines two techniques--mass spectrometry and protein charge ladders--to examine the relationship between the surface charge and hydrophobicity of a representative globular protein (bovine carbonic anhydrase II; BCA II) and its rate of amide hydrogen-deuterium (H/D) exchange. Mass spectrometric analysis indicated that the sequential acetylation of surface lysine-ε-NH3(+) groups--a type of modification that increases the net negative charge and hydrophobicity of the surface of BCA II without affecting its secondary or tertiary structure--resulted in a linear decrease in the aggregate rate of amide H/D exchange at pD 7.4, 15 °C. According to analysis with MS, the acetylation of each additional lysine generated between 1.4 and 0.9 additional hydrogens that are protected from H/D exchange during the 2 h exchange experiment at 15 °C, pD 7.4. NMR spectroscopy demonstrated that none of the hydrogen atoms which became protected upon acetylation were located on the side chain of the acetylated lysine residues (i.e., lys-ε-NHCOCH3) but were instead located on amide NHCO moieties in the backbone. The decrease in rate of exchange associated with acetylation paralleled a decrease in thermostability: the most slowly exchanging rungs of the charge ladder were the least thermostable (as measured by differential scanning calorimetry). This observation--that faster rates of exchange are associated with slower rates of denaturation--is contrary to the usual assumptions in protein chemistry. The fact that the rates of H/D exchange were similar for perbutyrated BCA II (e.g., [lys-ε-NHCO(CH2)2CH3]18) and peracetylated BCA II (e.g., [lys-ε-NHCOCH3]18) suggests that the electrostatic charge is more important than the hydrophobicity of surface groups in determining the rate of H/D exchange. These electrostatic effects on the kinetics of H/D exchange could complicate (or aid) the interpretation of experiments in which H/D exchange methods are used to probe the structural effects of non-isoelectric perturbations to proteins (i.e., phosphorylation, acetylation, or the binding of the protein to an oligonucleotide or to another charged ligand or protein).

Role of partial protein unfolding in alcohol-induced protein aggregation

Proteins aggregate in response to various stresses including changes in solvent conditions. Addition of alcohols has been recently shown to induce aggregation of disease-related as well as nondisease-related proteins. Here we probed the biophysical mechanisms underlying alcohol-induced protein aggregation, in particular the role of partial protein unfolding in aggregation. We have studied aggregation mechanisms due to benzyl alcohol which is used in numerous biochemical and biotechnological applications. We chose cytochrome c as a model protein, for the reason that various optical and structural probes are available to monitor its global and partial unfolding reactions. Benzyl alcohol induced the aggregation of cytochrome c in isothermal conditions and decreased the temperature at which the protein aggregates. However, benzyl alcohol did not perturb the overall native conformation of cytochrome c. Instead, it caused partial unfolding of a local protein region around the methionine residue at position 80. Site-specific optical probes, two-dimensional NMR titrations, and hydrogen exchange all support this conclusion. The protein aggregation temperature varied linearly with the melting temperature of the Met80 region. Stabilizing the Met80 region by heme iron reduction drastically decreased protein aggregation, which confirmed that the local unfolding of this region causes protein aggregation. These results indicate that a possible mechanism by which alcohols induce protein aggregation is through partial rather than complete unfolding of native proteins.