Dynamics Rationalize Proteolytic Susceptibility of the Major Birch Pollen Allergen Bet v 1

To stay or not to stay intact as an allergen: the endolysosomal degradation assay used as tool to analyze protein immunogenicity and T cell epitopes

Antigen uptake and processing of exogenous proteins is critical for adaptive immunity, particularly for T helper cell activation. Proteins undergo distinct proteolytic processing in endolysosomal compartments of antigen-presenting cells. The resulting peptides are presented on MHC class II molecules and specifically recognized by T cells. The in vitro endolysosomal degradation assay mimics antigen processing by incubating a protein of interest with a protease cocktail derived from the endolysosomal compartments of antigen presenting cells. The kinetics of protein degradation is monitored by gel electrophoresis and allows calculation of a protein's half-life and thus endolysosomal stability. Processed peptides are analyzed by mass spectrometry and abundant peptide clusters are shown to harbor T cell epitopes. The endolysosomal degradation assay has been widely used to study allergens, which are IgE-binding proteins involved in type I hypersensitivity. In this review article, we provide the first comprehensive overview of the endolysosomal degradation of 29 isoallergens and variants originating from the PR-10, Ole e 1-like, pectate lyase, defensin polyproline-linked, non-specific lipid transfer, mite group 1, 2, and 5, and tropomyosin protein families. The assay method is described in detail and suggestions for improved standardization and reproducibility are provided. The current hypothesis implies that proteins with high endolysosomal stability can induce an efficient immune response, whereas highly unstable proteins are degraded early during antigen processing and therefore not efficient for MHC II peptide presentation. To validate this concept, systematic analyses of high and low allergenic representatives of protein families should be investigated. In addition to purified molecules, allergen extracts should be degraded to analyze potential matrix effects and gastrointestinal proteolysis of food allergens. In conclusion, individual protein susceptibility and peptides obtained from the endolysosomal degradation assay are powerful tools for understanding protein immunogenicity and T cell reactivity. Systematic studies and linkage with in vivo sensitization data will allow the establishment of (machine-learning) tools to aid prediction of immunogenicity and allergenicity. The orthogonal method could in the future be used for risk assessment of novel foods and in the generation of protein-based immunotherapeutics.

pH-dependent structural diversity of profilin allergens determines thermal stability

The family of profilin allergens is a common class of proteins found in plants, viruses and various eukaryotes including mammals. Profilins are characterized by an evolutionary conserved structural fold, which is responsible for their cross-reactive nature of Immunoglobulin E (IgE) antibodies. Despite their high overall structural similarity, they exhibit substantial differences in their biophysical properties, such as thermal and pH stability. To understand the origin of these functional differences of Amb a 8, Art v 4 and Bet v 2, we performed constant pH molecular dynamics simulation in combination with Gaussian accelerated MD simulations. Depending on the respective protonation at different pH levels, we find distinct differences in conformational flexibility, which are consistent with experimentally determined melting temperatures. These variations in flexibility are accompanied by ensemble shifts in the conformational landscape and quantified and localized by residue-wise B-factors and dihedral entropies. These findings strengthen the link between flexibility of profilin allergens and their thermal stability. Thus, our results clearly show the importance of considering protonation dependent conformational ensembles in solution to elucidate biophysical differences between these structurally similar allergens.

Structural Alterations of Antigens at the Material Interface: An Early Decision Toolbox Facilitating Safe-by-Design Nanovaccine Development

Nanomaterials have found extensive interest in the development of novel vaccines, as adjuvants and/or carriers in vaccination platforms. Conjugation of protein antigens at the particle surface by non-covalent adsorption is the most widely used approach in licensed particulate vaccines. Hence, it is essential to understand proteins' structural integrity at the material interface in order to develop safe-by-design nanovaccines. In this study, we utilized two model proteins, the wild-type allergen Bet v 1 and its hypoallergenic fold variant (BM4), to compare SiO2 nanoparticles with Alhydrogel® as particulate systems. A set of biophysical and functional assays including circular dichroism spectroscopy and proteolytic degradation was used to examine the antigens' structural integrity at the material interface. Conjugation of both biomolecules to the particulate systems decreased their proteolytic stability. However, we observed qualitative and quantitative differences in antigen processing concomitant with differences in their fold stability. These changes further led to an alteration in IgE epitope recognition. Here, we propose a toolbox of biophysical and functional in vitro assays for the suitability assessment of nanomaterials in the early stages of vaccine development. These tools will aid in safe-by-design innovations and allow fine-tuning the properties of nanoparticle candidates to shape a specific immune response.

Energy penalties enhance flexible receptor docking in a model cavity

Protein flexibility remains a major challenge in library docking because of difficulties in sampling conformational ensembles with accurate probabilities. Here, we use the model cavity site of T4 lysozyme L99A to test flexible receptor docking with energy penalties from molecular dynamics (MD) simulations. Crystallography with larger and smaller ligands indicates that this cavity can adopt three major conformations: open, intermediate, and closed. Since smaller ligands typically bind better to the cavity site, we anticipate an energy penalty for the cavity opening. To estimate its magnitude, we calculate conformational preferences from MD simulations. We find that including a penalty term is essential for retrospective ligand enrichment; otherwise, high-energy states dominate the docking. We then prospectively docked a library of over 900,000 compounds for new molecules binding to each conformational state. Absent a penalty term, the open conformation dominated the docking results; inclusion of this term led to a balanced sampling of ligands against each state. High ranked molecules were experimentally tested by Tm upshift and X-ray crystallography. From 33 selected molecules, we identified 18 ligands and determined 13 crystal structures. Most interesting were those bound to the open cavity, where the buried site opens to bulk solvent. Here, highly unusual ligands for this cavity had been predicted, including large ligands with polar tails; these were confirmed both by binding and by crystallography. In docking, incorporating protein flexibility with thermodynamic weightings may thus access new ligand chemotypes. The MD approach to accessing and, crucially, weighting such alternative states may find general applicability.

Polarizable and non-polarizable force fields: Protein folding, unfolding, and misfolding

Molecular dynamics simulations are an invaluable tool to characterize the dynamic motions of proteins in atomistic detail. However, the accuracy of models derived from simulations inevitably relies on the quality of the underlying force field. Here, we present an evaluation of current non-polarizable and polarizable force fields (AMBER ff14SB, CHARMM 36m, GROMOS 54A7, and Drude 2013) based on the long-standing biophysical challenge of protein folding. We quantify the thermodynamics and kinetics of the β-hairpin formation using Markov state models of the fast-folding mini-protein CLN025. Furthermore, we study the (partial) folding dynamics of two more complex systems, a villin headpiece variant and a WW domain. Surprisingly, the polarizable force field in our set, Drude 2013, consistently leads to destabilization of the native state, regardless of the secondary structure element present. All non-polarizable force fields, on the other hand, stably characterize the native state ensembles in most cases even when starting from a partially unfolded conformation. Focusing on CLN025, we find that the conformational space captured with AMBER ff14SB and CHARMM 36m is comparable, but the ensembles from CHARMM 36m simulations are clearly shifted toward disordered conformations. While the AMBER ff14SB ensemble overstabilizes the native fold, CHARMM 36m and GROMOS 54A7 ensembles both agree remarkably well with experimental state populations. In addition, GROMOS 54A7 also reproduces experimental folding times most accurately. Our results further indicate an over-stabilization of helical structures with AMBER ff14SB. Nevertheless, the presented investigations strongly imply that reliable (un)folding dynamics of small proteins can be captured in feasible computational time with current additive force fields.

pH-Induced Local Unfolding of the Phl p 6 Pollen Allergen From cpH-MD

Susceptibility to endosomal degradation is a decisive contribution to a protein's immunogenicity. It is assumed that the processing kinetics of structured proteins are inherently linked to their probability of local unfolding. In this study, we quantify the impact of endosomal acidification on the conformational stability of the major timothy grass pollen allergen Phl p 6. We use state of the art sampling approaches in combination with constant pH MD techniques to profile pH-dependent local unfolding events in atomistic detail. Integrating our findings into the current view on type 1 allergic sensitization, we characterize local protein dynamics in the context of proteolytic degradation at neutral and acidic pH for the wild type protein and point mutants with varying proteolytic stability. We analyze extensive simulation data using Markov state models and retrieve highly reliable thermodynamic and kinetic information at varying pH levels. Thereby we capture the impact of endolysosomal acidification on the structure and dynamics of the Phl p 6 mutants. We find that upon protonation at lower pH values, the conformational flexibilities in key areas of the wild type protein, i.e., T-cell epitopes and early proteolytic cleavage sites, increase significantly. A decrease of the pH even leads to local unfolding in otherwise stable secondary structure elements, which is a prerequisite for proteolytic cleavage. This effect is even more pronounced in the destabilized mutant, while no unfolding was observed for the stabilized mutant. In summary, we report detailed structural models which rationalize the experimentally observed cleavage pattern during endosomal acidification.