Skip to main content

GP0.4 from bacteriophage T7: in silico characterisation of its structure and interaction with E. coli FtsZ



Proteins produced by bacteriophages can have potent antimicrobial activity. The study of phage-host interactions can therefore inform small molecule drug discovery by revealing and characterising new drug targets. Here we characterise in silico the predicted interaction of gene protein 0.4 (GP0.4) from the Escherichia coli (E. coli) phage T7 with E. coli filamenting temperature-sensitive mutant Z division protein (FtsZ). FtsZ is a tubulin homolog which plays a key role in bacterial cell division and that has been proposed as a drug target.


Using ab initio, fragment assembly structure modelling, we predicted the structure of GP0.4 with two programs. A structure similarity-based network was used to identify a U-shaped helix-turn-helix candidate fold as being favoured. ClusPro was used to dock this structure prediction to a homology model of E. coli FtsZ resulting in a favourable predicted interaction mode. Alternative docking methods supported the proposed mode which offered an immediate explanation for the anti-filamenting activity of GP0.4. Importantly, further strong support derived from a previously characterised insertion mutation, known to abolish GP0.4 activity, that is positioned in close proximity to the proposed GP0.4/FtsZ interface.


The mode of interaction predicted by bioinformatics techniques strongly suggests a mechanism through which GP0.4 inhibits FtsZ and further establishes the latter’s druggable intrafilament interface as a potential drug target.


Resistance to antibiotics has become a source of great public concern recently, compounded by the slow emergence of new antibiotics [1]. One possible solution may be phage therapy. Bacteriophages (referred to hereafter as phages) are among the simplest and most abundant types of microorganisms [2], and have been used therapeutically for close to a century. With the advent of antibiotics, the use of phages therapeutically has decreased in the western world. However, as resistance to antibiotics becomes a pressing issue, phages are once again being seen as a way to combat disease [3, 4].

There are several advantages to the therapeutic use of phages; they affect highly specific targets that minimise collateral damage [5], they can target bacterial strains resistant to antibiotics, and they can be used to supplement the effects of antibiotics [6]. Phages are, however, not without problems. For example, the delivery of phages is a major challenge, with the immune system presenting a large hurdle [3] and the use of phages as a delivery method requires accurate diagnosis of the disease-causing bacteria and therefore slows down treatment.

Study of phage-host interaction can also inform small molecule drug discovery by revealing new drug targets and pinpointing their weaknesses. Proteins in phages have naturally evolved to find effective methods to disrupt bacteria. Furthermore, in multi-protein complexes involved in processes such as cell regulation, small numbers of amino acids form hotspots which contribute most of the free energy during interactions [7]. By studying how phage proteins disrupt the protein–protein interfaces, we can identify potential hotspots and target them when designing drugs. Such drugs would dispense with the specificity limitation of the original phage since they can be designed to target broadly conserved bacterial mechanisms potentially rendering unnecessary the diagnosis of the specific pathogenic bacterium [8, 9].

Upon infection of bacteria, phages take over the host resources through the actions of proteins expressed early in the infection. One such protein from the Escherichia coli phage T7 is gene protein (GP) 0.4. GP0.4 directly inhibits the filamenting temperature-sensitive mutant Z division protein (FtsZ) of E. coli by preventing its assembly into protofilaments both in vivo and in vitro [10].

FtsZ is a tubulin homologue that plays a key role in the division of bacteria cells [1113]. Much like tubulin, purified FtsZ binds and hydrolyses GTP [13]. GTP binding induces the FtsZ to polymerise into one of two polar protofilament conformations; straight or gently curved [1115]. Between FtsZ monomers, an active site is formed which hydrolyses the GTP, and remains accessible to the GTP-rich cytoplasm [13, 15]. Therefore GTP binding can be rapidly restored allowing protofilaments to consist of mostly FtsZ-GTP subunits resistant to depolymerisation [13]. The FtsZ protofilaments associate laterally to form bundles or sheets [11, 13]. As the protofilaments bundle together they form the Z ring at the site of cytokinesis. Once assembled, the Z ring plays a crucial role in recruiting downstream proteins essential for cell division [12, 13]. Therefore the inhibition of FtsZ polymerisation prevents the division of the bacteria [10].

The exact mechanism through which GP0.4 interacts with FtsZ is unknown. However, a mutation in FtsZ that introduced a six nucleotide insertion (TCGGCG) overcame the toxicity of GP0.4 [10]. Here, using a suite of structural bioinformatics methods, we predict the structure of GP0.4 ab initio and determine a mode of interaction with FtsZ in accord with available data. The results suggest that the FtsZ protofilament interface is targeted in different ways by phage proteins for antimicrobial purposes or by the bacteria’s own proteins for regulatory purposes. These results add weight to the notion that this pocket is a druggable interface.


Sequence analysis

BLAST [16] was used with the GP0.4 and FtsZ sequences to identify homologous proteins in the UniProt database [17] and to search for any homologues already structurally characterised in the Protein Data Bank (PDB) [18]. The HHPred server was used to confirm that the structure of GP0.4 could not be inferred by distant homology to any known structure.

GP0.4 ab initio modelling

To elucidate a structure for GP0.4, ROSETTA 3.5 AbinitioRelax modelling [1924] was used. The AbinitioRelax application consists of two steps; the first step is a coarse-grained fragment-based search for conformation that uses a score function which favours protein-like features (Abinitio). The second step is an all atom refinement using the Rosetta full-atom force field (Relax). The Robetta server [22, 25] was used to generate the required fragments and the PSIPRED [26] secondary structure predictions used by ROSETTA.

Using ROSETTA, 10,000 ab initio models (or decoys) were produced for GP0.4 and clustered using default protocols on a Linux workstation. Representatives of the largest ten resulting clusters were considered as candidate fold predictions. This process was repeated for four homologues identified in the BLAST search. In addition to ROSETTA, GP0.4 and each homologous sequence was submitted to the QUARK ab initio server [27]. QUARK is an ab initio modelling method conceptually similar to ROSETTA. The server yielded a further ten models (representatives of the ten largest clusters) for each of the proteins.

Validation and comparison of GP0.4 models

GROMACS 5.0.1 molecular dynamics [2832] was used to test the stability of the models over a period of 5 ns using a cubic box filled with water as the solvent, chloride as the counterions, and the AMBER99SB-ILDN force field [33]. ProSA [34, 35] and QMEAN [3639] were used to obtain protein structure quality measurements.

The top ten models produced for each protein by ROSETTA and QUARK were submitted to the ProCKSI comparative server [40] to assess any structural similarity between them by producing a matrix of template modelling (TM) scores [41]. This matrix was visualised using CLuster ANalysis of Sequences (CLANS) [42] to cluster the models with an attraction value of >0.6.

FtsZ comparative modelling

ROSETTA comparative modelling (RosettaCM) [25] was used to make ten models of E. coli FtsZ based on ten homologues (2VAW, 2VXY, 4M8I, 1RQ2, 1OFU, 1W5F, 2VAP, 2R75, 3J4S and 4B45) obtained from HHPred [4345]. RosettaCM optimises an all-atom energy function over the conformational space defined by homologue structures to produce models with more accurate side chain and backbone conformations than previously available. In order to select the best model produced, the inbuilt RosettaCM scoring system was used in combination with the model quality assessment program ProSA [34, 35].


Four servers that performed well in the most recent CAPRI (Critical Assessment of PRediction of Interactions) competition [46] were used to predict how models of FtsZ and GP0.4 might interact—ClusPro [4750], Swarmdock [5153], Dock/PIERR [54] and GRAMM-X [55].

The final putative GP0.4 binding site was assessed for drugability with two complementary servers (DoGsiteScorer [56] and FTMap [5759]) and its conservation between FtsZ sequences quantified with ConSurf [6063].

Results and discussion

ab initio modelling of GP0.4

Since fragment-based ab initio modelling can benefit from consideration of homologous sequences, rather than focusing exclusively on the target, a BLAST [16] search with GP0.4 was done and identified four complete, non-redundant homologues in other phages (Fig. 1), none of which had been characterised.

Fig. 1

The alignment of GP0.4 (Uniprot accession: P03776), Yersinia phage YpsP-R (I6Q992), Salmonella phage Vi06 (E1XU80), and Citrobacter phage CR8 (W6PPJ7) using multiple sequence comparison by Log- Expectation (MUSCLE) [64, 65], and coloured using the Clustal colour scheme in Jalview [66]. The PSIPRED [67] secondary structure prediction for GP0.4 used in the ab initio modelling is shown below: red bars represent predicted α-helices and green arrows represent predicted β-strands

Fragment-based ab initio modelling was done locally using ROSETTA [1924] and remotely at the QUARK server [27] producing ten models (cluster centroids) each for GP0.4 itself (UniProt accession P03776), and homologues in Yersinia phage YpsP-R (I6Q992), Salmonella phage Vi06 (E1XU80), and Citrobacter phage CR8 (W6PPJ7). In order to look for consensus among the 100 model set, clustering was carried out based on pairwise TM scores obtained from the ProCKSI server [40]. TM scores are a well-established method to quantitatively measure the structural difference between two proteins [41] and range between 0 and 1: randomly chosen structures have an average score of 0.17 while a score of over 0.5 is considered as indicating that the two proteins have the same fold. Figure 2 shows GP0.4 and homologue models (nodes) linked by lines indicating a shared TM score of >0.6. This produced two clusters of note, one considerably larger than the other.

Fig. 2

Clustering by pairwise TM scores of models of GP0.4 (UniProt accession number: P03776), Yersinia phage YpsP-R (I6Q992), Salmonella phage Vi06 (E1XU80), and Citrobacter phage CR8 (W6PPJ7). Each point represents a model, coloured by origin blue (Rosetta) or red (QUARK). Shapes indicate the sequence modelled: square GP0.4 (P03776), triangle YpsP-R (I6Q992), circle Vi06 (E1XU80) and diamond CR8 (W6PPJ7). Lines indicate TM scores of >0.6. The figure was made using CLANS [42]

The larger cluster, containing ‘U’-shaped models was made up of QUARK predictions for GP0.4 and homologues from phages YpsP-R, Vi06 and CR8, along with Rosetta models for Vi06 and YpsP-R. The smaller and therefore less favoured cluster contained Rosetta models for GP0.4 and the CR8 homologue which were ‘V’-shaped. From each cluster, the GP0.4 model which made the most links to other structures was selected as a representative model for further analysis (Fig. 3). GROMACS [2832] molecular dynamics was used to test the stability of these models and found both stabilize at low structural deviation (<1.75 Å rmsd on Cα atoms) from the starting structure within 5 ns. Similarly, model quality assessment with ProSA [34, 35] and QMEAN [3639] did not favour one model or the other: within each cluster structures were found with ProSA scores as low as −4.36 and QMEAN scores as high as 0.734, within the range seen for experimental protein structures of a similar size. Furthermore, searches of the PDB using PDBeFOLD [68] and DALI [69] did not identify any significantly similar structures to either of the two GP0.4 models.

Fig. 3

Comparison of GP0.4 models from the major (a ‘U’-shaped) and minor (b ‘V’-shaped) clusters in Fig. 2. Each is coloured from N-(blue) to C-terminus (red) as a spectrum. Visualised in PyMOL Molecular Graphics System, Version 1.7.4 Schrödinger, LLC. PyMOL was also used for Figs.  4, 5, 6 and 7

FtsZ modelling

RosettaCM [25] was used to produce ten models of E.coli FtsZ using homologous structures from P. aeruginosa (1ofu, 2vaw), M. tuberculosis (1rq2), T. maritima (1w5f), A. aeolicus (2r75), M. jannaschii (2vap), B.subtilis (2vxy), B. thuringiensis (3j4 s), H. volcanii (4b45) and S. epidermidis (4m8i) as templates. FtsZ has a well conserved core domain followed by a variable C-terminal domain (Additional file 1: Figure S1). Indeed the ten models produced were very similar apart from that variable C-terminal domain which adopted a wide variety of poorly packed conformations. Previous studies of FtsZ–FtsZ interaction found that the C-terminal domain was not required for the assembly of protofilaments, but was essential for interaction with other membrane associated cell division proteins [13]. Since our focus here was on the inter-subunit interface of the protofilament targeted by GP0.4, the C-terminal region was eliminated and RosettaCM energy scores and ProSA [34, 35] scores were used to identify the most favoured model of the core domain.


Initial docking experiments between the GP0.4 ‘U’ and ‘V’ models and the FtsZ model were performed using the ClusPro server [4750]. The ‘U’ and the ‘V’ model were predicted to bind to two different locations using default (‘balanced’) settings. The docking site for the ‘U’ model was obtained from a cluster of 280/1000 models, and the docking site for the ‘V’ model was obtained from a cluster of 382/1000 models. The ‘V’ model bound to a region flanked by H6, H7 and S7, whereas the more favoured ‘U’ model bound to a region flanked by H1, H5 and H7 (Fig. 4).

Fig. 4

The top results when docking the ‘V’ form (a) and the ‘U’ form (b) of GP0.4 to the FtsZ model in ClusPro [4750]. The GP0.4 model is coloured blue while the FtsZ structure is coloured by secondary structure: red for α-helices, yellow for β-strands, green for coil. Regular secondary structure elements in FtsZ are labelled H for α-helices or S for β-strands. The α-helices (H) and β-strands (S) shown in the cartoon representation of FtsZ are labeled according to the 3D structures of Methanococcus jannaschi FtsZ

The results from ClusPro [4750] were compared to those obtained using Swarmdock [5153], Dock/PIERR [54] and GRAMM-X [55]. For the ‘U’ shaped model the top prediction from Dock/PIERR and the 6th highest prediction from Swarmdock were found to match the top prediction from ClusPro (Fig. 5). GRAMM-X also predicted that GP0.4 would bind to the same region of FtsZ but predicted a different docking pose (not shown).

Fig. 5

Cross-eyed stereo view of the consensus in docking of the ‘U’ form using ClusPro (Blue), Swarmdock (Magenta) and Dock/PIERR (Green). The FtsZ structure is coloured by secondary structure: red for alpha helix, yellow for beta sheet, green for coil

In contrast, there was little agreement between the docking server results for the ‘V’ shaped model. For example, there was no coincidence between any of the top ten predictions from ClusPro, Dock/PIERR or Swarmdock. The lack of a consensus in the ‘V’ shaped model docking results is further indirect evidence supporting the ‘U’-shaped model since consistency between docking results would be expected of reliable results. Further support for the docking predictions for the ‘U’-shaped model is obtained by considering the FtsZ protofilament. The ‘U’-shaped model inserts its well-conserved central helix-turn-helix motif into a cleft between helices 1, 5 and 7 on FtsZ thereby preventing formation of a further intrafilament interface (Fig. 6). The conserved FtsZ residues G21, N24, G47 and G107 contribute to the binding by making hydrogen bond interactions with GP0.4. Given that the consensus ‘U’-shaped model offered such a plausible explanation for the anti-filamenting activity of GP0.4, attention was then focussed exclusively on it, excluding the ‘V’-shaped complex from further consideration.

Fig. 6

a The interaction between an upper FtsZ (Red) monomer and a lower FtsZ (Green) monomer in a protofilament. b the locations at which GP0.4 (Blue) binds to the lower FtsZ (Green) monomer and MciZ (Purple) binds to the upper FtsZ (Red) monomer. These can be seen to target the same region of the interface

Mutant Ftsz modelling

Kiro et al. [10] identified a GP0.4 resistant variant of FtsZ, FtsZ9, where a 6 nucleotide insertion (TCGGCG) at position 18 resulted in a duplication of glycine–valine that allowed FtsZ9 to overcome the toxicity of GP0.4. It was reasoned that modelling this mutation might help identify how and where GP0.4 interacts with FtsZ. As shown in Fig. 7, the site of the mutation lies exactly at the base of the cleft to which GP0.4 is predicted to bind. This offers strong support to the binding model since the additional two residues encoded by FtsZ9 are thus perfectly placed to impeded GP0.4 binding and thereby confer resistance on the bacterium. The mutation also lies near the guanine nucleotide binding side of FtsZ, and the data of Kiro et al. [10] suggest that nucleotide binding is unaffected. However, we were unable to model a mutant FtsZ structure in which guanine nucleotides bound in the canonical fashion: conceivably, nucleotide binding is maintained in an altered fashion in FtsZ9.

Fig. 7

An enhanced view of the docked GP0.4. Interaction between GP0.4 (blue) and FtsZ residues N24 (2.8 Å), L178 (5.2 Å), and F182 (4.6 Å) are shown as lines. The orange circle indicates the location at where the mutation is inserted

Interestingly, MciZ, a developmental regulator of bacterial cell division was found to inhibit FtsZ polymerisation by targeting the same region of the intrafilament FtsZ interface as GP0.4 whilst binding to the opposite face of the FtsZ monomer [70], as shown in Fig. 6.

Characteristics of the predicted GP0.4 binding site

To help identify whether our putative GP0.4 binding pocket is a suitable candidate for small molecule drug design we used druggability servers; DoGsiteScorer [56] which looks for properties shared with known drug pockets and FTMap [5759] which identifies areas bound by small organic probe molecules. The druggability servers both identified our pocket as a suitable biological target for drug binding (Fig. 8).

Fig. 8

a The highest ranked drug pocket predicted by DoGsiteScorer (mauve) coincides with the putative GP0.4 binding site (GP0.4 surface shown as transparent grey). b The second highest ranked location of ligand binding predicted by FTMap (magenta) again coincides with the putative GP0.4 binding site (grey). In (a) and (b) FtsZ is coloured by secondary structure: red for alpha helix, yellow for beta sheet, green for coil. c ConSurf results indicating FtsZ amino acid conservation on a scale of blue to red, where blue is most conserved and red is least conserved

In order to assess whether the druggable predicted GP0.4 binding site is conserved among bacteria, we mapped sequence conservation of 150 FtsZ sequences onto the FtsZ structure using the ConSurf server. Figure 8 shows that the area the ‘U’ model bound to was well conserved. This provides evidence that a GP0.4 like protein or small molecule targeted to its binding site might be effective against FtsZ in other bacteria including pathogenic species: indeed, GP0.4 homologues were identified in phages for Yersinia and Salmonella supporting this idea.

The FtsZ filament interface targeted by GP0.4 has previously been highlighted as a possible target for small molecule targeting. However, there are concerns over targeting the GTP-binding site due to a risk of poisoning the eukaryotic homologue, tubulin [71]. Nevertheless, our results emphasise the druggable nature of a larger interface including the pocket targeted by GP0.4 They therefore encourage further efforts at exploiting the interface for small molecule drug design as well as offering possible routes forward for peptidomimetic inhibition.


Characterising the targets and inhibitory mechanisms of phage proteins is a valuable route to the discovery and validation of new potential drug targets. Here we bring a range of structural bioinformatics methods to bear on the interaction between phage GP0.4 and its bacterial target FtsZ. We provide evidence that GP0.4 adopts a ‘U’ shaped conformation that inserts into a cleft between helices 1, 5 and 7 on FtsZ. The hypothesis is strongly supported by data obtained for a GP0.4-resistant FtsZ mutation [10]. The presence of GP0.4 bound to this region, as shown in Fig. 6, would physically interfere with assembly of the FtsZ filament. The importance of this FtsZ–FtsZ interface was further demonstrated by MciZ, a regulatory protein that binds to the opposite side of the interface to inhibit bacterial cell division. The druggable nature of the broad intrafilament FtsZ interface should encourage future drug design.


  1. 1.

    Oyston PCF, Fox MA, Richards SJ, Clark GC. Novel peptide therapeutics for treatment of infections. J Med Microbiol. 2009;58(Pt 8):977–87.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Labrie SJ, Samson JE, Moineau S. Bacteriophage resistance mechanisms. Nat Rev Microbiol. 2010;8:317–27.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Lu TK, Koeris MS. The next generation of bacteriophage therapy. Curr Opin Microbiol. 2011;14:524–31.

    Article  PubMed  Google Scholar 

  4. 4.

    Pirnay J-P, De Vos D, Verbeken G, Merabishvili M, Chanishvili N, Vaneechoutte M, Zizi M, Laire G, Lavigne R, Huys I, Van den Mooter G, Buckling A, Debarbieux L, Pouillot F, Azeredo J, Kutter E, Dublanchet A, Górski A, Adamia R. The phage therapy paradigm: prêt-à-porter or sur-mesure? Pharm Res. 2011;28:934–7.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Keen EC. Phage therapy: concept to cure. Front Microbiol. 2012;3:238.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Kutateladze M, Adamia R. Bacteriophages as potential new therapeutics to replace or supplement antibiotics. Trends Biotechnol. 2010;28:591–5.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Bogan AA, Thorn KS. Anatomy of hot spots in protein interfaces. J Mol Biol. 1998;280:1–9.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Fischetti VA, Nelson D, Schuch R. Reinventing phage therapy: are the parts greater than the sum? Nat Biotechnol. 2006;24:1508–11.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Projan S. Phage-inspired antibiotics? Nat Biotechnol. 2004;22:167–8.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Kiro R, Molshanski-Mor S, Yosef I, Milam SL, Erickson HP, Qimron U. Gene product 0.4 increases bacteriophage T7 competitiveness by inhibiting host cell division. Proc Natl Acad Sci USA. 2013;110:19549–54.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Erickson HP, Taylor DW, Taylor KA, Bramhill D. Bacterial cell division protein FtsZ assembles into protofilament sheets and minirings, structural homologs of tubulin polymers. Proc Natl Acad Sci. 1996;93:519–23.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Loose M, Mitchison TJ. The bacterial cell division proteins FtsA and FtsZ self-organize into dynamic cytoskeletal patterns. Nat Cell Biol. 2014;16:38–46.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Margolin W. FtsZ and the division of prokaryotic cells and organelles. Nat Rev Mol Cell Biol. 2005;6:862–71.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Lu C, Reedy M, Erickson HP. Straight and curved conformations of FtsZ are regulated by GTP hydrolysis. J Bacteriol. 2000;182:164–70.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Mingorance J, Rivas G, Vélez M, Gómez-Puertas P, Vicente M. Strong FtsZ is with the force: mechanisms to constrict bacteria. Trends Microbiol. 2010;18:348–56.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    UniProt Consortium. UniProt: a hub for protein information. Nucleic Acids Res. 2015;43:204–12.

    Article  Google Scholar 

  18. 18.

    Berman HM. The protein data bank. Nucleic Acids Res. 2000;28:235–42.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Bonneau R, Tsai J, Ruczinski I, Chivian D, Rohl C, Strauss CE, Baker D. Rosetta in CASP4: progress in ab initio protein structure prediction. Proteins. 2001;5:119–26.

    Article  PubMed  Google Scholar 

  20. 20.

    Bonneau R, Strauss CEM, Rohl CA, Chivian D, Bradley P, Malmström L, Robertson T, Baker D. De novo prediction of three-dimensional structures for major protein families. J Mol Biol. 2002;322:65–78.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Bradley P, Misura KMS, Baker D. Toward high-resolution de novo structure prediction for small proteins. Science. 2005;309:1868–71.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Raman S, Vernon R, Thompson J, Tyka M, Sadreyev R, Pei J, Kim D, Kellogg E, DiMaio F, Lange O, Kinch L, Sheffler W, Kim B-H, Das R, Grishin NV, Baker D. Structure prediction for CASP8 with all-atom refinement using Rosetta. Proteins. 2009;77(Suppl 9):89–99.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Simons KT, Kooperberg C, Huang E, Baker D. Assembly of protein tertiary structures from fragments with similar local sequences using simulated annealing and Bayesian scoring functions. J Mol Biol. 1997;268:209–25.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Simons KT, Ruczinski I, Kooperberg C, Fox BA, Bystroff C, Baker D. Improved recognition of native-like protein structures using a combination of sequence-dependent and sequence-independent features of proteins. Proteins. 1999;34:82–95.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Song Y, DiMaio F, Wang RYR, Kim D, Miles C, Brunette T, Thompson J, Baker D. High-resolution comparative modeling with RosettaCM. Structure. 2013;21:1735–42.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Jones DT. Protein secondary structure prediction based on position-specific scoring matrices. J Mol Biol. 1999;292:195–202.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Xu D, Zhang Y. Ab initio protein structure assembly using continuous structure fragments and optimized knowledge-based force field. Proteins. 2012;80:1715–35.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Berendsen HJC, van der Spoel D, van Drunen R. GROMACS: a message-passing parallel molecular dynamics implementation. Comp Phys Comm. 1995;91:43–56.

    CAS  Article  Google Scholar 

  29. 29.

    Hess B, Kutzner C, van der Spoel D, Lindahl E. GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J Chem Theory Comput. 2008;4:435–47.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Lindahl E, Hess B, van der Spoel D. GROMACS 3.0: a package for molecular simulation and trajectory analysis. J Mol Model. 2001;7(8):306–17.

    CAS  Google Scholar 

  31. 31.

    Pandini A, Fornili A, Fraternali F, Kleinjung J. GSATools: analysis of allosteric communication and functional local motions using a structural alphabet. Bioinformatics. 2013;29:2053–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Van Der Spoel D, Lindahl E, Hess B, Groenhof G, Mark AE, Berendsen HJC. GROMACS: fast, flexible, and free. J Comput Chem. 2005;26:1701–18.

    Article  Google Scholar 

  33. 33.

    Lindorff-Larsen K, Piana S, Palmo K, Maragakis P, Klepeis JL, Dror RO, Shaw DE. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins. 2010;78:1950–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Wiederstein M, Sippl MJ. ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 2007;35:407–10.

    Article  Google Scholar 

  35. 35.

    Sippl MJ. Recognition of errors in three-dimensional structures of proteins. Proteins. 1993;17:355–62.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Benkert P, Biasini M, Schwede T. Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics. 2011;27:343–50.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Benkert P, Tosatto SCE, Schomburg D. QMEAN: a comprehensive scoring function for model quality assessment. Proteins. 2008;71:261–77.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Benkert P, Künzli M, Schwede T. QMEAN server for protein model quality estimation. Nucleic Acids Res. 2009;37:510–4.

    Article  Google Scholar 

  39. 39.

    Benkert P, Schwede T, Tosatto SC. QMEANclust: estimation of protein model quality by combining a composite scoring function with structural density information. BMC Struct Biol. 2009;9:35.

    Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Barthel D, Hirst JD, Błazewicz J, Burke EK, Krasnogor N. ProCKSI: a decision support system for protein (structure) comparison, knowledge, similarity and information. BMC Bioinformatics. 2007;8:416.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Zhang Y, Skolnick J. TM-align: a protein structure alignment algorithm based on the TM-score. Nucleic Acids Res. 2005;33:2302–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Frickey T, Lupas A. CLANS: a Java application for visualizing protein families based on pairwise similarity. Bioinformatics. 2004;20:3702–4.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Remmert M, Biegert A, Hauser A, Söding J. HHblits: lightning-fast iterative protein sequence searching by HMM-HMM alignment. Nat Methods. 2012;9:173–5.

    CAS  Article  Google Scholar 

  44. 44.

    Söding J, Biegert A, Lupas AN. The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res. 2005;33:244–8.

    Article  Google Scholar 

  45. 45.

    Söding J. Protein homology detection by HMM-HMM comparison. Bioinformatics. 2005;21:951–60.

    Article  PubMed  Google Scholar 

  46. 46.

    Lensink MF, Wodak SJ. Docking, scoring, and affinity prediction in CAPRI. Proteins. 2013;81:2082–95.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Comeau SR, Gatchell DW, Vajda S, Camacho CJ. ClusPro: an automated docking and discrimination method for the prediction of protein complexes. Bioinformatics. 2004;20:45–50.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Comeau SR, Gatchell DW, Vajda S, Camacho CJ. ClusPro: a fully automated algorithm for protein-protein docking. Nucleic Acids Res. 2004;32:96–9.

    Article  Google Scholar 

  49. 49.

    Kozakov D, Brenke R, Comeau SR, Vajda S. PIPER: an FFT-based protein docking program with pairwise potentials. Proteins. 2006;65:392–406.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Kozakov D, Beglov D, Bohnuud T, Mottarella SE, Xia B, Hall DR, Vajda S. How good is automated protein docking? Proteins. 2013;81:2159–66.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Torchala M, Moal IH, Chaleil RAG, Agius R, Bates PA. A Markov-chain model description of binding funnels to enhance the ranking of docked solutions. Proteins. 2013;81:2143–9.

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Torchala M, Bates PA. Predicting the structure of protein-protein complexes using the SwarmDock Web Server. Methods Mol Biol. 2014;1137:181–97.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Torchala M, Moal IH, Chaleil RAG, Fernandez-Recio J, Bates PA. SwarmDock: a server for flexible protein-protein docking. Bioinformatics. 2013;29:807–9.

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    Viswanath C, Ravikant DVS, Elber R. DOCK/PIERR: web server for structure prediction of protein-protein complexes. Methods Mol Biol. 2014;1137:199–207.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Tovchigrechko A, Vakser IA. GRAMM-X public web server for protein–protein docking. Nucleic Acids Res. 2006;34(2):W310–4.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  56. 56.

    Volkamer A, Kuhn D, Grombacher T, Rippmann F, Rarey M. Combining global and local measures for structure-based druggability predictions. J Chem Inf Model. 2012;52:360–72.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Brenke R, Kozakov D, Chuang G-Y, Beglov D, Hall D, Landon MR, Mattos C, Vajda S. Fragment-based identification of druggable “hot spots” of proteins using Fourier domain correlation techniques. Bioinformatics. 2009;25:621–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  58. 58.

    Kozakov D, Hall DR, Chuang G-Y, Cencic R, Brenke R, Grove LE, Beglov D, Pelletier J, Whitty A, Vajda S. Structural conservation of druggable hot spots in protein-protein interfaces. Proc Natl Acad Sci USA. 2011;108:13528–33.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Kozakov D, Grove LE, Hall DR, Bohnuud T, Mottarella SE, Luo L, Xia B, Beglov D, Vajda S. The FTMap family of web servers for determining and characterizing ligand-binding hot spots of proteins. Nat Protoc. 2015;10:733–55.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  60. 60.

    Ashkenazy H, Erez E, Martz E, Pupko T, Ben-Tal N. ConSurf 2010: calculating evolutionary conservation in sequence and structure of proteins and nucleic acids. Nucleic Acids Res. 2010;38:529–33.

    Article  Google Scholar 

  61. 61.

    Celniker G, Nimrod G, Ashkenazy H, Glaser F, Martz E, Mayrose I, Pupko T, Ben-Tal N. ConSurf: using evolutionary data to raise testable hypotheses about protein function. Isr J Chem. 2013;53:199–206.

    CAS  Article  Google Scholar 

  62. 62.

    Glaser F, Pupko T, Paz I, Bell RE, Bechor-Shental D, Martz E, Ben-Tal N. ConSurf: identification of functional regions in proteins by surface-mapping of phylogenetic information. Bioinformatics. 2003;19:163–4.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Landau M, Mayrose I, Rosenberg Y, Glaser F, Martz E, Pupko T, Ben-Tal N. ConSurf 2005: the projection of evolutionary conservation scores of residues on protein structures. Nucleic Acids Res. 2005;33:299–302.

    Article  Google Scholar 

  64. 64.

    Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004;5:113.

    Article  Google Scholar 

  65. 65.

    Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  66. 66.

    Waterhouse AM, Procter JB, Martin DMA, Clamp M, Barton GJ. Jalview version 2—a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009;25:1189–91.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  67. 67.

    Ferrè F, Clote P. DiANNA: a web server for disulfide connectivity prediction. Nucleic Acids Res. 2005;33:230–2.

    Article  Google Scholar 

  68. 68.

    Krissinel E, Henrick K. Secondary-structure matching (SSM), a new tool for fast protein structure alignment in three dimensions. Acta Crystallogr D Biol Crystallogr. 2004;60(12):2256–68.

    CAS  Article  PubMed  Google Scholar 

  69. 69.

    Holm L, Rosenström P. Dali server: conservation mapping in 3D. Nucleic Acids Res. 2010;38:545–9.

    Article  Google Scholar 

  70. 70.

    Bisson-Filho AW, Discola KF, Castellen P, Blasios V, Martins A, Sforça ML, Garcia W, Zeri ACM, Erickson HP, Dessen A, Gueiros-Filho FJ. FtsZ filament capping by MciZ, a developmental regulator of bacterial division. Proc Natl Acad Sci USA. 2015;112:E2130–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  71. 71.

    Schaffner-Barbero C, Martín-Fontecha M, Chacón P, Andreu JM. Targeting the assembly of bacterial cell division protein FtsZ with small molecules. ACS Chem Biol. 2012;7:269–77.

    CAS  Article  PubMed  Google Scholar 

Download references

Authors’ contributions

AJS carried out the experiments. DJR conceived and supervised the work. Both authors read and approved the final manuscript.


Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and material

Modelled structures are available on request from the authors.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.


None to report.

Author information



Corresponding author

Correspondence to Daniel J. Rigden.

Additional file

Additional file 1: Figure S1.

The alignment of FtsZ derived from P. aeruginosa, B. subtilis, S. aureus, and M. tuberculosis using MUltiple Sequence Comparison by Log-Expectation (muscle) [64, 65], and coloured using the clustal colour scheme. The alignment illustrates that E. coli FtsZ follows a known domain structure of FtsZ with a variable N-terminal segment (1–10), a highly conserved core region (10-316), a variable spacer (316–370) and finally a conserved C-terminal peptide (370–379) [13]. The figure was made with Jalview [66].

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Simpkin, A.J., Rigden, D.J. GP0.4 from bacteriophage T7: in silico characterisation of its structure and interaction with E. coli FtsZ. BMC Res Notes 9, 343 (2016).

Download citation


  • ab initio modelling
  • Homology modelling
  • Docking
  • FtsZ
  • Drug design