Comparative studies on ion-pair energetic,distribution among three domains of life: Archaea,eubacteria, and eukarya

Salt-bridges play a unique role in the structural and functional stability of proteins, especially under harsh environments. How these salt-bridges contribute to the overall thermodynamic stability of protein structure and function across different domains of life is elusive still date. To address the issue, statistical analyses on the energies of salt-bridges, involved in proteins' structure and function, are performed across three domains of life, that is, archaea, eubacteria, and eukarya. Results show that although the majority of salt-bridges are stable and conserved, yet the stability of archaeal proteins (∆∆G_net = −5.06 ± 3.8) is much more than that of eubacteria (∆∆G_net = −3.7 ± 2.9) and eukarya (∆∆G_net = −3.54 ± 3.1). Unlike earlier study with archaea, in eukarya and eubacteria, not all buried salt-bridge in our dataset are stable. Buried salt-bridges play surprising role in protein stability, whose variations are clearly observed among these domains. Greater desolvation penalty of buried salt-bridges is compensated by stable network of salt-bridges apart from equal contribution of bridge and background energy terms. On the basis proteins' secondary structure, topology, and evolution, our observation shows that salt-bridges when present closer to each other in sequence tend to form a greater number. Overall, our comparative study provides insight into the role of specific electrostatic interactions in proteins from different domains of life, which we hope, would be useful for protein engineering and bioinformatics study.     

Disordered domains and high surface charge confer hubs with the ability to interact with multiple proteins in interaction networks

We investigate the structural properties of hubs that enable them to interact with several partners in protein-protein interaction networks. We find that hubs have more observed and predicted disordered residues with fewer loops/coils, and more charged residues on the surface as compared to non-hubs. Smaller hubs have fewer disordered residues and more charged residues on the surface than larger hubs. We conclude that the global flexibility provided by disordered domains, and high surface charge are complementary factors that play a significant role in the binding ability of hubs.     

Computational protein design and large-scale assessment by I-TASSER structure assembly simulations

Protein design aims at designing new protein molecules of desired structure and functionality. One of the major obstacles to large-scale protein design are the extensive time and manpower requirements for experimental validation of designed sequences. Recent advances in protein structure prediction have provided potentials for an automated assessment of the designed sequences via folding simulations. We present a new protocol for protein design and validation. The sequence space is initially searched by Monte Carlo sampling guided by a public atomic potential, with candidate sequences selected by the clustering of sequence decoys. The designed sequences are then assessed by I-TASSER folding simulations, which generate full-length atomic structural models by the iterative assembly of threading fragments. The protocol is tested on 52 nonhomologous single-domain proteins, with an average sequence identity of 24% between the designed sequences and the native sequences. Despite this low sequence identity, three-dimensional models predicted for the first designed sequence have an RMSD of < 2 Å to the target structure in 62% of cases. This percentage increases to 77% if we consider the three-dimensional models from the top 10 designed sequences. Such a striking consistency between the target structure and the structural prediction from nonhomologous sequences, despite the fact that the design and folding algorithms adopt completely different force fields, indicates that the design algorithm captures the features essential to the global fold of the target. On average, the designed sequences have a free energy that is 0.39 kcal/(mol residue) lower than in the native sequences, potentially affording a greater stability to synthesized target folds.     

A Structural Portrait of the PDZ Domain Family

PDZ (PSD-95/Discs-large/ZO1) domains are interaction modules that typically bind to specific C-terminal sequences of partner proteins and assemble signaling complexes in multicellular organisms. We have analyzed the existing database of PDZ domain structures in the context of a specificity tree based on binding specificities defined by peptide-phage binding selections. We have identified 16 structures of PDZ domains in complex with high-affinity ligands and have elucidated four additional structures to assemble a structural database that covers most of the branches of the PDZ specificity tree. A detailed comparison of the structures reveals features that are responsible for the diverse specificities across the PDZ domain family. Specificity differences can be explained by differences in PDZ residues that are in contact with the peptide ligands, but these contacts involve both side-chain and main-chain interactions. Most PDZ domains bind peptides in a canonical conformation in which the ligand main chain adopts an extended β-strand conformation by interacting in an antiparallel fashion with a PDZ β-strand. However, a subset of PDZ domains bind peptides with a bent main-chain conformation and the specificities of these non-canonical domains could not be explained based on canonical structures. Our analysis provides a structural portrait of the PDZ domain family, which serves as a guide in understanding the structural basis for the diverse specificities across the family.     

Active-Site Monovalent Cations Revealed in a 1.55-Å-Resolution Hammerhead Ribozyme Structure

We have obtained a 1.55-Å crystal structure of a hammerhead ribozyme derived from Schistosoma mansoni under conditions that permit detailed observations of Na⁺ ion binding in the ribozyme's active site. At least two such Na⁺ ions are observed. The first Na⁺ ion binds to the N7 of G10.1 and the adjacent A9 phosphate in a manner identical with that previously observed for divalent cations. A second Na⁺ ion binds to the Hoogsteen face of G12, the general base in the hammerhead cleavage reaction, thereby potentially dissipating the negative charge of the catalytically active enolate form of the nucleotide base. A potential but more ambiguous third site bridges the A9 and scissile phosphates in a manner consistent with that of previous predictions. Hammerhead ribozymes have been observed to be active in the presence of high concentrations of monovalent cations, including Na⁺, but the mechanism by which monovalent cations substitute for divalent cations in hammerhead catalysis remains unclear. Our results enable us to suggest that Na⁺ directly and specifically substitutes for divalent cations in the hammerhead active site. The detailed geometry of the pre-catalytic active-site complex is also revealed with a new level of precision, thanks to the quality of the electron density maps obtained from what is currently the highest-resolution ribozyme structure in the Protein Data Bank.     

Travel depth, a new shape descriptor for macromolecules: application to ligand binding

Depth is a term frequently applied to the shape and surface of macromolecules, describing for example the grooves in DNA, the shape of an enzyme active site, or the binding site for a small molecule in a protein. Yet depth is a difficult property to define rigorously in a macromolecule, and few computational tools exist to quantify this notion, to visualize it, or analyze the results. We present our notion of travel depth, simply put the physical distance a solvent molecule would have to travel from a surface point to a suitably defined reference surface. To define the reference surface, we use the limiting form of the molecular surface with increasing probe size: the convex hull. We then present a fast, robust approximation algorithm to compute travel depth to every surface point. The travel depth is useful because it works for pockets of any size and complexity. It also works for two interesting special cases. First, it works on the grooves in DNA, which are unbounded in one direction. Second, it works on the case of tunnels, that is pockets that have no "bottom", but go through the entire macromolecule. Our algorithm makes it straightforward to quantify discussions of depth when analyzing structures. High-throughput analysis of macromolecule depth is also enabled by our algorithm. This is demonstrated by analyzing a database of protein-small molecule binding pockets, and the distribution of bound magnesium ions in RNA structures. These analyses show significant, but subtle effects of depth on ligand binding localization and strength.     

生物化学教学中实用素材的推荐

本文旨在与从事生物化学教学的教师分享几个实用素材,正在学习生物化学的学生亦可直接使用.素材包括多媒体类例如绪论导入歌曲"Biochemistry",ATP合酶工作原理的视频,RCSB PDB蛋白质结构数据库,RasMol分子图像观察软件等;实物类例如DNA双螺旋结构的纸质模型,tRNA一二三级结构的纸质模型等.侧重素材...     

Recovering the true targets of specific ligands by virtual screening of the protein data bank

The Protein Data Bank (PDB) has been processed to extract a screening protein library (sc-PDB) of 2148 entries. A knowledge-based detection algorithm has been applied to 18,000 PDB files to find regular expressions corresponding to either protein, ions, co-factors, solvent, or ligand atoms. The sc-PDB database comprises high-resolution X-ray structures of proteins for which (i) a well-defined active site exists, (ii) the bound-ligand is a small molecular weight molecule. The database has been screened by an inverse docking tool derived from the GOLD program to recover the known target of four unrelated ligands. Both the database and the inverse screening procedures are accurate enough to rank the true target of the four investigated ligands among the top 1% scorers, with 70-100 fold enrichment with respect to random screening. Applying the proposed screening procedure to a small-sized generic ligand was much less accurate suggesting that inverse screening shall be reserved to rather selective compounds.     

蛋白质结构与功能研究中的分子模拟技术

分子模拟技术为蛋白质的研究提供了一种崭新的手段,在理论上解决了结构预测和功能分析以及蛋白质工程实施方面所面临的难题。它在蛋白质的结构预测和模建工作中占有举足轻重的地位,实现了生物技术与计算机技术的完美结合。本文简要阐述了该技术的基本步骤和工作原理,并以目前应用最广的生物大分子领域的商品化分子模拟软件Accelrys公司基于Linux系统开发的InsightII为例,介绍了相关程序模块的功能和作用,同时结合该技术在蛋白质的结构预测和模建、结构与功能关系分析、分子设计等过程中的开发与应用,加以具体说明和展望。     

The role of RNA sequence and structure in RNA--protein interactions

We investigate the sequence and structural properties of RNA-protein interaction sites in 211 RNA-protein chain pairs, the largest set of RNA-protein complexes analyzed to date. Statistical analysis confirms and extends earlier analyses made on smaller data sets. There are 24.6% of hydrogen bonds between RNA and protein that are nucleobase specific, indicating the importance of both nucleobase-specific and -nonspecific interactions. While there is no significant difference between RNA base frequencies in protein-binding and non-binding regions, distinct preferences for RNA bases, RNA structural states, protein residues, and protein secondary structure emerge when nucleobase-specific and -nonspecific interactions are considered separately. Guanine nucleobase and unpaired RNA structural states are significantly preferred in nucleobase-specific interactions; however, nonspecific interactions disfavor guanine, while still favoring unpaired RNA structural states. The opposite preferences of nucleobase-specific and -nonspecific interactions for guanine may explain discrepancies between earlier studies with regard to base preferences in RNA-protein interaction regions. Preferences for amino acid residues differ significantly between nucleobase-specific and -nonspecific interactions, with nonspecific interactions showing the expected bias towards positively charged residues. Irregular protein structures are strongly favored in interactions with the protein backbone, whereas there is little preference for specific protein secondary structure in either nucleobase-specific interaction or -nonspecific interaction. Overall, this study shows strong preferences for both RNA bases and RNA structural states in protein-RNA interactions, indicating their mutual importance in protein recognition.     

Intelligent Access to Sequence and Structure Databases (IASSD) ? an interface for accessing information from major web databases

AvailabilityThe Jar file is freely available through e-mail from the corresponding author.     

A systematic analysis of the beta hairpin motif in the Protein Data Bank

The beta hairpin motif is a ubiquitous protein structural motif that can be found in molecules across the tree of life. This motif, which is also popular in synthetically designed proteins and peptides, is known for its stability and adaptability to broad functions. Here, we systematically probe all 49,000 unique beta hairpin substructures contained within the Protein Data Bank (PDB) to uncover key characteristics correlated with stable beta hairpin structure, including amino acid biases and enriched interstrand contacts. We find that position specific amino acid preferences, while seen throughout the beta hairpin structure, are most evident within the turn region, where they depend on subtle turn dynamics associated with turn length and secondary structure. We also establish a set of broad design principles, such as the inclusion of aspartic acid residues at a specific position and the careful consideration of desired secondary structure when selecting residues for the turn region, that can be applied to the generation of libraries encoding proteins or peptides containing beta hairpin structures.     

The limit of accuracy of protein modeling: influence of crystal packing on protein structure

The size of the protein database (PDB) makes it now feasible to arrive at statistical conclusions regarding structural effects of crystal packing. These effects are relevant for setting upper practical limits of accuracy on protein modeling. Proteins whose crystals have more than one molecule in the asymmetric unit or whose structures were determined at least twice by X-ray crystallography were paired and their differences analyzed. We demonstrate a clear influence of crystal environment on protein structure, including backbone conformations, hinge-like motions and side-chain conformations. The positions of surface water molecules tend to be variable in different crystal environments while those of ligands are not. Structures determined by independent groups vary more than structures determined by the same authors. The use of different refinement methods is a major source for this effect. Our pair-wise analysis derives a practical limit to the accuracy of protein modeling. For different crystal forms, the limit of accuracy (C(alpha), root-mean-square deviation (RMSD)) is approximately 0.8A for the entire protein, which includes approximately 0.3A due to crystal packing. For organized secondary elements, the upper limit of C(alpha) RMSD is 0.5-0.6A while for loops or protein surface it reaches 1.0A. Twenty percent of exposed side- chains exhibit different chi(1+2) conformations with approximately half of the effect also resulting from crystal packing. A web based tool for analysis and graphic presentation of surface areas of crystal contacts is available (http://ligin.weizmann.ac.il/cryco).     

Mapping of protein-protein interaction sites by the 'absence of interference' approach

Protein-protein interactions are critical to most biological processes, and locating protein-protein interfaces on protein structures is an important task in molecular biology. We developed a new experimental strategy called the ‘absence of interference’ approach to determine surface residues involved in protein-protein interaction of established yeast two-hybrid pairs of interacting proteins. One of the proteins is subjected to high-level randomization by error-prone PCR. The resulting library is selected by yeast two-hybrid system for interacting clones that are isolated and sequenced. The interaction region can be identified by an absence or depletion of mutations. For data analysis and presentation, we developed a Web interface that analyzes the mutational spectrum and displays the mutational frequency on the surface of the structure (or a structural model) of the randomized protein†. Additionally, this interface might be of use for the display of mutational distributions determined by other types of random mutagenesis experiments. We applied the approach to map the interface of the catalytic domain of the DNA methyltransferase Dnmt3a with its regulatory factor Dnmt3L. Dnmt3a was randomized with high mutational load. A total of 76 interacting clones were isolated and sequenced, and 648 mutations were identified. The mutational pattern allowed to identify a unique interaction region on the surface of Dnmt3a, which comprises about 500-600 Å². The results were confirmed by site-directed mutagenesis and structural analysis. The absence-of-interference approach will allow high-throughput mapping of protein interaction sites suitable for functional studies and protein docking.     

Facilities that make the PDB data collection more powerful

We describe a series of databases and tools that directly or indirectly support biomedical research on macromolecules, with focus on their applicability in protein structure bioinformatics research. DSSP, that determines secondary structures of proteins, has been updated to work well with extremely large structures in multiple formats. The PDBREPORT database that lists anomalies in protein structures has been remade to remove many small problems. These reports are now available as PDF‐formatted files with a computer‐readable summary. The VASE software has been added to analyze and visualize HSSP multiple sequence alignments for protein structures. The Lists collection of databases has been extended with a series of databases, most noticeably with a database that gives each protein structure a grade for usefulness in protein structure bioinformatics projects. The PDB‐REDO collection of reanalyzed and re‐refined protein structures that were solved by X‐ray crystallography has been improved by dealing better with sugar residues and with hydrogen bonds, and adding many missing surface loops. All academic software underlying these protein structure bioinformatics applications and databases are now publicly accessible, either directly from the authors or from the GitHub software repository.     

On the Structure of the Proton-Binding Site in the Fo Rotor of Chloroplast ATP Synthases

The recently reported crystal structures of the membrane-embedded proton-dependent c-ring rotors of a cyanobacterial F₁F_o ATP synthase and a chloroplast F₁F_o ATP synthase have provided new insights into the mechanism of this essential enzyme. While the overall features of these c-rings are similar, a discrepancy in the structure and hydrogen-bonding interaction network of the H⁺ sites suggests two distinct binding modes, potentially reflecting a mechanistic differentiation. Importantly, the conformation of the key glutamate side chain to which the proton binds is also altered. To investigate the nature of these differences, we use molecular dynamics simulations of both c-rings embedded in a phospholipid membrane. We observe that the structure of the c₁₅ ring from Spirulina platensis is unequivocally stable within the simulation time. By contrast, the proposed structure of the H⁺ site in the chloroplast c₁₄ ring changes rapidly and consistently into that reported for the c₁₅ ring, indicating that the latter represents a common binding mode. To assess this hypothesis, we have remodeled the c₁₄ ring by molecular replacement using the published structure factors. The resulting structure provides clear evidence in support of a common binding site conformation and is also considerably improved statistically. These findings, taken together with a sequence analysis of c-subunits in the ATP synthase family, indicate that the so-called proton-locked conformation observed in the c₁₅ ring may be a common characteristic not only of light-driven systems such as chloroplasts and cyanobacteria but also of a selection of other bacterial species.     

Homology‐based hydrogen bond information improves crystallographic structures in the PDB

The Protein Data Bank (PDB) is the global archive for structural information on macromolecules, and a popular resource for researchers, teachers, and students, amassing more than one million unique users each year. Crystallographic structure models in the PDB (more than 100,000 entries) are optimized against the crystal diffraction data and geometrical restraints. This process of crystallographic refinement typically ignored hydrogen bond (H‐bond) distances as a source of information. However, H‐bond restraints can improve structures at low resolution where diffraction data are limited. To improve low‐resolution structure refinement, we present methods for deriving H‐bond information either globally from well‐refined high‐resolution structures from the PDB‐REDO databank, or specifically from on‐the‐fly constructed sets of homologous high‐resolution structures. Refinement incorporating HOmology DErived Restraints (HODER), improves geometrical quality and the fit to the diffraction data for many low‐resolution structures. To make these improvements readily available to the general public, we applied our new algorithms to all crystallographic structures in the PDB: using massively parallel computing, we constructed a new instance of the PDB‐REDO databank ( https://pdb-redo.eu ). This resource is useful for researchers to gain insight on individual structures, on specific protein families (as we demonstrate with examples), and on general features of protein structure using data mining approaches on a uniformly treated dataset.     

Crystal structure of scytalidoglutamic peptidase with its first potent inhibitor provides insights into substrate specificity and catalysis

Scytalidoglutamic peptidase (SGP) from Scytalidium lignicolum is the founding member of the newly discovered\ family of peptidases, G1, so far found exclusively in fungi. The crystal structure of SGP revealed a previously undescribed fold for peptidases and a unique catalytic dyad of residues Gln53 and Glu136. Surprisingly, the beta-sandwich structure of SGP is strikingly similar to members of the carbohydrate-binding concanavalin A-like lectins/glucanases superfamily. By analogy with the active sites of aspartic peptidases, a mechanism employing nucleophillic attack by a water molecule activated by the general base functionality of Glu136 has been proposed. Here, we report the first crystal structures of SGP in complex with two transition state peptide analogs designed to mimic the tetrahedral intermediate of the proteolytic reaction. Of these two analogs, the one containing a central S-hydroxyl group is a potent sub-nanomolar inhibitor of SGP. The inhibitor binds non-covalently to the concave surface of the upper beta-sheet and enables delineation of the S4 to S3' substrate specificity pockets of the enzyme. Structural differences in these pockets account for the unique substrate preferences of SGP among peptidases having an acidic pH optimum. Inhibitor binding is accompanied by a structuring of the region comprising residues Tyr71-Gly80 from being mostly disordered in the apoenzyme and leading to positioning of crucial active site residues for establishing enzyme-inhibitor contacts. In addition, conformational rearrangements are seen in a disulfide bridged surface loop (Cys141-Cys148), which moves inwards, partially closing the open substrate binding cleft of the native enzyme. The non-hydrolysable scissile bond analog of the inhibitor is located in the active site forming close contacts with Gln53 and Glu136. The nucleophilic water molecule is displaced and a unique mode of binding is observed with the S-OH of the inhibitor occupying the oxyanion binding site of the proposed tetrahedral intermediate. Details of the enzyme-inhibitor interactions and mechanistic interpretations are discussed.     

Structural Analysis of Peroxide-Soaked MnSOD Crystals Reveals Side-On Binding of Peroxide to Active-Site Manganese

The superoxide dismutase (SOD) enzymes are important antioxidant agents that protect cells from reactive oxygen species. The SOD family is responsible for catalyzing the disproportionation of superoxide radical to oxygen and hydrogen peroxide. Manganese- and iron-containing SOD exhibit product inhibition whereas Cu/ZnSOD does not. Here, we report the crystal structure of Escherichia coli MnSOD with hydrogen peroxide cryotrapped in the active site. Crystallographic refinement to 1.55 Å and close inspection revealed electron density for hydrogen peroxide in three of the four active sites in the asymmetric unit. The hydrogen peroxide molecules are in the position opposite His26 that is normally assumed by water in the trigonal bipyramidal resting state of the enzyme. Hydrogen peroxide is present in active sites B, C, and D and is side-on coordinated to the active-site manganese. In chains B and D, the peroxide is oriented in the plane formed by manganese and ligands Asp167 and His26. In chain C, the peroxide is bound, making a 70° angle to the plane. Comparison of the peroxide-bound active site with the hydroxide-bound octahedral form shows a shifting of residue Tyr34 towards the active site when peroxide is bound. Comparison with peroxide-soaked Cu/ZnSOD indicates end-on binding of peroxide when the SOD does not exhibit inhibition by peroxide and side-on binding of peroxide in the product-inhibited state of MnSOD.