SWISS-MODEL Homology Modelling Report

Model Building Report

This document lists the results for the homology modelling project "B2UTZ3" submitted to SWISS-MODEL workspace on Oct. 26, 2017, 8:51 p.m..The submitted primary amino acid sequence is given in Table T1.

If you use any results in your research, please cite the relevant publications:

Marco Biasini; Stefan Bienert; Andrew Waterhouse; Konstantin Arnold; Gabriel Studer; Tobias Schmidt; Florian Kiefer; Tiziano Gallo Cassarino; Martino Bertoni; Lorenza Bordoli; Torsten Schwede. (2014). SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Research (1 July 2014) 42 (W1): W252-W258; doi: 10.1093/nar/gku340.
Arnold, K., Bordoli, L., Kopp, J. and Schwede, T. (2006) The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics, 22, 195-201.
Benkert, P., Biasini, M. and Schwede, T. (2011) Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics, 27, 343-350

Results

The SWISS-MODEL template library (SMTL version 2017-10-23, PDB release 2017-10-13) was searched with Blast (Altschul et al., 1997) and HHBlits (Remmert, et al., 2011) for evolutionary related structures matching the target sequence in Table T1. For details on the template search, see Materials and Methods. Overall 12 templates were found (Table T2).

Models

The following model was built (see Materials and Methods "Model Building"):

Model #01

FileBuilt withOligo-StateLigandsGMQEQMEAN
PDB ProMod3 Version 1.1.0. homo-dimer (matching prediction)
None
0.73-4.78
QMEAN-4.78
-3.63
All Atom-1.63
Solvation1.64
Torsion-4.26
TemplateSeq IdentityOligo-stateFound byMethodResolutionSeq SimilarityRangeCoverageDescription
5ckr.1.B41.47homo-dimerBLASTX-ray2.95Å0.40 16 - 353 0.96Phospho-N-acetylmuramoyl-pentapeptide-transferase
LigandAdded to ModelDescription
57M✕ - Binding site not conserved.
Muraymycin D2
57M✕ - Binding site not conserved.
Muraymycin D2

Target    MLYSLLYGYFNINLFQYLTFRAGLGFFIAFFLTLFLMPKFILWAKAKKANQPISSFV-----PSHQNKKDTPTMGGIVFV
5ckr.1.B --------WFAFNVLKYITFRSFTAVLIAFFLTLVLSPSFI--NRLRKIQRLFGGYVREYTPESHEVKKYTPTMGGIVIL

Target FATIVASVLCASLGNLYVLLGMIVLVGFSFVGFKDDYTKINQQNNAGMSVKMKFGM-LFVLSLVVSVLLSLKGLDTFLYA
5ckr.1.B IVVTLSTLLLMRWDIKYTWVVLLSFLSFGTIGFWDDYVKL--KNKKGISIKTKFLLQVLSASLISVLIYYWADIDTILYF

Target PFLKNPLFEMPTILAVGFWVLVFLSTSNAVNLTDGLDGLASVPSIFTLLSLSIFVYVAGNAEFSKYLLYPKVIDVGELFV
5ckr.1.B PFFKELYVDL-GVLYLPFAVFVIVGSANAVNLTDGLDGLAIGPAMTTATALGVVAYAVGHSKIAQYLNIPYVPYAGELTV

Target ISLALVGSLFGFLWYNCNPASVFMGDSGSLALGGFIAYNAIVSHNEILLVLMGLIFVVETLSVILQVGSYK-TRKKRLFL
5ckr.1.B FCFALVGAGLGFLWFNSFPAQMFMGDVGSLSIGASLATVALLTKSEFIFAVAAGVFVFETISVILQIIYFRWTGGKRLFK

Target MAPIHHHFEQKGWAENKVIVRFWIISMLSNLVALLSLKVR
5ckr.1.B RAPFHHHLELNGLPEPKIVVRMWIISILLAIIAISMLKLR




Materials and Methods

Template Search

Template search with Blast and HHBlits has been performed against the SWISS-MODEL template library (SMTL, last update: 2017-10-23, last included PDB release: 2017-10-13).

The target sequence was searched with BLAST (Altschul et al., 1997) against the primary amino acid sequence contained in the SMTL. A total of 4 templates were found.

An initial HHblits profile has been built using the procedure outlined in (Remmert, et al., 2011), followed by 1 iteration of HHblits against NR20. The obtained profile has then be searched against all profiles of the SMTL. A total of 18 templates were found.

Template Selection

For each identified template, the template's quality has been predicted from features of the target-template alignment. The templates with the highest quality have then been selected for model building.

Model Building

Models are built based on the target-template alignment using ProMod3. Coordinates which are conserved between the target and the template are copied from the template to the model. Insertions and deletions are remodelled using a fragment library. Side chains are then rebuilt. Finally, the geometry of the resulting model is regularized by using a force field. In case loop modelling with ProMod3 fails, an alternative model is built with PROMOD-II (Guex, et al., 1997).

Model Quality Estimation

The global and per-residue model quality has been assessed using the QMEAN scoring function (Benkert, et al., 2011) . For improved performance, weights of the individual QMEAN terms have been trained specifically for SWISS-MODEL.

Ligand Modelling

Ligands present in the template structure are transferred by homology to the model when the following criteria are met (Gallo -Casserino, to be published): (a) The ligands are annotated as biologically relevant in the template library, (b) the ligand is in contact with the model, (c) the ligand is not clashing with the protein, (d) the residues in contact with the ligand are conserved between the target and the template. If any of these four criteria is not satisfied, a certain ligand will not be included in the model. The model summary includes information on why and which ligand has not been included.

Oligomeric State Conservation

Homo-oligomeric structure of the target protein is predicted based on the analysis of pairwise interfaces of the identified template structures. For each relevant interface between polypeptide chains (interfaces with more than 10 residue-residue interactions), the QscoreOligomer (Mariani et al., 2011) is predicted from features such as similarity to target and frequency of observing this interface in the identified templates (Kiefer, Bertoni, Biasini, to be published). The prediction is performed with a random forest regressor using these features as input parameters to predict the probability of conservation for each interface. The QscoreOligomer of the whole complex is then calculated as the weight-averaged QscoreOligomer of the interfaces. The oligomeric state of the target is predicted to be the same as in the template when QscoreOligomer is predicted to be higher or equal to 0.5.

References

Altschul, S.F., Madden, T.L., Schaffer, A.A., Zhang, J., Zhang, Z., Miller, W. and Lipman, D.J. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res, 25, 3389-3402.
Remmert, M., Biegert, A., Hauser, A. and Soding, J. (2012) HHblits: lightning-fast iterative protein sequence searching by HMM-HMM alignment. Nat Methods, 9, 173-175.
Guex, N. and Peitsch, M.C. (1997) SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis, 18, 2714-2723.
Sali, A. and Blundell, T.L. (1993) Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol, 234, 779-815.
Benkert, P., Biasini, M. and Schwede, T. (2011) Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics, 27, 343-350.
Mariani, V., Kiefer, F., Schmidt, T., Haas, J. and Schwede, T. (2011) Assessment of template based protein structure predictions in CASP9. Proteins, 79 Suppl 10, 37-58.

Table T1:

Primary amino acid sequence for which templates were searched and models were built.

MLYSLLYGYFNINLFQYLTFRAGLGFFIAFFLTLFLMPKFILWAKAKKANQPISSFVPSHQNKKDTPTMGGIVFVFATIVASVLCASLGNLYVLLGMIVL
VGFSFVGFKDDYTKINQQNNAGMSVKMKFGMLFVLSLVVSVLLSLKGLDTFLYAPFLKNPLFEMPTILAVGFWVLVFLSTSNAVNLTDGLDGLASVPSIF
TLLSLSIFVYVAGNAEFSKYLLYPKVIDVGELFVISLALVGSLFGFLWYNCNPASVFMGDSGSLALGGFIAYNAIVSHNEILLVLMGLIFVVETLSVILQ
VGSYKTRKKRLFLMAPIHHHFEQKGWAENKVIVRFWIISMLSNLVALLSLKVR

Table T2:

TemplateSeq IdentityOligo-stateFound byMethodResolutionSeq SimilarityCoverageDescription
5ckr.1.B40.86homo-dimerHHblitsX-ray2.95Å0.400.99Phospho-N-acetylmuramoyl-pentapeptide-transferase
4j72.1.A40.86homo-dimerHHblitsX-ray3.30Å0.400.99Phospho-N-acetylmuramoyl-pentapeptide-transferase
4j72.1.B40.86homo-dimerHHblitsX-ray3.30Å0.400.99Phospho-N-acetylmuramoyl-pentapeptide-transferase
5ckr.1.B41.47homo-dimerBLASTX-ray2.95Å0.400.96Phospho-N-acetylmuramoyl-pentapeptide-transferase
4j72.1.A41.47homo-dimerBLASTX-ray3.30Å0.400.96Phospho-N-acetylmuramoyl-pentapeptide-transferase
4j72.1.B41.47homo-dimerBLASTX-ray3.30Å0.400.96Phospho-N-acetylmuramoyl-pentapeptide-transferase
5jnq.1.A36.83homo-dimerHHblitsX-ray2.60Å0.370.89Phospho-N-acetylmuramoyl-pentapeptide-transferase
5jnq.1.A40.68homo-dimerBLASTX-ray2.60Å0.400.84Phospho-N-acetylmuramoyl-pentapeptide-transferase
5lev.1.A16.17monomerHHblitsX-ray3.20Å0.290.86UDP-N-acetylglucosamine--dolichyl-phosphate N-acetylglucosaminephosphotransferase
2l8s.1.A20.00monomerHHblitsNMRNA0.340.04Integrin alpha-1
2kkz.1.A12.50monomerHHblitsNMRNA0.250.05Non-structural protein NS1
3few.1.A6.25monomerHHblitsX-ray2.45Å0.200.05Colicin S4