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.

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

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.

The beta hairpin structure within ribosomal protein S5 mediates interplay between domains II and IV and regulates HCV IRES function

Translation initiation in Hepatitis C Virus (HCV) is mediated by Internal Ribosome Entry Site (IRES), which is independent of cap-structure and uses a limited number of canonical initiation factors. During translation initiation IRES-40S complex formation depends on high affinity interaction of IRES with ribosomal proteins. Earlier, it has been shown that ribosomal protein S5 (RPS5) interacts with HCV IRES. Here, we have extensively characterized the HCV IRES-RPS5 interaction and demonstrated its role in IRES function. Computational modelling and RNA-protein interaction studies demonstrated that the beta hairpin structure within RPS5 is critically required for the binding with domains II and IV. Mutations disrupting IRES-RPS5 interaction drastically reduced the 80S complex formation and the corresponding IRES activity. Computational analysis and UV cross-linking experiments using various IRES-mutants revealed interplay between domains II and IV mediated by RPS5. In addition, present study demonstrated that RPS5 interaction is unique to HCV IRES and is not involved in 40S-3 ` UTR interaction. Further, partial silencing of RPS5 resulted in preferential inhibition of HCV RNA translation. However, global translation was marginally affected by partial silencing of RPS5. Taken together, results provide novel molecular insights into IRES-RPS5 interaction and unravel its functional significance in mediating internal initiation of translation.

Calcium binding properties of calcium dependent protein kinase 1 (CaCDPK1) from Cicer arietinum

Calcium plays a crucial role as a secondary messenger in all aspects of plant growth, development and survival. Calcium dependent protein kinases (CDPKs) are the major calcium decoders, which couple the changes in calcium level to an appropriate physiological response. The mechanism by which calcium regulates CDPK protein is not well understood. In this study, we investigated the interactions of Ca2+ ions with the CDPK1 isoform of Cicer arietinum (CaCDPK1) using a combination of biophysical tools. CaCDPK1 has four different EF hands as predicted by protein sequence analysis. The fluorescence emission spectrum of CaCDPK1 showed quenching with a 5 nm red shift upon addition of calcium, indicating conformational changes in the tertiary structure. The plot of changes in intensity against calcium concentrations showed a biphasic curve with binding constants of 1.29 mu M and 120 mu M indicating two kinds of binding sites. Isothermal calorimetric (ITC) titration with CaCl2 also showed a biphasic curve with two binding constants of 0.027 mu M and 1.7 mu M. Circular dichroism (CD) spectra showed two prominent peaks at 208 and 222 nm indicating that CaCDPK1 is a alpha-helical rich protein. Calcium binding further increased the alpha-helical content of CaCDPK1 from 75 to 81%. Addition of calcium to CaCDPK1 also increased fluorescence of 8-anilinonaphthalene-1-sulfonic acid (ANS) indicating exposure of hydrophobic surfaces. Thus, on the whole this study provides evidence for calcium induced conformational changes, exposure of hydrophobic surfaces and heterogeneity of EF hands in CaCDPK1. (C) 2015 Elsevier GmbH. All rights reserved.

Key challenges for the creation and maintenance of specialist protein resources

As the volume of data relating to proteins increases, researchers rely more and more on the analysis of published data, thus increasing the importance of good access to these data that vary from the supplemental material of individual articles, all the way to major reference databases with professional staff and long-term funding. Specialist protein resources fill an important middle ground, providing interactive web interfaces to their databases for a focused topic or family of proteins, using specialized approaches that are not feasible in the major reference databases. Many are labors of love, run by a single lab with little or no dedicated funding and there are many challenges to building and maintaining them. This perspective arose from a meeting of several specialist protein resources and major reference databases held at the Wellcome Trust Genome Campus (Cambridge, UK) on August 11 and 12, 2014. During this meeting some common key challenges involved in creating and maintaining such resources were discussed, along with various approaches to address them. In laying out these challenges, we aim to inform users about how these issues impact our resources and illustrate ways in which our working together could enhance their accuracy, currency, and overall value. Proteins 2015; 83:1005-1013. (c) 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

N-Terminal Disordered Domain of Saccharomyces cerevisiae Hop1 Protein Is Dispensable for DNA Binding, Bridging, and Synapsis of Double-Stranded DNA Molecules But Is Necessary for Spore Formation

The cytological architecture of the synaptonemal complex (SC), a meiosis-specific proteinaceous structure, is evolutionarily conserved among eukaryotes. However, little is known about the biochemical properties of SC components or the mechanisms underlying their roles in meiotic chromosome synapsis and recombination. Functional analysis of Saccharomyces cerevisiae Hop1, a key structural component of SC, has begun to reveal important insights into its function in interhomolog recombination. Previously, we showed that Hop1 is a structure-specific DNA-binding protein, exhibits higher binding affinity for the Holliday junction, and induces structural distortion at the core of the junction. Furthermore, Hop1 promotes DNA condensation and intra- and intermolecular synapsis between duplex DNA molecules. Here, we show that Hop1 possesses a modular domain organization, consisting of an intrinsically disordered N-terminal domain and a protease-resistant C-terminal domain (Hop1CTD). Furthermore, we found that Hop1CTD exhibits strong homotypic as well as heterotypic protein protein interactions, and its biochemical activities were similar to those of the full-length Hop1 protein. However, Hop1CTD failed to complement the meiotic recombination defects of the Delta hop1 strain, indicating that both N- and C-terminal domains of Hop1 are essential for meiosis and spore formation. Altogether, our findings reveal novel insights into the structure-function relationships of Hop1 and help to further our understanding of its role in meiotic chromosome synapsis and recombination.

A Multiantigenic DNA Vaccine That Induces Broad Hepatitis C Virus-Specific T-Cell Responses in Mice

There are 3 to 4 million new hepatitis C virus (HCV) infections annually around the world, but no vaccine is available. Robust T-cell mediated responses are necessary for effective clearance of the virus, and DNA vaccines result in a cell-mediated bias. Adjuvants are often required for effective vaccination, but during natural lytic viral infections damage-associated molecular patterns (DAMPs) are released, which act as natural adjuvants. Hence, a vaccine that induces cell necrosis and releases DAMPs will result in cell-mediated immunity (CMI), similar to that resulting from natural lytic viral infection. We have generated a DNA vaccine with the ability to elicit strong CMI against the HCV nonstructural (NS) proteins (3, 4A, 4B, and 5B) by encoding a cytolytic protein, perforin (PRF), and the antigens on a single plasmid. We examined the efficacy of the vaccines in C57BL/6 mice, as determined by gamma interferon enzyme-linked immunosorbent spot assay, cell proliferation studies, and intracellular cytokine production. Initially, we showed that encoding the NS4A protein in a vaccine which encoded only NS3 reduced the immunogenicity of NS3, whereas including PRF increased NS3 immunogenicity. In contrast, the inclusion of NS4A increased the immunogenicity of the NS3, NS4B, andNS5B proteins, when encoded in a DNA vaccine that also encoded PRF. Finally, vaccines that also encoded PRF elicited similar levels of CMI against each protein after vaccination with DNA encoding NS3, NS4A, NS4B, and NS5B compared to mice vaccinated with DNA encoding only NS3 or NS4B/5B. Thus, we have developed a promising ``multiantigen'' vaccine that elicits robust CMI. IMPORTANCE Since their development, vaccines have reduced the global burden of disease. One strategy for vaccine development is to use commercially viable DNA technology, which has the potential to generate robust immune responses. Hepatitis C virus causes chronic liver infection and is a leading cause of liver cancer. To date, no vaccine is currently available, and treatment is costly and often results in side effects, limiting the number of patients who are treated. Despite recent advances in treatment, prevention remains the key to efficient control and elimination of this virus. Here, we describe a novel DNA vaccine against hepatitis C virus that is capable of inducing robust cell-mediated immune responses in mice and is a promising vaccine candidate for humans.

Insights into the Functional Roles of N-Terminal and C-Terminal Domains of Helicobacter pylori DprA

DNA processing protein A (DprA) plays a crucial role in the process of natural transformation. This is accomplished through binding and subsequent protection of incoming foreign DNA during the process of internalization. DprA along with Single stranded DNA binding protein A (SsbA) acts as an accessory factor for RecA mediated DNA strand exchange. H. pylori DprA (HpDprA) is divided into an N-terminal domain and a C-terminal domain. In the present study, individual domains of HpDprA have been characterized for their ability to bind single stranded (ssDNA) and double stranded DNA (dsDNA). Oligomeric studies revealed that HpDprA possesses two sites for dimerization which enables HpDprA to form large and tightly packed complexes with ss and dsDNA. While the N-terminal domain was found to be sufficient for binding with ss or ds DNA, C-terminal domain has an important role in the assembly of poly-nucleoprotein complex. Using site directed mutagenesis approach, we show that a pocket comprising positively charged amino acids in the N-terminal domain has an important role in the binding of ss and dsDNA. Together, a functional cross talk between the two domains of HpDprA facilitating the binding and formation of higher order complex with DNA is discussed.

A carboxy terminal domain of the L protein of rinderpest virus possesses RNA triphosphatase activity - The first enzyme in the viral mRNA capping pathway

The large protein L of negative-sense RNA viruses is a multifunctional protein involved in transcription and replication of genomic RNA. It also possesses enzymatic activities involved in capping and methylation of viral mRNAs. The pathway for mRNA capping followed by the L protein of the viruses in the Morbillivirus genus has not been established, although it has been speculated that these viruses may follow the unconventional capping pathway as has been shown for some viruses of Rhabdoviridae family. We had earlier shown that the large protein L of Rinderpest virus expressed as recombinant L-P complex in insect cells as well as the ribonucleoprotein complex from purified virus possesses RNA triphosphatase (RTPase) and guanylyltransferase activities, in addition to RNA dependent RNA polymerase activity. In the present work, we demonstrate that RTPase as well as nucleoside triphosphatase (NTPase) activities are exhibited by a subdomain of the L protein in the C terminal region (a.a. 1640 1840). The RTPase activity depends absolutely on a divalent cation, either magnesium or manganese. Both the RTPase and NTPase activities of the protein show dual metal specificity. Two mutant proteins having alanine mutations in the glutamic acid residues in motif-A of the RTPase domain did not show RTPase activity, while exhibiting reduced NTPase activity suggesting overlapping active sites for the two enzymatic functions. The RTPase and NTPase activities of the L subdomain resemble those of the Vaccinia capping enzyme D1 and the baculovirus LEF4 proteins. (C) 2015 Elsevier Inc. All rights reserved.

Insights into stabilizing interactions in the distorted domain-swapped dimer of Salmonella typhimurium survival protein

The survival protein SurE from Salmonella typhimurium (StSurE) is a dimeric protein that functions as a phosphatase. SurE dimers are formed by the swapping of a loop with a pair of beta-strands and a C-terminal helix between two protomers. In a previous study, the Asp230 and His234 residues were mutated to Ala to abolish a hydrogen bond that was thought to be crucial for C-terminal helix swapping. These mutations led to functionally inactive and distorted dimers in which the two protomers were related by a rotation of 167 degrees. New salt bridges involving Glu112 were observed in the dimeric interface of the H234A and D230A/H234A mutants. To explore the role of these salt bridges in the stability of the distorted structure, E112A, E112A/D230A, E112A/H234A, E112A/D230A/H234A, R179L/H180A/H234A and E112A/R179L/H180A/H234A mutants were constructed. X-ray crystal structures of the E112A, E112A/H234A and E112A/D230A mutants could be determined. The dimeric structures of the E112A and E112A/H234A mutants were similar to that of native SurE, while the E112A/D230A mutant had a residual rotation of 11 degrees between the B chains upon superposition of the A chains of the mutant and native dimers. The native dimeric structure was nearly restored in the E112A/H234A mutant, suggesting that the new salt bridge observed in the H234A and D230A/H234A mutants was indeed responsible for the stability of their distorted structures. Catalytic activity was also restored in these mutants, implying that appropriate dimeric organization is necessary for the activity of SurE.