Protein structure prediction from inaccurate and sparse NMR data using an enhanced genetic algorithm

BackgroundProtein structure prediction (PSP) has important applications in different fields, such as drug design, disease prediction, and so on. In protein structure prediction, there are two important issues. The first one is the design of the structure model and the second one is the design of the optimization technology. Because of the complexity of the realistic protein structure, the structure model adopted in this paper is a simplified model, which is called off-lattice AB model. After the structure model is assumed, optimization technology is needed for searching the best conformation of a protein sequence based on the assumed structure model. However, PSP is an NP-hard problem even if the simplest model is assumed. Thus, many algorithms have been developed to solve the global optimization problem. In this paper, a hybrid algorithm, which combines genetic algorithm (GA) and tabu search (TS) algorithm, is developed to complete this task.ResultsIn order to develop an efficient optimization algorithm, several improved strategies are developed for the proposed genetic tabu search algorithm. The combined use of these strategies can improve the efficiency of the algorithm. In these strategies, tabu search introduced into the crossover and mutation operators can improve the local search capability, the adoption of variable population size strategy can maintain the diversity of the population, and the ranking selection strategy can improve the possibility of an individual with low energy value entering into next generation. Experiments are performed with Fibonacci sequences and real protein sequences. Experimental results show that the lowest energy obtained by the proposed GATS algorithm is lower than that obtained by previous methods.ConclusionsThe hybrid algorithm has the advantages from both genetic algorithm and tabu search algorithm. It makes use of the advantage of multiple search points in genetic algorithm, and can overcome poor hill-climbing capability in the conventional genetic algorithm by using the flexible memory functions of TS. Compared with some previous algorithms, GATS algorithm has better performance in global optimization and can predict 3D protein structure more effectively.

Chapter 1 of this thesis provides an overview of the structure and function of proteins as well as their investigation via NMR spectroscopy. Proteins are one of the most important building blocks of life. They are part of every organism and fulfill a variety of functions including, for example, the catalysis of many chemical reactions, which is the reason various processes in an organism become possible in the first place. Proteins are made up of amino acids, which are linked together via peptide bonds during the protein biosynthesis. The 22 proteinogenic amino acids result in an almost unlimited number of sequential combinations that give proteins their versatility. The sequence of amino acids in a protein is also called the primary structure. From this, locally limited structural motifs such as α-helices and β-sheets are formed, representing the secondary structure. Further intramolecular interactions create the tertiary structure of a protein, which can be seen as the overall three-dimensional structure of the protein. Furthermore, several individual proteins can assemble into a larger complex to form a quaternary structure. There are various methods for investigating the structure of proteins, including nuclear magnetic resonance (NMR) spectroscopy. In NMR spectroscopy, atomic nuclei with a non-zero nuclear spin are polarized by a constant magnetic field and then excited using radio frequency pulses. The emitted resonance frequencies are detected in the form of a free induction decay (FID) and converted into an NMR spectrum by a Fourier transformation. Many NMR experiments are available for researching proteins. By default, a two-dimensional 1H15N-HSQC spectrum is measured, which provides signals for each N-H bond in a protein meaning that one signal is obtained for each N-H in the backbone of the amino acids. Using 3D NMR experiments, these signals can then be assigned to the respective amino acids in the protein sequence. This makes it possible to investigate the dynamics of proteins as well as the formation of complexes with a ligand or another biomacromolecule if they interact. Further 3D NMR experiments also make it possible to assign the side-chains of the protein. With the help of the determined resonance assignments and spatial information from various NOESY spectra, a protein structure can then be calculated using a suitable computer program. Recently, the prediction of protein structures using AI (artificial intelligence)-supported programs has also become more prevalent. In particular, the artificial neural network AlphaFold has been extremely successful in predicting correct structures without any experimental data provided by the user. Chapter 2 of this thesis introduces the research field of small proteins. These are proteins with a length of usually less than 100 amino acids, although definitions in literature vary and upper limits of, for example, 50 or 200 amino acids exist as well. Various other names are also used to describe small proteins, such as μ-proteins or sORF-encoded peptides. The field of small proteins is still relatively new, as they were overlooked in research for many years. One reason for this was the lower limit of 100 codons that was set for the identification of open reading frames (ORFs) in genome projects. This served to exclude numerous biologically irrelevant ORFs from genome annotation. However, this meant that functional ORFs coding for small proteins were also ignored. Another reason was that experimental methods were not designed for the analysis of small proteins, meaning they could often not be detected at all. However, new analytical methods such as ribosome profiling and a generally increased focus on the detection of small proteins have led to their detection in all domains of life. In the research group of Prof. Dr. Schwalbe, small proteins of various prokaryotes were structurally characterized as part of a special priority program by the DFG. NMR spectroscopy was primarily used for this purpose, as it is the method of choice for the investigation of small proteins. The results of this research were published in several articles, two of which are part of this dissertation.

Background: The Mycobacterium tuberculosis complex (MTBC) bacteria include the slowly growing, host-associated bacteria Mycobacterium tuberculosis, Mycobacterium Bovis, Mycobacterium microti, Mycobacterium africanum, Mycobacterium pinnipedii. Aim: Comparative Functional Genomics Studies for understanding the Hypothetical Proteins in Mycobacterium tuberculosis variant microti 12. Objective: A computational genomics study was performed to understand the 247 hypothetical protein genes. Functional annotation of virtual proteins was performed on different servers to maximize confidence level. Methods: Sequence Retrieval. The whole genome sequences for the Mycobacterium tuberculosis micro variant 12 were retrieved from the KEGG database ( http://www.genome.jp/kegg/) and were used for screening 247 hypothetical proteins (Fig. 1 ). Functional Annotation and Sub-cellular localization. The Mycobacterium tuberculosis micro variant 12 hypothetical proteins were screened and sorted out from the genome and were individually analyzed for the presence of conserved functional domains by using computational biology tools like CDD-BLAST ( https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) ;Pfam ( http://pfam.xfam.org/ncbiseq/398365647); The subcellular localization of hypothetical proteins was determined by CELLO2GO ( http://cello.life.nctu.edu.tw). These web tools can search the defined conserved domains in the sequences available in the online servers or databases and assist in the classification of proteins in the appropriate families. Protein Structure Prediction. The in-silico structure predictions of the hypothetical protein sequences showing functional properties were carried out by using the PS2 Protein Structure Prediction Server ( http://www.ps2.life.nctu.edu.tw/). The online server helps to generate the 3D structures of the hypothetical proteins. The server accepts the sequences in FASTA format as a query to generate resultant proteins 3D structures. The structure determination is completely based on the conserved template regions detected during functional annotations. Protein-protein interaction through String database: The interaction of each hypothetical protein analyzed for functional characteristics was subjected to a protein-protein interaction server for the prediction of a possible functional role in interaction amongst the available known proteins ( https://string-db.org/). This information can help us to further validated the functional role of such hypothetical proteins and their possible role in the Mycobacterium Tuberculosis micro variant. Protein secondary structure prediction through JPred4: The secondary structure prediction of all the hypothetical proteins was determined through JPred4 ( http://www.compbio.dundee.ac.uk/jpred4/index.html) and served to identify the available secondary structures in the unknown hypothetical protein sequences. These further help us to understand the available templates in the uncharacterized protein sequences for the prediction of novel functions associated with these proteins. The predictions were further characterized by the Phyre2 server for structural modeling and prediction of templates based on comparative analysis based on conserved domains. Protein modeling, prediction, and analysis through Phyre2. The hypothetical proteins which were identified to have functional properties were further characterized by the Phyre2 server ( http://www.sbg.bio.ic.ac.uk/phyre2) for structural modeling and prediction of templates based on comparative analysis based on conserved domains. Results: A computational genomics study was performed to understand the 247 hypothetical protein genes Functional annotation of virtual proteins, and was performed on different servers to maximize confidence level. The functional prediction was performed by CDD-Blast and Pfam. The gene sequences of proteins have probably been successfully functionally annotated, characterized, and their subcellular localization and 3-D structural predictions have been predicted computationally. Online automated bioinformatics tools such as CDD-Blast, Pfam, CELLO2GO and PS2-Server were used for the structural and functional characterization of screened hypothetical proteins. The structure, function, and subcellular localization of a hypothetical protein from Mycobacterium tuberculosis variant microti 12 have been obtained and presented (Fig. 2 ). Also, the three-dimensional structure generated after using the template with the highest score was displayed as the template ID in the structure column of the respective hypothetical protein. However, as systems biology denies hypothetical protein functions, the structures of such proteins can be tested through biological processes and experiments, making them suitable for understanding their role in the life cycle, pathogenesis, and drug development. We can further explore these predictive possibilities in pharmaceuticals, and other clinically relevant studies. This study by HP helped find structure-function relationships in Mycobacterium tuberculosis variant microti 12 using a variety of bioinformatics tools. The string database made predictions about protein-protein interactions and the template helped us predict a hypothetical protein structure and even helped us find its 3D protein structure. Protein profiling can be performed on structures retrieved from these servers. This is useful for proteomics studies, including protein-protein interactions, protein expression of specific hypothetical proteins, and post-translational modifications of protein-coding genes. Further understanding of these hypothetical proteins can help us to know more about the Mycobacterium tuberculosis complex (MTBC) and may assist in Drugs and inhibitors against different pathogens within this complex. Conclusion: The all-inclusive bioinformatic study has helped to functionally elucidate 247 hypothetical proteins, which have resulted and made it easier to understand many functional proteins available in the Mycobacterium tuberculosis micro variant 12. The subcellular localization of the 247 sorted hypothetical proteins was also carried & which further helped us understand the localization of identified enzymes or proteins. We have successfully characterized the 247 unknown proteins of hypothetical protein sequences from Mycobacterium tuberculosis micro variant 12 to validate their structure and functions of the gene products. These predicted functions and three-dimensional structures may lead to establishing their role in the life cycle of the bacterium. This computationally generated data can also be further used for developing new protocols for new vaccines against Mycobacterium tuberculosis micro variant 12 that are essential for preventing infection, diseases, and transmission. This complete result of Hypothetical Protein is needed for further studies of the whole genomic of the Mycobacterium Tuberculosis micro variant 12 for their function interpretation which further help in the understanding of its functions as well as structure. Moreover, this interpretation would help us to study the evolution of Mycobacterium Tuberculosis micro variant 12 which further helps in the process of discovering the drugs to inhibit the causes of diseases.

We present a fast and simple protocol to obtain moderate-resolution backbone structures of helical proteins. This approach utilizes a combination of sparse backbone NMR data (residual dipolar couplings and paramagnetic relaxation enhancements) or EPR data with a residue-based force field and Monte Carlo/simulated annealing protocol to explore the folding energy landscape of helical proteins. By using only backbone NMR data, which are relatively easy to collect and analyze, and strategically placed spin relaxation probes, we show that it is possible to obtain protein structures with correct helical topology and backbone RMS deviations well below 4 Å. This approach offers promising alternatives for the structural determination of proteins in which nuclear Overha-user effect data are difficult or impossible to assign and produces initial models that will speed up the high-resolution structure determination by NMR spectroscopy.

Proteins are essential parts of our life and participate in virtually every process within a cell. The understanding of protein structures is vital to determine the function of a protein. Protein structure prediction (PSP) from amino acid sequence is one of the high focus problems in bioinformatics today. This is due to the fact that the biological function of the protein is determined by its three dimensional structure. Thus, protein structure prediction is a fundamental area of computational biology. Its importance is intensed by large amounts of sequence data coming from PDB (Protein Data Bank) and the fact that experimentally methods such as X-ray crystallography or Nuclear Magnetic Resonance (NMR)which are used to determining protein structures remains very expensive and time consuming. For minimizing the time, computational methods are used for protein folding and structure prediction problem. In this paper results of protein p53 are discussed.

One of the main research problems in Structural Bioinformatics is related to the prediction of three-dimensional structures (3-D) of polypeptides or proteins. The rate at which amino acid sequences are identified is increasing faster than the 3-D protein structure determination by experimental methods. Computational prediction methods have been developed during the last years, but the problem still remains challenging because of the complexity and high dimensionality of a protein conformational search space. In this article we present a hybrid genetic algorithm for the Protein Structure Prediction (PSP) Problem. A genetic algorithm is combined with a structured population, and it is hybridized with a path-relinking procedure that helps the algorithm to scape from local minima. We perform a set of experiments and show that the proposed hybrid genetic algorithm is effective in finding good quality solutions for the PSP Problem.

Structural genomics provides an important approach for characterizing and understanding systems biology. As a step toward better integrating protein three-dimensional (3D) structural information in cancer systems biology, we have constructed a Human Cancer Pathway Protein Interaction Network (HCPIN) by analysis of several classical cancer-associated signaling pathways and their physical protein-protein interactions. Many well known cancer-associated proteins play central roles as "hubs" or "bottlenecks" in the HCPIN. At least half of HCPIN proteins are either directly associated with or interact with multiple signaling pathways. Although some 45% of residues in these proteins are in sequence segments that meet criteria sufficient for approximate homology modeling (Basic Local Alignment Search Tool (BLAST) E-value <10(-6)), only approximately 20% of residues in these proteins are structurally covered using high accuracy homology modeling criteria (i.e. BLAST E-value <10(-6) and at least 80% sequence identity) or by actual experimental structures. The HCPIN Website provides a comprehensive description of this biomedically important multipathway network together with experimental and homology models of HCPIN proteins useful for cancer biology research. To complement and enrich cancer systems biology, the Northeast Structural Genomics Consortium is targeting >1000 human proteins and protein domains from the HCPIN for sample production and 3D structure determination. The long range goal of this effort is to provide a comprehensive 3D structure-function database for human cancer-associated proteins and protein complexes in the context of their interaction networks. The network-based target selection (BioNet) approach described here is an example of a general strategy for targeting co-functioning proteins by structural genomics projects.

Protein structure prediction constrained by solution X-ray scattering data and structural homology identification

Databases are the computerized platform where information is stored and can be retrieved easily by public users. Biological databases are the repositories of biological data. These biological data libraries contain facts and figures related to various disciplines of research including genomics, proteomics, microarray technology, metabolomics and phylogenetics. By using biological databases, a broad collection of essential biological information can be exploited ranging from function, structure and localization of gene, clinical consequences of mutation to similarity index among biological sequences and structures. Nowadays, different kinds of biological databases are available on the web. The present write up focuses on biological databases and bioinformatics tools for protein structure analysis. This review also aims to elaborate the searching schemes, available in different structural databases. The wide variety of different levels and types of information content related to 3D protein structures are available on web-based databases. Regarding the biological functions and 3D structures of various proteins, these databases provide a huge range of useful links, schematic diagrams as well as strategies for detailed analysis of proteins and other macromolecules structures. 3D structural illustration of proteins stored in structural databases is determined and visualized by X-ray crystallography, electron microscopy and NMR spectroscopy. On regular basis, a large number of protein structures are submitted by structural biologists, updated and curated by subject experts. Most familiar biological databases that store 3D protein and other macromolecules structures include, PDB, 3D Genomics, CATH, &amp; SCOP. These databases contain valuable information of overall protein structures, domains and motif structures, protein-protein complex systems and complex of protein with other biomolecules.

Covalent labeling mass spectrometry experiments are growing in popularity and provide important information regarding protein structure. Information obtained from these experiments correlates with residue solvent exposure within the protein in solution. However, it is impossible to determine protein structure from covalent labeling data alone. Incorporation of sparse covalent labeling data into the protein structure prediction software Rosetta has been shown to improve protein tertiary structure prediction. Here, covalent labeling techniques were analyzed computationally to provide insight into what labeling data is needed to optimize tertiary protein structure prediction in Rosetta. We have successfully implemented a new scoring functionality that provides improved predictions. We developed two new covalent labeling based score terms that use a "cone"-based neighbor count to quantify the relative solvent exposure of each amino acid. To test our method, we used a set of 20 proteins with structures deposited in the Protein Data Bank. Decoy model sets were generated for each of these 20 proteins, and the normalized covalent labeling score versus RMSD distributions were evaluated. On the basis of these distributions, we have determined an optimal subset of residues to use when performing covalent labeling experiments in order to maximize the structure prediction capabilities of the covalent labeling data. We also investigated how much false negative and false positive data can be tolerated without meaningfully impacting protein structure prediction. Using these new covalent labeling score terms, protein models were rescored and the resulting models improved by 3.9 Å RMSD on average. New models were also generated using Rosetta's AbinitioRelax program under the guidance of covalent labeling information, and improvement in model quality was observed.

Protein-Peptide Complex Prediction through Fragment Interaction Patterns

Predicting the 3D native conformation of a protein given the amino acid sequence is known as protein structure prediction (PSP) problem and is one of the greatest challenges of computational biology. The current work uses a real valued Genetic Algorithm (GA), a powerful variate of GA to simulate the PSP problem. This algorithm consists of basic evolutionary operators and a fitness vector. The fitness vector is designed by combining a set of knowledge based biophysical filters viz. persistence length, radius of gyration, packing fraction, hydrophobicity ratio and irregularity index of phi and psi. This vector converts all these biophysical measures into a real value by using specific weights or factors. The algorithm has been validated on six known globular proteins, with their length varying from 17-61 residues and total number of helices and strands in the range of 2-4. For all the test protein the algorithm converges rapidly and the converged structure shows a backbone RMSD (root mean square deviation) of 3-6A as compared to the native structure.

Article Free Access Share on Protein structure prediction and potential energy landscape analysis using continuous global minimization Authors: Ken A. Dill Department of Pharmaceutical Chemistry, University of California at San Francisco, San Francisco, CA Department of Pharmaceutical Chemistry, University of California at San Francisco, San Francisco, CAView Profile , Andrew T. Phillips Computer Science Department, United States Naval Academy, Annapolis, MD Computer Science Department, United States Naval Academy, Annapolis, MDView Profile , J. Ben Rosen Computer Science and Engineering Department, University of California at San Diego, San Diego, CA Computer Science and Engineering Department, University of California at San Diego, San Diego, CAView Profile Authors Info & Claims RECOMB '97: Proceedings of the first annual international conference on Computational molecular biologyJanuary 1997 Pages 109–117https://doi.org/10.1145/267521.267537Online:19 January 1997Publication History 2citation673DownloadsMetricsTotal Citations2Total Downloads673Last 12 Months2Last 6 weeks1 Get Citation AlertsNew Citation Alert added!This alert has been successfully added and will be sent to:You will be notified whenever a record that you have chosen has been cited.To manage your alert preferences, click on the button below.Manage my AlertsNew Citation Alert!Please log in to your account Save to BinderSave to BinderCreate a New BinderNameCancelCreateExport CitationPublisher SiteeReaderPDF

BackgroundAssociation and cooperation among structural domains play an important role in protein function and drug design. Despite remarkable advancements in highly accurate single-domain protein structure prediction through the collaborative efforts of the community using deep learning, challenges still exist in predicting multi-domain protein structures when the evolutionary signal for a given domain pair is weak or the protein structure is large.ResultsTo alleviate the above challenges, we proposed M-DeepAssembly, a protocol based on multi-objective protein conformation sampling algorithm for multi-domain protein structure prediction. Firstly, the inter-domain interactions and full-length sequence distance features are extracted through DeepAssembly and AlphaFold2, respectively. Secondly, subject to these features, we constructed a multi-objective energy model and designed a sampling algorithm for exploring and exploiting conformational space to generate ensembles. Finally, the output protein structure was selected from the ensembles using our in-house developed model quality assessment algorithm. On the test set of 164 multi-domain proteins, the results show that the average TM-score of M-DeepAssembly is 15.4% and 2.0% higher than AlphaFold2 and DeepAssembly, respectively. It is worth noting that there are models with higher accuracy in ensembles, achieving an improvement of 20.3% and 6.4% relative to the two baseline methods, although these models were not selected. Furthermore, when compared to the prediction results of AlphaFold2 for CASP15 multi-domain targets, M-DeepAssembly demonstrates certain performance advantages.ConclusionsM-DeepAssembly provides a distinctive multi-domain protein assembly algorithm, which can alleviate the current challenges of weak evolutionary signals and large structures to some extent by forming diverse ensembles using multi-objective protein conformation sampling algorithm. The proposed method contributes to exploring the functions of multi-domain proteins, especially providing new insights into targets with multiple conformational states.

In this study, we used the Hydrophobic-Polar (HP) two-dimensional square and three-dimensional cubic lattice models for the problem of Protein structure prediction (PSP). This kind of lattice reduces computational time and calculations, and the conformational space from 9n to 4n−2, even in this context, it is an appropriate benchmark problem for genetic algorithms. The contributions of the paper consist in: (1) implementation of a high-performing novel Genetic algorithm (GA); instead of considering only the self- avoiding walk (SAW) conformations approached in other work, we decided to allow any conformation to appear in the population at all stages of the proposed all conformations biomimetic Genetic Algorithm (ACGA). This increases the probability of achieving good native conformations (SAW ones), with the lowest energy that has the minimum number of contacts in the lattice, the one with the maximum topological neighbors. In addition to classical crossover and mutation operators, (2) we introduced specific translation operators for these two operations. We have proposed and implemented an HP Protein Visualizer tool which offers interpretability, a hybrid approach in that the visualizer give some insight to the algorithm, that analyse and optimise protein structures HP model. The program resulted based on performed research, provides a molecular modeling tool for studying protein folding using technologies such as React, Next.js, and Three.js for 3D rendering, and includes optimization algorithms to simulate protein folding.