Armadillo: domain boundary prediction by amino acid composition

The identification and annotation of protein domains provides a critical step in the accurate determination of molecular function. Both computational and experimental methods of protein structure determination may be deterred by large multi-domain proteins or flexible linker regions. Knowledge of domains and their boundaries may reduce the experimental cost of protein structure determination by allowing researchers to work on a set of smaller and possibly more successful alternatives. Current domain prediction methods often rely on sequence similarity to conserved domains and as such are poorly suited to detect domain structure in poorly conserved or orphan proteins. We present here a simple computational method to identify protein domain linkers and their boundaries from sequence information alone. Our domain predictor, Armadillo (http://armadillo.blueprint.org), uses any amino acid index to convert a protein sequence to a smoothed numeric profile from which domains and domain boundaries may be predicted. We derived an amino acid index called the domain linker propensity index (DLI) from the amino acid composition of domain linkers using a non-redundant structure dataset. The index indicates that Pro and Gly show a propensity for linker residues while small hydrophobic residues do not. Armadillo predicts domain linker boundaries from Z-score distributions and obtains 35% sensitivity with DLI in a two-domain, single-linker dataset (within +/-20 residues from linker). The combination of DLI and an entropy-based amino acid index increases the overall Armadillo sensitivity to 56% for two domain proteins. Moreover, Armadillo achieves 37% sensitivity for multi-domain proteins, surpassing most other prediction methods. Armadillo provides a simple, but effective method by which prediction of domain boundaries can be obtained with reasonable sensitivity. Armadillo should prove to be a valuable tool for rapidly delineating protein domains in poorly conserved proteins or those with no sequence neighbors. As a first-line predictor, domain meta-predictors could yield improved results with Armadillo predictions.     

The linear pentadecapeptide gramicidin is assembled by four multimodular nonribosomal peptide synthetases that comprise 16 modules with 56 catalytic domains

Linear gramicidin is a membrane channel forming pentadecapeptide that is produced via the nonribosomal pathway. It consists of 15 hydrophobic amino acids with alternating l- and d-configuration forming a beta-helix-like structure. It has an N-formylated valine and a C-terminal ethanolamine. Here we report cloning and sequencing of the entire biosynthetic gene cluster as well as initial biochemical analysis of a new reductase domain. The biosynthetic gene cluster was identified on two nonoverlapping fosmids and a 13-kilobase pair (kbp) interbridge fragment covering a region of 74 kbp. Four very large open reading frames, lgrA, lgrB, lgrC, and lgrD with 6.8, 15.5, 23.3, and 15.3 kbp, were identified and shown to encode nonribosomal peptide synthetases with two, four, six, and four modules, respectively. Within the 16 modules identified, seven epimerization domains in alternating positions were detected as well as a putative formylation domain fused to the first module LgrA and a putative reductase domain attached to the C-terminal module of LgrD. Analysis of the substrate specificity by phylogenetic studies using the residues of the substrate-binding pockets of all 16 adenylation domains revealed a good agreement of the substrate amino acids predicted with the sequence of linear gramicidin. Additional biochemical analysis of the three adenylation domains of modules 1, 2, and 3 confirmed the colinearity of this nonribosomal peptide synthetase assembly line. Module 16 was predicted to activate glycine, which would then, being the C-terminal residue of the peptide chain, be reduced by the adjacent reductase domain to give ethanolamine, thereby releasing the final product N-formyl-pentadecapeptide-ethanolamine. However, initial biochemical analysis of this reductase showed only a one-step reduction yielding the corresponding aldehyde in vitro.     

Novel Src homology 3 domain-binding motifs identified from proteomic screen of a Pro-rich region

The Src homology (SH) 3 domain has been shown recently to bind peptide sequences that lack the canonical PXXP motif. The diverse specificity in ligand recognition for a group of 15 SH3 domains has now been investigated using arrays of peptides derived from the proline-rich region of the SH2 domain-containing leukocyte protein of 76 kDa (SLP-76). A screen of the peptide arrays using individual or mixed SH3 domains has allowed the identification of a number of candidate SH3-binding peptides. Although some peptides contain the conventional PXXP motif, most are devoid of such a motif and are instead enriched in basic residues. Fluorescent polarization measurements using soluble peptides and purified SH3 domains demonstrated that several SH3 domains, including those from growth factor receptor-bound protein 2 (Grb2), NCK, and phospholipase C (PLC)-gamma1, bound with moderate affinities (10-100 microm) to a group of non-conventional peptides. Of particular interest, the PLC-gamma1 SH3 domain was found to associate with SLP-76 through at least three distinct sites, two of which bore a novel KKPP motif and the other contained the classic PXXP sequence. Intriguingly mutation of critical residues for the three sites not only affected binding of SLP-76 to the PLC-gamma1 SH3 domain but also to the Grb2 C-terminal SH3 domain, indicating that the binding sites in SLP-76 for the two SH3 domains are overlapped. Our studies suggest that the SH3 domain is an inherently promiscuous interaction module capable of binding to peptides that may or may not contain a PXXP motif. Furthermore the identification of numerous non-conventional SH3-binding peptides in SLP-76 implies that the global ligand pool for SH3 domains in a mammalian proteome may be significantly greater than previously acknowledged.     

The tyrocidine biosynthesis operon of Bacillus brevis: complete nucleotide sequence and biochemical characterization of functional internal adenylation domains.

The cyclic decapeptide antibiotic tyrocidine is produced by Bacillus brevis ATCC 8185 on an enzyme complex comprising three peptide synthetases, TycA, TycB, and TycC (tyrocidine synthetases 1, 2, and 3), via the nonribosomal pathway. However, previous molecular characterization of the tyrocidine synthetase-encoding operon was restricted to tycA, the gene that encodes the first one-module-bearing peptide synthetase. Here, we report the cloning and sequencing of the entire tyrocidine biosynthesis operon (39.5 kb) containing the tycA, tycB, and tycC genes. As deduced from the sequence data, TycB (404,562 Da) consists of three modules, including an epimerization domain, whereas TycC (723,577 Da) is composed of six modules and harbors a putative thioesterase domain at its C-terminal end. Each module incorporates one amino acid into the peptide product and can be further subdivided into domains responsible for substrate adenylation, thiolation, condensation, and epimerization (optional). We defined, cloned, and expressed in Escherichia coli five internal adenylation domains of TycB and TycC. Soluble His6-tagged proteins, ranging from 536 to 559 amino acids, were affinity purified and found to be active by amino acid-dependent ATP-PPi exchange assay. The detected amino acid specificities of the investigated domains manifested the colinear arrangement of the peptide product with the respective module in the corresponding peptide synthetases and explain the production of the four known naturally occurring tyrocidine variants. The Km values of the investigated adenylation domains for their amino acid substrates were found to be comparable to those published for undissected wild-type enzymes. These findings strongly support the functional integrities of single domains within multifunctional peptide synthetases. Directly downstream of the 3' end of the tycC gene, and probably transcribed in the tyrocidine operon, two tandem ABC transporters, which may be involved in conferring resistance against tyrocidine, and a putative thioesterase were found.     

Structure of antibody-bound peptides and retro-inverso analogues. A transferred nuclear Overhauser effect spectroscopy and molecular dynamics approach

The three-dimensional structures of the two L-peptides, H-CGGIRGERA-OH, called L(A), and H-CGGIRGERG-OH, called L(G), corresponding or close to the IRGERA sequence present in the C-terminal region (residues 130-135) of histone H3, and their retro-inverso analogues HO-mAreGriGGC-NH2, called RI(mA), and HO-mGreGriGGC-NH2, called RI(mG), have been studied by two-dimensional 1H NMR and molecular dynamics calculations in association with a monoclonal antibody generated against L(A). At 25 degrees C, the affinity constants of the monoclonal antibody with respect to RI(mA) and RI(mG) were 75- and 270-fold higher than those measured with the homologous L(A) and L(G) peptides, respectively. Due to the spontaneous epimerization of the mA malonic residue, RI(mA) gave rise to two sets of resonances. With regard to the NH amide region, one set was similar to that for RI(mG) while the second was similar to those for the parent L-peptides L(A) and L(G). The antibody-bound conformations of the two couples of L- and retro-inverso peptides have been analyzed using molecular modeling calculations based on the transferred NOE interproton distances. Folded structures appeared in both cases with a type II' beta-turn in the parent GGIR sequence and a type I' beta-turn in the retro-inverso reGr sequence.     

Fine-tuning of protein domain boundary by minimizing potential coiled coil regions

Structural determination of individual protein domains isolated from multidomain proteins is a common approach in the post-genomic era. Novel and thus uncharacterized domains liberated from intact proteins often self-associate due to incorrectly defined domain boundaries. Self-association results in missing signals, poor signal dispersion and a low signal-to-noise ratio in ¹H–¹⁵N HSQC spectra. We have found that a putative, non-canonical coiled coil region close to a domain boundary can cause transient hydrophobic self-association and monomer–dimer equilibrium in solution. Here we propose a rational method to predict putative coiled coil regions adjacent to the globular core domain using the program COILS. Except for the amino acid sequence, no preexisting knowledge concerning the domain is required. A small number of mutant proteins with a minimized coiled coil region have been rationally designed and tested. The engineered domains exhibit decreased self-association as assessed by ¹H–¹⁵N HSQC spectra with improved peak dispersion and sharper cross peaks. Two successful examples of isolating novel N-terminal domains from AAA-ATPases are demonstrated. Our method is useful for the experimental determination of domain boundaries suited for structural genomics studies. Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .     

EVEREST: automatic identification and classification of protein domains in all protein sequences

Background  

Proteins are comprised of one or several building blocks, known as domains. Such domains can be classified into families according to their evolutionary origin. Whereas sequencing technologies have advanced immensely in recent years, there are no matching computational methodologies for large-scale determination of protein domains and their boundaries. We provide and rigorously evaluate a novel set of domain families that is automatically generated from sequence data. Our domain family identification process, called EVEREST (EVolutionary Ensembles of REcurrent SegmenTs), begins by constructing a library of protein segments that emerge in an all vs. all pairwise sequence comparison. It then proceeds to cluster these segments into putative domain families. The selection of the best putative families is done using machine learning techniques. A statistical model is then created for each of the chosen families. This procedure is then iterated: the aforementioned statistical models are used to scan all protein sequences, to recreate a library of segments and to cluster them again.     

The influence of cholesterol on the interaction of HIV gp41 membrane proximal region-derived peptides with lipid bilayers

A small amino acid sequence (LWYIK) inside the HIV-1 gp41 ectodomain membrane proximal region (MPR) is commonly referred to as a cholesterol-binding domain. To further study this unique and peculiar property we have used fluorescence spectroscopy techniques to unravel the membrane interaction properties of three MPR-derived synthetic peptides: the membrane proximal region peptide-complete (MPRP-C) which corresponds to the complete MPR; the membrane proximal region peptide-short (MPRP-S), which corresponds to the last five MPR amino acid residues (the putative cholesterol-binding domain) and the membrane proximal region peptide-intermediate (MPRP-I), which corresponds to the MPRP-C peptide without the MPRP-S sequence. MPRP-C and MPRP-I membrane interaction is largely independent of the membrane phase. Membrane interaction of MPRP-S occurs for fluid phase membranes but not in gel phase membranes or cholesterol-containing bilayers. The gp41 ectodomain MPR may have a very specific function in viral fusion through the concerted and combined action of cholesterol-binding and non-cholesterol-binding domains (i.e. domains corresponding to MPRP-S and MPRP-I, respectively).     

Cloning, Expression, and Biological Function of a dTDP-Deoxyglucose Epimerase (gerF) Gene from Streptomyces sp. GERI-155

GERI-155 is a macrolide antibiotic containing two deoxyhexose molecules which has antimicrobial activities against gram-positive bacteria. The deoxyhexose biosynthetic gene cluster of GERI-155 from Streptomyces sp. GERI-155 genome has now been isolated. Four orf were identified and a putative orf, supposed to code for the dTDP-deoxyglucose epimerase gene, was designated as gerF. gerF was expressed in E. coli using recombinant expression vector pHJ3. The recombinant protein expressed in a soluble form. The enzyme was purified by Ni-affinity column using imidazole buffer as eluents. The molecular mass of the expressed protein correlated with the predicted mass (36,000 Da) deduced from the cloned gene sequence data. The purified enzyme produced maltol from dTDP-4-keto-6-deoxyglucose and it was confirmed that the expressed protein was dTDP-deoxyglucose epimerase catalyzing epimerization of C-3 and C-5 or C-3 of dTDP-4-keto-6-deoxyglucose.     

A beta-lysine adenylating enzyme and a beta-lysine binding protein involved in poly beta-lysine chain assembly in nourseothricin synthesis in Streptomyces noursei.

Nourseothricins (syn. Streptothricins), a group of nucleoside peptides produced by several streptomycete strains, contain a poly beta-lysine chain of variable length attached in amide linkage to the amino sugar moiety gulosamine of the nucleoside portion. We show that the nourseothricin-producing Streptomyces noursei contains an enzyme (NpsA) of an apparent M(r) 56,000 that specifically activates beta-lysine by adenylation but does not bind to it as a thioester. Cloning and sequencing of npsA from S. noursei including its flanking DNA regions revealed that it is closely linked to the nourseothricin resistance gene nat1 and some other genes on the chromosome possibly involved in nourseothricin biosynthesis. The deduced amino-acid sequence revealed that NpsA is a stand-alone adenylation domain with similarity to the adenylation domains of nonribosomal peptide synthetases (NRPS). Further analysis revealed that S. noursei contains a beta-lysine binding enzyme (NpsB) of about M(r) 64,100 which can be loaded by NpsA with beta-lysine as a thioester. Analysis of the deduced amino-acid sequence from the gene (npsB) of NpsB showed that it consists of two domains. The N-terminal domain of approximately 100 amino-acid residues has high similarity to PCP domains of NRPSs whereas the 450-amino-acid C-terminal domain has a high similarity to epimerization (E)-domains of NRPSs. Remarkably, in this E-domain the conserved H-H-motif is changed to H-Q, which suggests that either the domain is nonfunctional or has a specialized function. The presence of one single adenylating beta-lysine activating enzyme in nourseothricin-producing streptomycete and a separate binding protein suggests an iteratively operating NRPS-module catalyses synthesis of the poly beta-lysine chain.     

Characterization of the functional domains of Escherichia coli RNase II

RNase II is a single-stranded-specific 3'-exoribonuclease that degrades RNA generating 5'-mononucleotides. This enzyme is the prototype of an ubiquitous family of enzymes that are crucial in RNA metabolism and share a similar domain organization. By sequence prediction, three different domains have been assigned to the Escherichia coli RNase II: two RNA-binding domains at each end of the protein (CSD and S1), and a central RNB catalytic domain. In this work we have performed a functional characterization of these domains in order to address their role in the activity of RNase II. We have constructed a large set of RNase II truncated proteins and compared them to the wild-type regarding their exoribonucleolytic activity and RNA-binding ability. The dissociation constants were determined using different single- or double-stranded substrates. The results obtained revealed that S1 is the most important domain in the establishment of stable RNA-protein complexes, and its elimination results in a drastic reduction on RNA-binding ability. In addition, we also demonstrate that the N-terminal CSD plays a very specific role in RNase II, preventing a tight binding of the enzyme to single-stranded poly(A) chains. Moreover, the biochemical results obtained with RNB mutant that lacks both putative RNA-binding domains, revealed the presence of an additional region involved in RNA binding. Such region, was identified by sequence analysis and secondary structure prediction as a third putative RNA-binding domain located at the N-terminal part of RNB catalytic domain.     

Genes encoding synthetases of cyclic depsipeptides, anabaenopeptilides, in Anabaena strain 90

Anabaena strain 90 produces three hepatotoxic heptapeptides (microcystins), two seven-residue depsipeptides called anabaenopeptilide 90A and 90B, and three six-residue peptides called anabaenopeptins. The anabaenopeptilides belong to a group of cyanobacterial depsipeptides that share the structure of a six-amino-acid ring with a side-chain. Despite their similarity to known cyclic peptide toxins, no function has been assigned to the anabaenopeptilides. Degenerate oligonucleotide primers based on the conserved amino acid sequences of other peptide synthetases were used to amplify DNA from Anabaena 90, and the resulting polymerase chain reaction (PCR) products were used to identify a peptide synthetase gene cluster. Four genes encoding putative anabaenopeptilide synthetase domains were characterized. Three genes, apdA, apdB and apdD, contain two, four and one module, respectively, encoding a total of seven modules for activation and peptide bond formation of seven L-amino acids. Modules five and six also carry methyltransferase-like domains. Before the first module, there is a region similar in amino acid sequence to formyltransferases. A fourth gene (apdC), between modules six and seven, is similar in sequence to halogenase genes. Thus, the order of domains is co-linear with the positions of amino acid residues in the finished peptide. A mutant of Anabaena 90 was made by inserting a chloramphenicol resistance gene into the apdA gene. DNA amplification by PCR confirmed the insertion. Mass spectrometry analysis showed that anabaenopeptilides are not made in the mutant strain, but other peptides, such as microcystins and anabaenopeptins, are still produced by the mutant.     

Decoding protein-protein interactions through combinatorial chemistry: sequence specificity of SHP-1, SHP-2, and SHIP SH2 domains

A general, combinatorial library method for the rapid identification of high-affinity peptide ligands of protein modular domains is reported. The validity of this method has been demonstrated by determining the sequence specificity of four Src homology 2 (SH2) domains derived from protein tyrosine phosphatase SHP-1 and SHP-2 and inositol phosphatase SHIP. A phosphotyrosyl (pY) peptide library was screened against the SH2 domains, and the beads that carry high-affinity ligands of the SH2 domains were identified and peptides were sequenced by partial Edman degradation and mass spectrometry. The results reveal that the N-terminal SH2 domain of SHP-2 is capable of recognizing four different classes of pY peptides. Binding competition studies suggest that the four classes of pY peptides all bind to the same site on the SH2 domain surface. The C-terminal SH2 domains of SHP-1 and SHP-2 and the SHIP SH2 domain each bind to pY peptides of a single consensus sequence. Database searches using the consensus sequences identified most of the known as well as many potential interacting proteins of SHP-1 and/or SHP-2. Several proteins are found to bind to the SH2 domains of SHP-1 and SHP-2 through a new, nonclassical ITIM motif, (V/I/L)XpY(M/L/F)XP, which corresponds to the class IV peptides selected from the pY library. The combinatorial library method should be generally applicable to other protein domains.     

Portability of epimerization domain and role of peptidyl carrier protein on epimerization activity in nonribosomal peptide synthetases.

Incorporation of nonproteinogenic amino acids in small polypeptides synthesized by nonribosomal peptide synthetases (NRPS) significantly contributes to their biological activity. In these peptides, conversion of L-amino acids to the corresponding D-isomer is catalyzed by specialized NRPS modules that utilize an epimerization (E) domain. To understand the basis for the specific interaction of E domains with PCP domains (peptidyl carrier proteins, also described as T domains) and to investigate their substrate tolerance, we constructed a set of eight fusion proteins. The gene fragments encoding E and PCP-E domains of TycA (A-PCP-E), the one module tyrocidine synthetase A, were fused to different gene fragments encoding A and A-PCP domains, resulting in A/PCP-E and A-PCP/E types of fusion proteins (slash indicates site of fusion). We were able to show that the E domain of TycA, usually epimerizing Phe, does also accept the alternate substrates Trp, Ile, and Val, although with reduced efficiency. Interestingly, however, an epimerization activity was only observed in the case of fusion proteins where the PCP domain originates from modules containing an E domain. Sequence comparison revealed that such PCPs possess significant differences in the signature Ppant binding motif (CoreT: [GGDSI]), when compared to those carrier proteins, originating from ordinary C-A-PCP elongation modules (CoreT: [GGHSL]). By means of mutational analysis, we could show that epimerization activity is influenced by the nature of amino acid residues in proximity to the cofactor Ppant binding site. The aspartate residue in front of the invariant serine (Ppant binding site) especially seems to play an important role for the proper interaction between PCP and the E domain, as well as the presentation of the aminoacyl-S-Ppant substrate in the course of substrate epimerization. In conclusion, specialized PCP domains are needed for a productive interaction with E domains when constructing hybrid enzymes.     

The two homologous domains of human angiotensin I-converting enzyme are both catalytically active.

Molecular cloning of human endothelial angiotensin I-converting enzyme (kininase II; EC 3.4.15.1) (ACE) has recently shown that the enzyme contains two large homologous domains (called here the N and C domains), each bearing a putative active site, identified by sequence comparisons with the active sites of other zinc metallopeptidases. However, the previous experiments with zinc or competitive ACE inhibitors suggested a single active site in ACE. To establish whether both domains of ACE are enzymatically active, a series of ACE mutants, each containing only one intact domain, were constructed by deletion or point mutations of putative critical residues of the other domain, and expressed in heterologous Chinese hamster ovary cells. Both domains are enzymatically active and cleave the C-terminal dipeptide of hippuryl-His-Leu or angiotensin I. Moreover, both domains have an absolute zinc requirement for activity, are activated by chloride and are sensitive to competitive ACE inhibitors, and appear to function independently. However, the two domains display different catalytic constants and different patterns of chloride activation. At high chloride concentrations, the C domain hydrolyzes the two substrates tested faster than does the N domain. His-361,365 and His-959,963 are established as essential residues in the N and C domains, respectively, most likely involved in zinc binding, and Glu-362 in the N domain and Glu-960 in the C domain are essential catalytic residues. These observations provide strong evidence that ACE possesses two independent catalytic domains and suggest that they may have different functions.     

Identification and characterization of a novel class of extracellular poly(3-hydroxybutyrate) depolymerase from Bacillus sp. strain NRRL B-14911

The catalytic, linker, and denatured poly(3-hydroxybutyrate) (dPHB)-binding domains of bacterial extracellular PHB depolymerases (PhaZs) are classified into several different types. We now report a novel class of extracellular PHB depolymerase from Bacillus sp. strain NRRL B-14911. Its catalytic domain belongs to type 1, whereas its putative linker region neither possesses the sequence features of the three known types of linker domains nor exhibits significant amino acid sequence similarity to them. Instead, this putative linker region can be divided into two distinct linker domains of novel types: LD1 and LD2. LD1 shows significant amino acid sequence similarity to certain regions of a large group of PHB depolymerase-unrelated proteins. LD2 and its homologs are present in a small group of PhaZs. The remaining C-terminal portion of this PhaZ can be further divided into two distinct domains: SBD1 and SBD2. Each domain showed strong binding to dPHB, and there is no significant sequence similarity between them. Each domain neither possesses the sequence features of the two known types of dPHB-binding domains nor shows significant amino acid sequence similarity to them. These unique features indicate the presence of two novel and distinct types of dPHB-binding domains. Homologs of these novel domains also are present in the extracellular PhaZ of Bacillus megaterium and the putative extracellular PhaZs of Bacillus pseudofirmus and Bacillus sp. strain SG-1. The Bacillus sp. NRRL B-14911 PhaZ appears to be a representative of a novel class of extracellular PHB depolymerases.     

Accurate Prediction of Peptide Binding Sites on Protein Surfaces

Many important protein–protein interactions are mediated by the binding of a short peptide stretch in one protein to a large globular segment in another. Recent efforts have provided hundreds of examples of new peptides binding to proteins for which a three-dimensional structure is available (either known experimentally or readily modeled) but where no structure of the protein–peptide complex is known. To address this gap, we present an approach that can accurately predict peptide binding sites on protein surfaces. For peptides known to bind a particular protein, the method predicts binding sites with great accuracy, and the specificity of the approach means that it can also be used to predict whether or not a putative or predicted peptide partner will bind. We used known protein–peptide complexes to derive preferences, in the form of spatial position specific scoring matrices, which describe the binding-site environment in globular proteins for each type of amino acid in bound peptides. We then scan the surface of a putative binding protein for sites for each of the amino acids present in a peptide partner and search for combinations of high-scoring amino acid sites that satisfy constraints deduced from the peptide sequence. The method performed well in a benchmark and largely agreed with experimental data mapping binding sites for several recently discovered interactions mediated by peptides, including RG-rich proteins with SMN domains, Epstein-Barr virus LMP1 with TRADD domains, DBC1 with Sir2, and the Ago hook with Argonaute PIWI domain. The method, and associated statistics, is an excellent tool for predicting and studying binding sites for newly discovered peptides mediating critical events in biology.     

Utility of epimerization domains for the redesign of nonribosomal peptide synthetases

Many pharmacologically important agents are assembled on multimodular nonribosomal peptide synthetases (NRPSs) whose modules comprise a set of core domains with all essential catalytic functions necessary for the incorporation and modification of one building block. Very often, d-amino acids are found in such products which, with few exceptions, are generated by the action of NRPS integrated epimerization (E) domains that alter the stereochemistry of the corresponding peptidyl carrier protein (PCP) bound l-intermediate. In this study we present a quantitative investigation of substrate specificity of four different E domains (two 'peptidyl-' and two 'aminoacyl-'E domains) derived from different NRPSs towards PCP bound peptides. The respective PCP-E bidomain apo-proteins (TycB(3)-, FenD(2)-, TycA- and GrsA-PCP-E) were primed with various peptidyl-CoA precursors by utilizing the promiscuous phosphopantetheinyl transferase Sfp. PCP bound peptidyl-S-Ppant epimerization products were chemically cleaved and analyzed for their l/d-ratios by LCMS. We were able to show that all four E domains tolerate a broad variety of peptidyl-S-Ppant-substrates as evaluated by k(obs) values and final l/d-product equilibria determined for each reaction. The two C-terminal amino acids of the substrate seem to be recognized by 'peptidyl-'E domains. Interestingly, the 'aminoacyl-'E domains GrsA- and TycA-E were also able to convert the elongated intermediates. All four E domains accepted an N-methylated precursor as well and epimerized this substrate with high efficiency. Finally, we could demonstrate that the condensation (C) domain of TycB(1) is also able to process peptidyl substrates transferred by TycA. In conclusion, these findings are of great impact on future engineering attempts.