- Research article
- Open Access
Type I pyridoxal 5′-phosphate dependent enzymatic domains embedded within multimodular nonribosomal peptide synthetase and polyketide synthase assembly lines
© Milano et al.; licensee BioMed Central Ltd. 2013
- Received: 24 June 2013
- Accepted: 14 October 2013
- Published: 23 October 2013
Pyridoxal 5′-phosphate (PLP)-dependent enzymes of fold type I, the most studied structural class of the PLP-dependent enzyme superfamily, are known to exist as stand-alone homodimers or homotetramers. These enzymes have been found also embedded in multimodular and multidomain assembly lines involved in the biosynthesis of polyketides (PKS) and nonribosomal peptides (NRPS). The aim of this work is to provide a proteome-wide view of the distribution and characteristics of type I domains covalently integrated in these assemblies in prokaryotes.
An ad-hoc Hidden Markov profile was calculated using a sequence alignment derived from a multiple structural superposition of distantly related PLP-enzymes of fold type I. The profile was utilized to scan the sequence databank and to collect the proteins containing at least one type I domain linked to a component of an assembly line in bacterial genomes. The domains adjacent to a carrier protein were further investigated. Phylogenetic analysis suggested the presence of four PLP-dependent families: Aminotran_3, Beta_elim_lyase and Pyridoxal_deC, occurring mainly within mixed NRPS/PKS clusters, and Aminotran_1_2 found mainly in PKS clusters. Sequence similarity to the reference PLP enzymes with solved structures ranged from 24 to 42% identity. Homology models were built for each representative type I domain and molecular docking simulations with putative substrates were carried out. Prediction of the protein-protein interaction sites evidenced that the surface regions of the type I domains embedded within multienzyme assemblies were different from those of the self-standing enzymes; these structural features appear to be required for productive interactions with the adjacent domains in a multidomain context.
This work provides a systematic view of the occurrence of type I domain within NRPS and PKS assembly lines and it predicts their structural characteristics using computational methods. Comparison with the corresponding stand-alone enzymes highlighted the common and different traits related to various aspects of their structure-function relationship. Therefore, the results of this work, on one hand contribute to the understanding of the functional and structural diversity of the PLP-dependent type I enzymes and, on the other, pave the way to further studies aimed at their applications in combinatorial biosynthesis.
- Pyridoxal 5′-phosphate
- Fold type I
- Nonribosomal peptide synthetases
- Polyketide synthases
- Tailoring domains
- Hidden Markov models
- Homology modeling
- Protein-protein interaction
Pyridoxal 5′-phosphate (PLP), a derivative of Vitamin B6, is one of the most versatile organic cofactors in biology. In fact, PLP-dependent enzymes form a vast and complex group of proteins present in organisms belonging to all levels of the tree of life  and participate in a variety of reactions (Additional file 1: Scheme 1). In humans, for example, besides the classical role in transamination, many of these enzymes take part in the metabolism of neurotransmitters such as dopamine, serotonin, glycine, epinephrine, norepinephrine, d-serine, l-glutamate, γ-amino butyric acid and histamine . PLP-dependent enzymes have been classified in at least five evolutionarily unrelated families, characterized by specific three-dimensional folds [1, 3]. Many studies have been devoted to the elaboration of a rigorous classification of PLP enzymes, with the aim to identify their common structural features and to understand how different protein scaffolds can support similar substrate binding in the active sites . Among the different structural classes, the so-called fold type I  is the most populated in nature and the best characterized one. The subunit architecture of the fold includes one large and one small domain. The large domain contains a seven stranded β-sheet interacting with α-helices. The small domain at the C-terminal part of the chain folds as a three- or four-stranded β-sheet partly covered with helices. PLP-enzymes belonging to the fold type I, are known to exist as stand-alone proteins either homodimers or homotetramers. The active site is located in a crevice between the two domains at the subunit interface. The archetypal protein of this class is aspartate aminotransferase  which was the first PLP-dependent enzyme to be purified and crystallized. Since then, intense research work has been carried out to elucidate the details of its structural and functional properties.
Scrutiny of sequence data produced by genomic projects showed that fold type I domains can be found in multidomain frameworks, in prokaryotic systems, such as the transcriptional regulator MocR [7, 8] and also in multienzyme systems, polyketide synthases (PKS) and nonribosomal peptide synthetases (NRPS), involved in the biosynthesis of polyketides (PK), nonribosomal peptides (NRP) and hybrid PK/NRP secondary metabolites. A vast array of bioactive metabolites belonging to these classes includes also several important medicinal agents and biotechnologically relevant compounds [9, 10]. The canonical biosynthetic mechanisms of the two classes of structurally different secondary metabolites, polyketides (PK) and nonribosomal peptides (NRP), share several common features. Both NRPS and PKS type I systems require the participation of multienzyme complexes acting as assembly lines for the construction of polyketide or peptide chains by a sequence of condensation steps. In both systems each elongation step is catalyzed by a module containing the catalytic domains required for the insertion of a monomer into the growing chain. The first step in the biosynthetic pathway, is the ATP-dependent adenylation of the amino acid catalyzed by an adenylation domain (A) in the NRPS systems and the transfer of the acyl-group from an acyl-CoA onto the acyltransferase unit (AT) in the PKS systems. The monomer is then transferred to a carrier protein (CP) post-translationally primed with a phosphopantetheine arm, called thiolation domain (T) or peptidyl carrier protein (PCP) in NRPS and acyl carrier protein (ACP) in PKS multienzyme assemblies, respectively. In both systems, the carrier proteins mediate the transport of intermediates, linked by a thioester bond to the phosphopantetheine arm, along the assembly line. In NRPS systems, the key elongation reaction is the peptide bond formation by a nucleophilic attack of the α-amino group of the amino acid tethered to the downstream thiolation domain on the thioester bond of the intermediate tethered to the upstream peptidyl carrier protein, catalyzed by a condensation domain (C). In PKS assembly lines the elongation relies on the ketosynthase (KS) catalyzed carbon-carbon bond formation, by a Claisen condensation mechanism between the upstream acyl thioester and the downstream carbanionic acyl acceptor resulting from decarboxylation of malonyl- or methylmalonyl-ACP. The release of the product in both systems is usually catalyzed by a thioesterase (TE) domain located in the termination module and in most cases involves a macrocyclization; however different mechanisms have also been reported [9, 10]. In addition to these two distinct biosynthetic mechanisms, there is a large number of mixed clusters involved in the production of structurally complex compounds where both polyketide and peptide moieties can be recognized. The modular architecture and functional versatility makes possible the switching between NRPS and PKS assembly lines.
In NRPS, PKS and hybrid systems, besides the essential aforementioned catalytic domains, the presence of additional “tailoring domains”, which introduce structural modifications into the canonical building blocks and contribute to the amazing structural diversity in NRP and PK metabolites, is frequently observed. These domains can be encoded within the biosynthetic gene clusters, either fused to other catalytic domains or as self-standing domains . Their activities range from hydroxylation, halogenation, methylation, racemization, heterocyclization, lipidation and glycosilation. Several PLP-dependent proteins were also identified in multimodular biosynthetic machineries, some operating as self-standing domains, others incorporated in multidomain enzymes containing at least one carrier protein. Examples of PLP-dependent enzymes postulated to be involved in the formation of building blocks are found in the biosynthesis of peptidyl nucleoside antibiotic pacidamycin: PacE and PacS belonging to the type II fold, and PacT belonging to the type I fold . A PLP-dependent protein is involved in an interesting and very unusual chain releasing mechanism in the biosynthesis of the fungal polyketide mycotoxin fumonisin, namely by incorporation of two carbons and one amino group from alanine into the acyl chain . A stand-alone PLP dependent enzyme, MxcL, is hypothesized to participate in the final step of the biosynthesis of myxochelin B, the catecholate siderophore produced by Stigmatella aurantiaca Sg a15 , namely in the transamination of the aldehyde group present in the late biosynthetic intermediate. The insertion of an amino group into a polyketide biosynthetic precursor by transamination of the carbonyl function is a process operating also in the biosynthesis of antimicrobial polyamino antibiotics zeamines produced by Serratia plymuthica RVH1 . In particular, the gene zmn12 present in the complex biosynthetic cluster of this compound, encodes a protein containing a domain with a putative aminotransferase activity, homologous to the type I PLP–dependent glutamate-1-semialdehyde aminotransferase.
The functionally characterized PLP-dependent domains belonging to fold type I which operate in cis within mixed NRPS/PKS multienzyme systems are those involved in the biosynthesis of the potent antifungal cyclic lipopeptide mycosubtilin  and of the tripyrrolic metabolite prodigiosin . In the biosynthesis of mycosubilin , an aminotransferase domain (AMT) embedded within the PKS/NRPS hybrid enzyme MycA and located at the interface of the PKS and NRPS modules, catalyzes the incorporation of an amine group from the amine donor, Gln, into the protein-bound PLP and subsequently to the β-ketothioester tethered to the ACP domain of the polyketide moiety. A different role of a PLP domain was established in the formation of prodigiosin belonging to the family of tripyrrole red pigments prodiginines produced by Serratia and Streptomyces bacterial strains, which are attracting increasing interest because of their immunosuppressive, anticancer, antimicrobial, and antimalarial activities. The PLP-dependent domain, SerT, located on a module containing also two ACP domains, PigH, is predicted to generate a C2 fragment by decarboxylation of l-serine, which is then used for pyrrole B ring formation .
In this paper, we focus specifically on PLP-dependent domains of fold type I occurring covalently linked in multidomain frameworks related to NRPS and/or PKS-like assemblies in bacterial systems. Since the identification of these domains is relatively recent, we undertook an in silico analysis with the aim to contribute to the clarification of some aspects concerning their function, structural remodeling and relationship to the homologous, traditional PLP-dependent enzymes.
Construction of the hidden Markov model representative of fold type I PLP enzymes
List of structures of PLP-dependent type-I domains utilized for the calculation of the PLP_domain profile
PDB code a)
L-cysteine/L-cystine C-S lyase
Thermus thermophilus HB8
Cystathionine gamma lyase
Glutamate decarboxylase beta
Detection of sequences containing the type I PLP-domains through databank searches
At the completion of the databank searches, 206 sequences were collected using the criteria and the filtering procedure reported in Methods section; a detailed list is reported in Additional file 1: Table S1.
Phylogenic analysis of type I PLP domains embedded in the NRPS or PKS multienzyme assemblies
List of non redundant NRPS/PKS proteins containing type I PLP-dependent domains of the four families
Domain layout b)
Cluster type d)
Rhodobacter sp. SW2
Streptomyces sp. e14
Salinisphaera shabanensis E1L3A
Thiocapsa marina 5811
Thiorhodococcus drewsii AZ1
Streptomyces albus J1074
Lysobacter sp. ATCC 53042
The results of the annotation of the type I domain sequences through the B6 database  show heterogeneity with respect to the Pfam classification: the Aminotran_3 family indeed contains glutamate-1-semialdehyde 2,1-aminomutase, taurine--pyruvate aminotransferase, diamine aminotransferase, diaminobutyrate-2-oxoglutarate transaminase, ornithine--oxo-acid aminotransferase (Table 2). The Pfam families Aminotran_1_2 and Beta_elim_lyase appear homogeneous since they contain only 8-amino-7-oxononanoate synthase and tyrosine phenol lyase, respectively. The last and the less populated family Pyridoxal_deC contains, as expected, decarboxylases namely diaminobutyrate decarboxylase and glutamate decarboxylase.
Description of the domain architecture and sequence analysis of multidomain assemblies containing type I PLP enzymes
The organization and the identity of the domains contained in the parent sequences from which the type I domains were extracted, was determined through the script “Pfam_scan.pl” (see Methods). These results and those deriving from the antiSMASH [21, 22] analysis are summarized in Table 2 for the non-redundant subset and in Additional file 1: Table S1 for the entire set, respectively. Apparently, Aminotran_3 family represents the vast majority of the type I domains collected. They occur almost invariantly in mixed PKS/NRPS assemblies.
Regarding the family denoted by the Pfam tag Aminotran_1_2, it can be noted that many assemblies do not contain more than three domains. However, they are very likely involved in pathways related to NRPS or PKS because the corresponding coding sequences are often adjacent to those characteristic of such biosynthetic clusters (results not shown). For example, the coding sequence [UniProt: A0SZ00] (Additional file 1: Table S2) from Janthinobacterium lividum [GenBank:ABK64042], is located between the sequences [GenBank:ABK64060] and [GenBank:ABK64039], corresponding to a putative peptidyl carrier protein and a putative L-prolyl-AMP-ligase, respectively.
The Beta_elim_lyase domains occur in predicted mixed PKS/NRPS transacting clusters. Interestingly, the type I domains are incorporated in modules missing any A or AT domain. This situation is reminiscent of the cluster involved in the biosynthesis of the mixed NRP/PK metabolite leinamycin from Streptomyces atroolivaceus S-140  where the gene lnmJ encodes six PKS modules lacking the AT domains and a domain homologous to tyrosine phenol-lyases ([UniProt: Q8GGP2] in Table 2). It was also experimentally proved that the missing activities were provided by a discrete AT enzyme that loads the extender units in trans [23, 24]. Moreover, we found a similar example in the mixed NRPS/PKS gene cluster 6 predicted by antiSMASH  analysis of the genome of the bacteria Catenulispora acidiphila. In this cluster, to which belongs the sequence [UniProt:C7PXR3] shown in Table 2, the occurrence of a stand-alone A domain ([UniProt:C7PXP4]) is predicted.
Pyridoxal_deC domains are the rarest since they occur only in two instances of our set, one of which could not be annotated by antiSMASH .
Molecular modeling of the type I domains and docking of putative substrates
In the Aminotran_1_2 group, residues sandwiching the PLP cofactor, namely His137 and Ala211 (numbering refers to the template structure, Figure 3b), are invariant in all the homologous sequences considered (Figure 5b and Additional file 1: Figure S2b). Thr111 and Ser272 from the other subunit interact with the cofactor phosphate group and are conserved. His212 in Figure 5b is involved in the interaction with the O3 atom of PLP and is identically conserved in the homologous sequences. Regarding the binding of the substrate, it can be observed that His137 is involved in the binding of the substrate carboxylate . The same role is predicted in the docked complex between substrate Ser and the homology model (Additional file 1: Figure S3b). Likewise, Arg366 should also be mentioned among the residues involved in substrate binding in the serine palmitoyltransferase template. This residue is identically conserved and it is essential for the catalysis because it forms the key PLP:l-serine quinonoid intermediate that condenses with palmitoyl-CoA . Finally, according to the docked model of the PLP:l-serine quinonoid intermediate, the β-hydroxyl moiety of the substrate participates to a hydrogen bond with Ser242, residue conserved in all the selected members of the Aminotran_1_2 family (Additional file 1: Figure S2b).
Active site model of the type I domain representative of the Beta_elim_lyase family shows that some residues relevant for catalysis in the representative member of the group tyrosine phenol-lyase  are conserved. For example, Phe105 and Thr198, the residues sandwiching the PLP cofactor, are conserved (Figure 5c and Additional file 1: Figure S2c) although the latter residue is replaced by Ser in one case. The residues Arg386, Arg199 and Asn167, deemed to be involved in interaction with the carboxylic group of the reaction intermediate in tyrosine phenol-lyase, are conserved (Additional file 1: Figure S3c). On the contrary, residues that in the postulated mechanism of β-elimination reaction carry out the protonation of the substrate Cγ, are not present in the Beta_elim_lyase of type I embedded in multidomain context. In particular, Tyr53 of tyrosine phenol-lyase is replaced by an Arg residue in the model domain, while Arg381 that assists Tyr53 during the protonation is deleted in all the sequences reported in Additional file 1: Figure S2c, including the PLP domain involved in the biosynthesis of leinamycin [UniProt:Q8GGP2].
The Pyridoxal_deC group is the least populated; assessment of evolutionary conservation of residues possibly involved in catalysis is more difficult. His111 and Ala193 that, as in the Aminotran_1_2 family sandwich the PLP ring, are conserved. His222 and Ser71, ligands of the cofactor phosphate group, are also conserved. According to the docked model (Additional file 1: Figure S3d), the α-carboxylate group of the substrate is stabilized by an ion-pair interaction with Arg384. This residue is equivalent to Arg567 in glutamic acid decarboxylase, enzyme of the same family, that is responsible for the formation of a salt bridge to its substrate γ-aminobutyric acid . Similarly Tyr253, the catalytic residue of the Group II decarboxylases that performs the protonation of the Cα atom of the quinonoid intermediate , is conserved.
Prediction of protein-protein interaction sites on the homology models of the type I domains
Prediction of the presence of potential protein-protein binding sites was carried out for the models and the relative template structures in their dimeric forms. The results suggest that some type I domains covalently incorporated in multidomain contexts possess potential protein binding sites missing in the equivalent regions of their respective structural templates. These regions are in proximity of the active sites and contain residues largely conserved in the corresponding homologs. In particular, the Aminotran_3 domain displays a small potential region encompassing the high interaction-probability residues Phe25, Ile29, Lys30, Met32, Asp380 and Gly396 (positions are relative to the template numbering system in Figure 3a). The first four residues are located in the poorly conserved N-terminal region (Additional file 1: Figure S2a). Along with the other potentially interacting sites, they form a region around the active site mouth. The Aminotran_1_2 domain shows a surface region predicted as potential protein-protein interaction site that includes Asp333 and the sequence Ala353-Lys362 (the positions refer to the template numbering system in Figure 3b). This region is located on the rim of the active site. However, the potentially interacting surface of the model is, in this case, similar to that observed in the structural template (Additional file 1: Figure S4b). The putative interaction between the CP and the type I domains of polyketide synthase from Burkholderia thailandensis [UniProt:Q2T5Z2] have been predicted by protein-protein docking simulation using the ClusPro server . The homology model of the CP domain encompassed by the positions 932–998 of [UniProt: Q2T5Z2] has been calculated using the templates denoted by [PDB:2EHS], [PDB:1X3O], [PDB:2QNW] and [PDB:2JU2]. Although the results deriving from the docking experiments carried out with homology models should be considered with the great caution, it is interesting to note that the ten best complexes calculated by the ClusPro analysis suggest that the CP domain may interact with the surface regions of the PLP domain predicted as potential protein-protein interaction sites (Additional file 1: Figure S5). The Beta_elim_lyase model possesses a wide surface predicted as a potential interaction site, significantly larger than that predicted in the corresponding structural template (Additional file 1: Figure S4c). The interacting surface clusters into two patches located in the proximity of the active site. The first cluster is centered on the residues Arg343, His420, Gly423, Gly424, Pro431 and Tyr432. The second cluster incorporates the residues Leu104, Phe105, Pro106, Ile109 and Tyr110. The Pyridoxal_deC family, represented by the model of type I domain in the nonribosomal peptide synthetase module from Coxiella burnetii [UniProt:B6IZA3], displays a surface interaction propensity similar to the reference template structure of human cysteine sulfinic acid decarboxylase [PDB:2JIS] (Figure 3d). However, according to multiple sequence alignments between members of this type I domain and the other known members of the Group II decarboxylases family, the characteristic third N-terminal domain (N-domain) formed by three α-helices that fold upon dimerization , is lost in Pyridoxal_deC family upon incorporation in a multidomain context. It is therefore tempting to speculate that this structural deletion, which would result in an alternate entry into the active site, has evolved to accommodate the adjacent domains of the NRPS framework, e.g., a phosphopantetheine binding domain, and that α-decarboxylation could occur on the amino acid substrate tethered to the phosphopantetheine arm.
In this work we focused specifically on the type I PLP domains involved in putative bacterial NRPS and/or PKS assemblies as tailoring domains operating in cis. We collected a set of predicted type I PLP sequences incorporated in multidomain environments from eubacterial sources using the HMM search using the PLP_domain profile. Only those sequences containing a phosphopantetheine binding domain (CP) were considered. In most cases, the type I domains are placed downstream from a CP domain (Additional file 1: Table S1). As a test, the selection of sequences was carried out using as a criterion the association of the type I domains with adenylation, condensation or acyltransferase domains. In all cases, only subsets of the group of sequences retrieved with the original requirement were recovered, showing the correctness of the adopted procedure. The only exception to this “rule” was the β-ketoacyl synthase from Streptomyces violaceusniger [UniProt:G2P368]. This 1136 residue long sequence is atypical since it does not possess a recognizable phosphopantetheine binding domain and displays a segment, containing about 400 residues, in the central sequence region apparently lacking any relation to known Pfam domains. In a few cases [UniProt:G8X2R2] or [UniProt:Q82RP2], Pfam annotation does not report the presence of a CP domain. However, we were able to detect those domains using the “Pfam_scan.pl” annotation script. In some cases, use of a locally installed program provides an effective way of finely tuning the parameters.
The assemblies containing the type I domains we collected, were found only in the following eubacterial phyla: Acidobacteria, Actinobacteria, Bacteroidetes, Chloroflexi, Cyanobacteria, Firmicutes, Planctomycetes, Proteobacteria, and Verrucomicrobia. However, it should be considered that the observed distribution may be significantly biased by uneven species sampling.
Sequence alignments and phylogenetic analyses indicated that four groups can be distinguished within the collected set of the type I domains. The four groups share the identical conservation (with the few exceptions reported in Results section) of the key residues Asp and Lys which are involved in interaction with the cofactor pyridine nitrogen and in the formation of a Schiff base with the aldehyde group of PLP, respectively.
We further studied the structural features characterizing these type I protein subfamilies by multiple sequence alignment, homology modeling and search of the potential binding sites possibly involved in the interaction with protein partners.
The most numerous group is the Aminotran_3 family which is structurally related to the glutamate 1-semialdehyde-2,1-aminomutase . The proteins belonging to this group are predicted to occur mainly in mixed NRPS/PKS machineries. The aminotransferase domain embedded in the multidomain PKS/NRPS enzyme MycA from Bacillus subtilis involved in the synthesis of the cyclic lipopeptide mycosubtilin , is the example of the functionally characterized member of the group. The role of the aminotransferase domain is the insertion of the amino group into the polyketide biosynthetic intermediate; the amine source for the MycA enzyme was demonstrated to be glutamine . A similar function was proposed also for the aminotransferase domains encoded by the genes mxcL and zea12 [13, 14, 33] from myxochelin and zeamine biosynthetic clusters respectively, which belong to the same subfamily.
The multiple sequence alignment between a non-redundant set of sequences and the reference modeling template evidenced the structural features characteristic of these domains. As shown in Figure 3a, there is an insertion of about 12–14 residues at position 288–289 (numbering system of the structural template) that corresponds to the region 283–284 in Additional file 1: Figure S2a. The insertion is located on the surface of the domain (Figure 4a) and may be involved in modulating the interaction with the adjacent domains and/or in substrate recognition. Indeed, docking simulations suggests that a loop is involved in interaction with the CP domain (Additional file 1: Figure S5). In the surface area between the sequence positions 364–378 of the alignment in Figure 3a (corresponding to the region 364–378 in Additional file 1: Figure S2a), a deletion region is present. Further, the prediction of potential sites of protein-protein interaction assigns a significant potential to a region in the active site proximity, not visible in the equivalent position of the structural template (Additional file 1: Figure S4a). These observations, in particular the presence of a wide region of potential protein-protein interaction (Additional file 1: Figure S4b), suggest that this structural feature could be required for productive interaction with the adjacent domains in a multidomain context.
The Aminotran_1_2 group is structurally related to the Coenzyme A (CoA) family of the PLP-dependent enzymes of type I . They are predicted to be involved mainly in PKS pathways. Within the CoA subfamily, serine palmitoyltransferase is the most closely related enzyme. Sequence alignments and homology modeling indicate the conservation of the residues involved in substrate-CoA interaction and the absence of long insertion/deletions with respect to the structural template (Figure 3b). Exception is the sequence of the type I domain from the Anabaena circinalis protein [UniProt:B3EYK4] that shows two insertions, one of which 10-residue long (Additional file 1: Figure S2b). The only experimentally characterized member of this family is the PigH protein [UniProt: Q5W247] from Serratia marcescens , involved in the biosynthesis of the three-pyrrolic red pigment prodigiosin. PigH contains two CP domains followed by a type I PLP domain, SerT (as reported in Additional file 1: Table S1), predicted to catalyze the decarboxylation of l-serine and the formation of C2 fragment used in the formation of the pyrrole B ring of prodigiosin. Analysis of surface of the homology model of the representative type I domain of this family suggests the presence of an increased potential for protein interaction in the proximity of the active site mouth (Additional file 1: Figure S4b). However, the increase of the interaction potential with respect to the stand-alone counterpart is less evident than in the case of the Aminotran-3 family.
Extensive remodeling of the protein surface can be conjectured also for the Beta_elim_lyase family. In fact, as in Aminotran_3 group a wide region of potential protein-protein interaction (Additional file 1: Figure S4b) is observed. The sequence divergence from the model is evident at the N-terminal part of the domain. This family is indeed characterized by domains occurring mainly at the C-terminal edge of the multidomain module. Despite the conservation of many active site key residues observed in the tyrosine phenol-lyase enzyme, some side chains involved in the catalytic mechanisms are missing. In particular, as mentioned in the Results section, Tyr53 of tyrosine phenol-lyase is replaced by an Arg residue while Arg381 is missing. Although the relevance of such variations on the functionality of these domains cannot be presently assessed, it is worth mentioning that studies on classical β-eliminating lyases showed how the mutation of the residues corresponding to Arg381 and Tyr53 affected the activity of two proteins belonging to this family. The substitution of Arg381 with Ile and Val in tyrosine phenol-lyase caused a significant impairment in the activity towards Tyr, but not in case of other substrates; moreover the same substitutions were present in the wild type and fully active tryptophan indole-lyase enzymes . Similarly, the replacement of Tyr71 (corresponding to Tyr53) had different effects on the activity of the protein towards different substrates . These findings indicate that Tyr53 and Arg381 are not absolutely indispensable for the catalytic activity of all the PLP enzymes of this family. Regarding the prediction of protein interaction sites on the type I domains, the results suggest the presence of a wide area possibly involved in the interaction with other protein partners as observed in other aforementioned families.
The pyridoxal_deC family, the least numerous, is related to the decarboxylase family II  to which the enzyme glutamate decarboxylase belongs. However, the most similar reference structure found was cysteine sulfinic acid decarboxylase, an enzyme involved in hypotaurine biosynthesis  which functions as an autoantigen in human endocrine autoimmune diseases . Conservation of the residues essential for catalysis in glutamate decarboxylase, which is better characterized than cysteine sulfinic acid decarboxylase, suggests the possible preservation of the decarboxylase activity in this domain. The databank search showed that the pathogenic proteobacterium Coxiella burnetii, whose genome has been completely sequenced, possesses one of the two Pyridoxal_deC domains found during our databank searches. The other, showing high similarity with glutamate decarboxylase domains, was found in the cyanobacterium Lyngbya majuscula, in particular in the biosynthetic cluster of jamaicamide. Interestingly, it is embedded within the adenylation domain of the PKS/NRPS multidomain subunit JamL, but the precise function of the PLP-dependent domain in the biosynthesis of this metabolite has not been so far clarified .
The results of our work show that the domains belonging to the type I PLP dependent enzymes linked to a component of the multidomain frameworks related to NRPS and/or PKS-like assemblies are relatively rare but widespread among several bacterial phyla. These domains display conservation (except in a few cases in the Aminotran_3 family) of residues involved in cofactor binding and catalysis. However the prediction of protein-protein interaction sites suggests that the N- and C-terminal ends of the domain polypeptide chain display stronger sequence divergence with respect to the reference stand-alone structures (Additional file 1: Figure S2a). These regions are necessarily involved in connecting the linkers bridging the other domains in the multidomain subunits (Table 1). On the other hand, it must stressed that predictions of interacting sites are still very inaccurate and in this case they can be significantly biased by modeling inaccuracies especially related to the prediction of side chain solvent accessibility. Nevertheless, the differences between the prediction in the surface regions possibly involved in protein-protein interaction of the model and the templates are, at least in the case of Aminotran_3 and Aminotran_1_2 groups, very strong. These two families display the highest sequence similarity between the model and the template. Therefore, the strong differences observed can represent significant signals, while differences observed in the models of other families are less reliable. This hypothesis is supported by the results of the docking experiment carried out using the homology model of the PLP type I and CP domains of the Aminotran_3 family. Indeed, although the results of docking experiments carried out with homology models should be taken with the great caution, the ten best complexes indicate that the CP domain may interact with the surface regions of the PLP domain predicted as potential protein-protein interaction sites (Additional file 1: Figure S5).
In this context the quaternary architecture of multidomain assemblies incorporating a PLP-dependent type I domain should also be considered. Indeed, the PLP type I domains are dimers or tetramers because the proper formation of their catalytically competent active site requires the participation of residues from the adjacent subunits . This structural constraint fits well with the dimeric architecture of PKS systems. On the other hand, studies with individual NRPS domains showed that they were monomers ; however a dimeric structure was demonstrated in the multidomain synthetase VibF and a continuum of monomeric and dimeric oligomerization states in NRPS was proposed . The existence of a number of secondary metabolites of mixed PKS/NRPS origin shows that the two biosynthetic machineries are compatible and studies on multienzyme docking in hybrid megasynthetases indicated that NRPS subunits in mixed systems self-associate to interact with partner PKS homodimers .
This work offers a systematic view of the occurrence of the type I PLP-dependent enzymes within NRPS and PKS assembly lines and predicts their structural characteristics using in silico methods. The results of this research contribute to a deeper understanding of the functional and structural diversity of the PLP-enzyme family of fold type I and pave the way to further studies aimed at their applications in combinatorial biosynthesis. In fact, the success in the functioning of engineered biosynthetic clusters depends, to a great extent, on efficient molecular recognition between the single components.
Data sources and computational tools
All sequence data processed during the work were taken from the UniProt  release April, 2012 or Protein Data Banks . Most of the databank searches and analyses utilized the profile Hidden Markov Model (HMM) methodology as implemented in the package HMMER 3.0  or relied on the BLAST suite . Multiple sequence alignments were calculated either with the programs Clustal-W , MAFFT  or hmmalign ; sequence editing and alignment display relied on Jalview  or Seaview  editors.
Hidden Markov model of type I PLP-enzymes
A Hidden Markov Model (HMM) of the superfamily of PLP-dependent enzymes of fold type I was calculated with the HMMER v3.0 package . A profile HMM is normally calculated from a multiple alignment of a set of appropriately selected sequences belonging to the superfamily to be modeled. If the alignment is sufficiently accurate, the model should be able to recognize distantly related members of the same superfamily. Proteins belonging to the fold type I superfamily are characterized by sharing scant sequence identity, as low as 9% . This characteristic provides the opportunity to design a multiple sequence alignment able to encode the structural fingerprint shared by all the members, even very distant, of the superfamily. However, accurate calculation of a multiple sequence alignment containing distant sequences is intrinsically rather difficult and strongly error-prone. This inherent difficulty has been surmounted through use of the structure-based sequence alignments  implemented in Combinatorial Extension  tool.
The program “hmmsearch” of the HMMER v3.0 package scanned the UniProt bacterial subset (version April, 2012) using the PLP_domain query profile. The search output was filtered for the purpose of selecting multimodular sequences containing genuine PLP-domains of fold type I. As an initial and rough criterion, only hits embedded in sequences longer than 1000 residue and overlapping the query HMM profile for at least 110 residues were taken into consideration. The sequence length threshold assured that only enzymes within multidomain contexts were collected since the typical length of a stand-alone PLP enzyme is around 400 residues. The sequence coverage threshold was set taking into account that the total residues involved in the SCRs of the PLP_domain profile are 108. The limit should therefore assure presence of the SCR residues in the significant hits. Candidate sequences to be submitted to the following steps were collected;
The set of sequences selected at the end of step 1 was further processed by the scripts “hmmscan” in the HMMER package and “Pfam_scan.pl” available at the Pfam site . The scripts compare a sequence against the Pfam database to locate known domains. Only sequences containing at least one phosphopantetheine (PP)-binding domain combined with a type I PLP-domain were retained;
The sequences of the type I PLP domains surviving the selection step 2, were extracted from the parent sequences according to the position boundaries assigned to them by the Pfam models. The excised sequences were multiply aligned using the program hmmalign  and a new HMM profile was calculated. This profile is, by definition, specific for recognizing type I domains embedded in multidomain environments containing at least one PP-binding motif. A new search has been carried with this query profile;
Search output was filtered with the following criteria: hits embedded in sequences longer than 500 residues and overlapping the query profile at least 110 residues, were retained; domain labeling was again carried out and only sequences characterized by the co-presence of PP-binding and type I domains were considered. The minimum length was set taking into consideration that the sum of the lengths of a typical CP domain and a typical type I enzymes is about 500 residues.
Domain annotation and assembly characterization
The organization and the identity of the domains contained in the parent sequences from which the type I domains were extracted, were determined through annotation by means of the script Pfam_scan.pl. Gene cluster identity, substrate specificity and product prediction were carried out with the software pipeline antiSMASH v2 . The type I domains isolated from the parent sequences were also compared to the families contained in the B6 database and assigned to one of them .
Phylogenetic analyses were carried out with the software MEGA v5.1 . The Minimum Evolution method  implemented in the program was applied. Evolutionary distances were computed using the JTT matrix-based method  and all positions containing gaps and missing data were eliminated. A Bootstrap test with 1000 replicates was used to assess the predicted topology of the resulting trees. All the calculated trees were unrooted.
Protein structure analysis, modelling and docking
Protein structure superposition and inspection utilized the Combinatorial Extension  and PyMOL graphics program , respectively. Most of the data handling was carried out with Python or Perl scripts under the Linux environment or utilized routines of the EMBOSS suite . Homology modeling made use of Modeller v.9.9  and PyMod ; models were validated with standard programs such as ProsaII  and Procheck . Candidate templates for homology modeling were assessed through the Phyre v2.0 server for protein fold recognition . Molecular docking relied on Molegro Virtual Docker . Candidate molecules were prepared by adding explicit hydrogens, charges, and flexible torsions. The side-chains of the active-site residues were kept fixed during docking. A spherical energy grid with a 15 Å radius, centered on the carboxylate moiety of the aspartate residue interacting with the pyridine nitrogen atom of pyridoxal 5′-phosphate, and a grid resolution of 0.30 Å was used. Other parameters were set at their default values: scoring function, MolDock score; search algorithm: MolDock SE; number of runs: 10; maximum iterations: 1500; maximum population size, 50; maximum number of poses returned, 5; cluster similar poses with RMSD threshold: 1.00 Å. Only the highest scoring pose, according to the MolDock scheme, was kept. ESPript  produced figures of sequence alignments while PyMOL  was the program for protein structure analysis and figure design. Prediction of the potential presence of protein-protein interaction sites was carried out with the consensus method implemented in meta-PPISP at the web site http://pipe.scs.fsu.edu/meta-ppisp . Protein-protein docking was carried out using the ClusPro method available at the server http://cluspro.bu.edu .
This work has been partially funded by the Italian MIUR (Ministero dell’Istruzione, Università, Ricerca). The work will be submitted by TM in partial fulfillment of the requirements of the degree of “Dottorato di Ricerca in Biofisica” at Sapienza, Università di Roma. Authors are grateful to Dr. Maria Rosaria Fullone for critically reading the manuscript and helpful discussions.
- Schneider G, Kack H, Lindqvist Y: The manifold of vitamin B6 dependent enzymes. Structure 2000, 8: R1-R6.View ArticlePubMedGoogle Scholar
- Clayton PT: B6-responsive disorders: a model of vitamin dependency. J Inherit Metab Dis 2006, 29: 317–326.View ArticlePubMedGoogle Scholar
- Eliot AC, Kirsch JF: Pyridoxal phosphate enzymes: mechanistic, structural, and evolutionary considerations. Annu Rev Biochem 2004, 73: 383–415.View ArticlePubMedGoogle Scholar
- Denessiouk KA, Denesyuk AI, Lehtonen JV, Korpela T, Johnson MS: Common structural elements in the architecture of the cofactor-binding domains in unrelated families of pyridoxal phosphate-dependent enzymes. Proteins 1999, 35: 250–261.View ArticlePubMedGoogle Scholar
- Paiardini A, Bossa F, Pascarella S: Evolutionarily conserved regions and hydrophobic contacts at the superfamily level: the case of the fold-type I, pyridoxal-5′-phosphate-dependent enzymes. Protein Sci 2004, 13: 2992–3005.PubMed CentralView ArticlePubMedGoogle Scholar
- Kirsch JF, Eichele G, Ford GC, Vincent MG, Jansonius JN, Gehring H, Christen P: Mechanism of action of aspartate aminotransferase proposed on the basis of its spatial structure. J Mol Biol 1984, 174: 497–525.View ArticlePubMedGoogle Scholar
- Rigali S, Derouaux A, Giannotta F, Dusart J: Subdivision of the helix-turn-helix GntR family of bacterial regulators in the FadR, HutC, MocR, and YtrA subfamilies. J Biol Chem 2002, 277: 12507–12515.View ArticlePubMedGoogle Scholar
- Bramucci E, Milano T, Pascarella S: Genomic distribution and heterogeneity of MocR-like transcriptional factors containing a domain belonging to the superfamily of the pyridoxal-5′-phosphate dependent enzymes of fold type I. Biochem Biophys Res Commun 2011, 415: 88–93.View ArticlePubMedGoogle Scholar
- Fischbach MA, Walsh CT: Assembly-line enzymology for polyketide and nonribosomal peptide antibiotics: logic, machinery, and mechanisms. Chem Rev 2006, 106: 3468–3496.View ArticlePubMedGoogle Scholar
- Felnagle EA, Jackson EE, Chan YA, Podevels AM, Berti AD, McMahon MD, Thomas MG: Nonribosomal peptide synthetases involved in the production of medically relevant natural products. Mol Pharm 2008, 5: 191–211.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang W, Ostash B, Walsh CT: Identification of the biosynthetic gene cluster for the pacidamycin group of peptidyl nucleoside antibiotics. Proc Natl Acad Sci USA 2010, 107: 16828–16833.PubMed CentralView ArticlePubMedGoogle Scholar
- Gerber R, Lou L, Du L: A PLP-dependent polyketide chain releasing mechanism in the biosynthesis of mycotoxin fumonisins in Fusarium verticillioides . J Am Chem Soc 2009, 131: 3148–3149.View ArticlePubMedGoogle Scholar
- Silakowski B, Kunze B, Nordsiek G, Blocker H, Hofle G, Muller R: The myxochelin iron transport regulon of the myxobacterium Stigmatella aurantiaca Sg a15. Eur J Biochem 2000, 267: 6476–6485.View ArticlePubMedGoogle Scholar
- Masschelein J, Mattheus W, Gao LJ, Moons P, Van Houdt R, Uytterhoeven B, Lamberigts C, Lescrinier E, Rozenski J, Herdewijn P, Aertsen A, Michiels C, Lavigne R: A PKS/NRPS/FAS hybrid gene cluster from Serratia plymuthica RVH1 encoding the biosynthesis of three broad spectrum, zeamine-related antibiotics. PLoS One 2013, 8: e54143.PubMed CentralView ArticlePubMedGoogle Scholar
- Aron ZD, Dorrestein PC, Blackhall JR, Kelleher NL, Walsh CT: Characterization of a new tailoring domain in polyketide biogenesis: the amine transferase domain of MycA in the mycosubtilin gene cluster. J Am Chem Soc 2005, 127: 14986–14987.View ArticlePubMedGoogle Scholar
- Garneau-Tsodikova S, Dorrestein PC, Kelleher NL, Walsh CT: Protein assembly line components in prodigiosin biosynthesis: characterization of PigA, G, H, I, J. J Am Chem Soc 2006, 128: 12600–12601.View ArticlePubMedGoogle Scholar
- Paiardini A, Bossa F, Pascarella S: CAMPO, SCR_FIND and CHC_FIND: a suite of web tools for computational structural biology. Nucleic Acids Res 2005, 33: W50-W55.PubMed CentralView ArticlePubMedGoogle Scholar
- Olson SA: EMBOSS opens up sequence analysis. European Molecular Biology Open Software Suite. Brief Bioinform 2002, 3: 87–91.View ArticlePubMedGoogle Scholar
- Mehta PK, Christen P: The molecular evolution of pyridoxal-5′-phosphate-dependent enzymes. Adv Enzymol Relat Areas Mol Biol 2000, 74: 129–184.PubMedGoogle Scholar
- Percudani R, Peracchi A: The B6 database: a tool for the description and classification of vitamin B6-dependent enzymatic activities and of the corresponding protein families. BMC Bioinforma 2009, 10: 273.View ArticleGoogle Scholar
- Medema MH, Blin K, Cimermancic P, de Jager V, Zakrzewski P, Fischbach MA, Weber T, Takano E, Breitling R: antiSMASH: rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res 2011, 39: W339-W346.PubMed CentralView ArticlePubMedGoogle Scholar
- Blin K, Medema MH, Kazempour D, Fischbach MA, Breitling R, Takano E, Weber T: antiSMASH 2.0–a versatile platform for genome mining of secondary metabolite producers. Nucleic Acids Res 2013, 41: W204-W212.PubMed CentralView ArticlePubMedGoogle Scholar
- Tang GL, Cheng YQ, Shen B: Leinamycin biosynthesis revealing unprecedented architectural complexity for a hybrid polyketide synthase and nonribosomal peptide synthetase. Chem Biol 2004, 11: 33–45.View ArticlePubMedGoogle Scholar
- Cheng YQ, Tang GL, Shen B: Type I polyketide synthase requiring a discrete acyltransferase for polyketide biosynthesis. Proc Natl Acad Sci USA 2003, 100: 3149–3154.PubMed CentralView ArticlePubMedGoogle Scholar
- Schulze JO, Schubert WD, Moser J, Jahn D, Heinz DW: Evolutionary relationship between initial enzymes of tetrapyrrole biosynthesis. J Mol Biol 2006, 358: 1212–1220.View ArticlePubMedGoogle Scholar
- Ikushiro H, Islam MM, Okamoto A, Hoseki J, Murakawa T, Fujii S, Miyahara I, Hayashi H: Structural insights into the enzymatic mechanism of serine palmitoyltransferase from Sphingobacterium multivorum . J Biochem 2009, 146: 549–562.View ArticlePubMedGoogle Scholar
- Lowther J, Charmier G, Raman MC, Ikushiro H, Hayashi H, Campopiano DJ: Role of a conserved arginine residue during catalysis in serine palmitoyltransferase. FEBS Lett 2011, 585: 1729–1734.View ArticlePubMedGoogle Scholar
- Milic D, Demidkina TV, Faleev NG, Matkovic-Calogovic D, Antson AA: Insights into the catalytic mechanism of tyrosine phenol-lyase from X-ray structures of quinonoid intermediates. J Biol Chem 2008, 283: 29206–29214.PubMed CentralView ArticlePubMedGoogle Scholar
- Fenalti G, Law RH, Buckle AM, Langendorf C, Tuck K, Rosado CJ, Faux NG, Mahmood K, Hampe CS, Banga JP, et al.: GABA production by glutamic acid decarboxylase is regulated by a dynamic catalytic loop. Nat Struct Mol Biol 2007, 14: 280–286.View ArticlePubMedGoogle Scholar
- Kozakov D, Hall DR, Beglov D, Brenke R, Comeau SR, Shen Y, Li K, Zheng J, Vakili P, Paschalidis I, Vajda S: Achieving reliability and high accuracy in automated protein docking: ClusPro, PIPER, SDU, and stability analysis in CAPRI rounds 13–19. Proteins 2010, 78: 3124–3130.PubMed CentralView ArticlePubMedGoogle Scholar
- Giardina G, Montioli R, Gianni S, Cellini B, Paiardini A, Voltattorni CB, Cutruzzola F: Open conformation of human DOPA decarboxylase reveals the mechanism of PLP addition to Group II decarboxylases. Proc Natl Acad Sci USA 2011, 108: 20514–20519.PubMed CentralView ArticlePubMedGoogle Scholar
- Hennig M, Grimm B, Contestabile R, John RA, Jansonius JN: Crystal structure of glutamate-1-semialdehyde aminomutase: an alpha2-dimeric vitamin B6-dependent enzyme with asymmetry in structure and active site reactivity. Proc Natl Acad Sci USA 1997, 94: 4866–4871.PubMed CentralView ArticlePubMedGoogle Scholar
- Silakowski B, Nordsiek G, Kunze B, Blocker H, Muller R: Novel features in a combined polyketide synthase/non-ribosomal peptide synthetase: the myxalamid biosynthetic gene cluster of the myxobacterium Stigmatella aurantiaca Sga15. Chem Biol 2001, 8: 59–69.View ArticlePubMedGoogle Scholar
- Alexeev D, Alexeeva M, Baxter RL, Campopiano DJ, Webster SP, Sawyer L: The crystal structure of 8-amino-7-oxononanoate synthase: a bacterial PLP-dependent, acyl-CoA-condensing enzyme. J Mol Biol 1998, 284: 401–419.View ArticlePubMedGoogle Scholar
- Phillips RS, Demidkina TV, Faleev NG: Structure and mechanism of tryptophan indole-lyase and tyrosine phenol-lyase. Biochim Biophys Acta 2003, 1647: 167–172.View ArticlePubMedGoogle Scholar
- Liu P, Torrens-Spence MP, Ding H, Christensen BM, Li J: Mechanism of cysteine-dependent inactivation of aspartate/glutamate/cysteine sulfinic acid alpha-decarboxylases. Amino Acids 2013, 44: 391–404.View ArticlePubMedGoogle Scholar
- Skoldberg F, Rorsman F, Perheentupa J, Landin-Olsson M, Husebye ES, Gustafsson J, Kampe O: Analysis of antibody reactivity against cysteine sulfinic acid decarboxylase, a pyridoxal phosphate-dependent enzyme, in endocrine autoimmune disease. J Clin Endocrinol Metab 2004, 89: 1636–1640.View ArticlePubMedGoogle Scholar
- Edwards DJ, Marquez BL, Nogle LM, McPhail K, Goeger DE, Roberts MA, Gerwick WH: Structure and biosynthesis of the jamaicamides, new mixed polyketide-peptide neurotoxins from the marine cyanobacterium Lyngbya majuscula . Chem Biol 2004, 11: 817–833.View ArticlePubMedGoogle Scholar
- McPhalen CA, Vincent MG, Jansonius JN: X-ray structure refinement and comparison of three forms of mitochondrial aspartate aminotransferase. J Mol Biol 1992, 225: 495–517.View ArticlePubMedGoogle Scholar
- Sieber SA, Linne U, Hillson NJ, Roche E, Walsh CT, Marahiel MA: Evidence for a monomeric structure of nonribosomal peptide synthetases. Chem Biol 2002, 9: 997–1008.View ArticlePubMedGoogle Scholar
- Hillson NJ, Walsh CT: Dimeric structure of the six-domain VibF subunit of vibriobactin synthetase: mutant domain activity regain and ultracentrifugation studies. Biochemistry 2003, 42: 766–775.View ArticlePubMedGoogle Scholar
- Richter CD, Nietlispach D, Broadhurst RW, Weissman KJ: Multienzyme docking in hybrid megasynthetases. Nat Chem Biol 2008, 4: 75–81.View ArticlePubMedGoogle Scholar
- Uniprot Consortium: Update on activities at the Universal Protein Resource(UniProt) in 2013. Nucleic Acids Res 2013, 41(Database issue):D43–47.View ArticleGoogle Scholar
- Rose PW, Beran B, Bi C, Bluhm WF, Dimitropoulos D, Goodsell DS, Prlic A, Quesada M, Quinn GB, Westbrook JD, et al.: The RCSB Protein Data Bank: redesigned web site and web services. Nucleic Acids Res 2011, 39: D392-D401.PubMed CentralView ArticlePubMedGoogle Scholar
- Finn RD, Clements J, Eddy SR: HMMER web server: interactive sequence similarity searching. Nucleic Acids Res 2011, 39: W29-W37.PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer 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–3402.PubMed CentralView ArticlePubMedGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al.: Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23: 2947–2948.View ArticlePubMedGoogle Scholar
- Katoh K, Kuma K, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res 2005, 33: 511–518.PubMed CentralView ArticlePubMedGoogle Scholar
- Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ: Jalview Version 2–a multiple sequence alignment editor and analysis workbench. Bioinformatics 2009, 25: 1189–1191.PubMed CentralView ArticlePubMedGoogle Scholar
- Gouy M, Guindon S, Gascuel O: SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol 2010, 27: 221–224.View ArticlePubMedGoogle Scholar
- Shindyalov IN, Bourne PE: Protein structure alignment by incremental combinatorial extension (CE) of the optimal path. Protein Eng 1998, 11: 739–747.View ArticlePubMedGoogle Scholar
- Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, et al.: The Pfam protein families database. Nucleic Acids Res 2012, 40: D290-D301.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 24: 1596–1599.View ArticlePubMedGoogle Scholar
- Rzhetsky A, Nei M: Statistical properties of the ordinary least-squares, generalized least-squares, and minimum-evolution methods of phylogenetic inference. J Mol Evol 1992, 35: 367–375.View ArticlePubMedGoogle Scholar
- Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci 1992, 8: 275–282.PubMedGoogle Scholar
- Schrodinger LLC: The PyMOL Molecular Graphics System, Version 184.108.40.206. Portland, OR, USA; 2013.Google Scholar
- Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol 1993, 234: 779–815.View ArticlePubMedGoogle Scholar
- Bramucci E, Paiardini A, Bossa F, Pascarella S: PyMod: sequence similarity searches, multiple sequence-structure alignments, and homology modeling within PyMOL. BMC Bioinforma 2012, 13(4):S2.View ArticleGoogle Scholar
- 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: W407-W410.PubMed CentralView ArticlePubMedGoogle Scholar
- Laskowski RA, Rullmannn JA, MacArthur MW, Kaptein R, Thornton JM: AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR. J Biomol NMR 1996, 8: 477–486.View ArticlePubMedGoogle Scholar
- Bennett-Lovsey RM, Herbert AD, Sternberg MJ, Kelley LA: Exploring the extremes of sequence/structure space with ensemble fold recognition in the program Phyre. Proteins 2008, 70: 611–625.View ArticlePubMedGoogle Scholar
- Thomsen R, Christensen MH: MolDock: a new technique for high-accuracy molecular docking. J Med Chem 2006, 49: 3315–3321.View ArticlePubMedGoogle Scholar
- Gouet P, Courcelle E, Stuart DI, Metoz F: ESPript: analysis of multiple sequence alignments in PostScript. Bioinformatics 1999, 15: 305–308.View ArticlePubMedGoogle Scholar
- Qin S, Zhou HX: meta-PPISP: a meta web server for protein-protein interaction site prediction. Bioinformatics 2007, 23: 3386–3387.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.