Genetic Variations and Diseases in UniProtKB/Swiss‐Prot: The Ins and Outs of Expert Manual Curation

During the last few years, next‐generation sequencing (NGS) technologies have accelerated the detection of genetic variants resulting in the rapid discovery of new disease‐associated genes. However, the wealth of variation data made available by NGS alone is not sufficient to understand the mechanisms underlying disease pathogenesis and manifestation. Multidisciplinary approaches combining sequence and clinical data with prior biological knowledge are needed to unravel the role of genetic variants in human health and disease. In this context, it is crucial that these data are linked, organized, and made readily available through reliable online resources. The Swiss‐Prot section of the Universal Protein Knowledgebase (UniProtKB/Swiss‐Prot) provides the scientific community with a collection of information on protein functions, interactions, biological pathways, as well as human genetic diseases and variants, all manually reviewed by experts. In this article, we present an overview of the information content of UniProtKB/Swiss‐Prot to show how this knowledgebase can support researchers in the elucidation of the mechanisms leading from a molecular defect to a disease phenotype.     

Molecular docking of competitive phosphodiesterase inhibitors.

Mammalian phosphodiesterases types 3 and 4 (PDE3 and PDE4) hydrolyze cAMP and are essential for the regulation of this intracellular second messenger. These enzymes share structural and biochemical similarities, but each can be distinguished by its sensitivity to isoenzyme-specific, substrate-competitive inhibitors. We present a model configuration for the PDE4 substrate (cAMP) and a PDE4-specific inhibitor (rolipram) within the active site of the enzyme. The docked models were also used to examine the structural consequences of mutations that confer resistance to rolipram and other PDE4-specific inhibitors. The proposed rolipram-binding configuration is consistent with the substrate-competitive nature of inhibition and also provides a structural basis for the observed specificity of binding to the R- versus S-enantiomer. For mutations that render the enzyme rolipram-insensitive, there was generally an inverse relationship between the magnitude of the drug resistance and the distance of the altered residue from the predicted binding site. We observed a direct correlation between the net loss of protein residue interactions (van der Waals contacts and hydrogen bond interactions) and the degree of rolipram resistance. The positions of several drug sensitivity-determinant residues define a surface leading to the substrate- and drug-binding sites, suggesting a possible approach channel leading to the enzyme active site. The binding of other PDE4 inhibitors (high- and low-affinity) was also modeled and used to predict the involvement of residues that were not previously implicated in pharmacological interactions.     

The direct effects of tacrolimus and cyclosporin A on isolated human islets: A functional, survival and gene expression study

The use of immunosuppressive drugs in transplanted patients is associated with the development of diabetes, possibly due to β-cell toxicity. To better understand the mechanisms leading to post-transplant diabetes, we investigated the actions of prolonged exposure of isolated human islets to therapeutical levels of tacrolimus (Tac) or cyclosporin A (CsA). Islets were isolated from the pancreas of multiorgan donors by enzymatic digestion and density gradient centrifugation. Functional, survival and molecular studies were then performed after 4 days of incubation with therapeutical concentrations of Tac or CsA. Glucose-induced insulin secretion was significantly decreased in Tac, but not in CsA exposed islets, which was associated with a reduction of the amount of insulin granules as shown by electron microscopy. The percentage of apoptotic β-cells was higher in Tac than CsA exposed islets. Microarray experiments followed by Gene Set Enrichment Analysis revealed that gene expression was more markedly affected upon Tac treatment. In conclusion, Tac and CsA affect features of beta-cell differently, with several changes occurring at the molecular level.     

Localizing proteins in the cell from their phylogenetic profiles

We introduce a computational method for identifying subcellular locations of proteins from the phylogenetic distribution of the homologs of organellar proteins. This method is based on the observation that proteins localized to a given organelle by experiments tend to share a characteristic phylogenetic distribution of their homologs-a phylogenetic profile. Therefore any other protein can be localized by its phylogenetic profile. Application of this method to mitochondrial proteins reveals that nucleus-encoded proteins previously known to be destined for mitochondria fall into three groups: prokaryote-derived, eukaryote-derived, and organism-specific (i.e., found only in the organism under study). Prokaryote-derived mitochondrial proteins can be identified effectively by their phylogenetic profiles. In the yeast Saccharomyces cerevisiae, 361 nucleus-encoded mitochondrial proteins can be identified at 50% accuracy with 58% coverage. From these values and the proportion of conserved mitochondrial genes, it can be inferred that approximately 630 genes, or 10% of the nuclear genome, is devoted to mitochondrial function. In the worm Caenorhabditis elegans, we estimate that there are approximately 660 nucleus-encoded mitochondrial genes, or 4% of its genome, with approximately 400 of these genes contributed from the prokaryotic mitochondrial ancestor. The large fraction of organism-specific and eukaryote-derived genes suggests that mitochondria perform specialized roles absent from prokaryotic mitochondrial ancestors. We observe measurably distinct phylogenetic profiles among proteins from different subcellular compartments, allowing the general use of prokaryotic genomes in learning features of eukaryotic proteins.     

TBC1D7 Mutations are Associated with Intellectual Disability,Macrocrania, Patellar Dislocation,and Celiac Disease

TBC1D7 forms a complex with TSC1 and TSC2 that inhibits mTORC1 signaling and limits cell growth. Mutations in TBC1D7 were reported in a family with intellectual disability (ID) and macrocrania. Using exome sequencing, we identified two sisters homozygote for the novel c.17_20delAGAG, p.R7TfsX21 TBC1D7 truncating mutation. In addition to the already described macrocephaly and mild ID, they share osteoarticular defects, patella dislocation, behavioral abnormalities, psychosis, learning difficulties, celiac disease, prognathism, myopia, and astigmatism. Consistent with a loss‐of‐function of TBC1D7, the patient's cell lines show an increase in the phosphorylation of 4EBP1, a direct downstream target of mTORC1 and a delay in the initiation of the autophagy process. This second family allows enlarging the phenotypic spectrum associated with TBC1D7 mutations and defining a TBC1D7 syndrome. Our work reinforces the involvement of TBC1D7 in the regulation of mTORC1 pathways and suggests an altered control of autophagy as possible cause of this disease.     

The genome of the fire ant Solenopsis invicta

Ants have evolved very complex societies and are key ecosystem members. Some ants, such as the fire ant Solenopsis invicta, are also major pests. Here, we present a draft genome of S. invicta, assembled from Roche 454 and Illumina sequencing reads obtained from a focal haploid male and his brothers. We used comparative genomic methods to obtain insight into the unique features of the S. invicta genome. For example, we found that this genome harbors four adjacent copies of vitellogenin. A phylogenetic analysis revealed that an ancestral vitellogenin gene first underwent a duplication that was followed by possibly independent duplications of each of the daughter vitellogenins. The vitellogenin genes have undergone subfunctionalization with queen- and worker-specific expression, possibly reflecting differential selection acting on the queen and worker castes. Additionally, we identified more than 400 putative olfactory receptors of which at least 297 are intact. This represents the largest repertoire reported so far in insects. S. invicta also harbors an expansion of a specific family of lipid-processing genes, two putative orthologs to the transformer/feminizer sex differentiation gene, a functional DNA methylation system, and a single putative telomerase ortholog. EST data indicate that this S. invicta telomerase ortholog has at least four spliceforms that differ in their use of two sets of mutually exclusive exons. Some of these and other unique aspects of the fire ant genome are likely linked to the complex social behavior of this species.     

Animal Toxins: How is Complexity Represented in Databases?

Peptide toxins synthesized by venomous animals have been extensively studied in the last decades. To be useful to the scientific community, this knowledge has been stored, annotated and made easy to retrieve by several databases. The aim of this article is to present what type of information users can access from each database. ArachnoServer and ConoServer focus on spider toxins and cone snail toxins, respectively. UniProtKB, a generalist protein knowledgebase, has an animal toxin-dedicated annotation program that includes toxins from all venomous animals. Finally, the ATDB metadatabase compiles data and annotations from other databases and provides toxin ontology.