Template-based Modeling Methods Research Articles

Structural Perspective of NR4A Nuclear Receptor Family and Their Potential Endogenous Ligands.

There are 48 nuclear receptors in the human genome, and many members of this superfamily have been implicated in human diseases. The NR4A nuclear receptor family consisting of three members, NR4A1, NR4A2, and NR4A3 (formerly annotated as Nur77, Nurr1, and NOR1, respectively), are still orphan receptors but exert pathological effects on immune-related and neurological diseases. We previously reported that prostaglandin A1 (PGA1) and prostaglandin A2 (PGA2) are potent activators of NR4A3, which bind directly to the ligand-binding domain (LBD) of the receptor. Recently, the co-crystallographic structures of NR4A2-LBD bound to PGA1 and PGA2 were reported, followed by reports of the neuroprotective effects of these possible endogenous ligands in mouse models of Parkinson's disease. Based on these structures, we modeled the binding structures of the other two members (NR4A1 and NR4A3) with these potential endogenous ligands using a template-based modeling method, and reviewed the similarity and diversity of ligand-binding mechanisms in the nuclear receptor family.

Towards a Transferable Modeling Method of the Knee to Distinguish Between Future Healthy Joints from Osteoarthritic Joints: Data from the Osteoarthritis Initiative.

Computational models can be used to predict the onset and progression of knee osteoarthritis. Ensuring the transferability of these approaches among computational frameworks is urgent for their reliability. In this work, we assessed the transferability of a template-based modeling strategy, based on the finite element (FE) method, by implementing it on two different FE softwares and comparing their results and conclusions. For that, we simulated the knee joint cartilage biomechanics of 154 knees using healthy baseline conditions and predicted the degeneration that occurred after 8 years of follow-up. For comparisons, we grouped the knees using their Kellgren-Lawrence grade at the 8-year follow-up time and the simulated volume of cartilage tissue that exceeded age-dependent thresholds of maximum principal stress. We considered the medial compartment of the knee in the FE models and used ABAQUS and FEBio FE softwares for simulations. The two FE softwares detected different volumes of overstressed tissue in corresponding knee samples (p < 0.01). However, both programs correctly distinguished between the joints that remained healthy and those that developed severe osteoarthritis after the follow-up (AUC = 0.73). These results indicate that different software implementations of a template-based modeling method similarly classify future knee osteoarthritis grades, motivating further evaluations using simpler cartilage constitutive models and additional studies on the reproducibility of these modeling strategies.

MULTICOM2 open-source protein structure prediction system powered by deep learning and distance prediction

Protein structure prediction is an important problem in bioinformatics and has been studied for decades. However, there are still few open-source comprehensive protein structure prediction packages publicly available in the field. In this paper, we present our latest open-source protein tertiary structure prediction system—MULTICOM2, an integration of template-based modeling (TBM) and template-free modeling (FM) methods. The template-based modeling uses sequence alignment tools with deep multiple sequence alignments to search for structural templates, which are much faster and more accurate than MULTICOM1. The template-free (ab initio or de novo) modeling uses the inter-residue distances predicted by DeepDist to reconstruct tertiary structure models without using any known structure as template. In the blind CASP14 experiment, the average TM-score of the models predicted by our server predictor based on the MULTICOM2 system is 0.720 for 58 TBM (regular) domains and 0.514 for 38 FM and FM/TBM (hard) domains, indicating that MULTICOM2 is capable of predicting good tertiary structures across the board. It can predict the correct fold for 76 CASP14 domains (95% regular domains and 55% hard domains) if only one prediction is made for a domain. The success rate is increased to 3% for both regular and hard domains if five predictions are made per domain. Moreover, the prediction accuracy of the pure template-free structure modeling method on both TBM and FM targets is very close to the combination of template-based and template-free modeling methods. This demonstrates that the distance-based template-free modeling method powered by deep learning can largely replace the traditional template-based modeling method even on TBM targets that TBM methods used to dominate and therefore provides a uniform structure modeling approach to any protein. Finally, on the 38 CASP14 FM and FM/TBM hard domains, MULTICOM2 server predictors (MULTICOM-HYBRID, MULTICOM-DEEP, MULTICOM-DIST) were ranked among the top 20 automated server predictors in the CASP14 experiment. After combining multiple predictors from the same research group as one entry, MULTICOM-HYBRID was ranked no. 5. The source code of MULTICOM2 is freely available at https://github.com/multicom-toolbox/multicom/tree/multicom_v2.0.

Sequence alignment generation using intermediate sequence search for homology modeling

Protein tertiary structure is important information in various areas of biological research, however, the experimental cost associated with structure determination is high, and computational prediction methods have been developed to facilitate a more economical approach. Currently, template-based modeling methods are considered to be the most practical because the resulting predicted structures are often accurate, provided an appropriate template protein is available. During the first stage of template-based modeling, sensitive homology detection is essential for accurate structure prediction. However, sufficient structural models cannot always be obtained due to a lack of quality in the sequence alignment generated by a homology detection program. Therefore, an automated method that detects remote homologs accurately and generates appropriate alignments for accurate structure prediction is needed. In this paper, we propose an algorithm for suitable alignment generation using an intermediate sequence search for use with template-based modeling. We used intermediate sequence search for remote homology detection and intermediate sequences for alignment generation of remote homologs. We then evaluated the proposed method by comparing the sensitivity and selectivity of homology detection. Furthermore, based on the accuracy of the predicted structure model, we verify the accuracy of the alignments generated by our method. We demonstrate that our method generates more appropriate alignments for template-based modeling, especially for remote homologs. All source codes are available at https://github.com/shuichiro-makigaki/agora.

Utilizing Atlas-Based Modeling to Predict Knee Joint Cartilage Degeneration: Data from the Osteoarthritis Initiative

Currently, there are no clinically available tools or applications which could predict osteoarthritis development. Some computational models have been presented to simulate cartilage degeneration, but they are not clinically feasible due to time required to build subject-specific knee models. Therefore, the objective of this study was to develop a template-based modeling method for rapid prediction of knee joint cartilage degeneration. Knee joint models for 21 subjects were constructed with two different template approaches (multiple templates and one template) based on the MRI data. Geometries were also generated by manual segmentation. Evaluated volumes of cartilage degeneration for each subject, as assessed with the degeneration algorithm, were compared with experimentally observed 4 year follow up Kellgren-Lawrence (KL) grades. Furthermore, the effect of meniscus was tested by generating models with subject-specific meniscal supporting forces and those with the average meniscal supporting force from all models. All tested models were able to predict most severe cartilage degeneration to those subjects who had the highest KL grade after 4 year follow up. Surprisingly, in terms of statistical significance, the best result was obtained with one template approach and average meniscal support. This model was fully able to categorize all subjects to their experimentally defined groups (KL0, KL2 and KL3) based on the 4 year follow-up data. The results suggest that a template- or population-based approach, which is much faster than fully subject-specific, could be applied as a clinical prediction tool for osteoarthritis.

Protein alignment based on higher order conditional random fields for template-based modeling.

The query-template alignment of proteins is one of the most critical steps of template-based modeling methods used to predict the 3D structure of a query protein. This alignment can be interpreted as a temporal classification or structured prediction task and first order Conditional Random Fields have been proposed for protein alignment and proven to be rather successful. Some other popular structured prediction problems, such as speech or image classification, have gained from the use of higher order Conditional Random Fields due to the well known higher order correlations that exist between their labels and features. In this paper, we propose and describe the use of higher order Conditional Random Fields for query-template protein alignment. The experiments carried out on different public datasets validate our proposal, especially on distantly-related protein pairs which are the most difficult to align.

Performance of protein-structure predictions with the physics-based UNRES force field in CASP11.

Freely available on the web at http://www.unres.pl/ CONTACT: has5@cornell.edu.

Template-based protein-protein docking exploiting pairwise interfacial residue restraints.

Although many advanced and sophisticated ab initio approaches for modeling protein–protein complexes have been proposed in past decades, template-based modeling (TBM) remains the most accurate and widely used approach, given a reliable template is available. However, there are many different ways to exploit template information in the modeling process. Here, we systematically evaluate and benchmark a TBM method that uses conserved interfacial residue pairs as docking distance restraints [referred to as alpha carbon–alpha carbon (CA-CA)-guided docking]. We compare it with two other template-based protein–protein modeling approaches, including a conserved non-pairwise interfacial residue restrained docking approach [referred to as the ambiguous interaction restraint (AIR)-guided docking] and a simple superposition-based modeling approach. Our results show that, for most cases, the CA-CA-guided docking method outperforms both superposition with refinement and the AIR-guided docking method. We emphasize the superiority of the CA-CA-guided docking on cases with medium to large conformational changes, and interactions mediated through loops, tails or disordered regions. Our results also underscore the importance of a proper refinement of superimposition models to reduce steric clashes. In summary, we provide a benchmarked TBM protocol that uses conserved pairwise interface distance as restraints in generating realistic 3D protein–protein interaction models, when reliable templates are available. The described CA-CA-guided docking protocol is based on the HADDOCK platform, which allows users to incorporate additional prior knowledge of the target system to further improve the quality of the resulting models.

AIDA: ab initio domain assembly for automated multi-domain protein structure prediction and domain-domain interaction prediction.

Most proteins consist of multiple domains, independent structural and evolutionary units that are often reshuffled in genomic rearrangements to form new protein architectures. Template-based modeling methods can often detect homologous templates for individual domains, but templates that could be used to model the entire query protein are often not available. We have developed a fast docking algorithm ab initio domain assembly (AIDA) for assembling multi-domain protein structures, guided by the ab initio folding potential. This approach can be extended to discontinuous domains (i.e. domains with 'inserted' domains). When tested on experimentally solved structures of multi-domain proteins, the relative domain positions were accurately found among top 5000 models in 86% of cases. AIDA server can use domain assignments provided by the user or predict them from the provided sequence. The latter approach is particularly useful for automated protein structure prediction servers. The blind test consisting of 95 CASP10 targets shows that domain boundaries could be successfully determined for 97% of targets. The AIDA package as well as the benchmark sets used here are available for download at http://ffas.burnham.org/AIDA/. adam@sanfordburnham.org Supplementary data are available at Bioinformatics online.

Extended template-based modeling and evaluation method using consensus of binding mode of GPCRs for virtual screening.

G-protein-coupled receptors (GPCRs) are a pharmaceutically important protein family because they mediate numerous physiological functions. The crystal structures of several GPCR subtypes have been determined recently, encouraging efforts to apply structure-based virtual screening (SBVS) along with ligand-based virtual screening (LBVS) to improve the hit rate of active ligands from large chemical libraries. Three-dimensional models are also necessary for GPCR targets whose structures are unknown. Current challenges include the selection of structural templates from available structurally known GPCRs to use for accurate modeling and understanding the diversity of sites recognizing distinct ligands. We have developed and validated an extended template-based modeling and evaluation method for SBVS. Models were generated using a fragmental template procedure in addition to typical template-based modeling methods. The reliability of the models was evaluated using a virtual screening test with known active ligands and decoys and the consensus of the binding mode using the protein-ligand interaction fingerprint (PLIF) derived from the results of docking simulations. This novel workflow was applied to three targets with known structures (human dopamine receptor 3, human histamine H1 receptor, and human delta opioid receptor) and to a target with an unknown structure (human serotonin 2A receptor). In each case, model structures having high ligand selectivity with consensus binding mode were generated.

Improvement in Low-Homology Template-Based Modeling by Employing a Model Evaluation Method with Focus on Topology

Many template-based modeling (TBM) methods have been developed over the recent years that allow for protein structure prediction and for the study of structure-function relationships for proteins. One major problem all TBM algorithms face, however, is their unsatisfactory performance when proteins under consideration are low-homology. To improve the performance of TBM methods for such targets, a novel model evaluation method was developed here, and named MEFTop. Our novel method focuses on evaluating the topology by using two novel groups of features. These novel features included secondary structure element (SSE) contact information and 3-dimensional topology information. By combining MEFTop algorithm with FR-t5, a threading program developed by our group, we found that this modified TBM program, which was named FR-t5-M, exhibited significant improvements in predictive abilities for low-homology protein targets. We further showed that the MEFTop could be a generalized method to improve threading programs for low-homology protein targets. The softwares (FR-t5-M and MEFTop) are available to non-commercial users at our website: http://jianglab.ibp.ac.cn/lims/FRt5M/FRt5M.html.

Principles, Challenges and Advances in ab initio Protein Structure Prediction

The gap between known protein sequences and structures is increasing rapidly and experimental methods alone will not be able to fill in this gap. Therefore it is necessary to use computational methods to predict protein structures. Template based modeling methods could be used for sequences, which have detectable relationship with sequences of one or more experimentally determined protein structures. For predicting the structure of proteins, which does not share a detectable sequence relationship with experimental structures, ab initio protein structure prediction techniques must be used. The methods under ab initio protein structure prediction category aim to predict the structure of a protein from the sequence information alone, without any explicit use of previously known structures. These methods use thermodynamic principles and try to identify the native structure of a protein as the global minimum of a potential energy landscape. However, such methods are computationally complex and are extraordinarily challenging. There has been significant progress in the development of ab inito protein structure prediction methods over the past few years. This review describes the basic principles, the complexity, challenges and recent progresses of ab initio protein structure prediction.

GalaxyTBM: template-based modeling by building a reliable core and refining unreliable local regions

BackgroundProtein structures can be reliably predicted by template-based modeling (TBM) when experimental structures of homologous proteins are available. However, it is challenging to obtain structures more accurate than the single best templates by either combining information from multiple templates or by modeling regions that vary among templates or are not covered by any templates.ResultsWe introduce GalaxyTBM, a new TBM method in which the more reliable core region is modeled first from multiple templates and less reliable, variable local regions, such as loops or termini, are then detected and re-modeled by an ab initio method. This TBM method is based on “Seok-server,” which was tested in CASP9 and assessed to be amongst the top TBM servers. The accuracy of the initial core modeling is enhanced by focusing on more conserved regions in the multiple-template selection and multiple sequence alignment stages. Additional improvement is achieved by ab initio modeling of up to 3 unreliable local regions in the fixed framework of the core structure. Overall, GalaxyTBM reproduced the performance of Seok-server, with GalaxyTBM and Seok-server resulting in average GDT-TS of 68.1 and 68.4, respectively, when tested on 68 single-domain CASP9 TBM targets. For application to multi-domain proteins, GalaxyTBM must be combined with domain-splitting methods.ConclusionApplication of GalaxyTBM to CASP9 targets demonstrates that accurate protein structure prediction is possible by use of a multiple-template-based approach, and ab initio modeling of variable regions can further enhance the model quality.

Assessment of CASP7 predictions for template-based modeling targets

This manuscript presents the assessment of the template-based modeling category of the seventh Critical Assessment of Techniques for Protein Structure Prediction (CASP7). The accuracy of predicted protein models for 108 target domains was assessed based on a detailed comparison between the experimental and predicted structures. The assessment was performed using numerical measures for backbone and structural alignment accuracy, and by scoring correctly modeled hydrogen bond interactions in the predictions. Based on these criteria, our statistical analysis identified a number of groups whose predictions were on average significantly more accurate. Furthermore, the predictions for six target proteins were evaluated for the accuracy of their modeled cofactor binding sites. We also assessed the ability of predictors to improve over the best available single template structure, which showed that the best groups produced models closer to the target structure than the best single template for a significant number of targets. In addition, we assessed the accuracy of the error estimates (local confidence values) assigned to predictions on a per residue basis. Finally, we discuss some general conclusions about the state of the art of template-based modeling methods and their usefulness for practical applications.