Structural basis and prediction of substrate specificity in protein serine/threonine kinases

The large number of protein kinases makes it impractical to determine their specificities and substrates experimentally. Using the available crystal structures, molecular modeling, and sequence analyses of kinases and substrates, we developed a set of rules governing the binding of a heptapeptide substrate motif (surrounding the phosphorylation site) to the kinase and implemented these rules in a web-interfaced program for automated prediction of optimal substrate peptides, taking only the amino acid sequence of a protein kinase as input. We show the utility of the method by analyzing yeast cell cycle control and DNA damage checkpoint pathways. Our method is the only available predictive method generally applicable for identifying possible substrate proteins for protein serine/threonine kinases and helps in silico construction of signaling pathways. The accuracy of prediction is comparable to the accuracy of data from systematic large-scale experimental approaches.

Assessment of arginine 97 and lysine 72 as determinants of substrate specificity in cytochrome P450 2C9 (CYP2C9)

CYP2C9 is distinguished by a preference for substrates bearing a negative charge at physiological pH. Previous studies have suggested that CYP2C9 residues R97 and K72 may play roles in determining preference for anionic substrates by interaction at the active site or in the access channel. The aim of the present study was to assess the role of these two residues in determining substrate selectivity. R97 and K72 were substituted with negative, uncharged polar and hydrophobic residues using a degenerate polymerase chain reaction-directed strategy. Wild-type and mutant enzymes were expressed in bicistronic format with human cytochrome P450 reductase in Escherichia coli. Mutation of R97 led to a loss of holoenzyme expression for R97A, R97V, R97L, R97T, and R97E mutants. Low levels of hemoprotein were detected for R97Q, R97K, R97I, and R97P mutants. Significant apoenzyme was observed, suggesting that heme insertion or protein stability was compromised in R97 mutants. These observations are consistent with a structural role for R97 in addition to any role in substrate binding. By contrast, all K72 mutants examined (K72E, K72Q, K72V, and K72L) could be expressed as hemoprotein at levels comparable to wild-type. Type I binding spectra were obtained with wildtype and K72 mutants using diclofenac and ibuprofen. Mutation of K72 had little or no effect on the interaction with these substrates, arguing against a critical role in determining substrate specificity. Thus, neither residue appears to play a role in determining substrate specificity, but a structural role for R97 can be proposed consistent with recently published crystallographic data for CYP2C9 and CYP2C5.

Structural basis for the specificity of bipartite nuclear localization sequence binding by importin-α

Importin-alpha is the nuclear import receptor that recognizes cargo proteins carrying conventional basic monopartite and bipartite nuclear localization sequences (NLSs) and facilitates their transport into the nucleus. Bipartite NLSs contain two clusters of basic residues, connected by linkers of variable lengths. To determine the structural basis of the recognition of diverse bipartite NLSs by mammalian importin-alpha, we co-crystallized a non-autoinhibited mouse receptor protein with peptides corresponding to the NLSs from human retinoblastoma protein and Xenopus laevis phosphoprotein N1N2, containing diverse sequences and lengths of the linker. We show that the basic clusters interact analogously in both NLSs, but the linker sequences adopt different conformations, whereas both make specific contacts with the receptor. The available data allow us to draw general conclusions about the specificity of NLS binding by importin-alpha and facilitate an improved definition of the consensus sequence of a conventional basic/bipartite NLS (KRX10-12KRRK) that can be used to identify novel nuclear proteins.

Substrate specifity of protein kinases and computational prediction of substrates

To ensure signalling fidelity, kinases must act only on a defined subset of cellular targets. Appreciating the basis for this substrate specificity is essential for understanding the role of an individual protein kinase in a particular cellular process. The specificity in the cell is determined by a combination of peptide specificity of the kinase (the molecular recognition of the sequence surrounding the phosphorylation site), substrate recruitment and phosphatase activity. Peptide specificity plays a crucial role and depends on the complementarity between the kinase and the substrate and therefore on their three-dimensional structures. Methods for experimental identification of kinase substrates and characterization of specificity are expensive and laborious, therefore, computational approaches are being developed to reduce the amount of experimental work required in substrate identification. We discuss the structural basis of substrate specificity of protein kinases and review the experimental and computational methods used to obtain specificity information. (c) 2005 Elsevier B.V. All rights reserved.

Insights to substrate binding and processing by West Nile Virus NS3 protease through combined modeling, protease mutagenesis, and kinetic studies

West Nile Virus is becoming a widespread pathogen, infecting people on at least four continents with no effective treatment for these infections or many of their associated pathologies. A key enzyme that is essential for viral replication is the viral protease NS2B-NS3, which is highly conserved among all flaviviruses. Using a combination of molecular fitting of substrates to the active site of the crystal structure of NS3,site-directed enzyme and cofactor mutagenesis, and kinetic studies on proteolytic processing of panels of short peptide substrates, we have identified important enzyme-substrate interactions that define substrate specificity for NS3 protease. In addition to better understanding the involvement of S2, S3, and S4 enzyme residues in substrate binding, a residue within cofactor NS2B has been found to strongly influence the preference of flavivirus proteases for lysine or arginine at P2 in substrates. Optimization of tetrapeptide substrates for enhanced protease affinity and processing efficiency has also provided important clues for developing inhibitors of West Nile Virus infection.

Enzymatic cyclization of a potent Bowman-Birk protease inhibitor, sunflower trypsin inhibitor-1, and solution structure of an acyclic precursor peptide

The most potent known naturally occurring Bowman-Birk inhibitor, sunflower trypsin inhibitor-1 (SFTI-1), is a bicyclic 14-amino acid peptide from sunflower seeds comprising one disulfide bond and a cyclic backbone. At present, little is known about the cyclization mechanism of SFTI-1. We show here that an acyclic permutant of SFTI-1 open at its scissile bond, SFTI-1[ 6,5], also functions as an inhibitor of trypsin and that it can be enzymatically backbone-cyclized by incubation with bovine beta-trypsin. The resulting ratio of cyclic SFTI-1 to SFTI1[6,5] is similar to9:1 regardless of whether trypsin is incubated with SFTI-1[ 6,5] or SFTI-1. Enzymatic resynthesis of the scissile bond to form cyclic SFTI-1 is a novel mechanism of cyclization of SFTI-1[ 6,5]. Such a reaction could potentially occur on a trypsin affinity column as used in the original isolation procedure of SFTI-1. We therefore extracted SFTI-1 from sunflower seeds without a trypsin purification step and confirmed that the backbone of SFTI-1 is indeed naturally cyclic. Structural studies on SFTI-1[ 6,5] revealed high heterogeneity, and multiple species of SFTI-1[ 6,5] were identified. The main species closely resembles the structure of cyclic SFTI-1 with the broken binding loop able to rotate between a cis/trans geometry of the I7-P8 bond with the cis conformer being similar to the canonical binding loop conformation. The non-reactive loop adopts a beta-hairpin structure as in cyclic wild-type SFTI-1. Another species exhibits an isoaspartate residue at position 14 and provides implications for possible in vivo cyclization mechanisms.

Human tartrate-resistant acid phosphatase becomes an effective ATPase upon proteolytic activation

Proteolytic, cleavage in an exposed loop of human tartrate-resistant acid phosphatase (TRAcP) with trypsin leads to a significant increase in activity. At each pH value between 3.25 and 8.0 the cleaved enzyme is more active. Substrate specificity is also influenced by proteolysis. Only the cleaved form is able to hydrolyze unactivated substrates efficiently, and at pH > 6 cleaved TRAcP acquires a marked preference for ATP. The cleaved enzyme also has altered sensitivity to inhibitors. Interestingly, the magnitude and mode of inhibition by fluoride depends not only on the proteolytic state but also pH. The combined kinetic data imply a role of the loop residue D158 in catalysis in the cleaved enzyme. Notably, at low pH this residue may act as a proton donor for the leaving group. In this respect the mechanism of cleaved TRAcP resembles that of sweet potato purple acid phosphatase. (c) 2005 Elsevier Inc. Ail rights reserved.

Protein kinases associated with the yeast phosphoproteome

Background: Protein phosphorylation is an extremely important mechanism of cellular regulation. A large-scale study of phosphoproteins in a whole-cell lysate of Saccharomyces cerevisiae has previously identified 383 phosphorylation sites in 216 peptide sequences. However, the protein kinases responsible for the phosphorylation of the identified proteins have not previously been assigned. Results: We used Predikin in combination with other bioinformatic tools, to predict which of 116 unique protein kinases in yeast phosphorylates each experimentally determined site in the phosphoproteome. The prediction was based on the match between the phosphorylated 7-residue sequence and the predicted substrate specificity of each kinase, with the highest weight applied to the residues or positions that contribute most to the substrate specificity. We estimated the reliability of the predictions by performing a parallel prediction on phosphopeptides for which the kinase has been experimentally determined. Conclusion: The results reveal that the functions of the protein kinases and their predicted phosphoprotein substrates are often correlated, for example in endocytosis, cytokinesis, transcription, replication, carbohydrate metabolism and stress response. The predictions link phosphoproteins of unknown function with protein kinases with known functions and vice versa, suggesting functions for the uncharacterized proteins. The study indicates that the phosphoproteins and the associated protein kinases represented in our dataset have housekeeping cellular roles; certain kinases are not represented because they may only be activated during specific cellular responses. Our results demonstrate the utility of our previously reported protein kinase substrate prediction approach (Predikin) as a tool for establishing links between kinases and phosphoproteins that can subsequently be tested experimentally.

Molecular modelling of human CYP1B1 substrate interactions and investigation of allelic variant effects on metabolism

Molecular modelling of human CYP1B1 based on homology with the mammalian P450, CYP2C5, of known three-dimensional structure is reported. The enzyme model has been used to investigate the likely mode of binding for selected CYP1B1 substrates, particularly with regard to the possible effects of allelic variants of CYP1B1 on metabolism. In general, it appears that the CYP1B1 model is consistent with known substrate selectivity for the enzyme, and the sites of metabolism can be rationalized in terms of specific contacts with key amino acid residues within the CYP1B1 heme locus. Further-more, a mode of binding interaction for the inhibitor, a-naphthoflavone, is presented which accords with currently available information. The current paper shows that a combination of molecular modelling and experimental determinations on the substrate metabolism for CYP1B1 allelic variants can aid in the understanding of structure-function relationships within P450 enzymes. (C) 2003 Elsevier Science Ireland Ltd. All rights reserved.

Role of glutamic acid 216 in cytochrome P450 2D6 substrate binding and catalysis

Human cytochrome P450 (P450) 2D6 is an important enzyme involved in the metabolism of drugs, many of which are amines or contain other basic nitrogen atoms. Asp301 has generally been considered to be involved in electrostatic docking with the basic substrates, on the basis of previous modeling studies and site-directed mutagenesis. Substitution of Glu216 with a residue other than Asp strongly attenuated the binding of quinidine, bufuralol, and several other P450 2D6 ligands. Catalytic activity with the substrates bufuralol and 4-methoxyphenethylamine was strongly inhibited by neutral or basic mutations at Glu216 (>95%), to the same extent as the substitution of Asn at Asp301. Unlike the Asp301 mutants, the Gln216 mutant (E216Q) retained 40% enzyme efficiency with the substrate spirosulfonamide, devoid of basic nitrogen, suggesting that the substitutions at Glu216 affect binding of amine substrates more than other catalytic steps. Attempts to induce catalytic specificity toward new substrates by substitutions at Asp301 and Glu216 were unsuccessful. Collectively, the results provide evidence for electrostatic interaction of amine substrates with Glu216, and we propose that both of these acidic residues plus at least another residue(s) is (are) involved in binding the repertoire of P450 2D6 ligands.

Active site mutations and substrate inhibition in human sulfotransferase 1A1 and 1A3

Human SULT1A1 is primarily responsible for sulfonation of xenobiotics, including the activation of promutagens, and it has been implicated in several forms of cancer. Human SULT1A3 has been shown to be the major sulfotransferase that sulfonates dopamine. These two enzymes shares 93% amino acid sequence identity and have distinct but overlapping substrate preferences. The resolution of the crystal structures of these two enzymes has enabled us to elucidate the mechanisms controlling their substrate preferences and inhibition. The presence of two p-nitrophenol (pNP) molecules in the crystal structure of SULT1A1 was postulated to explain cooperativity at low and inhibition at high substrate concentrations, respectively. In SULT1A1, substrate inhibition occurs with pNP as the substrate but not with dopamine. For SULT1A3, substrate inhibition is found for dopamine but not with pNP. We investigated how substrate inhibition occurs in these two enzymes using molecular modeling, site-directed mutagenesis, and kinetic analysis. The results show that residue Phe-247 of SULT1A1, which interacts with both p-nitrophenol molecules in the active site, is important for substrate inhibition. Mutation of phenylalanine to leucine at this position in SULT1A1 results in substrate inhibition by dopamine. We also propose, based on modeling and kinetic studies, that substrate inhibition by dopamine in SULT1A3 is caused by binding of two dopamine molecules in the active site. © 2004 by The American Society for Biochemistry and Molecular Biology, Inc.

Elucidating the specificity of binding of sulfonylurea herbicides to acetohydroxyacid synthase

Acetohydroxyacid synthase (AHAS, EC 2.2.1.6) is the target for the sulfonylurea herbicides, which act as potent inhibitors of the enzyme. Chlorsulfuron (marketed as Glean) and sulforneturon methyl (marketed as Oust) are two commercially important members of this family of herbicides. Here we report crystal structures of yeast AHAS in complex with chlorsulfuron (at a resolution of 2.19 Angstrom), sulforneturon methyl (2.34 Angstrom), and two other sulfonylureas, metsulfuron methyl (2.29 Angstrom) and tribenuron methyl (2.58 Angstrom). The structures observed suggest why these inhibitors have different potencies and provide clues about the differential effects of mutations in the active site tunnel on various inhibitors. In all of the structures, the thiamin diphosphate cofactor is fragmented, possibly as the result of inhibitor binding. In addition to thiamin diphosphate, AHAS requires FAD for activity. Recently, it has been reported that reduction of FAD can occur as a minor side reaction due to reaction with the carbanion/enamine of the hydroxyethyl-ThDP intermediate that is formed midway through the catalytic cycle. Here we report that the isoalloxazine ring has a bent conformation that would account for its ability to accept electrons from the hydroxyethyl intermediate. Most sequence and mutation data suggest that yeast AHAS is a high-quality model for the plant enzyme.

Mesostigmatid mites (Acari : Mesostigmata) on rainforest tree trunks: Arboreal specialists, but substrate generalists?

Predatory mites (Acari: Mesostigmata) on tree trunks without significant epiphytic growth in a subtropical rainforest in Eastern Australia were assessed for habitat specificity (i.e. whether they are tree trunk specialists or occupying other habitats) and the influence of host tree and bark structure on their abundance, species richness and species composition. The trunks of nine tree species from eight plant families representing smooth, intermediate and rough bark textures were sampled using a knockdown insecticide spray. In total, 12 species or morphospecies of Mesostigmata (excluding Uropodina sensu stricto) were collected, most of which are undescribed. Comparison with collections from other habitats indicates that epicorticolous Mesostigmata are mainly represented by suspended soil dwellers (six species), secondarily by generalists (four species) and a bark specialist (one species). A typical ground-dwelling species was also found but was represented only by a single individual. In terms of abundance, 50.5% of individuals were suspended soil dwellers, 40.7% bark specialists, and 8.3% generalists. Host species and bark roughness had no significant effect on abundance or species richness. Furthermore, there was no clear effect on species composition. The distribution of the most frequently encountered species suggests that most mesostigmatid mites living on bark use many or most rainforest tree species, independent of bark roughness. These findings support the hypothesis that some epicorticolous Mesostigmata use tree trunks as 'highways' for dispersing between habitat patches, while others use it as a permanent habitat.