The conformational landscape of a serpin N-terminal subdomain facilitates folding and in-cell quality control.

Many multi-domain proteins including the serpin family of serine protease inhibitors contain non-sequential domains composed of regions that are far apart in sequence. Because proteins are translated vectorially from N- to C-terminus, such domains pose a particular challenge: how to balance the conformational lability necessary to form productive interactions between early and late translated regions while avoiding aggregation. This balance is mediated by the protein sequence properties and the interactions of the folding protein with the cellular quality control machinery. For serpins, particularly α1-antitrypsin (AAT), mutations often lead to polymer accumulation in cells and consequent disease suggesting that the lability/aggregation balance is especially precarious. Therefore, we investigated the properties of progressively longer AAT N-terminal fragments in solution and in cells. The N-terminal subdomain, residues 1-190 (AAT190), is monomeric in solution and efficiently degraded in cells. More ß-rich fragments, 1-290 and 1-323, form small oligomers in solution, but are still efficiently degraded, and even the polymerization promoting Siiyama (S53F) mutation did not significantly affect fragment degradation. In vitro, the AAT190 region is among the last regions incorporated into the final structure. Hydrogen-deuterium exchange mass spectrometry and enhanced sampling molecular dynamics simulations show that AAT190 has a broad, dynamic conformational ensemble that helps protect one particularly aggregation prone ß-strand from solvent. These AAT190 dynamics result in transient exposure of sequences that are buried in folded, full-length AAT, which may provide important recognition sites for the cellular quality control machinery and facilitate degradation and, under favorable conditions, reduce the likelihood of polymerization.

Long-range allostery mediates the regulation of plasminogen activator inhibitor-1 by cell adhesion factor vitronectin.

The serpin plasminogen activator inhibitor 1 (PAI-1) spontaneously undergoes a massive structural change from a metastable and active conformation, with a solvent-accessible reactive center loop (RCL), to a stable, inactive, or latent conformation, with the RCL inserted into the central ß-sheet. Physiologically, conversion to the latent state is regulated by the binding of vitronectin, which hinders the latency transition rate approximately twofold. The molecular mechanisms leading to this rate change are unclear. Here, we investigated the effects of vitronectin on the PAI-1 latency transition using all-atom path sampling simulations in explicit solvent. In simulated latency transitions of free PAI-1, the RCL is quite mobile as is the gate, the region that impedes RCL access to the central ß-sheet. This mobility allows the formation of a transient salt bridge that facilitates the transition; this finding rationalizes existing mutagenesis results. Vitronectin binding reduces RCL and gate mobility by allosterically rigidifying structural elements over 40 Å away from the binding site, thus blocking transition to the latent conformation. The effects of vitronectin are propagated by a network of dynamically correlated residues including a number of conserved sites that were previously identified as important for PAI-1 stability. Simulations also revealed a transient pocket populated only in the vitronectin-bound state, corresponding to a cryptic drug-binding site identified by crystallography. Overall, these results shed new light on PAI-1 latency transition regulation by vitronectin and illustrate the potential of path sampling simulations for understanding functional protein conformational changes and for facilitating drug discovery.

Characterization of interaction between blood coagulation factor VIII and LRP1 suggests dynamic binding by alternating complex contacts.

BACKGROUND: Deficiency in blood coagulation factor VIII (FVIII) results in life-threating bleeding (hemophilia A) treated by infusions of FVIII concentrates. To improve disease treatment, FVIII has been modified to increase its plasma half-life, which requires understanding mechanisms of FVIII catabolism. An important catabolic actor is hepatic low density lipoprotein receptor-related protein 1 (LRP1), which also regulates many other clinically significant processes. Previous studies showed complexity of FVIII site for binding LRP1. OBJECTIVES: To characterize binding sites between FVIII and LRP1 and suggest a model of the interaction. METHODS: A series of recombinant ligand-binding complement-type repeat (CR) fragments of LRP1 including mutated variants was generated in a baculovirus system and tested for FVIII interaction using surface plasmon resonance, tissue culture model, hydrogen-deuterium exchange mass spectrometry, and in silico. RESULTS: Multiple CR doublets within LRP1 clusters II and IV were identified as alternative FVIII-binding sites. These interactions follow the canonical binding mode providing major binding energy, and additional weak interactions are contributed by adjacent CR domains. A representative CR doublet was shown to have multiple contact sites on FVIII. CONCLUSIONS: FVIII and LRP1 interact via formation of multiple complex contacts involving both canonical and non-canonical binding combinations. We propose that FVIII-LRP1 interaction occurs via switching such alternative binding combinations in a dynamic mode, and that this mechanism is relevant to other ligand interactions of the low-density lipoprotein receptor family members including LRP1.

Modeling the native ensemble of PhuS using enhanced sampling MD and HDX-ensemble reweighting.

The cytoplasmic heme binding protein from Pseudomonas aeruginosa, PhuS, plays two essential roles in regulating heme uptake and iron homeostasis. First, PhuS shuttles exogenous heme to heme oxygenase (HemO) for degradation and iron release. Second, PhuS binds DNA and modulates the transcription of the prrF/H small RNAs (sRNAs) involved in the iron-sparing response. Heme binding to PhuS regulates this dual function, as the unliganded form binds DNA, whereas the heme-bound form binds HemO. Crystallographic studies revealed nearly identical structures for apo- and holo-PhuS, and yet numerous solution-based measurements indicate that heme binding is accompanied by large conformational rearrangements. In particular, hydrogen-deuterium exchange mass spectrometry (HDX-MS) of apo- versus holo-PhuS revealed large differences in deuterium uptake, notably in α-helices 6, 7, and 8 (α6,7,8), which contribute to the heme binding pocket. These helices were mostly labile in apo-PhuS but largely protected in holo-PhuS. In contrast, in silico-predicted deuterium uptake levels of α6,7,8 from molecular dynamics (MD) simulations of the apo- and holo-PhuS structures are highly similar, consistent only with the holo-PhuS HDX-MS data. To rationalize this discrepancy between crystal structures, simulations, and observed HDX-MS, we exploit a recently developed computational approach (HDXer) that fits the relative weights of conformational populations within an ensemble of structures to conform to a target set of HDX-MS data. Here, a combination of enhanced sampling MD, HDXer, and dimensionality reduction analysis reveals an apo-PhuS conformational landscape in which α6, 7, and 8 are significantly rearranged compared to the crystal structure, including a loss of secondary structure in α6 and the displacement of α7 toward the HemO binding interface. Circular dichroism analysis confirms the loss of secondary structure, and the extracted ensembles of apo-PhuS and of heme-transfer-impaired H212R mutant, are consistent with known heme binding and transfer properties. The proposed conformational landscape provides structural insights into the modulation by heme of the dual function of PhuS.

Interpreting Hydrogen-Deuterium Exchange Experiments with Molecular Simulations: Tutorials and Applications of the HDXer Ensemble Reweighting Software [Article v1.0].

Hydrogen-deuterium exchange (HDX) is a comprehensive yet detailed probe of protein structure and dynamics and, coupled to mass spectrometry, has become a powerful tool for investigating an increasingly large array of systems. Computer simulations are often used to help rationalize experimental observations of exchange, but interpretations have frequently been limited to simple, subjective correlations between microscopic dynamical fluctuations and the observed macroscopic exchange behavior. With this in mind, we previously developed the HDX ensemble reweighting approach and associated software, HDXer, to aid the objective interpretation of HDX data using molecular simulations. HDXer has two main functions; first, to compute H-D exchange rates that describe each structure in a candidate ensemble of protein structures, for example from molecular simulations, and second, to objectively reweight the conformational populations present in a candidate ensemble to conform to experimental exchange data. In this article, we first describe the HDXer approach, theory, and implementation. We then guide users through a suite of tutorials that demonstrate the practical aspects of preparing experimental data, computing HDX levels from molecular simulations, and performing ensemble reweighting analyses. Finally we provide a practical discussion of the capabilities and limitations of the HDXer methods including recommendations for a user's own analyses. Overall, this article is intended to provide an up-to-date, pedagogical counterpart to the software, which is freely available at https://github.com/Lucy-Forrest-Lab/HDXer.

Structure and dynamics of an α-fucosidase reveal a mechanism for highly efficient IgG transfucosylation.

Fucosylation is important for the function of many proteins with biotechnical and medical applications. Alpha-fucosidases comprise a large enzyme family that recognizes fucosylated substrates with diverse α-linkages on these proteins. Lactobacillus casei produces an α-fucosidase, called AlfC, with specificity towards α(1,6)-fucose, the only linkage found in human N-glycan core fucosylation. AlfC and certain point mutants thereof have been used to add and remove fucose from monoclonal antibody N-glycans, with significant impacts on their effector functions. Despite the potential uses for AlfC, little is known about its mechanism. Here, we present crystal structures of AlfC, combined with mutational and kinetic analyses, hydrogen-deuterium exchange mass spectrometry, molecular dynamic simulations, and transfucosylation experiments to define the molecular mechanisms of the activities of AlfC and its transfucosidase mutants. Our results indicate that AlfC creates an aromatic subsite adjacent to the active site that specifically accommodates GlcNAc in α(1,6)-linkages, suggest that enzymatic activity is controlled by distinct open and closed conformations of an active-site loop, with certain mutations shifting the equilibrium towards open conformations to promote transfucosylation over hydrolysis, and provide a potentially generalizable framework for the rational creation of AlfC transfucosidase mutants.

Successes and challenges in simulating the folding of large proteins.

Computational simulations of protein folding can be used to interpret experimental folding results, to design new folding experiments, and to test the effects of mutations and small molecules on folding. However, whereas major experimental and computational progress has been made in understanding how small proteins fold, research on larger, multidomain proteins, which comprise the majority of proteins, is less advanced. Specifically, large proteins often fold via long-lived partially folded intermediates, whose structures, potentially toxic oligomerization, and interactions with cellular chaperones remain poorly understood. Molecular dynamics based folding simulations that rely on knowledge of the native structure can provide critical, detailed information on folding free energy landscapes, intermediates, and pathways. Further, increases in computational power and methodological advances have made folding simulations of large proteins practical and valuable. Here, using serpins that inhibit proteases as an example, we review native-centric methods for simulating the folding of large proteins. These synergistic approaches range from Go and related structure-based models that can predict the effects of the native structure on folding to all-atom-based methods that include side-chain chemistry and can predict how disease-associated mutations may impact folding. The application of these computational approaches to serpins and other large proteins highlights the successes and limitations of current computational methods and underscores how computational results can be used to inform experiments. These powerful simulation approaches in combination with experiments can provide unique insights into how large proteins fold and misfold, expanding our ability to predict and manipulate protein folding.

A temperature-dependent conformational shift in p38α MAPK substrate-binding region associated with changes in substrate phosphorylation profile.

Febrile-range hyperthermia worsens and hypothermia mitigates lung injury, and temperature dependence of lung injury is blunted by inhibitors of p38 mitogen-activated protein kinase (MAPK). Of the two predominant p38 isoforms, p38α is proinflammatory and p38ß is cytoprotective. Here, we analyzed the temperature dependence of p38 MAPK activation, substrate interaction, and tertiary structure. Incubating HeLa cells at 39.5 °C stimulated modest p38 activation, but did not alter tumor necrosis factor-α (TNFα)-induced p38 activation. In in vitro kinase assays containing activated p38α and MAPK-activated kinase-2 (MK2), MK2 phosphorylation was 14.5-fold greater at 39.5 °C than at 33 °C. By comparison, we observed only 3.1- and 1.9-fold differences for activating transcription factor-2 (ATF2) and signal transducer and activator of transcription-1α (STAT1α) and a 7.7-fold difference for p38ß phosphorylation of MK2. The temperature dependence of p38α:substrate binding affinity, as measured by surface plasmon resonance, paralleled substrate phosphorylation. Hydrogen-deuterium exchange MS (HDX-MS) of p38α performed at 33, 37, and 39.5 °C indicated temperature-dependent conformational changes in an α helix near the common docking and glutamate:aspartate substrate-binding domains at the known binding site for MK2. In contrast, HDX-MS analysis of p38ß did not detect significant temperature-dependent conformational changes in this region. We observed no conformational changes in the catalytic domain of either isoform and no corresponding temperature dependence in the C-terminal p38α-interacting region of MK2. Because MK2 participates in the pathogenesis of lung injury, the observed changes in the structure and function of proinflammatory p38α may contribute to the temperature dependence of acute lung injury.

Recommendations for performing, interpreting and reporting hydrogen deuterium exchange mass spectrometry (HDX-MS) experiments.

Hydrogen deuterium exchange mass spectrometry (HDX-MS) is a powerful biophysical technique being increasingly applied to a wide variety of problems. As the HDX-MS community continues to grow, adoption of best practices in data collection, analysis, presentation and interpretation will greatly enhance the accessibility of this technique to nonspecialists. Here we provide recommendations arising from community discussions emerging out of the first International Conference on Hydrogen-Exchange Mass Spectrometry (IC-HDX; 2017). It is meant to represent both a consensus viewpoint and an opportunity to stimulate further additions and refinements as the field advances.

Interlaboratory Comparison of Hydrogen-Deuterium Exchange Mass Spectrometry Measurements of the Fab Fragment of NISTmAb.

Hydrogen-deuterium exchange mass spectrometry (HDX-MS) is an established, powerful tool for investigating protein-ligand interactions, protein folding, and protein dynamics. However, HDX-MS is still an emergent tool for quality control of biopharmaceuticals and for establishing dynamic similarity between a biosimilar and an innovator therapeutic. Because industry will conduct quality control and similarity measurements over a product lifetime and in multiple locations, an understanding of HDX-MS reproducibility is critical. To determine the reproducibility of continuous-labeling, bottom-up HDX-MS measurements, the present interlaboratory comparison project evaluated deuterium uptake data from the Fab fragment of NISTmAb reference material (PDB: 5K8A ) from 15 laboratories. Laboratories reported ∼89 800 centroid measurements for 430 proteolytic peptide sequences of the Fab fragment (∼78 900 centroids), giving ∼100% coverage, and ∼10 900 centroid measurements for 77 peptide sequences of the Fc fragment. Nearly half of peptide sequences are unique to the reporting laboratory, and only two sequences are reported by all laboratories. The majority of the laboratories (87%) exhibited centroid mass laboratory repeatability precisions of ⟨ sLab⟩ ≤ (0.15 ± 0.01) Da (1σx̅). All laboratories achieved ⟨sLab⟩ ≤ 0.4 Da. For immersions of protein at THDX = (3.6 to 25) °C and for D2O exchange times of tHDX = (30 s to 4 h) the reproducibility of back-exchange corrected, deuterium uptake measurements for the 15 laboratories is σreproducibility15 Laboratories( tHDX) = (9.0 ± 0.9) % (1σ). A nine laboratory cohort that immersed samples at THDX = 25 °C exhibited reproducibility of σreproducibility25C cohort( tHDX) = (6.5 ± 0.6) % for back-exchange corrected, deuterium uptake measurements.

Molecular Basis of Broad Spectrum N-Glycan Specificity and Processing of Therapeutic IgG Monoclonal Antibodies by Endoglycosidase S2.

Immunoglobulin G (IgG) glycosylation critically modulates antibody effector functions. Streptococcus pyogenes secretes a unique endo-ß-N-acetylglucosaminidase, EndoS2, which deglycosylates the conserved N-linked glycan at Asn297 on IgG Fc to eliminate its effector functions and evade the immune system. EndoS2 and specific point mutants have been used to chemoenzymatically synthesize antibodies with customizable glycosylation for gain of functions. EndoS2 is useful in these schemes because it accommodates a broad range of N-glycans, including high-mannose, complex, and hybrid types; however, its mechanism of substrate recognition is poorly understood. We present crystal structures of EndoS2 alone and bound to complex and high-mannose glycans; the broad N-glycan specificity is governed by critical loops that shape the binding site of EndoS2. Furthermore, hydrolytic experiments, domain-swap chimeras, and hydrogen-deuterium exchange mass spectrometry reveal the importance of the carbohydrate-binding module in the mechanism of IgG recognition by EndoS2, providing insights into engineering enzymes to catalyze customizable glycosylation reactions.

Monitoring dendrimer conformational transition using 19 F and 1 H2 O NMR.

The conformational transition of a fluorinated amphiphilic dendrimer is monitored by the 1 H signal from water, alongside the 19 F signal from the dendrimer. High-field NMR data (chemical shift δ, self-diffusion coefficient D, longitudinal relaxation rate R1 , and transverse relaxation rate R2 ) for both dendrimer (19 F) and water (1 H) match each other in detecting the conformational transition. Among all parameters for both nuclei, the water proton transverse-relaxation rate R2 (1 H2 O) displays the highest relative scale of change upon conformational transition of the dendrimer. Hydrogen/deuterium-exchange mass spectrometry reveals that the compact form of the dendrimer has slower proton exchange with water than the extended form. This result suggests that the sensitivity of R2 (1 H2 O) toward dendrimer conformation originates, at least partially, from the difference in proton exchange efficiency between different dendrimer conformations. Finally, we also demonstrated that this conformational transition could be conveniently monitored using a low-field benchtop NMR spectrometer via R2 (1 H2 O). The 1 H2 O signal thus offers a simple way to monitor structural changes of macromolecules using benchtop time-domain NMR.

Conformational transition of a non-associative fluorinated amphiphile in aqueous solution. II. Conformational transition vs. supramolecular assembly.

Unlike many known amphiphiles, the fluorinated amphiphilic dendrimer studied in this work demonstrated a concentration-dependent conformational transition rather than micellization or assembly. Hydrophobic and hydrophilic interactions with water were suggested as the most probable driving force of this transition. This assumption was consistent with the observed 19F chemical shift changes of the dendrimer compared to a known micelle-forming fluorinated amphiphile. Since water is an important factor in the process, trends of the concentration-dependent changes in water proton transverse relaxation rate served as an indicator of structural changes and/or supramolecular assembly. The conformational transition process was also confirmed by ion-mobility mass-spectrometry. We suggested that structural features, namely, steric hindrances, prevented the micellization/assembly of the dendrimer of this study. This conclusion might inform the approach to develop novel unconventional amphiphiles.

The Helicobacter pylori adhesin protein HopQ exploits the dimer interface of human CEACAMs to facilitate translocation of the oncoprotein CagA.

Helicobacter pylori infects half of the world's population, and strains that encode the cag type IV secretion system for injection of the oncoprotein CagA into host gastric epithelial cells are associated with elevated levels of cancer. CagA translocation into host cells is dependent on interactions between the H. pylori adhesin protein HopQ and human CEACAMs. Here, we present high-resolution structures of several HopQ-CEACAM complexes and CEACAMs in their monomeric and dimeric forms establishing that HopQ uses a coupled folding and binding mechanism to engage the canonical CEACAM dimerization interface for CEACAM recognition. By combining mutagenesis with biophysical and functional analyses, we show that the modes of CEACAM recognition by HopQ and CEACAMs themselves are starkly different. Our data describe precise molecular mechanisms by which microbes exploit host CEACAMs for infection and enable future development of novel oncoprotein translocation inhibitors and H. pylori-specific antimicrobial agents.

All-Atom Simulations Reveal How Single-Point Mutations Promote Serpin Misfolding.

Protein misfolding is implicated in many diseases, including serpinopathies. For the canonical inhibitory serpin α1-antitrypsin, mutations can result in protein deficiencies leading to lung disease, and misfolded mutants can accumulate in hepatocytes, leading to liver disease. Using all-atom simulations based on the recently developed bias functional algorithm, we elucidate how wild-type α1-antitrypsin folds and how the disease-associated S (Glu264Val) and Z (Glu342Lys) mutations lead to misfolding. The deleterious Z mutation disrupts folding at an early stage, whereas the relatively benign S mutant shows late-stage minor misfolding. A number of suppressor mutations ameliorate the effects of the Z mutation, and simulations on these mutants help to elucidate the relative roles of steric clashes and electrostatic interactions in Z misfolding. These results demonstrate a striking correlation between atomistic events and disease severity and shine light on the mechanisms driving chains away from their correct folding routes.

Remodeling KRAS.

Mutations in members of the RAS family of small GTPases have been associated with numerous human cancers. However, RAS family members are notoriously difficult to target. In this issue of Structure, Lu et al. (2017) examine the effects of two compounds with distinct chemical scaffolds on the structure and dynamics of an oncogenic KRAS mutant, thus highlighting the usefulness of HDX-MS for drug development.

Structure and Dynamics of FosA-Mediated Fosfomycin Resistance in Klebsiella pneumoniae and Escherichia coli.

Fosfomycin exhibits broad-spectrum antibacterial activity and is being reevaluated for the treatment of extensively drug-resistant pathogens. Its activity in Gram-negative organisms, however, can be compromised by expression of FosA, a metal-dependent transferase that catalyzes the conjugation of glutathione to fosfomycin, rendering the antibiotic inactive. In this study, we solved the crystal structures of two of the most clinically relevant FosA enzymes: plasmid-encoded FosA3 from Escherichia coli and chromosomally encoded FosA from Klebsiella pneumoniae (FosAKP). The structure, molecular dynamics, catalytic activity, and fosfomycin resistance of FosA3 and FosAKP were also compared to those of FosA from Pseudomonas aeruginosa (FosAPA), for which prior crystal structures exist. E. coli TOP10 transformants expressing FosA3 and FosAKP conferred significantly greater fosfomycin resistance (MIC, >1,024 µg/ml) than those expressing FosAPA (MIC, 16 µg/ml), which could be explained in part by the higher catalytic efficiencies of the FosA3 and FosAKP enzymes. Interestingly, these differences in enzyme activity could not be attributed to structural differences at their active sites. Instead, molecular dynamics simulations and hydrogen-deuterium exchange experiments with FosAKP revealed dynamic interconnectivity between its active sites and a loop structure that extends from the active site of each monomer and traverses the dimer interface. This dimer interface loop is longer and more extended in FosAKP and FosA3 than in FosAPA, and kinetic analyses of FosAKP and FosAPA loop-swapped chimeric enzymes highlighted its importance in FosA activity. Collectively, these data yield novel insights into fosfomycin resistance that could be leveraged to develop new strategies to inhibit FosA and potentiate fosfomycin activity.

IL-1 Family Cytokines Use Distinct Molecular Mechanisms to Signal through Their Shared Co-receptor.

Within the interleukin 1 (IL-1) cytokine family, IL-1 receptor accessory protein (IL-1RAcP) is the co-receptor for eight receptor-cytokine pairs, including those involving cytokines IL-1ß and IL-33. Unlike IL-1ß, IL-33 does not have a signaling complex that includes both its cognate receptor, ST2, and the shared co-receptor IL-1RAcP, which we now present here. Although the IL-1ß and IL-33 complexes shared structural features and engaged identical molecular surfaces of IL-1RAcP, these cytokines had starkly different strategies for co-receptor engagement and signal activation. Our data suggest that IL-1ß binds to IL-1RI to properly present the cytokine to IL-1RAcP, whereas IL-33 binds to ST2 in order to conformationally constrain the cognate receptor in an IL-1RAcP-receptive state. These findings indicate that members of the IL-1 family of cytokines use distinct molecular mechanisms to signal through their shared co-receptor, and they provide the foundation from which to design new therapies to target IL-33 signaling.

Common coding variant in SERPINA1 increases the risk for large artery stroke.

Large artery atherosclerotic stroke (LAS) shows substantial heritability not explained by previous genome-wide association studies. Here, we explore the role of coding variation in LAS by analyzing variants on the HumanExome BeadChip in a total of 3,127 cases and 9,778 controls from Europe, Australia, and South Asia. We report on a nonsynonymous single-nucleotide variant in serpin family A member 1 (SERPINA1) encoding alpha-1 antitrypsin [AAT; p.V213A; P = 5.99E-9, odds ratio (OR) = 1.22] and confirm histone deacetylase 9 (HDAC9) as a major risk gene for LAS with an association in the 3'-UTR (rs2023938; P = 7.76E-7, OR = 1.28). Using quantitative microscale thermophoresis, we show that M1 (A213) exhibits an almost twofold lower dissociation constant with its primary target human neutrophil elastase (NE) in lipoprotein-containing plasma, but not in lipid-free plasma. Hydrogen/deuterium exchange combined with mass spectrometry further revealed a significant difference in the global flexibility of the two variants. The observed stronger interaction with lipoproteins in plasma and reduced global flexibility of the Val-213 variant most likely improve its local availability and reduce the extent of proteolytic inactivation by other proteases in atherosclerotic plaques. Our results indicate that the interplay between AAT, NE, and lipoprotein particles is modulated by the gate region around position 213 in AAT, far away from the unaltered reactive center loop (357-360). Collectively, our findings point to a functionally relevant balance between lipoproteins, proteases, and AAT in atherosclerosis.