Solvent Sensitivity of Protein Unfolding: Dynamical Study of Chicken Villin Headpiece Subdomain in Water Ethanol Binary Mixture

We carry out a series of long atomistic molecular dynamics simulations to study the unfolding of a small protein, chicken villin headpiece (HP-36), in water-ethanol (EtOH) binary mixture. The prime objective of this work is to explore the sensitivity of protein unfolding dynamics toward increasing concentration of the cosolvent and unravel essential features of intermediates formed in search of a dynamical pathway toward unfolding. In water ethanol binary mixtures, HP-36 is found to unfold partially, under ambient conditions, that otherwise requires temperature as high as similar to 600 K to denature in pure aqueous solvent. However, an interesting course of pathway is observed to be followed in the process, guided by the formation of unique intermediates. The first step of unfolding is essentially the separation of the cluster formed by three hydrophobic (phenylalanine) residues, namely, Phe-7, Phe-11, and Phe-18, which constitute the hydrophobic core, thereby initiating melting of helix-2 of the protein. The initial steps are similar to temperature-induced unfolding as well as chemical unfolding using DMSO as cosolvent. Subsequent unfolding steps follow a unique path. As water-ethanol shows composition-dependent anomalies, so do the details of unfolding dynamics. With an increase in cosolvent concentration, different partially unfolded intermediates are found to be formed. This is reflected in a remarkable nonmonotonic composition dependence of several order parameters, including fraction of native contacts and protein-solvent interaction energy. The emergence of such partially unfolded states can be attributed to the preferential solvation of the hydrophobic residues by the ethyl groups of ethanol. We further quantify the local dynamics of unfolding by using a Marcus-type theory.

mulPBA: an efficient multiple protein structure alignment method based on a structural alphabet

The increasing number of available protein structures requires efficient tools for multiple structure comparison. Indeed, multiple structural alignments are essential for the analysis of function, evolution and architecture of protein structures. For this purpose, we proposed a new web server called multiple Protein Block Alignment (mulPBA). This server implements a method based on a structural alphabet to describe the backbone conformation of a protein chain in terms of dihedral angles. This sequence-like' representation enables the use of powerful sequence alignment methods for primary structure comparison, followed by an iterative refinement of the structural superposition. This approach yields alignments superior to most of the rigid-body alignment methods and highly comparable with the flexible structure comparison approaches. We implement this method in a web server designed to do multiple structure superimpositions from a set of structures given by the user. Outputs are given as both sequence alignment and superposed 3D structures visualized directly by static images generated by PyMol or through a Jmol applet allowing dynamic interaction. Multiple global quality measures are given. Relatedness between structures is indicated by a distance dendogram. Superimposed structures in PDB format can be also downloaded, and the results are quickly obtained. mulPBA server can be accessed at www.dsimb.inserm.fr/dsimb_tools/mulpba/.

Metal nanoparticles as better protein-aggregation prevention agents than chaperones

Prevention or suppression of protein aggregation is of great importance in the context of protein storage, transportation and delivery. Traditionally chaperones or other chemically active agents are used to stop or diffuse native protein aggregation. We have used gold nanoparticles to prevent thermal aggregation of alcohol dehydrogenase (ADH), a protein that maintains the alcohol level in the liver and stomach. A light-scattering assay has been used to investigate the effect of gold nanoparticles on thermal aggregation of ADH and the result of our study has been summarized in Fig. 1. The scattered light intensity from the solution containing ADH decreases when 45 nm gold nanoparticles are added prior to heating (thermal denaturation) the solution, which indicates prevention of aggregation. The aggregation of the protein is suppressed to the extent of 96% with picomolar concentration of 45 nm gold nanoparticles while micromolar amounts of other proteins and biological substances are necessary to achieve the same effect. The extent varies with the size and the concentration of the gold NPs for the same protein concentration.

Molecular Dynamics Perspective on the Protein Thermal Stability: A Case Study Using SAICAR Synthetase

The enzyme SAICAR synthetase ligates aspartate with CAIR (5'-phosphoribosyl-4-carboxy-5-aminoimidazole) forming SAICAR (5-amino-4-imidazole-N-succinocarboxamide ribonucleotide) in the presence of ATP. In continuation with our previous study on the thermostability of this enzyme in hyper-/thermophiles based on the structural aspects, here, we present the dynamic aspects that differentiate the mesophilic (E. coli, E. chaffeensis), thermophilic (G. kaustophilus), and hyperthermophilic (M. jannaschii, P. horikoshii) SAICAR synthetases by carrying out a total of 11 simulations. The five functional dimers from the above organisms were simulated using molecular dynamics for a period of 50 ns each at 300 K, 363 K, and an additional simulation at 333 K for the thermophilic protein. The basic features like root-mean-square deviations, root-mean-square fluctuations, surface accessibility, and radius of gyration revealed the instability of mesophiles at 363 K. Mean square displacements establish the reduced flexibility of hyper-/thermophiles at all temperatures. At the simulations time scale considered here, the long-distance networks are considerably affected in mesophilic structures at 363 K. In mesophiles, a comparatively higher number of short-lived (having less percent existence time) C alpha, hydrogen bonds, hydrophobic interactions are formed, and long-lived (with higher percentage existence time) contacts are lost. The number of time-averaged salt-bridges is at least 2-fold higher in hyperthermophiles at 363 K. The change in surface accessibility of salt-bridges at 363 K from 300 K is nearly doubled in mesophilic protein compared to proteins from other temperature classes.

Investigations of Ramachandran disallowed conformations in protein domain families

In peptide and protein structures, occurrence of (phi,psi.) angles in the disallowed region of the Ramachandran map almost always suggests local regions of error or poor accuracy. However, very rarely genuine disallowed conformations occur as noted in the current study in proteins of known structure available at ultra-high resolution (<= 1.2 (A) over circle). In the current work, extent of conservation of genuine disallowed conformations in homologous proteins of known structures has been analyzed. From a dataset of 124 protein domain families, with structure of at least one constituent member in each family available at a resolution of 1.2 (A) over circle or better, we have analyzed the conservation of 221 disallowed conformations. It is observed that the disallowed conformation is only moderately conservedin protein domain families. In the gross dataset no particular residue type adopting disallowed conformation elicit high conservation of residue type though there are alignment positions in the dataset with complete conservation of both the residue type and the disallowed conformation. Conserved disallowed conformation in protein domain families play biologically significant role in roughly 50% of the cases. The residues with the disallowed conformation or its flanking residues are often located within or around the functional site of the protein. (C) 2013 Elsevier B.V. All rights reserved.

Insulin like growth factor binding protein 4 promotes GBM progression and regulates key factors involved in EMT and invasion

Insulin like growth factor binding protein 4 (IGFBP4) regulates growth and development of tissues and organs by negatively regulating IGF signaling. Among most cancers, IGFBP4 has growth inhibitory role and reported as a down-regulated gene, except for renal cell carcinoma, wherein IGFBP4 promotes tumor progression. IGFBP4 expression has been shown to be higher in increasing grades of astrocytoma. However, the functional role of IGFBP4 in gliomas has not been explored. Surgical biopsies of 20 normal brain and 198 astrocytoma samples were analyzed for IGFBP4 expression by qRT-PCR. Highest expression of IGFBP4 mRNA was seen in GBM tumors compared to control brain tissues (median log2 of 2.035, p < 0.0001). Immunohistochemical analysis of 53 tissue samples revealed predominant nuclear staining of IGFBP4, seen maximally in GBMs when compared to DA and AA tumors (median LI = 29.12 +/- A 16.943, p < 0.001). Over expression of IGFBP4 in U343 glioma cells resulted in up-regulation of molecules involved in tumor growth, EMT and invasion such as pAkt, pErk, Vimentin, and N-cadherin and down-regulation of E-cadherin. Functionally, IGFBP4 over expression in these cells resulted in increased proliferation, migration and invasion as assessed by MTT, transwell migration, and Matrigel invasion assays. These findings were confirmed upon IGFBP4 knockdown in U251 glioma cells. Our data suggest a pro-tumorigenic role for IGFBP4 in glioma.

Transcriptional Repression of Tumor Suppressor CDC73, Encoding an RNA Polymerase II Interactor, by Wilms Tumor 1 Protein (WT1) Promotes Cell Proliferation IMPLICATION FOR CANCER THERAPEUTICS

The Wilms tumor 1 gene (WT1) can either repress or induce the expression of genes. Inconsistent with its tumor suppressor role, elevated WT1 levels have been observed in leukemia and solid tumors. WT1 has also been suggested to act as an oncogene by inducing the expression of MYC and BCL-2. However, these are only the correlational studies, and no functional study has been performed to date. Consistent with its tumor suppressor role, CDC73 binds to RNA polymerase II as part of a PAF1 transcriptional regulatory complex and causes transcriptional repression of oncogenes MYC and CCND1. It also represses beta-catenin-mediated transcription. Based on the reduced level of CDC73 in oral squamous cell carcinoma (OSCC) samples in the absence of loss-of-heterozygosity, promoter methylation, and mutations, we speculated that an inhibitory transcription factor is regulating its expression. The bioinformatics analysis predicted WT1 as an inhibitory transcription factor to regulate the CDC73 level. Our results showed that overexpression of WT1 decreased CDC73 levels and promoted proliferation of OSCC cells. ChIP and EMSA results demonstrated binding of WT1 to the CDC73 promoter. The 5-azacytidine treatment of OSCC cells led to an up-regulation of WT1 with a concomitant down-regulation of CDC73, further suggesting regulation of CDC73 by WT1. Exogenous CDC73 attenuated the protumorigenic activity of WT1 by apoptosis induction. An inverse correlation between expression levels of CDC73 and WT1 was observed in OSCC samples. These observations indicated that WT1 functions as an oncogene by repressing the expression of CDC73 in OSCC. We suggest that targeting WT1 could be a therapeutic strategy for cancer, including OSCC.

Filling-in Void and Sparse Regions in Protein Sequence Space by Protein-Like Artificial Sequences Enables Remarkable Enhancement in Remote Homology Detection Capability

Protein functional annotation relies on the identification of accurate relationships, sequence divergence being a key factor. This is especially evident when distant protein relationships are demonstrated only with three-dimensional structures. To address this challenge, we describe a computational approach to purposefully bridge gaps between related protein families through directed design of protein-like ``linker'' sequences. For this, we represented SCOP domain families, integrated with sequence homologues, as multiple profiles and performed HMM-HMM alignments between related domain families. Where convincing alignments were achieved, we applied a roulette wheel-based method to design 3,611,010 protein-like sequences corresponding to 374 SCOP folds. To analyze their ability to link proteins in homology searches, we used 3024 queries to search two databases, one containing only natural sequences and another one additionally containing designed sequences. Our results showed that augmented database searches showed up to 30% improvement in fold coverage for over 74% of the folds, with 52 folds achieving all theoretically possible connections. Although sequences could not be designed between some families, the availability of designed sequences between other families within the fold established the sequence continuum to demonstrate 373 difficult relationships. Ultimately, as a practical and realistic extension, we demonstrate that such protein-like sequences can be ``plugged-into'' routine and generic sequence database searches to empower not only remote homology detection but also fold recognition. Our richly statistically supported findings show that complementary searches in both databases will increase the effectiveness of sequence-based searches in recognizing all homologues sharing a common fold. (C) 2013 Elsevier Ltd. All rights reserved.

The N-terminal region containing the zinc finger domain of tobacco streak virus coat protein is essential for the formation of virus-like particles

Tobacco streak virus (TSV), a member of the genus Ilarvirus (family Bromoviridae), has a tripartite genome and forms quasi-isometric virions. All three viral capsids, encapsidating RNA 1, RNA 2 or RNA 3 and subgenomic RNA 4, are constituted of a single species of coat protein (CP). Formation of virus-like particles (VLPs) could be observed when the TSV CP gene was cloned and the recombinant CP (rCP) was expressed in E. coli. TSV VLPs were found to be stabilized by Zn2+ ions and could be disassembled in the presence of 500 mM CaCl2. Mutational analysis corroborated previous studies that showed that an N-terminal arginine-rich motif was crucial for RNA binding; however, the results presented here demonstrate that the presence of RNA is not a prerequisite for assembly of TSV VLPs. Instead, the N-terminal region containing the zinc finger domain preceding the arginine-rich motif is essential for assembly of these VLPs.