In silico identification and characterization of protein-ligand binding sites

Protein–ligand binding site prediction methods aim to predict, from amino acid sequence, protein–ligand interactions, putative ligands, and ligand binding site residues using either sequence information, structural information, or a combination of both. In silico characterization of protein–ligand interactions has become extremely important to help determine a protein’s functionality, as in vivo-based functional elucidation is unable to keep pace with the current growth of sequence databases. Additionally, in vitro biochemical functional elucidation is time-consuming, costly, and may not be feasible for large-scale analysis, such as drug discovery. Thus, in silico prediction of protein–ligand interactions must be utilized to aid in functional elucidation. Here, we briefly discuss protein function prediction, prediction of protein–ligand interactions, the Critical Assessment of Techniques for Protein Structure Prediction (CASP) and the Continuous Automated EvaluatiOn (CAMEO) competitions, along with their role in shaping the field. We also discuss, in detail, our cutting-edge web-server method, FunFOLD for the structurally informed prediction of protein–ligand interactions. Furthermore, we provide a step-by-step guide on using the FunFOLD web server and FunFOLD3 downloadable application, along with some real world examples, where the FunFOLD methods have been used to aid functional elucidation.

Analysis of protein networks in resting and collagen receptor (GPVI)-stimulated platelet sub-proteomes.

Proteomics approaches have made important contributions to the characterisation of platelet regulatory mechanisms. A common problem encountered with this method, however, is the masking of low-abundance (e.g. signalling) proteins in complex mixtures by highly abundant proteins. In this study, subcellular fractionation of washed human platelets either inactivated or stimulated with the glycoprotein (GP) VI collagen receptor agonist, collagen-related peptide, reduced the complexity of the platelet proteome. The majority of proteins identified by tandem mass spectrometry are involved in signalling. The effect of GPVI stimulation on levels of specific proteins in subcellular compartments was compared and analysed using in silico quantification, and protein associations were predicted using STRING (the search tool for recurring instances of neighbouring genes/proteins). Interestingly, we observed that some proteins that were previously unidentified in platelets including teneurin-1 and Van Gogh-like protein 1, translocated to the membrane upon GPVI stimulation. Newly identified proteins may be involved in GPVI signalling nodes of importance for haemostasis and thrombosis.

Activation of protein kinase cascades in the heart by hypertrophic G protein-coupled receptor agonists.

Cardiac myocyte hypertrophy involves changes in cell structure and alterations in protein expression regulated at both the transcriptional and translational levels. Hypertrophic G protein-coupled receptor (GPCR) agonists such as endothelin-(ET-1) and phenylephrine stimulate a number of protein kinase cascades in the heart. Mitogen-activated protein kinase (MAPK) cascades stimulated include the extracellularly regulated kinase cascade, the stress-activated protein kinase/c-Jun N-terminal kinase cascade, and the p38 MAPK cascade. All 3 pathways have been implicated in hypertrophy, but recent ex vivo evidence also suggests that there may be additional effects on cell survival. ET-1 and phenylephrine also stimulate the protein kinase B pathway, and this may be involved in the regulation of protein synthesis by these agonists. Thus, protein kinase-mediated signaling may be important in the regulation of the development of myocyte hypertrophy.

Regulation of protein kinase B and 4E-BP1 by oxidative stress in cardiac myocytes.

Stimulation of phosphatidylinositol 3'-kinase (PI3K) and protein kinase B (PKB) is implicated in the regulation of protein synthesis in various cells. One mechanism involves PI3K/PKB-dependent phosphorylation of 4E-BP1, which dissociates from eIF4E, allowing initiation of translation from the 7-methylGTP cap of mRNAs. We examined the effects of insulin and H(2)O(2) on this pathway in neonatal cardiac myocytes. Cardiac myocyte protein synthesis was increased by insulin, but was inhibited by H(2)O(2). PI3K inhibitors attenuated basal levels of protein synthesis and inhibited the insulin-induced increase in protein synthesis. Insulin or H(2)O(2) increased the phosphorylation (activation) of PKB through PI3K, but, whereas insulin induced a sustained response, the response to H(2)O(2) was transient. 4E-BP1 was phosphorylated in unstimulated cells, and 4E-BP1 phosphorylation was increased by insulin. H(2)O(2) stimulated dephosphorylation of 4E-BP1 by increasing protein phosphatase (PP1/PP2A) activity. This increased the association of 4E-BP1 with eIF4E, consistent with H(2)O(2) inhibition of protein synthesis. The effects of H(2)O(2) were sufficient to override the stimulation of protein synthesis and 4E-BP1 phosphorylation induced by insulin. These results indicate that PI3K and PKB are important regulators of protein synthesis in cardiac myocytes, but other factors, including phosphatase activity, modulate the overall response.

Analysis of protein-protein interactions in the feline calicivirus replication complex

Caliciviruses are a major cause of gastroenteritis in humans and cause a wide variety of other diseases in animals. Here, the characterization of protein-protein interactions between the individual proteins of Feline calicivirus (FCV), a model system for other members of the family Caliciviridae, is reported. Using the yeast two-hybrid system combined with a number of other approaches, it is demonstrated that the p32 protein (the picornavirus 2B analogue) of FCV interacts with p39 (2C), p30 (3A) and p76 (3CD). The FCV protease/RNA polymerase (ProPol) p76 was found to form homo-oligomers, as well as to interact with VPg and ORF2, the region encoding the major capsid protein VP1. A weak interaction was also observed between p76 and the minor capsid protein encoded by ORF3 (VP2). ORF2 protein was found to interact with VPg, p76 and VP2. The potential roles of the interactions in calicivirus replication are discussed.

Delineation and modelling of a nucleolar retention signal in the coronavirus nucleocapsid protein

Unlike nuclear localization signals, there is no obvious consensus sequence for the targeting of proteins to the nucleolus. The nucleolus is a dynamic subnuclear structure which is crucial to the normal operation of the eukaryotic cell. Studying nucleolar trafficking signals is problematic as many nucleolar retention signals (NoRSs) are part of classical nuclear localization signals (NLSs). In addition, there is no known consensus signal with which to inform a study. The avian infectious bronchitis virus (IBV), coronavirus nucleocapsid (N) protein, localizes to the cytoplasm and the nucleolus. Mutagenesis was used to delineate a novel eight amino acid motif that was necessary and sufficient for nucleolar retention of N protein and colocalize with nucleolin and fibrillarin. Additionally, a classical nuclear export signal (NES) functioned to direct N protein to the cytoplasm. Comparison of the coronavirus NoRSs with known cellular and other viral NoRSs revealed that these motifs have conserved arginine residues. Molecular modelling, using the solution structure of severe acute respiratory (SARS) coronavirus N-protein, revealed that this motif is available for interaction with cellular factors which may mediate nucleolar localization. We hypothesise that the N-protein uses these signals to traffic to and from the nucleolus and the cytoplasm.

Mechanisms of G protein activation via the D-2 dopamine receptor: evidence for persistent receptor/G protein interaction after agonist stimulation

Background and purpose: The aim of this report is to study mechanisms of G protein activation by agonists. Experimental approach: The association and dissociation of guanosine 5'-O-(3-[S-35] thio) triphosphate ([S-35] GTP gamma S) binding at G proteins in membranes of CHO cells stably transfected with the human dopamine D-2short receptor was studied in the presence of a range of agonists. Key results: Binding of [S-35] GTPgS was dissociable in the absence of agonist and dissociation was accelerated both in rate and extent by dopamine, an effect which was blocked by the dopamine D-2 receptor antagonist raclopride and by suramin, which inhibits receptor/G protein interaction. A range of agonists of varying efficacy increased the rate of dissociation of [S-35] GTPgS binding, with the more efficacious agonists resulting in faster dissociation. Agonists were able to dissociate about 70% of the pre-bound [S-35] GTPgS, leaving a component which may not be accessible to the agonist-bound receptor. The dissociable component of the [S-35] GTPgS binding was reduced with longer association times and increased [S-35] GTPgS concentrations. Conclusions and implications: These data are consistent with [S-35] GTPgS binding being initially to receptor-linked G proteins and then to G proteins which have separated from the agonist bound receptor. Under the conditions used typically for [S-35] GTPgS binding assays, therefore, much of the agonist-receptor complex remains in proximity to G proteins after they have been activated by agonist.

Proteogenomics and in silico structural and functional annotation of the barley powdery mildew Blumeria graminis f. sp. hordei

Blumeria graminis is an economically important obligate plant-pathogenic fungus, whose entire genome was recently sequenced and manually annotated using ab initio in silico predictions [7]. Employing large scale proteogenomic analysis we are now able to verify independently the existence of proteins predicted by 24% of open reading frame models. We compared the haustoria and sporulating hyphae proteomes and identified 71 proteins exclusively in haustoria, the feeding and effector-delivery organs of the pathogen. These proteins are ‘significantly smaller than the rest of the protein pool and predicted to be secreted. Most do not share any similarities with Swiss–Prot or Trembl entries nor possess any identifiable Pfam domains. We used a novel automated prediction pipeline to model the 3D structures of the proteins, identify putative ligand binding sites and predict regions of intrinsic disorder. This revealed that the protein set found exclusively in haustoria is significantly less disordered than the rest of the identified Blumeria proteins or random (and representative) protein sets generated from the yeast proteome. For most of the haustorial proteins with unknown functions no good templates could be found, from which to generate high quality models. Thus, these unknown proteins present potentially new protein folds that can be specific to the interaction of the pathogen with its host.

Role of intimin-Tir interactions and the Tir-cytoskeleton coupling protein in the colonization of calves and lambs by Escherichia coli O157:H7

Intimin facilitates intestinal colonization by enterohemorrhagic Escherichia coli O157:H7; however, the importance of intimin binding to its translocated receptor (Tir) as opposed to cellular coreceptors is unknown. The intimin-Tir interaction is needed for optimal actin assembly under adherent bacteria in vitro, a process which requires the Tir-cytoskeleton coupling protein (TccP/EspF(U)) in E. coli O157:H7. Here we report that E. coli O157:H7 tir mutants are at least as attenuated as isogenic eae mutants in calves and lambs, implying that the role of intimin in the colonization of reservoir hosts can be explained largely by its binding to Tir. Mutation of tccP uncoupled actin assembly from the intimin-Tir-mediated adherence of E. coli O157:H7 in vitro but did not impair intestinal colonization in calves and lambs, implying that pedestal formation may not be necessary for persistence. However, an E. coli O157:H7 tccP mutant induced typical attaching and effacing lesions in a bovine ligated ileal loop model of infection, suggesting that TccP-independent mechanisms of actin assembly may operate in vivo.

Effect of plasma insulin and branched-chain amino acids on skeletal muscle protein synthesis in fasted lambs

The increase in fractional rate of protein synthesis (K-s) in the skeletal muscle of growing rats during the transition from fasted to fed state has been explained by the synergistic action of a rise in plasma insulin and branched-chain amino acids (BCAA). Since growing lambs Also exhibit an increase in K-s with level of feed intake, the objective of the present study was to determine if this synergistic relationship between insulin and BCAA also occurs in ruminant animals. Six 30 kg fasted (72 h) lambs (8 months of age) received each of four treatments, which were based on continuous infusion into the jugular vein for 6 h of: (1) saline (155 mmol NaCl/l); (2) a mixture of BCAA (0.778 mumol leucine, 0.640 mumol isoleucine and 0.693 mumol valine/min.kg); (3) 18.7 mumol glucose/min.kg (to induce endogenous insulin secretion): (4) co-infusion of BCAA and glucose. Within each period all animals received the same isotope of phenylalanine, (Phe) as follows: (1) L-[1-C-13]Phe; (2) L-phenyl-[ring H-2(5)]-alanine; (3) L-[N-15]Phe; (4) L-[ring 2,6-H-3]Phe. Blood was sampled serially during infusions to measure plasma concentrations of insulin, glucose and amino acids, and plasma free Phe isotopic activity; biopsies were taken 6 h after the beginning of infusions to determine K-s in in. longissimus dorsi and vastus muscle. Compared with control (saline-infused) lambs, K-s was increased by an average of 40% at the end of glucose infusion, but this effect was not statistically significant in either of the muscles sampled. BCAA infusion, alone or in combination with glucose, also had no significant effect on K-s compared with control sheep. K-s was approximately 60% greater for vastus muscle than for m. longissimus dorsi (P<0.01), regardless of treatment. It is concluded that there are signals other than insulin and BCAA that are responsible for the feed-induced increase in K-s in muscle of growing ruminant animals.

Fas ligand-induced apoptosis is regulated by nitric oxide through the inhibition of fas receptor clustering and the nitrosylation of protein kinase Cepsilon

Apoptosis induced by the death-inducing ligand FasL (CD95L) is a major mechanism of cell death. Trophoblast cells express the Fas receptor yet survive in an environment that is rich in the ligand. We report that basal nitric oxide (NO) production is responsible for the resistance of trophoblasts to FasL-induced apoptosis. In this study we demonstrate that basal NO production resulted in the inhibition of receptor clustering following ligand binding. In addition NO also protected cells through the selective nitrosylation, and inhibition, of protein kinase Cepsilon (PKCepsilon) but not PKCalpha. In the absence of NO production PKCepsilon interacted with, and phosphorylated, the anti-apoptotic protein cFLIP. The interaction is predominantly with the short form of cFLIP and its phosphorylation reduces its recruitment to the death-inducing signaling complex (DISC) that is formed following binding of a death-inducing ligand to its receptor. Inhibition of cFLIP recruitment to the DISC leads to increased activation of caspase 8 and subsequently to apoptosis. Inhibition of PKCepsilon using siRNA significantly reversed the sensitivity to apoptosis induced by inhibition of NO synthesis suggesting that NO-mediated inhibition of PKCepsilon plays an important role in the regulation of Fas-induced apoptosis.

Role of Nrf2 in the regulation of CD36 and stress protein expression in murine macrophages - Activation by oxidatively modified LDL and 4-hydroxynonenal

CD36 is an important scavenger receptor mediating uptake of oxidized low- density lipoproteins ( oxLDLs) and plays a key role in foam cell formation and the pathogenesis of atherosclerosis. We report the first evidence that the transcription factor Nrf2 is expressed in vascular smooth muscle cells, and demonstrate that oxLDLs cause nuclear accumulation of Nrf2 in murine macrophages, resulting in the activation of genes encoding CD36 and the stress proteins A170, heme oxygenase- 1 ( HO- 1), and peroxiredoxin I ( Prx I). 4- Hydroxy- 2- nonenal ( HNE), derived from lipid peroxidation, was one of the most effective activators of Nrf2. Using Nrf2- deficient macrophages, we established that Nrf2 partially regulates CD36 expression in response to oxLDLs, HNE, or the electrophilic agent diethylmaleate. In murine aortic smooth muscle cells, expressing negligible levels of CD36, both moderately and highly oxidized LDL caused only limited Nrf2 translocation and negligible increases in A170, HO- 1, and Prx I expression. However, treatment of smooth muscle cells with HNE significantly enhanced nuclear accumulation of Nrf2 and increased A170, HO- 1, and Prx I protein levels. Because PPAR-gamma can be activated by oxLDLs and controls expression of CD36 in macrophages, our results implicate Nrf2 as a second important transcription factor involved in the induction of the scavenger receptor CD36 and antioxidant stress genes in atherosclerosis.

Automatic internal calibration in liquid chromatography/Fourier transform ion cyclotron resonance mass spectrometry of protein digests

Accurately measured peptide masses can be used for large-scale protein identification from bacterial whole-cell digests as an alternative to tandem mass spectrometry (MS/MS) provided mass measurement errors of a few parts-per-million (ppm) are obtained. Fourier transform ion cyclotron resonance (FTICR) mass spectrometry (MS) routinely achieves such mass accuracy either with internal calibration or by regulating the charge in the analyzer cell. We have developed a novel and automated method for internal calibration of liquid chromatography (LC)/FTICR data from whole-cell digests using peptides in the sample identified by concurrent MS/MS together with ambient polydimethyl-cyclosiloxanes as internal calibrants in the mass spectra. The method reduced mass measurement error from 4.3 +/- 3.7 ppm to 0.3 +/- 2.3 ppm in an E. coli LC/FTICR dataset of 1000 MS and MS/MS spectra and is applicable to all analyses of complex protein digests by FTICRMS. Copyright (c) 2006 John Wiley & Sons, Ltd.

Characterization of protein-polyphenol interactions

Dietary polyphenols have received attention for their biologically significant functions as antioxidants, anticarcinogens or antimutagens, which have led to their recognition as potential nutraceuticals. Polyphenols also characteristically possess a significant binding affinity for proteins, which can lead to the formation of soluble and insoluble protein-polyphenol complexes. Questions remain concerning whether and to what extent the protein-polyphenol interaction influences functionality. For example, is the formation of protein-polyphenol complexes an obstacle to the nutritional bioavailability of either species? This article discusses the development of suitable methodologies to investigate the physicochemical basis of protein-polyphenol interactions and the influence of structure-activity relationships on binding affinities. (C) 2003 Elsevier Ltd. All rights reserved.