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.

The long form of FLIP is an activator of caspase-8 at the Fas death-inducing signaling complex.

Death receptors, such as Fas and tumor necrosis factor-related apoptosis-inducing ligand receptors, recruit Fas-associated death domain and pro-caspase-8 homodimers, which are then autoproteolytically activated. Active caspase-8 is released into the cytoplasm, where it cleaves various proteins including pro-caspase-3, resulting in apoptosis. The cellular Fas-associated death domain-like interleukin-1-beta-converting enzyme-inhibitory protein long form (FLIP(L)), a structural homologue of caspase-8 lacking caspase activity because of several mutations in the active site, is a potent inhibitor of death receptor-induced apoptosis. FLIP(L) is proposed to block caspase-8 activity by forming a proteolytically inactive heterodimer with caspase-8. In contrast, we propose that FLIP(L)-bound caspase-8 is an active protease. Upon heterocomplex formation, a limited caspase-8 autoprocessing occurs resulting in the generation of the p43/41 and the p12 subunits. This partially processed form but also the non-cleaved FLIP(L)-caspase-8 heterocomplex are proteolytically active because they both bind synthetic substrates efficiently. Moreover, FLIP(L) expression favors receptor-interacting kinase (RIP) processing within the Fas-signaling complex. We propose that FLIP(L) inhibits caspase-8 release-dependent pro-apoptotic signals, whereas the single, membrane-restricted active site of the FLIP(L)-caspase-8 heterocomplex is proteolytically active and acts on local substrates such as RIP.

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.

Sequencing and characterizing the genome of Estrella lausannensis as an undergraduate project: training students and biological insights.

With the widespread availability of high-throughput sequencing technologies, sequencing projects have become pervasive in the molecular life sciences. The huge bulk of data generated daily must be analyzed further by biologists with skills in bioinformatics and by "embedded bioinformaticians," i.e., bioinformaticians integrated in wet lab research groups. Thus, students interested in molecular life sciences must be trained in the main steps of genomics: sequencing, assembly, annotation and analysis. To reach that goal, a practical course has been set up for master students at the University of Lausanne: the "Sequence a genome" class. At the beginning of the academic year, a few bacterial species whose genome is unknown are provided to the students, who sequence and assemble the genome(s) and perform manual annotation. Here, we report the progress of the first class from September 2010 to June 2011 and the results obtained by seven master students who specifically assembled and annotated the genome of Estrella lausannensis, an obligate intracellular bacterium related to Chlamydia. The draft genome of Estrella is composed of 29 scaffolds encompassing 2,819,825 bp that encode for 2233 putative proteins. Estrella also possesses a 9136 bp plasmid that encodes for 14 genes, among which we found an integrase and a toxin/antitoxin module. Like all other members of the Chlamydiales order, Estrella possesses a highly conserved type III secretion system, considered as a key virulence factor. The annotation of the Estrella genome also allowed the characterization of the metabolic abilities of this strictly intracellular bacterium. Altogether, the students provided the scientific community with the Estrella genome sequence and a preliminary understanding of the biology of this recently-discovered bacterial genus, while learning to use cutting-edge technologies for sequencing and to perform bioinformatics analyses.

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.

Development of a well-defined medium for the in vitro maturation of immature bovine cumulus-oocyte complexes.

PURPOSE: Our purpose was to develop a well-defined medium for the in vitro maturation (IVM) of immature bovine cumulus-oocyte complexes (COC). METHODS: The COC were cultured in the presence of three protein supplementations: fetal bovine serum (FBS), bovine serum albumin, and Synthetic Serum Substitute. The embryos obtained after in vitro fertilization of IVM oocytes were cocultured with Vero cells and their development to the morula and blastocyst stages was studied. RESULTS: When FBS was absent from the IVM medium, a significantly lower fertilization rate was observed, followed by a decrease in the percentage of embryos reaching the blastocyst stage. When FBS was replaced by a defined protein supplementation, the best results were obtained with Synthetic Serum Substitute. CONCLUSIONS: Adequate protein supplementation of the IVM medium optimizes the fertilization rate and the development of bovine IVM oocytes. The implication of these results in the human field is discussed.

A TARBP2-dependent miRNA expression profile underlies cancer stem cell properties and provides candidate therapeutic reagents in Ewing sarcoma.

We have recently demonstrated that human pediatric mesenchymal stem cells can be reprogrammed toward a Ewing sarcoma family tumor (ESFT) cancer stem cell (CSC) phenotype by mechanisms that implicate microRNAs (miRNAs). Here, we show that the miRNA profile of ESFT CSCs is shared by embryonic stem cells and CSCs from divergent tumor types. We also provide evidence that the miRNA profile of ESFT CSCs is the result of reversible disruption of TARBP2-dependent miRNA maturation. Restoration of TARBP2 activity and systemic delivery of synthetic forms of either of two of its targets, miRNA-143 or miRNA-145, inhibited ESFT CSC clonogenicity and tumor growth in vivo. Our observations suggest that CSC self-renewal and tumor maintenance may depend on deregulation of TARBP2-dependent miRNA expression.

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.

Splicing and the evolution of proteins in mammals.

It is often supposed that a protein's rate of evolution and its amino acid content are determined by the function and anatomy of the protein. Here we examine an alternative possibility, namely that the requirement to specify in the unprocessed RNA, in the vicinity of intron-exon boundaries, information necessary for removal of introns (e.g., exonic splice enhancers) affects both amino acid usage and rates of protein evolution. We find that the majority of amino acids show skewed usage near intron-exon boundaries, and that differences in the trends for the 2-fold and 4-fold blocks of both arginine and leucine show this to be owing to effects mediated at the nucleotide level. More specifically, there is a robust relationship between the extent to which an amino acid is preferred/avoided near boundaries and its enrichment/paucity in splice enhancers. As might then be expected, the rate of evolution is lowest near intron-exon boundaries, at least in part owing to splice enhancers, such that domains flanking intron-exon junctions evolve on average at under half the rate of exon centres from the same gene. In contrast, the rate of evolution of intronless retrogenes is highest near the domains where intron-exon junctions previously resided. The proportion of sequence near intron-exon boundaries is one of the stronger predictors of a protein's rate of evolution in mammals yet described. We conclude that after intron insertion selection favours modification of amino acid content near intron-exon junctions, so as to enable efficient intron removal, these changes then being subject to strong purifying selection even if nonoptimal for protein function. Thus there exists a strong force operating on protein evolution in mammals that is not explained directly in terms of the biology of the protein.

KnotProt: a database of proteins with knots and slipknots.

The protein topology database KnotProt, http://knotprot.cent.uw.edu.pl/, collects information about protein structures with open polypeptide chains forming knots or slipknots. The knotting complexity of the cataloged proteins is presented in the form of a matrix diagram that shows users the knot type of the entire polypeptide chain and of each of its subchains. The pattern visible in the matrix gives the knotting fingerprint of a given protein and permits users to determine, for example, the minimal length of the knotted regions (knot's core size) or the depth of a knot, i.e. how many amino acids can be removed from either end of the cataloged protein structure before converting it from a knot to a different type of knot. In addition, the database presents extensive information about the biological functions, families and fold types of proteins with non-trivial knotting. As an additional feature, the KnotProt database enables users to submit protein or polymer chains and generate their knotting fingerprints.

Human synthetic lethal inference as potential anti-cancer target gene detection

Background: Two genes are called synthetic lethal (SL) if mutation of either alone is not lethal, but mutation of both leads to death or a significant decrease in organism's fitness. The detection of SL gene pairs constitutes a promising alternative for anti-cancer therapy. As cancer cells exhibit a large number of mutations, the identification of these mutated genes' SL partners may provide specific anti-cancer drug candidates, with minor perturbations to the healthy cells. Since existent SL data is mainly restricted to yeast screenings, the road towards human SL candidates is limited to inference methods. Results: In the present work, we use phylogenetic analysis and database manipulation (BioGRID for interactions, Ensembl and NCBI for homology, Gene Ontology for GO attributes) in order to reconstruct the phylogenetically-inferred SL gene network for human. In addition, available data on cancer mutated genes (COSMIC and Cancer Gene Census databases) as well as on existent approved drugs (DrugBank database) supports our selection of cancer-therapy candidates.Conclusions: Our work provides a complementary alternative to the current methods for drug discovering and gene target identification in anti-cancer research. Novel SL screening analysis and the use of highly curated databases would contribute to improve the results of this methodology.