PROTEIN DISULFIDE ANALYSIS AND DESIGN

The significant contribution of naturally occurring disulfide bonds to protein stability has encouraged development of methods to engineer non-native disulfides in proteins. These have yielded mixed results. We summarize applications of the program MODIP for disulfide engineering. The program predicts sites in proteins where disulfides can be stably introduced. The program has also been used as an aid in conformational analysis of naturally occurring disulfides in a-helices, antiparallel and parallel beta-strands. Disulfides in a-helices occur only at N-termini, where the first cysteine residue is the N-cap residue of the helix. The disulfide occurs as a CXXC motif and can possess redox activity. In antiparallel beta-strands, disulfides occur exclusively at non-hydrogen bonded (NHB) registered pairs of antiparallel beta-sheets with only 1 known natural example occurring at a hydrogen bonded (HB) registered pair. Conformational analysis suggests that disulfides between HB residue pairs are under torsional strain. A similar analysis to characterize disulfides in parallel beta-strands was carried out. We observed that only 9 instances of cross-strand disulfides exist in a non-redundant dataset. Stereochemical analysis shows that while tbe chi(square) angles are similar to those of other disulfides, the chi(1) and chi(2) angles show more variation and that one of tbe strands is generally an edge strand.

Sparsely populated residue conformations in protein structures: Revisiting ``experimental'' Ramachandran maps

The Ramachandran map clearly delineates the regions of accessible conformational (phi-) space for amino acid residues in proteins. Experimental distributions of phi, values in high-resolution protein structures, reveal sparsely populated zones within fully allowed regions and distinct clusters in apparently disallowed regions. Conformational space has been divided into 14 distinct bins. Residues adopting these relatively rare conformations are presented and amino acid propensities for these regions are estimated. Inspection of specific examples in a completely arid, fully allowed region in the top left quadrant establishes that side-chain and backbone interactions may provide the energetic compensation necessary for populating this region of phi- space. Asn, Asp, and His residues showed the highest propensities in this region. The two distinct clusters in the bottom right quadrant which are formally disallowed on strict steric considerations correspond to the gamma turn (C7 axial) conformation (Bin 12) and the i + 1 position of Type II turns (Bin 13). Of the 516 non-Gly residues in Bin 13, 384 occupied the i + 1 position of Type II turns. Further examination of these turn segments revealed a high propensity to occur at the N-terminus of helices and as a tight turn in hairpins. The strand-helix motif with the Type II turn as a connecting element was also found in as many as 57 examples. Proteins 2014; 82:1101-1112. (c) 2013 Wiley Periodicals, Inc.

The Glycosylated Rv1860 Protein of Mycobacterium tuberculosis Inhibits Dendritic Cell Mediated TH1 and TH17 Polarization of T Cells and Abrogates Protective Immunity Conferred by BCG

We previously reported interferon gamma secretion by human CD4(+) and CD8(+) T cells in response to recombinant E. coli-expressed Rv1860 protein of Mycobacterium tuberculosis (MTB) as well as protection of guinea pigs against a challenge with virulent MTB following prime-boost immunization with DNA vaccine and poxvirus expressing Rv1860. In contrast, a Statens Serum Institute Mycobacterium bovis BCG (BCG-SSI) recombinant expressing MTB Rv1860 (BCG-TB1860) showed loss of protective ability compared to the parent BCG strain expressing the control GFP protein (BCG-GFP). Since Rv1860 is a secreted mannosylated protein of MTB and BCG, we investigated the effect of BCG-TB1860 on innate immunity. Relative to BCG-GFP, BCG-TB1860 effected a significant near total reduction both in secretion of cytokines IL-2, IL-12p40, IL-12p70, TNF-alpha, IL-6 and IL-10, and up regulation of co-stimulatory molecules MHC-II, CD40, CD54, CD80 and CD86 by infected bone marrow derived dendritic cells (BMDC), while leaving secreted levels of TGF-beta unchanged. These effects were mimicked by BCG-TB1860His which carried a 6-Histidine tag at the C-terminus of Rv1860, killed sonicated preparations of BCG-TB1860 and purified H37Rv-derived Rv1860 glycoprotein added to BCG-GFP, but not by E. coli-expressed recombinant Rv1860. Most importantly, BMDC exposed to BCG-TB1860 failed to polarize allogeneic as well as syngeneic T cells to secrete IFN-gamma and IL-17 relative to BCG-GFP. Splenocytes from mice infected with BCG-SSI showed significantly less proliferation and secretion of IL-2, IFN-gamma and IL-17, but secreted higher levels of IL-10 in response to in vitro restimulation with BCG-TB1860 compared to BCG-GFP. Spleens from mice infected with BCG-TB1860 also harboured significantly fewer DC expressing MHC-II, IL-12, IL-2 and TNF-alpha compared to mice infected with BCG-GFP. Glycoproteins of MTB, through their deleterious effects on DC may thus contribute to suppress the generation of a TH1- and TH17-dominated adaptive immune response that is vital for protection against tuberculosis.

Groundnut Bud Necrosis Virus Encoded NSm Associates with Membranes via Its C-Terminal Domain

Groundnut Bud Necrosis Virus (GBNV) is a tripartite ambisense RNA plant virus that belongs to serogroup IV of Tospovirus genus. Non-Structural protein-m (NSm), which functions as movement protein in tospoviruses, is encoded by the M RNA. In this communication, we demonstrate that despite the absence of any putative transmembrane domain, GBNV NSm associates with membranes when expressed in E. coli as well as in N. benthamiana. Incubation of refolded NSm with liposomes ranging in size from 200-250 nm resulted in changes in the secondary and tertiary structure of NSm. A similar behaviour was observed in the presence of anionic and zwitterionic detergents. Furthermore, the morphology of the liposomes was found to be modified in the presence of NSm. Deletion of coiled coil domain resulted in the inability of in planta expressed NSm to interact with membranes. Further, when the C-terminal coiled coil domain alone was expressed, it was found to be associated with membrane. These results demonstrate that NSm associates with membranes via the C-terminal coiled coil domain and such an association may be important for movement of viral RNA from cell to cell.

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.