Quantifying the impact of protein tertiary structure on molecular evolution

To investigate the evolutionary impact of protein structure, the experimentally determined tertiary structure and the protein-coding DNA sequence were collected for each of 1,195 genes. These genes were studied via a model of sequence change that explicitly incorporates effects on evolutionary rates due to protein tertiary structure. In the model, these effects act via the solvent accessibility environments and pairwise amino acid interactions that are induced by tertiary structure. To compare the hypotheses that structure does and does not have a strong influence on evolution, Bayes factors were estimated for each of the 1,195 sequences. Most of the Bayes factors strongly support the hypothesis that protein structure affects protein evolution. Furthermore, both solvent accessibility and pairwise interactions among amino acids are inferred to have important roles in protein evolution. Our results also indicate that the strength of the relationship between tertiary structure and evolution has a weak but real correlation to the annotation information in the Gene Ontology database. Although their influences on rates of evolution vary among protein families, we find that the mean impacts of solvent accessibility and pairwise interactions are about the same.     

Neuro-fuzzy structural classification of proteins for improved protein secondary structure prediction

Fourier transform infrared (FTIR) spectroscopy is a very flexible technique for characterization of protein secondary structure. Measurements can be carried out rapidly in a number of different environments based on only small quantities of proteins. For this technique to become more widely used for protein secondary structure characterization, however, further developments in methods to accurately quantify protein secondary structure are necessary. Here we propose a structural classification of proteins (SCOP) class specialized neural networks architecture combining an adaptive neuro-fuzzy inference system (ANFIS) with SCOP class specialized backpropagation neural networks for improved protein secondary structure prediction. Our study shows that proteins can be accurately classified into two main classes "all alpha proteins" and "all beta proteins" merely based on the amide I band maximum position of their FTIR spectra. ANFIS is employed to perform the classification task to demonstrate the potential of this architecture with moderately complex problems. Based on studies using a reference set of 17 proteins and an evaluation set of 4 proteins, improved predictions were achieved compared to a conventional neural network approach, where structure specialized neural networks are trained based on protein spectra of both "all alpha" and "all beta" proteins. The standard errors of prediction (SEPs) in % structure were improved by 4.05% for helix structure, by 5.91% for sheet structure, by 2.68% for turn structure, and by 2.15% for bend structure. For other structure, an increase of SEP by 2.43% was observed. Those results were confirmed by a "leave-one-out" run with the combined set of 21 FTIR spectra of proteins.     

Crowding and hydration effects on protein conformation: a study with sol-gel encapsulated proteins.

We are developing an experimental system for testing the effects of macromolecular crowding and molecular confinement on protein structure. In the present study, solvent effects on the secondary structure of two proteins were examined by circular dichroism following encapsulation in the hydrated pores of a silica glass matrix by the sol-gel method. Changes in the unfolded conformations of encapsulated apomyoglobin and reduced serum albumin were analyzed after equilibration with aqueous solutions of natural osmolytes, short-chain alcohols, polyethylene glycol, and a complete series of Hofmeister cations. In many instances, the alpha-helical content of the encapsulated protein was increased by addition of solutes at concentrations that have no effect on the protein in the absence of the glass. The results are discussed from the perspective of water structure. We argue that perturbed water at the silica interface causes an increase in the average free energy of the bulk water phase which, consequently, diminishes the strength of the hydrophobic effect inside the glass matrix and destabilizes the conformation of encapsulated proteins. We propose that solutes can increase the strength of the hydrophobic effect and influence folding equilibria without directly interacting with the protein. A hypothesis is provided for the apparent paradox that kosmotropic (strongly water binding) anions favor native protein structure, whereas chaotropic (weakly water binding) cations enhance native protein structure. The encapsulation results suggest that macromolecular crowding and molecular confinement are accompanied by hydration effects that may oppose or potentiate the stabilizing effects of excluded volume on protein structure, depending on the surface chemistry of the crowding agent and its influence on bulk water structure. In the crowded environment of a living cell, excluded volume effects, surface-induced water structure, and compatible solutes are expected to complement the dominant forces in protein folding.     

Phyletic relationships of protein structures based on spatial preference of residues

Summary A structure-based scoring matrix MDPRE was derived from amino acid spatial preferences in protein structures. Sequence alignment and evolutionary studies by using MDPRE matrix gave similar results as those from ordinary sequence and structure alignments. It is interesting that a matrix derived from structure data solely could give comparable alignment results, strongly indicating the intimate connection between protein sequences and structures. The branch order and length from this approach were close to those obtained by a structure comparison method. Thus, by applying this structure-based matrix, the trees obtained should reflect evolutionary characteristics of protein structure. This approach takes advantage over a direct structure comparison in that (1) only a sequence and MDPRE matrix are needed, making it simple and widely applicable (especially in the absence of 3-dimensional protein structure data); (2) an established algorithm for sequence alignment and tree building could be employed, providing opportunities for direct comparison between matrices from different methodologies. One of the most striking features of this method is its capability to detect protein structure homologies when the sequence identities are low. This was well reflected in the given examples of the alignment of dinucleotidebinding domains.     

The effect of backbone on the small-world properties of protein contact maps

In the last years, small-world behavior has been extensively described for proteins, when they are represented by the undirected graph defined by the inter-residue protein contacts. By adopting this representation it was possible to compute the average clustering coefficient (C) and characteristic path length (L) of protein structures, and their values were found to be similar to those of graphs characterized by small-world topology. In this comment, we analyze a large set of non-redundant protein structures (1753) and show that by randomly mimicking the protein collapse, the covalent structure of the protein chain significantly contributes to the small-world behavior of the inter-residue contact graphs. When protein graphs are generated, imposing constraints similar to those induced by the backbone connectivity, their characteristic path lengths and clustering coefficients are indistinguishable from those computed using the real contact maps showing that L and C values cannot be used for 'protein fingerprinting'. Moreover we verified that these results are independent of the selected protein representations, residue composition and protein secondary structures.     

Novel protein folds and their nonsequential structural analogs

Newly determined protein structures are classified to belong to a new fold, if the structures are sufficiently dissimilar from all other so far known protein structures. To analyze structural similarities of proteins, structure alignment tools are used. We demonstrate that the usage of nonsequential structure alignment tools, which neglect the polypeptide chain connectivity, can yield structure alignments with significant similarities between proteins of known three-dimensional structure and newly determined protein structures that possess a new fold. The recently introduced protein structure alignment tool, GANGSTA, is specialized to perform nonsequential alignments with proper assignment of the secondary structure types by focusing on helices and strands only. In the new version, GANGSTA+, the underlying algorithms were completely redesigned, yielding enhanced quality of structure alignments, offering alignment against a larger database of protein structures, and being more efficient. We applied DaliLite, TM-align, and GANGSTA+ on three protein crystal structures considered to be novel folds. Applying GANGSTA+ to these novel folds, we find proteins in the ASTRAL40 database, which possess significant structural similarities, albeit the alignments are nonsequential and in some cases involve secondary structure elements aligned in reverse orientation. A web server is available at http://agknapp.chemie.fu-berlin.de/gplus for pairwise alignment, visualization, and database comparison.     

The secondary structure of gap junctions. Influence of isolation methods and proteolysis.

Paired intercellular transmembrane channels, termed connexons, comprised of hexameric assemblies of gap junction protein, were isolated and purified from rat liver by exploiting their resistance to either Sarkosyl detergent solubilization or alkali extraction. The secondary structures of the gap junction proteins prepared by these methods were compared by circular dichroism (CD) spectroscopy. Both the spectra and the calculated net secondary structures of the proteins obtained by the two isolation methods were different. The protein isolated by the Sarkosyl treatment was found to be approximately 50% alpha-helical, while protein isolated by alkali extraction had a lower helix content (approximately 40%). In both types of preparations, however, the helical content of the gap junction protein was sufficiently large to be consistent with an all-helical model for the membrane-spanning parts of the structure. CD spectroscopy was also used to examine the effects of proteolytic digestion of the cytoplasmic domain on the net secondary structure of the detergent-treated gap junction protein. The membrane-bound fragments had a slightly higher proportion of their residues that were alpha-helical in nature, suggesting that the transmembrane and/or intra-gap domains are indeed enriched in this type of secondary structure. This information constrains the range of models which can be realistically proposed for the channel structure.     

Local interactions drive the formation of nonnative structure in the denatured state of human alpha-lactalbumin: a high resolution structural characterization of a peptide model in aqueous solution.

There are a small number of peptides derived from proteins that have a propensity to adopt structure in aqueous solution which is similar to the structure they possess in the parent protein. There are far fewer examples of protein fragments which adopt stable nonnative structures in isolation. Understanding how nonnative interactions are involved in protein folding is crucial to our understanding of the topic. Here we show that a small, 11 amino acid peptide corresponding to residues 101-111 of the protein alpha-lactalbumin is remarkably structured in isolation in aqueous solution. The peptide has been characterized by 1H NMR, and 170 ROE-derived constraints were used to calculate a structure. The calculations yielded a single, high-resolution structure for residues 101-107 that is nonnative in both the backbone and side-chain conformations. In the pH 6.5 crystal structure, residues 101-105 are in an irregular turn-like conformation and residues 106-111 form an alpha-helix. In the pH 4.2 crystal structure, residues 101-105 form an alpha-helix, and residues 106-111 form a loopike structure. Both of these structures are significantly different from the conformation adopted by our peptide. The structure in the peptide model is primarily the result of local side-chain interactions that force the backbone to adopt a nonnative 310/turn-like structure in residues 103-106. The structure in aqueous solution was compared to the structure in 30% trifluoroethanol (TFE), and clear differences were observed. In particular, one of the side-chain interactions, a hydrophobic cluster involving residues 101-105, is different in the two solvents and residues 107-111 are considerably more ordered in 30% TFE. The implications of the nonnative structure for the folding of alpha-lactalbumin is discussed.     

Denatured-state energy landscapes of a protein structural database reveal the energetic determinants of a framework model for folding

Position-specific denatured-state thermodynamics were determined for a database of human proteins by use of an ensemble-based model of protein structure. The results of modeling denatured protein in this manner reveal important sequence-dependent thermodynamic properties in the denatured ensembles as well as fundamental differences between the denatured and native ensembles in overall thermodynamic character. The generality and robustness of these results were validated by performing fold-recognition experiments, whereby sequences were matched with their respective folds based on amino acid propensities for the different energetic environments in the protein, as determined through cluster analysis. Correlation analysis between structure and energetic information revealed that sequence segments destined for β-sheet in the final native fold are energetically more predisposed to a broader repertoire of states than are sequence segments destined for α-helix. These results suggest that within the subensemble of mostly unstructured states, the energy landscapes are dominated by states in which parts of helices adopt structure, whereas structure formation for sequences destined for β-strand is far less probable. These results support a framework model of folding, which suggests that, in general, the denatured state has evolutionarily evolved to avoid low-energy conformations in sequences that ultimately adopt β-strand. Instead, the denatured state evolved so that sequence segments that ultimately adopt α-helix and coil will have a high intrinsic structure formation capability, thus serving as potential nucleation sites.     

Angle-distance image matching techniques for protein structure comparison

Proteins that contain similar structural elements often have analogous functions regardless of the degree of sequence similarity or structure connectivity in space. In general, protein structure comparison (PSC) provides a straightforward methodology for biologists to determine critical aspects of structure and function. Here, we developed a novel PSC technique based on angle-distance image (A-D image) transformation and matching, which is independent of sequence similarity and connectivity of secondary structure elements (SSEs). An A-D image is constructed by utilizing protein secondary structure information. According to various types of SSEs, the mutual SSE pairs of the query protein are classified into three different types of sub-images. Subsequently, corresponding sub-images between query and target protein structures are compared using modified cross-correlation approaches to identify the similarity of various patterns. Structural relationships among proteins are displayed by hierarchical clustering trees, which facilitate the establishment of the evolutionary relationships between structure and function of various proteins.Four standard testing datasets and one newly created dataset were used to evaluate the proposed method. The results demonstrate that proteins from these five datasets can be categorized in conformity with their spatial distribution of SSEs. Moreover, for proteins with low sequence identity that share high structure similarity, the proposed algorithms are an efficient and effective method for structural comparison.     

Using phosphorothioate-substituted RNA to investigate the thermodynamic role of phosphates in a sequence specific RNA-protein complex

Part of the binding affinity and specificity in RNA-protein complexes is often contributed by contacts between the protein and backbone phosphates that are held in position by the RNA structure. This study focuses on the well-characterized interaction between a dimer of the MS2 coat protein and a small RNA hairpin. Using a short oligoribonucleotide which contains all the necessary sequence elements required for tight protein binding, a single phosphorothioate linkage was introduced at 13 different positions. In each case, the R(P) and S(P) stereoisomers were separated and their affinities to the MS2 coat protein were determined. Comparison of these biochemical data with the crystal structure of the protein-hairpin complex indicates that introduction of a phosphorothioate only affects binding at sites where a protein-phosphate contact is observed in the crystal structure. This means that phosphorothioate-containing oligoribonucleotides should also be useful for mapping phosphate contacts in RNA-protein complexes for which no crystal structure is available.     

High-resolution structural validation of the computational redesign of human U1A protein

Achieving atomic-level resolution in the computational design of a protein structure remains a challenging problem despite recent progress. Rigorous experimental tests are needed to improve protein design algorithms, yet studies of the structure and dynamics of computationally designed proteins are very few. The NMR structure and backbone dynamics of a redesigned protein of 96 amino acids are compared here with the design target, human U1A protein. We demonstrate that the redesigned protein reproduces the target structure to within the uncertainty of the NMR coordinates, even as 65 out of 96 amino acids were simultaneously changed by purely computational methods. The dynamics of the backbone of the redesigned protein also mirror those of human U1A, suggesting that the protein design algorithm captures the shape of the potential energy landscape in addition to the local energy minimum.     

Evaluation of comparative protein modeling by MODELLER

We evaluate 3D models of human nucleoside diphosphate kinase, mouse cellular retinoic acid binding protein I, and human eosinophil neurotoxin that were calculated by MODELLER , a program for comparative protein modeling by satisfaction of spatial restraints. The models have good stereochemistry and are at least as similar to the crystallographic structures as the closest template structures. The largest errors occur in the regions that were not aligned correctly or where the template structures are not similar to the correct structure. These regions correspond predominantly to exposed loops, insertions of any length, and non-conserved side chains. When a template structure with more than 40% sequence identity to the target protein is available, the model is likely to have about 90% of the mainchain atoms modeled with an rms deviation from the X-ray structure of ≈ 1 Å, in large part because the templates are likely to be that similar to the X-ray structure of the target. This rms deviation is comparable to the overall differences between refined NMR and X-ray crystallography structures of the same protein. © 1995 Wiley-Liss, Inc.     

Design of lambda Cro fold: solution structure of a monomeric variant of the de novo protein

One of the classical DNA-binding proteins, bacteriophage lambda Cro, forms a homodimer with a unique fold of alpha-helices and beta-sheets. We have computationally designed an artificial sequence of 60 amino acid residues to stabilize the backbone tertiary structure of the lambda Cro dimer by simulated annealing using knowledge-based structure-sequence compatibility functions. The designed amino acid sequence has 25% identity with that of natural lambda Cro and preserves Phe58, which is important for formation of the stably folded structure of lambda Cro. The designed dimer protein and its monomeric variant, which was redesigned by the insertion of a beta-hairpin sequence at the C-terminal region to prevent dimerization, were synthesized and biochemically characterized to be well folded. The designed protein was monomeric under a wide range of protein concentrations and its solution structure was determined by NMR spectroscopy. The solved structure is similar to that of a monomeric variant of natural lambda Cro with a root-mean-square deviation of the polypeptide backbones at 2.1A and has a well-packed protein core. Thus, our knowledge-based functions provide approximate but essential relationships between amino acid sequences and protein structures, and are useful for finding novel sequences that are foldable into a given target structure.     

Visualizing transiently associated catalytic domains in assembly-line biosynthesis using cryo-electron microscopy

While cryo-electron microscopy (cryo-EM) has revolutionized the structure determination of supramolecular protein complexes that are refractory to structure determination by X-ray crystallography, structure determination by cryo-EM can nonetheless be complicated by excessive conformational flexibility or structural heterogeneity resulting from weak or transient protein–protein association. Since such transient complexes are often critical for function, specialized approaches must be employed for the determination of meaningful structure–function relationships. Here, we outline examples in which transient protein–protein interactions have been visualized successfully by cryo-EM in the biosynthesis of fatty acids, polyketides, and terpenes. These studies demonstrate the utility of chemical crosslinking to stabilize transient protein–protein complexes for cryo-EM structural analysis, as well as the use of partial signal subtraction and localized reconstruction to extract useful structural information out of cryo-EM data collected from inherently dynamic systems. While these approaches do not always yield atomic resolution insights on protein–protein interactions, they nonetheless enable direct experimental observation of complexes in assembly-line biosynthesis that would otherwise be too fleeting for structural analysis.     

Conformational Transitions upon Ligand Binding: Holo-Structure Prediction from Apo Conformations

Biological function of proteins is frequently associated with the formation of complexes with small-molecule ligands. Experimental structure determination of such complexes at atomic resolution, however, can be time-consuming and costly. Computational methods for structure prediction of protein/ligand complexes, particularly docking, are as yet restricted by their limited consideration of receptor flexibility, rendering them not applicable for predicting protein/ligand complexes if large conformational changes of the receptor upon ligand binding are involved. Accurate receptor models in the ligand-bound state (holo structures), however, are a prerequisite for successful structure-based drug design. Hence, if only an unbound (apo) structure is available distinct from the ligand-bound conformation, structure-based drug design is severely limited. We present a method to predict the structure of protein/ligand complexes based solely on the apo structure, the ligand and the radius of gyration of the holo structure. The method is applied to ten cases in which proteins undergo structural rearrangements of up to 7.1 Å backbone RMSD upon ligand binding. In all cases, receptor models within 1.6 Å backbone RMSD to the target were predicted and close-to-native ligand binding poses were obtained for 8 of 10 cases in the top-ranked complex models. A protocol is presented that is expected to enable structure modeling of protein/ligand complexes and structure-based drug design for cases where crystal structures of ligand-bound conformations are not available.     

Nascent structure in the kinase anchoring domain of microtubule-associated protein 2

Biological processes are often viewed as highly ordered interactions between well-folded protein domains. The specific interactions exhibited by certain highly abundant neuronal proteins such as microtubule-associated protein 2 (MAP2) and tau stand in stark contrast because these proteins do not show evidence of structure by standard biophysical assays, yet they do bind to specific targets. It is conceivable that there are regions of MAP2 and tau with propensity to form structural domains upon binding a target. To search for evidence of such regions, limited proteolysis experiments were carried out on MAP2c, the smallest MAP2 isoform. Increased protease resistance was observed around the binding site for the RII subunit of cAMP-dependent protein kinase. Protein constructs spanning this region were produced based on the long-lived tryptic fragments Ser44-Arg93 and Ile94-Arg182, and were probed for structure using spectroscopic methods. The results support the existence of regions of nascent structure in the N-terminal region of MAP2c, which are believed to contribute to its regulatory function.     

木薯SYP基因的序列特征及表达分析

植物突触融合蛋白(SYP)是一类与植物细胞内囊泡介导转运有关的蛋白。部分SYP基因与植物对生物和非生物胁迫的响应有关。该文利用生物信息学工具分析了木薯(Manihot esculenta)SYP基因及其蛋白结构、核苷酸多态性和系统进化关系, 并利用RT-PCR技术检测了木薯不同组织中SYP基因的表达。结果表明, 木薯SYP基因及其蛋白结构均具有明显的规律性和家族成员间的保守性; SYP基因的cDNA在基因间以及不同品种间具高度一致性, 核苷酸变异以同义替换为主。进化分析表明, 植物SYP基因可分为2个亚家族, 木薯SYP基因倾向于与蓖麻(Ricinus communis)SYP基因聚在进化树同一分支的末端。半定量RT-PCR分析表明, 5个木薯SYP家族成员具有组织特异性。上述研究结果为木薯SYP基因功能研究和功能单核苷酸标记的开发奠定了重要基础。     

Activated conformations of the ras-gene-encoded p21 protein. 1. An energy-refined structure for the normal p21 protein complexed with GDP.

A complete three-dimensional structure for the ras-gene-encoded p21 protein with Gly 12 and Gln 61, bound to GDP, has been constructed in four stages using the available alpha-carbon coordinates as deposited in the Brookhaven National Laboratories Protein Data Bank. No all-atom structure has been made available despite the fact that the first crystallographic structure for the p21 protein was reported almost four years ago. In the p21 protein, if amino acid substitutions are made at any one of a number of different positions in the amino acid sequence, the protein becomes permanently activated and causes malignant transformation of normal cells or, in some cell lines, differentiation and maturation. For example, all amino acids except Gly and Pro at position 12 result in an oncogenic protein; all amino acids except Gln, Glu and Pro at position 61 likewise cause malignant transformation of cells. We have constructed our all-atom structure of the non-oncogenic protein from the x-ray structure in order to determine how oncogenic amino acid substitutions affect the three-dimensional structure of this protein. In Stage 1 we generated a poly-alanine backbone (except at Gly and Pro residues) through the alpha-carbon structure, requiring the individual Ala, Pro or Gly residues to conform to standard amino acid geometry and to form trans-planar peptide bonds. Since no alpha-carbon coordinates for residues 60-65 have been determined, these residues were modeled by generating them in the extended conformation and then subjecting them to molecular dynamics using the computer application DISCOVER and energy minimization using DISCOVER and the ECEPP (Empirical Conformational Energies for Peptides Program). In Stage 2, the positions of residues that are homologous to corresponding residues of bacterial elongation factor Tu (EF-Tu) to which p21 bears an overall 40% sequence homology, were determined from their corresponding positions in a high-resolution structure of EF-Tu. Non-homologous loops were taken from the structure generated in Stage 1 and were placed between the appropriate homologous segments so as to connect them. In Stage 3, all bad contacts that occurred in this resulting structure were removed, and the coordinates of the alpha-carbon atoms were forced to superimpose as closely as possible on the corresponding atoms of the reference (x-ray) structure. Then the side chain positions of residues of the non-homologous loop regions were modeled using a combination of molecular dynamics and energy minimization using DISCOVER and ECEPP respectively.(ABSTRACT TRUNCATED AT 400 WORDS)