BackgroundIt had long been thought that a protein exhibits its specific function through its own specific 3D-structure under physiological conditions. However, subsequent research has shown that there are many proteins without specific 3D-structures under physiological conditions, so-called intrinsically disordered proteins (IDPs). This study presents a new technique for predicting intrinsically disordered regions in a protein, based on our average distance map (ADM) technique. The ADM technique was developed to predict compact regions or structural domains in a protein. In a protein containing partially disordered regions, a domain region is likely to be ordered, thus it is unlikely that a disordered region would be part of any domain. Therefore, the ADM technique is expected to also predict a disordered region between domains.ResultsThe results of our new technique are comparable to the top three performing techniques in the community-wide CASP10 experiment. We further discuss the case of p53, a tumor-suppressor protein, which is the most significant protein among cell cycle regulatory proteins. This protein exhibits a disordered character as a monomer but an ordered character when two p53s form a dimer.ConclusionOur technique can predict the location of an intrinsically disordered region in a protein with an accuracy comparable to the best techniques proposed so far. Furthermore, it can also predict a core region of IDPs forming definite 3D structures through interactions, such as dimerization. The technique in our study may also serve as a means of predicting a disordered region which would become an ordered structure when binding to another protein.

BackgroundAntibody, the primary effector molecule of the immune system, evolves after initial encounter with the antigen from a precursor form to a mature one to effectively deal with the antigen. Antibodies of a lineage diverge through antigen-directed isolated pathways of maturation to exhibit distinct recognition potential. In the context of evolution in immune recognition, diversity of antigen cannot be ignored. While there are reports on antibody lineage, structural perspective with respect to diverse recognition potential in a lineage has never been studied. Hence, it is crucial to evaluate how maturation leads to topological tailoring within a lineage enabling them to interact with significantly distinct antigens.ResultsA data-driven approach was undertaken for the study. Global experimental mouse and human antibody-antigen complex structures from PDB were compiled into a coherent database of germline-linked antibodies bound with distinct antigens. Structural analysis of all lineages showed variations in CDRs of both H and L chains. Observations of conformational adaptation made from analysis of static structures were further evaluated by characterizing dynamics of interaction in two lineages, mouse VH1–84 and human VH5–51. Sequence and structure analysis of the lineages explained that somatic mutations altered the geometries of individual antibodies with common structural constraints in some CDRs. Additionally, conformational landscape obtained from molecular dynamics simulations revealed that incoming pathogen led to further conformational divergence in the paratope (as observed across datasets) even while maintaining similar overall backbone topology. MM-GB/SA analysis showed binding energies to be in physiological range. Results of the study are coherent with experimental observations.ConclusionsThe findings of this study highlight basic structural principles shaping the molecular evolution of a lineage for significantly diverse antigens. Antibodies of a lineage follow different developmental pathways while preserving the imprint of the germline. From the study, it can be generalized that structural diversification of the paratope is an outcome of natural selection of a conformation from an available ensemble, which is further optimized for antigen interaction. The study establishes that starting from a common lineage, antibodies can mature to recognize a wide range of antigens. This hypothesis can be further tested and validated experimentally.

BackgroundSyphilis continues to be a major global health threat with 11 million new infections each year, and a global burden of 36 million cases. The causative agent of syphilis, Treponema pallidum subspecies pallidum, is a highly virulent bacterium, however the molecular mechanisms underlying T. pallidum pathogenesis remain to be definitively identified. This is due to the fact that T. pallidum is currently uncultivatable, inherently fragile and thus difficult to work with, and phylogenetically distinct with no conventional virulence factor homologs found in other pathogens. In fact, approximately 30% of its predicted protein-coding genes have no known orthologs or assigned functions. Here we employed a structural bioinformatics approach using Phyre2-based tertiary structure modeling to improve our understanding of T. pallidum protein function on a proteome-wide scale.ResultsPhyre2-based tertiary structure modeling generated high-confidence predictions for 80% of the T. pallidum proteome (780/978 predicted proteins). Tertiary structure modeling also inferred the same function as primary structure-based annotations from genome sequencing pipelines for 525/605 proteins (87%), which represents 54% (525/978) of all T. pallidum proteins. Of the 175 T. pallidum proteins modeled with high confidence that were not assigned functions in the previously annotated published proteome, 167 (95%) were able to be assigned predicted functions. Twenty-one of the 175 hypothetical proteins modeled with high confidence were also predicted to exhibit significant structural similarity with proteins experimentally confirmed to be required for virulence in other pathogens.ConclusionsPhyre2-based structural modeling is a powerful bioinformatics tool that has provided insight into the potential structure and function of the majority of T. pallidum proteins and helped validate the primary structure-based annotation of more than 50% of all T. pallidum proteins with high confidence. This work represents the first T. pallidum proteome-wide structural modeling study and is one of few studies to apply this approach for the functional annotation of a whole proteome.

BackgroundEarlier based on bioinformatics analyses, we had predicted the Mycobacterium tuberculosis (M.tb) proteins; Rv1555 and Rv1554, among the potential new tuberculosis drug targets. According to the ‘TB-drugome’ the Rv1555 protein is ‘druggable’ with sildenafil (Viagra), tadalafil (Cialis) and vardenafil (Levitra) drugs. In the present work, we intended to understand via computer modeling studies, how the above drugs are likely to inhibit the M.tb protein’s function.ResultsThe three-dimensional computer models for M.tb proteins; Rv1555 and Rv1554 constructed on the template of equivalent membrane anchor subunits of the homologous E.coli quinol fumarate reductase respiratory protein complex, followed by drug docking analyses, suggested that the binding of above drugs interferes with quinol binding sites. Also, we experimentally observed the in-vitro growth inhibition of E.coli bacteria containing the homologous M.tb protein sequences with sildenafil and tadalafil drugs.ConclusionsThe predicted binding sites of the drugs is likely to affect the above M.tb proteins function as quinol binding is known to be essential for electron transfer function during anaerobic respiration in the homologous E.coli protein complex. Therefore, sildenafil and related drugs currently used in the treatment of male erectile dysfunction targeting the human phosphodiesterase 5 enzyme may be evaluated for their plausible role as repurposed drugs to treat human tuberculosis.

BackgroundPost-crystallization dehydration methods, applying either vapor diffusion or humidity control devices, have been widely used to improve the diffraction quality of protein crystals. Despite the fact that RNA crystals tend to diffract poorly, there is a dearth of reports on the application of dehydration methods to improve the diffraction quality of RNA crystals.ResultsWe use dehydration techniques with a Free Mounting System (FMS, a humidity control device) to recover the poor diffraction quality of RNA crystals. These approaches were applied to RNA constructs that model various RNA-mediated repeat expansion disorders.ConclusionThe method we describe herein could serve as a general tool to improve diffraction quality of RNA crystals to facilitate structure determinations.Electronic supplementary materialThe online version of this article (doi:10.1186/s12900-016-0069-1) contains supplementary material, which is available to authorized users.

BackgroundWe comprehensively analyzed X-ray cocrystal structures of dipeptidyl peptidase IV (DPP-4) and its inhibitor to clarify whether DPP-4 alters its general or partial structure according to the inhibitor used and whether DPP-4 has a common rule for inhibitor binding.ResultsAll the main and side chains in the inhibitor binding area were minimally altered, except for a few side chains, despite binding to inhibitors of various shapes. Some residues (Arg125, Glu205, Glu206, Tyr662 and Asn710) in the area had binding modes to fix a specific atom of inhibitor to a particular spatial position in DPP-4. We found two specific water molecules that were common to 92 DPP-4 structures. The two water molecules were close to many inhibitors, and seemed to play two roles: maintaining the orientation of the Glu205 and Glu206 side chains through a network via the water molecules, and arranging the inhibitor appropriately at the S2 subsite.ConclusionsOur study based on high-quality resources may provide a necessary minimum consensus to help in the discovery of a novel DPP-4 inhibitor that is commercially useful.Electronic supplementary materialThe online version of this article (doi:10.1186/s12900-016-0062-8) contains supplementary material, which is available to authorized users.

BackgroundRNA ligases 2 are scarce and scattered across the tree of life. Two members of this family are well studied: the mitochondrial RNA editing ligase from the parasitic trypanosomes (Kinetoplastea), a promising drug target, and bacteriophage T4 RNA ligase 2, a workhorse in molecular biology. Here we report the identification of a divergent RNA ligase 2 (DpRNL) from Diplonema papillatum (Diplonemea), a member of the kinetoplastids’ sister group.MethodsWe identified DpRNL with methods based on sensitive hidden Markov Model. Then, using homology modeling and molecular dynamics simulations, we established a three dimensional structure model of DpRNL complexed with ATP and Mg2+.ResultsThe 3D model of Diplonema was compared with available crystal structures from Trypanosoma brucei, bacteriophage T4, and two archaeans. Interaction of DpRNL with ATP is predicted to involve double π-stacking, which has not been reported before in RNA ligases. This particular contact would shift the orientation of ATP and have considerable consequences on the interaction network of amino acids in the catalytic pocket. We postulate that certain canonical amino acids assume different functional roles in DpRNL compared to structurally homologous residues in other RNA ligases 2, a reassignment indicative of constructive neutral evolution. Finally, both structure comparison and phylogenetic analysis show that DpRNL is not specifically related to RNA ligases from trypanosomes, suggesting a unique adaptation of the latter for RNA editing, after the split of diplonemids and kinetoplastids.ConclusionHomology modeling and molecular dynamics simulations strongly suggest that DpRNL is an RNA ligase 2. The predicted innovative reshaping of DpRNL’s catalytic pocket is worthwhile to be tested experimentally.Electronic supplementary materialThe online version of this article (doi:10.1186/s12900-015-0046-0) contains supplementary material, which is available to authorized users.

BackgroundHelicobacter pylori is an important factor in the development of diseases such as ulcer and gastric cancer. This bacterium uses a periplasmic transporter, UreI, to deliver urea to the intracelullar space, where later it is transformed into ammonia by the cytoplasmic enzyme urease to survive the acidic condition of the human stomach. The UreI transporter presents a pH-dependent activity, where this pH-dependence remains unknown at a structural level. Althought the existance of several protonable residues in the periplasmic loops are related to the pH-dependent activity, we find interesting to have a clear view of the conformational changes involved in this phenomena through a molecular dynamic study.ResultsMolecular dynamic simulations of the UreI transporter at three different pH conditions were performed, revealing two main pH-dependent conformations, which we present as the open and close states. We find that salt bridges between the periplasmic loops are crucial interactions that stabilize these conformations. Besides, a cooperative behaviour exists between the six subunits of the system that is necessary to fulfill the activity of this transporter.ConclusionsWe found different pH-dependent conformations of the urea transporter UreI from Helicobacter pylori, which are related to salt-bridge interactions in the periplasmic regions. The behaviour of every channel in the system is not independent, given the existance of a cooperative behaviour through the formation of salt-bridges between the subunits of the hexameric system. We believe that our results will be related to the generation of new eradication therapies using this transporter as an attractive target, denoting that the knowledge of the possible pH-dependent conformations adopted for this transporter are important for the development of rational drug design approximations.Electronic supplementary materialThe online version of this article (doi:10.1186/s12900-015-0038-0) contains supplementary material, which is available to authorized users.

BackgroundInteractions between the orexin peptides and their cognate OX1 and OX2 receptors remain poorly characterized. Site-directed mutagenesis studies on orexin peptides and receptors have indicated amino acids important for ligand binding and receptor activation. However, a better understanding of specific pairwise interactions would benefit small molecule discovery.ResultsWe constructed a set of three-dimensional models of the orexin 1 receptor based on the 3D-structures of the orexin 2 receptor (released while this manuscript was under review), neurotensin receptor 1 and chemokine receptor CXCR4, conducted an exhaustive docking of orexin-A16–33 peptide fragment with ZDOCK and RDOCK, and analyzed a total of 4301 complexes through multidimensional scaling and clustering. The best docking poses reveal two alternative binding modes, where the C-terminus of the peptide lies deep in the binding pocket, on average about 5–6 Å above Tyr6.48 and close to Gln3.32. The binding modes differ in the about 100° rotation of the peptide; the peptide His26 faces either the receptor’s fifth transmembrane helix or the seventh helix. Both binding modes are well in line with previous mutation studies and partake in hydrogen bonding similar to suvorexant.ConclusionsWe present two binding modes for orexin-A into orexin 1 receptor, which help rationalize previous results from site-directed mutagenesis studies. The binding modes should serve small molecule discovery, and offer insights into the mechanism of receptor activation.Electronic supplementary materialThe online version of this article (doi:10.1186/s12900-015-0036-2) contains supplementary material, which is available to authorized users.

BackgroundDeinococcus radiodurans is an extremely radiation and desiccation resistant bacterium which can tolerate radiation doses up to 5,000 Grays without losing viability. We are studying the role of DNA repair and replication proteins for this unusual phenotype by a structural biology approach. The DNA polymerase III β subunit (β-clamp) acts as a sliding clamp on DNA, promoting the binding and processivity of many DNA-acting proteins, and here we report the crystal structure of D. radiodurans β-clamp (Drβ-clamp) at 2.0 Å resolution.ResultsThe sequence verification process revealed that at the time of the study the gene encoding Drβ-clamp was wrongly annotated in the genome database, encoding a protein of 393 instead of 362 amino acids. The short protein was successfully expressed, purified and used for crystallisation purposes in complex with Cy5-labeled DNA. The structure, which was obtained from blue crystals, shows a typical ring-shaped bacterial β-clamp formed of two monomers, each with three domains of identical topology, but with no visible DNA in electron density. A visualisation of the electrostatic surface potential reveals a highly negatively charged outer surface while the inner surface and the dimer forming interface have a more even charge distribution.ConclusionsThe structure of Drβ-clamp was determined to 2.0 Å resolution and shows an evenly distributed electrostatic surface charge on the DNA interacting side. We hypothesise that this charge distribution may facilitate efficient movement on encircled DNA and help ensure efficient DNA metabolism in D. radiodurans upon exposure to high doses of ionizing irradiation or desiccation.