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.

Enhancing Surface Coverage and Growth in Layer-by-Layer Assembly of Protein Nanoparticles

Thin films of bovine serum albumin (BSA) nanoparticles are fabricated via layer-by-layer assembly. The surface of BSA nanoparticles have two oppositely acting functional groups on the surface: amine (NH2) and carboxylate (COO-). The protonation and deprotonation of these functional groups at different pH vary the charge density on the particle surface, and entirely different growth can be observed by varying the nature of the complementary polymer and the pH of the particles. The complementary polymers used in this study are poly(dimethyldiallylammonium chloride) (PDDAC) and poly(acrylic acid) (PAA). The assembly of BSA nanoparticles based on electrostatic interaction with PDDAC suffers from the poor loading of the nanoparticles. The assembly with PAA aided by a hydrogen bonding interaction shows tremendous improvement in the growth of the assembly over PDDAC. Moreover, the pH of the BSA nanoparticles was observed to affect the loading of nanoparticles in the LbL assembly with PAA significantly.

Polyelectrolyte/silver nanocomposite multilayer films as multifunctional thin film platforms for remote activated protein and drug delivery

We demonstrate a nanoparticle loading protocol to develop a transparent, multifunctional polyelectrolyte multilayer film for externally activated drug and protein delivery. The composite film was designed by alternate adsorption of poly(allylamine hydrochloride) (PAH) and dextran sulfate (DS) on a glass substrate followed by nanoparticle synthesis through a polyol reduction method. The films showed a uniform distribution of spherical silver nanoparticles with an average diameter of 50 +/- 20 nm, which increased to 80 +/- 20 nm when the AgNO3 concentration was increased from 25 to 50 mM. The porous and supramolecular structure of the polyelectrolyte multilayer film was used to immobilize ciprofloxacin hydrochloride (CH) and bovine serum albumin (BSA) within the polymeric network of the film. When exposed to external triggers such as ultrasonication and laser light the loaded films were ruptured and released the loaded BSA and CH. The release of CH is faster than that of BSA due to a higher diffusion rate. Circular dichroism measurements confirmed that there was no significant change in the conformation of released BSA in comparison with native BSA. The fabricated films showed significant antibacterial activity against the bacterial pathogen Staphylococcus aureus. Applications envisioned for such drug-loaded films include drug and vaccine delivery through the transdermal route, antimicrobial or anti-inflammatory coatings on implants and drug-releasing coatings for stents. (C) 2013 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

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.

Evaluation of Brugia malayi sheath protein (Shp-1) as a diagnostic antigen for human lymphatic filariasis

Lymphatic filariasis is the second leading cause of permanent long-term disability globally and control of this disease needs efficient diagnostic methods. In this study, abundantly expressing microfilarial sheath protein (Shp-1) from Brugia malayi was characterized as a filarial diagnostic candidate using samples from different clinical population. Monoclonal antibodies were developed against E. coil expressed recombinant Shp-1 in order to assess its efficiency in filarial antigen detection assay system. Endemic Normal (EN, n = 170), asymptomatic microfilaeremics (MF, n = 65), symptomatic chronic pathology (CP, n = 45) and non endemic normal (NEN, n = 10) sera were analyzed by antigen capture enzyme-linked immunosorbent assay. Of the 290 individuals, all MF individuals (both brugian and bancroftian) were positive in this assay followed by CP and EN. When compared with SXP-1 and Og4C3 antigen assays, all assays detected Wb MF correctly, Bm MF was detected by Shp-1 and SXP-1 assays, and only Shp-1 was able to detect EN (12%) and CP (29%). Results showed that this assay may be useful for monitoring prior to mass drug administration. (c) 2014 Elsevier Inc. All rights reserved.

Specific Sequence of a Beta Turn in Human La Protein May Contribute to Species Specificity of Hepatitis C Virus

Human La protein is known to be an essential host factor for translation and replication of hepatitis C virus (HCV) RNA. Previously, we have demonstrated that residues responsible for interaction of human La protein with the HCV internal ribosomal entry site (IRES) around the initiator AUG within stem-loop IV form a beta-turn in the RNA recognition motif (RRM) structure. In this study, sequence alignment and mutagenesis suggest that the HCV RNA-interacting beta-turn is conserved only in humans and chimpanzees, the species primarily known to be infected by HCV. A 7-mer peptide corresponding to the HCV RNA-interacting region of human La inhibits HCV translation, whereas another peptide corresponding to the mouse La sequence was unable to do so. Furthermore, IRES-mediated translation was found to be significantly high in the presence of recombinant human La protein in vitro in rabbit reticulocyte lysate. We observed enhanced replication with HCV subgenomic and full-length replicons upon overexpression of either human La protein or a chimeric mouse La protein harboring a human La beta-turn sequence in mouse cells. Taken together, our results raise the possibility of creating an immunocompetent HCV mouse model using human-specific cell entry factors and a humanized form of La protein.