A functional proteomic method for the enrichment of peripheral membrane proteins reveals the collagen binding protein Hsp47 is exposed on the surface of activated human platelets

Platelets are small blood cells vital for hemostasis. Following vascular damage, platelets adhere to collagens and activate, forming a thrombus that plugs the wound and prevents blood loss. Stimulation of the platelet collagen receptor glycoprotein VI (GPVI) allows recruitment of proteins to receptor-proximal signaling complexes on the inner-leaflet of the plasma membrane. These proteins are often present at low concentrations; therefore, signaling-complex characterization using mass spectrometry is limited due to high sample complexity. We describe a method that facilitates detection of signaling proteins concentrated on membranes. Peripheral membrane proteins (reversibly associated with membranes) were eluted from human platelets with alkaline sodium carbonate. Liquid-phase isoelectric focusing and gel electrophoresis were used to identify proteins that changed in levels on membranes from GPVI-stimulated platelets. Immunoblot analysis verified protein recruitment to platelet membranes and subsequent protein phosphorylation was preserved. Hsp47, a collagen binding protein, was among the proteins identified and found to be exposed on the surface of GPVI-activated platelets. Inhibition of Hsp47 abolished platelet aggregation in response to collagen, while only partially reducing aggregation in response to other platelet agonists. We propose that Hsp47 may therefore play a role in hemostasis and thrombosis.

Transcriptomic analysis of rhizobium leguminosarum biovar viciae in symbiosis with host plants pisum sativum and Vicia cracca

Rhizobium leguminosarum bv. viciae forms nitrogen-fixing nodules on several legumes, including pea (Pisum sativum) and vetch (Vicia cracca), and has been widely used as a model to study nodule biochemistry. To understand the complex biochemical and developmental changes undergone by R. leguminosarum bv. viciae during bacteroid development, microarray experiments were first performed with cultured bacteria grown on a variety of carbon substrates (glucose, pyruvate, succinate, inositol, acetate, and acetoacetate) and then compared to bacteroids. Bacteroid metabolism is essentially that of dicarboxylate-grown cells (i.e., induction of dicarboxylate transport, gluconeogenesis and alanine synthesis, and repression of sugar utilization). The decarboxylating arm of the tricarboxylic acid cycle is highly induced, as is gamma-aminobutyrate metabolism, particularly in bacteroids from early (7-day) nodules. To investigate bacteroid development, gene expression in bacteroids was analyzed at 7, 15, and 21 days postinoculation of peas. This revealed that bacterial rRNA isolated from pea, but not vetch, is extensively processed in mature bacteroids. In early development (7 days), there were large changes in the expression of regulators, exported and cell surface molecules, multidrug exporters, and heat and cold shock proteins. fix genes were induced early but continued to increase in mature bacteroids, while nif genes were induced strongly in older bacteroids. Mutation of 37 genes that were strongly upregulated in mature bacteroids revealed that none were essential for nitrogen fixation. However, screening of 3,072 mini-Tn5 mutants on peas revealed previously uncharacterized genes essential for nitrogen fixation. These encoded a potential magnesium transporter, an AAA domain protein, and proteins involved in cytochrome synthesis.

Functional genomics in zebrafish permits rapid characterization of novel platelet membrane proteins

In this study, we demonstrate the suitability of the vertebrate Danio rerio (zebrafish) for functional screening of novel platelet genes in vivo by reverse genetics. Comparative transcript analysis of platelets and their precursor cell, the megakaryocyte, together with nucleated blood cell elements, endothelial cells, and erythroblasts, identified novel platelet membrane proteins with hitherto unknown roles in thrombus formation. We determined the phenotype induced by antisense morpholino oligonucleotide (MO)–based knockdown of 5 of these genes in a laser-induced arterial thrombosis model. To validate the model, the genes for platelet glycoprotein (GP) IIb and the coagulation protein factor VIII were targeted. MO-injected fish showed normal thrombus initiation but severely impaired thrombus growth, consistent with the mouse knockout phenotypes, and concomitant knockdown of both resulted in spontaneous bleeding. Knockdown of 4 of the 5 novel platelet proteins altered arterial thrombosis, as demonstrated by modified kinetics of thrombus initiation and/or development. We identified a putative role for BAMBI and LRRC32 in promotion and DCBLD2 and ESAM in inhibition of thrombus formation. We conclude that phenotypic analysis of MO-injected zebrafish is a fast and powerful method for initial screening of novel platelet proteins for function in thrombosis.

Rapid acidification and alkylation: Redox analysis of the MHC class I pathway

The technique of rapid acidification and alkylation can be used to characterise the redox status of oxidoreductases, and to determine numbers of free cysteine residues within substrate proteins. We have previously used this method to analyse interacting components of the MHC class I pathway, namely ERp57 and tapasin. Here, we have applied rapid acidification alkylation as a novel approach to analysing the redox status of MHC class I molecules. This analysis of the redox status of the MHC class I molecules HLA-A2 and HLA-B27, which is strongly associated with a group of inflammatory arthritic disorders referred to as Spondyloarthropathies, revealed structural and conformational information. We propose that this assay provides a useful tool in the study of in vivo MHC class I structure. (c) 2008 Elsevier B.V. All rights reserved.

Proteomic analysis of resting and thrombin-stimulated platelets reveals the translocation and functional relevance of HIP-55 in platelets

The platelet surface is a dynamic interface that changes rapidly in response to stimuli to coordinate the formation of thrombi at sites of vascular injury. Tight control is essential as loss of organisation may result in the inappropriate formation of thrombi (thrombosis) or excessive bleeding. In this paper we describe the comparative analysis of resting and thrombin-stimulated platelet membrane proteomes and associated proteins to identify proteins important to platelet function. Surface proteins were labelled using a biotin tag and isolated by NeurtrAvidin affinity chromatography. Liquid phase IEF and SDS-PAGE were used to separate proteins, and bands of increased intensity in the stimulated platelet fractions were digested and identified by FT-ICR mass spectrometry. Novel proteins were identified along with proteins known to be translocated to the platelet surface. Furthermore, many platelet proteins revealed changes in location associated with function, including G6B and Hip-55. HIP-55 is an SH3-binding protein important in T-cell receptor signalling. Further analysis of HIP-55 revealed that this adaptor protein becomes increasingly associated with both Syk and integrin beta 3 upon platelet activation. Analysis of HIP-55 deficient platelets revealed reduced fibrinogen binding upon thrombin stimulation, suggesting HIP-55 to be an important regulator of platelet function.

Intrinsic disorder prediction from the analysis of multiple protein fold recognition models

Motivation: Intrinsic protein disorder is functionally implicated in numerous biological roles and is, therefore, ubiquitous in proteins from all three kingdoms of life. Determining the disordered regions in proteins presents a challenge for experimental methods and so recently there has been much focus on the development of improved predictive methods. In this article, a novel technique for disorder prediction, called DISOclust, is described, which is based on the analysis of multiple protein fold recognition models. The DISOclust method is rigorously benchmarked against the top.ve methods from the CASP7 experiment. In addition, the optimal consensus of the tested methods is determined and the added value from each method is quantified. Results: The DISOclust method is shown to add the most value to a simple consensus of methods, even in the absence of target sequence homology to known structures. A simple consensus of methods that includes DISOclust can significantly outperform all of the previous individual methods tested.

Proteornics analysis unravels the functional repertoire of coronavirus nonstructural protein 3

Severe acute respiratory syndrome (SARS) coronavirus infection and growth are dependent on initiating signaling and enzyme actions upon viral entry into the host cell. Proteins packaged during virus assembly may subsequently form the first line of attack and host manipulation upon infection. A complete characterization of virion components is therefore important to understanding the dynamics of early stages of infection. Mass spectrometry and kinase profiling techniques identified nearly 200 incorporated host and viral proteins. We used published interaction data to identify hubs of connectivity with potential significance for virion formation. Surprisingly, the hub with the most potential connections was not the viral M protein but the nonstructurall protein 3 (nsp3), which is one of the novel virion components identified by mass spectrometry. Based on new experimental data and a bioinformatics analysis across the Coronaviridae, we propose a higher-resolution functional domain architecture for nsp3 that determines the interaction capacity of this protein. Using recombinant protein domains expressed in Escherichia coli, we identified two additional RNA-binding domains of nsp3. One of these domains is located within the previously described SARS-unique domain, and there is a nucleic acid chaperone-like domain located immediately downstream of the papain-like proteinase domain. We also identified a novel cysteine-coordinated metal ion-binding domain. Analyses of interdomain interactions and provisional functional annotation of the remaining, so-far-uncharacterized domains are presented. Overall, the ensemble of data surveyed here paint a more complete picture of nsp3 as a conserved component of the viral protein processing machinery, which is intimately associated with viral RNA in its role as a virion component.

Hepatitis C virus NS5A protein binds the SH3 domain of the Fyn tyrosine kinase with high affinity: mutagenic analysis of residues within the SH3 domain that contribute to the interaction

Background: The hepatitis C virus (HCV) non-structural 5A protein (NS5A) contains a highly conserved C-terminal polyproline motif with the consensus sequence Pro-X-X- Pro-X-Arg that is able to interact with the Src-homology 3 (SH3) domains of a variety of cellular proteins. Results: To understand this interaction in more detail we have expressed two N-terminally truncated forms of NS5A in E. coli and examined their interactions with the SH3 domain of the Src-family tyrosine kinase, Fyn. Surface plasmon resonance analysis revealed that NS5A binds to the Fyn SH3 domain with what can be considered a high affinity SH3 domain-ligand interaction (629 nM), and this binding did not require the presence of domain I of NS5A (amino acid residues 32-250). Mutagenic analysis of the Fyn SH3 domain demonstrated the requirement for an acidic cluster at the C-terminus of the RT-Src loop of the SH3 domain, as well as several highly conserved residues previously shown to participate in SH3 domain peptide binding. Conclusion: We conclude that the NS5A: Fyn SH3 domain interaction occurs via a canonical SH3 domain binding site and the high affinity of the interaction suggests that NS5A would be able to compete with cognate Fyn ligands within the infected cell.

Assembly rules for protein networks derived from phylogenetic-statistical analysis of whole genomes

Background: We report an analysis of a protein network of functionally linked proteins, identified from a phylogenetic statistical analysis of complete eukaryotic genomes. Phylogenetic methods identify pairs of proteins that co-evolve on a phylogenetic tree, and have been shown to have a high probability of correctly identifying known functional links. Results: The eukaryotic correlated evolution network we derive displays the familiar power law scaling of connectivity. We introduce the use of explicit phylogenetic methods to reconstruct the ancestral presence or absence of proteins at the interior nodes of a phylogeny of eukaryote species. We find that the connectivity distribution of proteins at the point they arise on the tree and join the network follows a power law, as does the connectivity distribution of proteins at the time they are lost from the network. Proteins resident in the network acquire connections over time, but we find no evidence that 'preferential attachment' - the phenomenon of newly acquired connections in the network being more likely to be made to proteins with large numbers of connections - influences the network structure. We derive a 'variable rate of attachment' model in which proteins vary in their propensity to form network interactions independently of how many connections they have or of the total number of connections in the network, and show how this model can produce apparent power-law scaling without preferential attachment. Conclusion: A few simple rules can explain the topological structure and evolutionary changes to protein-interaction networks: most change is concentrated in satellite proteins of low connectivity and small phenotypic effect, and proteins differ in their propensity to form attachments. Given these rules of assembly, power law scaled networks naturally emerge from simple principles of selection, yielding protein interaction networks that retain a high-degree of robustness on short time scales and evolvability on longer evolutionary time scales.

Analysing the ability to retain sidechain hydrogen-bonds in mutant proteins

Motivation: Hydrogen bonds are one of the most important inter-atomic interactions in biology. Previous experimental, theoretical and bioinformatics analyses have shown that the hydrogen bonding potential of amino acids is generally satisfied and that buried unsatisfied hydrogen-bond-capable residues are destabilizing. When studying mutant proteins, or introducing mutations to residues involved in hydrogen bonding, one needs to know whether a hydrogen bond can be maintained. Our aim, therefore, was to develop a rapid method to evaluate whether a sidechain can form a hydrogen-bond. Results: A novel knowledge-based approach was developed in which the conformations accessible to the residues involved are taken into account. Residues involved in hydrogen bonds in a set of high resolution crystal structures were analyzed and this analysis is then applied to a given protein. The program was applied to assess mutations in the tumour-suppressor protein, p53. This raised the number of distinct mutations identified as disrupting sidechain-sidechain hydrogen bonding from 181 in our previous analysis to 202 in this analysis.

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favourability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB) versus Thr(HB)-Asn(nHB). Significantly (p <= 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi(1) rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi(1) conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi(1) and chi(2) conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favourability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http:// www.rubic.rdg.ac.uk/betapairprefsparallel/). (c) 2005 Elsevier Ltd. All rights reserved.

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.

Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.

Thermodynamic analysis of ferrous ion binding to Escherichia coli ferritin EcFtnA

Iron oxidation in the bacterial ferritin EcFtnA from Escherichia coli shows marked differences from its homologue human H-chain ferritin (HuHF). While the amino acid residues that constitute the dinuclear center in these proteins are highly conserved, EcFtnA has a third iron-binding site (C site) in close proximity to the dinuclear center that is seemingly responsible for these differences. Here, we describe the first thermodynamic study of Fe2+ binding to EcFtnA and its variants to determine the location of the primary ferrous ion-binding sites on the protein and to better understand the role of the third C site in iron binding. Isothermal titration calorimetric analyses of the wild-type protein reveal the presence of two main classes of binding sites in the pH range of 6.5-7.5, ascribed to Fe2+ binding, first at the A and then the B sites. Site-directed mutagenesis of ligands in the A, B, or C sites affects the apparent Fe2+-binding stoichiometries at the unaltered sites. The data imply some degree of inter- and intrasubunit negative cooperative interaction between sites. Unlike HuHF where only the A site initially binds Fe2+, both A and B sites in EcFtnA bind Fe2+, implying a role for the C site in influencing the binding of Fe2+ at the B site of the di-iron center of EcFtnA. The ITC equations describing a binding model for three classes of independent binding sites are reported here for the first time.