Gotowa bibliografia na temat „Protein Structure Models”

The goal of the following thesis research was to develop a systematic approach for the refinement of low-resolution protein models, as a part of the protein structure prediction procedure. Significant progress has been made in the field of protein structure prediction and the contemporary methods are able to assemble correct topology for a large fraction of protein domains. But such approximate models are often not detailed enough for some important applications, including studies of reaction mechanisms, functional annotation, drug design or virtual ligand screening. The development of a method that could bring those structures closer to the native is then of great importance. The minimal requirements for a potential that can refine protein structures is the existence of a correlation between the energy with native similarity and the scoring of the native structure as being lowest in energy. Extensive tests of the contemporary all-atom physics-based force fields were conducted to assess their applicability for refinement. The tests revealed flatness of such potentials and enabled the identification of the key problems in the current approaches. Guided by these results, the optimization of the AMBER (ff03) force field was performed that aimed at creating a funnel shape of the potential, with the native structure at the global minimum. Such shape should facilitate the conformational search during refinement and drive it towards the native conformation. Adjusting the relative weights of particular energy components, and adding an explicit hydrogen bond potential significantly improved the average correlation coefficient of the energy with native similarity (from 0.25 for the original ff03 potential to 0.65 for the optimized force field). The fraction of proteins for which the native structure had lowest energy increased from 0.22 to 0.90. The new, optimized potential was subsequently used to refine protein models of various native-similarity. The test employed 47 proteins and 100 decoy structures per protein. When the lowest energy structure from each trajectory was compared with the starting decoy, we observed structural improvement for 70% of the models on average. Such an unprecedented result of a systematic refinement is extremely promising in the context of high-resolution structure prediction.

<p>Evolutionary algorithms (EAs) have been shown to be competent at solving complex, multimodal optimisation problems in applications where the search space is large and badly understood. EAs are therefore among the most promising classes of algorithms for solving the Protein Structure Prediction Problem (PSPP). The PSPP is how to derive the 3D-structure of a protein given only its sequence of amino acids. This dissertation defines, evaluates and shows limitations of simplified models for solving the PSPP. These simplified models are off-lattice extensions to the lattice HP model which has been proposed and is claimed to possess some of the properties of real protein folding such as the formation of a hydrophobic core. Lattice models usually model a protein at the amino acid level of detail, use simple energy calculations and are used mainly for search algorithm development. Off-lattice models usually model the protein at the atomic level of detail, use more complex energy calculations and may be used for comparison with real proteins. The idea is to combine the fast energy calculations of lattice models with the increased spatial possibilities of an off-lattice environment allowing for comparison with real protein structures. A hypothesis is presented which claims that a simplified off-lattice model which considers other amino acid properties apart from hydrophobicity will yield simulated structures with lower Root Mean Square Deviation (RMSD) to the native fold than a model only considering hydrophobicity. The hypothesis holds for four of five tested short proteins with a maximum of 46 residues. Best average RMSD for any model tested is above 6Å, i.e. too high for useful structure prediction and excludes significant resemblance between native and simulated structure. Hence, the tested models do not contain the necessary biological information to capture the complex interactions of real protein folding. It is also shown that the EA itself is competent and can produce near-native structures if given a suitable evaluation function. Hence, EAs are useful for eventually solving the PSPP.</p>

Bioinformatics is a fast-developing field, which makes use of computational methods to analyse and structure biological data. An important branch of bioinformatics is structure and function prediction of proteins, which is often based on finding relationships to already characterized proteins. It is known that two proteins with very similar sequences also share the same 3D structure. However, there are many proteins with similar structures that have no clear sequence similarity, which make it difficult to find these relationships. In this thesis, two methods for annotating protein domains are presented, one aiming at assigning the correct domain family or families to a protein sequence, and the other aiming at fold recognition. Both methods use hidden Markov models (HMMs) to find related proteins, and they both exploit the fact that structure is more conserved than sequence, but in two different ways. Most of the research presented in the thesis focuses on the structure-anchored HMMs, saHMMs. For each domain family, an saHMM is constructed from a multiple structure alignment of carefully selected representative domains, the saHMM-members. These saHMM-members are collected in the so called "midnight ASTRAL set", and are chosen so that all saHMM-members within the same family have mutual sequence identities below a threshold of about 20%. In order to construct the midnight ASTRAL set and the saHMMs, a pipe-line of software tools are developed. The saHMMs are shown to be able to detect the correct family relationships at very high accuracy, and perform better than the standard tool Pfam in assigning the correct domain families to new domain sequences. We also introduce the FI-score, which is used to measure the performance of the saHMMs, in order to select the optimal model for each domain family. The saHMMs are made available for searching through the FISH server, and can be used for assigning family relationships to protein sequences. The other approach presented in the thesis is secondary structure HMMs (ssHMMs). These HMMs are designed to use both the sequence and the predicted secondary structure of a query protein when scoring it against the model. A rigorous benchmark is used, which shows that HMMs made from multiple sequences result in better fold recognition than those based on single sequences. Adding secondary structure information to the HMMs improves the ability of fold recognition further, both when using true and predicted secondary structures for the query sequence.<br>Bioinformatik är ett område där datavetenskapliga och statistiska metoder används för att analysera och strukturera biologiska data. Ett viktigt område inom bioinformatiken försöker förutsäga vilken tredimensionell struktur och funktion ett protein har, utifrån dess aminosyrasekvens och/eller likheter med andra, redan karaktäriserade, proteiner. Det är känt att två proteiner med likande aminosyrasekvenser också har liknande tredimensionella strukturer. Att två proteiner har liknande strukturer behöver dock inte betyda att deras sekvenser är lika, vilket kan göra det svårt att hitta strukturella likheter utifrån ett proteins aminosyrasekvens. Den här avhandlingen beskriver två metoder för att hitta likheter mellan proteiner, den ena med fokus på att bestämma vilken familj av proteindomäner, med känd 3D-struktur, en given sekvens tillhör, medan den andra försöker förutsäga ett proteins veckning, d.v.s. ge en grov bild av proteinets struktur. Båda metoderna använder s.k. dolda Markov modeller (hidden Markov models, HMMer), en statistisk metod som bland annat kan användas för att beskriva proteinfamiljer. Med hjälp en HMM kan man förutsäga om en viss proteinsekvens tillhör den familj modellen representerar. Båda metoderna använder också strukturinformation för att öka modellernas förmåga att känna igen besläktade sekvenser, men på olika sätt. Det mesta av arbetet i avhandlingen handlar om strukturellt förankrade HMMer (structure-anchored HMMs, saHMMer). För att bygga saHMMerna används strukturbaserade sekvensöverlagringar, vilka genereras utifrån hur proteindomänerna kan läggas på varandra i rymden, snarare än utifrån vilka aminosyror som ingår i deras sekvenser. I varje proteinfamilj används bara ett särskilt, representativt urval av domäner. Dessa är valda så att då sekvenserna jämförs parvis, finns det inget par inom familjen med högre sekvensidentitet än ca 20%. Detta urval görs för att få så stor spridning som möjligt på sekvenserna inom familjen. En programvaruserie har utvecklats för att välja ut representanter för varje familj och sedan bygga saHMMer baserade på dessa. Det visar sig att saHMMerna kan hitta rätt familj till en hög andel av de testade sekvenserna, med nästan inga fel. De är också bättre än den ofta använda metoden Pfam på att hitta rätt familj till helt nya proteinsekvenser. saHMMerna finns tillgängliga genom FISH-servern, vilken alla kan använda via Internet för att hitta vilken familj ett intressant protein kan tillhöra. Den andra metoden som presenteras i avhandlingen är sekundärstruktur-HMMer, ssHMMer, vilka är byggda från vanliga multipla sekvensöverlagringar, men också från information om vilka sekundärstrukturer proteinsekvenserna i familjen har. När en proteinsekvens jämförs med ssHMMen används en förutsägelse om sekundärstrukturen, och den beräknade sannolikheten att sekvensen tillhör familjen kommer att baseras både på sekvensen av aminosyror och på sekundärstrukturen. Vid en jämförelse visar det sig att HMMer baserade på flera sekvenser är bättre än sådana baserade på endast en sekvens, när det gäller att hitta rätt veckning för en proteinsekvens. HMMerna blir ännu bättre om man också tar hänsyn till sekundärstrukturen, både då den riktiga sekundärstrukturen används och då man använder en teoretiskt förutsagd.<br>Jeanette Hargbo.

Here I investigate the protein structure refinement problem for homology-based protein structure models. The refinement problem has been identified as a major bottleneck in the structure prediction process and inhibits the goal of producing high-resolution experimental quality structures for target protein sequences. This thesis is composed of three investigations into aspects of template-based modelling and refinement. In the primary investigation, empirical evidence is provided to support the hypothesis that using multiple template-based structures to model a target sequence can improve the quality of the prediction over that obtained solely by using the single best prediction. A multi-objective genetic algorithm is used to optimize protein structure models by using the structural information from a set of predictions, guided by various objective functions. The effect of multi-objective optimization on model quality is examined. A benchmark of energy functions and model quality assessment methods is performed in the context of automated homology modelling to assess the ability of these methods at discriminating nearer-native structures from a set of predictions. These model quality assessment methods were unable to significantly improve the ranking of threading- based prediction methods though some model quality assessment methods improved model selection for methods which use sequence information alone. The results suggest that structural informational can provide valuable information for distinguishing better models where only sequence information has been used for modelling. The suitability of these energy functions for high-resolution refinement is discussed. Finally, a stochastic optimization algorithm is developed for refining homology-based protein structure models using evolutionary algorithms. This approach uses multiple structural model inputs, conformational sampling operators, and objective functions for guiding a search through conformational space. Single- and multi-objective genetic variants are applied to homology model predictions for 35 target proteins. The refinement results are discussed and the performance of both algorithmic variants compared and contrasted.

The aim of this work was to predict protein secondary structure using neural network models. Initially a Hopfield network was used but abandoned in favour of a layered network trained using the back propagation algorithm. In the early stages of this work an exploration of the many different approaches to this problem was undertaken. These included attempts to predict boundaries between secondary structures, the secondary structures of individual residues, and the secondary structures of sequences wholly within a particular secondary structure. Results indicated the latter to be the best approach to continue with. In addition two coding schemes were investigated: a coding scheme based on occurrences of pairs of residues and one based on the positions of residues. It was found that this positional coding scheme was the natural coding scheme for this problem. On segments of whole alpha-helix and whole non-alpha-helix 10 residues in length a prediction success of around 80% with a correlation coefficient of 0.52 was achieved with the positional coding scheme. On whole proteins, where predictions are made for individual residues, alpha-helix prediction drops to 73% with a correlation coefficient of 0.34. The relative predictability of alpha-helices of above and below average accessibility was also investigated showing that those of above average accessibility are more predictable than those with below average accessibility. The main body of this work concerns the apparent limit of predictability of alpha-helices. It was found that test set prediction did not depend on the number of hidden nodes. In fact, a single layer network performed as well as those with hidden nodes showing that the probolem is basically linearly separable. In addition, prediction success plateaus well below that of perfect prediction success. During training, test set prediction is seen to peak. The decrease in prediction success was found to be due to non-alpha-helix sequences that the network was unable to distinguish from real alpha-helix sequences. These regions of non-alpha-helix were shown to occur adjacent to actual alpha-helices with statistical significance. It is proposed that potential alpha-helices are disrupted by global constraints during the formation of tertiary structure. The effect of window size was also investigated as was beta-sheet prediction, but this was found to be limited by the small number of examples available with our approach. However, its distribution in the input space in relation to alpha-helix and coil was determined.

The three-dimensional structure of a protein, formed as a result of amino-acid sequences folding into compact domains, is regarded as a key factor in its biological function. How and why proteins fold into specific topologies, remain the key focus of scientific research in the field of biophysics. By stripping down complex reactions down to the most basic elements, biophysicists aim to develop simplified models for biological phenomena such as antibody discrimination, viral fusion or self-assembly. Focusing on small model peptide systems, rather than the full proteins from which they were derived, was hoped to result in accurate structural measurements and provide a more transparent comparison between simulation and experiment. The aim of this research was therefore to investigate how accurate these models were when compared against experiment. Furthermore, while breaking down the complex biological phenomena into simple models, there was also a conscious effort to ensure that the models were representative of real biological systems, and a major focus was therefore aimed at determining whether any meaningful biomedical insight may be extrapolated from such models. Peptides found in hormones (human chorionic gonadotropin, luteinizing hormone), viruses (HIV) and amyloid diseases (transthyretin) were selected in order to probe a variety of questions in relation to the aforementioned biological phenomena. Namely, how the primary sequence influenced the three-dimensional structure (and thus its biological function), how its environment could influence such a confirmation, and how these systems aggregated. This doctoral study has made use of a combination of computer simulations and experimental techniques to investigate a selection of biologically relevant peptides; utilising classical atomistic molecular dynamics (MD) simulations to characterise the free-energy landscapes of the chosen peptides, and compare these findings with the secondary structure content predicted by spectroscopic methods such as circular dichroism and infrared spectroscopy. The peptide systems studied within, were found to be characterised by rugged free-energy landscapes unlike their protein counterparts (defined by singular, deep minima). Furthermore, these landscapes were found to be highly plastic and sensitive to changes in the local environment.