Failures of inverse folding and threading with gapped alignment

To calculate the tertiary structure of a protein from its amino acid sequence, the thermodynamic approach requires a potential function of sequence and conformation that has its global minimum at the native conformation for many different proteins. Here we study the behavior of such functions for the simplest model system that still has some of the features of the protein folding problem, namely two-dimensional square lattice chain configurations involving two residue types. First we show that even the given contact potential, which by definition is used to identify the folding sequences and their unique native conformations, cannot always correctly select which sequences will fold to a given structure. Second, we demonstrate that the given contact potential is not always able to favor the native alignment of a native sequence on its own native conformation over other gapped alignments of different folding sequences onto that same conformation. Because of these shortcomings, even in this simple model system in which all conformations and all native sequences are known and determined directly by the given potential, we must reexamine our expectations for empirical potentials used for inverse folding and gapped alignment on more realistic representations of proteins. © 1996 Wiley-Liss, Inc.     

Amino acid sequence of a major globin from the sea cucumber Paracaudina chilensis

The sea cucumber Paracaudina chilensis (Echinodermata) contains three major globins I, II and III in coelomic cells. The complete amino acid sequence of globin I has been determined. It is composed of 157 amino acid residues, is acetylated at the N-terminus, and has a characteristic N-terminal extension of 9-10 residues when compared with vertebrate globins. The sequence of Paracaudina globin I showed slightly higher homology with human alpha globin (25%) rather than with the invertebrate Anadara alpha globin (22%). Paracaudina globin I also showed strong homology (59%) with globin D from another sea cucumber, Molpadia arenicola (Mauri, F.C. (1985) Ph.D. dissertation, University of Texas). The globin sequences from the phylum Echinodermata have an important position in the molecular evolution of the globins, because they are the invertebrate group most closely related to the vertebrates.     

Automated use of mutagenesis data in structure prediction

In the absence of experimental structural determination, numerous methods are available to indirectly predict or probe the structure of a target molecule. Genetic modification of a protein sequence is a powerful tool for identifying key residues involved in binding reactions or protein stability. Mutagenesis data is usually incorporated into the modeling process either through manual inspection of model compatibility with empirical data, or through the generation of geometric constraints linking sensitive residues to a binding interface. We present an approach derived from statistical studies of lattice models for introducing mutation information directly into the fitness score. The approach takes into account the phenotype of mutation (neutral or disruptive) and calculates the energy for a given structure over an ensemble of sequences. The structure prediction procedure searches for the optimal conformation where neutral sequences either have no impact or improve stability and disruptive sequences reduce stability relative to wild type. We examine three types of sequence ensembles: information from saturation mutagenesis, scanning mutagenesis, and homologous proteins. Incorporating multiple sequences into a statistical ensemble serves to energetically separate the native state and misfolded structures. As a result, the prediction of structure with a poor force field is sufficiently enhanced by mutational information to improve accuracy. Furthermore, by separating misfolded conformations from the target score, the ensemble energy serves to speed up conformational search algorithms such as Monte Carlo-based methods.     

Burial of nonpolar surface area and thermodynamic stabilization of globins as a function of chain elongation

Proteins are biosynthesized from N to C terminus before they depart from the ribosome and reach their bioactive state in the cell. At present, very little is known about the evolution of conformation and the free energy of the nascent protein with chain elongation. These parameters critically affect the extent of folding during ribosome‐assisted biosynthesis. Here, we address the impact of vectorial amino acid addition on the burial of nonpolar surface area and on the free energy of native‐like structure formation in the absence of the ribosomal machinery. We focus on computational predictions on proteins bearing the globin fold, which is known to encompass the 3/3, 2/2, and archaeal subclasses. We find that the burial of nonpolar surface increases progressively with chain elongation, leading to native‐like conformations upon addition of the last C‐terminal residues, corresponding to incorporation of the last two helices. Additionally, the predicted folding entropy for generating native‐like structures becomes less unfavorable at nearly complete chain lengths, suggesting a link between the late burial of nonpolar surface and water release. Finally, the predicted folding free energy takes a progressive favorable dip toward more negative values, as the chain gets longer. These results suggest that thermodynamic stabilization of the native structure of newly synthesized globins during translation in the cell is significantly enhanced as the chain elongates. This is especially true upon departure of the last C‐terminal residues from the ribosomal tunnel, which hosts ca., 30–40 amino acids. Hence, we propose that release from the ribosome is a crucial step in the life of single‐domain proteins in the cell. Proteins 2014; 82:2318–2331. © 2014 Wiley Periodicals, Inc.     

Globin gene family evolution and functional diversification in annelids

Globins are the most common type of oxygen-binding protein in annelids. In this paper, we show that circulating intracellular globin (Alvinella pompejana and Glycera dibranchiata), noncirculating intracellular globin (Arenicola marina myoglobin) and extracellular globin from various annelids share a similar gene structure, with two conserved introns at canonical positions B12.2 and G7.0. Despite sequence divergence between intracellular and extracellular globins, these data strongly suggest that these three globin types are derived from a common ancestral globin-like gene and evolved by duplication events leading to diversification of globin types and derived functions. A phylogenetic analysis shows a distinct evolutionary history of annelid extracellular hemoglobins with respect to intracellular annelid hemoglobins and mollusc and arthropod extracellular hemoglobins. In addition, dehaloperoxidase (DHP) from the annelid, Amphitrite ornata, surprisingly exhibits close phylogenetic relationships to some annelid intracellular globins. We have characterized the gene structure of A. ornata DHP to confirm assumptions about its homology with globins. It appears that it has the same intron position as in globin genes, suggesting a common ancestry with globins. In A. ornata, DHP may be a derived globin with an unusual enzymatic function.     

DockRank: Ranking docked conformations using partner‐specific sequence homology‐based protein interface prediction

Selecting near‐native conformations from the immense number of conformations generated by docking programs remains a major challenge in molecular docking. We introduce DockRank, a novel approach to scoring docked conformations based on the degree to which the interface residues of the docked conformation match a set of predicted interface residues. DockRank uses interface residues predicted by partner‐specific sequence homology‐based protein–protein interface predictor (PS‐HomPPI), which predicts the interface residues of a query protein with a specific interaction partner. We compared the performance of DockRank with several state‐of‐the‐art docking scoring functions using Success Rate (the percentage of cases that have at least one near‐native conformation among the top m conformations) and Hit Rate (the percentage of near‐native conformations that are included among the top m conformations). In cases where it is possible to obtain partner‐specific (PS) interface predictions from PS‐HomPPI, DockRank consistently outperforms both (i) ZRank and IRAD, two state‐of‐the‐art energy‐based scoring functions (improving Success Rate by up to 4‐fold); and (ii) Variants of DockRank that use predicted interface residues obtained from several protein interface predictors that do not take into account the binding partner in making interface predictions (improving success rate by up to 39‐fold). The latter result underscores the importance of using partner‐specific interface residues in scoring docked conformations. We show that DockRank, when used to re‐rank the conformations returned by ClusPro, improves upon the original ClusPro rankings in terms of both Success Rate and Hit Rate. DockRank is available as a server at http://einstein.cs.iastate.edu/DockRank/ . Proteins 2014; 82:250–267. © 2013 Wiley Periodicals, Inc.     

Loop propensity of the sequence YKGQP from staphylococcal nuclease: implications for the folding of nuclease.

Recently we performed molecular dynamics (MD) simulations on the folding of the hairpin peptide DTVKLMYKGQPMTFR from staphylococcal nuclease in explicit water. We found that the peptide folds into a hairpin conformation with native and nonnative hydrogen-bonding patterns. In all the folding events observed in the folding of the hairpin peptide, loop formation involving the region YKGQP was an important event. In order to trace the origins of the loop propensity of the sequence YKGQP, we performed MD simulations on the sequence starting from extended, polyproline II and native type I' turn conformations for a total simulation length of 300 ns, using the GROMOS96 force field under constant volume and temperature (NVT) conditions. The free-energy landscape of the peptide YKGQP shows minima corresponding to loop conformation with Tyr and Pro side-chain association, turn and extended conformational forms, with modest free-energy barriers separating the minima. To elucidate the role of Gly in facilitating loop formation, we also performed MD simulations of the mutated peptide YKAQP (Gly --> Ala mutation) under similar conditions starting from polyproline II conformation for 100 ns. Two minima corresponding to bend/turn and extended conformations were observed in the free-energy landscape for the peptide YKAQP. The free-energy barrier between the minima in the free-energy landscape of the peptide YKAQP was also modest. Loop conformation is largely sampled by the YKGQP peptide, while extended conformation is largely sampled by the YKAQP peptide. We also explain why the YKGQP sequence samples type II turn conformation in these simulations, whereas the sequence as part of the hairpin peptide DTVKLMYKGQPMTFR samples type I' turn conformation both in the X-ray crystal structure and in our earlier simulations on the folding of the hairpin peptide. We discuss the implications of our results to the folding of the staphylococcal nuclease.     

Axolotl hemoglobin: cDNA-derived amino acid sequences of two alpha globins and a beta globin from an adult Ambystoma mexicanum

Erythrocytes of the adult axolotl, Ambystoma mexicanum, have multiple hemoglobins. We separated and purified two kinds of hemoglobin, termed major hemoglobin (Hb M) and minor hemoglobin (Hb m), from a five-year-old male by hydrophobic interaction column chromatography on Alkyl Superose. The hemoglobins have two distinct alpha type globin polypeptides (alphaM and alpham) and a common beta globin polypeptide, all of which were purified in FPLC on a reversed-phase column after S-pyridylethylation. The complete amino acid sequences of the three globin chains were determined separately using nucleotide sequencing with the assistance of protein sequencing. The mature globin molecules were composed of 141 amino acid residues for alphaM globin, 143 for alpham globin and 146 for beta globin. Comparing primary structures of the five kinds of axolotl globins, including two previously established alpha type globins from the same species, with other known globins of amphibians and representatives of other vertebrates, we constructed phylogenetic trees for amphibian hemoglobins and tetrapod hemoglobins. The molecular trees indicated that alphaM, alpham, beta and the previously known alpha major globin were adult types of globins and the other known alpha globin was a larval type. The existence of two to four more globins in the axolotl erythrocyte is predicted.     

Identification of native protein folds amongst a large number of incorrect models. The calculation of low energy conformations from potentials of mean force

We present an approach that is able to detect native folds amongst a large number of non-native conformations. The method is based on the compilation of potentials of mean force of the interactions of the C beta atoms of all amino acid pairs from a database of known three-dimensional protein structures. These potentials are used to calculate the conformational energy of amino acid sequences in a number of different folds. For a substantial number of proteins we find that the conformational energy of the native state is lowest amongst the alternatives. Exceptions are proteins containing large prosthetic groups, Fe-S clusters or polypeptide chains that do not adopt globular folds. We discuss briefly potential applications in various fields of protein structural research.     

Genome-wide analysis of the globin gene family of C. elegans

The aim of our study was to annotate sequences for 35 putative globins from the nematode Caenorhabditis elegans. All these proteins are expressed, but seven of these differ from the gene predictions in Wormbase. The entire polypeptide sequences for 31 genes and the core globin domain of four proteins were confirmed or corrected. All core globin domains were aligned manually following a procedure that was designed to fit the putative sequences to the crystal structure based alignment of 56 known globin crystal structures. Neighbor-joining analysis of the resulting alignment showed that the majority of these globins are very divergent from each other, possibly suggesting a long evolutionary divergence. The surprisingly high number and low sequence conservation of putative globins in this small organism urges a detailed functional analysis.     

Recombinant protein sequences can trigger methylation of N-terminal amino acids in Escherichia coli.

Recombinant human hemoglobin rHb1.1 has been genetically engineered with the replacement of the wild-type valine residues at all N-termini with methionine, an Asn 108 Lys substitution on the beta globins, and a fusion of the two alpha globins with a glycine linker. When rHb1.1 was expressed in Escherichia coli, methylation of the N-terminal methionine of the alpha globin was discovered. Another mutant has been engineered with the alpha globin gene coding for N-terminal methionine followed by an insertion of alanine. Characterization of expressed hemoglobin from this variant revealed a methylated N-terminal alanine that occurred through two posttranslational events: initial excision of the N-terminal methionine, followed by methylation of alanine as the newly generated N-terminus. No methylation was observed for variants expressed with wild-type valine at the N-terminus of the alpha globin. The methylation of N-terminal amino acids was attributed to a specific protein sequence that can trigger methylation of proteins expressed in E. coli. Here we demonstrate that proline at position 4 in the protein sequence of alpha globin seems an essential part of that signaling. Although N-terminal methylation has been observed previously for native E. coli proteins with similar N-terminal sequences, methylation of the recombinant globins has allowed further delineation of the recognition sequence, and indicates that methylation of heterologous proteins can occur in E. coli.     

Loop modeling: Sampling, filtering, and scoring

We describe a fast and accurate protocol, LoopBuilder, for the prediction of loop conformations in proteins. The procedure includes extensive sampling of backbone conformations, side chain addition, the use of a statistical potential to select a subset of these conformations, and, finally, an energy minimization and ranking with an all-atom force field. We find that the Direct Tweak algorithm used in the previously developed LOOPY program is successful in generating an ensemble of conformations that on average are closer to the native conformation than those generated by other methods. An important feature of Direct Tweak is that it checks for interactions between the loop and the rest of the protein during the loop closure process. DFIRE is found to be a particularly effective statistical potential that can bias conformation space toward conformations that are close to the native structure. Its application as a filter prior to a full molecular mechanics energy minimization both improves prediction accuracy and offers a significant savings in computer time. Final scoring is based on the OPLS/SBG-NP force field implemented in the PLOP program. The approach is also shown to be quite successful in predicting loop conformations for cases where the native side chain conformations are assumed to be unknown, suggesting that it will prove effective in real homology modeling applications.     

Free-energy landscape of the villin headpiece in an all-atom force field

We investigate the landscape of the internal free-energy of the 36 amino acid villin headpiece with a modified basin hopping method in the all-atom force field PFF01, which was previously used to predictively fold several helical proteins with atomic resolution. We identify near native conformations of the protein as the global optimum of the force field. More than half of the twenty best simulations started from random initial conditions converge to the folding funnel of the native conformation, but several competing low-energy metastable conformations were observed. From 76,000 independently generated conformations we derived a decoy tree which illustrates the topological structure of the entire low-energy part of the free-energy landscape and characterizes the ensemble of metastable conformations. These emerge as similar in secondary content, but differ in tertiary arrangement.     

All-atom folding of the three-helix HIV accessory protein with an adaptive parallel tempering method

All-atom protein structure prediction from the amino acid sequence alone remains an important goal of biophysical chemistry. Recent progress in force field development and validation suggests that the PFF01 free-energy force field correctly predicts the native conformation of various helical proteins as the global optimum of its free-energy surface. Reproducible protein structure prediction requires the availability of efficient optimization methods to locate the global minima of such complex potentials. Here we investigate an adapted version of the parallel tempering method as an efficient parallel stochastic optimization method for protein structure prediction. Using this approach we report the reproducible all-atom folding of the three-helix 40 amino acid HIV accessory protein from random conformations to within 2.4 A backbone RMS deviation from the experimental structure with modest computational resources.     

Distinguishing native conformations of proteins from decoys with an effective free energy estimator based on the OPLS all-atom force field and the Surface Generalized Born solvent model

Protein decoy data sets provide a benchmark for testing scoring functions designed for fold recognition and protein homology modeling problems. It is commonly believed that statistical potentials based on reduced atomic models are better able to discriminate native-like from misfolded decoys than scoring functions based on more detailed molecular mechanics models. Recent benchmark tests on small data sets, however, suggest otherwise. In this work, we report the results of extensive decoy detection tests using an effective free energy function based on the OPLS all-atom (OPLS-AA) force field and the Surface Generalized Born (SGB) model for the solvent electrostatic effects. The OPLS-AA/SGB effective free energy is used as a scoring function to detect native protein folds among a total of 48,832 decoys for 32 different proteins from Park and Levitt's 4-state-reduced, Levitt's local-minima, Baker's ROSETTA all-atom, and Skolnick's decoy sets. Solvent electrostatic effects are included through the Surface Generalized Born (SGB) model. All structures are locally minimized without restraints. From an analysis of the individual energy components of the OPLS-AA/SGB energy function for the native and the best-ranked decoy, it is determined that a balance of the terms of the potential is responsible for the minimized energies that most successfully distinguish the native from the misfolded conformations. Different combinations of individual energy terms provide less discrimination than the total energy. The results are consistent with observations that all-atom molecular potentials coupled with intermediate level solvent dielectric models are competitive with knowledge-based potentials for decoy detection and protein modeling problems such as fold recognition and homology modeling.     

Inferring the conformation of RNA base pairs and triples from patterns of sequence variation.

The success of comparative analysis in resolving RNA secondary structure and numerous tertiary interactions relies on the presence of base covariations. Although the majority of base covariations in aligned sequences is associated to Watson-Crick base pairs, many involve non-canonical or restricted base pair exchanges (e.g. only G:C/A:U), reflecting more specific structural constraints. We have developed a computer program that determines potential base pairing conformations for a given set of paired nucleotides in a sequence alignment. This program (ISOPAIR) assumes that the base pair conformation is maintained through sequence variation without significantly affecting the path of the sugar-phosphate backbone. ISOPAIR identifies such 'isomorphic' structures for any set of input base pair or base triple sequences. The program was applied to base pairs and triples with known structures and sequence exchanges. In several instances, isomorphic structures were correctly identified with ISOPAIR. Thus, ISOPAIR is useful when assessing non-canonical base pair conformations in comparative analysis. ISOPAIR applications are limited to those cases where unusual base pair exchanges indeed reflect a non-canonical conformation.     

A differential geometric treatment of protein structure comparison

The technique of model-building a protein of known sequence but unknown tertiary structure from the structures of homologous proteins is probably so far the most reliable means of mapping from primary to tertiary structure. A key step towards the realization of the aim is to develop ways of aligning three-dimensional structures of homologus proteins, thereby deriving the rules useful for protein modelling. We have developed a generalized differential-geometric representation of protein local conformation for use in a protein comparison program which aligns protein sequences on the basis of their sequence and conformational knowledge. Because the differetial-geometric distance measure between local conformations is independent of the coordinate frame and remains chirality information, the comparison program is easily implemented, relatively rational and reasonably fast. The utility of this program for aligning closely and distantly related homologous proteins is demonstrated by multiple alignment of globins, serine proteinases and aspartic proteinase domains. Particularly, the method has reached the rational alignment between the mammalian and microbial serine proteinases as compared with many published alignment programs.     

A new protein folding recognition potential function

On the study of protein inverse folding problem, one goal is to find simple and efficient potential to evaluate the compatibility between structure and a given sequence. We present here a novo empirical mean force potential to address the importance of electrostatic interactions in protein inverse folding study. It is based on protein main chain polar fraction and constructed in a way similar with Sippl's from a database of 64 known independent three-dimensional protein structures. This potential was applied to recognize the protein native conformations among a conformation pool. Calculated results show that this potential is powerful in picking out native conformations, in addition it can also find structure similarity between proteins with low sequence similarity. The success of this new potential clearly shows the importance of electrostatic factors in protein inverse folding studies. © 1995 Wiley-Liss, Inc.