Skip to main content
Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
Nucleic Acids Res. 2018 Jul 2; 46(Web Server issue): W296–W303.
Published online 2018 May 21. doi: 10.1093/nar/gky427
PMCID: PMC6030848
PMID: 29788355

SWISS-MODEL: homology modelling of protein structures and complexes

Andrew Waterhouse,1,2 Martino Bertoni,1,2 Stefan Bienert,1,2 Gabriel Studer,1,2 Gerardo Tauriello,1,2 Rafal Gumienny,1,2 Florian T Heer,1,2 Tjaart A P de Beer,1,2 Christine Rempfer,1,2 Lorenza Bordoli,1,2 Rosalba Lepore,1,2 and Torsten Schwede1,2
Author information Article notes Copyright and License information PMC Disclaimer

Associated Data

Supplementary Materials

Abstract

Homology modelling has matured into an important technique in structural biology, significantly contributing to narrowing the gap between known protein sequences and experimentally determined structures. Fully automated workflows and servers simplify and streamline the homology modelling process, also allowing users without a specific computational expertise to generate reliable protein models and have easy access to modelling results, their visualization and interpretation. Here, we present an update to the SWISS-MODEL server, which pioneered the field of automated modelling 25 years ago and been continuously further developed. Recently, its functionality has been extended to the modelling of homo- and heteromeric complexes. Starting from the amino acid sequences of the interacting proteins, both the stoichiometry and the overall structure of the complex are inferred by homology modelling. Other major improvements include the implementation of a new modelling engine, ProMod3 and the introduction a new local model quality estimation method, QMEANDisCo. SWISS-MODEL is freely available at https://swissmodel.expasy.org.

INTRODUCTION

Three-dimensional structures of proteins provide valuable insights into their function on a molecular level and inform a broad spectrum of applications in life science research. Often, complexes of proteins are central to many cellular processes. A detailed description of their interactions and the overall quaternary structure is essential for a comprehensive understanding of biological systems, how protein complexes and networks operate and how we can modulate them (1,2). Given their biological relevance, it is not surprising that the number of large complexes deposited per year in the Protein Data Bank (PDB) is growing rapidly (3). A significant contribution to this trend originates from the continuous progress of structure determination technologies, including recent developments of Electron Microscopy (EM) based methods, which are particularly suited for large macromolecular assemblies (4). Still, compared to high-throughput methods for screening protein-protein interactions (i.e. yeast two-hybrid, affinity purification, phage-display etc.), the rate at which novel complex structures are determined experimentally is considerably lower. This uneven growth calls for computational methods to fill the gap.

Several approaches have been developed to address the computational prediction of protein-protein interactions (5). Co-evolution methods, based on correlated amino acid mutations in deep multiple sequence alignments (MSA), are efficiently used to identify interacting proteins based on sequence information alone (6,7). When the 3D structures of the binding partners are available, or can be reliably modelled, docking methods can be used to obtain a three-dimensional model of the complex based on geometric and physicochemical complementarity of the interacting molecules (8–11). Efficiently handling protein flexibility is still one of the major challenges in the development of effective docking simulation software; hence these methods are generally more accurate when little or no conformational change is required for binding. According to the community-wide experiment CAPRI (Critical Assessment of PRedicted Interactions (12)), considerable progress has been made in the field with the development of hybrid modelling strategies, that are able to incorporate available experimental information on the interaction (i.e. crosslinks, NMR, SAXS etc.) as constraints in the simulation of the docking process (13–15). Results from latest assessments show that significantly improved quality of models is obtained when multi-chain template information is available and used for modelling (16).

With more experimentally determined structures of protein complexes becoming available, it has been observed that interacting interfaces are often conserved among homologous complexes (17) and that templates are available for most of the known protein-protein interactions (18). These observations provided the rationale for comparative, or homology modelling, of protein complexes. Similar to comparative modelling of monomeric proteins, the information of a protein's quaternary structure is transferred by homology to another one, and a model of the complex is obtained based on the structures of the interacting homologs, i.e. interologs, as templates (19–21). The approach can be scaled to entire genomes and applied to binary as well as to higher-order protein assemblies (17,18,22,23). As highlighted by the introduction of the first assessment of protein assemblies in the recent CASP XII experiment (24), comparative modelling of protein complexes is receiving much attention and is expected to play a relevant role in the elucidation of the protein quaternary structure space.

SWISS-MODEL https://swissmodel.expasy.org was the first fully automated protein homology modelling server and has been continuously improved during the last 25 years (25–30). Its modelling functionality has been recently extended to include the modelling of homo- and heteromeric complexes, given the amino acid sequences of the interacting partners as starting point. Other recently introduced features include the development of a new modelling engine, ProMod3, with increased accuracy of the produced models, and an improved local model quality estimation method (QMEANDisCo) based on a novel version of QMEAN (31).

SWISS-MODEL currently generates ∼3000 models a day (∼2 models per minute), up from ∼1500 models a day in 2014 (30), making it one of the most widely used structure modelling servers worldwide. Its performance is continuously evaluated and compared with other state-of-the art servers in the field. To this aim, we are actively participating to the CAMEO project (Continuous Automated Model Evaluation, https://cameo3d.org) (32), a fully automated blind prediction assessment based on weekly pre-release of sequences from the PDB (33), allowing us to constantly monitor and improve the performance of the server.

MATERIALS AND METHODS

The modelling workflow

In comparative modelling, a 3D protein model of a target sequence is generated by extrapolating experimental information from an evolutionary related protein structure that serves as a template. In SWISS-MODEL, the default modelling workflow consists of the following main steps:

  1. Input data: The target protein can be provided as amino acid sequence, either in FASTA, Clustal format or as a plain text. Alternatively, a UniProtKB accession code (34) can be specified. If the target protein is heteromeric, i.e. it consists of different protein chains as subunits, amino acid sequences or UniProtKB accession codes must be specified for each subunit.
  2. Template search: Data provided in step 1 serve as a query to search for evolutionary related protein structures against the SWISS-MODEL template library SMTL (30). SWISS-MODEL performs this task by using two database search methods: BLAST (35,36), which is fast and sufficiently accurate for closely related templates, and HHblits (37), which adds sensitivity in case of remote homology.
  3. Template selection: When the template search is complete, templates are ranked according to expected quality of the resulting models, as estimated by Global Model Quality Estimate (GMQE) (30) and Quaternary Structure Quality Estimate (QSQE) (23). Top-ranked templates and alignments are compared to verify whether they represent alternative conformational states or cover different regions of the target protein. In this case, multiple templates are selected automatically and different models are built accordingly. To provide the user with the option to use alternative templates than those selected automatically, all templates are shown in a tabular form with a descriptive set of features. In addition, interactive graphical views facilitate the analysis and comparison of available templates in terms of their three-dimensional structures, sequence similarity and quaternary structure features.
  4. Model building: For each selected template, a 3D protein model is automatically generated by first transferring conserved atom coordinates as defined by the target-template alignment. Residue coordinates corresponding to insertions/deletions in the alignment are generated by loop modelling and a full-atom protein model is obtained by constructing the non-conserved amino acid side chains. SWISS-MODEL relies on the OpenStructure computational structural biology framework (38) and the ProMod3 modelling engine to perform this step. For more detailed information on model building we refer to a dedicated section in Results.
  5. Model quality estimation: To quantify modelling errors and give estimates on expected model accuracy, SWISS-MODEL relies on the QMEAN scoring function (31). QMEAN uses statistical potentials of mean force to generate global and per residue quality estimates. The local quality estimates are enhanced by pairwise distance constraints that represent ensemble information from all template structures found. For more information on quality estimation we refer to a dedicated section in Results.

SWISS-MODEL allows for further customization of steps 1 and 3. Expert users can directly upload custom target-template sequence alignments, template structures or DeepView project files (26) in separate input forms.

The SWISS-MODEL template library

The SWISS-MODEL Template Library (SMTL), available at https://swissmodel.expasy.org/templates/, is a curated template library, which is updated on a weekly basis according to the new PDB release (33). Every deposited PDB structure is automatically processed, annotated and indexed to support efficient querying of high quality structural data. SMTL entries are organized by quaternary structure assemblies, according to the ‘Biological Assembly’ annotation specified in the PDB. Biologically relevant ligands are annotated accordingly, as described in (30), and the annotation is then used by the modelling engine to determine whether a ligand is considered for inclusion into the final model. As of January 2018, the SMTL contains coordinates for a total of 92 474 unique protein sequences, mapping to 219 350 biological units, annotated as follows: 113 639 monomers, 71 555 homo-oligomers and 34 156 hetero-oligomers.

Integration with the SWISS-MODEL repository and cross-links to other services

The SWISS-MODEL Repository (39) (SMR, https://swissmodel.expasy.org/repository) is a database of automatically generated homology models for relevant model organisms and experimental structure information for all sequences in UniProtKB (34). Whenever a UniProtKB sequence is submitted to SWISS-MODEL, the generated model is automatically deposited into the SMR along with all data used to generate the model. Currently, the SMR contains 1 067 355 models from SWISS-MODEL and 129 416 structures from PDB with mapping to UniProtKB.

To facilitate exploration of available information on a given target protein, SWISS-MODEL provides cross-links to various other resources and databases. We include links to the RCSB (33), PDBsum (40), PDBe (41), CATH (42) and SwissDock (43). In addition, we also provide direct access to a specialised server for antibody modelling. The pre-screening of the target sequence has been extended in order to automatically identify whether an immunoglobulin sequence is present in the input. If a matching sequence signal is detected, data can be sent to the Prediction of Immunoglobulin Structure server PIGSPro (44–46) where the model of the antibody is generated according to the canonical structure method (47–49).

Documentation and technical implementation

An updated version of the documentation is provided to reflect the latest changes of the current SWISS-MODEL release. Tutorial pages and examples have been updated according to latest options and features available. Additionally, a video tutorial with a step-by-step guide on how to generate a model using SWISS-MODEL is available at https://swissmodel.expasy.org/docs/tutorial.

SWISS-MODEL is implemented in Python and Django with Javascript and jQuery for the front-end, and Python, C++ and OpenStructure (38) for the back-end functions. For visualization of protein structures, users can select between two interactive JavaScript/WebGL molecule viewers, PV (https://biasmv.github.io/pv/) and NGL (https://github.com/arose/ngl) (50). ProMod3 was developed using OpenStructure; its core is written in efficient C/C++ and a Python interface is provided for rapid prototyping.

RESULTS AND DISCUSSIONS

The ProMod3 modelling engine

The modelling engine is the heart of SWISS-MODEL. It builds an atomistic protein model given a template structure and a target-template sequence alignment. Until recently, the software package ProMod-II (26), using MODELLER (51) as a fall-back, was in use to perform this task. As of June 2016, the newly developed modelling engine ProMod3 is used exclusively. ProMod3 has been designed with the aim of providing rapid and flexible prototyping for future modelling developments in SWISS-MODEL.

Like its predecessors, ProMod3 extracts structural information from an aligned template structure in Cartesian space. Insertions and deletions, as defined by the sequence alignment, are resolved by first searching for viable candidate fragments in a structural database. This is a relevant modification, as ProMod-II mainly relied on ab-initio techniques to perform this step. Final candidates are selected using statistical potentials of mean force scoring methods. If no suitable fragments can be found, a conformational space search is performed using Monte Carlo sampling. Non-conserved side-chains are modelled using the 2010 backbone dependent rotamer library from the Dunbrack group (52). The optimal configuration of rotamers is estimated using the graph based TreePack algorithm (53) by minimizing the SCWRL4 energy function (54). As a final step, small structural distortions, unfavourable interactions or clashes introduced during the modelling process, are resolved by energy minimization. ProMod3 uses the OpenMM library (55) to perform the computations and the CHARMM27 force field (56) for parameterization.

A direct comparison between the previous and updated modelling engines has been performed in the context of the CAMEO experiment using 250 target proteins collected during the time range 20 October 2017–13 January 2018. For each target, a template search has been performed using HHblits against the SMTL at the time of the CAMEO submission. The best template, according to the HHblits e-value, and the corresponding target-template sequence alignment served as input for both engines. As shown in Figure ​Figure1,1, models generated with ProMod3 show significantly improved accuracy according to all-atom lDDT (Local Distance Difference Test) score (57), a superposition-free measure of the deviation of interatomic distances between model and native structures. The same also holds for other commonly used model quality metrics, i.e. GDT-HA (Global Distance Test High Accuracy score) (58) and TM-score (Template Modelling score) (59) (Supplementary Figure S1).

An external file that holds a picture, illustration, etc. Object name is gky427fig1.jpg

Performance comparison between ProMod-II and ProMod3 modelling engines. Performance is measured on a benchmark dataset of 250 targets collected during the CAMEO time range 20 October 2017–13 January 2018. For each target, the same template and target-template alignment were used as input for both modelling engines. Each data point represents the difference in model accuracy in terms of all-atom IDDT score. ProMod3 shows a statically significant improvement of 2.65 IDDT points on average (P-value = 1.1E–43) based on paired t-test.

Modelling the protein quaternary structure of homo- and hetero-oligomers

In SWISS-MODEL, we have recently introduced a new approach to model the stoichiometry and overall structure of protein complexes using the sequence of the interacting component as starting points (23). The method is based on a novel description of interface conservation as a function of evolutionary distance. The basic assumption is that biologically relevant interfaces are less free to vary than the rest of the protein surface (60,61). We capture such evolutionary constraints by measuring the ratio between interface and surface residue entropy distribution from multiple sequence alignments (MSA) of homologous proteins as a function of evolutionary distance. We employ this analysis, termed PPI fingerprint, to discriminate biologically relevant interfaces from crystal contacts, and for estimating the accuracy of models. Interface conservation analysis and geometric properties of protein complexes were used to train a supervised machine-learning algorithm, Support Vector Machines (SVM), to identify templates that maximize the estimated interface quality of the resulting model. The predicted interface quality, i.e. quaternary structure quality estimate (QSQE), is a score between 0 and 1 reflecting the expected accuracy of inter-chain contacts in a model based on a given alignment and template. Further details are provided in (23).

Model quality estimation

SWISS-MODEL provides quality estimates at several stages of the modelling process. Given a template structure and target-template alignment, the GMQE (30) and the QSQE (23) provide estimates of the expected quality of the resulting model at the tertiary and quaternary structure level. These estimates help the user identify optimal templates and are also utilized for the fully automated template selection procedure. Once models have been built, their quality is assessed by the QMEAN scoring function (31). QMEAN employs statistical potentials of mean force to generate quality estimates on a global and local scale. The latest version of QMEAN, QMEANDisCo, further enhances the accuracy of local quality estimates. It assesses the consistency of observed interatomic distances in the model with ensemble information extracted from experimentally determined protein structures that are homologues to the target sequence. To incorporate structural features, GMQE is updated after model building with the QMEAN global score and is then used for the model ranking. To facilitate interpretation of the obtained model quality estimates, the QMEAN global score is transformed to a Z-score, indicating whether the model scores as it would be expected from experimentally determined structures of similar size (31).

Performance comparison with other modelling servers

In order to provide objective assessments of modelling performance, SWISS-MODEL participates in the CAMEO project (https://cameo3d.org) (32). Taking some inspiration from CASP, CAMEO aims to provide a continuous, fully automated, assessment of predictions produced by various modelling servers using a common benchmark dataset of targets. CAMEO target sequences are obtained from the weekly pre-release of new PDB structures and submitted to participating methods at the same point in time. This ensures all servers have access to the same background information, i.e. same structures from PDB or protein sequences in UniProtKB, when running their predictions. Finally, in order to exclude trivial modelling cases, protein sequences exhibiting >85% sequence identity to available PDB structures are not considered in the CAMEO evaluation.

Based on the CAMEO results in the ‘3D Structure Prediction’ category, SWISS-MODEL is consistently ranked among the top-modelling servers for several crucial modelling aspects. Table ​Table11 shows the performance based on a benchmark dataset of 250 targets collected during the CAMEO time range 20 October 2017–13 January 2018. Full data on performance are provided as supplementary materials (Supplementary Tables S1–S7). Notably, SWISS-MODEL has the lowest response time to generate models and excels at model quality for binding sites (IDDT-BS), for high-quality models (lDDT-easy) and for quaternary structure prediction (QS-score). SWISS-MODEL is optimized for comparative protein modelling cases, where high-quality models can be generated and used in a variety of practical research applications (62). For difficult remote homology or de novo modelling targets, other methods perform better in the CAMEO assessment (63–65). It is worth mentioning that among the participating servers, only SWISS-MODEL and Robetta provided results for oligomeric targets. Therefore quaternary structure predictions were assessed on a common subset of oligomeric proteins where both methods returned a model, for a total of 32 targets. Finally, based on the assessment of model confidence, SWISS-MODEL significantly outperforms other modelling servers in providing accurate local confidence estimates of the returned models.

Table 1.

Performance comparison in the context of the CAMEO continuous evaluation platform
ServerResponse time (hh:mm:ss) (N = 168)lDDT total (N = 168)lDDT easy (N = 37)lDDT medium (N = 90)lDDT hard (N = 41)lDDT BS (N = 69)QS-Score (N = 32)Model confidence (N = 168)
SWISS-MODEL00:15:4866.2286.0169.7140.6770.8863.950.85
HHpredB01:16:15*65.9582.10*69.6843.1871.47-0.79*
NaiveBLAST01:20:27*58.93*82.86*64.20*25.76*63.88*-0.68*
PRIMO02:12:08*60.26*84.51*65.07*27.82*67.30*-0.67*
SPARKS-X02:35:21*63.14*80.06*65.57*42.5367.76*-0.54*
RaptorX06:28:57*69.15*83.35*72.10*49.88*68.85-0.65*
IntFOLD4-TS32:47:59*68.41*83.76*70.8849.11*71.65-0.84
Robetta37:00:07*71.60*85.1774.00*54.08*67.48*60.200.81*

Performance is measured based on a benchmark dataset of 250 targets collected during the CAMEO time range 20 October 2017–13 January 2018. Results from SWISS-MODEL and seven other modelling servers were collected from CAMEO and the performance evaluated on a common subset of targets where all compared servers returned a model. Each column indicates average performance values in terms of Response Time, model accuracy (IDDT, QS-score) and self-assessment of model quality (Model Confidence). lDDT evaluation has further been split according to CAMEOs definition of target difficulty; per column subset sizes are shown in brackets. Asterisks indicate a statistically significant difference (P-value < 0.05) compared to SWISS-MODEL based on paired t-test.

Case study: Modelling the Ferredoxin/Ferredoxin-NADP(+) Reductase complex

To illustrate the new features of SWISS-MODEL, we describe here the modelling of the hetero-dimeric complex formed by Ferredoxin-NADP(+) Reductase (FNR) and its physiological electron donor Ferredoxin (Fd). In higher plants, these proteins are part of the electron transport chain of thylakoid membranes where they catalyse the last step of NADP+ reduction. In non-photosynthetic tissues, i.e. roots, the reaction operates in the opposite direction and is mediated by the tissue specific isoforms of the enzymes (66). Crystal structures of the leaf electron transfer complex FNR:Fd have been reported from Zea mays, providing structural details of the protein–protein interactions involved in electron transfer during light dependent reactions of photosynthesis (67). Only recently, a three-dimensional structure of the FNR:Fd complex formed by the root isotypes has been determined (68). To illustrate the modelling of the root FNR:Fd complex, the native structure has been removed from the SMTL and is used only for validation and visualization purposes.

The amino acid sequences of root FNR (UniProtKB: B4G043) and root Fd (UniProtKB: P27788) from Z. mays were submitted to SWISS-MODEL. Results of the template search are shown in Figure ​Figure2A,2A, where templates are clustered and displayed in a decision tree according to their quaternary structure features: oligomeric state, stoichiometry, topology and interface similarity. Each leaf of the tree corresponds to a template and target-template alignment (based on HHblits, BLAST or both); templates are labelled with their SMTL ID; bars indicate sequence identity and coverage to the target (darker shades of blue indicate higher sequence identity). Three homologous template complexes could be identified, with similar coverage and identity to the target sequences: FNR sequence coverage between 75 and 78%, sequence identity between 50 and 52%; Fd sequence coverage 63–64% and sequence identity between 66 and 70%. Based on clustering results, available templates have the same oligomeric state, stoichiometry and topology. In terms of structural interface similarity, on the other hand, they form three different clusters. As shown in the PPI fingerprint plot (Figure ​(Figure2B),2B), two templates (SMTL ID: 1gaq.1 (67) and SMTL ID: 1ewy.1 (69)) display a similar interface conservation pattern, typical to that observed for biologically relevant interfaces (23). Instead, 3w5u.1 shows a different PPI fingerprint curve, with conservation score constantly close to zero, as typically observed for crystallization artefacts (23,70). Notably, the quality of the model interface based on this template is expected to be very low (QSQE = 0.13). Indeed, after inspecting the structure and the corresponding study (71), we could confirm that the interface observed in the 3w5u.1 biounit is the result of a cross-link experiment and, as it can be appreciated in Figure ​Figure2C,2C, does not correspond to the biologically functional interface. A stronger conservation signal is visible for template SMTL ID: 1gaq.1 (green line in Figure ​Figure2B),2B), which according to both tertiary and quaternary structure quality estimates (GMQE = 0.62; QSQE = 0.54) is considered the best template among the available options; hence it is selected for modelling. The resulting model is shown in Figure ​Figure2D,2D, where it has been superimposed onto the experimental structure of the complex (PDB ID: 5H5J (68), shown in light gray). Results of local quality assessment can be visualized onto the model where the colour gradient, from blue to red, indicates high to low quality as measured according to all-atom IDDT score. As it can be observed, three-dimensional structures of both FNR and Fd are modelled with good accuracy (Cα-RMSD: FNR = 2.8Å; Fd = 1.6 Å; IDDT: FNR = 77, Fd = 74 Å). Only small regions, mostly found on surface loops or terminal tails, show lower quality compared to the rest of the protein structures. The relative arrangement of the two proteins in the modelled complex is similar to that observed in the native structure, but the Fd subunit has a different orientation. The QS-score, which expresses the fraction of shared interface contacts between model and native complex, is 0.52. This value is higher than what we observe when comparing biologically functional FNR:Fd complexes, i.e. the templates and the native structure of the target protein, where the pairwise QS-scores range from 0.10 to 0.48, probably due to the elusive localization of the Fd moiety in crystals structures (69). Notably, the interface quality of the model, i.e. QS-score = 0.52, agrees very well with its estimated accuracy, i.e. QSQE = 0.54. The same correspondence is also observed between the local model quality estimated by QMEANDisCo and that measured according to all-atom IDDT score (Pearson correlation on the full complex = 0.80; Supplementary Figure S2). Finally, the FNR:Fd complex model was compared to that obtained using our previous modelling engine, ProMod-II. An average improvement of 2.5 IDDT points per chain is obtained with ProMod3. This is consistent with the results of our performance evaluation on a benchmark dataset (Figure ​(Figure11).

An external file that holds a picture, illustration, etc. Object name is gky427fig2.jpg

Modelling example of the Ferrodoxin/Ferredoxin-NADP(+) Reductase hetero dimeric complex. (A) Decision tree of templates clustered according to their quaternary structure features: oligomeric state, stoichiometry, topology and interface similarity. Three different clusters are formed based on interface similarity between templates. (B) PPI fingerprint analysis of available template structures. The ratio between interface and surface residue entropy (interface conservation, y-axis) is reported as a function of evolutionary distance (sequence identity, x-axis). Templates corresponding to SMTL ID: 1ewy.1 (in blue) and SMTL ID: 1gaq.1 (in green) show the typical conservation pattern observed for biologically relevant interfaces, with stronger conservation signal in the sequence identity range between 40 and 60%. Considering also remote homologs (below 40% sequence identity), only the interface in template SMTL ID: 1gaq.1 is deemed as conserved (interface/surface conservation ratio below zero). Template corresponding to SMTL ID: 3w5u.1 (in red) displays an interface/surface conservation ratio close to zero, as observed in crystal contacts/artefacts. (C) Structure superposition of available templates. Each template is coloured according to same colouring scheme of Figure ​Figure2A2A and B. Templates corresponding to SMTL ID: 1ewy.1 (in blue) and 1gaq.1 (in green) show similar arrangement of FNR and Fd in the complex. Template SMTL ID: 3w5u.1 (red) shows a different localization of the Fd moiety. Cross-linked cysteines are shown in sticks. (D) Structure superposition between model and native structure of the root FNR:Fd complex. The model is coloured according to its local quality using a colour gradient from blue (high quality) to red (low quality) as measured by all-atom IDDT score. The native structure of the complex is shown in light gray.

CONCLUSIONS

Computational structural modelling methods have established themselves as a valuable complement to experimental structural biology efforts towards increasing our understanding of the protein universe and of its properties. In this endeavour, comparative modelling techniques have matured into fully automated pipelines, providing easy access to reliable 3D models and broadening the spectrum of users and applications of protein models. SWISS-MODEL pioneered the field of fully automated comparative modelling servers 25 years ago and it has been continuously developed and improved since then.

With the new version of SWISS-MODEL presented here, we aimed at extending the scope of automated homology modelling to address the modelling of protein assemblies by efficiently using the information on quaternary structures available in the PDB. The success of this approach clearly depends on the availability of homologous complexes that can be used as templates for modelling. As such, ongoing structural biology efforts leading to structures of macromolecular complexes being determined at unprecedented speed are tremendously beneficial for making our approach increasingly applicable and effective. An important aspect is the ability to handle ambiguous or conflicting information present in available structural data, which is crucial for the development of stable and fully automated pipelines. Here, we showed how our PPI fingerprint analysis and model quality estimates could provide additional criteria to improve the automatic identification of templates, which in turn results into more accurate models and a biologically meaningful representation of their oligomeric state. Finally, we introduced an improved modelling engine and increased the precision of model quality estimates, leading to more accurate models and realistic error estimates at the same time.

Supplementary Material

Supplementary Data

ACKNOWLEDGEMENTS

We would like to thank all SWISS-MODEL users for their continuous feedbacks and participation to our user surveys, which motivate and inspire the continuous development of SWISS-MODEL. We thank sciCORE at the University of Basel for providing computational resources and system administration support.

SUPPLEMENTARY DATA

Supplementary Data are available at NAR Online.

FUNDING

SIB Swiss Institute of Bioinformatics. Swiss Foundation for Excellence and Talent in Biomedical Research (to G.S.). Funding for open access charge: SIB Swiss Institute of Bioinformatics. 'Fellowship for Excellence’ international PhD program of the Biozentrum Basel (to M.B.).

Conflict of interest statement. None declared.

REFERENCES

1. Fuller J.C., Burgoyne N.J., Jackson R.M.. Predicting druggable binding sites at the protein-protein interface. Drug Discov. Today. 2009; 14:155–161. [PubMed] [Google Scholar]
2. Nim S., Jeon J., Corbi-Verge C., Seo M.H., Ivarsson Y., Moffat J., Tarasova N., Kim P.M.. Pooled screening for antiproliferative inhibitors of protein-protein interactions. Nat. Chem. Biol. 2016; 12:275–281. [PMC free article] [PubMed] [Google Scholar]
3. Dutta S., Berman H.M.. Large macromolecular complexes in the Protein Data Bank: a status report. Structure. 2005; 13:381–388. [PubMed] [Google Scholar]
4. Marsh J.A., Teichmann S.A.. Structure, dynamics, assembly, and evolution of protein complexes. Annu. Rev. Biochem. 2015; 84:551–575. [PubMed] [Google Scholar]
5. Tramontano A. The computational prediction of protein assemblies. Curr. Opin. Struct. Biol. 2017; 46:170–175. [PubMed] [Google Scholar]
6. Pazos F., Helmer-Citterich M., Ausiello G., Valencia A.. Correlated mutations contain information about protein-protein interaction. J. Mol. Biol. 1997; 271:511–523. [PubMed] [Google Scholar]
7. Hopf T.A., Scharfe C.P., Rodrigues J.P., Green A.G., Kohlbacher O., Sander C., Bonvin A.M., Marks D.S.. Sequence co-evolution gives 3D contacts and structures of protein complexes. eLife. 2014; 3:e03430. [PMC free article] [PubMed] [Google Scholar]
8. Morris G.M., Lim-Wilby M.. Molecular docking. Methods Mol. Biol. 2008; 443:365–382. [PubMed] [Google Scholar]
9. Chaudhury S., Berrondo M., Weitzner B.D., Muthu P., Bergman H., Gray J.J.. Benchmarking and analysis of protein docking performance in Rosetta v3.2. PLoS One. 2011; 6:e22477. [PMC free article] [PubMed] [Google Scholar]
10. Kurkcuoglu Z., Koukos P.I., Citro N., Trellet M.E., Rodrigues J., Moreira I.S., Roel-Touris J., Melquiond A.S.J., Geng C., Schaarschmidt J. et al.. Performance of HADDOCK and a simple contact-based protein-ligand binding affinity predictor in the D3R Grand Challenge 2. J. Computer-aided Mol. Des. 2018; 32:175–185. [PMC free article] [PubMed] [Google Scholar]
11. Peterson L.X., Togawa Y., Esquivel-Rodriguez J., Terashi G., Christoffer C., Roy A., Shin W.H., Kihara D.. Modeling the assembly order of multimeric heteroprotein complexes. PLoS Comput. Biol. 2018; 14:e1005937. [PMC free article] [PubMed] [Google Scholar]
12. Janin J. Assessing predictions of protein-protein interaction: the CAPRI experiment. Protein Sci. 2005; 14:278–283. [PMC free article] [PubMed] [Google Scholar]
13. Janin J. The targets of CAPRI rounds 20–27. Proteins. 2013; 81:2075–2081. [PubMed] [Google Scholar]
14. Rodrigues J.P., Karaca E., Bonvin A.M.. Information-driven structural modelling of protein-protein interactions. Methods Mol. Biol. 2015; 1215:399–424. [PubMed] [Google Scholar]
15. Geng C., Narasimhan S., Rodrigues J.P., Bonvin A.M.. Information-driven, ensemble flexible peptide docking using HADDOCK. Methods Mol. Biol. 2017; 1561:109–138. [PubMed] [Google Scholar]
16. Peterson L.X., Shin W.H., Kim H., Kihara D.. Improved performance in CAPRI round 37 using LZerD docking and template-based modeling with combined scoring functions. Proteins. 2017; 86(Suppl. 1):311–320. [PMC free article] [PubMed] [Google Scholar]
17. Zhang Q.C., Petrey D., Norel R., Honig B.H.. Protein interface conservation across structure space. Proc. Natl. Acad. Sci. U.S.A. 2010; 107:10896–10901. [PMC free article] [PubMed] [Google Scholar]
18. Kundrotas P.J., Zhu Z., Janin J., Vakser I.A.. Templates are available to model nearly all complexes of structurally characterized proteins. Proc. Natl. Acad. Sci. U.S.A. 2012; 109:9438–9441. [PMC free article] [PubMed] [Google Scholar]
19. Matthews L.R., Vaglio P., Reboul J., Ge H., Davis B.P., Garrels J., Vincent S., Vidal M.. Identification of potential interaction networks using sequence-based searches for conserved protein-protein interactions or “interologs”. Genome Res. 2001; 11:2120–2126. [PMC free article] [PubMed] [Google Scholar]
20. Aloy P., Ceulemans H., Stark A., Russell R.B.. The relationship between sequence and interaction divergence in proteins. J. Mol. Biol. 2003; 332:989–998. [PubMed] [Google Scholar]
21. Szilagyi A., Zhang Y.. Template-based structure modeling of protein-protein interactions. Curr. Opin. Struct. Biol. 2014; 24:10–23. [PMC free article] [PubMed] [Google Scholar]
22. Yu H., Luscombe N.M., Lu H.X., Zhu X., Xia Y., Han J.D., Bertin N., Chung S., Vidal M., Gerstein M.. Annotation transfer between genomes: protein-protein interologs and protein-DNA regulogs. Genome Res. 2004; 14:1107–1118. [PMC free article] [PubMed] [Google Scholar]
23. Bertoni M., Kiefer F., Biasini M., Bordoli L., Schwede T.. Modeling protein quaternary structure of homo- and hetero-oligomers beyond binary interactions by homology. Scientific Rep. 2017; 7:10480. [PMC free article] [PubMed] [Google Scholar]
24. Lafita A., Bliven S., Kryshtafovych A., Bertoni M., Monastyrskyy B., Duarte J.M., Schwede T., Capitani G.. Assessment of protein assembly prediction in CASP12. Proteins. 2018; 86(Suppl. 1):247–256. [PMC free article] [PubMed] [Google Scholar]
25. Peitsch M.C. ProMod and Swiss-Model: Internet-based tools for automated comparative protein modelling. Biochem. Soc. Trans. 1996; 24:274–279. [PubMed] [Google Scholar]
26. Guex N., Peitsch M.C.. SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis. 1997; 18:2714–2723. [PubMed] [Google Scholar]
27. Schwede T., Kopp J., Guex N., Peitsch M.C.. SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res. 2003; 31:3381–3385. [PMC free article] [PubMed] [Google Scholar]
28. Arnold K., Bordoli L., Kopp J., Schwede T.. The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006; 22:195–201. [PubMed] [Google Scholar]
29. Bordoli L., Schwede T.. Automated protein structure modeling with SWISS-MODEL Workspace and the Protein Model Portal. Methods Mol. Biol. 2012; 857:107–136. [PMC free article] [PubMed] [Google Scholar]
30. Biasini M., Bienert S., Waterhouse A., Arnold K., Studer G., Schmidt T., Kiefer F., Gallo Cassarino T., Bertoni M., Bordoli L. et al.. SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res. 2014; 42:W252–W258. [PMC free article] [PubMed] [Google Scholar]
31. Benkert P., Biasini M., Schwede T.. Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics. 2011; 27:343–350. [PMC free article] [PubMed] [Google Scholar]
32. Haas J., Barbato A., Behringer D., Studer G., Roth S., Bertoni M., Mostaguir K., Gumienny R., Schwede T.. Continuous automated model evaluation (CAMEO) complementing the critical assessment of structure prediction in CASP12. Proteins. 2018; 86(Suppl. 1):387–398. [PMC free article] [PubMed] [Google Scholar]
33. Berman H.M., Battistuz T., Bhat T.N., Bluhm W.F., Bourne P.E., Burkhardt K., Feng Z., Gilliland G.L., Iype L., Jain S. et al.. The Protein Data Bank. Acta Crystallogr. D, Biol. Crystallogr. 2002; 58:899–907. [PubMed] [Google Scholar]
34. The UniProt, C. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2017; 45:D158–D169. [PMC free article] [PubMed] [Google Scholar]
35. Altschul S.F., Madden T.L., Schaffer A.A., Zhang J., Zhang Z., Miller W., Lipman D.J.. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997; 25:3389–3402. [PMC free article] [PubMed] [Google Scholar]
36. Camacho C., Coulouris G., Avagyan V., Ma N., Papadopoulos J., Bealer K., Madden T.L.. BLAST+: architecture and applications. BMC Bioinformatics. 2009; 10:421. [PMC free article] [PubMed] [Google Scholar]
37. Remmert M., Biegert A., Hauser A., Soding J.. HHblits: lightning-fast iterative protein sequence searching by HMM-HMM alignment. Nat. Methods. 2011; 9:173–175. [PubMed] [Google Scholar]
38. Biasini M., Schmidt T., Bienert S., Mariani V., Studer G., Haas J., Johner N., Schenk A.D., Philippsen A., Schwede T.. OpenStructure: an integrated software framework for computational structural biology. Acta Crystallogr. D, Biol. Crystallogr. 2013; 69:701–709. [PMC free article] [PubMed] [Google Scholar]
39. Bienert S., Waterhouse A., de Beer T.A., Tauriello G., Studer G., Bordoli L., Schwede T.. The SWISS-MODEL Repository—new features and functionality. Nucleic Acids Res. 2017; 45:D313–D319. [PMC free article] [PubMed] [Google Scholar]
40. de Beer T.A., Berka K., Thornton J.M., Laskowski R.A.. PDBsum additions. Nucleic Acids Res. 2014; 42:D292–D296. [PMC free article] [PubMed] [Google Scholar]
41. Velankar S., van Ginkel G., Alhroub Y., Battle G.M., Berrisford J.M., Conroy M.J., Dana J.M., Gore S.P., Gutmanas A., Haslam P. et al.. PDBe: improved accessibility of macromolecular structure data from PDB and EMDB. Nucleic Acids Res. 2016; 44:D385–D395. [PMC free article] [PubMed] [Google Scholar]
42. Dawson N.L., Lewis T.E., Das S., Lees J.G., Lee D., Ashford P., Orengo C.A., Sillitoe I.. CATH: an expanded resource to predict protein function through structure and sequence. Nucleic Acids Res. 2017; 45:D289–D295. [PMC free article] [PubMed] [Google Scholar]
43. Grosdidier A., Zoete V., Michielin O.. SwissDock, a protein-small molecule docking web service based on EADock DSS. Nucleic Acids Res. 2011; 39:W270–W277. [PMC free article] [PubMed] [Google Scholar]
44. Marcatili P., Olimpieri P.P., Chailyan A., Tramontano A.. Antibody modeling using the prediction of immunoglobulin structure (PIGS) web server [corrected]. Nat. Protoc. 2014; 9:2771–2783. [PubMed] [Google Scholar]
45. Messih M.A., Lepore R., Marcatili P., Tramontano A.. Improving the accuracy of the structure prediction of the third hypervariable loop of the heavy chains of antibodies. Bioinformatics. 2014; 30:2733–2740. [PMC free article] [PubMed] [Google Scholar]
46. Lepore R., Olimpieri P.P., Messih M.A., Tramontano A.. PIGSPro: prediction of immunoGlobulin structures v2. Nucleic Acids Res. 2017; 45:W17–W23. [PMC free article] [PubMed] [Google Scholar]
47. Chothia C., Lesk A.M.. Canonical structures for the hypervariable regions of immunoglobulins. J. Mol. Biol. 1987; 196:901–917. [PubMed] [Google Scholar]
48. Morea V., Tramontano A., Rustici M., Chothia C., Lesk A.M.. Conformations of the third hypervariable region in the VH domain of immunoglobulins. J. Mol. Biol. 1998; 275:269–294. [PubMed] [Google Scholar]
49. Tramontano A., Chothia C., Lesk A.M.. Framework residue 71 is a major determinant of the position and conformation of the second hypervariable region in the VH domains of immunoglobulins. J. Mol. Biol. 1990; 215:175–182. [PubMed] [Google Scholar]
50. Rose A.S., Hildebrand P.W.. NGL Viewer: a web application for molecular visualization. Nucleic Acids Res. 2015; 43:W576–W579. [PMC free article] [PubMed] [Google Scholar]
51. Sali A., Blundell T.L.. Comparative protein modelling by satisfaction of spatial restraints. J. Mol. Biol. 1993; 234:779–815. [PubMed] [Google Scholar]
52. Shapovalov M.V., Dunbrack R.L. Jr. A smoothed backbone-dependent rotamer library for proteins derived from adaptive kernel density estimates and regressions. Structure. 2011; 19:844–858. [PMC free article] [PubMed] [Google Scholar]
53. Xu J. Miyano S, Mesirov J, Kasif S, Istrail S, Pevzner PA, Waterman M. Rapid protein Side-Chain packing via tree decomposition. Research in Computational Molecular Biology: 9th Annual International Conference, RECOMB 2005, Cambridge, MA, USA, May 14–18, 2005. Proceedings. 2005; Berlin, Heidelberg: Springer; 423–439. [Google Scholar]
54. Krivov G.G., Shapovalov M.V., Dunbrack R.L. Jr. Improved prediction of protein side-chain conformations with SCWRL4. Proteins. 2009; 77:778–795. [PMC free article] [PubMed] [Google Scholar]
55. Eastman P., Swails J., Chodera J.D., McGibbon R.T., Zhao Y., Beauchamp K.A., Wang L.P., Simmonett A.C., Harrigan M.P., Stern C.D. et al.. OpenMM 7: rapid development of high performance algorithms for molecular dynamics. PLoS Comput. Biol. 2017; 13:e1005659. [PMC free article] [PubMed] [Google Scholar]
56. Mackerell A.D. Jr, Feig M., Brooks C.L. 3rd.. Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations. J. Comput. Chem. 2004; 25:1400–1415. [PubMed] [Google Scholar]
57. Mariani V., Biasini M., Barbato A., Schwede T.. lDDT: a local superposition-free score for comparing protein structures and models using distance difference tests. Bioinformatics. 2013; 29:2722–2728. [PMC free article] [PubMed] [Google Scholar]
58. Read R.J., Chavali G.. Assessment of CASP7 predictions in the high accuracy template-based modeling category. Proteins. 2007; 69(Suppl. 8):27–37. [PubMed] [Google Scholar]
59. Xu J., Zhang Y.. How significant is a protein structure similarity with TM-score = 0.5. Bioinformatics. 2010; 26:889–895. [PMC free article] [PubMed] [Google Scholar]
60. Elcock A.H., McCammon J.A.. Identification of protein oligomerization states by analysis of interface conservation. Proc. Natl. Acad. Sci. U.S.A. 2001; 98:2990–2994. [PMC free article] [PubMed] [Google Scholar]
61. Capra J.A., Singh M.. Predicting functionally important residues from sequence conservation. Bioinformatics. 2007; 23:1875–1882. [PubMed] [Google Scholar]
62. Schwede T., Sali A., Honig B., Levitt M., Berman H.M., Jones D., Brenner S.E., Burley S.K., Das R., Dokholyan N.V. et al.. Outcome of a workshop on applications of protein models in biomedical research. Structure. 2009; 17:151–159. [PMC free article] [PubMed] [Google Scholar]
63. Kim D.E., Chivian D., Baker D.. Protein structure prediction and analysis using the Robetta server. Nucleic Acids Res. 2004; 32:W526–W531. [PMC free article] [PubMed] [Google Scholar]
64. Yang Y., Faraggi E., Zhao H., Zhou Y.. Improving protein fold recognition and template-based modeling by employing probabilistic-based matching between predicted one-dimensional structural properties of query and corresponding native properties of templates. Bioinformatics. 2011; 27:2076–2082. [PMC free article] [PubMed] [Google Scholar]
65. Kallberg M., Wang H., Wang S., Peng J., Wang Z., Lu H., Xu J.. Template-based protein structure modeling using the RaptorX web server. Nat. Protoc. 2012; 7:1511–1522. [PMC free article] [PubMed] [Google Scholar]
66. Aliverti A., Faber R., Finnerty C.M., Ferioli C., Pandini V., Negri A., Karplus P.A., Zanetti G.. Biochemical and crystallographic characterization of ferredoxin-NADP(+) reductase from nonphotosynthetic tissues. Biochemistry. 2001; 40:14501–14508. [PubMed] [Google Scholar]
67. Kurisu G., Kusunoki M., Katoh E., Yamazaki T., Teshima K., Onda Y., Kimata-Ariga Y., Hase T.. Structure of the electron transfer complex between ferredoxin and ferredoxin-NADP(+) reductase. Nat. Struct. Biol. 2001; 8:117–121. [PubMed] [Google Scholar]
68. Shinohara F., Kurisu G., Hanke G., Bowsher C., Hase T., Kimata-Ariga Y.. Structural basis for the isotype-specific interactions of ferredoxin and ferredoxin: NADP(+) oxidoreductase: an evolutionary switch between photosynthetic and heterotrophic assimilation. Photosynth. Res. 2017; 134:281–289. [PubMed] [Google Scholar]
69. Morales R., Kachalova G., Vellieux F., Charon M.H., Frey M.. Crystallographic studies of the interaction between the ferredoxin-NADP+ reductase and ferredoxin from the cyanobacterium Anabaena: looking for the elusive ferredoxin molecule. Acta Crystallogr. D, Biol. Crystallogr. 2000; 56:1408–1412. [PubMed] [Google Scholar]
70. Duarte J.M., Srebniak A., Scharer M.A., Capitani G.. Protein interface classification by evolutionary analysis. BMC Bioinformatics. 2012; 13:334. [PMC free article] [PubMed] [Google Scholar]
71. Kimata-Ariga Y., Kubota-Kawai H., Lee Y.H., Muraki N., Ikegami T., Kurisu G., Hase T.. Concentration-dependent oligomerization of cross-linked complexes between ferredoxin and ferredoxin-NADP+ reductase. Biochem. Biophys. Res. Commun. 2013; 434:867–872. [PubMed] [Google Scholar]
Articles from Nucleic Acids Research are provided here courtesy of Oxford University Press
-