The alteration of structural network upon transient association between proteins studied using graph theory

Proteins such as enzymes perform their function by predominant non-covalent bond interactions between transiently interacting units. There is an impact on the overall structural topology of the protein, albeit transient nature of such interactions, that enable proteins to deactivate or activate. This aspect of the alteration of the structural topology is studied by employing protein structural networks, which are node-edge representative models of protein structure, reported as a robust tool for capturing interactions between residues. Several methods have been optimized to collect meaningful, functionally relevant information by studying alteration of structural networks. In this article, different methods of comparing protein structural networks are employed, along with spectral decomposition of graphs to study the subtle impact of protein-protein interactions. A detailed analysis of the structural network of interacting partners is performed across a dataset of around 900 pairs of bound complexes and corresponding unbound protein structures. The variation in network parameters at, around, and far away from the interface are analyzed. Finally, we present interesting case studies, where an allosteric mechanism of structural impact is understood from communication-path detection methods. The results of this analysis are beneficial in understanding protein stability, for future engineering, and docking studies.

Autonomous face mask detection using single shot multibox detector, and ResNet-50 with identity retrieval through face matching using deep siamese neural network

The COVID-19 pandemic poses a global health challenge. The World Health Organization states that face masks are proven to be effective, especially in public areas. Real-time monitoring of face masks is challenging and exhaustive for humans. To reduce human effort and to provide an enforcement mechanism, an autonomous system has been proposed to detect non-masked people and retrieve their identity using computer vision. The proposed method introduces a novel and efficient method that involves fine-tuning the pre-trained ResNet-50 model with a new head layer for classification between masked and non-masked people. The classifier is trained using adaptive momentum optimization algorithm with decaying learning rate and binary cross-entropy loss. Data augmentation and dropout regularization are employed to achieve best convergence. During real-time application of our classifier on videos, a Caffe face detector model based on Single Shot MultiBox Detector is used to extract the face regions of interest from each frame, on which the trained classifier is applied for detecting the non-masked people. The faces of these people are then captured, which is passed on to a deep siamese neural network, based on VGG-Face model for face matching. The captured faces are compared with the reference images from the database, by extracting the features and calculating cosine distance. If the faces match, the details of that person are retrieved from the database and displayed on the web application. The proposed method has secured best results where the trained classifier has achieved 99.74% accuracy, and the identity retrieval model achieved 98.24% accuracy.

Computational tools to study RNA-protein complexes

RNA is the key player in many cellular processes such as signal transduction, replication, transport, cell division, transcription, and translation. These diverse functions are accomplished through interactions of RNA with proteins. However, protein–RNA interactions are still poorly derstood in contrast to protein–protein and protein–DNA interactions. This knowledge gap can be attributed to the limited availability of protein-RNA structures along with the experimental difficulties in studying these complexes. Recent progress in computational resources has expanded the number of tools available for studying protein-RNA interactions at various molecular levels. These include tools for predicting interacting residues from primary sequences, modelling of protein-RNA complexes, predicting hotspots in these complexes and insights into derstanding in the dynamics of their interactions. Each of these tools has its strengths and limitations, which makes it significant to select an optimal approach for the question of interest. Here we present a mini review of computational tools to study different aspects of protein-RNA interactions, with focus on overall application, development of the field and the future perspectives.

Automated Multimodal Fusion Based Hyperparameter Tuned Deep Learning Model for Brain Tumor Diagnosis

As medical image processing research has progressed, image fusion has emerged as a realistic solution, automatically extracting relevant data from many images before fusing them into a single, unified image. Medical imaging techniques, such as Computed Tomography (CT), Magnetic Resonance Imaging (MRI), etc., play a crucial role in the diagnosis and classification of brain tumors (BT). A single imaging technique is not sufficient for correct diagnosis of the disease. In case the scans are ambiguous, it can lead doctors to incorrect diagnoses, which can be unsafe to the patient. The solution to this problem is fusing images from different scans containing complementary information to generate accurate images with minimum uncertainty. This research presents a novel method for the automated identification and classification of brain tumors using multi-modal deep learning (AMDL-BTDC). The proposed AMDL-BTDC model initially performs image pre-processing using bilateral filtering (BF) technique. Next, feature vectors are generated using a pair of pre-trained deep learning models called EfficientNet and SqueezeNet. Slime Mold Algorithm is used to acquire the DL models’ optimal hyperparameter settings (SMA). In the end, an autoencoder (AE) model is used for BT classification once features have been fused. The suggested model’s superior performance over other techniques under diverse measures was validated by extensive testing on the benchmark medical imaging dataset.

Sequence Divergence and Functional Specializations of the Ancient Spliceosomal SF3b: Implications in Flexibility and Adaptations of the Multi-Protein Complex

Multi-protein assemblies are complex molecular systems that perform highly sophisticated biochemical functions in an orchestrated manner. They are subject to changes that are governed by the evolution of individual components. We performed a comparative analysis of the ancient and functionally conserved spliceosomal SF3b complex, to recognize molecular signatures that contribute to sequence divergence and functional specializations. For this, we recognized homologous sequences of individual SF3b proteins distributed across 10 supergroups of eukaryotes and identified all seven protein components of the complex in 578 eukaryotic species. Using sequence and structural analysis, we establish that proteins occurring on the surface of the SF3b complex harbor more sequence variation than the proteins that lie in the core. Further, we show through protein interface conservation patterns that the extent of conservation varies considerably between interacting partners. When we analyze phylogenetic distributions of individual components of the complex, we find that protein partners that are known to form independent subcomplexes are observed to share similar profiles, reaffirming the link between differential conservation of interface regions and their inter-dependence. When we extend our analysis to individual protein components of the complex, we find taxa-specific variability in molecular signatures of the proteins. These trends are discussed in the context of proline-rich motifs of SF3b4, functional and drug binding sites of SF3b1. Further, we report key protein-protein interactions between SF3b1 and SF3b6 whose presence is observed to be lineage-specific across eukaryotes. Together, our studies show the association of protein location within the complex and subcomplex formation patterns with the sequence conservation of SF3b proteins. In addition, our study underscores evolutionarily flexible elements that appear to confer adaptive features in individual components of the multi-protein SF3b complexes and may contribute to its functional adaptability.

Identification, functional characterization, assembly and structure of ToxIN type III toxin-antitoxin complex from E. coli

Toxin–antitoxin (TA) systems are proposed to play crucial roles in bacterial growth under stress conditions such as phage infection. The type III TA systems consist of a protein toxin whose activity is inhibited by a noncoding RNA antitoxin. The toxin is an endoribonuclease, while the antitoxin consists of multiple repeats of RNA. The toxin assembles with the individual antitoxin repeats into a cyclic complex in which the antitoxin forms a pseudoknot structure. While structure and functions of some type III TA systems are characterized, the complex assembly process is not well understood. Using bioinformatics analysis, we have identified type III TA systems belonging to the ToxIN family across different Escherichia coli strains and found them to be clustered into at least five distinct clusters. Furthermore, we report a 2.097 Å resolution crystal structure of the first E. coli ToxIN complex that revealed the overall assembly of the protein-RNA complex. Isothermal titration calorimetry experiments showed that toxin forms a high-affinity complex with antitoxin RNA resulting from two independent (5′ and 3′ sides of RNA) RNA binding sites on the protein. These results further our understanding of the assembly of type III TA complexes in bacteria.

Profiles of Natural and Designed Protein-Like Sequences Effectively Bridge Protein Sequence Gaps: Implications in Distant Homology Detection

Sequence-based approaches are fundamental to guide experimental investigations in obtaining structural and/or functional insights into uncharacterized protein families. Powerful profile-based sequence search methods rely on a sequence space continuum to identify non-trivial relationships through homology detection. The computational design of protein-like sequences that serve as "artificial linkers" is useful in identifying relationships between distant members of a structural fold. Such sequences act as intermediates and guide homology searches between distantly related proteins. Here, we describe an approach that represents natural intermediate sequences and designed protein-like sequences as HMM (Hidden Markov Models) profiles, to improve the sensitivity of existing search methods. Searches made within the "Profile database" were shown to recognize the parent structural fold for 90% of the search queries at query coverage better than 60%. For 1040 protein families with no available structure, fold associations were made through searches in the database of natural and designed sequence profiles. Most of the associations were made with the Alpha-alpha superhelix, Transmembrane beta-barrels, TIM barrel, and Immunoglobulin-like beta-sandwich folds. For 11 domain families of unknown functions, we provide confident fold associations using the profiles of designed sequences and a consensus from other fold recognition methods. For two DUFs (Domain families of Unknown Functions), we performed detailed functional annotation through comparisons with characterized templates of families of known function.

Comparative Analysis of Structural and Dynamical Features of Ribosome Upon Association With mRNA Reveals Potential Role of Ribosomal Proteins

Ribosomes play a critical role in maintaining cellular proteostasis. The binding of messenger RNA (mRNA) to the ribosome regulates kinetics of protein synthesis. To generate an understanding of the structural, mechanistic, and dynamical features of mRNA recognition in the ribosome, we have analysed mRNA-protein interactions through a structural comparison of the ribosomal complex in the presence and absence of mRNA. To do so, we compared the 3-Dimensional (3D) structures of components of the two assembly structures and analysed their structural differences because of mRNA binding, using elastic network models and structural network-based analysis. We observe that the head region of 30S ribosomal subunit undergoes structural displacement and subunit rearrangement to accommodate incoming mRNA. We find that these changes are observed in proteins that lie far from the mRNA-protein interface, implying allostery. Further, through perturbation response scanning, we show that the proteins S13, S19, and S20 act as universal sensors that are sensitive to changes in the inter protein network, upon binding of 30S complex with mRNA and other initiation factors. Our study highlights the significance of mRNA binding in the ribosome complex and identifies putative allosteric sites corresponding to alterations in structure and/or dynamics, in regions away from mRNA binding sites in the complex. Overall, our work provides fresh insights into mRNA association with the ribosome, highlighting changes in the interactions and dynamics of the ribosome assembly because of the binding.

Rewards of divergence in sequences, 3-D structures and dynamics of yeast and human spliceosome SF3b complexes

The evolution of homologous and functionally equivalent multiprotein assemblies is intriguing considering sequence divergence of constituent proteins. Here, we studied the implications of protein sequence divergence on the structure, dynamics and function of homologous yeast and human SF3b spliceosomal subcomplexes. Human and yeast SF3b comprise of 7 and 6 proteins respectively, with all yeast proteins homologous to their human counterparts at moderate sequence identity. SF3b6, an additional component in the human SF3b, interacts with the N-terminal extension of SF3b1 while the yeast homologue Hsh155 lacks the equivalent region. Through detailed homology studies, we show that SF3b6 is absent not only in yeast but in multiple lineages of eukaryotes implying that it is critical in specific organisms. We probed for the potential role of SF3b6 in the spliceosome assembled form through structural and flexibility analyses. By analysing normal modes derived from anisotropic network models of SF3b1, we demonstrate that when SF3b1 is bound to SF3b6, similarities in the magnitude of residue motions (0.86) and inter-residue correlated motions (0.94) with Hsh155 are significantly higher than when SF3b1 is considered in isolation (0.21 and 0.89 respectively). We observed that SF3b6 promotes functionally relevant 'open-to-close' transition in SF3b1 by enhancing concerted residue motions. Such motions are found to occur in the Hsh155 without SF3b6. The presence of SF3b6 influences motions of 16 residues that interact with U2 snRNA/branchpoint duplex and supports the participation of its interface residues in long-range communication in the SF3b1. These results advocate that SF3b6 potentially acts as an allosteric regulator of SF3b1 for BPS selection and might play a role in alternative splicing. Furthermore, we observe variability in the relative orientation of SF3b4 and in the local structure of three β-propeller domains of SF3b3 with reference to their yeast counterparts. Such differences influence the inter-protein interactions of SF3b between these two organisms. Together, our findings highlight features of SF3b evolution and suggests that the human SF3b may have evolved sophisticated mechanisms to fine tune its molecular function.

Signatures of conserved and unique molecular features in Afrotheria

Afrotheria is a clade of African-origin species with striking dissimilarities in appearance and habitat. In this study, we compared whole proteome sequences of six Afrotherian species to obtain a broad viewpoint of their underlying molecular make-up, to recognize potentially unique proteomic signatures. We find that 62% of the proteomes studied here, predominantly involved in metabolism, are orthologous, while the number of homologous proteins between individual species is as high as 99.5%. Further, we find that among Afrotheria, L. africana has several orphan proteins with 112 proteins showing < 30% sequence identity with their homologues. Rigorous sequence searches and complementary approaches were employed to annotate 156 uncharacterized protein sequences and 28 species-specific proteins. For 122 proteins we predicted potential functional roles, 43 of which we associated with protein- and nucleic-acid binding roles. Further, we analysed domain content and variations in their combinations within Afrotheria and identified 141 unique functional domain architectures, highlighting proteins with potential for specialized functions. Finally, we discuss the potential relevance of highly represented protein families such as MAGE-B2, olfactory receptor and ribosomal proteins in L. africana and E. edwardii, respectively. Taken together, our study reports the first comparative study of the Afrotherian proteomes and highlights salient molecular features.

An insight into the skin glands, dermal scales and secretions of the caecilian amphibian Ichthyophis beddomei

The caecilian amphibians are richly endowed with cutaneous glands, which produce secretory materials that facilitate survival in the hostile subterranean environment. Although India has a fairly abundant distribution of caecilians, there are only very few studies on their skin and secretion. In this background, the skin of Ichthyophis beddomei from the Western Ghats of Kerala, India, was subjected to light and electron microscopic analyses. There are two types of dermal glands, mucous and granular. The mucous gland has a lumen, which is packed with a mucous. The mucous-producing cells are located around the lumen. In the granular gland, a lumen is absent; the bloated secretory cells, filling the gland, are densely packed with granules of different sizes which are elegantly revealed in TEM. There is a lining of myo-epithelial cells in the peripheral regions of the glands. Small flat disk-like dermal scales, dense with squamulae, are embedded in pockets in the dermis, distributed among the cutaneous glands. 1-4 scales of various sizes are present in each scale pocket. Scanning electron microscopic observation of the skin surface revealed numerous glandular openings. The skin gland secretions, exuded through the pores, contain fatty acids, alcohols, steroid, hydrocarbons, terpene, aldehyde and a few unknown compounds.

Artificial protein sequences enable recognition of vicinal and distant protein functional relationships

High divergence in protein sequences makes the detection of distant protein relationships through homology-based approaches challenging. Grouping protein sequences into families, through similarities in either sequence or 3-D structure, facilitates in the improved recognition of protein relationships. In addition, strategically designed protein-like sequences have been shown to bridge distant structural domain families by serving as artificial linkers. In this study, we have augmented a search database of known protein domain families with such designed sequences, with the intention of providing functional clues to domain families of unknown structure. When assessed using representative query sequences from each family, we obtain a success rate of 94% in protein domain families of known structure. Further, we demonstrate that the augmented search space enabled fold recognition for 582 families with no structural information available a priori. Additionally, we were able to provide reliable functional relationships for 610 orphan families. We discuss the application of our method in predicting functional roles through select examples for DUF4922, DUF5131, and DUF5085. Our approach also detects new associations between families that were previously not known to be related, as demonstrated through new sub-groups of the RNA polymerase domain among three distinct RNA viruses. Taken together, designed sequences-augmented search databases direct the detection of meaningful relationships between distant protein families. In turn, they enable fold recognition and offer reliable pointers to potential functional sites that may be probed further through direct mutagenesis studies.

Molecular and Structural Basis of Cross-Reactivity in M. tuberculosis Toxin-Antitoxin Systems

Mycobacterium tuberculosis genome encodes over 80 toxin-antitoxin (TA) systems. While each toxin interacts with its cognate antitoxin, the abundance of TA systems presents an opportunity for potential non-cognate interactions. TA systems mediate manifold interactions to manage pathogenicity and stress response network of the cell and non-cognate interactions may play vital roles as well. To address if non-cognate and heterologous interactions are feasible and to understand the structural basis of their interactions, we have performed comprehensive computational analyses on the available 3D structures and generated structural models of paralogous M. tuberculosis VapBC and MazEF TA systems. For a majority of the TA systems, we show that non-cognate toxin-antitoxin interactions are structurally incompatible except for complexes like VapBC15 and VapBC11, which show similar interfaces and potential for cross-reactivity. For TA systems which have been experimentally shown earlier to disfavor non-cognate interactions, we demonstrate that they are structurally and stereo-chemically incompatible. For selected TA systems, our detailed structural analysis identifies specificity conferring residues. Thus, our work improves the current understanding of TA interfaces and generates a hypothesis based on congenial binding site, geometric complementarity, and chemical nature of interfaces. Overall, our work offers a structure-based explanation for non-cognate toxin-antitoxin interactions in M. tuberculosis.

Bioinformatic and mutational studies of related toxin-antitoxin pairs in Mycobacterium tuberculosis predict and identify key functional residues

Mycobacterium tuberculosis possesses an unusually large representation of type II toxin-antitoxin (TA) systems, whose functions and targets are mostly unknown. To better understand the basis of their unique expansion and to probe putative functional similarities among these systems, here we computationally and experimentally investigated their sequence relationships. Bioinformatic and phylogenetic investigations revealed that 51 sequences of the VapBC toxin family group into paralogous sub-clusters. On the basis of conserved sequence fingerprints within paralogues, we predicted functional residues and residues at the putative TA interface that are useful to evaluate TA interactions. Substitution of these likely functional residues abolished the toxin's growth-inhibitory activity. Furthermore, conducting similarity searches in 101 mycobacterial and ∼4500 other prokaryotic genomes, we assessed the relative conservation of the M. tuberculosis TA systems and found that most TA orthologues are well-conserved among the members of the M. tuberculosis complex, which cause tuberculosis in animal hosts. We found that soil-inhabiting, free-living Actinobacteria also harbor as many as 12 TA pairs. Finally, we identified five novel putative TA modules in M. tuberculosis. For one of them, we demonstrate that overexpression of the putative toxin, Rv2514c, induces bacteriostasis and that co-expression of the cognate antitoxin Rv2515c restores bacterial growth. Taken together, our findings reveal that toxin sequences are more closely related than antitoxin sequences in M. tuberculosis Furthermore, the identification of additional TA systems reported here expands the known repertoire of TA systems in M. tuberculosis.

Mycobacterium tuberculosis Rv0366c-Rv0367c encodes a non-canonical PezAT-like toxin-antitoxin pair

Toxin-antitoxin (TA) systems are ubiquitously existing addiction modules with essential roles in bacterial persistence and virulence. The genome of Mycobacterium tuberculosis encodes approximately 79 TA systems. Through computational and experimental investigations, we report for the first time that Rv0366c-Rv0367c is a non-canonical PezAT-like toxin-antitoxin system in M. tuberculosis. Homology searches with known PezT homologues revealed that residues implicated in nucleotide, antitoxin-binding and catalysis are conserved in Rv0366c. Unlike canonical PezA antitoxins, the N-terminal of Rv0367c is predicted to adopt the ribbon-helix-helix (RHH) motif for deoxyribonucleic acid (DNA) recognition. Further, the modelled complex predicts that the interactions between PezT and PezA involve conserved residues. We performed a large-scale search in sequences encoded in 101 mycobacterial and 4500 prokaryotic genomes and show that such an atypical PezAT organization is conserved in 20 other mycobacterial organisms and in families of class Actinobacteria. We also demonstrate that overexpression of Rv0366c induces bacteriostasis and this growth defect could be restored upon co-expression of cognate antitoxin, Rv0367c. Further, we also observed that inducible expression of Rv0366c in Mycobacterium smegmatis results in decreased cell-length and enhanced tolerance against a front-line tuberculosis (TB) drug, ethambutol. Taken together, we have identified and functionally characterized a novel non-canonical TA system from M. tuberculosis.

Use of designed sequences in protein structure recognition

Background

Knowledge of the protein structure is a pre-requisite for improved understanding of molecular function. The gap in the sequence-structure space has increased in the post-genomic era. Grouping related protein sequences into families can aid in narrowing the gap. In the Pfam database, structure description is provided for part or full-length proteins of 7726 families. For the remaining 52% of the families, information on 3-D structure is not yet available. We use the computationally designed sequences that are intermediately related to two protein domain families, which are already known to share the same fold. These strategically designed sequences enable detection of distant relationships and here, we have employed them for the purpose of structure recognition of protein families of yet unknown structure.

Results

We first measured the success rate of our approach using a dataset of protein families of known fold and achieved a success rate of 88%. Next, for 1392 families of yet unknown structure, we made structural assignments for part/full length of the proteins. Fold association for 423 domains of unknown function (DUFs) are provided as a step towards functional annotation.

Conclusion

The results indicate that knowledge-based filling of gaps in protein sequence space is a lucrative approach for structure recognition. Such sequences assist in traversal through protein sequence space and effectively function as ‘linkers’, where natural linkers between distant proteins are unavailable.

Reviewers

This article was reviewed by Oliviero Carugo, Christine Orengo and Srikrishna Subramanian.

Electronic supplementary material

The online version of this article (10.1186/s13062-018-0209-6) contains supplementary material, which is available to authorized users.

Immunohistochemistry in Primary B-cell Lymphoma of Tibia Masquerading as Chronic Osteomyelitis

The diagnosis between chronic osteomyelitis, Ewing sarcoma and lymphoma often is being confusing in many occasions. As the latter two conditions are malignant, early diagnosis and interventions are crucial. We present a 28 year old male with features of chronic osteomyelitis of right tibia 2 years back then changed through Ewing sarcoma to B-cell lymphoma in the histological diagnosis. This case report highlights the difficulties that arise in diagnosing primary bone lymphomas which may masquerade as chronic osteomyelitis and hence the need for immunohistochemistry in chronic recurrent osteomyelitis.

Protein sequence design and its applications

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.

De-DUFing the DUFs: Deciphering distant evolutionary relationships of Domains of Unknown Function using sensitive homology detection methods

Background

In the post-genomic era where sequences are being determined at a rapid rate, we are highly reliant on computational methods for their tentative biochemical characterization. The Pfam database currently contains 3,786 families corresponding to “Domains of Unknown Function” (DUF) or “Uncharacterized Protein Family” (UPF), of which 3,087 families have no reported three-dimensional structure, constituting almost one-fourth of the known protein families in search for both structure and function.

Results

We applied a ‘computational structural genomics’ approach using five state-of-the-art remote similarity detection methods to detect the relationship between uncharacterized DUFs and domain families of known structures. The association with a structural domain family could serve as a start point in elucidating the function of a DUF. Amongst these five methods, searches in SCOP-NrichD database have been applied for the first time. Predictions were classified into high, medium and low- confidence based on the consensus of results from various approaches and also annotated with enzyme and Gene ontology terms. 614 uncharacterized DUFs could be associated with a known structural domain, of which high confidence predictions, involving at least four methods, were made for 54 families. These structure-function relationships for the 614 DUF families can be accessed on-line at http://proline.biochem.iisc.ernet.in/RHD_DUFS/. For potential enzymes in this set, we assessed their compatibility with the associated fold and performed detailed structural and functional annotation by examining alignments and extent of conservation of functional residues. Detailed discussion is provided for interesting assignments for DUF3050, DUF1636, DUF1572, DUF2092 and DUF659.

Conclusions

This study provides insights into the structure and potential function for nearly 20 % of the DUFs. Use of different computational approaches enables us to reliably recognize distant relationships, especially when they converge to a common assignment because the methods are often complementary. We observe that while pointers to the structural domain can offer the right clues to the function of a protein, recognition of its precise functional role is still ‘non-trivial’ with many DUF domains conserving only some of the critical residues. It is not clear whether these are functional vestiges or instances involving alternate substrates and interacting partners.

Reviewers

This article was reviewed by Drs Eugene Koonin, Frank Eisenhaber and Srikrishna Subramanian.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0069-2) contains supplementary material, which is available to authorized users.

CHEXVIS: a tool for molecular channel extraction and visualization

Background

Understanding channel structures that lead to active sites or traverse the molecule is important in the study of molecular functions such as ion, ligand, and small molecule transport. Efficient methods for extracting, storing, and analyzing protein channels are required to support such studies. Further, there is a need for an integrated framework that supports computation of the channels, interactive exploration of their structure, and detailed visual analysis of their properties.

Results

We describe a method for molecular channel extraction based on the alpha complex representation. The method computes geometrically feasible channels, stores both the volume occupied by the channel and its centerline in a unified representation, and reports significant channels. The representation also supports efficient computation of channel profiles that help understand channel properties. We describe methods for effective visualization of the channels and their profiles. These methods and the visual analysis framework are implemented in a software tool, ChExVis. We apply the method on a number of known channel containing proteins to extract pore features. Results from these experiments on several proteins show that ChExVis performance is comparable to, and in some cases, better than existing channel extraction techniques. Using several case studies, we demonstrate how ChExVis can be used to study channels, extract their properties and gain insights into molecular function.

Conclusion

ChExVis supports the visual exploration of multiple channels together with their geometric and physico-chemical properties thereby enabling the understanding of the basic biology of transport through protein channels. The ChExVis web-server is freely available at http://vgl.serc.iisc.ernet.in/chexvis/. The web-server is supported on all modern browsers with latest Java plug-in.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0545-9) contains supplementary material, which is available to authorized users.

Nodular Erythema Elevatum Diutinum Mimicking Kaposi's Sarcoma in a Human Immunodeficiency Virus Infected Patient

Erythema elevatum diutinum (EED) has been emerging as a specific Human Immunodeficiency Virus (HIV) associated dermatosis in recent times. It is an extremely rare chronic disease of unknown origin and part of the spectrum of leukocytoclastic vasculitis. We describe a case of EED simulating Kaposi's sarcoma in a 52-year-old HIV infected female patient with no previous opportunistic infections and CD4+ count of 164/mm(3). Therapy with oral dapsone (100 mg/day) for two weeks resulted in resolution of some lesions.

Mapping post-translational modifications of mammalian testicular specific histone variant TH2B in tetraploid and haploid germ cells and their implications on the dynamics of nucleosome structure

Histones regulate a variety of chromatin templated events by their post-translational modifications (PTMs). Although there are extensive reports on the PTMs of canonical histones, the information on the histone variants remains very scanty. Here, we report the identification of different PTMs, such as acetylation, methylation, and phosphorylation of a major mammalian histone variant TH2B. Our mass spectrometric analysis has led to the identification of both conserved and unique modifications across tetraploid spermatocytes and haploid spermatids. We have also computationally derived the 3-dimensional model of a TH2B containing nucleosome in order to study the spatial orientation of the PTMs identified and their effect on nucleosome stability and DNA binding potential. From our nucleosome model, it is evident that substitution of specific amino acid residues in TH2B results in both differential histone-DNA and histone-histone contacts. Furthermore, we have also observed that acetylation on the N-terminal tail of TH2B weakens the interactions with the DNA. These results provide direct evidence that, similar to somatic H2B, the testis specific histone TH2B also undergoes multiple PTMs, suggesting the possibility of chromatin regulation by such covalent modifications in mammalian male germ cells.

NrichD database: sequence databases enriched with computationally designed protein-like sequences aid in remote homology detection

NrichD (http://proline.biochem.iisc.ernet.in/NRICHD/) is a database of computationally designed protein-like sequences, augmented into natural sequence databases that can perform hops in protein sequence space to assist in the detection of remote relationships. Establishing protein relationships in the absence of structural evidence or natural ‘intermediately related sequences’ is a challenging task. Recently, we have demonstrated that the computational design of artificial intermediary sequences/linkers is an effective approach to fill naturally occurring voids in protein sequence space. Through a large-scale assessment we have demonstrated that such sequences can be plugged into commonly employed search databases to improve the performance of routinely used sequence search methods in detecting remote relationships. Since it is anticipated that such data sets will be employed to establish protein relationships, two databases that have already captured these relationships at the structural and functional domain level, namely, the SCOP database and the Pfam database, have been ‘enriched’ with these artificial intermediary sequences. NrichD database currently contains 3 611 010 artificial sequences that have been generated between 27 882 pairs of families from 374 SCOP folds. The data sets are freely available for download. Additional features include the design of artificial sequences between any two protein families of interest to the user.

Filling-in void and sparse regions in protein sequence space by protein-like artificial sequences enables remarkable enhancement in remote homology detection capability

Protein functional annotation relies on the identification of accurate relationships, sequence divergence being a key factor. This is especially evident when distant protein relationships are demonstrated only with three-dimensional structures. To address this challenge, we describe a computational approach to purposefully bridge gaps between related protein families through directed design of protein-like "linker" sequences. For this, we represented SCOP domain families, integrated with sequence homologues, as multiple profiles and performed HMM-HMM alignments between related domain families. Where convincing alignments were achieved, we applied a roulette wheel-based method to design 3,611,010 protein-like sequences corresponding to 374 SCOP folds. To analyze their ability to link proteins in homology searches, we used 3024 queries to search two databases, one containing only natural sequences and another one additionally containing designed sequences. Our results showed that augmented database searches showed up to 30% improvement in fold coverage for over 74% of the folds, with 52 folds achieving all theoretically possible connections. Although sequences could not be designed between some families, the availability of designed sequences between other families within the fold established the sequence continuum to demonstrate 373 difficult relationships. Ultimately, as a practical and realistic extension, we demonstrate that such protein-like sequences can be "plugged-into" routine and generic sequence database searches to empower not only remote homology detection but also fold recognition. Our richly statistically supported findings show that complementary searches in both databases will increase the effectiveness of sequence-based searches in recognizing all homologues sharing a common fold.

Cascaded walks in protein sequence space: use of artificial sequences in remote homology detection between natural proteins

Over the past two decades, many ingenious efforts have been made in protein remote homology detection. Because homologous proteins often diversify extensively in sequence, it is challenging to demonstrate such relatedness through entirely sequence-driven searches. Here, we describe a computational method for the generation of 'protein-like' sequences that serves to bridge gaps in protein sequence space. Sequence profile information, as embodied in a position-specific scoring matrix of multiply aligned sequences of bona fide family members, serves as the starting point in this algorithm. The observed amino acid propensity and the selection of a random number dictate the selection of a residue for each position in the sequence. In a systematic manner, and by applying a 'roulette-wheel' selection approach at each position, we generate parent family-like sequences and thus facilitate an enlargement of sequence space around the family. When generated for a large number of families, we demonstrate that they expand the utility of natural intermediately related sequences in linking distant proteins. In 91% of the assessed examples, inclusion of designed sequences improved fold coverage by 5-10% over searches made in their absence. Furthermore, with several examples from proteins adopting folds such as TIM, globin, lipocalin and others, we demonstrate that the success of including designed sequences in a database positively sensitized methods such as PSI-BLAST and Cascade PSI-BLAST and is a promising opportunity for enormously improved remote homology recognition using sequence information alone.

Perspectives of plant-associated microbes in heavy metal phytoremediation

"Phytoremediation" know-how to do-how is rapidly expanding and is being commercialized by harnessing the phyto-microbial diversity. This technology employs biodiversity to remove/contain pollutants from the air, soil and water. In recent years, there has been a considerable knowledge explosion in understanding plant-microbes-heavy metals interactions. Novel applications of plant-associated microbes have opened up promising areas of research in the field of phytoremediation technology. Various metabolites (e.g., 1-aminocyclopropane-1-carboxylic acid deaminase, indole-3-acetic acid, siderophores, organic acids, etc.) produced by plant-associated microbes (e.g., plant growth promoting bacteria, mycorrhizae) have been proposed to be involved in many biogeochemical processes operating in the rhizosphere. The salient functions include nutrient acquisition, cell elongation, metal detoxification and alleviation of biotic/abiotic stress in plants. Rhizosphere microbes accelerate metal mobility, or immobilization. Plants and associated microbes release inorganic and organic compounds possessing acidifying, chelating and/or reductive power. These functions are implicated to play an essential role in plant metal uptake. Overall the plant-associated beneficial microbes enhance the efficiency of phytoremediation process directly by altering the metal accumulation in plant tissues and indirectly by promoting the shoot and root biomass production. The present work aims to provide a comprehensive review of some of the promising processes mediated by plant-associated microbes and to illustrate how such processes influence heavy metal uptake through various biogeochemical processes including translocation, transformation, chelation, immobilization, solubilization, precipitation, volatilization and complexation of heavy metals ultimately facilitating phytoremediation.

Pharmacognostical Standardization of Tephrosia purpurea Pers Root

Wild Indigo or Purple Tephrosia or fish poison occurs throughout the Indian subcontinent. It is widely used in the treatment of inflammation, diabetes, rheumatism, asthma, diarrhoea and many other ailments. But so far the pharmacognostic standardization has not been reported for its proper identification. Hence the present study is a pharmacognosy work carried out for the root part. This may help in the identification of the plant species. A thin transverse section, powder microscopy, measurement of the dimensions of cell structures, fluorescence analysis and physico chemical parameters were conducted for the root. From the TS, the secondary xylem fibres and vessels were found to be the tissues of diagnostic importance. The xylem vessels were of two types: narrow and long; broad and short. The important characters in the powdered microscopy were vessel elements, fibres and xylem parenchyma cells. The different fluorescent light shades were obtained under short and long UV light for both powder as well as the extracts of the root. The proximate analysis values were also obtained in a satisfactory way. Combining all these data a suitable root profile for plant can be constructed which may help in the identification of quality of the plant part.

Ozonation of reactive dye, Procion red MX-5B catalyzed by metal ions

Ozonation of aqueous solution of Procion red MX-5B, a commercial azo reactive dye was investigated. The effect of various experimental variables on decolorization and degradation of the dye has been studied. Color removal was faster than organic removal. Dye mineralization was evident by the generation of sulphate, chloride, nitrate, oxalic and acetic acid during ozonation. The decolorization and organic removal followed pseudo-first-order kinetics. Among the metal ions studied, manganese catalyzed ozonation showed better decolorization and organic removal. Fourier transform infrared (FTIR) spectroscopy results indicated the disappearance of bands for aromatics and appearance of bands due to carboxylic acids in the ozonated samples. Acetyl benzoic acid, diethyl phthalate and phthalic anhydride were identified as intermediates by gas chromatography mass spectra (GC/MS). Based on these findings a tentative pathway for the degradation MX-5B was postulated.

Pulmonary strongyloidiasis

In patients with abnormalities of cell-mediated immunity, strongyloides hyperinfection syndrome may occur producing pulmonary infection that may manifest as asthma, chronic bronchitis, haemoptysis, eosinophilia and pulmonary infiltrates. We report a case of an uncontrolled asthma patient who presented with fever and haemoptysis. She had no evidence of immunosuppression and CT chest showed a lesion suspicious of malignancy, sputum cytology showed strongyloides larvae and the patient had complete recovery with treatment.

Length variations amongst protein domain superfamilies and consequences on structure and function

Background

Related protein domains of a superfamily can be specified by proteins of diverse lengths. The structural and functional implications of indels in a domain scaffold have been examined.

Methodology

In this study, domain superfamilies with large length variations (more than 30% difference from average domain size, referred as ‘length-deviant’ superfamilies and ‘length-rigid’ domain superfamilies (<10% length difference from average domain size) were analyzed for the functional impact of such structural differences. Our delineated dataset, derived from an objective algorithm, enables us to address indel roles in the presence of peculiar structural repeats, functional variation, protein-protein interactions and to examine ‘domain contexts’ of proteins tolerant to large length variations. Amongst the top-10 length-deviant superfamilies analyzed, we found that 80% of length-deviant superfamilies possess distant internal structural repeats and nearly half of them acquired diverse biological functions. In general, length-deviant superfamilies have higher chance, than length-rigid superfamilies, to be engaged in internal structural repeats. We also found that ∼40% of length-deviant domains exist as multi-domain proteins involving interactions with domains from the same or other superfamilies. Indels, in diverse domain superfamilies, were found to participate in the accretion of structural and functional features amongst related domains. With specific examples, we discuss how indels are involved directly or indirectly in the generation of oligomerization interfaces, introduction of substrate specificity, regulation of protein function and stability.

Conclusions

Our data suggests a multitude of roles for indels that are specialized for domain members of different domain superfamilies. These specialist roles that we observe and trends in the extent of length variation could influence decision making in modeling of new superfamily members. Likewise, the observed limits of length variation, specific for each domain superfamily would be particularly relevant in the choice of alignment length search filters commonly applied in protein sequence analysis.

Determination of kinetic constants of hybrid textile wastewater treatment system

The present study is related to treatment of textile wastewater in microaerophilic-aerobic hybrid reactor. The study showed the effectiveness of biological treatment of wastewater involving appropriate microorganism and suitable reactors. COD and color were reduced to 82-94%, and 99% respectively for textile wastewater. The reactor was operated at highest loading of 16.4 g COD g l(-1)d(-1) and obtained 80% COD and 72% color removal. Biokinetic models were applied to data obtained from experimental studies in continuously operated hybrid reactor. Treatment efficiencies of the reactor were investigated at different hydraulic retention times (2.3-9.1d) and organic loading rates (2.6-16.4 g COD l(-1)d(-1)). Second-order and a Stover-Kincannon models were best fitted to the hybrid column reactor. The second-order substrate removal rate constant (k(2(S))) was found as 41.44 d(-1) for hybrid reactor. Applying the modified Stover-Kincannon model to the hybrid reactor, the maximum removal rate constant (U(max)) and saturation value constant (K(B)) were found to be 212 g l(-1)d(-1) and 22.89 g l(-1)d(-1), respectively.

CUSP: an algorithm to distinguish structurally conserved and unconserved regions in protein domain alignments and its application in the study of large length variations

BACKGROUND

Distantly related proteins adopt and retain similar structural scaffolds despite length variations that could be as much as two-fold in some protein superfamilies. In this paper, we describe an analysis of indel regions that accommodate length variations amongst related proteins. We have developed an algorithm CUSP, to examine multi-membered PASS2 superfamily alignments to identify indel regions in an automated manner. Further, we have used the method to characterize the length, structural type and biochemical features of indels in related protein domains.

RESULTS

CUSP, examines protein domain structural alignments to distinguish regions of conserved structure common to related proteins from structurally unconserved regions that vary in length and type of structure. On a non-redundant dataset of 353 domain superfamily alignments from PASS2, we find that 'length- deviant' protein superfamilies show > 30% length variation from their average domain length. 60% of additional lengths that occur in indels are short-length structures (< 5 residues) while 6% of indels are > 15 residues in length. Structural types in indels also show class-specific trends.

CONCLUSION

The extent of length variation varies across different superfamilies and indels show class-specific trends for preferred lengths and structural types. Such indels of different lengths even within a single protein domain superfamily could have structural and functional consequences that drive their selection, underlying their importance in similarity detection and computational modelling. The availability of systematic algorithms, like CUSP, should enable decision making in a domain superfamily-specific manner.

Molecular beacons: a real-time polymerase chain reaction assay for detecting Escherichia coli from fresh produce and water

Molecular beacons (MBs) are oligonucleotide probes that fluoresce upon hybridization. The development of a real-time polymerase chain reaction (PCR) assay to detect the presence of Escherichia coli using these fluorogenic reporter molecules is reported. MBs were designed to recognize a 19-bp region of the uid A gene, coding for an enzyme beta-glucuronidase. The specificity of the MB-based PCR assay was evaluated for various E. coli strains as well as bacteria species that are present in nature. The capability of the assay to detect E. coli in drinking water and produce was demonstrated. Positive detection of E. coli was demonstrated when >10(1) CFU mL(-1) (colony forming unit) was present in the water samples and fresh produce after 18 h of enrichment. These assays could be carried out entirely in sealed PCR tubes, enabling rapid and semiautomated detection of E. coli in food and environmental samples.

Microbially-influenced degradation of solidified/stabilized metal waste

In the present study, a refined microbially-influenced degradation method was used to evaluate the stability of a solidified synthetic waste containing chromium salt, cement and fly ash in two different proportions. The experimental samples showed evidence of microbial growth by leaching of sulfate. Chromium leached by Thiobacillus thiooxidans from the experimental samples 'C1' (10.26% CrCl3 .6H2O; 89.74% cement) and 'FC1' (10.26% CrCl3 .6H2O; 10% fly ash; 79.74% cement), after 30 days of exposure was 14.53 mg/g and 9.53 mg/g, respectively. The corresponding concentration of chromium in the leachate was 0.189 mg/l and 0.124 mg/l, respectively, which was lower than the toxicity characteristic leaching procedure (TCLP), regulatory limit (5 mg/l). Replacement of cement by 10% fly ash in FC1 restricted the leaching of chromium more effectively. Model equations based on two shrinking core models namely, acid dissolution and bulk diffusion model, were used to analyze the kinetics of microbial degradation. Of the two approaches, the bulk diffusion model fit the data better than the acid dissolution model as indicated by the correlation coefficients of >0.97.

The antigen binding sites of various hCG monoclonal antibodies show homology to different domains of LH receptor

The common feature of receptors and antibodies against the ligand is that both display very specific, high affinity binding towards the ligand. Therefore, it can be hypothesized that the paratope of antibodies may exhibit homology with distinct domains of the receptor. By locating the hormone epitopes and determining the structure of the paratopes, it should be possible to identify the contact points between the ligand and the receptor. This hypothesis has been tested using hCG monoclonal antibodies (MAbs) recognizing different epitopes and having different effects on hormone binding and response. The beta subunit and heterodimer specific antibodies inhibited both hormone binding and response, while the alpha subunit specific antibodies inhibited response without affecting binding. The single chain fragment variables (ScFvs) produced from these antibodies also retained the properties of the parent antibodies. The amino acid sequences of the ScFvs exhibited homology to different regions of the receptor; the beta subunit specific antibody being homologous to the concave surface of the leucine rich repeats (LRR) of the receptor, particularly the concave surface of the LRRs, while the heterodimer specific antibody showed homology to the hinge region. The alpha subunit specific antibody showed homology to the transmembrane domain of the receptor. The exact locations of the epitopes of the monoclonal antibodies in the hormone molecule have also been identified. The data presented here also support the model of glycoprotein hormone-receptor interaction in which the hormone binds to the extracellular domain through the beta subunit and then the alpha subunit is brought in contact with the transmembrane domain leading to signal transduction.

Cascade PSI-BLAST web server: a remote homology search tool for relating protein domains

Owing to high evolutionary divergence, it is not always possible to identify distantly related protein domains by sequence search techniques. Intermediate sequences possess sequence features of more than one protein and facilitate detection of remotely related proteins. We have demonstrated recently the employment of Cascade PSI-BLAST where we perform PSI-BLAST for many ‘generations’, initiating searches from new homologues as well. Such a rigorous propagation through generations of PSI-BLAST employs effectively the role of intermediates in detecting distant similarities between proteins. This approach has been tested on a large number of folds and its performance in detecting superfamily level relationships is ∼35% better than simple PSI-BLAST searches. We present a web server for this search method that permits users to perform Cascade PSI-BLAST searches against the Pfam, SCOP and SwissProt databases. The URL for this server is .

Recent trends in remote homology detection: an Indian Medley

The development of remote homology detection methods is a challenging area in Bioinformatics. Sequence analysis-based approaches that address this problem have employed the use of profiles, templates and Hidden Markov Models (HMMs). These methods often face limitations due to poor sequence similarities and non-uniform sequence dispersion in protein sequence space. Search procedures are often asymmetrical due to over or under-representation of some protein families and outliers often remain undetected. Intermediate sequences that share high similarities with more than one protein can help overcome such problems. Methods such as MulPSSM and Cascade PSI-BLAST that employ intermediate sequences achieve better coverage of members in searches. Others employ peptide modules or conserved patterns of motifs or residues and are effective in overcoming dependencies on high sequence similarity to establish homology by using conserved patterns in searches. We review some of these recent methods developed in India in the recent past.