De novo inference of protein function from coarse-grained dynamics

Inference of molecular function of proteins is the fundamental task in the quest for understanding cellular processes. The task is getting increasingly difficult with thousands of new proteins discovered each day. The difficulty arises primarily due to lack of high-throughput experimental technique for assessing protein molecular function, a lacunae that computational approaches are trying hard to fill. The latter too faces a major bottleneck in absence of clear evidence based on evolutionary information. Here we propose a de novo approach to annotate protein molecular function through structural dynamics match for a pair of segments from two dissimilar proteins, which may share even <10% sequence identity. To screen these matches, corresponding 1 mu s coarse-grained (CG) molecular dynamics trajectories were used to compute normalized root-mean-square-fluctuation graphs and select mobile segments, which were, thereafter, matched for all pairs using unweighted three-dimensional autocorrelation vectors. Our in-house custom-built forcefield (FF), extensively validated against dynamics information obtained from experimental nuclear magnetic resonance data, was specifically used to generate the CG dynamics trajectories. The test for correspondence of dynamics-signature of protein segments and function revealed 87% true positive rate and 93.5% true negative rate, on a dataset of 60 experimentally validated proteins, including moonlighting proteins and those with novel functional motifs. A random test against 315 unique fold/function proteins for a negative test gave >99% true recall. A blind prediction on a novel protein appears consistent with additional evidences retrieved therein. This is the first proof-of-principle of generalized use of structural dynamics for inferring protein molecular function leveraging our custom-made CG FF, useful to all. (C) 2014 Wiley Periodicals, Inc.

Direct regulation of topoisomerase activity by a nucleoid-associated protein

The topological homeostasis of bacterial chromosomes is maintained by the balance between compaction and the topological organization of genomes. Two classes of proteins play major roles in chromosome organization: the nucleoid-associated proteins (NAPs) and topoisomerases. The NAPs bind DNA to compact the chromosome, whereas topoisomerases catalytically remove or introduce supercoils into the genome. We demonstrate that HU, a major NAP of Mycobacterium tuberculosis specifically stimulates the DNA relaxation ability of mycobacterial topoisomerase I (TopoI) at lower concentrations but interferes at higher concentrations. A direct physical interaction between M. tuberculosis HU (MtHU) and TopoI is necessary for enhancing enzyme activity both in vitro and in vivo. The interaction is between the amino terminal domain of MtHU and the carboxyl terminal domain of TopoI. Binding of MtHU did not affect the two catalytic trans-esterification steps but enhanced the DNA strand passage, requisite for the completion of DNA relaxation, a new mechanism for the regulation of topoisomerase activity. An interaction-deficient mutant of MtHU was compromised in enhancing the strand passage activity. The species-specific physical and functional cooperation between MtHU and TopoI may be the key to achieve the DNA relaxation levels needed to maintain the optimal superhelical density of mycobacterial genomes.

Charge density distribution and electrostatic interactions of ethionamide: an inhibitor of the enoyl acyl carrier protein reductase (inhA) enzyme of Mycobacterium tuberculosis

An experimental charge density analysis of an anti-TB drug ethionamide was carried out from high resolution X-ray diffraction at 100 K to understand its charge density distribution and electrostatic properties. The experimental results were validated from periodic theoretical charge density calculations performed using CRYSTAL09 at the B3LYP/6-31G** level of theory. The electron density rho(bcp)(r) and the Laplacian of electron density del(2)(rho bcp)(r) of the molecule calculated from both the methods display the charge density distribution of the ethionamide molecule in the crystal field. The electrostatic potential map shows a large electropositive region around the pyridine ring and a large electronegative region at the vicinity of the thiol atom. The calculated experimental dipole moment is 10.6D, which is higher than the value calculated from theory (8.2D). The topological properties of C-H center dot center dot center dot S, N-H center dot center dot center dot N and N-H center dot center dot center dot S hydrogen bonds were calculated, revealing their strength. The charge density analysis of the ethionamide molecule determined from both the experiment and theory gives the topological and electrostatic properties of the molecule, which allows to precisely understand the nature of intra and intermolecular interactions.

Sensitivity of polarization fluctuations to the nature of protein-water interactions: Study of biological water in four different protein-water systems

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (similar to 80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of similar to 2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water < dM(W) (0)delta M-W (t) > is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (similar to 50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work. (C) 2014 AIP Publishing LLC.

An overview of recent advances in structural bioinformatics of protein-protein interactions and a guide to their principles

Rich data bearing on the structural and evolutionary principles of protein protein interactions are paving the way to a better understanding of the regulation of function in the cell. This is particularly the case when these interactions are considered in the framework of key pathways. Knowledge of the interactions may provide insights into the mechanisms of crucial `driver' mutations in oncogenesis. They also provide the foundation toward the design of protein protein interfaces and inhibitors that can abrogate their formation or enhance them. The main features to learn from known 3-D structures of protein protein complexes and the extensive literature which analyzes them computationally and experimentally include the interaction details which permit undertaking structure-based drug discovery, the evolution of complexes and their interactions, the consequences of alterations such as post-translational modifications, ligand binding, disease causing mutations, host pathogen interactions, oligomerization, aggregation and the roles of disorder, dynamics, allostery and more to the protein and the cell. This review highlights some of the recent advances in these areas, including design, inhibition and prediction of protein protein complexes. The field is broad, and much work has been carried out in these areas, making it challenging to cover it in its entirety. Much of this is due to the fast increase in the number of molecules whose structures have been determined experimentally and the vast increase in computational power. Here we provide a concise overview. (C) 2014 Elsevier Ltd. All rights reserved.

Plasmodium berghei glycine cleavage system T-protein is non-essential for parasite survival in vertebrate and invertebrate hosts

T-protein, an aminomethyltransferase, represents one of the four components of glycine cleavage system (GCS) and catalyzes the transfer of methylene group from H-protein intermediate to tetrahydrofolate (THF) forming N-5, N-10-methylene THF (CH2-THF) with the release of ammonia. The malaria parasite genome encodes T-, H- and L-proteins, but not P-protein which is a glycine decarboxylase generating the aminomethylene group. A putative GCS has been considered to be functional in the parasite mitochondrion despite the absence of a detectable P-protein homologue. In the present study, the mitochondrial localization of T-protein in the malaria parasite was confirmed by immunofluorescence and its essentiality in the entire parasite life cycle was studied by targeting the T-protein locus in Plasmodium berghei (Pb). PbT knock out parasites did not show any growth defect in asexual, sexual and liver stages indicating that the T-protein is dispensable for parasite survival in vertebrate and invertebrate hosts. The absence of P-protein homologue and the non-essentiality of T protein suggest the possible redundancy of GCS activity in the malaria parasite. Nevertheless, the H- and L-proteins of GCS could be essential for malaria parasite because of their involvement in alpha-lcetoacid dehydrogenase reactions. (C) 2014 Elsevier B.V. All rights reserved.

From Workstations to Workbenches: Towards Predicting Physicochemically Viable Protein-Protein Interactions Across a Host and a Pathogen

The understanding of protein-protein interactions is indispensable in comprehending most of the biological processes in a cell. Small-scale experiments as well as large-scale high-throughput techniques over the past few decades have facilitated identification and analysis of protein-protein interactions which form the basis of much of our knowledge on functional and regulatory aspects of proteins. However, such rich catalog of interaction data should be used with caution when establishing protein-protein interactions in silico, as the high-throughput datasets are prone to false positives. Numerous computational means developed to pursue genome-wide studies on protein-protein interactions at times overlook the mechanistic and molecular details, thus questioning the reliability of predicted protein-protein interactions. We review the development, advantages, and shortcomings of varied approaches and demonstrate that by providing a structural viewpoint in terms of shape complementarity and interaction energies at protein-protein interfaces coupled with information on expression and localization of proteins homologous to an interacting pair, it is possible to assess the credibility of predicted interactions in biological context. With a focus on human pathogen Mycobacterium tuberculosis H37Rv, we show that such scrupulous use of details at the molecular level can predict physicochemically viable protein-protein interactions across host and pathogen. Such predicted interactions have the potential to provide molecular basis of probable mechanisms of pathogenesis and hence open up ways to explore their usefulness as targets in the light of drug discovery. (c) 2014 IUBMB Life, 66(11):759-774, 2014

Mycobacterium tuberculosis Cell Division Protein, FtsE, is an ATPase in Dimeric Form

FtsE is one of the earliest cell division proteins that assembles along with FtsX at the mid-cell site during cell division in Escherichia coli. Both these proteins are highly conserved across diverse bacterial genera and are predicted to constitute an ABC transporter type complex, in which FtsE is predicted to bind ATP and hydrolyse it, and FtsX is predicted to be an integral membrane protein. We had earlier reported that the MtFtsE of the human pathogen, Mycobacterium tuberculosis, binds ATP and interacts with MtFtsX on the cell membrane of M. tuberculosis and E. coli. In this study, we demonstrate that MtFtsE is an ATPase, the active form of which is a dimer, wherein the participating monomers are held together by non-covalent interactions, with the Cys84 of each monomer present at the dimer interface. Under oxidising environment, the dimer gets stabilised by the formation of Cys84-Cys84 disulphide bond. While the recombinant MtFtsE forms a dimer on the membrane of E. coli, the native MtFtsE seems to be in a different conformation in the M. tuberculosis membrane. Although disulphide bridges were not observed on the cytoplasmic side (reducing environment) of the membrane, the two participating monomers could be isolated as dimers held together by non-covalent interactions. Taken together, these findings show that MtFtsE is an ATPase in the non-covalent dimer form, with the Cys84 of each monomer present in the reduced form at the dimer interface, without participating in the dimerisation or the catalytic activity of the protein.

Chaperone protein HYPK interacts with the first 17 amino acid region of Huntingtin and modulates mutant HTT-mediated aggregation and cytotoxicity

Huntington's disease is a polyglutamine expansion disorder, characterized by mutant HTT-mediated aggregate formation and cytotoxicity. Many reports suggests roles of N-terminal 17 amino acid domain of HTT (HTT-N17) towards subcellular localization, aggregate formation and subsequent pathogenicity induced by N-terminal HTT harboring polyQ stretch in pathogenic range. HYPK is a HTT-interacting chaperone which can reduce N-terminal mutant HTT-mediated aggregate formation and cytotoxicity in neuronal cell lines. However, how HYPK interacts with N-terminal fragment of HTT remained unknown. Here we report that specific interaction of HYPK with HTT-N17 is crucial for the chaperone activity of HYPK. Deletion of HTT-N17 leads to formation of tinier, SDS-soluble nuclear aggregates formed by N-terminal mutant HTT. The increased cytotoxicity imparted by these tiny aggregates might be contributed due to loss of interaction with HYPK. (C) 2014 Elsevier Inc. All rights reserved.

Sensitive detection of C-reactive protein using optical fiber Bragg gratings

An accurate and highly sensitive sensor platform has been demonstrated for the detection of C-reactive protein (CRP) using optical fiber Bragg gratings (FBGs). The CRP detection has been carried out by monitoring the shift in Bragg wavelength (Delta lambda(B)) of an etched FBG (eFBG) coated with an anti-CRP antibody (aCRP)-graphene oxide (GO) complex. The complex is characterized by Fourier transform infrared spectroscopy, X-ray photoelectron spectroscopy and atomic force microscopy. A limit of detection of 0.01 mg/L has been achieved with a linear range of detection from 0.01 mg/L to 100 mg/L which includes clinical range of CRP. The eFBG sensor coated with only aCRP (without GO) show much less sensitivity than that of aCRP-GO complex coated eFBG. The eFBG sensors show high specificity to CRP even in the presence of other interfering factors such as urea, creatinine and glucose. The affinity constant of similar to 1.1 x 10(10) M-1 has been extracted from the data of normalized shift (Delta lambda(B)/lambda(B)) as a function of CRP concentration. (C) 2014 Elsevier B.V. All rights reserved.

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

NrichD ( ext-link-type=''uri'' xlink:href=''http://proline.biochem.iisc.ernet.in/NRICHD/'' xlink:type=''simple''>http://proline.biochem.iisc.ernet.in/NRICHD/)< /named-content> 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.

Demonstration of helicase activity in the nonstructural protein, NSs, of the negative-sense RNA virus, Groundnut bud necrosis virus

The nonstructural protein NSs, encoded by the S RNA of groundnut bud necrosis virus (GBNV) (genus Tospovirus, family Bunyaviridae) has earlier been shown to possess nucleic-acid-stimulated NTPase and 50 a phosphatase activity. ATP hydrolysis is an essential function of a true helicase. Therefore, NSs was tested for DNA helicase activity. The results demonstrated that GBNV NSs possesses bidirectional DNA helicase activity. An alanine mutation in the Walker A motif (K189A rNSs) decreased DNA helicase activity substantially, whereas a mutation in the Walker B motif resulted in a marginal decrease in this activity. The parallel loss of the helicase and ATPase activity in the K189A mutant confirms that NSs acts as a non-canonical DNA helicase. Furthermore, both the wild-type and K189A NSs could function as RNA silencing suppressors, demonstrating that the suppressor activity of NSs is independent of its helicase or ATPase activity. This is the first report of a true helicase from a negative-sense RNA virus.

A Gaussian network model study suggests that structural fluctuations are higher for inactive states than active states of protein kinases

We performed Gaussian network model based normal mode analysis of 3-dimensional structures of multiple active and inactive forms of protein kinases. In 14 different kinases, a more number of residues (1095) show higher structural fluctuations in inactive states than those in active states (525), suggesting that, in general, mobility of inactive states is higher than active states. This statistically significant difference is consistent with higher crystallographic B-factors and conformational energies for inactive than active states, suggesting lower stability of inactive forms. Only a small number of inactive conformations with the DFG motif in the ``in'' state were found to have fluctuation magnitudes comparable to the active conformation. Therefore our study reports for the first time, intrinsic higher structural fluctuation for almost all inactive conformations compared to the active forms. Regions with higher fluctuations in the inactive states are often localized to the aC-helix, aG-helix and activation loop which are involved in the regulation and/or in structural transitions between active and inactive states. Further analysis of 476 kinase structures involved in interactions with another domain/protein showed that many of the regions with higher inactive-state fluctuation correspond to contact interfaces. We also performed extensive GNM analysis of (i) insulin receptor kinase bound to another protein and (ii) holo and apo forms of active and inactive conformations followed by multi-factor analysis of variance. We conclude that binding of small molecules or other domains/proteins reduce the extent of fluctuation irrespective of active or inactive forms. Finally, we show that the perceived fluctuations serve as a useful input to predict the functional state of a kinase.

Typical and atypical domain combinations in human protein kinases: functions, disease causing mutations and conservation in other primates

Ser/Thr and Tyr protein kinases orchestrate many signalling pathways and hence loss in this balance leads to many disease phenotypes. Due to their high abundance, diversity and importance, efforts have been made in the past to classify kinases and annotate their functions at both gross and fine levels. These kinases are conventionally classified into subfamilies based on the sequences of catalytic domains. Usually the domain architecture of a full-length kinase is consistent with the subfamily classification made based on the sequence of kinase domain. Important contributions of modular domains to the overall function of the kinase are well known. Recently occurrence of two kinds of outlier kinases-''Hybrid'' and ``Rogue'' has been reported. These show considerable deviations in their domain architectures from the typical domain architecture known for the classical kinase subfamilies. This article provides an overview of the different subfamilies of human kinases and the role of non-kinase domains in functions and diseases. Importantly this article provides analysis of hybrid and rogue kinases encoded in the human genome and highlights their conservation in closely related primate species. These kinases are examples of elegant rewiring to bring about subtle functional differences compared to canonical variants.