Inferring protein function is an integral part of genome annotation and analysis. This process is usually performed in silico, and most in silico inferences are based on sequence homology approaches, which can fail when in presence of divergent sequences. However, because protein structures and their biological roles are intertwined, protein function can also be inferred by searching for structural homology. Many excellent tools have been released in recent years with regards to protein structure prediction, structural homology searches and protein visualization. Unfortunately, these tools are disconnected from each other and often use a web server-based approach that is ill-suited to high-throughput genome-wide analyses. To help assist genome annotation, we built a structural homology-based pipeline called 3DFI (for tridimensional functional inference) leveraging some of the best structural homology tools. This pipeline was built with simplicity of use in mind and enables genome-wide structural homology inferences.

Viruses acquire host genes via horizontal transfer and can express them to manipulate host biology during infections. Some homologs retain sequence identity, but evolutionary divergence can obscure host origins. We use structural modeling to compare vaccinia virus proteins with metazoan proteomes. We identify vaccinia A47L as a homolog of gasdermins, the executioners of pyroptosis. An X-ray crystal structure of A47 confirms this homology, and cell-based assays reveal that A47 interferes with caspase function. We also identify vaccinia C1L as the product of a cryptic gene fusion event coupling a Bcl-2-related fold with a pyrin domain. C1 associates with components of the inflammasome, a cytosolic innate immune sensor involved in pyroptosis, yet paradoxically enhances inflammasome activity, suggesting differential modulation during infections. Our findings demonstrate the increasing power of structural homology screens to reveal proteins with unique combinations of domains that viruses capture from host genes and combine in unique ways.

Tea (Camellia sinensis (L). O. Kuntze) is known as the oldest, mild stimulating caffeine containing non-alcoholic beverage. One of the major threats in south Asian tea industry is the blister blight leaf disease (BB), caused by the fungus Exobasidium vexans Masse. SSR DNA marker EST SSR 073 is used as a molecular marker to tag blister blight disease resistance trait of tea. The amino acid sequences were derived from cDNA sequences related to EST SSR 073 of BB susceptible (TRI 2023) and BB resistant (TRI 2043) cultivars. An attempt has been made to understand the structural characteristics and variations of EST SSR 073 locus that may reveal the factors influencing the BB resistance of tea with multiple bioinformatics tools such as ORF finder, ExPasy ProtParam tools, modeler V 9.17, Rampage server, UCSF-Chimera, and HADDOCK docking server.

Symbiosis is a complex genetic regulatory biological evolution which is highly specific pertaining to plant species and microbial strains. Biological nitrogen fixation in legumes is a functional combination of nodulation by nod genes and regulation by nif, fix genes. Three rhizobial strains (Rhizobium leguminosarum, Bradyrhizobium japonicum, and Mesorhizobium ciceri) that we considered for in silico analysis of nif A are proved to be the best isolates with respect to N2 fixing for ground nut, chick pea and soya bean (in vitro) out of 47 forest soil samples. An attempt has been made to understand the structural characteristics and variations of nif genes that may reveal the factors influencing the nitrogen fixation. The primary, secondary and tertiary structure of nif A protein was analyzed by using multiple bioinformatics tools such as chou-Fasman, GOR, ExPasy ProtParam tools, Prosa -web. Literature shows that the homology modeling of nif A protein have not been explored yet which insisted the immediate development for better understanding of nif A structure and its influence on biological nitrogen fixation. In the present predicted 3D structure, the nif A protein was analyzed by three different software tools (Phyre2, Swiss model, Modeller) and validated accordingly which can be considered as an acceptable model. However further in silico studies are suggested to determine the specific factors responsible for nitrogen fixing in the present three rhizobial strains.

GPR6 is an orphan G protein-coupled receptor that has been associated with the cannabinoid family because of its recognition of a sub-set of cannabinoid ligands. The high abundance of GPR6 in the central nervous system, along with high constitutive activity and a link to several neurodegenerative diseases make GPR6 a promising biological target. In fact, diverse research groups have demonstrated that GPR6 represents a possible target for the treatment of neurodegenerative disorders such as Parkinson's disease, Alzheimer's disease, and Huntington's disease. Several patents have claimed the use of a wide range of pyrazine derivatives as GPR6 inverse agonists for the treatment of Parkinson's disease symptoms and other dyskinesia syndromes. However, the full pharmacological importance of GPR6 has not yet been fully explored due to the lack of high potency, readily available ligands targeting GPR6. The long-term goal of the present study is to develop such ligands. In this paper, we describe our initial steps towards this goal. A human GPR6 homology model was constructed using a suite of computational techniques. This model permitted the identification of unique GPR6 structural features and the exploration of the GPR6 binding crevice. A subset of patented pyrazine analogs were docked in the resultant GPR6 inactive state model to validate the model, rationalize the structure-activity relationships from the reported patents and identify the key residues in the binding crevice for ligand recognition. We will take this structural knowledge into the next phase of GPR6 project, in which scaffold hopping will be used to design new GPR6 ligands.

The adaptive potential of plant fungal pathogens is largely governed by the gene content of a species, consisting of core and accessory genes across the pathogen isolate repertoire. To approximate the complete gene repertoire of a globally significant crop fungal pathogen, a pan genomic analysis was undertaken for Pyrenophora tritici-repentis (Ptr), the causal agent of tan (or yellow) spot disease in wheat. In this study, 15 new Ptr genomes were sequenced, assembled and annotated, including isolates from three races not previously sequenced. Together with 11 previously published Ptr genomes, a pangenome for 26 Ptr isolates from Australia, Europe, North Africa and America, representing nearly all known races, revealed a conserved core-gene content of 57 % and presents a new Ptr resource for searching natural homologues (orthologues not acquired by horizontal transfer from another species) using remote protein structural homology. Here, we identify for the first time a non-synonymous mutation in the Ptr necrotrophic effector gene ToxB, multiple copies of the inactive toxb within an isolate, a distant natural Pyrenophora homologue of a known Parastagonopora nodorum necrotrophic effector (SnTox3), and clear genomic break points for the ToxA effector horizontal transfer region. This comprehensive genomic analysis of Ptr races includes nine isolates sequenced via long read technologies. Accordingly, these resources provide a more complete representation of the species, and serve as a resource to monitor variations potentially involved in pathogenicity.

The Duck interferon gamma (DuIFN-gamma) cDNA was cloned from a phytohaemaglutinin-stimulated duck spleen cDNA library screened using a chicken IFN-gamma (ChIFN-gamma) cDNA probe. The DuIFN-gamma cDNA is 1392 nt long and shows 99% and 80% sequence identity with another cloned DuIFN-gamma cDNA, and with ChIFN-gamma cDNA, respectively. The cDNA contains a 495 bp ORF that encodes a putative 164 amino acid (AA) protein that shares 67% identity with ChIFN-gamma, but only 30-35% identity with mammalian IFN-gamma. The predicted three-dimensional (3D) structures of DuIFN-gamma and ChIFN-gamma are similar when analysed by comparative protein modelling. Culture supernatant collected from COS cells transfected with DuIFN-gamma cDNA was able to activate nitrite secretion from a chicken macrophage cell line (HD11) in a dose-dependent fashion. This activity could not be neutralised by an anti-ChIFN-gamma monoclonal antibody (Mab 85) that was able to neutralise the activity of ChIFN-gamma. Recombinant DuIFN-gamma (rDuIFN-gamma) protein was expressed in E. coli as an N-terminally His-tagged protein and was purified on a nickel affinity column. The eluted protein, which was detected as a approximately 18 kDa band with a purity of >90%, was also detected by Western blot using the anti-ChIFN-gamma monoclonal antibody (Mab 9.1). The rDuIFN-gamma was shown to activate nitrite secretion by HD11 cells in a dose-dependent fashion with a specific activity that was approximately 16-fold lower than a rChIFN-gamma control. Two rabbit antisera raised against rDuIFN-gamma were able to neutralise COS cell-expressed DuIFN-gamma activity; one of these also neutralised ChIFN-gamma activity. These findings indicate that DuIFN-gamma shares structural and functional identity with ChIFN-gamma, which is consistent with our previous results which demonstrated cross reactivity with other lymphokines from the two species.

Alkyl hydroperoxidase activity provides an important antioxidant defense for bacterial cells. The catalytic mechanism requires two peroxidases, AhpC and AhpD, where AhpD plays the role of an essential adaptor protein.

Invertebrates rely heavily on immune-like molecules with highly diversified variability so as to counteract infections. However, the mechanisms and the relationship between this variability and functionalities are not well understood. Here, we showed that the C-terminal domain of hemocyanin (HMC) from shrimp Litopenaeus vannamei contained an evolutionary conserved domain with highly variable genetic sequence, which is structurally homologous to immunoglobulin (Ig). This domain is responsible for recognizing and binding to bacteria or red blood cells, initiating agglutination and hemolysis. Furthermore, when HMC is separated into three fractions using anti-human IgM, IgG, or IgA, the subpopulation, which reacted with anti-human IgM (HMC-M), showed the most significant antimicrobial activity. The high potency of HMC-M is a consequence of glycosylation, as it contains high abundance of α-d-mannose relative to α-d-glucose and N-acetyl-d-galactosamine. Thus, the removal of these glycans abolished the antimicrobial activity of HMC-M. Our results present a comprehensive investigation of the role of HMC in fighting against infections through genetic variability and epigenetic modification.

Inside-out activation of integrins is mediated via the binding of talin and kindlin to integrin β-subunit cytoplasmic tails. The kindlin FERM domain is interrupted by a pleckstrin homology (PH) domain within its F2 subdomain. Here, we present data confirming the importance of the kindlin-1 PH domain for integrin activation and its x-ray crystal structure at a resolution of 2.1 Å revealing a C-terminal second α-helix integral to the domain but found only in the kindlin protein family. An isoform-specific salt bridge occludes the canonical phosphoinositide binding site, but molecular dynamics simulations display transient switching to an alternative open conformer. Molecular docking reveals that the opening of the pocket would enable potential ligands to bind within it. Although lipid overlay assays suggested the PH domain binds inositol monophosphates, surface plasmon resonance demonstrated weak affinities for inositol 3,4,5-triphosphate (Ins(3,4,5)P(3); K(D) ∼100 μM) and no monophosphate binding. Removing the salt bridge by site-directed mutagenesis increases the PH domain affinity for Ins(3,4,5)P(3) as measured by surface plasmon resonance and enables it to bind PtdIns(3,5)P(2) on a dot-blot. Structural comparison with other PH domains suggests that the phosphate binding pocket in the kindlin-1 PH domain is more occluded than in kindlins-2 and -3 due to its salt bridge. In addition, the apparent affinity for Ins(3,4,5)P(3) is affected by the presence of PO(4) ions in the buffer. We suggest the physiological ligand of the kindlin-1 PH domain is most likely not an inositol phosphate but another phosphorylated species.

With many genomes now sequenced, computational annotation methods to characterize genes and proteins from their sequence are increasingly important. The BioSapiens Network has developed tools to address all stages of this process, and here we review progress in the automated prediction of protein function based on protein sequence and structure.

Tospovirus is a plant-infecting genus within the family Bunyaviridae, which also includes four animal-infecting genera: Hantavirus, Nairovirus, Phlebovirus and Orthobunyavirus. Compared to these members, the structures of Tospovirus proteins still are poorly understood. Despite multiple studies have attempted to identify candidate N protein regions involved in RNA binding and protein multimerization for tospovirus using yeast two-hybrid systems (Y2HS) and site-directed mutagenesis, the tospovirus ribonucleocapsids (RNPs) remains largely uncharacterized at the molecular level and the lack of structural information prevents detailed insight into these interactions.

Viruses acquire host genes via horizontal gene transfer and can express them to manipulate host biology during infections. Some viral and host homologs retain sequence identity, but evolutionary divergence can obscure host origins. We used structural modeling to compare vaccinia virus proteins with metazoan proteomes. We identified vaccinia A47L as a homolog of gasdermins, the executioners of pyroptosis. An X-ray crystal structure of A47 confirmed this homology and cell-based assays revealed that A47 inhibits pyroptosis. We also identified vaccinia C1L as the product of a cryptic gene fusion event coupling a Bcl-2 related fold with a pyrin domain. C1 associates with components of the inflammasome, a cytosolic innate immune sensor involved in pyroptosis, yet paradoxically enhances inflammasome activity, suggesting a benefit to poxvirus replication in some circumstances. Our findings demonstrate the potential of structural homology screens to reveal genes that viruses capture from hosts and repurpose to benefit viral fitness.

The large (L) proteins of non-segmented negative stranded (NNS) RNA viruses contain the core RNA dependent RNA polymerase activity for RNA replication and transcription as well as the activities for polyadenylating and capping the mRNA transcripts and for methylating the cap structures. There is currently no structural information available for these large multi-functional proteins. Phylogenetic analyses have led to the division of the L protein primary structure into six functional domains of high conservation that are linked by variable regions. The studies in this report investigate the role of specific amino acids within domain VI of the VSV L protein, which contains a 2'-O-ribose methyltransferase (MTase) domain. We generated a structural homology model of residues 1644-1842 within domain VI based on the crystal structure determined for the known 2'-O-ribose MTase of E. coli, RrmJ. The information generated by this homology model directed us to residues structurally important for MTase activity and SAM binding. Selected residues were analyzed by site-specific mutagenesis and the mutant L proteins were assayed for their effects on RNA synthesis and cap methylation. The goal of this study was to functionally test the model in order to gain insight into the structural constraints of this region of the L protein. The data presented here revealed specific mutations that affect transcription, replication, and 5' cap methylation, many of which resulted in polymerases temperature sensitive for RNA synthesis.

Different bioinformatics methods illuminate different aspects of protein function, from specific catalytic activities to broad functional categories. Here, a triple-pronged approach to predict function for a domain of unknown function, DUF2086, is applied. Distant homology to characterised enzymes and conservation of key residues suggest an oxygenase function. Modelling indicates that the substrate is most likely a nucleic acid. Finally, genomic context analysis linking DUF2086 to DNA repair, leads to a predicted activity of oxidative demethylation of damaged bases in DNA. The newly assigned activity is sporadically present in phyla not containing near relatives of the similarly active repair protein AlkB.

Aquaporins (AQPs) facilitate the passage of small neutral polar molecules across membranes of the cell. In animals there are four distinct AQP subfamilies, whereof AQP8 homologues constitute one of the smallest subfamilies with just one member in man. AQP8 conducts water, ammonia, urea, glycerol and H2O2 through various membranes of animal cells. This passive channel has been connected to a number of phenomena, such as volume change of mitochondria, ammonia neurotoxicity, and mitochondrial dysfunction related to oxidative stress. Currently, there is no experimentally determined structure of an AQP8, hence the structural understanding of this subfamily is limited. The recently solved structure of the plant AQP, AtTIP2;1, which has structural and functional features in common with AQP8s, has opened up for construction of homology models that are likely to be more accurate than previous models.

The Msh4-Msh5 protein complex in eukaryotes is involved in stabilizing Holliday junctions and its progenitors to facilitate crossing over during Meiosis I. These functions of the Msh4-Msh5 complex are essential for proper chromosomal segregation during the first meiotic division. The Msh4/5 proteins are homologous to the bacterial mismatch repair protein MutS and other MutS homologs (Msh2, Msh3, Msh6). Saccharomyces cerevisiae msh4/5 point mutants were identified recently that show two fold reduction in crossing over, compared to wild-type without affecting chromosome segregation. Three distinct classes of msh4/5 point mutations could be sorted based on their meiotic phenotypes. These include msh4/5 mutations that have a) crossover and viability defects similar to msh4/5 null mutants; b) intermediate defects in crossing over and viability and c) defects only in crossing over. The absence of a crystal structure for the Msh4-Msh5 complex has hindered an understanding of the structural aspects of Msh4-Msh5 function as well as molecular explanation for the meiotic defects observed in msh4/5 mutations. To address this problem, we generated a structural model of the S. cerevisiae Msh4-Msh5 complex using homology modeling. Further, structural analysis tailored with evolutionary information is used to predict sites with potentially critical roles in Msh4-Msh5 complex formation, DNA binding and to explain asymmetry within the Msh4-Msh5 complex. We also provide a structural rationale for the meiotic defects observed in the msh4/5 point mutations. The mutations are likely to affect stability of the Msh4/5 proteins and/or interactions with DNA. The Msh4-Msh5 model will facilitate the design and interpretation of new mutational data as well as structural studies of this important complex involved in meiotic chromosome segregation.

Huntington's disease (HD) is an inherited progressive neurodegenerative disorder caused by a pathological CAG trinucleotide repeat expansion in the large multi-exon gene, huntingtin (HTT). Although multiple pathogenic mechanisms have been proposed for HD, there is increasing interest in the RNA processing of the HTT gene. In mammals, most multi-exon genes are alternatively spliced; however, few alternative transcripts have been described for HTT. Given the numerous protein bands detected in mouse and human brain tissue by Western blotting using anti-huntingtin antibodies, we examined whether alternative splicing of HTT may account for some of these fragments. Using RT-PCR in mouse brain, we detected two novel splice variants of Htt that lacked the 111-bp exon 29 (Htt∆ex29) or retained a 57-bp portion of intron 28 (Htt(+57)in28) via use of a cryptic splice site. The alternative transcripts were present in wild-type and homozygous Hdh(Q150/Q150) mouse brain at all ages and in all brain regions and peripheral tissues studied. Differential splicing of Htt∆ex29 was found in the cerebellum of Hdh(Q150/Q150) mice with a significant reduction in transcript levels in mutant animals. In human brain, we detected similar splice variants lacking exons 28 and 29. The ability of alternatively spliced transcripts to encode different protein isoforms with individual functions in the cell, combined with the known role of splicing in disease, renders these novel transcripts of interest in the context of HD pathogenesis.

Protein-engineering methods have been exploited to produce a surrogate system for the extracellular neurotransmitter-binding site of a heteromeric human ligand-gated ion channel, the glycine receptor. This approach circumvents two major issues: the inherent experimental difficulties in working with a membrane-bound ion channel and the complication that a heteromeric assembly is necessary to create a key, physiologically relevant binding site. Residues that form the orthosteric site in a highly stable ortholog, acetylcholine-binding protein, were selected for substitution. Recombinant proteins were prepared and characterized in stepwise fashion exploiting a range of biophysical techniques, including X-ray crystallography, married to the use of selected chemical probes. The decision making and development of the surrogate, which is termed a glycine-binding protein, are described, and comparisons are provided with wild-type and homomeric systems that establish features of molecular recognition in the binding site and the confidence that the system is suited for use in early-stage drug discovery targeting a heteromeric α/β glycine receptor.

Neurotransmitters are released at presynaptic active zones (AZs). In the fly Drosophila, monoclonal antibody (MAB) nc82 specifically labels AZs. We employ nc82 to identify Bruchpilot protein (BRP) as a previously unknown AZ component. BRP shows homology to human AZ protein ELKS/CAST/ERC, which binds RIM1 in a complex with Bassoon and Munc13-1. The C terminus of BRP displays structural similarities to multifunctional cytoskeletal proteins. During development, transcription of the bruchpilot locus (brp) coincides with neuronal differentiation. Panneural reduction of BRP expression by RNAi constructs permits a first functional characterization of this large AZ protein: larvae show reduced evoked but normal spontaneous transmission at neuromuscular junctions. In adults, we observe loss of T bars at active zones, absence of synaptic components in electroretinogram, locomotor inactivity, and unstable flight (hence "bruchpilot"-crash pilot). We propose that BRP is critical for intact AZ structure and normal-evoked neurotransmitter release at chemical synapses of Drosophila.