Network properties of protein-decoy structures

Convergence of the vast sequence space of proteins into a highly restricted fold/conformational space suggests a simple yet unique underlying mechanism of protein folding that has been the subject of much debate in the last several decades. One of the major challenges related to the understanding of protein folding or in silico protein structure prediction is the discrimination of non-native structures/decoys from the native structure. Applications of knowledge-based potentials to attain this goal have been extensively reported in the literature. Also, scoring functions based on accessible surface area and amino acid neighbourhood considerations were used in discriminating the decoys from native structures. In this article, we have explored the potential of protein structure network (PSN) parameters to validate the native proteins against a large number of decoy structures generated by diverse methods. We are guided by two principles: (a) the PSNs capture the local properties from a global perspective and (b) inclusion of non-covalent interactions, at all-atom level, including the side-chain atoms, in the network construction accommodates the sequence dependent features. Several network parameters such as the size of the largest cluster, community size, clustering coefficient are evaluated and scored on the basis of the rank of the native structures and the Z-scores. The network analysis of decoy structures highlights the importance of the global properties contributing to the uniqueness of native structures. The analysis also exhibits that the network parameters can be used as metrics to identify the native structures and filter out non-native structures/decoys in a large number of data-sets; thus also has a potential to be used in the protein `structure prediction' problem.

A touch of history and a peep into the future of the lipid-quinone known as coenzyme Q and ubiquinone

Coenzyme Q (ubiquinone), a fully substituted benzoquinone with polyprenyl side chain, participates in many cellular redox activities. Paradoxically it was discovered only in 1957, albeit being ubiquitous. It required a person, F. L. Crane, a place, Enzyme Institute, Madison, USA, and a time when D. E. Green was directing vigorous research on mitochondria. Located at the transition of 2-electron flavoproteins and 1-electron cytochrome carriers, it facilitates electron transfer through the elegant Q-cycle in mitochondria to reduce O-2 to H2O, and to H2O2, now a significant signal-transducing agent, as a minor activity in shunt pathway (animals) and alternative oxidase (plants). The ability to form Q-radical by losing an electron and a proton was ingeniously used by Mitchell to explain the formation of the proton gradient, considered the core of energy transduction, and also in acidification in vacuoles. Known to be a mobile membrane constituent (microsomes, plasma membrane and Golgi apparatus), allowing it to reach multiple sites, coenzyme Q is expected to have other activities. Coenzyme Q protects circulating lipoproteins being a better lipid antioxidant than even vitamin E. Binding to proteins such as QPS, QPN, QPC and uncoupling protein in mitochondria, QA and QB in the reaction centre in R. sphaeroides, and disulfide bond-forming protein in E. coli (possibly also in Golgi), coenzyme Q acquires selective functions. A characteristic of orally dosed coenzyme Q is its exclusive absorption into the liver, but not the other tissues. This enrichment of Q is accompanied by significant decrease of blood pressure and of serum cholesterol. Inhibition of formation of mevalonate, the common precursor in the branched isoprene pathway, by the minor product, coenzyme Q, decreases the major product, cholesterol. Relaxation of contracted arterial smooth muscle by a side-chain truncated product of coenzyme Q explains its effect of decreasing blood pressure. Extensive clinical studies carried out on oral supplements of coenzyine Q, initially by K. Folkers and Y. Yamamura and followed many others, revealed a large number of beneficial effects, significantly in cardiovascular diseases. Such a variety of effects by this lipid quinone cannot depend on redox activity alone. The fat-soluble vitamins (A, D, E and K) that bear structural relationship with coenzyme Q are known to be active in their polar forms. A vignette of modified forms of coenzyme Q taking active role in its multiple effects is emerging.

Photoresist functionalisation method for high-density protein microarrays using photolithography

Since the last decade, there is a growing need for patterned biomolecules for various applications ranging from diagnostic devices to enabling fundamental biological studies with high throughput. Protein arrays facilitate the study of protein-protein, protein-drug or protein-DNA interactions as well as highly multiplexed immunosensors based on antibody-antigen recognition. Protein microarrays are typically fabricated using piezoelectric inkjet printing with resolution limit of similar to 70-100 mu m limiting the array density. A considerable amount of research has been done on patterning biomolecules using customised biocompatible photoresists. Here, a simple photolithographic process for fabricating protein microarrays on a commercially available diazo-naphthoquinone-novolac-positive tone photoresist functionalised with 3-aminopropyltriethoxysilane is presented. The authors demonstrate that proteins immobilised using this procedure retain their activity and therefore form functional microarrays with the array density limited only by the resolution of lithography, which is more than an order of magnitude compared with inkjet printing. The process described here may be useful in the integration of conventional semiconductor manufacturing processes with biomaterials relevant for the creation of next-generation bio-chips.

Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X-ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra-molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. © 2012 Wiley Periodicals, Inc.

Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X-ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra-molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. (c) 2012 Wiley Periodicals, Inc.

Phospholipid Mediated Activation of Calcium Dependent Protein Kinase 1 (CaCDPK1) from Chickpea: A New Paradigm of Regulation

Phospholipids, the major structural components of membranes, can also have functions in regulating signaling pathways in plants under biotic and abiotic stress. The effects of adding phospholipids on the activity of stress-induced calcium dependent protein kinase (CaCDPK1) from chickpea are reported here. Both autophosphorylation as well as phosphorylation of the added substrate were enhanced specifically by phosphatidylcholine and to a lesser extent by phosphatidic acid, but not by phosphatidylethanolamine. Diacylgylerol, the neutral lipid known to activate mammalian PKC, stimulated CaCDPK1 but at higher concentrations. Increase in V-max of the enzyme activity by these phospholipids significantly decreased the K-m indicating that phospholipids enhance the affinity towards its substrate. In the absence of calcium, addition of phospholipids had no effect on the negligible activity of the enzyme. Intrinsic fluorescence intensity of the CaCDPK1 protein was quenched on adding PA and PC. Higher binding affinity was found with PC (K-1/2 = 114 nM) compared to PA (K-1/2 = 335 nM). We also found that the concentration of PA increased in chickpea plants under salt stress. The stimulation by PA and PC suggests regulation of CaCDPK1 by these phospholipids during stress response.

Structural and molecular basis of interaction of HCV non-structural protein 5A with human casein kinase 1 alpha and PKR

Background: Interaction of non-structural protein 5A (NS5A) of Hepatitis C virus (HCV) with human kinases namely, casein kinase 1 alpha (ck1 alpha) and protein kinase R (PKR) have different functional implications such as regulation of viral replication and evasion of interferon induced immune response respectively. Understanding the structural and molecular basis of interactions of the viral protein with two different human kinases can be useful in developing strategies for treatment against HCV. Results: Serine 232 of NS5A is known to be phosphorylated by human ck1 alpha. A structural model of NS5A peptide containing phosphoacceptor residue Serine 232 bound to ck1 alpha has been generated using the known 3-D structures of kinase-peptide complexes. The substrate interacting residues in ck1 alpha has been identified from the model and these are found to be conserved well in the ck1 family. ck1 alpha - substrate peptide complex has also been used to understand the structural basis of association between ck1 alpha and its other viral stress induced substrate, tumour suppressor p53 transactivation domain which has a crystal structure available. Interaction of NS5A with another human kinase PKR is primarily genotype specific. NS5A from genotype 1b has been shown to interact and inhibit PKR whereas NS5A from genotype 2a/3a are unable to bind and inhibit PKR efficiently. This is one of the main reasons for the varied response to interferon therapy in HCV patients across different genotypes. Using PKR crystal structure, sequence alignment and evolutionary trace analysis some of the critical residues responsible for the interaction of NS5A 1b with PKR have been identified. Conclusions: The substrate interacting residues in ck1 alpha have been identified using the structural model of kinase substrate peptide. The PKR interacting NS5A 1b residues have also been predicted using PKR crystal structure, NS5A sequence analysis along with known experimental results. Functional significance and nature of interaction of interferon sensitivity determining region and variable region 3 of NS5A in different genotypes with PKR which was experimentally shown are also supported by the findings of evolutionary trace analysis. Designing inhibitors to prevent this interaction could enable the HCV genotype 1 infected patients respond well to interferon therapy.

Site-specific N-Linked Glycosylation of Receptor Guanylyl Cyclase C Regulates Ligand Binding, Ligand-mediated Activation and Interaction with Vesicular Integral Membrane Protein 36, VIP36

Guanylyl cyclase C (GC-C) is a multidomain, membrane-associated receptor guanylyl cyclase. GC-C is primarily expressed in the gastrointestinal tract, where it mediates fluid-ion homeostasis, intestinal inflammation, and cell proliferation in a cGMP-dependent manner, following activation by its ligands guanylin, uroguanylin, or the heat-stable enterotoxin peptide (ST). GC-C is also expressed in neurons, where it plays a role in satiation and attention deficiency/hyperactive behavior. GC-C is glycosylated in the extracellular domain, and differentially glycosylated forms that are resident in the endoplasmic reticulum (130 kDa) and the plasma membrane (145 kDa) bind the ST peptide with equal affinity. When glycosylation of human GC-C was prevented, either by pharmacological intervention or by mutation of all of the 10 predicted glycosylation sites, ST binding and surface localization was abolished. Systematic mutagenesis of each of the 10 sites of glycosylation in GC-C, either singly or in combination, identified two sites that were critical for ligand binding and two that regulated ST-mediated activation. We also show that GC-C is the first identified receptor client of the lectin chaperone vesicular integral membrane protein, VIP36. Interaction with VIP36 is dependent on glycosylation at the same sites that allow GC-C to fold and bind ligand. Because glycosylation of proteins is altered in many diseases and in a tissue-dependent manner, the activity and/or glycan-mediated interactions of GC-C may have a crucial role to play in its functions in different cell types.

Dramatic Structural Changes Resulting from the Loss of a Crucial Hydrogen Bond in the Hinge Region Involved in C-Terminal Helix Swapping in SurE: A Survival Protein from Salmonella typhimurium

Domain swapping is an interesting feature of some oligomeric proteins in which each protomer of the oligomer provides an identical surface for exclusive interaction with a segment or domain belonging to another protomer. Here we report results of mutagenesis experiments on the structure of C-terminal helix swapped dimer of a stationary phase survival protein from Salmonella typhimurium (StSurE). Wild type StSurE is a dimer in which a large helical segment at the C-terminus and a tetramerization loop comprising two beta strands are swapped between the protomers. Key residues in StSurE that might promote C-terminal helix swapping were identified by sequence and structural comparisons. Three mutants in which the helix swapping is likely to be avoided were constructed and expressed in E. coli. Three-dimensional X-ray crystal structures of the mutants H234A and D230A/H234A could be determined at 2.1 angstrom and 2.35 angstrom resolutions, respectively. Contrary to expectations, helix swapping was mostly retained in both the mutants. The loss of the crucial D230 OD2- H234 NE2 hydrogen bond (2.89 angstrom in the wild type structure) in the hinge region was compensated by new inter and intra-chain interactions. However, the two fold molecular symmetry was lost and there were large conformational changes throughout the polypeptide. In spite of these changes, the dimeric structure and an approximate tetrameric organization were retained, probably due to the interactions involving the tetramerization loop. Mutants were mostly functionally inactive, highlighting the importance of precise inter-subunit interactions for the symmetry and function of StSurE.

Network properties of decoys and CASP predicted models: a comparison with native protein structures

Protein structure space is believed to consist of a finite set of discrete folds, unlike the protein sequence space which is astronomically large, indicating that proteins from the available sequence space are likely to adopt one of the many folds already observed. In spite of extensive sequence-structure correlation data, protein structure prediction still remains an open question with researchers having tried different approaches (experimental as well as computational). One of the challenges of protein structure prediction is to identify the native protein structures from a milieu of decoys/models. In this work, a rigorous investigation of Protein Structure Networks (PSNs) has been performed to detect native structures from decoys/ models. Ninety four parameters obtained from network studies have been optimally combined with Support Vector Machines (SVM) to derive a general metric to distinguish decoys/models from the native protein structures with an accuracy of 94.11%. Recently, for the first time in the literature we had shown that PSN has the capability to distinguish native proteins from decoys. A major difference between the present work and the previous study is to explore the transition profiles at different strengths of non-covalent interactions and SVM has indeed identified this as an important parameter. Additionally, the SVM trained algorithm is also applied to the recent CASP10 predicted models. The novelty of the network approach is that it is based on general network properties of native protein structures and that a given model can be assessed independent of any reference structure. Thus, the approach presented in this paper can be valuable in validating the predicted structures. A web-server has been developed for this purpose and is freely available at http://vishgraph.mbu.iisc.ernet.in/GraProStr/PSN-QA.html.

Algal biofuel: bountiful lipid from Chlorococcum sp proliferating in municipal wastewater

Algae biofuel have emerged as viable renewable energy sources and are the potential alternatives to fossil-based fuels in recent times. Algae have the potential to generate significant quantities of commercially viable biofuel apart from treating wastewater. Three algal species, viz. Chlorococcum sp., Microcystis sp. and Phormidium sp. proliferating in wastewater ponds were isolated and cultured in the laboratory myxotrophically under similar wastewater conditions. Chlorococcum sp. attained a mean biomass productivity of 0.09 g. I(-1)d(-1) with the maximum `biomass density of 1.33 g I-1 and comparatively higher lipid content of 30.55% (w/w) on the ninth day of the culture experiment. Under similar conditions Microcystis sp. and Phormidium sp. attained mean biomass productivities of 0.058 and 0.063 g I-1 d(-1) with a total lipid content of 8.88% and 18.66% respectively. Biochemical composition (carbohydrates, proteins, lipids and phosphates) variations and lipid accumulation studies were performed by comparison of the ratios of carbohydrate to protein; lipid to protein (L/P) and lipid to phosphates using attenuated total reflectance-Fourier transform infrared spectroscopy which showed higher L/P ratio during the stationary phase of algal growth. Composition analysis of fatty acid methyl ester has been performed using gas chromatography and mass spectrometry. Chlorococcum sp. with higher productivity and faster growth rate has higher lipid content with about 67% of saturated fatty acid dominated by palmitate (36.3%) followed by an unsaturate as linoleate (14%) and has proved to be an economical and viable feedstock for biofuel production compared to the other wastewater-grown species.

Redox regulation of protein tyrosine phosphatase 1B (PTP1B): Importance of steric and electronic effects on the unusual cyclization of the sulfenic acid intermediate to a sulfenyl amide

The redox regulation of protein tyrosine phosphatase 1B (PTP1B) via the unusual transformation of its sulfenic acid (PTP1B-SOH) to a cyclic sulfenyl amide intermediate is studied by using small molecule chemical models. These studies suggest that the sulfenic acids derived from the H2O2-mediated reactions o-amido thiophenols do not efficiently cyclize to sulfenyl amides and the sulfenic acids produced in situ can be trapped by using methyl iodide. Theoretical calculations suggest that the most stable conformer of such sulfenic acids are stabilized by n(O) -> sigma* (S-OH) orbital interactions, which force the -OH group to adopt a position trans to the S center dot center dot center dot O interaction, leading to an almost linear arrangement of the O center dot center dot center dot S-O moiety and this may be the reason for the slow cyclization of such sulfenic acids to their corresponding sulfenyl amides. On the other hand, additional substituents at the 6-position of o-amido phenylsulfenic acids that can induce steric environment and alter the electronic properties around the sulfenic acid moiety by S center dot center dot center dot N or S center dot center dot center dot O nonbonded interactions destabilize the sulfenic acids by inducing strain in the molecule. This may lead to efficient the cyclization of such sulfenic acids. This model study suggests that the amino acid residues in the close proximity of the sulfenic acid moiety in PTP1B may play an important role in the cyclization of PTP1B-SOH to produce the corresponding sulfenyl amide.

DNA and protein targeting 1,2,4-triazole based water soluble dinickel(II) complexes enhances antiproliferation and lactate dehydrogenase inhibition

Water soluble dinickel(II) complexes Ni-2(L)(2)(1-2)](NO3)(4) (1-2), where L1-2 are triazole based dinucleating ligands, were synthesized and characterized. The DNA binding, protein binding, DNA hydrolysis and anticancer properties were investigated. The interactions of complexes 1 and 2 with calf thymus DNA were studied by spectroscopic techniques, including absorption and fluorescence spectroscopy. The DNA binding constant values of the complexes 1 and 2 were found to be 2.36 x 10(5) and 4.87 x 10(5) M-1 and the binding affinities are in the following order: 2 > 1. Both the dinickel(II) complexes 1 and 2, promoted the hydrolytic cleavage of plasmid pBR322 DNA under both anaerobic and aerobic conditions. Kinetic data for DNA hydrolysis promoted by 1 and 2 under physiological conditions give the observed rate constants (k(obs)) of 5.05 +/- 0.2 and 5.65 +/- 0.1 h(-1), respectively, which shows 10(8)-fold rate acceleration over the uncatalyzed reaction of ds-DNA. Meanwhile, the interactions of the complex with BSA have also been studied by spectroscopy. Both the complexes 1 and 2 display strong binding propensity and the binding constant (K-b), number of binding sites (n) were obtained are 0.71 x 10(6) 1.47] and 5.62 x 10(6) 1.98] M-1, respectively. The complexes 1 and 2 also promoted the apoptosis against human carcinoma (HeLa, and BeWo) cancer cells. Cytotoxicity of the complexes was further confirmed by lactate dehydrogenase enzyme level in cancer cell lysate and content media. (c) 2013 Elsevier Ltd. All rights reserved.

An automated approach to network features of protein structure ensembles

Network theory applied to protein structures provides insights into numerous problems of biological relevance. The explosion in structural data available from PDB and simulations establishes a need to introduce a standalone-efficient program that assembles network concepts/parameters under one hood in an automated manner. Herein, we discuss the development/application of an exhaustive, user-friendly, standalone program package named PSN-Ensemble, which can handle structural ensembles generated through molecular dynamics (MD) simulation/NMR studies or from multiple X-ray structures. The novelty in network construction lies in the explicit consideration of side-chain interactions among amino acids. The program evaluates network parameters dealing with topological organization and long-range allosteric communication. The introduction of a flexible weighing scheme in terms of residue pairwise cross-correlation/interaction energy in PSN-Ensemble brings in dynamical/chemical knowledge into the network representation. Also, the results are mapped on a graphical display of the structure, allowing an easy access of network analysis to a general biological community. The potential of PSN-Ensemble toward examining structural ensemble is exemplified using MD trajectories of an ubiquitin-conjugating enzyme (UbcH5b). Furthermore, insights derived from network parameters evaluated using PSN-Ensemble for single-static structures of active/inactive states of 2-adrenergic receptor and the ternary tRNA complexes of tyrosyl tRNA synthetases (from organisms across kingdoms) are discussed. PSN-Ensemble is freely available from http://vishgraph.mbu.iisc.ernet.in/PSN-Ensemble/psn_index.html.