Improvement of 3D protein models using multiple templates guided by single-template model quality assessment

Motivation: Modelling the 3D structures of proteins can often be enhanced if more than one fold template is used during the modelling process. However, in many cases, this may also result in poorer model quality for a given target or alignment method. There is a need for modelling protocols that can both consistently and significantly improve 3D models and provide an indication of when models might not benefit from the use of multiple target-template alignments. Here, we investigate the use of both global and local model quality prediction scores produced by ModFOLDclust2, to improve the selection of target-template alignments for the construction of multiple-template models. Additionally, we evaluate clustering the resulting population of multi- and single-template models for the improvement of our IntFOLD-TS tertiary structure prediction method. Results: We find that using accurate local model quality scores to guide alignment selection is the most consistent way to significantly improve models for each of the sequence to structure alignment methods tested. In addition, using accurate global model quality for re-ranking alignments, prior to selection, further improves the majority of multi-template modelling methods tested. Furthermore, subsequent clustering of the resulting population of multiple-template models significantly improves the quality of selected models compared with the previous version of our tertiary structure prediction method, IntFOLD-TS.

Towards data warehousing and mining of protein unfolding simulation data

OBJECTIVES: The prediction of protein structure and the precise understanding of protein folding and unfolding processes remains one of the greatest challenges in structural biology and bioinformatics. Computer simulations based on molecular dynamics (MD) are at the forefront of the effort to gain a deeper understanding of these complex processes. Currently, these MD simulations are usually on the order of tens of nanoseconds, generate a large amount of conformational data and are computationally expensive. More and more groups run such simulations and generate a myriad of data, which raises new challenges in managing and analyzing these data. Because the vast range of proteins researchers want to study and simulate, the computational effort needed to generate data, the large data volumes involved, and the different types of analyses scientists need to perform, it is desirable to provide a public repository allowing researchers to pool and share protein unfolding data. METHODS: To adequately organize, manage, and analyze the data generated by unfolding simulation studies, we designed a data warehouse system that is embedded in a grid environment to facilitate the seamless sharing of available computer resources and thus enable many groups to share complex molecular dynamics simulations on a more regular basis. RESULTS: To gain insight into the conformational fluctuations and stability of the monomeric forms of the amyloidogenic protein transthyretin (TTR), molecular dynamics unfolding simulations of the monomer of human TTR have been conducted. Trajectory data and meta-data of the wild-type (WT) protein and the highly amyloidogenic variant L55P-TTR represent the test case for the data warehouse. CONCLUSIONS: Web and grid services, especially pre-defined data mining services that can run on or 'near' the data repository of the data warehouse, are likely to play a pivotal role in the analysis of molecular dynamics unfolding data.

Refolded dengue virus type 2 NS1 protein expressed in Escherichia coli preserves structural and immunological properties of the native protein

The dengue virus NS1 protein has been shown to be a protective antigen under different experimental conditions but the recombinant protein produced in bacterial expression systems is usually not soluble and loses structural and immunological features of the native viral protein In the present study, experimental conditions leading to purification and refolding of the recombinant dengue virus type 2 (DENV-2) NS1 protein expressed in Escherichia coil are described The refolded recombinant protein was recovered as heat-stable soluble dimers with preserved structural features, as demonstrated by spectroscopic methods In addition, antibodies against epitopes of the NS1 protein expressed in eukaryotic cells recognized the refolded protein expressed in E coli but not the denatured form or the same protein submitted to a different refolding condition Collectively, the results demonstrate that the recombinant NS1 protein preserved important conformation and antigenic determinants of the native virus protein and represents a valuable reagent either for the development of vaccines or for diagnostic methods. (C) 2010 Elsevier B V All rights reserved

Structure, Dynamics, and RNA Interaction Analysis of the Human SBDS Protein

Shwachman-Bodian-Diamond syndrome is an autosomal recessive genetic syndrome with pleiotropic phenotypes, including pancreatic deficiencies, bone marrow dysfunctions with increased risk of myelodysplasia or leukemia, and skeletal abnormalities. This syndrome has been associated with mutations in the SBDS gene, which encodes a conserved protein showing orthologs in Archaea and eukaryotes. The Shwachman-Bodian-Diamond syndrome pleiotropic phenotypes may be an indication of different cell type requirements for a fully functional SBDS protein. RNA-binding activity has been predicted for archaeal and yeast SBDS orthologs, with the latter also being implicated in ribosome biogenesis. However, full-length SBDS orthologs function in a species-specific manner, indicating that the knowledge obtained from model systems may be of limited use in understanding major unresolved issues regarding SBDS function, namely, the effect of mutations in human SBDS on its biochemical function and the specificity of RNA interaction. We determined the solution structure and backbone dynamics of the human SBDS protein and describe its RNA binding site using NMR spectroscopy. Similarly to the crystal structures of Archaea, the overall structure of human SBDS comprises three well-folded domains. However, significant conformational exchange was observed in NMR dynamics experiments for the flexible linker between the N-terminal domain and the central domain, and these experiments also reflect the relative motions of the domains. RNA titrations monitored by heteronuclear correlation experiments and chemical shift mapping analysis identified a classic RNA binding site at the N-terminal FYSH (fungal, Yhr087wp, Shwachman) domain that concentrates most of the mutations described for the human SBDS. (C) 2010 Elsevier Ltd. All rights reserved.

Molecular modeling databases: A new way in the search of protein targets for drug development

DBMODELING is a relational database of annotated comparative protein structure models and their metabolic, pathway characterization. It is focused on enzymes identified in the genomes of Mycobacterium tuberculosis and Xylella fastidiosa. The main goal of the present database is to provide structural models to be used in docking simulations and drug design. However, since the accuracy of structural models is highly dependent on sequence identity between template and target, it is necessary to make clear to the user that only models which show high structural quality should be used in such efforts. Molecular modeling of these genomes generated a database, in which all structural models were built using alignments presenting more than 30% of sequence identity, generating models with medium and high accuracy. All models in the database are publicly accessible at http://www.biocristalografia.df.ibilce.unesp.br/tools. DBMODELING user interface provides users friendly menus, so that all information can be printed in one stop from any web browser. Furthermore, DBMODELING also provides a docking interface, which allows the user to carry out geometric docking simulation, against the molecular models available in the database. There are three other important homology model databases: MODBASE, SWISSMODEL, and GTOP. The main applications of these databases are described in the present article. © 2007 Bentham Science Publishers Ltd.

Relationship between global structural parameters and Enzyme Commission hierarchy: Implications for function prediction

In protein databases there is a substantial number of proteins structurally determined but without function annotation. Understanding the relationship between function and structure can be useful to predict function on a large scale. We have analyzed the similarities in global physicochemical parameters for a set of enzymes which were classified according to the four Enzyme Commission (EC) hierarchical levels. Using relevance theory we introduced a distance between proteins in the space of physicochemical characteristics. This was done by minimizing a cost function of the metric tensor built to reflect the EC classification system. Using an unsupervised clustering method on a set of 1025 enzymes, we obtained no relevant clustering formation compatible with EC classification. The distance distributions between enzymes from the same EC group and from different EC groups were compared by histograms. Such analysis was also performed using sequence alignment similarity as a distance. Our results suggest that global structure parameters are not sufficient to segregate enzymes according to EC hierarchy. This indicates that features essential for function are rather local than global. Consequently, methods for predicting function based on global attributes should not obtain high accuracy in main EC classes prediction without relying on similarities between enzymes from training and validation datasets. Furthermore, these results are consistent with a substantial number of studies suggesting that function evolves fundamentally by recruitment, i.e., a same protein motif or fold can be used to perform different enzymatic functions and a few specific amino acids (AAs) are actually responsible for enzyme activity. These essential amino acids should belong to active sites and an effective method for predicting function should be able to recognize them. (C) 2012 Elsevier Ltd. All rights reserved.

Characterization of novel StAR (steroidogenic acute regulatory protein) mutations causing non-classic lipoid adrenal hyperplasia

Context Steroidogenic acute regulatory protein (StAR) is crucial for transport of cholesterol to mitochondria where biosynthesis of steroids is initiated. Loss of StAR function causes lipoid congenital adrenal hyperplasia (LCAH). Objective StAR gene mutations causing partial loss of function manifest atypical and may be mistaken as familial glucocorticoid deficiency. Only a few mutations have been reported. Design To report clinical, biochemical, genetic, protein structure and functional data on two novel StAR mutations, and to compare them with published literature. Setting Collaboration between the University Children's Hospital Bern, Switzerland, and the CIBERER, Hospital Vall d'Hebron, Autonomous University, Barcelona, Spain. Patients Two subjects of a non-consanguineous Caucasian family were studied. The 46,XX phenotypic normal female was diagnosed with adrenal insufficiency at the age of 10 months, had normal pubertal development and still has no signs of hypergonodatropic hypogonadism at 32 years of age. Her 46,XY brother was born with normal male external genitalia and was diagnosed with adrenal insufficiency at 14 months. Puberty was normal and no signs of hypergonadotropic hypogonadism are present at 29 years of age. Results StAR gene analysis revealed two novel compound heterozygote mutations T44HfsX3 and G221S. T44HfsX3 is a loss-of-function StAR mutation. G221S retains partial activity (~30%) and is therefore responsible for a milder, non-classic phenotype. G221S is located in the cholesterol binding pocket and seems to alter binding/release of cholesterol. Conclusions StAR mutations located in the cholesterol binding pocket (V187M, R188C, R192C, G221D/S) seem to cause non-classic lipoid CAH. Accuracy of genotype-phenotype prediction by in vitro testing may vary with the assays employed.

X-ray structure of the human hyperplastic discs protein: An ortholog of the C-terminal domain of poly(A)-binding protein

The poly(A)-binding protein (PABP) recognizes the 3′ mRNA poly(A) tail and plays an essential role in eukaryotic translation initiation and mRNA stabilization/degradation. PABP is a modular protein, with four N-terminal RNA-binding domains and an extensive C terminus. The C-terminal region of PABP is essential for normal growth in yeast and has been implicated in mediating PABP homo-oligomerization and protein–protein interactions. A small, proteolytically stable, highly conserved domain has been identified within this C-terminal segment. Remarkably, this domain is also present in the hyperplastic discs protein (HYD) family of ubiquitin ligases. To better understand the function of this conserved region, an x-ray structure of the PABP-like segment of the human HYD protein has been determined at 1.04-Å resolution. The conserved domain adopts a novel fold resembling a right-handed supercoil of four α-helices. Sequence profile searches and comparative protein structure modeling identified a small ORF from the Arabidopsis thaliana genome that encodes a structurally similar but distantly related PABP/HYD domain. Phylogenetic analysis of the experimentally determined (HYD) and homology modeled (PABP) protein surfaces revealed a conserved feature that may be responsible for binding to a PABP interacting protein, Paip1, and other shared interaction partners.

Top-down free-energy minimization on protein potential energy landscapes

The hierarchical properties of potential energy landscapes have been used to gain insight into thermodynamic and kinetic properties of protein ensembles. It also may be possible to use them to direct computational searches for thermodynamically stable macroscopic states, i.e., computational protein folding. To this end, we have developed a top-down search procedure in which conformation space is recursively dissected according to the intrinsic hierarchical structure of a landscape's effective-energy barriers. This procedure generates an inverted tree similar to the disconnectivity graphs generated by local minima-clustering methods, but it fundamentally differs in the manner in which the portion of the tree that is to be computationally explored is selected. A key ingredient is a branch-selection algorithm that takes advantage of statistically predictive properties of the landscape to guide searches down the tree branches that are most likely to lead to the physically relevant macroscopic states. Using the computational folding of a β-hairpin-forming peptide as an example, we show that such predictive properties indeed exist and can be used for structure prediction by free-energy global minimization.

Fourier Transform Infrared Microspectroscopy Detects Changes in Protein Secondary Structure Associated with Desiccation Tolerance in Developing Maize Embryos1

Isolated immature maize (Zea mays L.) embryos have been shown to acquire tolerance to rapid drying between 22 and 25 d after pollination (DAP) and to slow drying from 18 DAP onward. To investigate adaptations in protein profile in association with the acquisition of desiccation tolerance in isolated, immature maize embryos, we applied in situ Fourier transform infrared microspectroscopy. In fresh, viable, 20- and 25-DAP embryo axes, the shapes of the different amide-I bands were identical, and this was maintained after flash drying. On rapid drying, the 20-DAP axes had a reduced relative proportion of α-helical protein structure and lost viability. Rapidly dried 25-DAP embryos germinated (74%) and had a protein profile similar to the fresh control axes. On slow drying, the α-helical contribution in both the 20- and 25-DAP embryo axes increased compared with that in the fresh control axes, and survival of desiccation was high. The protein profile in dry, mature axes resembled that after slow drying of the immature axes. Rapid drying resulted in an almost complete loss of membrane integrity in the 20-DAP embryo axes and much less so in the 25-DAP axes. After slow drying, low plasma membrane permeability ensued in both the 20- and 25-DAP axes. We conclude that slow drying of excised, immature embryos leads to an increased proportion of α-helical protein structures in their axes, which coincides with additional tolerance of desiccation stress.

West nile virus core protein: Tetramer structure and ribbon formation

We have determined the crystal structure of the core (C) protein from the Kunjin subtype of West Nile virus (WNV), closely related to the NY99 strain of WNV, currently a major health threat in the U.S. WNV is a member of the Flaviviridae family of enveloped RNA viruses that contains many important human pathogens. The C protein is associated with the RNA genome and forms the internal core which is surrounded by the envelope in the virion. The C protein structure contains four a. helices and forms dimers that are organized into tetramers. The tetramers form extended filamentous ribbons resembling the stacked alpha helices seen in HEAT protein structures.

Wurst: a protein threading server with a structural scoring function, sequence profiles and optimized substitution matrices

Wurst is a protein threading program with an emphasis on high quality sequence to structure alignments (http://www.zbh.uni-hamburg.de/wurst). Submitted sequences are aligned to each of about 3000 templates with a conventional dynamic programming algorithm, but using a score function with sophisticated structure and sequence terms. The structure terms are a log-odds probability of sequence to structure fragment compatibility, obtained from a Bayesian classification procedure. A simplex optimization was used to optimize the sequence-based terms for the goal of alignment and model quality and to balance the sequence and structural contributions against each other. Both sequence and structural terms operate with sequence profiles.

The 2.1 angstrom crystal structure of the far-red fluorescent protein HcRed: Inherent conformational flexibility of the chromophore

We have determined the crystal structure of HcRed, a far-red fluorescent protein isolated from Heteractis crispa, to 2.1 resolution. HcRed was observed to form a dimer, in contrast to the monomeric form of green fluorescent protein (GFP) or the tetrameric forms of the GFP-like proteins (eqFP611, Rtms5 and DsRed). Unlike the well-defined chromophore conformation observed in GFP and the GFP-like proteins, the HcRed chromophore was observed to be considerably mobile. Within the HcRed structure, the cyclic tripeptide chromophore, Glu64-Tyr65-Gly66, was observed to adopt both a cis coplanar and a tran. non-coplanar conformation. As a result of these two con formations, the hydroxyphenyl moiety of the chromophore makes distinct interactions within the interior of the b-can. These data together with a quantum chemical model of the chromophore, suggest the cis coplanar conformation to be consistent with the fluorescent properties of HcRed, and the trans non-coplanar conformation to be consistent with non-fluorescent properties of hcCP, the chromoprotein parent of HcRed. Moreover, within the GFP-like family, it appears that where conformational freedom is permissible then flexibility in the chromophore conformation is possible. 2005 Elsevier Ltd. All rights reserved.

Predicting the solvent accessibility of transmembrane residues from protein sequence

In this study, we propose a novel method to predict the solvent accessible surface areas of transmembrane residues. For both transmembrane alpha-helix and beta-barrel residues, the correlation coefficients between the predicted and observed accessible surface areas are around 0.65. On the basis of predicted accessible surface areas, residues exposed to the lipid environment or buried inside a protein can be identified by using certain cutoff thresholds. We have extensively examined our approach based on different definitions of accessible surface areas and a variety of sets of control parameters. Given that experimentally determining the structures of membrane proteins is very difficult and membrane proteins are actually abundant in nature, our approach is useful for theoretically modeling membrane protein tertiary structures, particularly for modeling the assembly of transmembrane domains. This approach can be used to annotate the membrane proteins in proteomes to provide extra structural and functional information.

The Effect Of Celastrol, A Triterpene With Antitumorigenic Activity, On Conformational And Functional Aspects Of The Human 90kda Heat Shock Protein Hsp90α, A Chaperone Implicated In The Stabilization Of The Tumor Phenotype.

Hsp90 is a molecular chaperone essential for cell viability in eukaryotes that is associated with the maturation of proteins involved in important cell functions and implicated in the stabilization of the tumor phenotype of various cancers, making this chaperone a notably interesting therapeutic target. Celastrol is a plant-derived pentacyclic triterpenoid compound with potent antioxidant, anti-inflammatory and anticancer activities; however, celastrol's action mode is still elusive. In this work, we investigated the effect of celastrol on the conformational and functional aspects of Hsp90α. Interestingly, celastrol appeared to target Hsp90α directly as the compound induced the oligomerization of the chaperone via the C-terminal domain as demonstrated by experiments using a deletion mutant. The nature of the oligomers was investigated by biophysical tools demonstrating that a two-fold excess of celastrol induced the formation of a decameric Hsp90α bound throughout the C-terminal domain. When bound, celastrol destabilized the C-terminal domain. Surprisingly, standard chaperone functional investigations demonstrated that neither the in vitro chaperone activity of protecting against aggregation nor the ability to bind a TPR co-chaperone, which binds to the C-terminus of Hsp90α, were affected by celastrol. Celastrol interferes with specific biological functions of Hsp90α. Our results suggest a model in which celastrol binds directly to the C-terminal domain of Hsp90α causing oligomerization. However, the ability to protect against protein aggregation (supported by our results) and to bind to TPR co-chaperones are not affected by celastrol. Therefore celastrol may act primarily by inducing specific oligomerization that affects some, but not all, of the functions of Hsp90α. To the best of our knowledge, this study is the first work to use multiple probes to investigate the effect that celastrol has on the stability and oligomerization of Hsp90α and on the binding of this chaperone to Tom70. This work provides a novel mechanism by which celastrol binds Hsp90α.