Template-based Modeling Research Articles

Decrypting the effect of sulfate stress on the cyanobacterium Anabaena sp. PCC 7120 using physiological, proteome and transcript analysis in conjunction with bioinformatics

Sulfate limitation affects growth and various metabolic processes in cyanobacteria. However, there is little understanding of the underlying mechanism of these processes. Hence, the study investigated the effects of sulfate limitation in the cyanobacterium, Anabaena sp. PCC 7120. Results showed that sulfate stress affected photopigment contents in the test cyanobacterium delineating reduced growth. Proteome profile showed that differential expression profiles of different proteins involved in the various metabolic pathways suggested sulfur-dependent pathways were altered due to sulfate limitation. The most damaging effect was observed on the proteins belonging to the photosynthesis pathway, affirming the relationship between sulfur and carbon metabolism. Expression of proteins involved in the pathways such as nitrogen, protein synthesis, chaperone category, amino acid biosynthesis, and redox-responsive proteins were down-regulated. Moreover, forty-seven hypothetical proteins were also influenced by sulfate limitation. Molecular chaperones act as a physical barrier to prevent the aggregation of housekeeping proteins. Down expression of these proteins might affect the sustainability of the other proteins. Besides, transcript levels of the selected proteins were in concordance with physiological as well as proteome data. Additionally, template-based modeling and ab initio full atomic relaxation were performed via bioinformatic analysis through different tools to gain deeper insight into the structural and functional attributes of hypothetical proteins. Based on study, it could be concluded that mostly the proteins were downregulated which may lead to PCD in the cyanobacterium and also suggested that sulfur is a central hub and crucial for various sulfur-dependent pathways.

Protein Homology Modeling in the Low Sequence Similarity Regime

Predicting the 3-D structure of a protein from its sequence based on a template protein structure is still one of the most exact modeling techniques present today. However, template-based modeling is heavily dependent on the selection of a single template structure and the sequence alignment between target and template. Mainly when the target and template sequence identity is low, the error from the alignment introduces larger errors to the model structure. An iterative method to correct such alignment mistakes is used in this study with a benchmark set from CASP in the extremely low sequence-identity regime. This is a protocol developed and tested before and it evaluates the alignment quality by building rough 3-D models for each alignment. Then by using a genetic algorithm it iteratively creates a new set of alignments. Since the method evaluates models, not sequence alignments, structural features are automatically incorporated into the alignment protocol. In the current study, models from structural alignment have been built by Modeller program to show the maximum possible quality of the model that can be obtained from that template structure with the iterative modeling protocol. Then the results and correctly aligned segments from the iterative modeling protocol are analyzed. Finally, it has been shown that if a good local fragment assessment scoring function is developed, the correctly aligned segments exist in the pool of alignments created by the protocol. Thus, the improvement of modeling in the low sequence identity regime is conceivable.

IBio-GATS-A Semi-Automated Workflow for Structural Modelling of Insect Odorant Receptors.

Insects utilize seven transmembrane (7TM) odorant receptor (iOR) proteins, with an inverted topology compared to G-protein coupled receptors (GPCRs), to detect chemical cues in the environment. For pest biocontrol, chemical attractants are used to trap insect pests. However, with the influx of invasive insect pests, novel odorants are urgently needed, specifically designed to match 3D iOR structures. Experimental structural determination of these membrane receptors remains challenging and only four experimental iOR structures from two evolutionarily distant organisms have been solved. Template-based modelling (TBM) is a complementary approach, to generate model structures, selecting templates based on sequence identity. As the iOR family is highly divergent, a different template selection approach than sequence identity is needed. Bio-GATS template selection for GPCRs, based on hydrophobicity correspondence, has been morphed into iBio-GATS, for template selection from available experimental iOR structures. This easy-to-use semi-automated workflow has been extended to generate high-quality models from any iOR sequence from the selected template, using Python and shell scripting. This workflow was successfully validated on Apocrypta bakeri Orco and Machilis hrabei OR5 structures. iBio-GATS models generated for the fruit fly iOR, OR59b and Orco, yielded functional ligand binding results concordant with experimental mutagenesis findings, compared to AlphaFold2 models.

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.

RNA threading with secondary structure and sequence profile.

RNA threading aims to identify remote homologies for template-based modeling of RNA 3D structure. Existing RNA alignment methods primarily rely on secondary structure alignment. They are often time- and memory-consuming, limiting large-scale applications. In addition, the accuracy is far from satisfactory. Using RNA secondary structure and sequence profile, we developed a novel RNA threading algorithm, named RNAthreader. To enhance the alignment process and minimize memory usage, a novel approach has been introduced to simplify RNA secondary structures into compact diagrams. RNAthreader employs a two-step methodology. Initially, integer programming and dynamic programming are combined to create an initial alignment for the simplified diagram. Subsequently, the final alignment is obtained using dynamic programming, taking into account the initial alignment derived from the previous step. The benchmark test on 80 RNAs illustrates that RNAthreader generates more accurate alignments than other methods, especially for RNAs with pseudoknots. Another benchmark, involving 30 RNAs from the RNA-Puzzles experiments, exhibits that the models constructed using RNAthreader templates have a lower average RMSD than those created by alternative methods. Remarkably, RNAthreader takes less than two hours to complete alignments with ∼5000 RNAs, which is 3-40 times faster than other methods. These compelling results suggest that RNAthreader is a promising algorithm for RNA template detection. https://yanglab.qd.sdu.edu.cn/RNAthreader.

CoDock-Ligand: combined template-based docking and CNN-based scoring in ligand binding prediction

For ligand binding prediction, it is crucial for molecular docking programs to integrate template-based modeling with a precise scoring function. Here, we proposed the CoDock-Ligand docking method that combines template-based modeling and the GNINA scoring function, a Convolutional Neural Network-based scoring function, for the ligand binding prediction in CASP15. Among the 21 targets, we obtained successful predictions in top 5 submissions for 14 targets and partially successful predictions for 4 targets. In particular, for the most complicated target, H1114, which contains 56 metal cofactors and small molecules, our docking method successfully predicted the binding of most ligands. Analysis of the failed systems showed that the predicted receptor protein presented conformational changes in the backbone and side chains of the binding site residues, which may cause large structural deviations in the ligand binding prediction. In summary, our hybrid docking scheme was efficiently adapted to the ligand binding prediction challenges in CASP15.

Utilization of AlphaFold models for drug discovery: Feasibility and challenges. Histone deacetylase 11 as a case study

Histone deacetylase 11 (HDAC11), an enzyme that cleaves acyl groups from acylated lysine residues, is the sole member of class IV of HDAC family with no reported crystal structure so far. The catalytic domain of HDAC11 shares low sequence identity with other HDAC isoforms which complicates the conventional template-based homology modeling. AlphaFold is a neural network machine learning approach for predicting the 3D structures of proteins with atomic accuracy even in absence of similar structures. However, the structures predicted by AlphaFold are missing small molecules as ligands and cofactors. In our study, we first optimized the HDAC11 AlphaFold model by adding the catalytic zinc ion followed by assessment of the usability of the model by docking of the selective inhibitor FT895. Minimization of the optimized model in presence of transplanted inhibitors, which have been described as HDAC11 inhibitors, was performed. Four complexes were generated and proved to be stable using three replicas of 50 ns MD simulations and were successfully utilized for docking of the selective inhibitors FT895, MIR002 and SIS17. For SIS17, The most reasonable pose was selected based on structural comparison between HDAC6, HDAC8 and the HDAC11 optimized AlphaFold model. The manually optimized HDAC11 model is thus able to explain the binding behavior of known HDAC11 inhibitors and can be used for further structure-based optimization.

Improved computational epitope profiling using structural models identifies a broader diversity of antibodies that bind to the same epitope.

The function of an antibody is intrinsically linked to the epitope it engages. Clonal clustering methods, based on sequence identity, are commonly used to group antibodies that will bind to the same epitope. However, such methods neglect the fact that antibodies with highly diverse sequences can exhibit similar binding site geometries and engage common epitopes. In a previous study, we described SPACE1, a method that structurally clustered antibodies in order to predict their epitopes. This methodology was limited by the inaccuracies and incomplete coverage of template-based modeling. In addition, it was only benchmarked at the level of domain-consistency on one virus class. Here, we present SPACE2, which uses the latest machine learning-based structure prediction technology combined with a novel clustering protocol, and benchmark it on binding data that have epitope-level resolution. On six diverse sets of antigen-specific antibodies, we demonstrate that SPACE2 accurately clusters antibodies that engage common epitopes and achieves far higher dataset coverage than clonal clustering and SPACE1. Furthermore, we show that the functionally consistent structural clusters identified by SPACE2 are even more diverse in sequence, genetic lineage, and species origin than those found by SPACE1. These results reiterate that structural data improve our ability to identify antibodies that bind to the same epitope, adding information to sequence-based methods, especially in datasets of antibodies from diverse sources. SPACE2 is openly available on GitHub (https://github.com/oxpig/SPACE2).

A novel homozygous PIGO mutation associated with severe infantile epileptic encephalopathy, profound developmental delay and psychomotor retardation: structural and 3D modelling investigations and genotype-phenotype correlation.

The PIGO gene encodes the GPI-ethanolamine phosphate transferase 3, which is crucial for the final synthetic step of the glycosylphosphatidylinositol-anchor serving to attach various proteins to their cell surface. These proteins are intrinsic for normal neuronal and embryonic development. In the current research work, a clinical investigation was conducted on a patient from a consanguineous family suffering from epileptic encephalopathy, characterized by severe seizures, developmental delay, hypotonia, ataxia and hyperphosphatasia. Molecular analysis was performed using Whole Exome Sequencing (WES). The molecular investigation revealed a novel homozygous variant c.1132C > T in thePIGOgene, in which a highly conserved Leucine was changed to a Phenylalanine (p.L378F). To investigate the impact of the non-synonymous mutation, a 3D structural model of the PIGO protein was generated using the AlphaFold protein structure database as a resource for template-based tertiary structure modeling. A structural analysis by applying some bioinformatic tools on both variants 378L and 378F models predicted the pathogenicity of the non-synonymous mutation and its potential functional and structural effects on PIGO protein. We also discussed the phenotypic and genotypic variability associated with the PIGO deficiency. To our best knowledge, this is the first report of a patient diagnosed with infantile epileptic encephalopathy showing a high elevation of serum alkaline phosphatase level. Our findings, therefore, widen the genotype and phenotype spectrum of GPI-anchor deficiencies and broaden the cohort of patients withPIGOassociated epileptic encephalopathy with an elevated serum alkaline phosphatase level.

Large-scale template-based structural modeling of T-cell receptors with known antigen specificity reveals complementarity features.

T-cell receptor (TCR) recognition of foreign peptides presented by the major histocompatibility complex (MHC) initiates the adaptive immune response against pathogens. While a large number of TCR sequences specific to different antigenic peptides are known to date, the structural data describing the conformation and contacting residues for TCR-peptide-MHC complexes is relatively limited. In the present study we aim to extend and analyze the set of available structures by performing highly accurate template-based modeling of these complexes using TCR sequences with known specificity. Identification of CDR3 sequences and their further clustering, based on available spatial structures, V- and J-genes of corresponding T-cell receptors, and epitopes, was performed using the VDJdb database. Modeling of the selected CDR3 loops was conducted using a stepwise introduction of single amino acid substitutions to the template PDB structures, followed by optimization of the TCR-peptide-MHC contacting interface using the Rosetta package applications. Statistical analysis and recursive feature elimination procedures were carried out on computed energy values and properties of contacting amino acid residues between CDR3 loops and peptides, using R. Using the set of 29 complex templates (including a template with SARS-CoV-2 antigen) and 732 specificity records, we built a database of 1585 model structures carrying substitutions in either TCRα or TCRβ chains with some models representing the result of different mutation pathways for the same final structure. This database allowed us to analyze features of amino acid contacts in TCR - peptide interfaces that govern antigen recognition preferences and interpret these interactions in terms of physicochemical properties of interacting residues. Our results provide a methodology for creating high-quality TCR-peptide-MHC models for antigens of interest that can be utilized to predict TCR specificity.

A Structural Refinement Technique for Protein-RNA Complexes Based on a Combination of AI-based Modeling and Flexible Docking: A Study of Musashi-1 Protein

Abstract An efficient structural refinement technique for protein-RNA complexes is proposed based on a combination of AI-based modeling and flexible docking. Specifically, an enhanced sampling method called parallel cascade selection molecular dynamics (PaCS-MD) was extended to include flexible docking to construct protein-RNA complexes from those obtained by AI-based modeling (AlphaFold2). With the present technique, the conformational sampling of flexible RNA regions is accelerated by PaCS-MD, enabling one to construct plausible models for protein-RNA complexes. For demonstration, PaCS-MD constructed several protein-RNA complexes of the RNA-binding Musashi-1 (MSI1) family of proteins, which were validated by comparing a group of crucial residues for RNA-binding with experimental complexes. Our analyses suggest that PaCS-MD improves the quality of complex modeling compared to the standard protocol based on template-based modeling (Phyre2). Furthermore, PaCS-MD could also be a beneficial technique for constructing complexes of non-native RNA-binding to proteins.

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.

Model templates: transdisciplinary application and entanglement

The omnipresence of the same basic equations, function forms, algorithms, and quantitative methods is one of the most spectacular characteristics of contemporary modeling practice. Recently, the emergence of the discussion of templates and template transfer has addressed this striking cross-disciplinary reach of certain mathematical forms and computational algorithms. In this paper, we develop a notion of a model template, consisting of its mathematical structure, ontology, prototypical properties and behaviors, focal conceptualizations, and the paradigmatic questions it addresses. We apply this notion to three widely disseminated and powerful model templates: the Sherrington-Kirkpatrick model of spin glasses, scale-free networks, and the Kuramoto model of synchronization. We argue that what appears to be an interdisciplinary model transfer between different domains turns out, from a broader perspective, to be the application of transdisciplinary model templates across a multitude of domains. We also point out a further feature of template-based modeling that so far has not been discussed: template entanglement. Such entanglement enhances and makes manifest the conceptual side of model templates.

To split or not to split: CASP15 targets and their processing into tertiary structure evaluation units.

Processing of CASP15 targets into evaluation units (EUs) and assigning them to evolutionary-based prediction classes is presented in this study. The targets were first split into structural domains based on compactness and similarity to other proteins. Models were then evaluated against these domains and their combinations. The domains were joined into larger EUs if predictors' performance on the combined units was similar to that on individual domains. Alternatively, if most predictors performed better on the individual domains, then they were retained as EUs. As a result, 112 evaluation units were created from 77 tertiary structure prediction targets. The EUs were assigned to four prediction classes roughly corresponding to target difficulty categories in previous CASPs: TBM (template-based modeling, easy or hard), FM (free modeling), and the TBM/FM overlap category. More than a third of CASP15 EUs were attributed to the historically most challenging FM class, where homology or structural analogy to proteins of known fold cannot be detected.

The ClusPro AbEMap web server for the prediction of antibody epitopes.

Antibodies play an important role in the immune system by binding to molecules called antigens at their respective epitopes. These interfaces or epitopes are structural entities determined by the interactions between an antibody and an antigen, making them ideal systems to analyze by using docking programs. Since the advent of high-throughput antibody sequencing, the ability to perform epitope mapping using only the sequence of the antibody has become a high priority. ClusPro, a leading protein-protein docking server, together with its template-based modeling version, ClusPro-TBM, have been re-purposed to map epitopes for specific antibody-antigen interactions by using the Antibody Epitope Mapping server (AbEMap). ClusPro-AbEMap offers three different modes for users depending on the information available on the antibody as follows: (i) X-ray structure, (ii) computational/predicted model of the structure or (iii) only the amino acid sequence. The AbEMap server presents a likelihood score for each antigen residue of being part of the epitope. We provide detailed information on the server's capabilities for the three options and discuss how to obtain the best results. In light of the recent introduction of AlphaFold2 (AF2), we also show how one of the modes allows users to use their AF2-generated antibody models as input. The protocol describes the relative advantages of the server compared to other epitope-mapping tools, its limitations and potential areas of improvement. The server may take 45-90 min depending on the size of the proteins.

Template-Based Modelling of the Structure of Fungal Effector Proteins

The discovery of new fungal effector proteins is necessary to enable the screening of cultivars for disease resistance. Sequence-based bioinformatics methods have been used for this purpose, but only a limited number of functional effector proteins have been successfully predicted and subsequently validated experimentally. A significant obstacle is that many fungal effector proteins discovered so far lack sequence similarity or conserved sequence motifs. The availability of experimentally determined three-dimensional (3D) structures of a number of effector proteins has recently highlighted structural similarities amongst groups of sequence-dissimilar fungal effectors, enabling the search for similar structural folds amongst effector sequence candidates. We have applied template-based modelling to predict the 3D structures of candidate effector sequences obtained from bioinformatics predictions and the PHI-BASE database. Structural matches were found not only with ToxA- and MAX-like effector candidates but also with non-fungal effector-like proteins—including plant defensins and animal venoms—suggesting the broad conservation of ancestral structural folds amongst cytotoxic peptides from a diverse range of distant species. Accurate modelling of fungal effectors were achieved using RaptorX. The utility of predicted structures of effector proteins lies in the prediction of their interactions with plant receptors through molecular docking, which will improve the understanding of effector–plant interactions.

PrePPI: A Structure Informed Proteome-wide Database of Protein–Protein Interactions

We present an updated version of the Predicting Protein-Protein Interactions (PrePPI) webserver which predicts PPIs on a proteome-wide scale. PrePPI combines structural and non-structural evidence within a Bayesian framework to compute a likelihood ratio (LR) for essentially every possible pair of proteins in a proteome; the current database is for the human interactome. The structural modeling (SM) component is derived from template-based modeling and its application on a proteome-wide scale is enabled by a unique scoring function used to evaluate a putative complex. The updated version of PrePPI leverages AlphaFold structures that are parsed into individual domains. As has been demonstrated in earlier applications, PrePPI performs extremely well as measured by receiver operating characteristic curves derived from testing on E. coli and human protein–protein interaction (PPI) databases. A PrePPI database of ∼1.3 million human PPIs can be queried with a webserver application that comprises multiple functionalities for examining query proteins, template complexes, 3D models for predicted complexes, and related features (https://honiglab.c2b2.columbia.edu/PrePPI). PrePPI is a state-of-the-art resource that offers an unprecedented structure-informed view of the human interactome.

Computer-Aided Formalization of Internal Consistent Product Family Models

The product line approach enables the rapid building of software products by reusing core assets. How to manage the commonality and variability of product family is the core problem of product line engineering. The feature modeling approach is one of the main approaches to solve this problem. However, this approach does not define the internal functions of features. Some researchers use formal descriptions to define the function of features, but they do not consider the connection between feature specification and feature model. These product family modeling approaches result in product family models that suffer from inconsistent descriptions, which are difficult to resolve manually. To this end, this paper proposes a computer-aided approach to construct an internally consistent product family model. The method performs an inconsistency check during the formalization of each feature in the given initial feature model, pushing the evolution of the feature model and feature specification. The creation and evolution of feature models are achieved through a template-based feature modeling approach. The method is implemented by means of templates, each of which contains two parts of knowledge: the elements required for a class of feature model, and a way to convert the given elements into a formal description. An inconsistency check is achieved through a rule-based detection approach. Each rule contains two parts: a class of conditions that the feature model and the feature specification satisfy, and the actions to resolve inconsistencies once the conditions are satisfied. We have developed a tool. The tool guides the user to build an internally consistent product family model that contains the feature model and feature specifications. Experiments have shown that our approach can provide effective automation support for the process of formal product family modeling.

Benchmarking TriadAb using targets from the second antibody modeling assessment.

Computational modeling and design of antibodies has become an integral part of today's research and development in antibody therapeutics. Here we describe the Triad Antibody Homology Modeling (TriadAb) package, a functionality of the Triad protein design platform that predicts the structure of any heavy and light chain sequences of an antibody Fv domain using template-based modeling. To gauge the performance of TriadAb, we benchmarked against the results of the Second Antibody Modeling Assessment (AMA-II). On average, TriadAb produced main-chain carbonyl root-mean-square deviations between models and experimentally determined structures at 1.10Å, 1.45Å, 1.41Å, 3.04Å, 1.47Å, 1.27Å, 1.63Å in the framework and the six complementarity-determining regions (H1, H2, H3, L1, L2, L3), respectively. The inaugural results are comparable to those reported in AMA-II, corroborating with our internal bench-based experiences that models generated using TriadAb are sufficiently accurate and useful for antibody engineering using the sequence design capabilities provided by Triad.

Improving Protein-Ligand Interaction Modeling with cryo-EM Data, Templates, and Deep Learning in 2021 Ligand Model Challenge.

Elucidating protein-ligand interaction is crucial for studying the function of proteins and compounds in an organism and critical for drug discovery and design. The problem of protein-ligand interaction is traditionally tackled by molecular docking and simulation, which is based on physical forces and statistical potentials and cannot effectively leverage cryo-EM data and existing protein structural information in the protein-ligand modeling process. In this work, we developed a deep learning bioinformatics pipeline (DeepProLigand) to predict protein-ligand interactions from cryo-EM density maps of proteins and ligands. DeepProLigand first uses a deep learning method to predict the structure of proteins from cryo-EM maps, which is averaged with a reference (template) structure of the proteins to produce a combined structure to add ligands. The ligands are then identified and added into the structure to generate a protein-ligand complex structure, which is further refined. The method based on the deep learning prediction and template-based modeling was blindly tested in the 2021 EMDataResource Ligand Challenge and was ranked first in fitting ligands to cryo-EM density maps. These results demonstrate that the deep learning bioinformatics approach is a promising direction for modeling protein-ligand interactions on cryo-EM data using prior structural information.