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.

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.

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.

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.

Binding of an RNA trafficking response element to heterogeneous nuclear ribonucleoproteins A1 and A2

Heterogeneous nuclear ribonucleoprotein (hnRNP) A2 binds a 21-nucleotide myelin basic protein mRNA response element, the A2RE, and A2RE-like sequences in other localized mRNAs, and is a trans-acting factor in oligodendrocyte cytoplasmic RNA trafficking. Recombinant human hnRNPs A1 and A2 were used in a biosensor to explore interactions with A2RE and the cognate oligodeoxyribonucleotide. Both proteins have a single site that bound oligonucleotides with markedly different sequences but did not bind in the presence of heparin. Both also possess a second, specific site that bound only A2RE and was unaffected by heparin, hnRNP A2 bound A2RE in the latter site with a K-d near 50 nM, whereas the K-d for hnRNP A1 was above 10 muM. UV cross-linking assays led to a similar conclusion. Mutant A2RE sequences, that in earlier qualitative studies appeared not to bind hnRNP A2 or support RNA trafficking in oligodendrocytes, had dissociation constants above 5 muM for this protein. The two concatenated RNA recognition motifs (RRMs), but not the individual RRMs, mimicked the binding behavior of hnRNP A2. These data highlight the specificity of the interaction of A2RE with these hnRNPs and suggest that the sequence-specific A2RE-binding site on hnRNP A2 is formed by both RRMs acting in cis.

Characterization of metal ion-induced [3H]inositol hexakisphosphate binding to rat cerebellar membranes

The binding of [3H]inositol hexakisphosphate ([3H] InsP6) to rat cerebellar membranes has been characterized with the objective of establishing the role, if any, of a membrane protein receptor. In the presence of EDTA, we have previously identified an InsP6-binding site with a capacity of approximately 20 pmol/mg protein (Hawkins, P. T., Reynolds, D. J. M., Poyner, D. R., and Hanley, M. R. (1990) Biochem. Biophys. Res. Commun. 167, 819-827). However, in the presence of 1 mM Mg2+, the capacity of [3H]InsP6 binding to membranes was increased approximately 9-fold. This enhancing effect of Mg2+ was reversed by addition of 10 microM of several cation chelators, suggesting that the increased binding required trace quantities of other metal cations. This is supported by experiments where it was possible to saturate binding by addition of excess membranes, despite not significantly depleting radioligand, pointing to removal of some other factor. Removal of endogenous cations from the binding assay by pretreatment with chelex resin also prevents the Mg(2+)-induced potentiation. Consideration of the specificity of the chelators able to abolish this potentiation suggested involvement of Fe3+ or Al3+. Both these ions (but not several others) were able to increase [3H]InsP6 binding to chelex-pretreated membranes at concentrations of 1 microM. It is possible to demonstrate synergy between Fe3+ and Mg2+ under these conditions. We propose that [3H]InsP6 may interact with membranes through non-protein recognition possibly via phospholipids, in a manner dependent upon trace metals. The implications of this for InsP6 biology are considered.

Identification of a novel inositol phosphate recognition site: specific [3H]inositol hexakisphosphate binding to brain regions and cerebellar membranes

[3H]Inositol hexakisphosphate (InsP6) binds with a heterogeneous distribution to frozen sections of unfixed rat brain and is displaced by unlabelled InsP6. The pattern of binding correlates with binding to neuronal cell bodies. [3H]InsP6 binding to cerebellar membranes has been further characterised, is reversible, and saturable, and exhibits high specificity for inositol polyphosphates. The IC50 for competition by unlabelled InsP6 is approximately 100nM, whereas inositol 1,3,4,5,6 pentakisphosphate (Ins(13456)P5), inositol 1,3,4,5 tetrakisphosphate (Ins(1345)P4), and inositol 1,4,5 trisphosphate (Ins(145)P3) bind with an affinity at least one order of magnitude lower. [3H]InsP6 binding is clearly distinct from previously characterised Ins(145)P3 (ref. 1, 2) and Ins(1345)P4 (ref. 3) binding, both in terms of pharmacology and brain distribution.

Oxidative modification of a specific apolipoprotein B lysine residue confers altered receptor specificity on LDL

During inflammation, many cell types release reactive oxygen species (ROS) via the respiratory burst. These ROS are potent oxidants of LDL and its major protein, apolipoprotein B. Whilst native LDL is taken up by endothelial cells via a feedback controlled receptor-regulated process, oxidative modification of LDL renders it a ligand for many scavenger receptors. Scavenger receptors include CD-36, LOX-1 and the prototypic macrophage SR A I/II, all of which are variably expressed. Uncontrolled uptake of oxidised LDL is implicated in the pathogenesis of atherosclerosis. In addition, oxidised LDL increases CCR2 protein and mRNA expression on monocytes, and thus may contribute to monocyte retention and perpetuation in inflammatory, unstable atherosclerotic lesions. However, little data are available on the effects of specific minor modifications to apolipoprotein B. In order to identify the sequence specificity and nature of oxidative modifications which confer altered properties on LDL, we have investigated the effects of modified peptides (which correspond to the putative LDLR binding domain) on LDL uptake by HUVECs and U937 monocytes.

Discovering properties of new DNA-binding activity of proteins

Protein-DNA interactions are an essential feature in the genetic activities of life, and the ability to predict and manipulate such interactions has applications in a wide range of fields. This Thesis presents the methods of modelling the properties of protein-DNA interactions. In particular, it investigates the methods of visualising and predicting the specificity of DNA-binding Cys2His2 zinc finger interaction. The Cys2His2 zinc finger proteins interact via their individual fingers to base pair subsites on the target DNA. Four key residue positions on the a- helix of the zinc fingers make non-covalent interactions with the DNA with sequence specificity. Mutating these key residues generates combinatorial possibilities that could potentially bind to any DNA segment of interest. Many attempts have been made to predict the binding interaction using structural and chemical information, but with only limited success. The most important contribution of the thesis is that the developed model allows for the binding properties of a given protein-DNA binding to be visualised in relation to other protein-DNA combinations without having to explicitly physically model the specific protein molecule and specific DNA sequence. To prove this, various databases were generated, including a synthetic database which includes all possible combinations of the DNA-binding Cys2His2 zinc finger interactions. NeuroScale, a topographic visualisation technique, is exploited to represent the geometric structures of the protein-DNA interactions by measuring dissimilarity between the data points. In order to verify the effect of visualisation on understanding the binding properties of the DNA-binding Cys2His2 zinc finger interaction, various prediction models are constructed by using both the high dimensional original data and the represented data in low dimensional feature space. Finally, novel data sets are studied through the selected visualisation models based on the experimental DNA-zinc finger protein database. The result of the NeuroScale projection shows that different dissimilarity representations give distinctive structural groupings, but clustering in biologically-interesting ways. This method can be used to forecast the physiochemical properties of the novel proteins which may be beneficial for therapeutic purposes involving genome targeting in general.

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background - MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions. Results - A large dataset comprising MHC-peptide structural complexes was created by re-modelling pre-determined x-ray crystallographic structures. Static energetic analysis, following energy minimisation, was performed on the dataset in order to characterise interactions between bound peptides and the MHC Class I molecule, partitioning the interactions within the groove into van der Waals, electrostatic and total non-bonded energy contributions. Conclusion - The QSAR techniques of Genetic Function Approximation (GFA) and Genetic Partial Least Squares (G/PLS) algorithms were used to identify key interactions between the two molecules by comparing the calculated energy values with experimentally-determined BL50 data. Although the peptide termini binding interactions help ensure the stability of the MHC Class I-peptide complex, the central region of the peptide is also important in defining the specificity of the interaction. As thermodynamic studies indicate that peptide association and dissociation may be driven entropically, it may be necessary to incorporate entropic contributions into future calculations.

Optimisation of a droplet digital PCR assay for the diagnosis of Schistosoma japonicum infection: A duplex approach with DNA binding dye chemistry

Schistosomiasis is a chronically debilitating helminth infection with a significant socio-economic and public health impact. Accurate diagnostics play a pivotal role in achieving current schistosomiasis control and elimination goals. However, many of the current diagnostic procedures, which rely on detection of schistosome eggs, have major limitations including lack of accuracy and the inability to detect pre-patent infections. DNA-based detection methods provide a viable alternative to the current tests commonly used for schistosomiasis diagnosis. Here we describe the optimisation of a novel droplet digital PCR (ddPCR) duplex assay for the diagnosis of Schistosoma japonicum infection which provides improved detection sensitivity and specificity. The assay involves the amplification of two specific and abundant target gene sequences in S. japonicum; a retrotransposon (SjR2) and a portion of a mitochondrial gene (nad1). The assay detected target sequences in different sources of schistosome DNA isolated from adult worms, schistosomules and eggs, and exhibits a high level of specificity, thereby representing an ideal tool for the detection of low levels of parasite DNA in different clinical samples including parasite cell free DNA in the host circulation and other bodily fluids. Moreover, being quantitative, the assay can be used to determine parasite infection intensity and, could provide an important tool for the detection of low intensity infections in low prevalence schistosomiasis-endemic areas.

Receptor activity-modifying protein-directed G protein signaling specificity for the calcitonin gene-related peptide family of receptors

The calcitonin gene-related peptide (CGRP) family of G protein- coupled receptors (GPCRs) is formed through the association of the calcitonin receptor-like receptor (CLR) and one of three receptor activity-modifying proteins (RAMPs). Binding of one of the three peptide ligands, CGRP, adrenomedullin (AM), and intermedin/adrenomedullin 2 (AM2), is well known to result in aGαs-mediated increase in cAMP. Here we used modified yeast strains that couple receptor activation to cell growth, via chimeric yeast/Gα subunits, and HEK-293 cells to characterize the effect of different RAMP and ligand combinations on this pathway. We not only demonstrate functional couplings to both Gαs and Gαq but also identify a Gαi component to CLR signaling in both yeast and HEK-293 cells, which is absent in HEK-293S cells. We show that the CGRP family of receptors displays both ligand- and RAMPdependent signaling bias among the Gαs, Gαi, and Gαq/11 pathways. The results are discussed in the context of RAMP interactions probed through molecular modeling and molecular dynamics simulations of the RAMP-GPCR-G protein complexes. This study further highlights the importance of RAMPs to CLR pharmacology and to bias in general, as well as identifying the importance of choosing an appropriate model system for the study of GPCR pharmacology.

Comparison between the 1% and 2% ibopamine provocative test in primary open-angle glaucoma patients: sensitivity, specificity and tolerability

PURPOSE: To compare intraocular pressure (IOP) rise in normal individuals and primary open-angle glaucoma patients and the safety and efficacy of ibopamine eye drops in different concentrations as a provocative test for glaucoma. METHODS: Glaucoma patients underwent (same eye) the ibopamine provocative test with two concentrations, 1% and 2%, in a random sequence at least 3 weeks apart, but not more than 3 months. The normal individuals were randomly submitted to one of the concentrations of ibopamine (1% and 2%). The test was considered positive if there was an IOP rise greater than 3 or 4 mmHg at 30 or 45 minutes to test which subset of the test has the best sensitivity (Se)/specificity (Sp). RESULTS: There was no statistically significant difference in any of the IOP measurements, comparing 1% with 2% ibopamine. The IOP was significantly higher at 30 and 45 minutes with both concentrations (p<0.001). The best sensitivity/specificity ratio was achieved with the cutoff point set as greater than 3 mmHg at 45 minutes with 2% ibopamine (area under the ROC curve: 0.864, Se: 84.6%; Sp:73.3%). All patients described a slight burning after ibopamine's instillation. CONCLUSION: 2% ibopamine is recommended as a provocative test for glaucoma. Because both concentrations have similar ability to rise IOP, 1% ibopamine may be used to treat ocular hypotony.

Sensitivity and specificity of machine learning classifiers for glaucoma diagnosis using Spectral Domain OCT and standard automated perimetry

PURPOSE: To evaluate the sensitivity and specificity of machine learning classifiers (MLCs) for glaucoma diagnosis using Spectral Domain OCT (SD-OCT) and standard automated perimetry (SAP). METHODS: Observational cross-sectional study. Sixty two glaucoma patients and 48 healthy individuals were included. All patients underwent a complete ophthalmologic examination, achromatic standard automated perimetry (SAP) and retinal nerve fiber layer (RNFL) imaging with SD-OCT (Cirrus HD-OCT; Carl Zeiss Meditec Inc., Dublin, California). Receiver operating characteristic (ROC) curves were obtained for all SD-OCT parameters and global indices of SAP. Subsequently, the following MLCs were tested using parameters from the SD-OCT and SAP: Bagging (BAG), Naive-Bayes (NB), Multilayer Perceptron (MLP), Radial Basis Function (RBF), Random Forest (RAN), Ensemble Selection (ENS), Classification Tree (CTREE), Ada Boost M1(ADA),Support Vector Machine Linear (SVML) and Support Vector Machine Gaussian (SVMG). Areas under the receiver operating characteristic curves (aROC) obtained for isolated SAP and OCT parameters were compared with MLCs using OCT+SAP data. RESULTS: Combining OCT and SAP data, MLCs' aROCs varied from 0.777(CTREE) to 0.946 (RAN).The best OCT+SAP aROC obtained with RAN (0.946) was significantly larger the best single OCT parameter (p<0.05), but was not significantly different from the aROC obtained with the best single SAP parameter (p=0.19). CONCLUSION: Machine learning classifiers trained on OCT and SAP data can successfully discriminate between healthy and glaucomatous eyes. The combination of OCT and SAP measurements improved the diagnostic accuracy compared with OCT data alone.

Essential regulatory elements for NHE3 gene transcription in renal proximal tubule cells

The objectives of the present study were to identify the cis-elements of the promoter absolutely required for the efficient rat NHE3 gene transcription and to locate positive and negative regulatory elements in the 5’-flanking sequence (5’FS), which might modulate the gene expression in proximal tubules, and to compare this result to those reported for intestinal cell lines. We analyzed the promoter activity of different 5’FS segments of the rat NHE3 gene, in the OKP renal proximal tubule cell line by measuring the activity of the reporter gene luciferase. Because the segment spanning the first 157 bp of 5’FS was the most active it was studied in more detail by sequential deletions, point mutations, and gel shift assays. The essential elements for gene transcription are in the region -85 to -33, where we can identify consensual binding sites for Sp1 and EGR-1, which are relevant to NHE3 gene basal transcription. Although a low level of transcription is still possible when the first 25 bp of the 5’FS are used as promoter, efficient transcription only occurs with 44 bp of 5’FS. There are negative regulatory elements in the segments spanning -1196 to -889 and -467 to -152, and positive enhancers between -889 and -479 bp of 5’FS. Transcription factors in the OKP cell nuclear extract efficiently bound to DNA elements of rat NHE3 promoter as demonstrated by gel shift assays, suggesting a high level of similarity between transcription factors of both species, including Sp1 and EGR-1.