A lock-and-key model for protein-protein interactions

MOTIVATION: Protein-protein interaction networks are one of the major post-genomic data sources available to molecular biologists. They provide a comprehensive view of the global interaction structure of an organism's proteome, as well as detailed information on specific interactions. Here we suggest a physical model of protein interactions that can be used to extract additional information at an intermediate level: It enables us to identify proteins which share biological interaction motifs, and also to identify potentially missing or spurious interactions. RESULTS: Our new graph model explains observed interactions between proteins by an underlying interaction of complementary binding domains (lock-and-key model). This leads to a novel graph-theoretical algorithm to identify bipartite subgraphs within protein-protein interaction networks where the underlying data are taken from yeast two-hybrid experimental results. By testing on synthetic data, we demonstrate that under certain modelling assumptions, the algorithm will return correct domain information about each protein in the network. Tests on data from various model organisms show that the local and global patterns predicted by the model are indeed found in experimental data. Using functional and protein structure annotations, we show that bipartite subnetworks can be identified that correspond to biologically relevant interaction motifs. Some of these are novel and we discuss an example involving SH3 domains from the Saccharomyces cerevisiae interactome. AVAILABILITY: The algorithm (in Matlab format) is available (see http://www.maths.strath.ac.uk/~aas96106/lock_key.html).     

A second gene for peroxisomal HMG-CoA reductase? A genomic reassessment

HMG-CoA reductase (HMGCR) catalyzes the conversion of HMG-CoA to mevalonate, the rate-limiting step of eukaryotic isoprenoid biosynthesis, and is the main target of cholesterol-lowering drugs. The classical form of the enzyme is a transmembrane-protein anchored to the endoplasmic reticulum. However, during the last years several lines of evidence pointed to the existence of a second isoform of HMGCR localized in peroxisomes, where mevalonate is converted further to farnesyl diphosphate. This finding is relevant for our understanding of the complex regulation and compartmentalization of the cholesterogenic pathway. Here we review experimental evidence suggesting that the peroxisomal activity might be due to a second HMGCR gene in mammals. We then present a comprehensive analysis of completely sequenced eukaryotic genomes, as well as the human and mouse genome drafts. Our results provide evidence for a large number of independent duplications of HMGCR in all eukaryotic kingdoms, but not for a second gene in mammals. We conclude that the peroxisomal HMGCR activity in mammals is due to alternative targeting of the ER enzyme to peroxisomes by an as yet uncharacterized mechanism.     

Surface display of antibodies

To screen antibody libraries that contain many millions of different clones, a selection system is required with an efficiency comparable to that of the immune system. This can be achieved by displaying antibodies on the surface of microorganisms containing the antibody's gene, analogous to the expression of the IgM antigen receptor on the surface of unactivated B-lymphocytes. Specific clones can then be selected using immobilized antigens. The minor coat protein of filamentous phages, pIII, which initiates the infection of E.coli by binding to their F-pili, and the major coat protein, pVIII, have been used as carriers for displaying antibodies on the phage surface. Recombinant antibodies have also been targeted to the cell surface of bacteria by fusing them with outer membrane components derived from lipoproteins, OmpA and an IgA protease. However, only the pIII system has been routinely used for screening antibody libraries. Here we describe the various antibody surface display systems and the screening of antibody libraries generated from the gene repertoire of lymphocytes and by gene synthesis. Finally, we have made a short comparison of the bacterial production of Fabs versus single chain antibodies (scFv).     

Measurement of the adhesive strength of a gingival cell line to agar

A new technique is described for measuring the adhesive strength of a gingival cell line to an agar substratum by the modification of the original "blister" test for adhesives. A cell monolayer was developed on a Petri dish with a hole in the center of the growing surface, overlayed with agar, and the system pressurized to debond the cells from the agar surface. Pressure changes were measured by a capacitance pressure transducer the output of which was measured by a strip-chart recorder. The modulus (E) of the agar overlay was determined and used in the calculation of the adhesive-bond strength (gammaalpha). The gammaalpha yield for the gingival cell line (cell-agar debond) was 48.8 ergs per cm2, and for the control (no cells) (agar-polystyrene debond) was 30.0 ergs per cm2.     

Structure-based phylogenetic analysis of short-chain alcohol dehydrogenases and reclassification of the 17beta-hydroxysteroid dehydrogenase family.

Short-chain alcohol dehydrogenases (SCAD) constitute a large and diverse family of ancient origin. Several of its members play an important role in human physiology and disease, especially in the metabolism of steroid substrates (e.g., prostaglandins, estrogens, androgens, and corticosteroids). Their involvement in common human disorders such as endocrine-related cancer, osteoporosis, and Alzheimer disease makes them an important candidate for drug targets. Recent phylogenetic analysis of SCAD is incomplete and does not allow any conclusions on very ancient divergences or on a functional characterization of novel proteins within this complex family. We have developed a 3D structure-based approach to establish the deep-branching pattern within the SCAD family. In this approach, pairwise superpositions of X-ray structures were used to calculate similarity scores as an input for a tree-building algorithm. The resulting phylogeny was validated by comparison with the results of sequence-based algorithms and biochemical data. It was possible to use the 3D data as a template for the reliable determination of the phylogenetic position of novel proteins as a first step toward functional predictions. We were able to discern new patterns in the phylogenetic relationships of the SCAD family, including a basal dichotomy of the 17beta-hydroxysteroid dehydrogenases (17beta-HSDs). These data provide an important contribution toward the development of type-specific inhibitors for 17beta-HSDs for the treatment and prevention of disease. Our structure-based phylogenetic approach can also be applied to increase the reliability of evolutionary reconstructions in other large protein families.     

Inhibition of 17beta-hydroxysteroid dehydrogenases by phytoestrogens: comparison with other steroid metabolizing enzymes

Effects of phytoestrogens on human health have been reported for decades. These include not only beneficial action in cancer prevention but also endocrine disruption in males. Since then many molecular mechanisms underlying these effects have been identified. Targets of phytoestrogens comprise steroid receptors, steroid metabolising enzymes, elements of signal transduction and apoptosis pathways, and even the DNA processing machinery. Understanding the specific versus pleiotropic effects of selected phytoestrogens will be crucial for their biomedical application. This review will concentrate on the influence of phytoestrogens on 17beta-hydroxysteroid dehydrogenases from a comparative perspective with other steroid metabolizing enzymes.     

Mass appeal: metabolite identification in mass spectrometry-focused untargeted metabolomics

¹H nuclear magnetic resonance (¹H NMR)-based metabolomics was utilized to elucidate the earthworm sub-lethal toxicity after exposure to the persistent environmental contaminant phenanthrene. Earthworms were exposed to 0.05, 0.2 and 0.4 mg/cm² of phenanthrene [which correspond to 1/32nd to 1/4th of the 48-h LC₅₀ (concentration that causes 50 % mortality), respectively] via contact tests over 1, 2 and 3 days of dermal contact. ¹H NMR-based metabolomic analysis of the polar and non-polar fractions of the earthworm tissue extracts revealed heightened Eisenia fetida toxic responses with both longer exposure times and higher phenanthrene concentrations. Principal component analysis (PCA) of the polar fraction showed significant separation between control and exposed earthworms along PC1 for all phenanthrene concentrations on each day. The PCA of the non-polar fraction showed significant separation between the controls and exposed earthworms for only the first day of exposure. These results suggested that alanine, glutamate, maltose, and fatty acids were potential indicators of phenanthrene exposure. Interruption in energy production due to a deactivation of the succinate dehydrogenase enzyme in the Krebs cycle was also postulated in exposed earthworms. Cross-validated partial least squares-regression models showed that the polar metabolic profile of E. fetida was weakly but significantly correlated to phenanthrene exposure concentrations after day 1 and day 2 of exposure. Overall, this study indicates that with longer exposures, contact time becomes more important than concentration in discriminating between control and exposed earthworms. This study also shows that NMR-based metabolomics has promise as a powerful ecotoxicological tool for elucidating the mode of toxicity of contaminants.     

Handling Uncertainty in Dynamic Models: The Pentose Phosphate Pathway in Trypanosoma brucei

Dynamic models of metabolism can be useful in identifying potential drug targets, especially in unicellular organisms. A model of glycolysis in the causative agent of human African trypanosomiasis, Trypanosoma brucei, has already shown the utility of this approach. Here we add the pentose phosphate pathway (PPP) of T. brucei to the glycolytic model. The PPP is localized to both the cytosol and the glycosome and adding it to the glycolytic model without further adjustments leads to a draining of the essential bound-phosphate moiety within the glycosome. This phosphate “leak” must be resolved for the model to be a reasonable representation of parasite physiology. Two main types of theoretical solution to the problem could be identified: (i) including additional enzymatic reactions in the glycosome, or (ii) adding a mechanism to transfer bound phosphates between cytosol and glycosome. One example of the first type of solution would be the presence of a glycosomal ribokinase to regenerate ATP from ribose 5-phosphate and ADP. Experimental characterization of ribokinase in T. brucei showed that very low enzyme levels are sufficient for parasite survival, indicating that other mechanisms are required in controlling the phosphate leak. Examples of the second type would involve the presence of an ATP:ADP exchanger or recently described permeability pores in the glycosomal membrane, although the current absence of identified genes encoding such molecules impedes experimental testing by genetic manipulation. Confronted with this uncertainty, we present a modeling strategy that identifies robust predictions in the context of incomplete system characterization. We illustrate this strategy by exploring the mechanism underlying the essential function of one of the PPP enzymes, and validate it by confirming the model predictions experimentally.