Hydrogen Exchange Protection Factors Research Articles

Dynamics and thermal stability of the bypass polymerase, DinB homolog (Dbh).

The DinB homolog polymerase (Dbh) is a member of the Y-family of translesion DNA polymerases that can synthesize using a damaged DNA template. Since Dbh comes from the thermophilic archaeon Sulfolobus acidocaldarius, it is capable of functioning over a wide range of temperatures. Existing X-ray structures were determined at temperatures where the protein is least active. Here we use NMR and circular dichroism to understand how the structure and dynamics of Dbh are affected by temperature (2°C-65°C) and metal ion binding in solution. We measured hydrogen exchange protection factors, temperature coefficients, and chemical shift perturbations with and without magnesium and manganese. We report on regions of the protein that become more dynamic as the temperature is increased toward the functional temperature. Hydrogen exchange protection factors and temperature coefficients reveal that both the thumb and finger domains are very dynamic relative to the palm and little-finger (LF) domains. These trends remain true at high temperature with dynamics increasing as temperatures increase from 35°C to 50°C. Notably, NMR spectra show that the Dbh tertiary structure cold denatures beginning at 25°C and increases in denaturation as the temperature is lowered to 5°C with little change observed by CD. Above 35°C, chemical shift perturbation analysis in the presence and absence of magnesium and manganese reveals three ion binding sites, without DNA bound. In contrast, these bound metals are not apparent in any Dbh crystal structures of the protein without DNA. Two ion binding sites are confirmed to be near the active site, as reported in other Y-family polymerases, and we report a novel ion binding site in the LF domain. Thus, the solution-state structure of the Dbh polymerase is distinct from that of the solid-state structures and shows an unusually high cold denaturation temperature.

Proteins' Evolution upon Point Mutations.

As the reader must be already aware, state-of-the-art protein folding prediction methods have reached a smashing success in their goal of accurately determining the three-dimensional structures of proteins. Yet, a solution to simple problems such as the effects of protein point mutations on their (i) native conformation; (ii) marginal stability; (iii) ensemble of high-energy nativelike conformations; and (iv) metamorphism propensity and, hence, their evolvability, remains as an unsolved problem. As a plausible solution to the latter, some properties of the amide hydrogen-deuterium exchange, a highly sensitive probe of the structure, stability, and folding of proteins, are assessed from a new perspective. The preliminary results indicate that the protein marginal stability change upon point mutations provides the necessary and sufficient information to estimate, through a Boltzmann factor, the evolution of the amide hydrogen exchange protection factors and, consequently, that of the ensemble of folded conformations coexisting with the native state. This work contributes to our general understanding of the effects of point mutations on proteins and may spur significant progress in our efforts to develop methods to determine the appearance of new folds and functions accurately.

Structural Ensemble of the Insulin Monomer.

Experimental evidence suggests that monomeric insulin exhibits significant conformational heterogeneity, and modifications of apparently disordered regions affect both biological activity and the longevity of pharmaceutical formulations, presumably through receptor binding and fibrillation/degradation, respectively. However, a microscopic understanding of conformational heterogeneity has been lacking. Here, we integrate all-atom molecular dynamics simulations with an analysis pipeline to investigate the structural ensemble of human insulin monomers. We find that 60% of the structures present at least one of the following elements of disorder: melting of the A-chain N-terminal helix, detachment of the B-chain N-terminus, and detachment of the B-chain C-terminus. We also observe partial melting and extension of the B-chain helix and significant conformational heterogeneity in the region containing the B-chain β-turn. We then estimate hydrogen-exchange protection factors for the sampled ensemble and find them in line with experimental results for KP-insulin, although the simulations underestimate the importance of unfolded states. Our results help explain the ready exchange of specific amide sites that appear to be protected in crystal structures. Finally, we discuss the implications for insulin function and stability.

Remarkable Rigidity of the Single α-Helical Domain of Myosin-VI As Revealed by NMR Spectroscopy.

Although the α-helix has long been recognized as an all-important element of secondary structure, it generally requires stabilization by tertiary interactions with other parts of a protein’s structure. Highly charged single α-helical (SAH) domains, consisting of a high percentage (>75%) of Arg, Lys, and Glu residues, are exceptions to this rule but have been difficult to characterize structurally. Our study focuses on the 68-residue medial tail domain of myosin-VI, which is found to contain a highly ordered α-helical structure extending from Glu-6 to Lys-63. High hydrogen exchange protection factors (15–150), small (ca. 4 Hz) 3JHNHα couplings, and a near-perfect fit to an ideal model α-helix for its residual dipolar couplings (RDCs), measured in a filamentous phage medium, support the high regularity of this helix. Remarkably, the hydrogen exchange rates are far more homogeneous than the protection factors derived from them, suggesting that for these transiently broken helices the intrinsic exchange rates derived from the amino acid sequence are not appropriate reference values. 15N relaxation data indicate a very high degree of rotational diffusion anisotropy (D∥/D⊥ ≈ 7.6), consistent with the hydrodynamic behavior predicted for such a long, nearly straight α-helix. Alignment of the helix by a paramagnetic lanthanide ion attached to its N-terminal region shows a decrease in alignment as the distance from the tagging site increases. This decrease yields a precise measure for the persistence length of 224 ± 10 Å at 20 °C, supporting the idea that the role of the SAH helix is to act as an extension of the myosin-VI lever arm.

Monitoring Hydrogen Exchange During Protein Folding by Fast Pressure Jump NMR Spectroscopy.

A method is introduced that permits direct observation of the rates at which backbone amide hydrogens become protected from solvent exchange after rapidly dropping the hydrostatic pressure inside the NMR sample cell from denaturing (2.5 kbar) to native (1 bar) conditions. The method is demonstrated for a pressure-sensitized ubiquitin variant that contains two Val to Ala mutations. Increased protection against hydrogen exchange with solvent is monitored as a function of time during the folding process. Results for 53 backbone amides show narrow clustering with protection occurring with a time constant of ca. 85 ms, but slower protection is observed around a reverse turn near the C-terminus of the protein. Remarkably, the native NMR spectrum returns with this slower time constant of ca. 150 ms, indicating that the almost fully folded protein retains molten globule characteristics with severe NMR line broadening until the final hydrogen bonds are formed. Prior to crossing the transition state barrier, hydrogen exchange protection factors are close to unity, but with slightly elevated values in the β1-β2 hairpin, previously shown to be already lowly populated in the urea-denatured state.

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.

Predicting the Predictive Power of IDP Ensembles

The function of intrinsically disordered proteins may be interpreted in terms of their structural ensembles. The article by Schwalbe and colleagues in this issue of Structure combines NMR and SAXS constraints to generate structural ensembles that unveil important functional and pathological features.

Long-Time Scale Fluctuations of Human Prion Protein Determined by Restrained MD Simulations

Cellular prion protein (PrP(C)) has the ability to trigger transmissible lethal diseases after in vivo maturation into a toxic amyloidogenic misfolded form (PrP(Sc)). Here, we use hydrogen exchange protection factors in restrained molecular dynamics simulations to characterize long-time scale fluctuations in human PrP(C). We find that the regions of residues 138-141 and 183-192 form new β-strands in several exchange-competent structures. Moreover, these structural changes are associated with the disruption of native contacts that when tethered prevent fibril formation. Our findings illustrate the structural plasticity of PrP(C) and are valuable for understanding the conversion of PrP(C) to PrP(Sc).

Prediction of Native-State Hydrogen Exchange from Perfectly Funneled Energy Landscapes

Simulations based on perfectly funneled energy landscapes often capture many of the kinetic features of protein folding. We examined whether simulations based on funneled energy functions can also describe fluctuations in native-state protein ensembles. We quantitatively compared the site-specific local stability determined from structure-based folding simulations, with hydrogen exchange protection factors measured experimentally for ubiquitin, chymotrypsin inhibitor 2, and staphylococcal nuclease. Different structural definitions for the open and closed states based on the number of native contacts for each residue, as well as the hydrogen-bonding state, or a combination of both criteria were evaluated. The predicted exchange patterns agree with the experiments under native conditions, indicating that protein topology indeed has a dominant effect on the exchange kinetics. Insights into the simplest mechanistic interpretation of the amide exchange process were thus obtained.

Structure and Disorder in an Unfolded State under Nondenaturing Conditions from Ensemble Models Consistent with a Large Number of Experimental Restraints

Obtaining detailed structural models of disordered states of proteins under nondenaturing conditions is important for a better understanding of both functional intrinsically disordered proteins and unfolded states of folded proteins. Extensive experimental characterization of the drk N-terminal SH3 domain unfolded state has shown that, although it appears to be highly disordered, it possesses significant nonrandom secondary and tertiary structure. In our previous attempts to generate structural models of the unfolded state using the program ENSEMBLE, we were limited by insufficient experimental restraints and conformational sampling. In this study, we have vastly expanded our experimental restraint set to include 1H– 15N residual dipolar couplings, small-angle X-ray scattering measurements, nitroxide paramagnetic relaxation enhancements, O 2-induced 13C paramagnetic shifts, hydrogen-exchange protection factors, and 15N R 2 data, in addition to the previously used nuclear Overhauser effects, amino terminal Cu 2+–Ni 2+ binding paramagnetic relaxation enhancements, J-couplings, chemical shifts, hydrodynamic radius, and solvent accessibility restraints. We have also implemented a new ensemble calculation methodology that uses iterative conformational sampling and seeks to calculate the simplest possible ensemble models. As a result, we can now generate ensembles that are consistent with much larger experimental data sets than was previously possible. Although highly heterogeneous and having broad molecular size distributions, the calculated drk N-terminal SH3 domain unfolded-state ensembles have very different properties than expected for random or statistical coils and possess significant nonnative α-helical structure and both native-like and nonnative tertiary structure.

Redox-Dependent Dynamics of a Dual Thioredoxin Fold Protein: Evolution of Specialized Folds

An enzyme system protecting bacteria from oxidative stress includes the flavoprotein AhpF and the peroxiredoxin AhpC. The N-terminal domain of AhpF (NTD), with two fused thioredoxin (Trx) folds, belongs to the hyperthermophilic protein disulfide oxidoreductase family. The NTD is distinct in that it contains a redox active a fold with a CxxC sequence and a redox inactive b fold that has lost the CxxC motif. Here we characterize the stability, the (15)N backbone relaxation, and the hydrogen-deuterium exchange properties of reduced [NTD-(SH)(2)] and oxidized (NTD-S(2)) NTD from Salmonella typhimurium. While both NTD-(SH)(2) and NTD-S(2) exhibit similar equilibrium unfolding transitions and order parameters, R(ex) relaxation terms are quite distinct with considerably more intermediate time scale motions in NTD-S(2). Hydrogen exchange protection factors show that the slowly exchanging core corresponds to residues in the b fold in both NTD-(SH)(2) and NTD-S(2). Interestingly, folded-state dynamic fluctuations in the catalytic a fold are significantly increased for residues in NTD-S(2) compared to NTD-(SH)(2). Taken together, these data demonstrate that oxidation of the active site disulfide does not significantly increase stability but results in a dramatic increase in conformational heterogeneity in residues primarily in the redox active a fold. Differences in dynamics between the two folds of the NTD suggest that each evolved a specialized function which, in the a fold, couples redox state to internal motions which may enhance catalysis and specificity and, in the b fold, provides a redox insensitive stable core.

Structural Interpretation of Hydrogen Exchange Protection Factors in Proteins: Characterization of the Native State Fluctuations of CI2

Protection factors obtained from equilibrium hydrogen exchange experiments are an important source of structural information on both native and nonnative states of proteins. We present a method for determining ensembles of protein structures by using hydrogen exchange data as restraints in molecular dynamics simulations in conjunction with an empirical force-field. The method is applied to determine the ensemble of structures representing the native state of chymotrypsin inhibitor 2 (CI2), including the rare, large fluctuations responsible for hydrogen exchange.

Determination of an ensemble of structures representing the intermediate state of the bacterial immunity protein Im7

We present a detailed structural characterization of the intermediate state populated during the folding and unfolding of the bacterial immunity protein Im7. We achieve this result by incorporating a variety of experimental data available for this species in molecular dynamics simulations. First, we define the structure of the exchange-competent intermediate state of Im7 by using equilibrium hydrogen-exchange protection factors. Second, we use this ensemble to predict Phi-values and compare the results with the experimentally determined Phi-values of the kinetic refolding intermediate. Third, we predict chemical-shift measurements and compare them with the measured chemical shifts of a mutational variant of Im7 for which the kinetic folding intermediate is the most stable state populated at equilibrium. Remarkably, we found that the properties of the latter two species are predicted with high accuracy from the exchange-competent intermediate that we determined, suggesting that these three states are characterized by a similar architecture in which helices I, II, and IV are aligned in a native-like, but reorganized, manner. Furthermore, the structural ensemble that we obtained enabled us to rationalize the results of tryptophan fluorescence experiments in the WT protein and a series of mutational variants. The results show that the integration of diverse sets of experimental data at relatively low structural resolution is a powerful approach that can provide insights into the structural organization of this conformationally heterogeneous three-helix intermediate with unprecedented detail and highlight the importance of both native and non-native interactions in stabilizing its structure.

Slow Protein Conformational Dynamics from Multiple Experimental Structures: The Helix/Sheet Transition of Arc Repressor

Conformational transitions underlie the function of many biomolecular systems. Resolving intermediate structural changes, however, is challenging for both experiments and all-atom simulations because the duration of transitions is short relative to the lifetime of the stable species. Simplified descriptions based on a single experimental structure, such as elastic network models or Gō models, are not immediately applicable. Here, we develop a general method that combines multiple coarse-grained models to capture slow conformational transitions. Individually, each model describes one of the experimental structures; together, they approximate the complete energy surface. We demonstrate the method for the helix-to-sheet transition in Arc repressor N11L. We find that the transition involves the partial unfolding of the switch region, and rapid refolding into the alternate structure. Transient local unfolding is consistent with the low hydrogen exchange protection factors of the switch region. Also in agreement with experiment, the isomerization occurs independently of the global folding/dimerization transition.

Contribution of the multi-turn segment in the reversible thermal stability of hyperthermophile rubredoxin: NMR thermal chemical exchange analysis of sequence hybrids

Pyrococcus furiosus ( Pf) rubredoxin is the most thermostable protein characterized to date. Reflecting the complications arising from irreversible denaturation of this protein, predictions of which structural regions confer differential thermal stability have utilized kinetic stability measurements, hydrogen exchange protection factors, long range hydrogen bond NMR spin couplings, and molecular dynamics simulations, and have primarily implicated the three-stranded β-sheet and the adjacent metal binding site. Herein, NMR chemical exchange experiments demonstrate reversible two-state unfolding at the thermal transition temperature ( T m) for hybrids of Pf and the mesophile Clostridium pasteurianum ( Cp) rubredoxins which interchange residues 14–33, the so-called multi-turn segment. This complementary pair of hybrid rubredoxins exhibits largely additive incremental thermal stabilizations vs. the parental proteins. Both stabilization free energy measurements as well as incremental T m values indicate that a minimum of 37% of the total differential thermal stability resides in this multi-turn segment. Such a proportionality between ΔΔ G and incremental T m values is predicted for hybrid pairs exhibiting thermodynamic additivity in which the differential stability is predominantly enthalpic.

Identification of Novel Interactions in HIV-1 Capsid Protein Assembly by High-resolution Mass Spectrometry

The pleomorphic nature of the immature and mature HIV-1 virions has made it difficult to characterize intersubunit interactions using traditional approaches. While the structures of isolated domains are known, the challenge is to identify intersubunit interactions and thereby pack these domains into supramolecular structures. Using high-resolution mass spectrometry, we have measured the amide hydrogen exchange protection factors for the soluble capsid protein (CA) and CA assembled in vitro. Comparison of the protection factors as well as chemical crosslinking experiments has led to a map of the subunit/subunit interfaces in the assembled tubes. This analysis provides direct biochemical evidence for the homotypic N domain and C domain interactions proposed from cryo-electron microscopy image reconstruction of CA tubes. Most significantly, we have identified a previously unrecognized intersubunit N domain–C domain interaction. The detection of this interaction reconciles previously discrepant biophysical and genetic data.

Accessing the global minimum conformation of stefin A dimer by annealing under partially denaturing conditions

Stefin A folds as a monomer under strongly native conditions. We have observed that under partially denaturing conditions in the temperature range from 74 to 93 °C it folds into a dimer, while it is monomeric above the melting temperature of 95 °C. Below 74 °C the dimer is trapped and it does not dissociate. The dimer is a folded and structured protein as judged by CD and NMR, nevertheless it is no more functional as an inhibitor of cysteine proteases. The monomer-dimer transition proceeds at a slow rate and the activation energy of dimerization at 99 kcal/mol is comparable to the unfolding enthalpy. A large and negative dimerization enthalpy of −111(±8) kcal/mol was calculated from the temperature dependence of the dissociation constant. An irreversible pretransition at 10–15 deg. below the global unfolding temperature has been observed previously by DSC and can now be assigned to the monomer-dimer transition. Backbone resonances of all the dimer residues were assigned using 15N isotopically enriched protein. The dimer is symmetric and the chemical shift differences between the monomer and dimer are localized around the tripartite hydrophobic wedge, which otherwise interacts with cysteine proteases. Hydrogen exchange protection factors of the residues affected by dimer formation are higher in the dimer than in the monomer. The monomer to dimer transition is accompanied by a rapid exchange of all of the amide protons which are protected in the dimer, indicating that the transition state is unfolded to a large extent. Our results demonstrate that the native monomeric state of stefin A is actually metastable but is favored by the kinetics of folding. The substantial energy barrier which separates the monomer from the more stable dimer traps each state under native conditions.

Identification of cooperative folding units in a set of native proteins.

Cooperative unfolding penalties are calculated by statistically evaluating an ensemble of denatured states derived from native structures. The ensemble of denatured states is determined by dividing the native protein into short contiguous segments and defining all possible combinations of native, i.e., interacting, and non-native, i.e., non-interacting, segments. We use a novel knowledge-based scoring function, derived from a set of non-homologous proteins in the Protein Data Bank, to describe the interactions among residues. This procedure is used for the structural identification of cooperative folding cores for four globular proteins: bovine pancreatic trypsin inhibitor, horse heart cytochrome c, French bean plastocyanin, and staphylococcal nuclease. The theoretical folding units are shown to correspond to regions that exhibit enhanced stability against denaturation as determined from experimental hydrogen exchange protection factors. Using a sequence similarity score for related sequences, we show that, in addition to residues necessary for enzymatic function, those amino acids comprising structurally important folding cores are also preferentially conserved during evolution. This implies that the identified folding cores may be part of an array of fundamental structural folding units.

Structural characterisation and comparison of the native and A-states of equine lysozyme

Native state 1H NMR resonance assignments for 125 of the 129 residues of equine lysozyme have enabled measurement of the hydrogen exchange kinetics for over 60 backbone amide and three tryptophan indole hydrogen atoms in the native state. Native holo equine lysozyme hydrogen exchange protection factors are as large as 10 6, the most protected residues being located in elements of secondary structure. High exchange protection in the domain interface correlates with the binding of Ca 2+ in this region. Equine lysozyme differs from most non-Ca 2+ binding lysozymes in forming a highly populated partially folded state at low pH. The protein in this A-state at pH 2.0 has been found to bind 1-anilino-naphthalene-8-sulphonate with the enhancement of fluorescent intensity and blue shift in the spectral maximum characteristic of molten globules. NMR spectra indicate that the A-state is globally much less ordered than native equine lysozyme but does not contain significant regions of random coil structure. The amides most protected against hydrogen exchange in the A-state (protection factors up to 10 2 at 5°C) correspond to residues of three of the four α-helices of the native state; the side-chains of these residues form a hydrophobic cluster that includes five aromatic residues. Circular dichroism and tryptophan fluorescence indicate that these residues are substantially more constrained than similar residues in “classical” molten globules. Taken together, the data suggest a model for the A-state of equine lysozyme in which a more ordered core is surrounded by a less ordered but still compact polypeptide chain.

Predicting the equilibrium protein folding pathway: structure-based analysis of staphylococcal nuclease.

The equilibrium folding pathway of staphylococcal nuclease (SNase) has been approximated using a statistical thermodynamic formalism that utilizes the high-resolution structure of the native state as a template to generate a large ensemble of partially folded states. Close to 400,000 different states ranging from the native to the completely unfolded states were included in the analysis. The probability of each state was estimated using an empirical structural parametrization of the folding energetics. It is shown that this formalism predicts accurately the stability of the protein, the cooperativity of the folding/unfolding transition observed by differential scanning calorimetry (DSC) or urea denaturation and the thermodynamic parameters for unfolding. More importantly, this formalism provides a quantitative account of the experimental hydrogen exchange protection factors measured under native conditions for SNase. These results suggest that the computer-generated distribution of states approximates well the ensemble of conformations existing in solution. Furthermore, this formalism represents the first model capable of quantitatively predicting within a unified framework the probability distribution of states seen under native conditions and its change upon unfolding.