Hydrogen-deuterium Exchange NMR Research Articles

Protein Structure Prediction from NMR Hydrogen-Deuterium Exchange Data.

Amide hydrogen-deuterium exchange (HDX) has long been used to determine regional flexibility and binding sites in proteins; however, the data are too sparse for full structural characterization. Experiments that measure HDX rates, such as HDX-NMR, have far higher throughput compared to structure determination via X-ray crystallography, cryo-EM, or a full suite of NMR experiments. Data from HDX-NMR experiments encode information on the protein structure, making HDX a prime candidate to be supplemented by computational algorithms for protein structure prediction. We have developed a methodology to incorporate HDX-NMR data into ab initio protein structure prediction using the Rosetta software framework to predict structures based on experimental agreement. To demonstrate the efficacy of our algorithm, we examined 38 proteins with HDX-NMR data available, comparing the predicted model with and without the incorporation of HDX data into scoring. The root-mean-square deviation (rmsd, a measure of the average atomic distance between superimposed models) of the predicted model improved by 1.42 Å on average after incorporating the HDX-NMR data into scoring. The average rmsd improvement for the proteins where the selected model rmsd changed after incorporating HDX data was 3.63 Å, including one improvement of more than 11 Å and seven proteins improving by greater than 4 Å, with 12/15 proteins improving overall. Additionally, for independent verification, two proteins that were not part of the original benchmark were scored including HDX data, with a dramatic improvement of the selected model rmsd of nearly 9 Å for one of the proteins. Moreover, we have developed a confidence metric allowing us to successfully identify near-native models in the absence of a native structure. Improvement in model selection with a strong confidence measure demonstrates that protein structure prediction with HDX-NMR is a powerful tool which can be performed with minimal additional computational strain and expense.

Tryptophan-mediated Dimerization of the TssL Transmembrane Anchor Is Required for Type VI Secretion System Activity

The type VI secretion system (T6SS) is a multiprotein complex used by bacteria to deliver effectors into target cells. The T6SS comprises a bacteriophage-like contractile tail structure anchored to the cell envelope by a membrane complex constituted of the TssJ outer-membrane lipoprotein and the TssL and TssM inner-membrane proteins. TssJ establishes contact with the periplasmic domain of TssM whereas the transmembrane segments of TssM and its cytoplasmic domain interact with TssL. TssL protrudes in the cytoplasm but is anchored by a C-terminal transmembrane helix (TMH). Here, we show that TssL TMH dimerization is required for the stability of the protein and for T6SS function. Using the TOXCAT assay and point mutations of the 23 residues of the TssL TMH, we identified Thr194 and Trp199 as necessary for TssL TMH dimerization. NMR hydrogen–deuterium exchange experiments demonstrated the existence of a dimer with the presence of Trp185 and Trp199 at the interface. A structural model based on molecular dynamic simulations shows that TssL TMH dimer formation involves π–π interactions resulting from the packing of the two Trp199 rings at the C-terminus and of the six aromatic rings of Tyr184, Trp185 and Trp188 at the N-terminus of the TMH.

Formation of Annular Protofibrillar Assembly by Cysteine Tripeptide: Unraveling the Interactions with NMR, FTIR, and Molecular Dynamics.

Both hydrogen-bonding and hydrophobic interactions play a significant role in molecular assembly, including self-assembly of proteins and peptides. In this study, we report the formation of annular protofibrillar structure (diameter ∼500 nm) made of a newly synthesized s-benzyl-protected cysteine tripeptide, which was primarily stabilized by hydrogen-bonding and hydrophobic interactions. Atomic force microscopy and field emission scanning electron microscopy analyses found small oligomers (diameter ∼60 nm) to bigger annular (outer diameter ∼300 nm; inner diameter, 100 nm) and protofibrillar structures after 1-2 days of incubation. Rotating-frame Overhauser spectroscopic (ROESY) analysis revealed the presence of several nonbonded proton-proton interactions among the residues, such as amide protons with methylene group, aromatic protons with tertiary butyl group, and methylene protons with tertiary butyl group. These added significant stability to bring the peptides closer to form a well-ordered assembled structure. Hydrogen-deuterium exchange NMR measurement further suggested that two individual amide protons among the three amide groups were strongly engaged with the adjacent tripeptide via H-bond interaction. However, the remaining amide proton was found to be exposed to solvent and remained noninteracting with other tripeptide molecules. In addition to chemical shift values, a significant change in amide bond vibrations of the tripeptide was found due to the formation of the self-assembled structure. The amide I mode of vibrations involving two amide linkages appeared at 1641 and 1695 cm-1 in the solid state. However, in the assembled state, the stretching band at 1695 cm-1 became broad and slightly shifted to ∼1689 cm-1. On the contrary, the band at 1641 cm-1 shifted to 1659 cm-1 and indicated that the -C═O bond associated with this vibration became stronger in the assembled state. These changes in Fourier transform infrared spectroscopy frequency clearly indicated changes in the amide backbone conformation and the associated hydrogen-bonding pattern due to the formation of the assembled structure. In addition to hydrogen bonding, molecular dynamics simulation indicated that the number of π-π interactions also increased with increasing number of tripeptides participated in the self-assembly process. Combined results envisaged a cross β-sheet assembly unit consisting of four intermolecular hydrogen bonds. Such noncovalent peptide assemblies glued by hydrogen-bonding and other weak forces may be useful in developing nanocapsule and related materials.

Conformational Flexibility Differentiates Naturally Occurring Bet v 1 Isoforms.

The protein Bet v 1 represents the main cause for allergic reactions to birch pollen in Europe and North America. Structurally homologous isoforms of Bet v 1 can have different properties regarding allergic sensitization and Th2 polarization, most likely due to differential susceptibility to proteolytic cleavage. Using NMR relaxation experiments and molecular dynamics simulations, we demonstrate that the initial proteolytic cleavage sites in two naturally occurring Bet v 1 isoforms, Bet v 1.0101 (Bet v 1a) and Bet v 1.0102 (Bet v 1d), are conformationally flexible. Inaccessible cleavage sites in helices and strands are highly flexible on the microsecond-millisecond time scale, whereas those located in loops display faster nanosecond-microsecond flexibility. The data consistently show that Bet v 1.0102 is more flexible and conformationally heterogeneous than Bet v 1.0101. Moreover, NMR hydrogen-deuterium exchange measurements reveal that the backbone amides in Bet v 1.0102 are significantly more solvent exposed, in agreement with this isoform’s higher susceptibility to proteolytic cleavage. The differential conformational flexibility of Bet v 1 isoforms, along with the transient exposure of inaccessible sites to the protein surface, may be linked to proteolytic susceptibility, representing a potential structure-based rationale for the observed differences in Th2 polarization and allergic sensitization.

Quenched hydrogen-deuterium exchange NMR of a disease-relevant Aβ(1-42) amyloid polymorph.

Alzheimer’s disease is associated with the aggregation into amyloid fibrils of Aβ(1–42) and Aβ(1–40) peptides. Interestingly, these fibrils often do not obtain one single structure but rather show different morphologies, so-called polymorphs. Here, we compare quenched hydrogen-deuterium (H/D) exchange of a disease-relevant Aβ(1–42) fibril for which the 3D structure has been determined by solid-state NMR with H/D exchange previously determined on another structural polymorph. This comparison reveals secondary structural differences between the two polymorphs suggesting that the two polymorphisms can be classified as segmental polymorphs.

Protein knotting through concatenation significantly reduces folding stability.

Concatenation by covalent linkage of two protomers of an intertwined all-helical HP0242 homodimer from Helicobacter pylori results in the first example of an engineered knotted protein. While concatenation does not affect the native structure according to X-ray crystallography, the folding kinetics is substantially slower compared to the parent homodimer. Using NMR hydrogen-deuterium exchange analysis, we showed here that concatenation destabilises significantly the knotted structure in solution, with some regions close to the covalent linkage being destabilised by as much as 5 kcal mol−1. Structural mapping of chemical shift perturbations induced by concatenation revealed a pattern that is similar to the effect induced by concentrated chaotrophic agent. Our results suggested that the design strategy of protein knotting by concatenation may be thermodynamically unfavourable due to covalent constrains imposed on the flexible fraying ends of the template structure, leading to rugged free energy landscape with increased propensity to form off-pathway folding intermediates.

Intrinsic local destabilization of the C-terminus predisposes integrin α1 I domain to a conformational switch induced by collagen binding.

Integrin-collagen interactions play a critical role in a myriad of cellular functions that include immune response, and cell development and differentiation, yet their mechanism of binding is poorly understood. There is increasing evidence that conformational flexibility assumes a central role in the molecular mechanisms of protein-protein interactions and here we employ NMR hydrogen-deuterium exchange (HDX) experiments to explore the impact of slower timescale dynamic events. To gain insight into the mechanisms underlying collagen-induced conformational switches, we have undertaken a comparative study between the wild type integrin α1 I and a gain-of-function E317A mutant. NMR HDX results suggest a relationship between regions exhibiting a reduced local stability in the unbound I domain and those that undergo significant conformational changes upon binding. Specifically, the αC and α7 helices within the C-terminus are at the center of such major perturbations and present reduced local stabilities in the unbound state relative to other structural elements. Complementary isothermal titration calorimetry experiments have been performed to derive complete thermodynamic binding profiles for association of the collagen-like triple-helical peptide with wild type α1 I and E317A mutant. The differential energetics observed for E317A are consistent with the HDX experiments and support a model in which intrinsically destabilized regions predispose conformational rearrangement in the integrin I domain. This study highlights the importance of exploring different timescales to delineate allosteric and binding events.

Phosphorylation Interferes with Maturation of Amyloid-β Fibrillar Structure in the N Terminus

Neurodegeneration is characterized by the ubiquitous presence of modifications in protein deposits. Despite their potential significance in the initiation and progression of neurodegenerative diseases, the effects of posttranslational modifications on the molecular properties of protein aggregates are largely unknown. Here, we study the Alzheimer disease-related amyloid-β (Aβ) peptide and investigate how phosphorylation at serine 8 affects the structure of Aβ aggregates. Serine 8 is shown to be located in a region of high conformational flexibility in monomeric Aβ, which upon phosphorylation undergoes changes in local conformational dynamics. Using hydrogen-deuterium exchange NMR and fluorescence quenching techniques, we demonstrate that Aβ phosphorylation at serine 8 causes structural changes in the N-terminal region of Aβ aggregates in favor of less compact conformations. Structural changes induced by serine 8 phosphorylation can provide a mechanistic link between phosphorylation and other biological events that involve the N-terminal region of Aβ aggregates. Our data therefore support an important role of posttranslational modifications in the structural polymorphism of amyloid aggregates and their modulatory effect on neurodegeneration.

The Knotted Protein UCH-L1 Exhibits Partially Unfolded Forms under Native Conditions that Share Common Structural Features with Its Kinetic Folding Intermediates

The human ubiquitin C-terminal hydrolase, UCH-L1, is an abundant neuronal deubiquitinase that is associated with Parkinson’s disease. It contains a complex Gordian knot topology formed by the polypeptide chain alone. Using a combination of fluorescence-based kinetic measurements, we show that UCH-L1 has two distinct kinetic folding intermediates that are transiently populated on parallel pathways between the denatured and native states. NMR hydrogen-deuterium exchange (HDX) experiments indicate the presence of partially unfolded forms (PUFs) of UCH-L1 under native conditions. HDX measurements as a function of urea concentration were used to establish the structure of the PUFs and pulse-labelled HDX NMR was used to show that the PUFs and the folding intermediates are likely the same species. In both cases, a similar stable core encompassing most of the central β-sheet is highly structured and α-helix 3, which is partially formed, packs against it. In contrast to the stable β-sheet core, the peripheral α-helices display significant local fluctuations leading to rapid exchange. The results also suggest that the main difference between the two kinetic intermediates is structure and packing of α-helices 3 and 7 and the degree of structure in β-strand 5. Together, the fluorescence and NMR results establish that UCH-L1 neither folds through a continuum of pathways nor by a single discrete pathway. Its folding is complex, the β-sheet core forms early and is present in both intermediate states, and the rate-limiting step which is likely to involve the threading of the chain to form the 52-knot occurs late on the folding pathway.

Unraveling the folding mechanism of the smallest knotted protein, MJ0366.

Understanding the mechanism by which polypeptide chains thread themselves into topologically knotted structures has emerged to be a challenging subject not least because of the additional complexity associated with the spontaneous and efficient knotting and folding events. While recent theoretical calculations have made significant progress in establishing the atomistic folding pathways for a number of knotted proteins, experimental data on the folding stabilities and kinetic pathways of knotted proteins has been sparse. Using MJ0366 from Methanocaldococcus jannaschii, the smallest knotted protein known to date, as a model system, we set out to systematically investigate its folding equilibrium, kinetics, and internal dynamics under native and chemically denatured states. NMR hydrogen-deuterium exchange analysis indicates that the knotted region is the most stable structural element within the novel fold. Additionally, (15)N spin relaxation analysis reveals the presence of residual structures in urea-denatured MJ0366. Despite the apparent two-state equilibrium unfolding behavior during chemical denaturation, the kinetic unfolding pathway of MJ0366 involves the dissociation of the homodimeric native state into a native-like monomeric intermediate followed by unfolding into a denatured state. Our results provide comprehensive structural information regarding the folding dynamics and kinetic pathways of MJ0366, whose small size is ideal for converging experimental and theoretical findings to better understand the underlying principles of the folding of knotted proteins.

Effects of Ligand Binding on the Rigidity and Mobility of Proteins: An Experimental and Computational Approach

Cyclophilin A is an enzyme which plays a role in the of proteins. It also binds to and aids the function of the immunosuppressant drug cyclosporin A as well as binding to the HIV-1 capsid protein.We expand upon a method recently used in our study of HIV-1 protease and use the computational tools FIRST and FRODA to model the flexibility and motion of cyclophilin A in the presence and absence of its principal ligands.In order to verify the simulation process, we have also conducted hydrogen-deuterium exchange NMR (HDX) experiments, the results of which are also presented.The computational procedure is far quicker than more in-depth approaches such as molecular dynamics, yet it yields results which are comparable with experimental data, on a timescale of CPU-minutes.We use the computational tools FIRST and FRODA to model the flexibility and motion of cyclophilin A in order to predict the results of hydrogen-deuterium exchange NMR (HDX) experiments.The rigidity analysis software FIRST can be used to predict the folding cores of proteins identified as slowly exchanging residues in HDX. This prediction is improved using the protein mobility software FRODA.We are using these methods to investigate the effect of ligand binding on cyclophilin A computationally and experimentally.

Deciphering Protein Stability in Cells

Deciphering Protein Stability in Cells

Aβ Peptide Fibrillar Architectures Controlled by Conformational Constraints of the Monomer

Anomalous self-assembly of the Aβ peptide into fibrillar amyloid deposits is strongly correlated with the development of Alzheimer's disease. Aβ fibril extension follows a template guided “dock and lock” mechanism where polymerisation is catalysed by the fibrillar ends. Using surface plasmon resonance (SPR) and quenched hydrogen-deuterium exchange NMR (H/D-exchange NMR), we have analysed the fibrillar structure and polymerisation properties of both the highly aggregation prone Aβ1–40 Glu22Gly (Aβ40Arc) and wild type Aβ1–40 (Aβ40WT). The solvent protection patterns from H/D exchange experiments suggest very similar structures of the fibrillar forms. However, through cross-seeding experiments monitored by SPR, we found that the monomeric form of Aβ40WT is significantly impaired to acquire the fibrillar architecture of Aβ40Arc. A detailed characterisation demonstrated that Aβ40WT has a restricted ability to dock and isomerise with high binding affinity onto Aβ40Arc fibrils. These results have general implications for the process of fibril assembly, where the rate of polymerisation, and consequently the architecture of the formed fibrils, is restricted by conformational constraints of the monomers. Interestingly, we also found that the kinetic rate of fibril formation rather than the thermodynamically lowest energy state determines the overall fibrillar structure.

β-Hairpins with native-like and non-native hydrogen bonding patterns could form during the refolding of staphylococcal nuclease

Refolding of staphylococcal nuclease has been studied recently by hydrogen–deuterium exchange and NMR spectroscopy. These studies infer that β-hairpin formed by strand 2 and strand 3 connected by reverse turn forms early during the refolding of nuclease. Typically, hydrogen–deuterium exchange NMR techniques are usually carried out on a time scale of milliseconds whereas β-hairpins are known to fold on a much shorter time scale. It follows that in the experiments, the hydrogen–deuterium exchange protection patterns could be arising from a significant population of fully formed hairpins. In order to demonstrate it is the fully formed hairpins which gives rise to the hydrogen–deuterium exchange protection patterns, we have considered molecular dynamics simulation of the peptide 21DTVKLM YKGQPMTFR 35 from staphylococcal nuclease corresponding to the β-hairpin region, using GROMOS96 force field under NVT conditions. Starting from unfolded conformational states, the peptide folds into hairpin conformations with native-like and non-native hydrogen bonding patterns. Subsequent to folding, equilibrium conditions prevail. The computed protection factors and atom depth values, at equilibrium, of the various amide protons agree qualitatively with experimental observations. A collection of molecules following the trajectories observed in the simulations can account for experimental observations. These simulations provide a molecular picture of the formed hairpins and their conformational features during the refolding experiments on nuclease, monitored by hydrogen–deuterium exchange.

Conformational features of disulfide intact and reduced forms of hen egg white lysozyme in aqueous solution in the presence of trifluoroethanol (TFE): Implications for protein folding intermediates

The conformational features of reduced and disulfide intact hen egg white lysozyme in aqueous 2,2,2-trifluoroethanol (TFE) solutions have been examined using circular dichroism and fluorescence spectroscopies. We find that non-native like and molten globule like states are induced in reduced lysozyme, with non-native like and native like α helicities and with reduced accessibility of tryptophans to collisional quencher acrylamide, in the presence of TFE as judged from circular dichroism and fluorescence experiments. We correlate our results to kinetic hydrogen–deuterium exchange NMR results of the refolding of lysozyme available in the literature. We discuss the implications of our results for disulfide bond formation and protein folding intermediates.

Antibody-Protein Complexes

The X-ray crystallographic structures of five antibody-protein antigen complexes are briefly described. The three-dimensional structure of the antibody-antigen interface is compared to the results obtained by two other techniques: analysis of mutant binding data and NMR hydrogen-deuterium exchange.