Gene ontology improves template selection in comparative protein docking

Structural characterization of protein-protein interactions is essential for our ability to study life processes at the molecular level. Computational modeling of protein complexes (protein docking) is important as the source of their structure and as a way to understand the principles of protein interaction. Rapidly evolving comparative docking approaches utilize target/template similarity metrics, which are often based on the protein structure. Although the structural similarity, generally, yields good performance, other characteristics of the interacting proteins (eg, function, biological process, and localization) may improve the prediction quality, especially in the case of weak target/template structural similarity. For the ranking of a pool of models for each target, we tested scoring functions that quantify similarity of Gene Ontology (GO) terms assigned to target and template proteins in three ontology domains—biological process, molecular function, and cellular component (GO-score). The scoring functions were tested in docking of bound, unbound, and modeled proteins. The results indicate that the combined structural and GO-terms functions improve the scoring, especially in the twilight zone of structural similarity, typical for protein models of limited accuracy.     

Detection of unrealistic molecular environments in protein structures based on expected electron densities

Understanding the relationship between protein structure and biological function is a central theme in structural biology. Advances are severely hampered by errors in experimentally determined protein structures. Detection and correction of such errors is therefore of utmost importance. Electron densities in molecular structures obey certain rules which depend on the molecular environment. Here we present and discuss a new approach that relates electron densities computed from a structural model to densities expected from prior observations on identical or closely related molecular environments. Strong deviations of computed from expected densities reveal unrealistic molecular structures. Most importantly, structure analysis and error detection are independent of experimental data and hence may be applied to any structural model. The comparison to state-of-the-art methods reveals that our approach is able to identify errors that formerly remained undetected. The new technique, called RefDens, is accessible as a public web service at .     

Disulfide bonds convert small heat shock protein Hsp16.3 from a chaperone to a non-chaperone: implications for the evolution of cysteine in molecular chaperones

Molecular chaperones mainly function in assisting newly synthesized polypeptide folding and protect non-native proteins from aggregation, with known structural features such as the ability of spontaneous folding/refolding and high conformational flexibility. In this report, we verified the assumption that the lack of disulfide bonds in molecular chaperones is a prerequisite for such unique structural features. Using small heat shock protein (one sub-class of chaperones) Hsp16.3 as a model system, our results show the following: (1) Cysteine-free Hsp16.3 wild type protein can efficiently exhibit chaperone activity and spontaneously refold/reassemble with high conformational flexibility. (2) Whereas Hsp16.3 G89C mutant with inter-subunit disulfide bonds formed seems to lose the nature of chaperone proteins, i.e., under stress conditions, it neither acts as molecular chaperone nor spontaneously refolds/reassembles. Structural analysis indicated that the mutant exists as an unstable molten globule-like state, which incorrectly exposes hydrophobic surfaces and irreversibly tends to form aggregates that can be suppressed by the other molecular chaperone (alpha-crystallin). By contrast, reduction of disulfide bond in the Hsp16.3 G89C mutant can significantly recover its character as a molecular chaperone. In light of these results, we propose that disulfide bonds could severely disturb the structure/function of molecular chaperones like Hsp16.3. Our results might not only provide insights into understanding the structural basis of chaperone upon binding substrates, but also explain the observation that the occurrence of cysteine in molecular chaperones is much lower than that in other protein families, subsequently being helpful to understand the evolution of protein family.     

The role of protein structure in genomics

The genome projects produce an enormous amount of sequence data that needs to be annotated in terms of molecular structure and biological function. These tasks have triggered additional initiatives like structural genomics. The intention is to determine as many protein structures as possible, in the most efficient way, and to exploit the solved structures for the assignment of biological function to hypothetical proteins. We discuss the impact of these developments on protein classification, gene function prediction, and protein structure prediction.     

Adaptability in protein structures: structural dynamics and implications in ligand design

The basic framework of understanding the mechanisms of protein functions is achieved from the knowledge of their structures which can model the molecular recognition. Recent advancement in the structural biology has revealed that in spite of the availability of the structural data, it is nontrivial to predict the mechanism of the molecular recognition which progresses via situation-dependent structural adaptation. The mutual selectivity of protein–protein and protein–ligand interactions often depends on the modulations of conformations empowered by their inherent flexibility, which in turn regulates the function. The mechanism of a protein’s function, which used to be explained by the ideas of ‘lock and key’ has evolved today as the concept of ‘induced fit’ as well as the ‘population shift’ models. It is felt that the ‘dynamics’ is an essential feature to take into account for understanding the mechanism of protein’s function. The design principles of therapeutic molecules suffer from the problems of plasticity of the receptors whose binding conformations are accurately not predictable from the prior knowledge of a template structure. On the other hand, flexibility of the receptors provides the opportunity to improve the binding affinity of a ligand by suitable substitution that will maximize the binding by modulating the receptors surface. In this paper, we discuss with example how the protein’s flexibility is correlated with its functions in various systems, revealing the importance of its understanding and for making applications. We also highlight the methodological challenges to investigate it computationally and to account for the flexible nature of the molecules in drug design.     

Force spectroscopy with single bio-molecules

For many biological molecules, force is an important functional and structural parameter. With the rapidly growing knowledge about the relation between structure, function, and force, single-molecule force spectroscopy has become a versatile analytical tool for the structural and functional investigation of single bio-molecules in their native environments. Within the past year, detailed insights into binding potentials of receptor ligand pairs, protein folding pathways, molecular motors, DNA mechanics and the functioning of DNA-binding agents (such as proteins and drugs), as well as the function of molecular motors, have been obtained.     

Molecular dynamics simulations and oxidation rates of methionine residues of granulocyte colony-stimulating factor at different pH values

To understand the connection between the conformation of a protein molecule and the oxidation of its methionine residues, we measured the rates of oxidation of methionine residues by H(2)O(2) in granulocyte colony-stimulating factor (G-CSF) as a function of pH and also studied the structural properties of this protein as a function of pH via molecular dynamics simulations. We found that each of the four methionine groups in G-CSF have significant and different rates of oxidation as a function of pH. Moreover, Met(1), in the unstructured N-terminal region, has a rate of oxidation as low as half that of free methionine. The structural properties of G-CSF as a function of pH are evaluated in terms of properties such as hydrogen bonding, deviations from X-ray structure, helical/helical packing, and the atomic covariance fluctuation matrix of alpha-carbons. We found that dynamics (structural fluctuations) are essential in explaining oxidation and that a static picture, such as that resulting from X-ray data, fails in this regard. Moreover, the simulation results also indicate that the solvent-accessible area, traditionally used to measure solvent accessibility of a protein site, of the sulfur atom of methionine residues does not correlate well with the rate of oxidation. Instead, we identified a structural property, average two-shell water coordination number, that correlates well with measured oxidation rates.     

基于结构的单绕蛋白聚类图构建与分析

蛋白质分子进化规律研究是分子进化研究的重点,对揭示生命起源与进化机制有重要意义。本文对已知空间结构及物种信息的单绕蛋白,利用结构比对信息,构建了不同层次单绕样本系统聚类图。分析发现:功能相似蛋白存在明显聚集现象,同一超家族样本基本聚在一个大支中,同一家族样本集中在所属超家族下的小支中,功能约束下单绕样本聚类图与物种进化图有较好对应关系。结果表明:单绕蛋白的结构演化反映了蛋白质功能的约束,特定功能单绕样本的结构差异具有种属特异性,结构演化包含了物种进化信息。     

A biological cosmos of parallel universes: Does protein structural plasticity facilitate evolution?

While Darwin pictured organismal evolution as "descent with modification" more than 150 years ago, a detailed reconstruction of the basic evolutionary transitions at the molecular level is only emerging now. In particular, the evolution of today's protein structures and their concurrent functions has remained largely mysterious, as the destruction of these structures by mutation seems far easier than their construction. While the accumulation of genomic and structural data has indicated that proteins are related via common ancestors, naturally occurring protein structures are often considered to be evolutionarily robust, thus leaving open the question of how protein structures can be remodelled while selective pressure forces them to function. New information on the proteome, however, increasingly explains the nature of local and global conformational diversity in protein evolution, which allows the acquisition of novel functions via molecular transition forms containing ancestral and novel structures in dynamic equilibrium. Such structural plasticity may permit the evolution of new protein folds and help account for both the origins of new biological functions and the nature of molecular defects.     

Solid-state NMR characterization of conformational plasticity within the transmembrane domain of the influenza A M2 proton channel

Membrane protein function within the membrane interstices is achieved by mechanisms that are not typically available to water-soluble proteins. The whole balance of molecular interactions that stabilize three-dimensional structure in the membrane environment is different from that in an aqueous environment. As a result interhelical interactions are often dominated by non-specific van der Waals interactions permitting dynamics and conformational heterogeneity in these interfaces. Here, solid-state NMR data of the transmembrane domain of the M2 protein from influenza A virus are used to exemplify such conformational plasticity in a tetrameric helical bundle. Such data lead to very high resolution structural restraints that can identify both subtle and substantial structural differences associated with various states of the protein. Spectra from samples using two different preparation protocols, samples prepared in the presence and absence of amantadine, and spectra as a function of pH are used to illustrate conformational plasticity.     

Solid-state NMR characterization of conformational plasticity within the transmembrane domain of the influenza A M2 proton channel

Membrane protein function within the membrane interstices is achieved by mechanisms that are not typically available to water-soluble proteins. The whole balance of molecular interactions that stabilize three-dimensional structure in the membrane environment is different from that in an aqueous environment. As a result interhelical interactions are often dominated by non-specific van der Waals interactions permitting dynamics and conformational heterogeneity in these interfaces. Here, solid-state NMR data of the transmembrane domain of the M2 protein from influenza A virus are used to exemplify such conformational plasticity in a tetrameric helical bundle. Such data lead to very high resolution structural restraints that can identify both subtle and substantial structural differences associated with various states of the protein. Spectra from samples using two different preparation protocols, samples prepared in the presence and absence of amantadine, and spectra as a function of pH are used to illustrate conformational plasticity.     

Evidence of Functional Protein Dynamics from X-Ray Crystallographic Ensembles

It is widely recognized that representing a protein as a single static conformation is inadequate to describe the dynamics essential to the performance of its biological function. We contrast the amino acid displacements below and above the protein dynamical transition temperature, T_D∼215K, of hen egg white lysozyme using X-ray crystallography ensembles that are analyzed by molecular dynamics simulations as a function of temperature. We show that measuring structural variations across an ensemble of X-ray derived models captures the activation of conformational states that are of functional importance just above T_D, and they remain virtually identical to structural motions measured at 300K. Our results highlight the ability to observe functional structural variations across an ensemble of X-ray crystallographic data, and that residue fluctuations measured in MD simulations at room temperature are in quantitative agreement with the experimental observable.     

A Proline-Tryptophan Turn in the Intrinsically Disordered Domain 2 of NS5A Protein Is Essential for Hepatitis C Virus RNA Replication

Hepatitis C virus (HCV) nonstructural protein 5A (NS5A) and its interaction with the human chaperone cyclophilin A are both targets for highly potent and promising antiviral drugs that are in the late stages of clinical development. Despite its high interest in regards to the development of drugs to counteract the worldwide HCV burden, NS5A is still an enigmatic multifunctional protein poorly characterized at the molecular level. NS5A is required for HCV RNA replication and is involved in viral particle formation and regulation of host pathways. Thus far, no enzymatic activity or precise molecular function has been ascribed to NS5A that is composed of a highly structured domain 1 (D1), as well as two intrinsically disordered domains 2 (D2) and 3 (D3), representing half of the protein. Here, we identify a short structural motif in the disordered NS5A-D2 and report its NMR structure. We show that this structural motif, a minimal Pro³¹⁴–Trp³¹⁶ turn, is essential for HCV RNA replication, and its disruption alters the subcellular distribution of NS5A. We demonstrate that this Pro-Trp turn is required for proper interaction with the host cyclophilin A and influences its peptidyl-prolyl cis/trans isomerase activity on residue Pro³¹⁴ of NS5A-D2. This work provides a molecular basis for further understanding of the function of the intrinsically disordered domain 2 of HCV NS5A protein. In addition, our work highlights how very small structural motifs present in intrinsically disordered proteins can exert a specific function.     

Crystal structure of a hypothetical protein,TM841 of Thermotoga maritima,reveals its function as a fatty acid-binding protein

We determined the three-dimensional (3D) crystal structure of protein TM841, a protein product from a hypothetical open-reading frame in the genome of the hyperthermophile bacterium Thermotoga maritima, to 2.0 A resolution. The protein belongs to a large protein family, DegV or COG1307 of unknown function. The 35 kDa protein consists of two separate domains, with low-level structural resemblance to domains from other proteins with known 3D structures. These structural homologies, however, provided no clues for the function of TM841. But the electron density maps revealed clear density for a bound fatty-acid molecule in a pocket between the two protein domains. The structure indicates that TM841 has the molecular function of fatty-acid binding and may play a role in the cellular functions of fatty acid transport or metabolism.     

Insights from molecular modeling and dynamics simulation of pathogen resistance (R) protein from brinjal

Resistance (R) protein recognizes molecular signature of pathogen infection and activates downstream hypersensitive response signalling in plants. R protein works as a molecular switch for pathogen defence signalling and represent one of the largest plant gene family. Hence, understanding molecular structure and function of R proteins has been of paramount importance for plant biologists. The present study is aimed at predicting structure of R proteins signalling domains (CC-NBS) by creating a homology model, refining and optimising the model by molecular dynamics simulation and comparing ADP and ATP binding. Based on sequence similarity with proteins of known structures, CC-NBS domains were initially modelled using CED- 4 (cell death abnormality protein) and APAF-1 (apoptotic protease activating factor) as multiple templates. The final CC-NBS structural model was built and optimized by molecular dynamic simulation for 5 nanoseconds (ns). Docking of ADP and ATP at active site shows that both ligand bind specifically with same residues and with minor difference (1 Kcal/mol) in binding energy. Sharing of binding site by ADP and ATP and low difference in their binding site makes CC-NBS suitable for working as molecular switch. Furthermore, structural superimposition elucidate that CC-NBS and CARD (caspase recruitment domains) domain of CED-4 have low RMSD value of 0.9 A° Availability of 3D structural model for both CC and NBS domains will . help in getting deeper insight in these pathogen defence genes.     

Computer-based screening of functional conformers of proteins

A long-standing goal in biology is to establish the link between function, structure, and dynamics of proteins. Considering that protein function at the molecular level is understood by the ability of proteins to bind to other molecules, the limited structural data of proteins in association with other bio-molecules represents a major hurdle to understanding protein function at the structural level. Recent reports show that protein function can be linked to protein structure and dynamics through network centrality analysis, suggesting that the structures of proteins bound to natural ligands may be inferred computationally. In the present work, a new method is described to discriminate protein conformations relevant to the specific recognition of a ligand. The method relies on a scoring system that matches critical residues with central residues in different structures of a given protein. Central residues are the most traversed residues with the same frequency in networks derived from protein structures. We tested our method in a set of 24 different proteins and more than 260,000 structures of these in the absence of a ligand or bound to it. To illustrate the usefulness of our method in the study of the structure/dynamics/function relationship of proteins, we analyzed mutants of the yeast TATA-binding protein with impaired DNA binding. Our results indicate that critical residues for an interaction are preferentially found as central residues of protein structures in complex with a ligand. Thus, our scoring system effectively distinguishes protein conformations relevant to the function of interest.     

Interaction-site prediction for protein complexes: a critical assessment

MOTIVATION: Proteins function through interactions with other proteins and biomolecules. Protein-protein interfaces hold key information toward molecular understanding of protein function. In the past few years, there have been intensive efforts in developing methods for predicting protein interface residues. A review that presents the current status of interface prediction and an overview of its applications and project future developments is in order. SUMMARY: Interface prediction methods rely on a wide range of sequence, structural and physical attributes that distinguish interface residues from non-interface surface residues. The input data are manipulated into either a numerical value or a probability representing the potential for a residue to be inside a protein interface. Predictions are now satisfactory for complex-forming proteins that are well represented in the Protein Data Bank, but less so for under-represented ones. Future developments will be directed at tackling problems such as building structural models for multi-component structural complexes.     

Molecular cloning and characterization of spiggin. An androgen-regulated extraorganismal adhesive with structural similarities to von Willebrand Factor-related proteins

One of the most definitive examples of a vertebrate extraorganismal structural protein can be found in three-spined sticklebacks (Gasterosteus aculeatus). In the breeding male the kidney hypertrophies and synthesizes an adhesive protein called "spiggin," which is secreted into the urinary bladder from where it is employed as a structural thread for nest building. This paper describes the first molecular characterization of spiggin and demonstrates that this adhesive is a protein complex assembled from a potential of three distinct subunits (alpha, beta, and gamma). These subunits arise by alternative splicing, and 11-ketoandrogens induce their expression in stickleback kidneys. Analysis of the predicted amino acid sequence of each subunit reveals a modular organization whose structural elements display a similarity to the multimerization domains found within von Willebrand Factor-related proteins. These results implicate that spiggin utilizes a conserved multimerization mechanism for the formation of a viscous agglutinate from its constituent subunits in the urinary bladders of male sticklebacks. This novel extraorganismal structural protein is therefore ideally suited to its function as an adhesive thread.     

The mammalian autophagy initiator complex contains 2 HORMA domain proteins

ATG101 is an essential component of the ULK complex responsible for initiating cellular autophagy in mammalian cells; its 3-dimensional structure and molecular function, however, are currently unclear. Here we present the X-ray structure of human ATG101. The protein displays an open HORMA domain fold. Both structural properties and biophysical evidence indicate that ATG101 is locked in this conformation, in contrast to the prototypical HORMA domain protein MAD2. Moreover, we discuss a potential mode of dimerization with ATG13 as a fundamental aspect of ATG101 function.