Functional diversification of duplicate genes through subcellular adaptation of encoded proteins.

BACKGROUND: Gene duplication is the primary source of new genes with novel or altered functions. It is known that duplicates may obtain these new functional roles by evolving divergent expression patterns and/or protein functions after the duplication event. Here, using yeast (Saccharomyces cerevisiae) as a model organism, we investigate a previously little considered mode for the functional diversification of duplicate genes: subcellular adaptation of encoded proteins. RESULTS: We show that for 24-37% of duplicate gene pairs derived from the S. cerevisiae whole-genome duplication event, the two members of the pair encode proteins that localize to distinct subcellular compartments. The propensity of yeast duplicate genes to evolve new localization patterns depends to a large extent on the biological function of their progenitor genes. Proteins involved in processes with a wider subcellular distribution (for example, catabolism) frequently evolved new protein localization patterns after duplication, whereas duplicate proteins limited to a smaller number of organelles (for example, highly expressed biosynthesis/housekeeping proteins with a slow rate of evolution) rarely relocate within the cell. Paralogous proteins evolved divergent localization patterns by partitioning of ancestral localizations ('sublocalization'), but probably more frequently by relocalization to new compartments ('neolocalization'). We show that such subcellular reprogramming may occur through selectively driven substitutions in protein targeting sequences. Notably, our data also reveal that relocated proteins functionally adapted to their new subcellular environments and evolved new functional roles through changes of their physico-chemical properties, expression levels, and interaction partners. CONCLUSION: We conclude that protein subcellular adaptation represents a common mechanism for the functional diversification of duplicate genes.

Patterns of evolution of host proteins involved in retroviral pathogenesis.

BACKGROUND: Evolutionary analysis may serve as a useful approach to identify and characterize host defense and viral proteins involved in genetic conflicts. We analyzed patterns of coding sequence evolution of genes with known (TRIM5alpha and APOBEC3G) or suspected (TRIM19/PML) roles in virus restriction, or in viral pathogenesis (PPIA, encoding Cyclophilin A), in the same set of human and non-human primate species. RESULTS AND CONCLUSION: This analysis revealed previously unidentified clusters of positively selected sites in APOBEC3G and TRIM5alpha that may delineate new virus-interaction domains. In contrast, our evolutionary analyses suggest that PPIA is not under diversifying selection in primates, consistent with the interaction of Cyclophilin A being limited to the HIV-1M/SIVcpz lineage. The strong sequence conservation of the TRIM19/PML sequences among primates suggests that this gene does not play a role in antiretroviral defense.

Suppression of short interfering RNA-mediated gene silencing by the structural proteins of hepatitis C virus.

Viruses have evolved strategies to overcome the antiviral effects of the host at different levels. Besides specific defence mechanisms, the host responds to viral infection via the interferon pathway and also by RNA interference (RNAi). However, several viruses have been identified that suppress RNAi. We addressed the question of whether hepatitis C virus (HCV) suppresses RNAi, using cell lines constitutively expressing green fluorescent protein (GFP) and inducibly expressing HCV proteins. It was found that short interfering RNA-mediated GFP gene silencing was inhibited when the entire HCV polyprotein was expressed. Further studies showed that HCV structural proteins, and in particular envelope protein 2 (E2), were responsible for this inhibition. Co-precipitation assays demonstrated that E2 bound to Argonaute-2 (Ago-2), a member of the RNA-induced silencing complex, RISC. Thus, HCV E2 that interacts with Ago-2 is able to suppress RNAi.

Modified SH2 domain to phototrap and identify phosphotyrosine proteins from subcellular sites within cells.

Spatial regulation of tyrosine phosphorylation is important for many aspects of cell biology. However, phosphotyrosine accounts for less than 1% of all phosphorylated substrates, and it is typically a very transient event in vivo. These factors complicate the identification of key tyrosine kinase substrates, especially in the context of their extraordinary spatial organization. Here, we describe an approach to identify tyrosine kinase substrates based on their subcellular distribution from within cells. This method uses an unnatural amino acid-modified Src homology 2 (SH2) domain that is expressed within cells and can covalently trap phosphotyrosine proteins on exposure to light. This SH2 domain-based photoprobe was targeted to cellular structures, such as the actin cytoskeleton, mitochondria, and cellular membranes, to capture tyrosine kinase substrates unique to each cellular region. We demonstrate that RhoA, one of the proteins associated with actin, can be phosphorylated on two tyrosine residues within the switch regions, suggesting that phosphorylation of these residues might modulate RhoA signaling to the actin cytoskeleton. We conclude that expression of SH2 domains within cellular compartments that are capable of covalent phototrapping can reveal the spatial organization of tyrosine kinase substrates that are likely to be important for the regulation of subcellular structures.

A-kinase anchoring proteins: scaffolding proteins in the heart.

The pleiotropic cyclic nucleotide cAMP is the primary second messenger responsible for autonomic regulation of cardiac inotropy, chronotropy, and lusitropy. Under conditions of prolonged catecholaminergic stimulation, cAMP also contributes to the induction of both cardiac myocyte hypertrophy and apoptosis. The formation of localized, multiprotein complexes that contain different combinations of cAMP effectors and regulatory enzymes provides the architectural infrastructure for the specialization of the cAMP signaling network. Scaffolds that bind protein kinase A are called "A-kinase anchoring proteins" (AKAPs). In this review, we discuss recent advances in our understanding of how PKA is compartmentalized within the cardiac myocyte by AKAPs and how AKAP complexes modulate cardiac function in both health and disease.

Untangling the evolution of Rab G proteins: implications of a comprehensive genomic analysis.

BACKGROUND: Membrane-bound organelles are a defining feature of eukaryotic cells, and play a central role in most of their fundamental processes. The Rab G proteins are the single largest family of proteins that participate in the traffic between organelles, with 66 Rabs encoded in the human genome. Rabs direct the organelle-specific recruitment of vesicle tethering factors, motor proteins, and regulators of membrane traffic. Each organelle or vesicle class is typically associated with one or more Rab, with the Rabs present in a particular cell reflecting that cell's complement of organelles and trafficking routes. RESULTS: Through iterative use of hidden Markov models and tree building, we classified Rabs across the eukaryotic kingdom to provide the most comprehensive view of Rab evolution obtained to date. A strikingly large repertoire of at least 20 Rabs appears to have been present in the last eukaryotic common ancestor (LECA), consistent with the 'complexity early' view of eukaryotic evolution. We were able to place these Rabs into six supergroups, giving a deep view into eukaryotic prehistory. CONCLUSIONS: Tracing the fate of the LECA Rabs revealed extensive losses with many extant eukaryotes having fewer Rabs, and none having the full complement. We found that other Rabs have expanded and diversified, including a large expansion at the dawn of metazoans, which could be followed to provide an account of the evolutionary history of all human Rabs. Some Rab changes could be correlated with differences in cellular organization, and the relative lack of variation in other families of membrane-traffic proteins suggests that it is the changes in Rabs that primarily underlies the variation in organelles between species and cell types.

Use of an in-house approach to study the three-dimensional structures of various outer membrane proteins: structure of the alcaligin outer membrane transporter FauA from Bordetella pertussis.

Bordetella pertussis is the bacterial agent of whooping cough in humans. Under iron-limiting conditions, it produces the siderophore alcaligin. Released to the extracellular environment, alcaligin chelates iron, which is then taken up as a ferric alcaligin complex via the FauA outer membrane transporter. FauA belongs to a family of TonB-dependent outer membrane transporters that function using energy derived from the proton motive force. Using an in-house protocol for membrane-protein expression, purification and crystallization, FauA was crystallized in its apo form together with three other TonB-dependent transporters from different organisms. Here, the protocol used to study FauA is described and its three-dimensional structure determined at 2.3 A resolution is discussed.

Interaction of Fas(Apo-1/CD95) with proteins implicated in the ubiquitination pathway.

Fas(Apo-1/CD95), a receptor belonging to the tumor necrosis factor receptor family, induces apoptosis when triggered by Fas ligand. Upon its activation, the cytoplasmic domain of Fas binds several proteins which transmit the death signal. We used the yeast two-hybrid screen to isolate Fas-associated proteins. Here we report that the ubiquitin-conjugating enzyme UBC9 binds to Fas at the interface between the death domain and the membrane-proximal region of Fas. This interaction is also seen in vivo. UBC9 transiently expressed in HeLa cells bound to the co-expressed cytoplasmic segment of Fas. FAF1, a Fas-associated protein that potentiates apoptosis (Chu et al. (1996) Proc. Natl. Acad. Sci. USA 92, 11894-11898), was found to contain sequences similar to ubiquitin. These results suggest that proteins related to the ubiquitination pathway may modulate the Fas signaling pathway.

Tunable reporter signal production in feedback-uncoupled arsenic bioreporters.

Escherichia coli-based bioreporters for arsenic detection are typically based on the natural feedback loop that controls ars operon transcription. Feedback loops are known to show a wide range linear response to the detriment of the overall amplification of the incoming signal. While being a favourable feature in controlling arsenic detoxification for the cell, a feedback loop is not necessarily the most optimal for obtaining highest sensitivity and response in a designed cellular reporter for arsenic detection. Here we systematically explore the effects of uncoupling the topology of arsenic sensing circuitry on the developed reporter signal as a function of arsenite concentration input. A model was developed to describe relative ArsR and GFP levels in feedback and uncoupled circuitry, which was used to explore new ArsR-based synthetic circuits. The expression of arsR was then placed under the control of a series of constitutive promoters, which differed in promoter strength, and which could be further modulated by TetR repression. Expression of the reporter gene was maintained under the ArsR-controlled Pars promoter. ArsR expression in the systems was measured by using ArsR-mCherry fusion proteins. We find that stronger constitutive ArsR production decreases arsenite-dependent EGFP output from Pars and vice versa. This leads to a tunable series of arsenite-dependent EGFP outputs in a variety of systematically characterized circuitries. The higher expression levels and sensitivities of the response curves in the uncoupled circuits may be useful for improving field-test assays using arsenic bioreporters.

Biophysical characterization of tubulin tyrosine ligase-like proteins on microtubule dynamics

Report for the scientific sojourn carried out at the Cell Biology and Biophysics Unit from the National Institutes of Health, from 2010 to 2012.

Distributed all-atom simulations of intrinsically unstructured proteins

The Computational Biophysics Group at the Universitat Pompeu Fabra (GRIB-UPF) hosts two unique computational resources dedicated to the execution of large scale molecular dynamics (MD) simulations: (a) the ACMD molecular-dynamics software, used on standard personal computers with graphical processing units (GPUs); and (b) the GPUGRID. net computing network, supported by users distributed worldwide that volunteer GPUs for biomedical research. We leveraged these resources and developed studies, protocols and open-source software to elucidate energetics and pathways of a number of biomolecular systems, with a special focus on flexible proteins with many degrees of freedom. First, we characterized ion permeation through the bactericidal model protein Gramicidin A conducting one of the largest studies to date with the steered MD biasing methodology. Next, we addressed an open problem in structural biology, the determination of drug-protein association kinetics; we reconstructed the binding free energy, association, and dissaciociation rates of a drug like model system through a spatial decomposition and a Makov-chain analysis. The work was published in the Proceedings of the National Academy of Sciences and become one of the few landmark papers elucidating a ligand-binding pathway. Furthermore, we investigated the unstructured Kinase Inducible Domain (KID), a 28-peptide central to signalling and transcriptional response; the kinetics of this challenging system was modelled with a Markovian approach in collaboration with Frank Noe’s group at the Freie University of Berlin. The impact of the funding includes three peer-reviewed publication on high-impact journals; three more papers under review; four MD analysis components, released as open-source software; MD protocols; didactic material, and code for the hosting group.

Docking to heme proteins.

In silico screening has become a valuable tool in drug design, but some drug targets represent real challenges for docking algorithms. This is especially true for metalloproteins, whose interactions with ligands are difficult to parametrize. Our docking algorithm, EADock, is based on the CHARMM force field, which assures a physically sound scoring function and a good transferability to a wide range of systems, but also exhibits difficulties in case of some metalloproteins. Here, we consider the therapeutically important case of heme proteins featuring an iron core at the active site. Using a standard docking protocol, where the iron-ligand interaction is underestimated, we obtained a success rate of 28% for a test set of 50 heme-containing complexes with iron-ligand contact. By introducing Morse-like metal binding potentials (MMBP), which are fitted to reproduce density functional theory calculations, we are able to increase the success rate to 62%. The remaining failures are mainly due to specific ligand-water interactions in the X-ray structures. Testing of the MMBP on a second data set of non iron binders (14 cases) demonstrates that they do not introduce a spurious bias towards metal binding, which suggests that they may reliably be used also for cross-docking studies.