Characterization of all possible single-nucleotide change caused amino acid substitutions in the kinase domain of Bruton tyrosine kinase

PON-Fold: Prediction of Substitutions Affecting Protein Folding Rate

Most proteins fold into characteristic three-dimensional structures. The rate of folding and unfolding varies widely and can be affected by variations in proteins. We developed a novel machine-learning-based method for the prediction of the folding rate effects of amino acid substitutions in two-state folding proteins. We collected a data set of experimentally defined folding rates for variants and used them to train a gradient boosting algorithm starting with 1161 features. Two predictors were designed. The three-class classifier had, in blind tests, specificity and sensitivity ranging from 0.324 to 0.419 and from 0.256 to 0.451, respectively. The other tool was a regression predictor that showed a Pearson correlation coefficient of 0.525. The error measures, mean absolute error and mean squared error, were 0.581 and 0.603, respectively. One of the previously presented tools could be used for comparison with the blind test data set, our method called PON-Fold showed superior performance on all used measures. The applicability of the tool was tested by predicting all possible substitutions in a protein domain. Predictions for different conformations of proteins, open and closed forms of a protein kinase, and apo and holo forms of an enzyme indicated that the choice of the structure had a large impact on the outcome. PON-Fold is freely available.

Site-directed modification of multifunctional lignocellulose-degrading enzymes of straw based on homologous modeling

Studying the straw lignocellulose strengthening mechanism during simultaneous degradation has important practical significance for improving resource utilization and reducing environmental pollution. In this paper, the degradation ability of four straw lignocellulose-degrading enzymes was evaluated by molecular docking and molecular dynamics. Using the significantly binds to straw lignocellulose-degrading enzyme as a template, a multifunctional lignocellulose-degrading enzyme 3CBH-1KS5-4XQD-1B85 was constructed based on amino acid recombination and homologous modeling. Five efficient degrading enzymes (3CBH-1, 3CBH-2, 3CBH-3, 3CBH-4, and 3CBH-5) were designed by site-directed mutagenesis of 3CBH-1KS5-4XQD-1B85 amino acid at position 346. Molecular dynamics showed that the degradation ability of 3CBH-1 was significant and it was 1.45 times higher than 3CBH-1KS5-4XQD-1B85. Moreover, the mechanism of enhanced degradability and the stability of the enzymes were explored. With the aid of Taguchi experiments, the suitable external environment for degrading straw was determined. In the presence of inhibitors (organic acids and phenolic compounds), the binding energy of 3CBH-1 (238.46 ± 30.96 kJ/mol) is 36.42% higher than that of 3CBH-1KS5-4XQD-1B85 (174.79 ± 20.35 kJ/mol) without external environmental stimulation. Based on homology modeling, this paper constructed a site-directed mutagenesis scheme of multifunctional enzymes, and the aim was to obtain multifunctional and efficient straw lignocellulose-degrading enzymes through protein engineering, which provided a feasible scheme for straw biodegradation.

Research-based flow cytometry assays for pathogenic assessment in the human B-cell biology of gene variants revealed in the diagnosis of inborn errors of immunity: a Bruton's tyrosine kinase case-study

Introduction

Inborn errors of immunity (IEI) are an expanding group of rare diseases whose field has been boosted by next-generation sequencing (NGS), revealing several new entities, accelerating routine diagnoses, expanding the number of atypical presentations and generating uncertainties regarding the pathogenic relevance of several novel variants.

Methods

Research laboratories that diagnose and provide support for IEI require accurate, reproducible and sustainable phenotypic, cellular and molecular functional assays to explore the pathogenic consequences of human leukocyte gene variants and contribute to their assessment. We have implemented a set of advanced flow cytometry-based assays to better dissect human B-cell biology in a translational research laboratory. We illustrate the utility of these techniques for the in-depth characterization of a novel (c.1685G>A, p.R562Q) de novo gene variant predicted as probably pathogenic but with no previous insights into the protein and cellular effects, located in the tyrosine kinase domain of the Bruton's tyrosine kinase (BTK) gene, in an apparently healthy 14-year-old male patient referred to our clinic for an incidental finding of low immunoglobulin (Ig) M levels with no history of recurrent infections.

Results and discussion

A phenotypic analysis of bone marrow (BM) revealed a slightly high percentage of pre-B-I subset in BM, with no blockage at this stage, as typically observed in classical X-linked agammaglobulinemia (XLA) patients. The phenotypic analysis in peripheral blood also revealed reduced absolute numbers of B cells, all pre-germinal center maturation stages, together with reduced but detectable numbers of different memory and plasma cell isotypes. The R562Q variant allows Btk expression and normal activation of anti-IgM-induced phosphorylation of Y551 but diminished autophosphorylation at Y223 after anti IgM and CXCL12 stimulation. Lastly, we explored the potential impact of the variant protein for downstream Btk signaling in B cells. Within the canonical nuclear factor kappa B (NF-κB) activation pathway, normal IκBα degradation occurs after CD40L stimulation in patient and control cells. In contrast, disturbed IκBα degradation and reduced calcium ion (Ca2⁺) influx occurs on anti-IgM stimulation in the patient's B cells, suggesting an enzymatic impairment of the mutated tyrosine kinase domain.

Individual Genetic Heterogeneity

Genetic variation has been widely covered in literature, however, not from the perspective of an individual in any species. Here, a synthesis of genetic concepts and variations relevant for individual genetic constitution is provided. All the different levels of genetic information and variation are covered, ranging from whether an organism is unmixed or hybrid, has variations in genome, chromosomes, and more locally in DNA regions, to epigenetic variants or alterations in selfish genetic elements. Genetic constitution and heterogeneity of microbiota are highly relevant for health and wellbeing of an individual. Mutation rates vary widely for variation types, e.g., due to the sequence context. Genetic information guides numerous aspects in organisms. Types of inheritance, whether Mendelian or non-Mendelian, zygosity, sexual reproduction, and sex determination are covered. Functions of DNA and functional effects of variations are introduced, along with mechanism that reduce and modulate functional effects, including TARAR countermeasures and intraindividual genetic conflict. TARAR countermeasures for tolerance, avoidance, repair, attenuation, and resistance are essential for life, integrity of genetic information, and gene expression. The genetic composition, effects of variations, and their expression are considered also in diseases and personalized medicine. The text synthesizes knowledge and insight on individual genetic heterogeneity and organizes and systematizes the central concepts.

Targeted RNAseq Improves Clinical Diagnosis of Very Early-Onset Pediatric Immune Dysregulation

Despite increased use of whole exome sequencing (WES) for the clinical analysis of rare disease, overall diagnostic yield for most disorders hovers around 30%. Previous studies of mRNA have succeeded in increasing diagnoses for clearly defined disorders of monogenic inheritance. We asked if targeted RNA sequencing could provide similar benefits for primary immunodeficiencies (PIDs) and very early-onset inflammatory bowel disease (VEOIBD), both of which are difficult to diagnose due to high heterogeneity and variable severity. We performed targeted RNA sequencing of a panel of 260 immune-related genes for a cohort of 13 patients (seven suspected PID cases and six VEOIBD) and analyzed variants, splicing, and exon usage. Exonic variants were identified in seven cases, some of which had been previously prioritized by exome sequencing. For four cases, allele specific expression or lack thereof provided additional insights into possible disease mechanisms. In addition, we identified five instances of aberrant splicing associated with four variants. Three of these variants had been previously classified as benign in ClinVar based on population frequency. Digenic or oligogenic inheritance is suggested for at least two patients. In addition to validating the use of targeted RNA sequencing, our results show that rare disease research will benefit from incorporating contributing genetic factors into the diagnostic approach.

Identification of four novel mutations in BTK from six Chinese families with X-linked agammaglobulinemia

BACKGROUND

The defect of Bruton's tyrosine kinase (BTK) gene resulted in X-linked agammaglobulinemia (XLA), which is characterized by recurrent bacterial infections, immunodeficiency with low B-cell numbers and immunoglobulin. Diagnosis of XLA depends on clinical phenotype and genetic testing.

METHODS

Six unrelated Chinese families with high suspicion of XLA were enrolled in this study. Potential pathogenic variants were detected and validated by Whole Exome Sequencing (WES) and Sanger Sequencing. Western blot, Quantitative PCR (qPCR) analysis and immunofluorescence analysis were used to evaluate the preliminary function of candidate BTK variants.

RESULTS

A total of six variants were identified, four of which were not reported before. The novel missense mutation(c.1900T>G) and deletion(c.897delG) were found that the mutant protein and mRNA expression levels have fallen by Western Blot and qPCR identification. We also constructed minigene expression vector to determine the deletion (c.1751-6_1755delttctagGGGTT) resulting a 35bp skipping in exon 18. Meanwhile, the break point of gross deletion (Exon2-5) discovered based on WES was confirmed to be located at site ChX:101367539_101376531 through qPCR and Gap-PCR.

CONCLUSION

This study makes definitive diagnosis for 6 families with suspected XLA and further expands the spectrum of BTK mutations, providing new information for the diagnosis of the disease.

Structure-Function Relationships of Covalent and Non-Covalent BTK Inhibitors

Low-molecular weight chemical compounds have a longstanding history as drugs. Target specificity and binding efficiency represent major obstacles for small molecules to become clinically relevant. Protein kinases are attractive cellular targets; however, they are challenging because they present one of the largest protein families and share structural similarities. Bruton tyrosine kinase (BTK), a cytoplasmic protein tyrosine kinase, has received much attention as a promising target for the treatment of B-cell malignancies and more recently autoimmune and inflammatory diseases. Here we describe the structural properties and binding modes of small-molecule BTK inhibitors, including irreversible and reversible inhibitors. Covalently binding compounds, such as ibrutinib, acalabrutinib and zanubrutinib, are discussed along with non-covalent inhibitors fenebrutinib and RN486. The focus of this review is on structure-function relationships.

Enhanced plant-microbe remediation of PCBs in soil using enzyme modification technique combined with molecular docking and molecular dynamics

The study on the enhanced mechanisms of the enzymes involved in plant absorption, plant degradation, and microbial mineralization in the remediation of soils contaminated with polychlorinated biphenyls (PCBs) is of great significance for the application of plant-microbe combined remediation technique in PCB-contaminated soils. The present study first used a combination of molecular docking and molecular dynamics methods to calculate the effects of the plant absorption enzyme, plant degradation enzyme, and microbial mineralization enzyme on the PCBs in the soil environment. A multifunctional plant degradation enzyme was constructed with three functional roles of absorption, degradation, and mineralization using amino acid sequence recombination and site-directed mutagenesis to modify the template of plant degradation enzyme. Finally, using the Taguchi experimental design-assisted molecular dynamics simulation method, the suitable external environmental conditions of plant-microbe combined remediation of the PCB-contaminated soil were determined. In total, six multifunctional plant degradation enzymes were designed, which exhibited a significantly improved efficiency of PCB degradation. In comparison to the complex of plant absorption enzyme, plant degradation enzyme, and microorganism mineralization enzyme (6QIM-3GZX-1B85), the six multifunctional plant degradation enzymes exhibited significantly higher efficiency (2.10-2.38 times) in degrading the PCBs, with a maximum of 2.69 times under suitable external environmental conditions.

Primary Immunodeficiencies in India: Molecular Diagnosis and the Role of Next-Generation Sequencing

Primary immunodeficiency diseases (PIDs) are a group of clinically and genetically heterogeneous disorders showing ethnic and geographic diversities. Next-generation sequencing (NGS) is a comprehensive tool to diagnose PID. Although PID is common in India, data on the genetic spectrum of PIDs are limited due to financial restrictions. The study aims to characterize the clinical and genetic spectrum of PID patients in India and highlight the importance of a cost-effective targeted gene panel sequencing approach for PID in a resource-limited setting. The study includes 229 patients with clinical and laboratory features suggestive of PIDs. Mutation analysis was done by Sanger sequencing and NGS targeting a customized panel of genes. Pathogenic variants were identified in 97 patients involving 42 different genes with BTK and IL12RB1 being the most common mutated genes. Autosomal recessive and X-linked recessive inheritance were seen in 51.6% and 23.7% of patients. Mendelian susceptibility to mycobacterial diseases (MSMD) and IL12RB1 mutations was more common in our population compared to the Western world and the Middle East. Two patients with hypomorphic RAG1 mutations and one female with skewed CYBB mutation were also identified. Another 40 patients had variants classified as variants of uncertain significance (VUS). The study shows that targeted NGS is an effective diagnostic strategy for PIDs in countries with limited diagnostic resources. Molecular diagnosis of PID helps in genetic counseling and to make therapeutic decisions including the need for a stem cell transplantation.

Variation benchmark datasets: update, criteria, quality and applications

Development of new computational methods and testing their performance has to be carried out using experimental data. Only in comparison to existing knowledge can method performance be assessed. For that purpose, benchmark datasets with known and verified outcome are needed. High-quality benchmark datasets are valuable and may be difficult, laborious and time consuming to generate. VariBench and VariSNP are the two existing databases for sharing variation benchmark datasets used mainly for variation interpretation. They have been used for training and benchmarking predictors for various types of variations and their effects. VariBench was updated with 419 new datasets from 109 papers containing altogether 329 014 152 variants; however, there is plenty of redundancy between the datasets. VariBench is freely available at http://structure.bmc.lu.se/VariBench/. The contents of the datasets vary depending on information in the original source. The available datasets have been categorized into 20 groups and subgroups. There are datasets for insertions and deletions, substitutions in coding and non-coding region, structure mapped, synonymous and benign variants. Effect-specific datasets include DNA regulatory elements, RNA splicing, and protein property for aggregation, binding free energy, disorder and stability. Then there are several datasets for molecule-specific and disease-specific applications, as well as one dataset for variation phenotype effects. Variants are often described at three molecular levels (DNA, RNA and protein) and sometimes also at the protein structural level including relevant cross references and variant descriptions. The updated VariBench facilitates development and testing of new methods and comparison of obtained performances to previously published methods. We compared the performance of the pathogenicity/tolerance predictor PON-P2 to several benchmark studies, and show that such comparisons are feasible and useful, however, there may be limitations due to lack of provided details and shared data.

Database URL: http://structure.bmc.lu.se/VariBench

Functional effects of protein variants

Genetic and other variations frequently affect protein functions. Scientific articles can contain confusing descriptions about which function or property is affected, and in many cases the statements are pure speculation without any experimental evidence. To clarify functional effects of protein variations of genetic or non-genetic origin, a systematic conceptualisation and framework are introduced. This framework describes protein functional effects on abundance, activity, specificity and affinity, along with countermeasures, which allow cells, tissues and organisms to tolerate, avoid, repair, attenuate or resist (TARAR) the effects. Effects on abundance discussed include gene dosage, restricted expression, mis-localisation and degradation. Enzymopathies, effects on kinetics, allostery and regulation of protein activity are subtopics for the effects of variants on activity. Variation outcomes on specificity and affinity comprise promiscuity, specificity, affinity and moonlighting. TARAR mechanisms redress variations with active and passive processes including chaperones, redundancy, robustness, canalisation and metabolic and signalling rewiring. A framework for pragmatic protein function analysis and presentation is introduced. All of the mechanisms and effects are described along with representative examples, most often in relation to diseases. In addition, protein function is discussed from evolutionary point of view. Application of the presented framework facilitates unambiguous, detailed and specific description of functional effects and their systematic study.

Impact of amino acid substitution in the kinase domain of Bruton tyrosine kinase and its association with X-linked agammaglobulinemia

X-linked agammaglobulinemia (XLA) is a rare disease that affects the immune system, characterized by a serial development of bacterial infection from the onset of infantile age. Bruton tyrosine kinase (BTK) is a non-receptor cytoplasmic kinase that plays a crucial role in the B-lymphocyte maturation. The altered expression, mutation and/or structural variations of BTK are responsible for causing XLA. Here, we have performed extensive sequence and structure analyses of BTK to find deleterious variations and their pathogenic association with XLA. First, we screened the pathogenic variations in the BTK from a pool of publicly available resources, and their pathogenicity/tolerance and stability predictions were carried out. Finally, two pathogenic variations (E589G and M630K) were studied in detail and subjected to all-atom molecular dynamics simulation for 200 ns. Intramolecular hydrogen bonds (H-bonds), secondary structure, and principal component analysis revealed significant conformational changes in variants that support the structural basis of BTK dysfunction in XLA. The free energy landscape analysis revealed the presence of multiple energy minima, suggests that E589G brings a large destabilization and consequently unfolding behavior compared to M630K. Overall, our study suggests that amino acid substitutions, E589G, and M630K, significantly alter the structural conformation and stability of BTK.

Identification of biochemically neutral positions in liver pyruvate kinase

Understanding how each residue position contributes to protein function has been a long-standing goal in protein science. Substitution studies have historically focused on conserved protein positions. However, substitutions of nonconserved positions can also modify function. Indeed, we recently identified nonconserved positions that have large substitution effects in human liver pyruvate kinase (hLPYK), including altered allosteric coupling. To facilitate a comparison of which characteristics determine when a nonconserved position does vs does not contribute to function, the goal of the current work was to identify neutral positions in hLPYK. However, existing hLPYK data showed that three features commonly associated with neutral positions-high sequence entropy, high surface exposure, and alanine scanning-lacked the sensitivity needed to guide experimental studies. We used multiple evolutionary patterns identified in a sequence alignment of the PYK family to identify which positions were least patterned, reasoning that these were most likely to be neutral. Nine positions were tested with a total of 117 amino acid substitutions. Although exploring all potential functions is not feasible for any protein, five parameters associated with substrate/effector affinities and allosteric coupling were measured for hLPYK variants. For each position, the aggregate functional outcomes of all variants were used to quantify a "neutrality" score. Three positions showed perfect neutral scores for all five parameters. Furthermore, the nine positions showed larger neutral scores than 17 positions located near allosteric binding sites. Thus, our strategy successfully enriched the dataset for positions with neutral and modest substitutions.

New insights into the pathogenicity of non-synonymous variants through multi-level analysis

Precise classification of non-synonymous single nucleotide variants (SNVs) is a fundamental goal of clinical genetics. Next-generation sequencing technology is effective for establishing the basis of genetic diseases. However, identification of variants that are causal for genetic diseases remains a challenge. We analyzed human non-synonymous SNVs from a multilevel perspective to characterize pathogenicity. We showed that computational tools, though each having its own strength and weakness, tend to be overly dependent on the degree of conservation. For the mutations at non-degenerate sites, the amino acid sites of pathogenic substitutions show a distinct distribution in the classes of protein domains compared with the sites of benign substitutions. Overlooked disease susceptibility of genes explains in part the failures of computational tools. The more pathogenic sites observed, the more likely the gene is expressed in a high abundance or in a high tissue-specific manner, and have a high node degree of protein-protein interaction. The destroyed functions due to some false-negative mutations may arise because of a reprieve from the epigenetic repressed state which shouldn't happen in multiple biological conditions, instead of the defective protein. Our work adds more to our knowledge of non-synonymous SNVs' pathogenicity, thus will benefit the field of clinical genetics.

Genetic screening of male patients with primary hypogammaglobulinemia can guide diagnosis and clinical management

The precise diagnosis of an immunodeficiency is sometimes difficult to assess, especially due to the large spectrum of phenotypic variation reported among patients. Common variable immunodeficiency disorders (CVID) do not have, for a large part, an identified genetic cause. The identification of a causal genetic mutation is important to confirm, or in some cases correct, the diagnosis. We screened >150 male patients with hypogammaglobulinemia for mutations in three genes involved in pediatric X-linked primary immunoglobulin deficiency: CD40LG, SH2D1A and BTK. The SH2D1A screening allowed to reclassify two individuals with an initial CVID presentation as XLP after mutations identification. All these mutations were associated with a lack of protein expression. In addition, 4 patients with a primary diagnosis of CVID and one with a primary IgG subclass deficiency were requalified as XLA after identifying BTK mutations. Interestingly, two out of these 5 patients carried a damaging coding BTK mutation associated with a lower, but detectable, BTK expression in monocytes, suggesting that a dysfunctional protein explains the disease phenotype in these patients. In conclusion, our results advocate to include SH2D1A and BTK in newly developed targeted NGS genetic testing, to contribute to providing the most appropriate medical treatment and genetic counselling.

Genome-Wide Identification and Characterization of Tyrosine Kinases in the Silkworm, Bombyx mori

The tyrosine kinases (TKs) are important parts of metazoan signaling pathways and play significant roles in cell growth, development, apoptosis and disease. Genome-wide characterization of TKs has been conducted in many metazoans, however, systematic information about this family in Lepidoptera is still lacking. We retrieved 33 TK-encoding genes in silkworm and classified them into 25 subfamilies by sequence analysis, without members in AXL, FRK, PDGFR, STYK1 and TIE subfamilies. Although domain sequences in each subfamily are conserved, TKs in vertebrates tend to be remarkably conserved and stable. Our results of phylogenetic analysis supported the previous conclusion for the second major expansion of TK family. Gene-Ontology (GO) analysis revealed that a higher proportion of BmTKs played roles in binding, catalysis, signal transduction, metabolism, biological regulation and response to stimulus, compared to all silkworm genes annotated in GO. Moreover, the expression profile analysis of BmTKs among multiple tissues and developmental stages demonstrated that many genes exhibited stage-specific and/or sex-related expression during embryogenesis, molting and metamorphosis, and that 8 BmTKs presented tissue-specific high expression. Our study provides systematic description of silkworm tyrosine kinases, and may also provide further insights into metazoan TKs and assist future studies addressing their functions.

PON-SC - program for identifying steric clashes caused by amino acid substitutions

Background

Amino acid substitutions due to DNA nucleotide replacements are frequently disease-causing because of affecting functionally important sites. If the substituting amino acid does not fit into the protein, it causes structural alterations that are often harmful. Clashes of amino acids cause local or global structural changes. Testing structural compatibility of variations has been difficult due to the lack of a dedicated method that could handle vast amounts of variation data produced by next generation sequencing technologies.

Results

We developed a method, PON-SC, for detecting protein structural clashes due to amino acid substitutions. The method utilizes side chain rotamer library and tests whether any of the common rotamers can be fitted into the protein structure. The tool was tested both with variants that cause and do not cause clashes and found to have accuracy of 0.71 over five test datasets.

Conclusions

We developed a fast method for residue side chain clash detection. The method provides in addition to the prediction also visualization of the variant in three dimensional structure.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1947-7) contains supplementary material, which is available to authorized users.

Identification and characterization of tyrosine kinases in anole lizard indicate the conserved tyrosine kinase repertoire in vertebrates

The tyrosine kinases (TKs) play principal roles in regulation of multicellular aspects of the organism and are implicated in many cancer types and congenital disorders. The anole lizard has recently been introduced as a model organism for laboratory-based studies of organismal function and field studies of ecology and evolution. However, the TK family of anole lizard has not been systematically identified and characterized yet. In this study, we identified 82 TK-encoding genes in the anole lizard genome and classified them into 28 subfamilies through phylogenetic analysis, with no member from ROS and STYK1 subfamilies identified. Although TK domain sequences and domain organization in each subfamily were conserved, the total number of TKs in different species was much variable. In addition, extensive evolutionary analysis in metazoans indicated that TK repertoire in vertebrates tends to be remarkably stable. Phylogenetic analysis of Eph subfamily indicated that the divergence of EphA and EphB occurred prior to the whole genome duplication (WGD) but after the split of Urochordates and vertebrates. Moreover, the expression pattern analysis of lizard TK genes among 9 different tissues showed that 14 TK genes exhibited tissue-specific expression and 6 TK genes were widely expressed. Comparative analysis of TK expression suggested that the tissue specifically expressed genes showed different expression pattern but the widely expressed genes showed similar pattern between anole lizard and human. These results may provide insights into the evolutionary diversification of animal TK genes and would aid future studies on TK protein regulation of key growth and developmental processes.

PON-P and PON-P2 predictor performance in CAGI challenges: Lessons learned

Computational tools are widely used for ranking and prioritizing variants for characterizing their disease relevance. Since numerous tools have been developed, they have to be properly assessed before being applied. Critical Assessment of Genome Interpretation (CAGI) experiments have significantly contributed toward the assessment of prediction methods for various tasks. Within and outside the CAGI, we have addressed several questions that facilitate development and assessment of variation interpretation tools. These areas include collection and distribution of benchmark datasets, their use for systematic large-scale method assessment, and the development of guidelines for reporting methods and their performance. For us, CAGI has provided a chance to experiment with new ideas, test the application areas of our methods, and network with other prediction method developers. In this article, we discuss our experiences and lessons learned from the various CAGI challenges. We describe our approaches, their performance, and impact of CAGI on our research. Finally, we discuss some of the possibilities that CAGI experiments have opened up and make some suggestions for future experiments.

From identification of the BTK kinase to effective management of leukemia

BTK is a cytoplasmic protein-tyrosine kinase, whose corresponding gene was isolated in the early 1990s. BTK was initially identified by positional cloning of the gene causing X-linked agammaglobulinemia and independently in a search for new kinases. Given the phenotype of affected patients, namely lack of B-lymphocytes and plasma cells with the ensuing inability to mount humoral immune responses, BTK inhibitors were anticipated to have beneficial effects on antibody-mediated pathologies, such as autoimmunity. In contrast to, for example, the SRC-family of cytoplasmic kinases, there was no obvious way in which structural alterations would yield constitutively active forms of BTK, and such mutations were also not found in leukemias or lymphomas. In 2007, the first efficient inhibitor, ibrutinib, was reported and soon became approved both in the United States and in Europe for the treatment of three B-cell malignancies, mantle cell lymphoma, chronic lymphocytic leukemia and Waldenström's macroglobulinemia. Over the past few years, additional inhibitors have been developed, with acalabrutinib being more selective, and recently demonstrating fewer clinical adverse effects. The antitumor mechanism is also not related to mutations in BTK. Instead tumor residency in lymphoid organs is inhibited, making these drugs highly versatile. BTK is one of the only 10 human kinases that carry a cysteine in the adenosine triphosphate-binding cleft. As this allows for covalent, irreversible inhibitor binding, it provides these compounds with a highly advantageous character. This quality may be crucial and bodes well for the future of BTK-modifying medicines, which have been estimated to reach annual multi-billion dollar sales in the future.

KinMutRF: a random forest classifier of sequence variants in the human protein kinase superfamily

Background

The association between aberrant signal processing by protein kinases and human diseases such as cancer was established long time ago. However, understanding the link between sequence variants in the protein kinase superfamily and the mechanistic complex traits at the molecular level remains challenging: cells tolerate most genomic alterations and only a minor fraction disrupt molecular function sufficiently and drive disease.

Results

KinMutRF is a novel random-forest method to automatically identify pathogenic variants in human kinases. Twenty six decision trees implemented as a random forest ponder a battery of features that characterize the variants: a) at the gene level, including membership to a Kinbase group and Gene Ontology terms; b) at the PFAM domain level; and c) at the residue level, the types of amino acids involved, changes in biochemical properties, functional annotations from UniProt, Phospho.ELM and FireDB. KinMutRF identifies disease-associated variants satisfactorily (Acc: 0.88, Prec:0.82, Rec:0.75, F-score:0.78, MCC:0.68) when trained and cross-validated with the 3689 human kinase variants from UniProt that have been annotated as neutral or pathogenic. All unclassified variants were excluded from the training set. Furthermore, KinMutRF is discussed with respect to two independent kinase-specific sets of mutations no included in the training and testing, Kin-Driver (643 variants) and Pon-BTK (1495 variants). Moreover, we provide predictions for the 848 protein kinase variants in UniProt that remained unclassified.

A public implementation of KinMutRF, including documentation and examples, is available online (http://kinmut2.bioinfo.cnio.es). The source code for local installation is released under a GPL version 3 license, and can be downloaded from https://github.com/Rbbt-Workflows/KinMut2.

Conclusions

KinMutRF is capable of classifying kinase variation with good performance. Predictions by KinMutRF compare favorably in a benchmark with other state-of-the-art methods (i.e. SIFT, Polyphen-2, MutationAssesor, MutationTaster, LRT, CADD, FATHMM, and VEST). Kinase-specific features rank as the most elucidatory in terms of information gain and are likely the improvement in prediction performance. This advocates for the development of family-specific classifiers able to exploit the discriminatory power of features unique to individual protein families.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2723-1) contains supplementary material, which is available to authorized users.

Substitution scanning identifies a novel, catalytically active ibrutinib-resistant BTK cysteine 481 to threonine (C481T) variant

Irreversible Bruton tyrosine kinase (BTK) inhibitors, ibrutinib and acalabrutinib have demonstrated remarkable clinical responses in multiple B-cell malignancies. Acquired resistance has been identified in a sub-population of patients in which mutations affecting BTK predominantly substitute cysteine 481 in the kinase domain for catalytically active serine, thereby ablating covalent binding of inhibitors. Activating substitutions in the BTK substrate phospholipase Cγ2 (PLCγ2) instead confers resistance independent of BTK. Herein, we generated all six possible amino acid substitutions due to single nucleotide alterations for the cysteine 481 codon, in addition to threonine, requiring two nucleotide substitutions, and performed functional analysis. Replacement by arginine, phenylalanine, tryptophan or tyrosine completely inactivated the catalytic activity, whereas substitution with glycine caused severe impairment. BTK with threonine replacement was catalytically active, similar to substitution with serine. We identify three potential ibrutinib resistance scenarios for cysteine 481 replacement: (1) Serine, being catalytically active and therefore predominating among patients. (2) Threonine, also being catalytically active, but predicted to be scarce, because two nucleotide changes are needed. (3) As BTK variants replaced with other residues are catalytically inactive, they presumably need compensatory mutations, therefore being very scarce. Glycine and tryptophan variants were not yet reported but likely also provide resistance.

Next Generation Sequencing Data Analysis in Primary Immunodeficiency Disorders - Future Directions

Primary immunodeficiency diseases (PIDs) comprise a group of highly heterogeneous immune system diseases and around 300 forms of PID have been described to date. Next Generation Sequencing (NGS) has recently become an increasingly used approach for gene identification and molecular diagnosis of human diseases. Herein we summarize the practical considerations for the interpretation of NGS data and the techniques for searching disease-related PID genes, and suggest future directions for research in this field.

Computational Analysis of the Binding Specificities of PH Domains

Pleckstrin homology (PH) domains share low sequence identities but extremely conserved structures. They have been found in many proteins for cellular signal-dependent membrane targeting by binding inositol phosphates to perform different physiological functions. In order to understand the sequence-structure relationship and binding specificities of PH domains, quantum mechanical (QM) calculations and sequence-based combined with structure-based binding analysis were employed in our research. In the structural aspect, the binding specificities were shown to correlate with the hydropathy characteristics of PH domains and electrostatic properties of the bound inositol phosphates. By comparing these structure properties with sequence-based profiles of physicochemical properties, PH domains can be classified into four functional subgroups according to their binding specificities and affinities to inositol phosphates. The method not only provides a simple and practical paradigm to predict binding specificities for functional genomic research but also gives new insight into the understanding of the basis of diseases with respect to PH domain structures.

wKinMut-2: Identification and Interpretation of Pathogenic Variants in Human Protein Kinases

Most genomic alterations are tolerated while only a minor fraction disrupts molecular function sufficiently to drive disease. Protein kinases play a central biological function and the functional consequences of their variants are abundantly characterized. However, this heterogeneous information is often scattered across different sources, which makes the integrative analysis complex and laborious. wKinMut-2 constitutes a solution to facilitate the interpretation of the consequences of human protein kinase variation. Nine methods predict their pathogenicity, including a kinase-specific random forest approach. To understand the biological mechanisms causative of human diseases and cancer, information from pertinent reference knowledge bases and the literature is automatically mined, digested, and homogenized. Variants are visualized in their structural contexts and residues affecting catalytic and drug binding are identified. Known protein-protein interactions are reported. Altogether, this information is intended to assist the generation of new working hypothesis to be corroborated with ulterior experimental work. The wKinMut-2 system, along with a user manual and examples, is freely accessible at http://kinmut2.bioinfo.cnio.es, the code for local installations can be downloaded from https://github.com/Rbbt-Workflows/KinMut2.

Classification of Amino Acid Substitutions in Mismatch Repair Proteins Using PON-MMR2

Variations in mismatch repair (MMR) system genes are causative of Lynch syndrome and other cancers. Thousands of variants have been identified in MMR genes, but the clinical relevance is known for only a small proportion. Recently, the InSiGHT group classified 2,360 MMR variants into five classes. One-third of variants, majority of which is nonsynonymous variants, remain to be of uncertain clinical relevance. Computational tools can be used to prioritize variants for disease relevance investigations. Previously, we classified 248 MMR variants as likely pathogenic and likely benign using PON-MMR. We have developed a novel tool, PON-MMR2, which is trained on a larger and more reliable dataset. In performance comparison, PON-MMR2 outperforms both generic tolerance prediction methods as well as methods optimized for MMR variants. It achieves accuracy and MCC of 0.89 and 0.78, respectively, in cross-validation and 0.86 and 0.69, respectively, on an independent test dataset. We classified 354 class 3 variants in InSiGHT database as well as all possible amino acid substitutions in four MMR proteins. Likely harmful variants mainly appear in the protein core, whereas likely benign variants are on the surface. PON-MMR2 is a highly reliable tool to prioritize variants for functional analysis. It is freely available at http://structure.bmc.lu.se/PON-MMR2/.

Harmful somatic amino acid substitutions affect key pathways in cancers

Background

Cancer is characterized by the accumulation of large numbers of genetic variations and alterations of multiple biological phenomena. Cancer genomics has largely focused on the identification of such genetic alterations and the genes containing them, known as ‘cancer genes’. However, the non-functional somatic variations out-number functional variations and remain as a major challenge. Recurrent somatic variations are thought to be cancer drivers but they are present in only a small fraction of patients.

Methods

We performed an extensive analysis of amino acid substitutions (AASs) from 6,861 cancer samples (whole genome or exome sequences) classified into 30 cancer types and performed pathway enrichment analysis. We also studied the overlap between the cancers based on proteins containing harmful AASs and pathways affected by them.

Results

We found that only a fraction of AASs (39.88 %) are harmful even in known cancer genes. In addition, we found that proteins containing harmful AASs in cancers are often centrally located in protein interaction networks. Based on the proteins containing harmful AASs, we identified significantly affected pathways in 28 cancer types and indicate that proteins containing harmful AASs can affect pathways despite the frequency of AASs in them. Our cross-cancer overlap analysis showed that it would be more beneficial to identify affected pathways in cancers rather than individual genes and variations.

Conclusion

Pathways affected by harmful AASs reveal key processes involved in cancer development. Our approach filters out the putative benign AASs thus reducing the list of cancer variations allowing reliable identification of affected pathways. The pathways identified in individual cancer and overlap between cancer types open avenues for further experimental research and for developing targeted therapies and interventions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-015-0125-x) contains supplementary material, which is available to authorized users.