Search

The world is currently facing the COVID-19 pandemic, for which mild symptoms include fever and dry cough. In severe cases, it could lead to pneumonia and ultimately death in some instances. Moreover, the causative pathogen is highly contagious and there are no drugs or vaccines for it yet. The pathogen, SARS-CoV-2, is one of the human coronaviruses which was identified to infect humans first in December 2019. SARS-CoV-2 shares evolutionary relationship to other highly pathogenic viruses such as Severe Acute Respiratory Syndrome (SARS) and Middle East respiratory syndrome (MERS). We have exploited this similarity to model a target non-structural protein, NSP1, since it is implicated in the regulation of host gene expression by the virus and hijacking of host machinery. We next interrogated the capacity to repurpose around 2300 FDA-approved drugs and more than 3,00,000 small molecules of natural origin towards drug identification through virtual screening and molecular dynamics. Interestingly, we observed simple molecules like lactose, previously known anti-virals and few secondary metabolites of plants as promising hits. These herbal plants are already practiced in Ayurveda over centuries to treat respiratory problems and inflammation. Disclaimer: we would not like to recommend uptake of these small molecules for suspect COVID patients until it is approved by competent national or international authorities.

Background: Studying plasticity in gene expression in natural systems is crucial, for predicting and managing the effects of climate change on plant species. To understand the contribution of gene expression level variations to abiotic stress compensation in a Himalaya plant (Primula sikkimensis), we carried out a transplant experiment within (Ambient), and beyond (Below Ambient and Above Ambient) the altitudinal range limit of species. We sequenced nine transcriptomes (three each from each altitudinal range condition) using Illumina sequencing technology. We compared the fitness variation of transplants among three transplant conditions., Results: A large number of significantly differentially expressed genes (DEGs) between below ambient versus ambient (109) and above ambient versus ambient (85) were identified. Transcripts involved in plant growth and development were mostly up-regulated in below ambient conditions. Transcripts involved in signalling, defence, and membrane transport were mostly up-regulated in above ambient condition. Pathway analysis revealed that most of the genes involved in metabolic processes, secondary metabolism, and flavonoid biosynthesis were differentially expressed in below ambient conditions, whereas most of the genes involved in photosynthesis and plant hormone signalling were differentially expressed in above ambient conditions. In addition, we observed higher reproductive fitness in transplant individuals at below ambient condition compared to above ambient conditions; contrary to what we expect from the cold adaptive P. sikkimensis plants., Conclusions: We reveal P. sikkimensis's capacity for rapid adaptation to climate change through transcriptome variation, which may facilitate the phenotypic plasticity observed in morphological and life history traits. The genes and pathways identified provide a genetic resource for understanding the temperature stress (both the hot and cold stress) tolerance mechanism of P. sikkimensis in their natural environment.

Homologous domains embedded in multidomain proteins of different domain architectures (DA) may exhibit subtle, but important, differences in their structure and function. Here, we consider two multidomain proteins, Arf nucleotide binding site opener (ARNO) and G protein-coupled receptor kinase 2 (GRK2), which have very different DAs, but both contain pleckstrin homology (PH) domains. We analyzed the roles of residues selectively conserved in these subfamilies of PH domains from ARNO and GRK2 proteins. DA-specific residues in PH domain are found to contribute to structural and functional specialization of ARNO and GRK2 in terms of (a) specific intra- and interprotein interactions; (b) specificity for phospholipids; and (c) participation in conformational excursions, leading to various functional forms. Our approach can also be applied to subfamilies of other protein families to identify subfamily-specific residues and their specialized roles., (© 2019 The Authors. Published by FEBS Press and John Wiley & Sons Ltd.)

Substrate channeling is a mechanism for the internal transfer of hydrophobic, unstable or toxic intermediates from the active site of one enzyme to another. Such transfer has previously been described to be mediated by a hydrophobic tunnel, the use of electrostatic highways or pivoting and by conformational changes. The enzyme PaaZ is used by many bacteria to degrade environmental pollutants. PaaZ is a bifunctional enzyme that catalyzes the ring opening of oxepin-CoA and converts it to 3-oxo-5,6-dehydrosuberyl-CoA. Here we report the structures of PaaZ determined by electron cryomicroscopy with and without bound ligands. The structures reveal that three domain-swapped dimers of the enzyme form a trilobed structure. A combination of small-angle X-ray scattering (SAXS), computational studies, mutagenesis and microbial growth experiments suggests that the key intermediate is transferred from one active site to the other by a mechanism of electrostatic pivoting of the CoA moiety, mediated by a set of conserved positively charged residues.

Plants are essential facilitators of human life on planet earth. Plants play a critical functional role in mediating the quality of air, availability of food and the sustainability of agricultural resources. However, plants are in constant interaction with its environment and often hampered by various types of stresses like biotic and abiotic ones. Biotic stress is a significant reason for crop-loss and causes yield loss in the range of 31-42%, post-harvest loss due to biotic stress is in the range of 6-20%, and abiotic stress causes 6-20% of the crop damage. Recognizing the molecular factors driving plant stress-related events, and developing molecular strategies to aid plants to tolerate, resist or adapt to biotic and abiotic stress are critical for sustainable agriculture practice. In this review, we discuss how recent advances in bioinformatics, plant genomics, and data science could help to improve our understanding of plant stress biology and improve the scale of global food production. We present various areas of scientific and technological advances, such as increased availability of genomics data through whole genome sequencing that require attention. We also discuss emerging techniques including CRISPR-Cas9 based genome engineering systems to develop plant varieties that can handle combinatorial stress signals. Growing trend of converging multiple omics technologies and availability of accurate, multi-scale models of plant stress through the study of orthologs and synteny studies, would improve our knowledge of how plants perceive, respond, and manage stress to thrive as resilient crop species and thus help to reduce global food crisis., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Domain-domain interactions in multi-domain proteins play an important role in the combined function of individual domains for the overall biological activity of the protein. The functions of the tethered domains are often coupled and hence, limited numbers of domain architectures with defined folds are known in nature. Therefore, it is of interest to document the available fold-fold combinations and their preference in multi-domain proteins. Hence, we analyzed all multi-domain proteins with known structures in the protein databank and observed that only about 860 fold-fold combinations are present among them. Analyses of multi-domain proteins represented in sequence database result in recognition of 29,860 fold-fold combinations and it accounts for only 2.8% of the theoretically possible 1,036,080 (1439C2) fold-fold combinations. The observed preference for fold-fold combinations in multi-domain proteins is interesting in the context of multiple functions through structural adaptation by gene fusion.

Chemical signaling is ubiquitous and employs a variety of receptor types to detect the cacophony of molecules relevant for each living organism. Insects, our most diverse taxon, have evolved unique olfactory receptors with as little as 10% sequence identity between receptor types. We have identified a promiscuous volatile, 2-methyltetrahydro-3-furanone (coffee furanone), that elicits chemosensory and behavioral activity across multiple insect orders and receptors. In vivo and in vitro physiology showed that coffee furanone was detected by roughly 80% of the recorded neurons expressing the insect-specific olfactory receptor complex in the antenna of Drosophila melanogaster , at concentrations similar to other known, and less promiscuous, ligands. Neurons expressing specialized receptors, other chemoreceptor types, or mutants lacking the complex entirely did not respond to this compound. This indicates that coffee furanone is a promiscuous ligand for the insect olfactory receptor complex itself and did not induce non-specific cellular responses. In addition, we present homology modeling and docking studies with selected olfactory receptors that suggest conserved interaction regions for both coffee furanone and known ligands. Apart from its physiological activity, this known food additive elicits a behavioral response for several insects, including mosquitoes, flies, and cockroaches. A broad-scale behaviorally active molecule non-toxic to humans thus has significant implications for health and agriculture. Coffee furanone serves as a unique tool to unlock molecular, physiological, and behavioral relationships across this diverse receptor family and animal taxa.

Reversible phosphorylation maintained by protein kinases and phosphatases is an integral part of intracellular signalling, and phosphorylation on tyrosine is extensively utilised in higher eukaryotes. Tyrosine phosphatases are enzymes that not only scavenge phosphotyrosine but are also involved in wide range of signalling pathways. As a result, mutations in these enzymes have been implicated in the pathogenesis of several diseases like cancer, autoimmune disorders, and muscle-related diseases. The genes that harbour phosphatase domain also display diversity in co-existing domains suggesting the recruitment of the catalytic machinery in diverse pathways. We have examined the current draft of the human genome, using a combination of 3 sequence search methods and validations, and identified 101 genes encoding tyrosine phosphatase-containing gene products, agreeing with previous reports. Such gene products adopt 37 unique domain architectures (DAs), including few new ones and harbouring few co-existing domains that have not been reported before. This semi-automated computational approach for detection of gene products belonging to a particular superfamily can now be easily applied at whole genome level on other mammalian genomes and for other protein domains as well., Competing Interests: Declaration of conflicting interests:The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Olfactory receptors are important transmembrane proteins that enable organisms to perceive odours and react to them. Structural understanding of insect olfactory receptors is scarce. In this review, we discuss different transmembrane helix prediction methods, consensus methods, topology prediction methods which can enable topology prediction of these proteins. We discuss the current success rates by applying the algorithms on few G-protein coupled receptors of known structure and olfactory receptor sequences and outstanding challenges. Finally, we discuss the impact of topology prediction on biology and modeling of ORs., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Flavivirus causes arthropod-borne severe diseases that sometimes lead to the death. The Flavivirus species including Dengue virus, Zika virus and yellow fever virus are transmitted by the bite of Aedes mosquitoes. All these viral species target the people living in their respective endemic zone causing a high mortality rate. Recent studies show that immune factors present in the Ae. aegypti saliva is the hidden culprit promoting blood meal collection, suppressing host immune molecules and promoting disease establishment. This study was designed to develop a subunit vaccine using Aedes mosquito salivary proteins targeting the aforementioned Flaviviruses. Subunit vaccine was designed very precisely by combining the immunogenic B-cell epitope with CTL and HTL epitopes and also suitable adjuvant and linkers. Immunogenicity, allergenicity and physiochemical characterization were also performed for scientific validation. Molecular docking and molecular dynamics simulations studies were carried out to confirm the stable affinity between the vaccine protein (3D) and TLR3 receptor. At last, in silico cloning was executed to get the subunit vaccine restriction clone into pET28a vectro to express it in microbial expression system. Additionally, this study warrants the experimental evaluation for the validation purposes., (Copyright © 2018 Elsevier B.V. All rights reserved.)

Background: Protein-protein interactions are crucial for normal biological processes and to regulate cellular reactions that affect gene expression and function. Several previous studies have emphasized the roles of residues at the interface of protein-protein complexes in conferring stability and specificity to the complex. Interface residues in a protein are well known for their interactions with sidechain and main chain atoms with the interacting protein. However, the extent of intra-protein interactions involving interface residues in a protein-protein complex and their relative contribution in comparison to inter-protein interactions are not clearly understood. This paper probes this feature using a dataset of protein-protein complexes of known 3-D structure., Results: We have analysed a dataset of 45 transient protein-protein complex structures with at least one of the interacting proteins with a known structure available also in the unbound form. We observe that a large proportion of interface residues (1608 out of 2137 interface residues, 75%) are involved in intra and inter-protein interactions simultaneously. The amino acid propensities of such interfacial residues involved in bifurcated interactions are found to be highly similar to the general propensities to occur at protein-protein interfaces. Finally, we observe that a majority (83%) of intra-protein interactions of interface residues with bifurcated interactions, are also observed in the protein uncomplexed form., Conclusions: We have shown, to the best of our knowledge for the first time, that a vast majority of the protein-protein interface residues are involved in extensive intra-protein interactions apart from inter-protein interactions. For a majority of such interface residues the microenvironment in the tertiary structure is pre-formed and retained upon complex formation with its cognate partner during transient interactions., Reviewers: This article was reviewed by Arumay Pal and Mallur Madhusudhan.

Computational approaches to high-throughput data are gaining importance because of explosion of sequences in the post-genomic era. This explosion of sequence data creates a huge gap among the domains of sequence structure and function, since the experimental techniques to determine the structure and function are very expensive, time taking, and laborious in nature. Therefore, there is an urgent need to emphasize on the development of computational approaches in the field of biological systems. Engagement of proteins in quaternary arrangements, such as domain swapping, might be relevant for higher compatibility of such genes at stress conditions. In this study, the capacity to engage in domain swapping was predicted from mere sequence information in the whole genome of holy Basil ( Ocimum tenuiflorum ), which is well known to be an anti-stress agent. Approximately, one-fourth of the proteins of O tenuiflorum are predicted to undergo three-dimensional (3D)-domain swapping. Furthermore, function annotation was carried out on all the predicted domain-swap sequences from the O tenuiflorum and Arabidopsis thaliana for their distribution in different Pfam protein families and gene ontology (GO) terms. These domain-swapped protein sequences are associated with many Pfam protein families with a wide range of GO annotation terms. A comparative analysis of domain-swap-predicted sequences in O tenuiflorum with gene products in A thaliana reveals that around 26% (2522 sequences) are close homologues across the 2 genomes. Functional annotation of predicted domain-swapped sequences infers that predicted domain-swap sequences are involved in diverse molecular functions, such as in gene regulation of abiotic stress conditions and adaptation to different environmental niches. Finally, the positively predicted sequences of A thaliana and O tenuiflorum were also examined for their presence in stress regulome, as recorded in our STIFDB database, to check the involvement of these proteins in different abiotic stresses., Competing Interests: Declaration of conflicting interests:The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

The number of protein structures is increasing due to the individual initiatives and rapid development of structure determination techniques. Structure-based sequence alignments of distantly related proteins enable the investigation of structural, evolutionary and functional relationships between proteins and their domains leading to their common evolutionary origin. Protein Alignments organized as Structural Superfamilies (PASS2) is a database that provides such alignments of members of protein domain superfamilies of known structure and with less than 40% sequence identity. PASS2 has been continuously updated in accordance to Structural Classification of Proteins (SCOP), and now Structural Classification of Proteins - extended (SCOPe). The current update directly corresponds to SCOPe 2.06, dealing with 2006 domain superfamilies of known structure and about 14 000 domains. Alignments have been augmented by features such as hidden Markov models, highly conserved residues, structural motifs and gene ontology terms, which are available for download. In this update, we introduce the concepts of 'extreme structural outliers' and 'split superfamilies' as well., (© The Author(s) 2019. Published by Oxford University Press.)

Since proteins evolve by divergent evolution, proteins with distant homology to each other may or may not bear similar functions. Improved computational approaches are required to recognize distant homologues that are functionally similar. One of the methods of assigning function to sequences is to use profiles derived from sequences of known structure. We describe an update of the Genomic Distribution of protein structural domain Superfamilies (GenDiS) database, namely GenDiS+, which provides a projection of SCOP superfamily members on the sequence space (NR database, NCBI). The sequences are validated using structure-based sequence alignment profiles and domain and full-length sequence alignments. GenDiS+ is a `tour de force' for detecting homologues within around 160 000 taxonomic identifiers, starting from nearly 11 000 domains of known structure. Features, like full-sequence alignment and phylogeny, domain sequence alignment and phylogeny, list of associated structural and sequence domains with strength of interactions, links to databases like Pfam, UniProt and ModBase and list of sequences with a PDB structure, are provided., (© The Author(s) 2019. Published by Oxford University Press.)

[This corrects the article DOI: 10.1371/journal.pone.0100554.].

Domains are the basic building blocks of proteins which can combine to give rise to different domain architectures. Annotation of domains in a sequence is the first step towards understanding the biological function. Since there are a limited number of folds and evolutionarily related proteins have a similar structure, function can be inferred through remote homology. Computational sequence searches were performed for remote homologues on genomes of around ∼160 000 different organisms, starting from nearly 11 000 superfamily queries of known structure. Case studies revealed that most of the associated domains are involved in the same biological process. Using all the proteins predicted to have at least one structural domain, a coverage of 61% of Pfam families was achieved which is higher than the existing methods (43.36% by SIFTS). Taxonomic analysis of the proteins revealed 493 superfamilies in all the major kingdoms of life and a few lateral gene transfers between viruses and cellular organisms. The distribution of remote homologues across different classes, folds and superfamilies was studied and reveals that sequences are unequally distributed across structural classes. Finally, domain architectures were computed for the homologues and these data were compiled for each superfamily and organism.

Extensive research performed on Toll-like receptor (TLR) signaling has identified residues in the Toll/interleukin-1 receptor (TIR) domains that are essential for its proper functioning. Among these residues, those in BB loop are particularly significant as single amino acid mutations in this region can cause drastic changes in downstream signaling. However, while the effect of these mutations on the function is well studied (like the P681H mutation in TLR2, the A795P mutation in TLR3, and the P714H mutation in TLR4), their influence on the dynamics and inter-residue networks is not well understood. The effects of local perturbations induced by these mutations could propagate throughout the TIR domain, influencing interactions with other TIR domain-containing proteins. The identification of these subtle changes in inter-residue interactions can provide new insights and structural rationale for how single-point mutations cause drastic changes in TIR-TIR interactions. We employed molecular dynamics simulations and protein structure network (PSN) analyses to investigate the structural transitions with special emphasis on TLR2 and TLR3. Our results reveal that phosphorylation of the Tyr 759 residue in the TIR domain of TLR3 introduces rigidity to its BB loop. Subtle differences in the intra BB loop hydrogen bonding network between TLR3 and TLR2 are also observed. The PSN analyses indicate that the TIR domain is highly connected and pinpoints key differences in the inter-residue interactions between the wild-type and mutant TIR domains, suggesting that TIR domain structure is prone to allosteric effects, consistent with the current view of the influence of allostery on TLR signaling., (© 2018 Wiley Periodicals, Inc.)

Glycidyl ethers and their vicinal diols are important building blocks in the organic synthesis of anti-cancer and anti-obesity drugs. Ylehd, an epoxide hydrolase from tropical marine yeast Yarrowia lipolytica , was explored for its enantioselective properties by kinetic, thermodynamic and in silico studies. Kinetic resolution of racemic phenyl glycidyl ether (PGE) yielded ( S )-epoxide while for benzyl glycidyl ether (BGE) ( R )-epoxide was obtained, with vicinal diols of the opposite configuration. Amongst the enantiomers of PGE and BGE, the ( S )-selective conversion of benzyl glycidyl ether to its corresponding diol, ( S )-3-benzyloxy-1,2-propanediol while retaining ( R )-BGE was most favourable with 95% ee in 20 min. Enantioselective conversion of specific enantiomer of BGE to its corresponding diols was attributed to the favourable kinetic and thermodynamic parameters as well as to the number and proximity of water molecules near the base H325 in the active site pocket. The easily available and highly active Ylehd could be a potential biocatalyst for large scale preparation of pharmaceutically relevant chiral ( R )-benzyl glycidyl ether and ( S )-3-benzyloxy-1,2-propanediol., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

The Toll-like receptors (TLRs) are critical components of the innate immune system due to their ability to detect conserved pathogen-associated molecular patterns, present in bacteria, viruses, and other microorganisms. Ligand detection by TLRs leads to a signaling cascade, mediated by interactions among TIR domains present in the receptors, the bridging adaptors and sorting adaptors. The BB loop is a highly conserved region present in the TIR domain and is crucial for mediating interactions among TIR domain-containing proteins. Mutations in the BB loop of the Toll-like receptors, such as the A795P mutation in TLR3 and the P712H mutation (Lps d mutation) in TLR4, have been reported to disrupt or alter downstream signaling. While the phenotypic effect of these mutations is known, the underlying effect of these mutations on the structure, dynamics and interactions with other TIR domain-containing proteins is not well understood. Here, we have attempted to investigate the effect of the BB loop mutations on the dimer form of TLRs, using TLR2 and TLR3 as case studies. Our results based on molecular dynamics simulations, protein-protein interaction analyses and protein structure network analyses highlight significant differences between the dimer interfaces of the wild-type and mutant forms and provide a logical reasoning for the effect of these mutations on adaptor binding to TLRs. Furthermore, it also leads us to propose a hypothesis for the differential requirement of signaling and bridging adaptors by TLRs. This could aid in further understanding of the mechanisms governing such signaling pathways., (© 2018 Wiley Periodicals, Inc.)

Background: The eukaryotic voltage-gated sodium channel(e-Nav) is a large asymmetric transmembrane protein with important functions concerning neurological function. No structure has been resolved at high resolution for this protein., Methods: A homology model of the transmembrane and extracellular regions of an Anopheles gambiae para-like channel with emphasis on the pore entrance has been constructed, based upon the templates provided by a prokaryotic sodium channel and a potassium two-pore channel. The latter provides a template for the extracellular regions, which are located above the entrance to the pore, which is likely to open at a side of a dome formed by these loops., Results: A model created with this arrangement shows a structure similar to low-resolution cryoelectron microscope images of a related structure. The pore entrance also shows favorable electrostatic interface., Conclusion: Residues responsible for the negative charge around the pore have been traced in phylogeny to highlight their importance. This model is intended for the study of pore-blocking toxins., (Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.org.)

Point mutations in genes encoding sarcomeric proteins are the leading cause of inherited primary cardiomyopathies. Among them are mutations in the TNNT2 gene that encodes cardiac troponin T (TnT). These mutations are clustered in the tropomyosin (Tm) binding region of TnT, TNT1 (residues 80-180). To understand the mechanistic changes caused by pathogenic mutations in the TNT1 region, six hypertrophic cardiomyopathy (HCM) and two dilated cardiomyopathy (DCM) mutants were studied by biochemical approaches. Binding assays in the absence and presence of actin revealed changes in the affinity of some, but not all, TnT mutants for Tm relative to WT TnT. HCM mutants were hypersensitive and DCM mutants were hyposensitive to Ca 2+ in regulated actomyosin ATPase activities. To gain better insight into the disease mechanism, we modeled the structure of TNT1 and its interactions with Tm. The stability predictions made by the model correlated well with the affinity changes observed in vitro of TnT mutants for Tm. The changes in Ca 2+ sensitivity showed a strong correlation with the changes in binding affinity. We suggest the primary reason by which these TNNT2 mutations between residues 92 and 144 cause cardiomyopathy is by changing the affinity of TnT for Tm within the TNT1 region., Competing Interests: Conflict of interest statement: J.A.S. is a founder of and owns shares in MyoKardia, Inc., a biotechnology company developing therapeutics to treat the underlying biomechanical defects leading to hypertrophic cardiomyopathy or dilated cardiomyopathy. J.A.S. is also a founder of Cytokinetics and MyoKardia and a member of their scientific advisory boards.

Olfactory/odorant receptors (ORs) probably govern eusocial behaviour in honey bees through detection of cuticular hydrocarbons (CHCs) and queen mandibular gland pheromones (QMP). CHCs are involved in nest-mate recognition whereas QMP acts as sex pheromone for drones and as retinue pheromone for female workers. Further studies on the effect of eusociality on the evolution of ORs are hindered by the non-availability of comprehensive OR sets of solitary species. We report complete OR repertoires from two solitary bees Dufourea novaeangliae (112 ORs) and Habropoda laboriosa (151 ORs). We classify these ORs into 34 phylogenetic clades/subfamilies. Differences in the OR sets of solitary and eusocial bees are observed in individual subfamilies like subfamily 9-exon (putative CHC receptors) and L (contains putative QMP receptor group). A subfamily (H) including putative floral scent receptors is expanded in the generalist honey bees only, but not in the specialists. On the contrary, subfamily J is expanded in all bees irrespective of their degree of social complexity or food preferences. Finally, we show species-lineage specific and OR-subfamily specific differences in the putative cis-regulatory DNA motifs of the ORs from six hymenopteran species. Out of these, [A/G]CGCAAGCG[C/T] is a candidate master transcription factor binding site for multiple olfactory genes.

Background: Pathogenic bacteria have evolved various strategies to counteract host defences. They are also exposed to environments that are undergoing constant changes. Hence, in order to survive, bacteria must adapt themselves to the changing environmental conditions by performing regulations at the transcriptional and/or post-transcriptional levels. Roles of RNA-binding proteins (RBPs) as virulence factors have been very well studied. Here, we have used a sequence search-based method to compare and contrast the proteomes of 16 pathogenic and three non-pathogenic E. coli strains as well as to obtain a global picture of the RBP landscape (RBPome) in E. coli., Results: Our results show that there are no significant differences in the percentage of RBPs encoded by the pathogenic and the non-pathogenic E. coli strains. The differences in the types of Pfam domains as well as Pfam RNA-binding domains, encoded by these two classes of E. coli strains, are also insignificant. The complete and distinct RBPome of E. coli has been established by studying all known E. coli strains till date. We have also identified RBPs that are exclusive to pathogenic strains, and most of them can be exploited as drug targets since they appear to be non-homologous to their human host proteins. Many of these pathogen-specific proteins were uncharacterised and their identities could be resolved on the basis of sequence homology searches with known proteins. Detailed structural modelling, molecular dynamics simulations and sequence comparisons have been pursued for selected examples to understand differences in stability and RNA-binding., Conclusions: The approach used in this paper to cross-compare proteomes of pathogenic and non-pathogenic strains may also be extended to other bacterial or even eukaryotic proteomes to understand interesting differences in their RBPomes. The pathogen-specific RBPs reported in this study, may also be taken up further for clinical trials and/or experimental validations.

Remarkable features that are achieved in a protein-protein complex to precise levels are stability and specificity. Deviation from the normal levels of specificity and stability, which is often caused by mutations, could result in disease conditions. Chemical nature, 3-D arrangement and dynamics of interface residues code for both specificity and stability. This article reviews roles of interfacial residues in transient protein-protein complexes. It is proposed that aside from hotspot residues conferring stability to the complex, a small set of 'rigid' residues at the interface that maintain conformation between complexed and uncomplexed forms, play a major role in conferring specificity. Exceptionally, 'super hotspot' residues, which confer both stability and specificity, are attractive sites for interaction with small molecule inhibitors., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Background: TRIF is a key protein in antiviral innate immunity, operating downstream of TLRs. TRIF activation leads to the production of interferon-β and pro-inflammatory cytokines. There is evidence from experiments to suggest that the N-terminal domain of TRIF binds to its TIR domain to avoid constitutive activation. However, no structure of a complex between the N-terminal domain and the TIR domain exists till date. The disordered nature of the region connecting the N-terminal domain and the TIR domain compounds the issue of elucidating the mechanism of autoinhibition of TRIF. In this study, we have employed an integrative approach consisting of mutual information analysis, docking, molecular dynamics simulations and residue network analysis, in combination with existing experimental data to provide a glimpse of TRIF in its autoinhibited state., Results: Our extensive docking approach reveals that the N-terminal domain binds to the BB loop-B helix region of the TIR domain, consistent with experimental observations. Long length molecular dynamics simulations of 1 microsecond performed on the docked model highlights residues participating in hydrogen bonding and hydrophobic interactions at the interface. A pair of residues present in the vicinity of the interface is also predicted by mutual information analysis, to co-evolve. Residues mediating long-range interactions within the TIR domain of TRIF were identified using residue network analysis., Conclusions: Based on the results of the modelling and residue network analysis, we propose that the N-terminal domain binds to the BB loop region of the TIR domain, thereby preventing its homodimersation. The binding of TRIF to TLR3 or TRAM could induce a slight conformational change, causing the interactions between the N-terminal domain and TIR domain to disrupt, thereby exposing the BB loop and rendering it amenable for higher-order oligomerisation., Reviewers: This article was reviewed by Michael Gromiha, Srikrishna Subramaniam and Peter Bond (nominated by Chandra Verma).

Background: RNA-binding proteins (RBPs) interact with their cognate RNA(s) to form large biomolecular assemblies. They are versatile in their functionality and are involved in a myriad of processes inside the cell. RBPs with similar structural features and common biological functions are grouped together into families and superfamilies. It will be useful to obtain an early understanding and association of RNA-binding property of sequences of gene products. Here, we report a web server, RStrucFam, to predict the structure, type of cognate RNA(s) and function(s) of proteins, where possible, from mere sequence information., Results: The web server employs Hidden Markov Model scan (hmmscan) to enable association to a back-end database of structural and sequence families. The database (HMMRBP) comprises of 437 HMMs of RBP families of known structure that have been generated using structure-based sequence alignments and 746 sequence-centric RBP family HMMs. The input protein sequence is associated with structural or sequence domain families, if structure or sequence signatures exist. In case of association of the protein with a family of known structures, output features like, multiple structure-based sequence alignment (MSSA) of the query with all others members of that family is provided. Further, cognate RNA partner(s) for that protein, Gene Ontology (GO) annotations, if any and a homology model of the protein can be obtained. The users can also browse through the database for details pertaining to each family, protein or RNA and their related information based on keyword search or RNA motif search., Conclusions: RStrucFam is a web server that exploits structurally conserved features of RBPs, derived from known family members and imprinted in mathematical profiles, to predict putative RBPs from sequence information. Proteins that fail to associate with such structure-centric families are further queried against the sequence-centric RBP family HMMs in the HMMRBP database. Further, all other essential information pertaining to an RBP, like overall function annotations, are provided. The web server can be accessed at the following link: http://caps.ncbs.res.in/rstrucfam .

We developed a computational pipeline for homology based identification of the complete repertoire of olfactory receptor (OR) genes in the Asian honey bee species, Apis florea Apis florea is phylogenetically the most basal honey bee species and also the most distant sister species to the Western honey bee Apis mellifera, for which all OR genes had been identified before. Using our pipeline, we identified 180 OR genes in A. florea, which is very similar to the number of ORs identified in A. mellifera (177 ORs). Many characteristics of the ORs including gene structure, synteny of tandemly repeated ORs and basic phylogenetic clustering are highly conserved. The composite phylogenetic tree of A. florea and A. mellifera ORs could be divided into 21 clades which are in harmony with the existing Hymenopteran tree. However, we found a few nonorthologous OR relationships between both species as well as independent pseudogenization of ORs suggesting separate evolutionary changes. Particularly, a subgroup of the OR gene clade XI, which had been hypothesized to code cuticular hydrocarbon receptors showed a high number of species-specific ORs RNAseq analysis detected a total number of 145 OR transcripts in male and 162 in female antennae. Most of the OR genes were highly expressed on the female antennae. However, we detected five distinct male-biased OR genes, out of which three genes (AfOr11, AfOr18, AfOr170P) were shown to be male-biased in A. mellifera, too, thus corroborating a behavioral function in sex-pheromone communication., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: The study of survival and communication of pathogenic bacteria is important to combat diseases caused by such micro-organisms. Bacterial cells communicate with each other using a density-dependent cell-cell communication process called Quorum Sensing (QS). LuxS protein is an important member of interspecies quorum-sensing system, involved in the biosynthesis of Autoinducer-2 (AI-2), and has been identified as a drug target. Despite the above mentioned significance, their evolution has not been fully studied, particularly from a structural perspective., Results: Search for LuxS in the non-redundant database of protein sequences yielded 3106 sequences. Phylogenetic analysis of these sequences revealed grouping of sequences into five distinct clusters belonging to different phyla and according to their habitat. A majority of the neighbouring genes of LuxS have been found to be hypothetical proteins. However, gene synteny analyses in different bacterial genomes reveal the presence of few interesting gene neighbours. Moreover, LuxS gene was found to be a component of an operon in only six out of 36 genomes. Analysis of conserved motifs in representative LuxS sequences of different clusters revealed the presence of conserved motifs common to sequences of all the clusters as well as motifs unique to each cluster. Homology modelling of LuxS protein sequences of each cluster revealed few structural features unique to protein of each cluster. Analyses of surface electrostatic potentials of the homology models of each cluster showed the interactions that are common to all the clusters, as well as cluster-specific potentials and therefore interacting partners, which may be unique to each cluster., Conclusions: LuxS protein evolved early during the course of bacterial evolution, but has diverged into five subtypes. Analysis of sequence motifs and homology models of representative members reveal cluster-specific structural properties of LuxS. Further, it is also shown that LuxS protein may be involved in various protein-protein or protein-RNA interactions, which may regulate the activity of LuxS proteins in bacteria.

Accurate assessment of genetic variation in human DNA sequencing studies remains a nontrivial challenge in clinical genomics and genome informatics. Ascribing functional roles and/or clinical significances to single nucleotide variants identified from a next-generation sequencing study is an important step in genome interpretation. Experimental characterization of all the observed functional variants is yet impractical; thus, the prediction of functional and/or regulatory impacts of the various mutations using in silico approaches is an important step toward the identification of functionally significant or clinically actionable variants. The relationships between genotypes and the expressed phenotypes are multilayered and biologically complex; such relationships present numerous challenges and at the same time offer various opportunities for the design of in silico variant assessment strategies. Over the past decade, many bioinformatics algorithms have been developed to predict functional consequences of single nucleotide variants in the protein coding regions. In this review, we provide an overview of the bioinformatics resources for the prediction, annotation and visualization of coding single nucleotide variants. We discuss the currently available approaches and major challenges from the perspective of protein sequence, structure, function and interactions that require consideration when interpreting the impact of putatively functional variants. We also discuss the relevance of incorporating integrated workflows for predicting the biomedical impact of the functionally important variations encoded in a genome, exome or transcriptome. Finally, we propose a framework to classify variant assessment approaches and strategies for incorporation of variant assessment within electronic health records., (© The Author 2015. Published by Oxford University Press. For Permissions, please email: journals.permissions@oup.com.)

The Transformer2 (Tra2) proteins in humans are homologues of the Drosophila Tra2 protein. One of the two RNA-binding paralogs, Tra2β, has been very well-studied over the past decade, but not much is known about Tra2α. It was very recently shown that the two proteins demonstrate the phenomenon of paralog compensation. Here, we provide a structural basis for this genetic backup circuit, using molecular modelling and dynamics studies. We show that the two proteins display similar binding specificities, but differential affinities to a short GAA-rich RNA stretch. Starting from the 6-nucleotide RNA in the solution structure, close to 4000 virtual mutations were modelled on RNA and the domain-RNA interactions were studied after energy minimisation to convergence. Separately, another known 13-nucleotide stretch was docked and the domain-RNA interactions were observed through a 100-ns dynamics trajectory. We have also demonstrated the 'compensatory' mechanism at the level of domains in one of the domain repeat-containing RNA-binding proteins.

Myosins are actin-based motor proteins involved in many cellular movements. It is interesting to study the evolutionary patterns and the functional attributes of various types of myosins. Computational search algorithms were performed to identify putative myosin members by phylogenetic analysis, sequence motifs, and coexisting domains. This study is aimed at understanding the distribution and the likely biological functions of myosins encoded in various taxa and available eukaryotic genomes. We report here a phylogenetic analysis of around 4,064 myosin motor domains, built entirely from complete or near-complete myosin repertoires incorporating many unclassified, uncharacterized sequences and new myosin classes, with emphasis on myosins from Fungi, Haptophyta, and other Stramenopiles, Alveolates, and Rhizaria (SAR). The identification of large classes of myosins in Oomycetes, Cellular slime molds, Choanoflagellates, Pelagophytes, Eustigmatophyceae, Fonticula, Eucoccidiorida, and Apicomplexans with novel myosin motif variants that are conserved and thus presumably functional extends our knowledge of this important family of motor proteins. This work provides insights into the distribution and probable function of myosins including newly identified myosin classes.

3D-domain swapping is one of the mechanisms of protein oligomerization and the proteins exhibiting this phenomenon have many biological functions. These proteins, which undergo domain swapping, have acquired much attention owing to their involvement in human diseases, such as conformational diseases, amyloidosis, serpinopathies, proteionopathies etc. Early realisation of proteins in the whole human genome that retain tendency to domain swap will enable many aspects of disease control management. Predictive models were developed by using machine learning approaches with an average accuracy of 78% (85.6% of sensitivity, 87.5% of specificity and an MCC value of 0.72) to predict putative domain swapping in protein sequences. These models were applied to many complete genomes with special emphasis on the human genome. Nearly 44% of the protein sequences in the human genome were predicted positive for domain swapping. Enrichment analysis was performed on the positively predicted sequences from human genome for their domain distribution, disease association and functional importance based on Gene Ontology (GO). Enrichment analysis was also performed to infer a better understanding of the functional importance of these sequences. Finally, we developed hinge region prediction, in the given putative domain swapped sequence, by using important physicochemical properties of amino acids.

Protein-protein interactions play a very important role in the process of cellular functionality. Intricate details about the interactions between the proteins in a macromolecular assembly are important to understand the function and significance of protein complexes. We are reporting about a database of protein-protein interactions in huge macromolecular assemblies (PIMADb) that records the intrinsic details of 189,532 interchain interactions in 40,049 complexes from the Protein Data Bank. These details include the results of the quantification and analysis of all the interactions in the complex. The availability of interprotomer interaction networks can enable the design of point mutation experiments. PIMADb can be accessed from the URL: http://caps.ncbs.res.in/pimadb.

Insertions/deletions are common evolutionary tools employed to alter the structural and functional repertoire of protein domains. An insert situated proximal to the active site or ligand binding site frequently impacts protein function; however, the effect of distal indels on protein activity and/or stability are often not studied. In this paper, we have investigated a distal insert, which influences the function and stability of a unique DNA polymerase, called terminal deoxynucleotidyl transferase (TdT). TdT (EC:2.7.7.31) is a monomeric 58 kDa protein belonging to family X of eukaryotic DNA polymerases and known for its role in V(D)J recombination as well as in non-homologous end-joining (NHEJ) pathways. Two murine isoforms of TdT, with a length difference of twenty residues and having different biochemical properties, have been studied. All-atom molecular dynamics simulations at different temperatures and interaction network analyses were performed on the short and long-length isoforms. We observed conformational changes in the regions distal to the insert position (thumb subdomain) in the longer isoform, which indirectly affects the activity and stability of the enzyme through a mediating loop (Loop1). A structural rationale could be provided to explain the reduced polymerization rate as well as increased thermosensitivity of the longer isoform caused by peripherally located length variations within a DNA polymerase. These observations increase our understanding of the roles of length variants in introducing functional diversity in protein families in general.

Precise functioning and fine-tuning of Toll-like receptor 4 (TLR4) signaling is a critical requirement for the smooth functioning of the innate immune system, since aberrant TLR4 activation causes excessive production of pro-inflammatory cytokines and interferons. This can result in life threatening conditions such as septic shock and other inflammatory disorders. The TRIF-related adaptor molecule (TRAM) adaptor protein is unique to the TLR4 signaling pathway and abrogation of TRAM-mediated TLR4 signaling is a promising strategy for developing therapeutics aimed at disrupting TRAM interactions with other components of the TLR4 signaling complex. The VIPER motif from the vaccinia virus-producing protein, A46 has been reported to disrupt TRAM-TLR4 interactions. We have exploited this information, in combination with homology modeling and docking approaches, to identify a potential binding site on TRAM lined by the BB loop and αC helix. Virtual screening of commercially available small molecules targeting the binding site enabled to short-list 12 small molecules to abrogate TRAM-mediated TLR4 signaling. Molecular dynamics and molecular mechanics calculations have been performed for the analysis of these receptor-ligand interactions.

Differentially evolved responses to various stress conditions in plants are controlled by complex regulatory circuits of transcriptional activators, and repressors, such as transcription factors (TFs). To understand the general and condition-specific activities of the TFs and their regulatory relationships with the target genes (TGs), we have used a homogeneous stress gene expression dataset generated on ten natural ecotypes of the model plant Arabidopsis thaliana, during five single and six combined stress conditions. Knowledge-based profiles of binding sites for 25 stress-responsive TF families (187 TFs) were generated and tested for their enrichment in the regulatory regions of the associated TGs. Condition-dependent regulatory sub-networks have shed light on the differential utilization of the underlying network topology, by stress-specific regulators and multifunctional regulators. The multifunctional regulators maintain the core stress response processes while the transient regulators confer the specificity to certain conditions. Clustering patterns of transcription factor binding sites (TFBS) have reflected the combinatorial nature of transcriptional regulation, and suggested the putative role of the homotypic clusters of TFBS towards maintaining transcriptional robustness against cis-regulatory mutations to facilitate the preservation of stress response processes. The Gene Ontology enrichment analysis of the TGs reflected sequential regulation of stress response mechanisms in plants., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Design of proteins has far-reaching potentials in diverse areas that span repurposing of the protein scaffold for reactions and substrates that they were not naturally meant for, to catching a glimpse of the ephemeral proteins that nature might have sampled during evolution. These non-natural proteins, either in synthesized or virtual form have opened the scope for the design of entities that not only rival their natural counterparts but also offer a chance to visualize the protein space continuum that might help to relate proteins and understand their associations. Here, we review the recent advances in protein engineering and design, in multiple areas, with a view to drawing attention to their future potential., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

The modern sequencing technologies are generating large volumes of information at the transcriptome and genome level. Translation of this information into a biological meaning is far behind the race due to which a significant portion of proteins discovered remain as proteins of unknown function (PUFs). Attempts to uncover the functional significance of PUFs are limited due to lack of easy and high throughput functional annotation tools. Here, we report an approach to assign putative functions to PUFs, identified in the transcriptome of mulberry, a perennial tree commonly cultivated as host of silkworm. We utilized the mulberry PUFs generated from leaf tissues exposed to drought stress at whole plant level. A sequence and structure based computational analysis predicted the probable function of the PUFs. For rapid and easy annotation of PUFs, we developed an automated pipeline by integrating diverse bioinformatics tools, designated as PUFs Annotation Server (PUFAS), which also provides a web service API (Application Programming Interface) for a large-scale analysis up to a genome. The expression analysis of three selected PUFs annotated by the pipeline revealed abiotic stress responsiveness of the genes, and hence their potential role in stress acclimation pathways. The automated pipeline developed here could be extended to assign functions to PUFs from any organism in general. PUFAS web server is available at http://caps.ncbs.res.in/pufas/ and the web service is accessible at http://capservices.ncbs.res.in/help/pufas.

Motivation: In the post-genomic era, automatic annotation of protein sequences using computational homology-based methods is highly desirable. However, often protein sequences diverge to an extent where detection of homology and automatic annotation transfer is not straightforward. Sophisticated approaches to detect such distant relationships are needed. We propose a new approach to identify deep evolutionary relationships of proteins to overcome shortcomings of the available methods., Results: We have developed a method to identify remote homologues more effectively from any protein sequence database by using several cascading events with Hidden Markov Models (C-HMM). We have implemented clustering of hits and profile generation of hit clusters to effectively reduce the computational timings of the cascaded sequence searches. Our C-HMM approach could cover 94, 83 and 40% coverage at family, superfamily and fold levels, respectively, when applied on diverse protein folds. We have compared C-HMM with various remote homology detection methods and discuss the trade-offs between coverage and false positives., Availability and Implementation: A standalone package implemented in Java along with a detailed documentation can be downloaded from https://github.com/RSLabNCBS/C-HMM SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online., Contact: mini@ncbs.res.in., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

RNA-binding proteins (RBPs) are involved in various post-transcriptional gene regulatory processes and are also functionally important members of the ribosome and the spliceosome. However, RBPs and their interactions with RNA are less well-studied in comparison to DNA-binding proteins. We have classified the existing RBP structures, available in complexes with RNA and RNA/DNA hybrids, into different structural families and created Hidden Markov Models (HMMs). These structure-centric family HMMs, along with the sequence-centric family HMMs, were used as a primary database to systematically search the human proteome for the presence of putative RBPs. We have found more than 2600 gene products with RBP signatures in humans, of which around 28% are likely to bind to RNA but not DNA, whereas 9% might bind to both RNA and DNA. 11% of them do not contain an explicit functional annotation yet. Nearly 30% of the putative RBPs are exclusively nuclear, 15% have known disease associations and around 30% are enzymes. Around 40% of the proteins identified in this study are novel and have not been reported by recent large-scale studies on human RBPs.

Protein-protein interactions are essential for the basic biological machinery of the cell. This is important for processes like protein synthesis, enzyme kinetics, molecular assembly and signal transduction. A high number of macromolecular structural complexes are known due to recent advances in structure determination techniques. Therefore, it is of interest to develop an interactive tool to objectively analyze large protein complexes. Hence, we describe the development and utility of a web enabled application named 'Protein-Protein Interaction in Macro-molecular Assembly' (PIMA) for the analysis of large protein assemblies. The intricate details of physical interactions amongst protein subunits in a large complex are presented as simple user preferred interactive network diagrams.

Structure-based sequence alignment is an essential step in assessing and analysing the relationship of distantly related proteins. PASS2 is a database that records such alignments for protein domain superfamilies and has been constantly updated periodically. This update of the PASS2 version, named as PASS2.5, directly corresponds to the SCOPe 2.04 release. All SCOPe structural domains that share less than 40% sequence identity, as defined by the ASTRAL compendium of protein structures, are included. The current version includes 1977 superfamilies and has been assembled utilizing the structure-based sequence alignment protocol. Such an alignment is obtained initially through MATT, followed by a refinement through the COMPARER program. The JOY program has been used for structural annotations of such alignments. In this update, we have automated the protocol and focused on inclusion of new features such as mapping of GO terms, absolutely conserved residues among the domains in a superfamily and inclusion of PDBs, that are absent in SCOPe 2.04, using the HMM profiles from the alignments of the superfamily members and are provided as a separate list. We have also implemented a more user-friendly manner of data presentation and options for downloading more features. PASS2.5 version is available at http://caps.ncbs.res.in/pass2/., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The reelin gene is conserved across many vertebrate species, including humans. The protein product of this gene plays several important roles in early brain development and regulation of neural network plasticity of a matured brain structure. With an extended structure of 3461 amino acid sequences, consisting of eight reelin repeats, the human reelin sequence stands out as an exceptional model for evolutionary studies. In this study, sequence analysis of the human reelin and its homologues and reelin sequences from 104 other species is described in detail. Interesting sequence conservation patterns of individual repeats have been highlighted. Sequence phylogeny of the reelin sequences indicates a pattern similar to the evolution of the species, thereby serving as a highly conserved family for evolutionary purposes. Multiple sequence alignment of different reelin domain repeats, derived from homologues, suggests specific functions for individual repeats and high sequence conservation across reelin repeats from different organisms, albeit with few unusual domain architectures. A three-dimensional structural model of the full-length human reelin is now available that provides clues on residues at the dimer interface.

Background: Modeling protein-protein interactions (PPIs) using docking algorithms is useful for understanding biomolecular interactions and mechanisms. Typically, a docking algorithm generates a large number of docking poses, and it is often challenging to select the best native-like pose. A further challenge is to recognize key residues, termed as hotspots, at protein-protein interfaces, which contribute more in stabilizing a protein-protein interface., Results: We had earlier developed a computer algorithm, called PPCheck, which ascribes pseudoenergies to measure the strength of PPIs. Native-like poses could be successfully identified in 27 out of 30 test cases, when applied on a separate set of decoys that were generated using FRODOCK. PPCheck, along with conservation and accessibility scores, was able to differentiate 'native-like and non-native-like poses from 1883 decoys of Critical Assessment of Prediction of Interactions (CAPRI) targets with an accuracy of 60%. PPCheck was trained on a 10-fold mixed dataset and tested on a 10-fold mixed test set for hotspot prediction. We obtain an accuracy of 72%, which is in par with other methods, and a sensitivity of 59%, which is better than most existing methods available for hotspot prediction that uses similar datasets. Other relevant tests suggest that PPCheck can also be reliably used to identify conserved residues in a protein and to perform computational alanine scanning., Conclusions: PPCheck webserver can be successfully used to differentiate native-like and non-native-like docking poses, as generated by docking algorithms. The webserver can also be a convenient platform for calculating residue conservation, for performing computational alanine scanning, and for predicting protein-protein interface hotspots. While PPCheck can differentiate the generated decoys into native-like and non-native-like decoys with a fairly good accuracy, the results improve dramatically when features like conservation and accessibility are included. The method can be successfully used in ranking/scoring the decoys, as obtained from docking algorithms.

Background: Krishna Tulsi, a member of Lamiaceae family, is a herb well known for its spiritual, religious and medicinal importance in India. The common name of this plant is 'Tulsi' (or 'Tulasi' or 'Thulasi') and is considered sacred by Hindus. We present the draft genome of Ocimum tenuiflurum L (subtype Krishna Tulsi) in this report. The paired-end and mate-pair sequence libraries were generated for the whole genome sequenced with the Illumina Hiseq 1000, resulting in an assembled genome of 374 Mb, with a genome coverage of 61 % (612 Mb estimated genome size). We have also studied transcriptomes (RNA-Seq) of two subtypes of O. tenuiflorum, Krishna and Rama Tulsi and report the relative expression of genes in both the varieties., Results: The pathways leading to the production of medicinally-important specialized metabolites have been studied in detail, in relation to similar pathways in Arabidopsis thaliana and other plants. Expression levels of anthocyanin biosynthesis-related genes in leaf samples of Krishna Tulsi were observed to be relatively high, explaining the purple colouration of Krishna Tulsi leaves. The expression of six important genes identified from genome data were validated by performing q-RT-PCR in different tissues of five different species, which shows the high extent of urosolic acid-producing genes in young leaves of the Rama subtype. In addition, the presence of eugenol and ursolic acid, implied as potential drugs in the cure of many diseases including cancer was confirmed using mass spectrometry., Conclusions: The availability of the whole genome of O.tenuiflorum and our sequence analysis suggests that small amino acid changes at the functional sites of genes involved in metabolite synthesis pathways confer special medicinal properties to this herb.

Olfactory receptors (ORs) belong to the class A G-Protein Coupled Receptor superfamily of proteins. Unlike G-Protein Coupled Receptors, ORs exhibit a combinatorial response to odors/ligands. ORs display an affinity towards a range of odor molecules rather than binding to a specific set of ligands and conversely a single odorant molecule may bind to a number of olfactory receptors with varying affinities. The diversity in odor recognition is linked to the highly variable transmembrane domains of these receptors. The purpose of this study is to decode the odor-olfactory receptor interactions using in silico docking studies. In this study, a ligand (odor molecules) dataset of 125 molecules was used to carry out in silico docking using the GLIDE docking tool (SCHRODINGER Inc Pvt LTD). Previous studies, with smaller datasets of ligands, have shown that orthologous olfactory receptors respond to similarly-tuned ligands, but are dramatically different in their efficacy and potency. Ligand docking results were applied on homologous pairs (with varying sequence identity) of ORs from human and mouse genomes and ligand binding residues and the ligand profile differed among such related olfactory receptor sequences. This study revealed that homologous sequences with high sequence identity need not bind to the same/ similar ligand with a given affinity. A ligand profile has been obtained for each of the 20 receptors in this analysis which will be useful for expression and mutation studies on these receptors.

During 11-12 August 2014, a Protein Bioinformatics and Community Resources Retreat was held at the Wellcome Trust Genome Campus in Hinxton, UK. This meeting brought together the principal investigators of several specialized protein resources (such as CAZy, TCDB and MEROPS) as well as those from protein databases from the large Bioinformatics centres (including UniProt and RefSeq). The retreat was divided into five sessions: (1) key challenges, (2) the databases represented, (3) best practices for maintenance and curation, (4) information flow to and from large data centers and (5) communication and funding. An important outcome of this meeting was the creation of a Specialist Protein Resource Network that we believe will improve coordination of the activities of its member resources. We invite further protein database resources to join the network and continue the dialogue.

We have developed an integrated database for Mycobacterium tuberculosis H37Rv (Mtb) that collates information on protein sequences, domain assignments, functional annotation and 3D structural information along with protein-protein and protein-small molecule interactions. SInCRe (Structural Interactome Computational Resource) is developed out of CamBan (Cambridge and Bangalore) collaboration. The motivation for development of this database is to provide an integrated platform to allow easily access and interpretation of data and results obtained by all the groups in CamBan in the field of Mtb informatics. In-house algorithms and databases developed independently by various academic groups in CamBan are used to generate Mtb-specific datasets and are integrated in this database to provide a structural dimension to studies on tuberculosis. The SInCRe database readily provides information on identification of functional domains, genome-scale modelling of structures of Mtb proteins and characterization of the small-molecule binding sites within Mtb. The resource also provides structure-based function annotation, information on small-molecule binders including FDA (Food and Drug Administration)-approved drugs, protein-protein interactions (PPIs) and natural compounds that bind to pathogen proteins potentially and result in weakening or elimination of host-pathogen protein-protein interactions. Together they provide prerequisites for identification of off-target binding., (© The Author(s) 2015. Published by Oxford University Press.)

Domain swapping mechanism is a specialised mode of oligomerization of proteins in which part of a protein is exchanged in a non-covalent manner between constituent subunits. This mechanism is highly affected by several physiological conditions. Here, we present a detailed analysis ofthe effect of pH on different regions of the domain swapped oligomer by considering examples which are known to be sensitive to pH in transiting from monomeric to domain-swapped dimeric form. The energetic calculations were performed using a specialized method which considers changes in pH and subsequent changes in the interactions between subunits. This analysis provides definitive hints about the pH-dependence switch from monomer to domain-swapped oligomer and the steps that may be involved in the swapping mechanism.