Identification and semi-quantification of protein allergens in complex mixtures using proteomic and AllerCatPro 2.0 bioinformatic analyses: a proof-of-concept investigation.

The demand for botanicals and natural substances in consumer products has increased in recent years. These substances usually contain proteins and these, in turn, can pose a risk for immunoglobulin E (IgE)-mediated sensitization and allergy. However, no method has yet been accepted or validated for assessment of potential allergenic hazards in such materials. In the studies here, a dual proteomic-bioinformatic approach is proposed to evaluate holistically allergenic hazards in complex mixtures of plants, insects, or animal proteins. Twelve commercial preparations of source materials (plant products, dust mite extract, and preparations of animal dander) known to contain allergenic proteins were analyzed by label-free proteomic analyses to identify and semi-quantify proteins. These were then evaluated by bioinformatics using AllerCatPro 2.0 (https://allercatpro.bii.a-star.edu.sg/) to predict no, weak, or strong evidence for allergenicity and similarity to source-specific allergens. In total, 4,586 protein sequences were identified in the 12 source materials combined. Of these, 1,665 sequences were predicted with weak or strong evidence for allergenic potential. This first-tier approach provided top-level information about the occurrence and abundance of proteins and potential allergens. With regards to source-specific allergens, 129 allergens were identified. The sum of the relative abundance of these allergens ranged from 0.8% (lamb's quarters) to 63% (olive pollen). It is proposed here that this dual proteomic-bioinformatic approach has the potential to provide detailed information on the presence and relative abundance of allergens, and can play an important role in identifying potential allergenic hazards in complex protein mixtures for the purposes of safety assessments.

Application of AllerCatPro 2.0 for protein safety assessments of consumer products.

Foreign proteins are potentially immunogenic, and a proportion of these are able to induce immune responses that result in allergic sensitization. Subsequent exposure of sensitized subjects to the inducing protein can provoke a variety of allergic reactions that may be severe, or even fatal. It has therefore been recognized for some time that it is important to determine a priori whether a given protein has the potential to induce allergic responses in exposed subjects. For example, the need to assess whether transgene products expressed in genetically engineered crop plants have allergenic properties. This is not necessarily a straightforward exercise (as discussed elsewhere in this edition), but the task becomes even more challenging when there is a need to conduct an overall allergenicity safety assessment of complex mixtures of proteins in botanicals or other natural sources that are to be used in consumer products. This paper describes a new paradigm for the allergenicity safety assessment of proteins that is based on the use of AllerCatPro 2.0, a new version of a previously described web application model developed for the characterization of the allergenic potential of proteins. Operational aspects of AllerCatPro 2.0 are described with emphasis on the application of new features that provide improvements in the predictions of allergenic properties such as the identification of proteins with high allergenic concern. Furthermore, the paper provides a description of strategies of how AllerCatPro 2.0 can best be deployed as a screening tool for identifying suitable proteins as ingredients in consumer products as well as a tool, in conjunction with label-free proteomic analysis, for identifying and semiquantifying protein allergens in complex materials. Lastly, the paper discusses the steps that are recommended for formal allergenicity safety assessment of novel consumer products which contain proteins, including consideration and integration of predicted consumer exposure metrics. The article therefore provides a holistic perspective of the processes through which effective protein safety assessments can be made of potential allergenic hazards and risks associated with exposure to proteins in consumer products, with a particular focus on the use of AllerCatPro 2.0 for this purpose.

Comparison of protein extraction protocols and allergen mapping from black soldier fly Hermetia illucens.

Exploration of important insect proteins - including allergens - and proteomes can be limited by protein extraction buffer selection and the complexity of the proteome. Herein, LC-MS/MS-based proteomics experiments were used to assess the protein extraction efficiencies for a suite of extraction buffers and the effect of ingredient processing on proteome and allergen detection. Discovery proteomics revealed that SDS-based buffer yields the maximum number of protein groups from three types of BSF samples. Bioinformatic analysis revealed that buffer composition and ingredient processing could influence allergen detection. Upon applying multi-level filtering criteria, 33 putative allergens were detected by comparing the detected BSF proteins to sequences from public allergen protein databases. A targeted LC-MRM-MS assay was developed for the pan-allergen tropomyosin and used to assess the influence of buffer composition and ingredient processing using peptide abundance measurements. SIGNIFICANCE: We demonstrated that the selection of protein extraction buffer and the processing method could influence protein yield and cross-reactive allergen detection from processed and un-processed black soldier fly (BSF) samples. In total, 33 putative allergens were detected by comparing the detected BSF proteins to sequences from public allergen protein databases. An LC-MRM-MS assay was developed for tropomyosin, indicating the importance of buffer selection and processing conditions to reduce BSF samples' allergenicity.

Corrigendum: A Web Server for GPCR-GPCR Interaction Pair Prediction.

[This corrects the article DOI: 10.3389/fendo.2022.825195.].

AllerCatPro 2.0: a web server for predicting protein allergenicity potential.

Proteins in food and personal care products can pose a risk for an immediate immunoglobulin E (IgE)-mediated allergic response. Bioinformatic tools can assist to predict and investigate the allergenic potential of proteins. Here we present AllerCatPro 2.0, a web server that can be used to predict protein allergenicity potential with better accuracy than other computational methods and new features that help assessors making informed decisions. AllerCatPro 2.0 predicts the similarity between input proteins using both their amino acid sequences and predicted 3D structures towards the most comprehensive datasets of reliable proteins associated with allergenicity. These datasets currently include 4979 protein allergens, 162 low allergenic proteins, and 165 autoimmune allergens with manual expert curation from the databases of WHO/International Union of Immunological Societies (IUIS), Comprehensive Protein Allergen Resource (COMPARE), Food Allergy Research and Resource Program (FARRP), UniProtKB and Allergome. Various examples of profilins, autoimmune allergens, low allergenic proteins, very large proteins, and nucleotide input sequences showcase the utility of AllerCatPro 2.0 for predicting protein allergenicity potential. The AllerCatPro 2.0 web server is freely accessible at https://allercatpro.bii.a-star.edu.sg.

A Web Server for GPCR-GPCR Interaction Pair Prediction.

The GGIP web server (https://protein.b.dendai.ac.jp/GGIP/) provides a web application for GPCR-GPCR interaction pair prediction by a support vector machine. The server accepts two sequences in the FASTA format. It responds with a prediction that the input GPCR sequence pair either interacts or not. GPCRs predicted to interact with the monomers constituting the pair are also shown when query sequences are human GPCRs. The server is simple to use. A pair of amino acid sequences in the FASTA format is pasted into the text area, a PDB ID for a template structure is selected, and then the 'Execute' button is clicked. The server quickly responds with a prediction result. The major advantage of this server is that it employs the GGIP software, which is presently the only method for predicting GPCR-interaction pairs. Our web server is freely available with no login requirement. In this article, we introduce some application examples of GGIP for disease-associated mutation analysis.

Global spectrum of population-specific common missense variation in cytochrome P450 pharmacogenes.

Next-generation sequencing technology has afforded the discovery of many novel variants that are of significance to inheritable pharmacogenomics (PGx) traits but a large proportion of them have unknown consequences. These include missense variants resulting in single amino acid substitutions in cytochrome P450 (CYP) proteins that can impair enzyme function, leading to altered drug efficacy and toxicity. While most unknown variants are rare, an overlooked minority are variants that are collectively rare but enriched in specific populations. Here, we analyzed sequence variation data in 141,456 individuals from across eight study populations in gnomAD for 38 CYP genes to identify such variants in addition to common variants. By further comparison with data from two PGx-specific databases (PharmVar and PharmGKB) and ClinVar, we identified 234 missense variants in 35 CYP genes, of which 107 were unknown to these databases. Most unknown variants (n = 83) were population-specific common variants and several (n = 7) were found in important CYP pharmacogenes (CYP2D6, CYP4F2, and CYP2C19). Overall, 29% (n = 31) of 107 unknown variants were predicted to affect CYP enzyme function although further biochemical characterization is necessary. These variants may elucidate part of the unexplained interpopulation differences observed in drug response.

Protein extraction protocols for optimal proteome measurement and arginine kinase quantitation from cricket Acheta domesticus for food safety assessment.

Insects have been consumed by people for millennia and have recently been proposed as a complementary, sustainable source of protein to feed the world's growing population. Insects and crustaceans both belong to the arthropod family. Crustacean (shellfish) allergies are common and potentially severe; hence, the cross-reactivity of the immune system with insect proteins is a potential health concern. Herein, LC-MS/MS was used to explore the proteome of whole, roasted whole and roasted powdered cricket products. Eight protein extraction protocols were compared using the total number of protein and distinct peptide identifications. Within these data, 20 putative allergens were identified, of which three were arginine kinase (AK) proteoforms. Subsequently, a multiple reaction monitoring MS assay was developed for the AK proteoforms and applied to a subset of extracts. This targeted assay demonstrated that allergen abundance/detectability varies according to the extraction method as well as the food processing method.

SVEP1 as a Genetic Modifier of TEK-Related Primary Congenital Glaucoma.

Purpose: Affecting children by age 3, primary congenital glaucoma (PCG) can cause debilitating vision loss by the developmental impairment of aqueous drainage resulting in high intraocular pressure (IOP), globe enlargement, and optic neuropathy. TEK haploinsufficiency accounts for 5% of PCG in diverse populations, with low penetrance explained by variable dysgenesis of Schlemm's canal (SC) in mice. We report eight families with TEK-related PCG, and provide evidence for SVEP1 as a disease modifier in family 8 with a higher penetrance and severity. Methods: Exome sequencing identified coding/splice site variants with an allele frequency less than 0.0001 (gnomAD). TEK variant effects were assayed in construct-transfected HEK293 cells via detection of autophosphorylated (active) TEK protein. An enucleated eye from an affected member of family 8 was examined via histology. SVEP1 expression in developing outflow tissues was detected by immunofluorescent staining of 7-day mouse anterior segments. SVEP1 stimulation of TEK expression in human umbilical vascular endothelial cells (HUVECs) was measured by TaqMan quantitative PCR. Results: Heterozygous TEK loss-of-function alleles were identified in eight PCG families, with parent-child disease transmission observed in two pedigrees. Family 8 exhibited greater disease penetrance and severity, histology revealed absence of SC in one eye, and SVEP1:p.R997C was identified in four of the five affected individuals. During SC development, SVEP1 is secreted by surrounding tissues. SVEP1:p.R997C abrogates stimulation of TEK expression by HUVECs. Conclusions: We provide further evidence for PCG caused by TEK haploinsufficiency, affirm autosomal dominant inheritance in two pedigrees, and propose SVEP1 as a modifier of TEK expression during SC development, affecting disease penetrance and severity.

Gluconobacter dominates the gut microbiome of the Asian palm civet Paradoxurus hermaphroditus that produces kopi luwak.

Coffee beans derived from feces of the civet cat are used to brew coffee known as kopi luwak (the Indonesian words for coffee and palm civet, respectively), which is one of the most expensive coffees in the world owing to its limited supply and strong market demand. Recent metabolomics studies have revealed that kopi luwak metabolites differ from metabolites found in other coffee beans. To produce kopi luwak, coffee beans are first eaten by civet cats. It has been proposed that fermentation inside the civet cat digestive tract may contribute to the distinctively smooth flavor of kopi luwak, but the biological basis has not been determined. Therefore, we characterized the microbiome of civet cat feces using 16S rRNA gene sequences to determine the bacterial taxa that may influence fermentation processes related to kopi luwak. Moreover, we compared this fecal microbiome with that of 14 other animals, revealing that Gluconobacter is a genus that is, uniquely found in feces of the civet cat. We also found that Gluconobacter species have a large number of cell motility genes, which may encode flagellar proteins allowing colonization of the civet gut. In addition, genes encoding enzymes involved in the metabolism of hydrogen sulfide and sulfur-containing amino acids were over-represented in Gluconobacter. These genes may contribute to the fermentation of coffee beans in the digestive tract of civet cats.

AllerCatPro-prediction of protein allergenicity potential from the protein sequence.

MOTIVATION: Due to the risk of inducing an immediate Type I (IgE-mediated) allergic response, proteins intended for use in consumer products must be investigated for their allergenic potential before introduction into the marketplace. The FAO/WHO guidelines for computational assessment of allergenic potential of proteins based on short peptide hits and linear sequence window identity thresholds misclassify many proteins as allergens. RESULTS: We developed AllerCatPro which predicts the allergenic potential of proteins based on similarity of their 3D protein structure as well as their amino acid sequence compared with a data set of known protein allergens comprising of 4180 unique allergenic protein sequences derived from the union of the major databases Food Allergy Research and Resource Program, Comprehensive Protein Allergen Resource, WHO/International Union of Immunological Societies, UniProtKB and Allergome. We extended the hexamer hit rule by removing peptides with high probability of random occurrence measured by sequence entropy as well as requiring 3 or more hexamer hits consistent with natural linear epitope patterns in known allergens. This is complemented with a Gluten-like repeat pattern detection. We also switched from a linear sequence window similarity to a B-cell epitope-like 3D surface similarity window which became possible through extensive 3D structure modeling covering the majority (74%) of allergens. In case no structure similarity is found, the decision workflow reverts to the old linear sequence window rule. The overall accuracy of AllerCatPro is 84% compared with other current methods which range from 51 to 73%. Both the FAO/WHO rules and AllerCatPro achieve highest sensitivity but AllerCatPro provides a 37-fold increase in specificity. AVAILABILITY AND IMPLEMENTATION: https://allercatpro.bii.a-star.edu.sg/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Discovering novel SNPs that are correlated with patient outcome in a Singaporean cancer patient cohort treated with gemcitabine-based chemotherapy.

BACKGROUND: Single Nucleotide Polymorphisms (SNPs) can influence patient outcome such as drug response and toxicity after drug intervention. The purpose of this study is to develop a systematic pathway approach to accurately and efficiently predict novel non-synonymous SNPs (nsSNPs) that could be causative to gemcitabine-based chemotherapy treatment outcome in Singaporean non-small cell lung cancer (NSCLC) patients. METHODS: Using a pathway approach that incorporates comprehensive protein-protein interaction data to systematically extend the gemcitabine pharmacologic pathway, we identified 77 related nsSNPs, common in the Singaporean population. After that, we used five computational criteria to prioritize the SNPs based on their importance for protein function. We specifically selected and screened six candidate SNPs in a patient cohort with NSCLC treated with gemcitabine-based chemotherapy. RESULT: We performed survival analysis followed by hematologic toxicity analyses and found that three of six candidate SNPs are significantly correlated with the patient outcome (P < 0.05) i.e. ABCG2 Q141K (rs2231142), SLC29A3 S158F (rs780668) and POLR2A N764K (rs2228130). CONCLUSIONS: Our computational SNP candidate enrichment workflow approach was able to identify several high confidence biomarkers predictive for personalized drug treatment outcome while providing a rationale for a molecular mechanism of the SNP effect. TRIAL REGISTRATION: NCT00695994. Registered 10 June, 2008 'retrospectively registered'.

Genome-Wide Analysis of Protein-Coding Variants in Leprosy.

Although genome-wide association studies have greatly advanced our understanding of the contribution of common noncoding variants to leprosy susceptibility, protein-coding variants have not been systematically investigated. We carried out a three-stage genome-wide association study of protein-coding variants in Han Chinese, of whom were 7,048 leprosy patients and 14,398 were healthy control subjects. Seven coding variants of exome-wide significance were discovered, including two rare variants: rs145562243 in NCKIPSD (P = 1.71 × 10-9, odds ratio [OR] = 4.35) and rs149308743 in CARD9 (P = 2.09 × 10-8, OR = 4.75); three low-frequency variants: rs76418789 in IL23R (P = 1.03 × 10-10, OR = 1.36), rs146466242 in FLG (P = 3.39 × 10-12, OR = 1.45), and rs55882956 in TYK2 (P = 1.04 × 10-6, OR = 1.30); and two common variants: rs780668 in SLC29A3 (P = 2.17 × 10-9, OR = 1.14) and rs181206 in IL27 (P = 1.08 × 10-7, OR = 0.83). Discovered protein-coding variants, particularly low-frequency and rare ones, showed involvement of skin barrier and endocytosis/phagocytosis/autophagy, in addition to known innate and adaptive immunity, in the pathogenesis of leprosy, highlighting the merits of protein-coding variant studies for complex diseases.

Angiopoietin receptor TEK mutations underlie primary congenital glaucoma with variable expressivity.

Primary congenital glaucoma (PCG) is a devastating eye disease and an important cause of childhood blindness worldwide. In PCG, defects in the anterior chamber aqueous humor outflow structures of the eye result in elevated intraocular pressure (IOP); however, the genes and molecular mechanisms involved in the etiology of these defects have not been fully characterized. Previously, we observed PCG-like phenotypes in transgenic mice that lack functional angiopoietin-TEK signaling. Herein, we identified rare TEK variants in 10 of 189 unrelated PCG families and demonstrated that each mutation results in haploinsufficiency due to protein loss of function. Multiple cellular mechanisms were responsible for the loss of protein function resulting from individual TEK variants, including an absence of normal protein production, protein aggregate formation, enhanced proteasomal degradation, altered subcellular localization, and reduced responsiveness to ligand stimulation. Further, in mice, hemizygosity for Tek led to the formation of severely hypomorphic Schlemm's canal and trabecular meshwork, as well as elevated IOP, demonstrating that anterior chamber vascular development is sensitive to Tek gene dosage and the resulting decrease in angiopoietin-TEK signaling. Collectively, these results identify TEK mutations in patients with PCG that likely underlie disease and are transmitted in an autosomal dominant pattern with variable expressivity.

GGIP: Structure and sequence-based GPCR-GPCR interaction pair predictor.

G Protein-Coupled Receptors (GPCRs) are important pharmaceutical targets. More than 30% of currently marketed pharmaceutical medicines target GPCRs. Numerous studies have reported that GPCRs function not only as monomers but also as homo- or hetero-dimers or higher-order molecular complexes. Many GPCRs exert a wide variety of molecular functions by forming specific combinations of GPCR subtypes. In addition, some GPCRs are reportedly associated with diseases. GPCR oligomerization is now recognized as an important event in various biological phenomena, and many researchers are investigating this subject. We have developed a support vector machine (SVM)-based method to predict interacting pairs for GPCR oligomerization, by integrating the structure and sequence information of GPCRs. The performance of our method was evaluated by the Receiver Operating Characteristic (ROC) curve. The corresponding area under the curve was 0.938. As far as we know, this is the only prediction method for interacting pairs among GPCRs. Our method could accelerate the analyses of these interactions, and contribute to the elucidation of the global structures of the GPCR networks in membranes. Proteins 2016; 84:1224-1233. © 2016 Wiley Periodicals, Inc.

HPMV: human protein mutation viewer - relating sequence mutations to protein sequence architecture and function changes.

Next-generation sequencing advances are rapidly expanding the number of human mutations to be analyzed for causative roles in genetic disorders. Our Human Protein Mutation Viewer (HPMV) is intended to explore the biomolecular mechanistic significance of non-synonymous human mutations in protein-coding genomic regions. The tool helps to assess whether protein mutations affect the occurrence of sequence-architectural features (globular domains, targeting signals, post-translational modification sites, etc.). As input, HPMV accepts protein mutations - as UniProt accessions with mutations (e.g. HGVS nomenclature), genome coordinates, or FASTA sequences. As output, HPMV provides an interactive cartoon showing the mutations in relation to elements of the sequence architecture. A large variety of protein sequence architectural features were selected for their particular relevance to mutation interpretation. Clicking a sequence feature in the cartoon expands a tree view of additional information including multiple sequence alignments of conserved domains and a simple 3D viewer mapping the mutation to known PDB structures, if available. The cartoon is also correlated with a multiple sequence alignment of similar sequences from other organisms. In cases where a mutation is likely to have a straightforward interpretation (e.g. a point mutation disrupting a well-understood targeting signal), this interpretation is suggested. The interactive cartoon can be downloaded as standalone viewer in Java jar format to be saved and viewed later with only a standard Java runtime environment. The HPMV website is: http://hpmv.bii.a-star.edu.sg/ .

Novel SNP improves differential survivability and mortality in non-small cell lung cancer patients.

BACKGROUND: Non-small cell lung cancer (NSCLC) is a major cause of cancer-related death worldwide due to poor patient prognosis and clinical outcome. Here, we studied the genetic variations underlying NSCLC pathogenesis based on their association to patient outcome after gemcitabine therapy. RESULTS: Bioinformatics analysis was used to investigate possible effects of POLA2 G583R (POLA2+1747 GG/GA, dbSNP ID: rs487989) in terms of protein function. Using biostatistics, POLA2+1747 GG/GA (rs487989, POLA2 G583R) was identified as strongly associated with mortality rate and survival time among NSCLC patients. It was also shown that POLA2+1747 GG/GA is functionally significant for protein localization via green fluorescent protein (GFP)-tagging and confocal laser scanning microscopy analysis. The single nucleotide polymorphism (SNP) causes DNA polymerase alpha subunit B to localize in the cytoplasm instead of the nucleus. This inhibits DNA replication in cancer cells and confers a protective effect in individuals with this SNP. CONCLUSIONS: The results suggest that POLA2+1747 GG/GA may be used as a prognostic biomarker of patient outcome in NSCLC pathogenesis.

Mutations in SCO2 are associated with autosomal-dominant high-grade myopia.

Myopia, or near-sightedness, is an ocular refractive error of unfocused image quality in front of the retinal plane. Individuals with high-grade myopia (dioptric power greater than -6.00) are predisposed to ocular morbidities such as glaucoma, retinal detachment, and myopic maculopathy. Nonsyndromic, high-grade myopia is highly heritable, and to date multiple gene loci have been reported. We performed exome sequencing in 4 individuals from an 11-member family of European descent from the United States. Affected individuals had a mean dioptric spherical equivalent of -22.00 sphere. A premature stop codon mutation c.157C>T (p.Gln53*) cosegregating with disease was discovered within SCO2 that maps to chromosome 22q13.33. Subsequent analyses identified three additional mutations in three highly myopic unrelated individuals (c.341G>A, c.418G>A, and c.776C>T). To determine differential gene expression in a developmental mouse model, we induced myopia by applying a -15.00D lens over one eye. Messenger RNA levels of SCO2 were significantly downregulated in myopic mouse retinae. Immunohistochemistry in mouse eyes confirmed SCO2 protein localization in retina, retinal pigment epithelium, and sclera. SCO2 encodes for a copper homeostasis protein influential in mitochondrial cytochrome c oxidase activity. Copper deficiencies have been linked with photoreceptor loss and myopia with increased scleral wall elasticity. Retinal thinning has been reported with an SC02 variant. Human mutation identification with support from an induced myopic animal provides biological insights of myopic development.

Whole exome sequencing identifies a mutation for a novel form of corneal intraepithelial dyskeratosis.

BACKGROUND: Corneal intraepithelial dyskeratosis is an extremely rare condition. The classical form, affecting Native American Haliwa-Saponi tribe members, is called hereditary benign intraepithelial dyskeratosis (HBID). Herein, we present a new form of corneal intraepithelial dyskeratosis for which we identified the causative gene by using deep sequencing technology. METHODS AND RESULTS: A seven member Caucasian French family with two corneal intraepithelial dyskeratosis affected individuals (6-year-old proband and his mother) was ascertained. The proband presented with bilateral complete corneal opacification and dyskeratosis. Palmoplantar hyperkeratosis and laryngeal dyskeratosis were associated with the phenotype. Histopathology studies of cornea and vocal cord biopsies showed dyskeratotic keratinisation. Quantitative PCR ruled out 4q35 duplication, classically described in HBID cases. Next generation sequencing with mean coverage of 50× using the Illumina Hi Seq and whole exome capture processing was performed. Sequence reads were aligned, and screened for single nucleotide variants and insertion/deletion calls. In-house pipeline filtering analyses and comparisons with available databases were performed. A novel missense mutation M77T was discovered for the gene NLRP1 which maps to chromosome 17p13.2. This was a de novo mutation in the proband's mother, following segregation in the family, and not found in 738 control DNA samples. NLRP1 expression was determined in adult corneal epithelium. The amino acid change was found to destabilise significantly the protein structure. CONCLUSIONS: We describe a new corneal intraepithelial dyskeratosis and how we identified its causative gene. The NLRP1 gene product is implicated in inflammation, autoimmune disorders, and caspase mediated apoptosis. NLRP1 polymorphisms are associated with various diseases.