Insights into the Fold Organization of TIM Barrel from Interaction Energy Based Structure Networks

There are many well-known examples of proteins with low sequence similarity, adopting the same structural fold. This aspect of sequence-structure relationship has been extensively studied both experimentally and theoretically, however with limited success. Most of the studies consider remote homology or ``sequence conservation'' as the basis for their understanding. Recently ``interaction energy'' based network formalism (Protein Energy Networks (PENs)) was developed to understand the determinants of protein structures. In this paper we have used these PENs to investigate the common non-covalent interactions and their collective features which stabilize the TIM barrel fold. We have also developed a method of aligning PENs in order to understand the spatial conservation of interactions in the fold. We have identified key common interactions responsible for the conservation of the TIM fold, despite high sequence dissimilarity. For instance, the central beta barrel of the TIM fold is stabilized by long-range high energy electrostatic interactions and low-energy contiguous vdW interactions in certain families. The other interfaces like the helix-sheet or the helix-helix seem to be devoid of any high energy conserved interactions. Conserved interactions in the loop regions around the catalytic site of the TIM fold have also been identified, pointing out their significance in both structural and functional evolution. Based on these investigations, we have developed a novel network based phylogenetic analysis for remote homologues, which can perform better than sequence based phylogeny. Such an analysis is more meaningful from both structural and functional evolutionary perspective. We believe that the information obtained through the ``interaction conservation'' viewpoint and the subsequently developed method of structure network alignment, can shed new light in the fields of fold organization and de novo computational protein design.

Local Structural Differences in Homologous Proteins: Specificities in Different SCOP Classes

The constant increase in the number of solved protein structures is of great help in understanding the basic principles behind protein folding and evolution. 3-D structural knowledge is valuable in designing and developing methods for comparison, modelling and prediction of protein structures. These approaches for structure analysis can be directly implicated in studying protein function and for drug design. The backbone of a protein structure favours certain local conformations which include alpha-helices, beta-strands and turns. Libraries of limited number of local conformations (Structural Alphabets) were developed in the past to obtain a useful categorization of backbone conformation. Protein Block (PB) is one such Structural Alphabet that gave a reasonable structure approximation of 0.42 angstrom. In this study, we use PB description of local structures to analyse conformations that are preferred sites for structural variations and insertions, among group of related folds. This knowledge can be utilized in improving tools for structure comparison that work by analysing local structure similarities. Conformational differences between homologous proteins are known to occur often in the regions comprising turns and loops. Interestingly, these differences are found to have specific preferences depending upon the structural classes of proteins. Such class-specific preferences are mainly seen in the all-beta class with changes involving short helical conformations and hairpin turns. A test carried out on a benchmark dataset also indicates that the use of knowledge on the class specific variations can improve the performance of a PB based structure comparison approach. The preference for the indel sites also seem to be confined to a few backbone conformations involving beta-turns and helix C-caps. These are mainly associated with short loops joining the regular secondary structures that mediate a reversal in the chain direction. Rare beta-turns of type I' and II' are also identified as preferred sites for insertions.

Crowd Sourcing a New Paradigm for Interactome Driven Drug Target Identification in Mycobacterium tuberculosis

A decade since the availability of Mycobacterium tuberculosis (Mtb) genome sequence, no promising drug has seen the light of the day. This not only indicates the challenges in discovering new drugs but also suggests a gap in our current understanding of Mtb biology. We attempt to bridge this gap by carrying out extensive re-annotation and constructing a systems level protein interaction map of Mtb with an objective of finding novel drug target candidates. Towards this, we synergized crowd sourcing and social networking methods through an initiative `Connect to Decode' (C2D) to generate the first and largest manually curated interactome of Mtb termed `3interactome pathway' (IPW), encompassing a total of 1434 proteins connected through 2575 functional relationships. Interactions leading to gene regulation, signal transduction, metabolism, structural complex formation have been catalogued. In the process, we have functionally annotated 87% of the Mtb genome in context of gene products. We further combine IPW with STRING based network to report central proteins, which may be assessed as potential drug targets for development of drugs with least possible side effects. The fact that five of the 17 predicted drug targets are already experimentally validated either genetically or biochemically lends credence to our unique approach.

Methods for Detecting Early Warnings of Critical Transitions in Time Series Illustrated Using Simulated Ecological Data

Many dynamical systems, including lakes, organisms, ocean circulation patterns, or financial markets, are now thought to have tipping points where critical transitions to a contrasting state can happen. Because critical transitions can occur unexpectedly and are difficult to manage, there is a need for methods that can be used to identify when a critical transition is approaching. Recent theory shows that we can identify the proximity of a system to a critical transition using a variety of so-called `early warning signals', and successful empirical examples suggest a potential for practical applicability. However, while the range of proposed methods for predicting critical transitions is rapidly expanding, opinions on their practical use differ widely, and there is no comparative study that tests the limitations of the different methods to identify approaching critical transitions using time-series data. Here, we summarize a range of currently available early warning methods and apply them to two simulated time series that are typical of systems undergoing a critical transition. In addition to a methodological guide, our work offers a practical toolbox that may be used in a wide range of fields to help detect early warning signals of critical transitions in time series data.

Analysis of Conformational Variation in Macromolecular Structural Models

Experimental conditions or the presence of interacting components can lead to variations in the structural models of macromolecules. However, the role of these factors in conformational selection is often omitted by in silico methods to extract dynamic information from protein structural models. Structures of small peptides, considered building blocks for larger macromolecular structural models, can substantially differ in the context of a larger protein. This limitation is more evident in the case of modeling large multi-subunit macromolecular complexes using structures of the individual protein components. Here we report an analysis of variations in structural models of proteins with high sequence similarity. These models were analyzed for sequence features of the protein, the role of scaffolding segments including interacting proteins or affinity tags and the chemical components in the experimental conditions. Conformational features in these structural models could be rationalized by conformational selection events, perhaps induced by experimental conditions. This analysis was performed on a non-redundant dataset of protein structures from different SCOP classes. The sequence-conformation correlations that we note here suggest additional features that could be incorporated by in silico methods to extract dynamic information from protein structural models.

The Hinge Region of Human Thyroid-Stimulating Hormone (TSH) Receptor Operates as a Tunable Switch between Hormone Binding and Receptor Activation

The mechanism by which the hinge regions of glycoprotein hormone receptors couple hormone binding to activation of downstream effecters is not clearly understood. In the present study, agonistic (311.62) and antagonistic (311.87) monoclonal antibodies (MAbs) directed against the TSH receptor extracellular domain were used to elucidate role of the hinge region in receptor activation. MAb 311.62 which identifies the LRR/Cb-2 junction (aa 265-275), increased the affinity of TSHR for the hormone while concomitantly decreasing its efficacy, whereas MAb 311.87 recognizing LRR 7-9 (aa 201-259) acted as a non-competitive inhibitor of Thyroid stimulating hormone (TSH) binding. Binding of MAbs was sensitive to the conformational changes caused by the activating and inactivating mutations and exhibited differential effects on hormone binding and response of these mutants. By studying the effects of these MAbs on truncation and chimeric mutants of thyroid stimulating hormone receptor (TSHR), this study confirms the tethered inverse agonistic role played by the hinge region and maps the interactions between TSHR hinge region and exoloops responsible for maintenance of the receptor in its basal state. Mechanistic studies on the antibody-receptor interactions suggest that MAb 311.87 is an allosteric insurmountable antagonist and inhibits initiation of the hormone induced conformational changes in the hinge region, whereas MAb 311.62 acts as a partial agonist that recognizes a conformational epitope critical for coupling of hormone binding to receptor activation. The hinge region, probably in close proximity with the alpha-subunit in the hormone-receptor complex, acts as a tunable switch between hormone binding and receptor activation.

Comparative Transcriptomics of H. pylori Strains AM5, SS1 and Their hpyAVIBM Deletion Mutants: Possible Roles of Cytosine Methylation

Helicobacter pylori is an important human pathogen and one of the most successful chronic colonizers of the human body. H. pylori uses diverse mechanisms to modulate its interaction with the host in order to promote chronic infection and overcome host immune response. Restriction-modification genes are a major part of strain-specific genes present in H. pylori. The role of N-6 -adenine methylation in bacterial gene regulation and virulence is well established but not much is known about the effect of C-5 -cytosine methylation on gene expression in prokaryotes. In this study, it was observed by microarray analysis and RT-PCR, that deletion of an orphan C-5 -cytosine methyltransferase, hpyAVIBM in H. pylori strains AM5and SS1 has a significant effect on the expression of number of genes belonging to motility, adhesion and virulence. AM Delta DhpyAVIBM mutant strain has a different LPS profile and is able to induce high IL-8 production compared to wild-type. hpyAVIBM from strain 26695 is able to complement mutant SS1 and AM5 strains. This study highlights a possible significance of cytosine methylation in the physiology of H. pylori.

Coalescent Method in Conjunction with Niche Modeling Reveals Cryptic Diversity among Centipedes in the Western Ghats of South India

Background: There has been growing interest in integrative taxonomy that uses data from multiple disciplines for species delimitation. Typically, in such studies, monophyly is taken as a proxy for taxonomic distinctiveness and these units are treated as potential species. However, monophyly could arise due to stochastic processes. Thus here, we have employed a recently developed tool based on coalescent approach to ascertain the taxonomic distinctiveness of various monophyletic units. Subsequently, the species status of these taxonomic units was further tested using corroborative evidence from morphology and ecology. This inter-disciplinary approach was implemented on endemic centipedes of the genus Digitipes (Attems 1930) from the Western Ghats (WG) biodiversity hotspot of India. The species of the genus Digitipes are morphologically conserved, despite their ancient late Cretaceous origin. Principal Findings: Our coalescent analysis based on mitochondrial dataset indicated the presence of nine putative species. The integrative approach, which includes nuclear, morphology, and climate datasets supported distinctiveness of eight putative species, of which three represent described species and five were new species. Among the five new species, three were morphologically cryptic species, emphasizing the effectiveness of this approach in discovering cryptic diversity in less explored areas of the tropics like the WG. In addition, species pairs showed variable divergence along the molecular, morphological and climate axes. Conclusions: A multidisciplinary approach illustrated here is successful in discovering cryptic diversity with an indication that the current estimates of invertebrate species richness for the WG might have been underestimated. Additionally, the importance of measuring multiple secondary properties of species while defining species boundaries was highlighted given variable divergence of each species pair across the disciplines.

An Antitubulin Agent BCFMT Inhibits Proliferation of Cancer Cells and Induces Cell Death by Inhibiting Microtubule Dynamics

Using cell based screening assay, we identified a novel anti-tubulin agent (Z)-5-((5-(4-bromo-3-chlorophenyl)furan-2-yl)methylene)-2-thioxothiazoli din-4-one (BCFMT) that inhibited proliferation of human cervical carcinoma (HeLa) (IC50, 7.2 +/- 1.8 mu M), human breast adenocarcinoma (MCF-7) (IC50, 10.0 +/- 0.5 mu M), highly metastatic breast adenocarcinoma (MDA-MB-231) (IC50, 6.0 +/- 1 mu M), cisplatin-resistant human ovarian carcinoma (A2780-cis) (IC50, 5.8 +/- 0.3 mu M) and multi-drug resistant mouse mammary tumor (EMT6/AR1) (IC50, 6.5 +/- 1 mu M) cells. Using several complimentary strategies, BCFMT was found to inhibit cancer cell proliferation at G2/M phase of the cell cycle apparently by targeting microtubules. In addition, BCFMT strongly suppressed the dynamics of individual microtubules in live MCF-7 cells. At its half maximal proliferation inhibitory concentration (10 mu M), BCFMT reduced the rates of growing and shortening phases of microtubules in MCF-7 cells by 37 and 40%, respectively. Further, it increased the time microtubules spent in the pause (neither growing nor shortening detectably) state by 135% and reduced the dynamicity (dimer exchange per unit time) of microtubules by 70%. In vitro, BCFMT bound to tubulin with a dissociation constant of 8.3 +/- 1.8 mu M, inhibited tubulin assembly and suppressed GTPase activity of microtubules. BCFMT competitively inhibited the binding of BODIPY FL-vinblastine to tubulin with an inhibitory concentration (K-i) of 5.2 +/- 1.5 mu M suggesting that it binds to tubulin at the vinblastine site. In cultured cells, BCFMT-treatment depolymerized interphase microtubules, perturbed the spindle organization and accumulated checkpoint proteins (BubR1 and Mad2) at the kinetochores. BCFMT-treated MCF-7 cells showed enhanced nuclear accumulation of p53 and its downstream p21, which consequently activated apoptosis in these cells. The results suggested that BCFMT inhibits proliferation of several types of cancer cells including drug resistance cells by suppressing microtubule dynamics and indicated that the compound may have chemotherapeutic potential.

Novel Levamisole Derivative Induces Extrinsic Pathway of Apoptosis in Cancer Cells and Inhibits Tumor Progression in Mice

Background: Levamisole, an imidazo(2,1-b) thiazole derivative, has been reported to be a potential antitumor agent. In the present study, we have investigated the mechanism of action of one of the recently identified analogues, 4a (2-benzyl-6-(4'-fluorophenyl)-5-thiocyanato-imidazo2,1-b]1,3,4]thi adiazole). Materials and Methods: ROS production and expression of various apoptotic proteins were measured following 4a treatment in leukemia cell lines. Tumor animal models were used to evaluate the effect of 4a in comparison with Levamisole on progression of breast adenocarcinoma and survival. Immunohistochemistry and western blotting studies were performed to understand the mechanism of 4a action both ex vivo and in vivo. Results: We have determined the IC50 value of 4a in many leukemic and breast cancer cell lines and found CEM cells most sensitive (IC50 5 mu M). Results showed that 4a treatment leads to the accumulation of ROS. Western blot analysis showed upregulation of pro-apoptotic proteins t-BID and BAX, upon treatment with 4a. Besides, dose-dependent activation of p53 along with FAS, FAS-L, and cleavage of CASPASE-8 suggest that it induces death receptor mediated apoptotic pathway in CEM cells. More importantly, we observed a reduction in tumor growth and significant increase in survival upon oral administration of 4a (20 mg/kg, six doses) in mice. In comparison, 4a was found to be more potent than its parental analogue Levamisole based on both ex vivo and in vivo studies. Further, immunohistochemistry and western blotting studies indicate that 4a treatment led to abrogation of tumor cell proliferation and activation of apoptosis by the extrinsic pathway even in animal models. Conclusion: Thus, our results suggest that 4a could be used as a potent chemotherapeutic agent.

An Exported Heat Shock Protein 40 Associates with Pathogenesis-Related Knobs in Plasmodium falciparum Infected Erythrocytes

Cell surface structures termed knobs are one of the most important pathogenesis related protein complexes deployed by the malaria parasite Plasmodium falciparum at the surface of the infected erythrocyte. Despite their relevance to the disease, their structure, mechanisms of traffic and their process of assembly remain poorly understood. In this study, we have explored the possible role of a parasite-encoded Hsp40 class of chaperone, namely PFB0090c/PF3D7_0201800 (KAHsp40) in protein trafficking in the infected erythrocyte. We found the gene coding for PF3D7_0201800 to be located in a chromosomal cluster together with knob components KAHRP and PfEMP3. Like the knob components, KAHsp40 too showed the presence of PEXEL motif required for transport to the erythrocyte compartment. Indeed, sub-cellular fractionation and immunofluorescence analysis (IFA) showed KAHsp40 to be exported in the erythrocyte cytoplasm in a stage dependent manner localizing as punctuate spots in the erythrocyte periphery, distinctly from Maurer's cleft, in structures which could be the reminiscent of knobs. Double IFA analysis revealed co-localization of PF3D7_0201800 with the markers of knobs (KAHRP, PfEMP1 and PfEMP3) and components of the PEXEL translocon (Hsp101, PTEX150). KAHsp40 was also found to be in a complex with KAHRP, PfEMP3 and Hsp101 as confirmed by co-immunoprecipitation assay. Our results suggest potential involvement of a parasite encoded Hsp40 in chaperoning knob assembly in the erythrocyte compartment.

Distinct and Contrasting Transcription Initiation Patterns at Mycobacterium tuberculosis Promoters

Although sequencing of Mycobacterium tuberculosis genome lead to better understanding of transcription units and gene functions, interactions occurring during transcription initiation between RNA polymerase and promoters is yet to be elucidated. Different stages of transcription initiation include promoter specific binding of RNAP, isomerization, abortive initiation and promoter clearance. We have now analyzed these events with four promoters of M. tuberculosis viz. P-gyrB1, P-gyrR, P-rrnPCL1 and P-metU. The promoters differed from each other in their rates of open complex formation, decay, promoter clearance and abortive transcription. The equilibrium binding and kinetic studies of various steps revealed distinct rate limiting events for each of the promoter, which also differed markedly in their characteristics from the respective promoters of Mycobacterium smegmatis. Surprisingly, the transcription at gyr promoter was enhanced in the presence of initiating nucleotides and decreased in the presence of alarmone, pppGpp, a pattern typically seen with rRNA promoters studied so far. The gyr promoter of M. smegmatis, on the other hand, was not subjected to pppGpp mediated regulation. The marked differences in the transcription initiation pathway seen with rrn and gyr promoters of M. smegmatis and M. tuberculosis suggest that such species specific differences in the regulation of expression of the crucial housekeeping genes could be one of the key determinants contributing to the differences in growth rate and lifestyle of the two organisms. Moreover, the distinct rate limiting steps during transcription initiation of each one of the promoters studied point at variations in their intracellular regulation.

Crystallographic Study of Novel Transthyretin Ligands Exhibiting Negative-Cooperativity between Two Thyroxine Binding Sites

Background: Transthyretin (TTR) is a homotetrameric serum and cerebrospinal fluid protein that transports thyroxine (T4) and retinol by binding to retinol binding protein. Rate-limiting tetramer dissociation and rapid monomer misfolding and disassembly of TTR lead to amyloid fibril formation in different tissues causing various amyloid diseases. Based on the current understanding of the pathogenesis of TTR amyloidosis, it is considered that the inhibition of amyloid fibril formation by stabilization of TTR in native tetrameric form is a viable approach for the treatment of TTR amyloidosis. Methodology and Principal Findings: We have examined interactions of the wtTTR with a series of compounds containing various substitutions at biphenyl ether skeleton and a novel compound, previously evaluated for binding and inhibiting tetramer dissociation, by x-ray crystallographic approach. High resolution crystal structures of five ligands in complex with wtTTR provided snapshots of negatively cooperative binding of ligands in two T4 binding sites besides characterizing their binding orientations, conformations, and interactions with binding site residues. In all complexes, the ligand has better fit and more potent interactions in first T4 site i.e. (AC site) than the second T4 site (BD site). Together, these results suggest that AC site is a preferred ligand binding site and retention of ordered water molecules between the dimer interfaces further stabilizes the tetramer by bridging a hydrogen bond interaction between Ser117 and its symmetric copy. Conclusion: Novel biphenyl ether based compounds exhibit negative-cooperativity while binding to two T4 sites which suggests that binding of only single ligand molecule is sufficient to inhibit the TTR tetramer dissociation.

Differential Activities of the Two Closely Related Withanolides, Withaferin A and Withanone: Bioinformatics and Experimental Evidences

Background and Purpose: Withanolides are naturally occurring chemical compounds. They are secondary metabolites produced via oxidation of steroids and structurally consist of a steroid-backbone bound to a lactone or its derivatives. They are known to protect plants against herbivores and have medicinal value including anti-inflammation, anti-cancer, adaptogenic and anti-oxidant effects. Withaferin A (Wi-A) and Withanone (Wi-N) are two structurally similar withanolides isolated from Withania somnifera, also known as Ashwagandha in Indian Ayurvedic medicine. Ashwagandha alcoholic leaf extract (i-Extract), rich in Wi-N, was shown to kill cancer cells selectively. Furthermore, the two closely related purified phytochemicals, Wi-A and Wi-N, showed differential activity in normal and cancer human cells in vitro and in vivo. We had earlier identified several genes involved in cytotoxicity of i-Extract in human cancer cells by loss-of-function assays using either siRNA or randomized ribozyme library. Methodology/Principal Findings: In the present study, we have employed bioinformatics tools on four genes, i.e., mortalin, p53, p21 and Nrf2, identified by loss-of-function screenings. We examined the docking efficacy of Wi-N and Wi-A to each of the four targets and found that the two closely related phytochemicals have differential binding properties to the selected cellular targets that can potentially instigate differential molecular effects. We validated these findings by undertaking parallel experiments on specific gene responses to either Wi-N or Wi-A in human normal and cancer cells. We demonstrate that Wi-A that binds strongly to the selected targets acts as a strong cytotoxic agent both for normal and cancer cells. Wi-N, on the other hand, has a weak binding to the targets; it showed milder cytotoxicity towards cancer cells and was safe for normal cells. The present molecular docking analyses and experimental evidence revealed important insights to the use of Wi-A and Wi-N for cancer treatment and development of new anti-cancer phytochemical cocktails.

Crystal Structure of a Monomeric Thiolase-Like Protein Type 1 (TLP1) from Mycobacterium smegmatis

An analysis of the Mycobacterium smegmatis genome suggests that it codes for several thiolases and thiolase-like proteins. Thiolases are an important family of enzymes that are involved in fatty acid metabolism. They occur as either dimers or tetramers. Thiolases catalyze the Claisen condensation of two acetyl-Coenzyme A molecules in the synthetic direction and the thiolytic cleavage of 3-ketoacyl-Coenzyme A molecules in the degradative direction. Some of the M. smegmatis genes have been annotated as thiolases of the poorly characterized SCP2-thiolase subfamily. The mammalian SCP2-thiolase consists of an N-terminal thiolase domain followed by an additional C-terminal domain called sterol carrier protein-2 or SCP2. The M. smegmatis protein selected in the present study, referred to here as the thiolase-like protein type 1 (MsTLP1), has been biochemically and structurally characterized. Unlike classical thiolases, MsTLP1 is a monomer in solution. Its structure has been determined at 2.7 angstrom resolution by the single wavelength anomalous dispersion method. The structure of the protomer confirms that the N-terminal domain has the thiolase fold. An extra C-terminal domain is indeed observed. Interestingly, it consists of six beta-strands forming an anti-parallel beta-barrel which is completely different from the expected SCP2-fold. Detailed sequence and structural comparisons with thiolases show that the residues known to be essential for catalysis are not conserved in MsTLP1. Consistent with this observation, activity measurements show that MsTLP1 does not catalyze the thiolase reaction. This is the first structural report of a monomeric thiolase-like protein from any organism. These studies show that MsTLP1 belongs to a new group of thiolase related proteins of unknown function.