Rational design and purification of human Bruton's tyrosine kinase SH3-SH2 protein for structure-function studies

Bruton's tyrosine kinase (Btk) is a cytoplasmic protein tyrosine kinase consisting of N-terminal pleckstrin homology (PH) domain followed by Tec homology (TH) domain, Src homology 3 and 2 (SH3 and SH2) domains, and a C-terminal kinase domain. Mutations in the human BTK gene cause the severe immunodeficiency disease X-linked agammaglobulinemia (XLA). The structural and functional basis of several XLA-causing mutations remains unknown, since only the structures of the PH and SH3 domains of human Btk are currently available. In this study, we overexpressed and purified a protein consisting of the SH3 and SH2 domains of human Btk for biochemical and structural analysis. The purified protein was only partially soluble and had a tendency to dimerize, which made it unsuitable for further studies. To overcome the problems of low solubility and dimerization, subdomain interactions were engineered without altering the function of the protein.     

Six X-linked agammaglobulinemia-causing missense mutations in the Src homology 2 domain of Bruton's tyrosine kinase: phosphotyrosine-binding and circular dichroism analysis

Src homology 2 (SH2) domains recognize phosphotyrosine (pY)-containing sequences and thereby mediate their association to ligands. Bruton's tyrosine kinase (Btk) is a cytoplasmic protein tyrosine kinase, in which mutations cause a hereditary immunodeficiency disease, X-linked agammaglobulinemia (XLA). Mutations have been found in all Btk domains, including SH2. We have analyzed the structural and functional effects of six disease-related amino acid substitutions in the SH2 domain: G302E, R307G, Y334S, L358F, Y361C, and H362Q. Also, we present a novel Btk SH2 missense mutation, H362R, leading to classical XLA. Based on circular dichroism analysis, the conformation of five of the XLA mutants studied differs from the native Btk SH2 domain, while mutant R307G is structurally identical. The binding of XLA mutation-containing SH2 domains to pY-Sepharose was reduced, varying between 1 and 13% of that for the native SH2 domain. The solubility of all the mutated proteins was remarkably reduced. SH2 domain mutations were divided into three categories: 1) Functional mutations, which affect residues presumably participating directly in pY binding (R307G); 2) structural mutations that, via conformational change, not only impair pY binding, but severely derange the structure of the SH2 domain and possibly interfere with the overall conformation of the Btk molecule (G302E, Y334S, L358F, and H362Q); and 3) structural-functional mutations, which contain features from both categories above (Y361C).     

Insights into the conformational flexibility of Bruton's tyrosine kinase from multiple ligand complex structures

Bruton's tyrosine kinase (BTK) plays a key role in B cell receptor signaling and is considered a promising drug target for lymphoma and inflammatory diseases. We have determined the X-ray crystal structures of BTK kinase domain in complex with six inhibitors from distinct chemical classes. Five different BTK protein conformations are stabilized by the bound inhibitors, providing insights into the structural flexibility of the Gly-rich loop, helix C, the DFG sequence, and activation loop. The conformational changes occur independent of activation loop phosphorylation and do not correlate with the structurally unchanged WEI motif in the Src homology 2-kinase domain linker. Two novel activation loop conformations and an atypical DFG conformation are observed representing unique inactive states of BTK. Two regions within the activation loop are shown to structurally transform between 3(10)- and α-helices, one of which collapses into the adenosine-5'-triphosphate binding pocket. The first crystal structure of a Tec kinase family member in the pharmacologically important DFG-out conformation and bound to a type II kinase inhibitor is described. The different protein conformations observed provide insights into the structural flexibility of BTK, the molecular basis of its regulation, and the structure-based design of specific inhibitors.     

Bruton tyrosine kinase-like protein, BtkSD, is present in the marine sponge Suberites domuncula

Sponges, the simplest and most ancient phylum of Metazoa, encode in their genome complex and highly sophisticated proteins that evolved together with multicellularity and are found only in metazoan animals. We report here the finding of a Bruton tyrosine kinase (BTK)-like protein in the marine sponge Suberites domuncula (Demospongiae). The nucleotide sequence of one sponge cDNA predicts a 700-aa-long protein, which contains all of the characteristic domains for the Tec family of protein tyrosine kinases (PTKs). The highest homology (38% identity, 55% overall similarity) was found with human BTK and TEC PTKs. Sponge PTK was therefore named BtkSD. Human BTK is involved in the maturation of B cells and mutations in the BTK gene cause X-linked agammaglobulinemia. Kinases from the Tec family are not present in Caenorhabditis elegans and, until now, they were found only in insects and higher animal taxa. Our finding implies that the BTK/TEC genes are of a very ancient origin.     

In silico analysis of Single Nucleotide Polymorphisms (SNPs) in human BRAF gene

BRAF gene mutations are frequently seen in both inherited and somatic diseases. However, the harmful mutations for BRAF gene have not been predicted in silico. Owing to the importance of BRAF gene in cell division, differentiation and secretion processes, the functional analysis was carried out to explore the possible association between genetic mutations and phenotypic variations. Genomic analysis of BRAF was initiated with SIFT followed by PolyPhen and SNPs&GO servers to retrieve the 85 deleterious non-synonymous SNPs (nsSNPs) from dbSNP. A total of 5 mutations i.e. c.406T>G (S136A), c.1446G>T (R462I), c.1556 A>G (K499E), c.1860 T>A (V600E) and c.2352 C>T (P764L) that are found to exert benign effects on the BRAF protein structure and function were chosen for further analysis. Protein structural analysis with these amino acid variants was performed by using I-Mutant, FOLD-X, HOPE, NetSurfP, Swiss PDB viewer, Chimera and NOMAD-Ref servers to check their solvent accessibility, molecular dynamics and energy minimization calculations. Our in silico analysis suggested that S136A and P764L variants of BRAF could directly or indirectly destabilize the amino acid interactions and hydrogen bond networks thus explain the functional deviations of protein to some extent. Screening for BRAF, S136A and P764Lvariants may be useful for disease molecular diagnosis and also to design the molecular inhibitors of BRAF pathways.     

SARS coronavirus E protein in phospholipid bilayers: an x-ray study

We investigated the structure of the hydrophobic domain of the severe acute respiratory syndrome E protein in model lipid membranes by x-ray reflectivity and x-ray scattering. In particular, we used x-ray reflectivity to study the location of an iodine-labeled residue within the lipid bilayer. The label imposes spatial constraints on the protein topology. Experimental data taken as a function of protein/lipid ratio P/L and different swelling states support the hairpin conformation of severe acute respiratory syndrome E protein reported previously. Changes in the bilayer thickness and acyl-chain ordering are presented as a function of P/L, and discussed in view of different structural models.     

Identification of the hydrophobic ligand binding pocket of the S1P1 receptor

Sphingosine 1-phosphate (S1P), a naturally occurring sphingolipid mediator and also a second messenger with growth factor-like actions in almost every cell type, is an endogenous ligand of five G protein-coupled receptors (GPCRs) in the endothelial differentiation gene family. The lack of GPCR crystal structures sets serious limitations to rational drug design and in silico searches for subtype-selective ligands. Here we report on the experimental validation of a computational model of the ligand binding pocket of the S1P1 GPCR surrounding the aliphatic portion of S1P. The extensive mutagenesis-based validation confirmed 18 residues lining the hydrophobic ligand binding pocket, which, combined with the previously validated three head group-interacting residues, now complete the mapping of the S1P ligand recognition site. We identified six mutants (L3.43G/L3.44G, L3.43E/L3.44E, L5.52A, F5.48G, V6.40L, and F6.44G) that maintained wild type [32P]S1P binding with abolished ligand-dependent activation by S1P. These data suggest a role for these amino acids in the conformational transition of S1P1 to its activated state. Three aromatic mutations (F5.48Y, F6.44G, and W6.48A) result in differential activation, by S1P or SEW2871, indicating that structural differences between the two agonists can partially compensate for differences in the amino acid side chain. The now validated ligand binding pocket provided us with a pharmacophore model, which was used for in silico screening of the NCI, National Institutes of Health, Developmental Therapeutics chemical library, leading to the identification of two novel nonlipid agonists of S1P1.     

BTKbase, mutation database for X-linked agammaglobulinemia (XLA).

X-linked agammaglobulinemia (XLA) is an immunodeficiency caused by mutations in the gene coding for Bruton's agammaglobulinemia tyrosine kinase (BTK). A database (BTKbase) of BTK mutations has been compiled and the recent update lists 225 entries from 189 unrelated families showing 148 unique molecular events. Each patient is given a unique patient identity number (PIN). Information is included regarding the phenotype including symptoms. Mutations in all the five domains of BTK have been noticed to cause the disease, the most common event being missense mutations. The mutations appear almost uniformly throughout the molecule and frequently affect CpG sites forming arginine residues. A decreased frequency of missense mutations was found in the TH, SH3 and upper lobe of the kinase domain. The putative structural implications of all the missense mutations are given in the database.     

TEC family kinases in health and disease--loss-of-function of BTK and ITK and the gain-of-function fusions ITK-SYK and BTK-SYK

The TEC family is ancient and constitutes the second largest family of cytoplasmic tyrosine kinases. In 1993, loss-of-function mutations in the BTK gene were reported as the cause of X-linked agammaglobulinemia. Of all the existing 90 tyrosine kinases in humans, Bruton's tyrosine kinase (BTK) is the kinase for which most mutations have been identified. These experiments of nature collectively provide a form of mutation scanning with direct implications for the several hundred endogenous signaling proteins carrying domains also found in BTK. In 2009, an inactivating mutation in the ITK gene was shown to cause susceptibility to lethal Epstein-Barr virus infection. Both kinases represent interesting targets for inhibition: in the case of BTK, as an immunosuppressant, whereas there is evidence that the inhibition of inducible T-cell kinase (ITK) could influence the infectivity of HIV and also have anti-inflammatory activity. Since 2006, several patients carrying a fusion protein, originating from a translocation joining genes encoding the kinases ITK and spleen tyrosine kinase (SYK), have been shown to develop T-cell lymphoma. We review these disease processes and also describe the role of the N-terminal pleckstrin homology-Tec homology (PH-TH) domain doublet of BTK and ITK in the downstream intracellular signaling of such fusion proteins.     

BTKbase, mutation database for X-linked agammaglobulinemia (XLA).

X-linked agammaglobulinemia (XLA) is an immunodeficiency caused by mutations in the gene coding for Bruton's agammaglobulinemia tyrosine kinase (BTK). A database (BTKbase) of BTK mutations has been compiled and the recent update lists 463 mutation entries from 406 unrelated families showing 303 unique molecular events. In addition to mutations, the database also lists variants or polymorphisms. Each patient is given a unique patient identity number (PIN). Information is included regarding the phenotype including symptoms. Mutations in all the five domains of BTK have been noticed to cause the disease, the most common event being missense mutations. The mutations appear almost uniformly throughout the molecule and frequently affect CpG sites that code for arginine residues. The putative structural implications of all the missense mutations are given in the database. The improved version of the registry having a number of new features is available at http://www. helsinki.fi/science/signal/btkbase.html     

A new mutation that predicted a drastic alteration of the BTK protein function

X-linked agammaglobulinemia (XLA) or Bruton disease is a relatively rare constitutionally immune disorder due to a genetic mutation of BTK (Bruton tyrosine kinase) gene which encodes for BTK protein. BTK is a signal-transducing protein expressed in hematopoietic lineages. The genetic disorder is responsible for B cell lymphocytes' maturation arrest. The humoral immunodeficiency caused by BTK mutation is linked with recurrent bacterial and viral infections. Genetic investigations of the prepositus as well as the other members of the family are necessary to characterize a mutation in BTK gene to confirm the diagnosis and reveal a hereditary transmission or de novo mutation.     

Identification of novel drug-resistant EGFR mutant inhibitors by in silico screening using comprehensive assessments of protein structures

EGFR is a target protein for the treatment of non small cell lung cancer (NSCLC). The mutations associated with the activation of EGFR kinase activity, such as L858R and G719S, destabilize the inactive conformation of EGFR and are closely linked with the development of NSCLC. The additional T790M mutation reportedly causes drug resistance against the commercially available EGFR inhibitors, gefitinib and erlotinib. In this study, we searched for novel G719S/T790M EGFR inhibitors by a new in silico screening strategy, using two datasets. The results of in silico screening using protein-ligand docking are affected by the selection of 3D structure of the target protein. As the first strategy, we chose the 3D structures for in silico screening by test dockings using the G719S/T790M crystal structure, its molecular dynamics snapshots, and known inhibitors of the drug-resistant EGFR. In the second strategy, we selected the 3D structures by test dockings using all of the EGFR structures, regardless of the mutations, and all of the known EGFR inhibitors. Using each of the 3D structures selected by the strategies, 1000 compounds were chosen from the 71,588 compounds. Kinase assays identified 15 G719S/T790M EGFR inhibitors, including two compounds with novel scaffolds. Analyses of their structure-activity relationships revealed that interactions with the mutated Met790 residue specifically increase the inhibitory activity against G719S/T790M EGFR.     

BTKbase, mutation database for X-linked agammaglobulinemia (XLA)

X-linked agammaglobulinemia (XLA) is an immunodeficiency caused by mutations in the gene coding for Bruton's agammaglobulinemia tyrosine kinase (BTK). A database (BTKbase) of BTK mutations has been compiled and the recent update lists 368 entries from 318 unrelated families showing 228 unique molecular events. In addition to mutations the database lists also some polymorphisms and site-directed mutations. Each patient is given a unique patient identity number (PIN). Information is provided regarding the phenotype including symptoms. Mutations in all the five domains of BTK have been noticed to cause the disease, the most common event being missense mutations. The mutations appear almost uniformly throughout the molecule and frequently affect CpG sites forming arginine residues. These hot spots have generally pyrimidines 5'and purines 3'to the mutated cytosine. A decreased frequency of missense mutations was found in the TH, SH3 and the upper lobe of the kinase domain. The putative structural implications of all the missense mutations are given in the database showing 228 unique molecular events, including a novel missense mutation causing an R28C substitution as previously seen in the Xid mouse.     

Role of the PHTH module in protein substrate recognition by Bruton's agammaglobulinemia tyrosine kinase

Defects in Bruton's tyrosine kinase (Btk) are responsible for X chromosome-linked agammaglobulinemia in patients. Mutations in each of the structural domains of Btk have been detected in patients, yet a mechanistic explanation for most of these mutant phenotypes is lacking. To understand the possible role of the unique pleckstrin homology and Tec homology (PHTH) module of Btk, we have compared the enzymatic properties of full-length Btk and a Btk mutant lacking the PHTH module (BtkDeltaPHTH). Here we show that Btk and BtkDeltaPHTH have similar basal catalytic activity but very different abilities to recognize protein substrates. Furthermore, the catalytic domain of Btk is inactive, in contrast to the catalytic domain of the prototypical Src tyrosine kinase that retains full catalytic ability. These data suggest that the PHTH module plays an important role in protein substrate recognition, that Btk and Src likely have different interdomain organizations and regulations, and that alterations in substrate recognition might play a role in X chromosome-linked agammaglobulinemia.     

Activating mutations of N-WASP alter Shigella pathogenesis

The pathogenesis of Shigella requires binding to the host protein N-WASP. To examine the roles of structural conformation and phospho-regulation of N-WASP during Shigella pathogenesis, mutant N-WASP constructs predicted to result in a constitutively open conformation (L229P and L232P) or either a phospho-mimicking (Y253E) or phospho-disruptive (Y253F) structure were constructed. Pyrene actin assays demonstrated that the N-WASP L229P and L232P constructs are constitutively active. Despite the increase in actin polymerization seen in vitro, cell lines expressing N-WASP L229P and L232P supported shorter actin tails when infected with Shigella.Shigella actin tails were unchanged in cells expressing N-WASP phospho-regulation mutant proteins. Shigella invasion, intracellular, and intercellular motility were not altered in cells expressing N-WASP L229P or L232P. However, plaque numbers were increased in cells expressing N-WASP L229P and L232P. These data demonstrate that N-WASP structural conformation is an important regulator of Shigella pathogenesis in distinct segments of its lifecycle.     

X-linked agammaglobulinemia (XLA): A genetic tyrosine kinase (Btk) disease

X-linked agammaglobulinemia is a heritable immunodeficiency disease caused by a differentiation abnormality, resulting in the virtual absence of B Iymphocytes and plasma cells. The affected gene encodes a cytoplasmic protein tyrosine kinase, Bruton's agammaglobulinemia tyrosine kinase, designated Btk. Btk and the other family members, Tec, Itk and Bmx, contain five regions, four of which are common structural and functional modules that are found in other signaling proteins. Mutations affect all domains of the gene, but amino acid substitutions seem to be confined to certain regions. More than 150 unique mutations have been identified and are collected in a mutation database, BTKbase. Here we discuss the three-dimensional structural implications of such mutations and their putative functional role. Of special interest are mutations affecting the pleckstrin homology domain, as Btk is the only disease-associated protein so far reported to carry mutations in this particular module.     

Novel adiponectin variants identified in type 2 diabetic patients reveal multimerization and secretion defects

ADIPOQ, encoding adiponectin, is a candidate gene for type 2 diabetes (T2D) identified by genome-wide linkage analyses with supporting evidence showing the protein function in sensitizing insulin actions. In an endeavor to characterize candidate genes causing T2D in Thai patients, we identified 10 novel ADIPOQ variations, several of which were non-synonymous variations observed only in the patients. To examine the impact of these non-synonymous variations on adiponectin structure and biochemical characteristics, we conducted a structural analysis of the wild-type and variant proteins by in silico modeling and further characterized biochemical properties of the variants with predicted structural abnormalities from the modeling by molecular and biochemical studies. The recombinant plasmids containing wild-type and variant ADIPOQ cDNAs derived from the variations identified by our study (R55H, R112H, and R131H) and previous work (G90S and R112C) were constructed and transiently expressed and co-expressed in cultured HEK293T cells to investigate their oligomerization, interaction, and secretion. We found that the novel R55H variant impaired protein multimerization but it did not exert the effect over the co-expressed wild-type protein while novel R131H variant impaired protein secretion and also affected the co-expressed wild-type protein in a dominant negative fashion. The R131H variant could traffic from the endoplasmic reticulum to the Golgi, trans-Golgi network, and early endosome but could not be secreted. The R131H variant was likely to be degraded through the lysosomal system and inhibition of its degradation rescued the variant protein from secretion defect. We have shown the possibility of using in silico modeling for predicting the effect of amino acid substitution on adiponectin oligomerization. This is also the first report that demonstrates a dominant negative effect of the R131H variant on protein secretion and the possibility of using protein degradation inhibitors as therapeutic agents in the patients carrying adiponectin variants with secretion defect.     

Structural effects of Parkinson's disease linked DJ-1 mutations

Mutations in the protein DJ-1 are associated with familial forms of Parkinson's disease, indicating that DJ-1 may be involved in pathways related to the etiology of this disorder. Here we have used solution state NMR and circular dichroism spectroscopies to evaluate the extent of structural perturbations associated with five different Parkinson's disease linked DJ-1mutations: L166P, E64D, M26I, A104T, and D149A. Comparison of the data with those obtained for the wild-type protein shows that the L166P mutation leads to severe and global destabilization and unfolding of the protein structure, while the structure of the E64D mutation, as expected, is nearly unperturbed. Interestingly, the remaining three mutants all show different degrees of structural perturbation, which are accompanied by a reduction in the thermodynamic stability of the protein. The observed structural and thermodynamic differences are likely to underlie any functional variations between these mutants and the wild type, which in turn are likely responsible for the pathogenicity of these mutations.     

Rational design and synthesis of a novel anti-leukemic agent targeting Bruton's tyrosine kinase (BTK), LFM-A13 [alpha-cyano-beta-hydroxy-beta-methyl-N-(2, 5-dibromophenyl)propenamide

In a systematic effort to design potent inhibitors of the anti-apoptotic tyrosine kinase BTK (Bruton's tyrosine kinase) as anti-leukemic agents with apoptosis-promoting and chemosensitizing properties, we have constructed a three-dimensional homology model of the BTK kinase domain. Our modeling studies revealed a distinct rectangular binding pocket near the hinge region of the BTK kinase domain with Leu460, Tyr476, Arg525, and Asp539 residues occupying the corners of the rectangle. The dimensions of this rectangle are approximately 18 x 8 x 9 x 17 A, and the thickness of the pocket is approximately 7 A. Advanced docking procedures were employed for the rational design of leflunomide metabolite (LFM) analogs with a high likelihood to bind favorably to the catalytic site within the kinase domain of BTK. The lead compound LFM-A13, for which we calculated a Ki value of 1.4 microM, inhibited human BTK in vitro with an IC50 value of 17.2 +/- 0.8 microM. Similarly, LFM-A13 inhibited recombinant BTK expressed in a baculovirus expression vector system with an IC50 value of 2.5 microM. The energetically favorable position of LFM-A13 in the binding pocket is such that its aromatic ring is close to Tyr476, and its substituent group is sandwiched between residues Arg525 and Asp539. In addition, LFM-A13 is capable of favorable hydrogen bonding interactions with BTK via Asp539 and Arg525 residues. Besides its remarkable potency in BTK kinase assays, LFM-A13 was also discovered to be a highly specific inhibitor of BTK. Even at concentrations as high as 100 micrograms/ml (approximately 278 microM), this novel inhibitor did not affect the enzymatic activity of other protein tyrosine kinases, including JAK1, JAK3, HCK, epidermal growth factor receptor kinase, and insulin receptor kinase. In accordance with the anti-apoptotic function of BTK, treatment of BTK+ B-lineage leukemic cells with LFM-A13 enhanced their sensitivity to ceramide- or vincristine-induced apoptosis. To our knowledge, LFM-A13 is the first BTK-specific tyrosine kinase inhibitor and the first anti-leukemic agent targeting BTK.