RESUMO
Protein structure prediction has a number of important ad hoc similarity measures for evaluating predictions, but would benefit from a measure that is able to provide a common framework for a broad range of comparisons. Here we show that a mutual information-like measure can provide a comprehensive framework for evaluating protein structure prediction of all types. We discuss the concept of information, its application to secondary structure, and the obstacle to applying it to 3D structure. On the basis of the insights from the secondary structure case, we present an approach to work around the 3D difficulties, and develop a method to measure the mutual information provided by a 3D structure prediction. We integrate the evaluation of all types of protein structure prediction into a single framework, and compare the amount of information provided by various prediction methods, including secondary structure prediction. Within this broadened framework, the idea that structure is better preserved than sequence during evolution is evaluated quantitatively for the globin family. A nearly perfect sequence match in the globin family corresponds to about 300 bits of information, whereas a nearly perfect structural match for the same two proteins corresponds to about 2500 bits of information, where bits of information describes the probability of obtaining a match of similar closeness by chance. Mutual information provides both a theoretical basis for evaluating structure similarity and an explanatory surround for existing similarity measures.
Assuntos
Teoria da Informação , Conformação Proteica , Bases de Dados de Proteínas , Entropia , Proteínas/química , Análise de Sequência de Proteína , Relação Estrutura-AtividadeRESUMO
Template based protein structure prediction (commonly referred to as homology or comparative modeling) uses knowledge of solved structures to model a protein sequence's native or true fold. First, a parent structure is found and then a template structure is built by mapping the target sequence onto the parent structure. This putative structure is refined using a combination of backbone moves, side-chain packing, and loop modeling. Template based protein structure prediction has always held great promise to produce atomically accurate models close to the native conformation based on two major assumptions. First, similar sequences exhibit similar protein folds. Second, soluble proteins populate a discrete fold space with many representatives already solved in our Protein Data Bank (PDB). Ironically, beginning so close to the native structure is also the primary source of problems confronting this method and is the reason for the lack of progress in this category of structure prediction. In this review, the general concepts and procedures for template based structure prediction are outlined based on the following topics: sequence alignment, parent structure selection, template structure building, refinement, evaluation, and final structure selection. Then, a description of established software and algorithms is provided where the advantages and limitations of the different methods will be pointed out. This is followed by a discussion of the developments in template based structure prediction up to the 7th Critical Assessment of Structure Prediction meeting. Lastly, we will address the increased difficulty in improving templates that start so close to the native structure, and discuss the improvements needed in this field.
Assuntos
Biologia Computacional/métodos , Proteínas/química , Homologia Estrutural de Proteína , Algoritmos , Sequência de Aminoácidos , Animais , Humanos , Internet , Modelos Moleculares , Dados de Sequência MolecularRESUMO
We describe an information-theory-based measure of the quality of secondary structure prediction (RELINFO). RELINFO has a simple yet intuitive interpretation: it represents the factor by which secondary structure choice at a residue has been restricted by a prediction scheme. As an alternative interpretation of secondary structure prediction, RELINFO complements currently used methods by providing an information-based view as to why a prediction succeeds and fails. To demonstrate this score's capabilities, we applied RELINFO to an analysis of a large set of secondary structure predictions obtained from the first five rounds of the Critical Assessment of Structure Prediction (CASP) experiment. RELINFO is compared with two other common measures: percent correct (Q3) and secondary structure overlap (SOV). While the correlation between Q3 and RELINFO is approximately 0.85, RELINFO avoids certain disadvantages of Q3, including overestimating the quality of a prediction. The correlation between SOV and RELINFO is approximately 0.75. The valuable SOV measure unfortunately suffers from a saturation problem, and perhaps has unfairly given the general impression that secondary structure prediction has reached its limit since SOV hasn't improved much over the recent rounds of CASP. Although not a replacement for SOV, RELINFO has greater dispersion. Over the five rounds of CASP assessed here, RELINFO shows that predictions targets have been more difficult in successive CASP experiments, yet the predictions quality has continued to improve measurably over each round. In terms of information, the secondary structure prediction quality has almost doubled from CASP1 to CASP5. Therefore, as a different perspective of accuracy, RELINFO can help to improve prediction of protein secondary structure by providing a measure of difficulty as well as final quality of a prediction.
Assuntos
Biologia Computacional/métodos , Estrutura Secundária de Proteína , Proteínas/química , Algoritmos , SoftwareRESUMO
Most current template-based structure prediction methods concentrate on finding the correct backbone conformation and then packing sidechains within that backbone. Our packing-based method derives distance constraints from conserved relative packing groups (RPGs). In our refinement approach, the RPGs provide a level of resolution that restrains global topology while allowing conformational sampling. In this study, we test our template-based structure prediction method using 51 prediction units from CASP7 experiments. RPG-based constraints are able to substantially improve approximately two-thirds of starting templates. Upon deeper investigation, we find that true positive spatial constraints, especially those non-local in sequence, derived from the RPGs were important to building nearer native models. Surprisingly, the fraction of incorrect or false positive constraints does not strongly influence the quality of the final candidate. This result indicates that our RPG-based true positive constraints sample the self-consistent, cooperative interactions of the native structure. The lack of such reinforcing cooperativity explains the weaker effect of false positive constraints. Generally, these findings are encouraging indications that RPGs will improve template-based structure prediction.
Assuntos
Caspase 7/química , Modelos Moleculares , Homologia Estrutural de Proteína , Algoritmos , Simulação por Computador , Humanos , Conformação ProteicaRESUMO
In a fine-grained computational analysis of protein structure, we investigated the relationships between a residue's backbone conformations and its side-chain packing as well as conformations. To produce continuous distributions in high resolution, we ran molecular dynamics simulations over a set of protein folds (dynameome). In effect, the dynameome dataset samples not only the states well represented in the PDB but also the known states that are not well represented in the structural database. In our analysis, we characterized the mutual influence among the backbone phi,psi angles with the first side-chain torsion angles (chi(1)) and the volumes occupied by the side-chains. The dependencies of these relationships on side-chain environment and amino acids are further explored. We found that residue volumes exhibit dependency on backbone 2 degrees structure conformation: side-chains pack more densely in extended beta-sheet than in alpha-helical structures. As expected, residue volumes on the protein surface were larger than those in the interior. The first side-chain torsion angles are found to be dependent on the backbone conformations in agreement with previous studies, but the dynameome dataset provides higher resolution of rotamer preferences based on the backbone conformation. All three gauche(-), gauche(+), and trans rotamers show different patterns of phi,psi dependency, and variations in chi(1) value are skewed from their canonical values to relieve the steric strains. By demonstrating the utility of dynameomic modeling on the native state ensemble, this study reveals details of the interplay among backbone conformations, residue volumes and side-chain conformations.