The Lipase Engineering Database: a navigation and analysis tool for protein families

The Lipase Engineering Database (LED) (http://www.led.uni-stuttgart.de) integrates information on sequence, structure, and function of lipases, esterases, and related proteins. Sequence data on 806 protein entries are assigned to 38 homologous families, which are grouped into 16 superfamilies with no global sequence similarity between each other. For each family, multisequence alignments are provided with functionally relevant residues annotated. Pre-calculated phylogenetic trees allow navigation inside superfamilies. Experimental structures of 45 proteins are superposed and consistently annotated. The LED has been applied to systematically analyze sequence-structure-function relationships of this vast and diverse enzyme class. It is a useful tool to identify functionally relevant residues apart from the active site residues, and to design mutants with desired substrate specificity.     

The triterpene cyclase protein family: a systematic analysis

Triterpene cyclases catalyze a broad range of cyclization reactions to form polycyclic triterpenes. Triterpene cyclases that convert squalene to hopene are named squalene-hopene cyclases (SHC) and triterpene cyclases that convert oxidosqualene are named oxidosqualene cyclases (OSC). Many sequences have been published, but there is only one structure available for each of SHCs and OSCs. Although they catalyze a similar reaction, the sequence similarity between SHCs and OSCs is low. A family classification based on phylogenetic analysis revealed 20 homologous families which are grouped into two superfamilies, SHCs and OSCs. Based on this family assignment, the Triterpene Cyclase Engineering Database (TTCED) was established. It integrates available information on sequence and structure of 639 triterpene cyclases as well as on structurally and functionally relevant amino acids. Family specific multiple sequence alignments were generated to identify the functionally relevant residues. Based on sequence alignments, conserved residues in SHCs and OSCs were analyzed and compared to experimentally confirmed mutational data. Functional schematic models of the central cavities of OSCs and SHCs were derived from structure comparison and sequence conservation analysis. These models demonstrate the high similarity of the substrate binding cavity of SHCs and OSCs and the equivalences of the respective residues. The TTCED is a novel source for comprehensive information on the triterpene cyclase family, including a compilation of previously described mutational data. The schematic models present the conservation analysis in a readily available fashion and facilitate the correlation of residues to a specific function or substrate interaction.     

Conformational analysis and clustering of short and medium size loops connecting regular secondary structures: a database for modeling and prediction.

Loops are regions of nonrepetitive conformation connecting regular secondary structures. We identified 2,024 loops of one to eight residues in length, with acceptable main-chain bond lengths and peptide bond angles, from a database of 223 protein and protein-domain structures. Each loop is characterized by its sequence, main-chain conformation, and relative disposition of its bounding secondary structures as described by the separation between the tips of their axes and the angle between them. Loops, grouped according to their length and type of their bounding secondary structures, were superposed and clustered into 161 conformational classes, corresponding to 63% of all loops. Of these, 109 (51% of the loops) were populated by at least four nonhomologous loops or four loops sharing a low sequence identity. Another 52 classes, including 12% of the loops, were populated by at least three loops of low sequence similarity from three or fewer nonhomologous groups. Loop class suprafamilies resulting from variations in the termini of secondary structures are discussed in this article. Most previously described loop conformations were found among the classes. New classes included a 2:4 type IV hairpin, a helix-capping loop, and a loop that mediates dinucleotide-binding. The relative disposition of bounding secondary structures varies among loop classes, with some classes such as beta-hairpins being very restrictive. For each class, sequence preferences as key residues were identified; those most frequently at these conserved positions than in proteins were Gly, Asp, Pro, Phe, and Cys. Most of these residues are involved in stabilizing loop conformation, often through a positive phi conformation or secondary structure capping. Identification of helix-capping residues and beta-breakers among the highly conserved positions supported our decision to group loops according to their bounding secondary structures. Several of the identified loop classes were associated with specific functions, and all of the member loops had the same function; key residues were conserved for this purpose, as is the case for the parvalbumin-like calcium-binding loops. A significant number, but not all, of the member loops of other loop classes had the same function, as is the case for the helix-turn-helix DNA-binding loops. This article provides a systematic and coherent conformational classification of loops, covering a broad range of lengths and all four combinations of bounding secondary structure types, and supplies a useful basis for modelling of loop conformations where the bounding secondary structures are known or reliably predicted.     

Crystal structure of the CN-hydrolase SA0302 from the pathogenic bacterium Staphylococcus aureus belonging to the Nit and NitFhit Branch of the nitrilase superfamily

The nitrilases include a variety of enzymes with functional specificities of nitrilase, amidase, and hydrolase reactions. The crystal structure of the uncharacterized protein SA0302 from the pathogenic microorganism Staphylococcus aureus is solved at 1.7?Å resolution. The protein contains 261 amino acids and presents a four-layer αββα sandwich with a chain topology similar to that of a few known CN-hydrolase folds. In the crystal, the proteins are arranged as dimers whose monomers are related by a pseudo twofold rotation symmetry axis. Analysis of the sequences and structures of CN-hydrolases with known 3D structures shows that SA0302 definitely is a member of Branch 10 (Nit and NitFhit) of the nitrilase superfamily. Enzyme activities and substrate specificities of members of this branch are not yet characterized, in contrast to those of the members of Branches 1–9. Although the sequence identities between Branch 10 members are rather low, less than 30%, five conserved regions are common in this subfamily. Three of them contain functionally important catalytic residues, and the two other newly characterized ones are associated with crucial intramolecular and intermolecular interactions. Sequence homology of the area near the active site shows clearly that the catalytic triad of SA0302 is Glu41-Lys110-Cys146. We suggest also that the active site includes a fourth residue, the closely located Glu119. Despite an extensive similarity with other Nit-family structural folds, SA0302 displays an important difference. Protein loop 111–122, which follows the catalytic Lys110, is reduced to half the number of amino acids found in other Nit-family members. This leaves the active site fully accessible to solvent and substrates. We have identified conservative sequence motifs around the three core catalytic residues, which are inherent solely to Branch 10 of the nitrilase superfamily. On the basis of these new sequence fingerprints, 10 previously uncharacterized proteins also could be assigned to this hydrolase subfamily.

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:19     

Crystal structure of the YdjC-family protein TTHB029 from Thermus thermophilus HB8: structural relationship with peptidoglycan N-acetylglucosamine deacetylase

The YdjC-family protein is widely distributed, from human to bacteria, but so far no three-dimensional structure and functional analysis of this family of proteins has been reported. We determined the three-dimensional structure of the YdjC homolog TTHB029 at a resolution of 2.9 Å. The overall structure of the monomer consists of (βα)-barrel fold forming a homodimer. Asp21, His60, and His127 residues coordinate to Mg²⁺ as a possible active site. TTHB029 shows structural similarity to the peptidoglycan N-acetylglucosamine deacetylase from Streptococcus pneumoniae (SpPgdA). The active site groove of SpPgdA includes the Zn²⁺ coordinated to Asp276, His326, and His330. Despite the low sequence identity, metal-binding residues of Asp-His-His were conserved among the two enzymes. There were definitive differences, however, in that one of the histidines of the metal-binding site was substituted for the other histidine located on the other loop. Moreover, these important metal-binding residues and the residues of the presumed active site are fully conserved in YdjC-family protein.     

Sequence variation in the β7–β8 loop of bacterial class A sortase enzymes alters substrate selectivity

Gram-positive bacteria contain sortase enzymes on their cell surfaces that catalyze transpeptidation reactions critical for proper cellular function. In vitro, sortases are used in sortase-mediated ligation (SML) reactions for a variety of protein engineering applications. Historically, sortase A from Staphylococcus aureus (saSrtA) has been the enzyme of choice to catalyze SML reactions. However, the stringent specificity of saSrtA for the LPXTG sequence motif limits its uses. Here, we describe the impact on substrate selectivity of a structurally conserved loop with a high degree of sequence variability in all classes of sortases. We investigate the contribution of this β7–β8 loop by designing and testing chimeric sortase enzymes. Our chimeras utilize natural sequence variation of class A sortases from eight species engineered into the SrtA sequence from Streptococcus pneumoniae. While some of these chimeric enzymes mimic the activity and selectivity of the WT protein from which the loop sequence was derived (e.g., that of saSrtA), others results in chimeric Streptococcus pneumoniae SrtA enzymes that are able to accommodate a range of residues in the final position of the substrate motif (LPXTX). Using mutagenesis, structural comparisons, and sequence analyses, we identify three interactions facilitated by β7–β8 loop residues that appear to be broadly conserved or converged upon in class A sortase enzymes. These studies provide the foundation for a deeper understanding of sortase target selectivity and can expand the sortase toolbox for future SML applications.     

Construction of a dictionary of sequence motifs that characterize groups of related proteins.

An automatic procedure is proposed to identify, from the protein sequence database, conserved amino acid patterns (or sequence motifs) that are exclusive to a group of functionally related proteins. This procedure is applied to the PIR database and a dictionary of sequence motifs that relate to specific superfamilies constructed. The motifs have a practical relevance in identifying the membership of specific superfamilies without the need to perform sequence database searches in 20% of newly determined sequences. The sequence motifs identified represent functionally important sites on protein molecules. When multiple blocks exist in a single motif they are often close together in the 3-D structure. Furthermore, occasionally these motif blocks were found to be split by introns when the correlation with exon structures was examined.     

Analysis of the substrate specificity loop of the HAD superfamily cap domain

The haloacid dehalogenase (HAD) superfamily includes a variety of enzymes that catalyze the cleavage of substrate C-Cl, P-C, and P-OP bonds via nucleophilic substitution pathways. All members possess the alpha/beta core domain, and many also possess a small cap domain. The active site of the core domain is formed by four loops (corresponding to sequence motifs 1-4), which position substrate and cofactor-binding residues as well as the catalytic groups that mediate the "core" chemistry. The cap domain is responsible for the diversification of chemistry within the family. A tight beta-turn in the helix-loop-helix motif of the cap domain contains a stringently conserved Gly (within sequence motif 5), flanked by residues whose side chains contribute to the catalytic site formed at the domain-domain interface. To define the role of the conserved Gly in the structure and function of the cap domain loop of the HAD superfamily members phosphonoacetaldehyde hydrolase and beta-phosphoglucomutase, the Gly was mutated to Pro, Val, or Ala. The catalytic activity was severely reduced in each mutant. To examine the impact of Gly substitution on loop 5 conformation, the X-ray crystal structure of the Gly50Pro phosphonoacetaldehyde hydrolase mutant was determined. The altered backbone conformation at position 50 had a dramatic effect on the spatial disposition of the side chains of neighboring residues. Lys53, the Schiff Base forming lysine, had rotated out of the catalytic site and the side chain of Leu52 had moved to fill its place. On the basis of these studies, it was concluded that the flexibility afforded by the conserved Gly is critical to the function of loop 5 and that it is a marker by which the cap domain substrate specificity loop can be identified within the amino acid sequence of HAD family members.     

Mass spectrometry and site-directed mutagenesis identify several autophosphorylated residues required for the activity of PrkC,a Ser/Thr kinase from Bacillus subtilis

We have shown recently that PrkC, which is involved in developmental processes in Bacillus subtilis, is a Ser/Thr kinase with features of the receptor kinase family of eukaryotic Hanks kinases. In this study, we expressed and purified from Escherichia coli the cytoplasmic domain of PrkC containing the kinase and a short juxtamembrane region. This fragment, which we designate PrkCc, undergoes autophosphorylation in E.coli. PrkCc is further autophosphorylated in vitro, apparently through a trans-kinase, intermolecular reaction. PrkC also displays kinase activity with myelin basic protein. Using high mass accuracy electrospray tandem mass spectrometry (LC-MS/MS) and nanoelectrospray tandem mass spectrometry, we identified seven phosphorylated threonine and one serine residue in PrkCc. All the corresponding residues were replaced by systematic site-directed mutagenesis and the purified mutant proteins were tested for in vitro kinase activity. Single and multiple replacement of four threonine residues, clustered between residues 162 and 167 in a putative activation loop, substantially reduced kinase activity and the effect was clearly additive. Replacement of the other three threonine residues, clustered between residues 290 and 320, had relatively little effect on activity. In contrast, substitution of Ser214, which is conserved in closely related receptor kinase-like bacterial proteins, independently affected activity and may represent a novel regulatory mechanism. When projected onto a 3D structure of PrkC modelled on the structure of known Hanks kinases, the first cluster of phospho-threonine residues falls precisely in the activation loop, controlling the access of substrate and ATP to the catalytic site of many eukaryotic receptor kinases, whereas the second cluster is located in the juxtamembrane region. These results indicate that regulation of PrkC kinase activity (and presumably autophosphorylation) includes a conserved activation loop mechanism. The juxtamembrane phospho-threonine residues may be essential, for example for the recruitment of other proteins necessary for a PrkC signalling cascade or for coupling to other signalling pathways. This is the first structure-function analysis of a bacterial receptor-like kinase of the Hanks family.     

Comparison of topological clustering within protein networks using edge metrics that evaluate full sequence,full structure,and active site microenvironment similarity

The development of accurate protein function annotation methods has emerged as a major unsolved biological problem. Protein similarity networks, one approach to function annotation via annotation transfer, group proteins into similarity-based clusters. An underlying assumption is that the edge metric used to identify such clusters correlates with functional information. In this contribution, this assumption is evaluated by observing topologies in similarity networks using three different edge metrics: sequence (BLAST), structure (TM-Align), and active site similarity (active site profiling, implemented in DASP). Network topologies for four well-studied protein superfamilies (enolase, peroxiredoxin (Prx), glutathione transferase (GST), and crotonase) were compared with curated functional hierarchies and structure. As expected, network topology differs, depending on edge metric; comparison of topologies provides valuable information on structure/function relationships. Subnetworks based on active site similarity correlate with known functional hierarchies at a single edge threshold more often than sequence- or structure-based networks. Sequence- and structure-based networks are useful for identifying sequence and domain similarities and differences; therefore, it is important to consider the clustering goal before deciding appropriate edge metric. Further, conserved active site residues identified in enolase and GST active site subnetworks correspond with published functionally important residues. Extension of this analysis yields predictions of functionally determinant residues for GST subgroups. These results support the hypothesis that active site similarity-based networks reveal clusters that share functional details and lay the foundation for capturing functionally relevant hierarchies using an approach that is both automatable and can deliver greater precision in function annotation than current similarity-based methods.     

The Short-chain Dehydrogenase/Reductase Engineering Database (SDRED): A classification and analysis system for a highly diverse enzyme family

The Short-chain Dehydrogenases/Reductases Engineering Database (SDRED) covers one of the largest known protein families (168 150 proteins). Assignment to the superfamilies of Classical and Extended SDRs was achieved by global sequence similarity and by identification of family-specific sequence motifs. Two standard numbering schemes were established for Classical and Extended SDRs that allow for the determination of conserved amino acid residues, such as cofactor specificity determining positions or superfamily specific sequence motifs. The comprehensive sequence dataset of the SDRED facilitates the refinement of family-specific sequence motifs. The glycine-rich motifs for Classical and Extended SDRs were refined to improve the precision of superfamily classification. In each superfamily, the majority of sequences formed a tightly connected sequence network and belonged to a large homologous family. Despite their different sequence motifs and their different sequence length, the two sequence networks of Classical and Extended SDRs are not separate, but connected by edges at a threshold of 40% sequence similarity, indicating that all SDRs belong to a large, connected network. The SDRED is accessible at https://sdred.biocatnet.de/.     

Amino acid substitution analysis of E. coli thymidylate synthase: The study of a highly conserved region at the N-terminus

Amino acid substitution analysis within a highly conserved region of Escherichia coli thymidylate synthase (TS), using suppression of amber mutations by tRNA suppressors, has yielded a bank of 124 new mutationally altered TS proteins. These mutant proteins have been used to study the structure-function relationship of the Escherichia coli TS protein at the N-terminus corresponding to residues 20 through 35. This region contains a block of amino acids whose sequence has been well conserved among other known TS proteins from various organisms. Positions 20 through 25 contain a surface loop structure and positions 26 through 35 encompass a β-strand. We find that residues surrounding a β-bulge structure within the β-strand are particularly sensitive to amino acid substitution, suggesting that this structure is maintained by a highly ordered packing arrangement. Three residues in the surface loop that are present at the base of the substrate binding pocket are also sensitive to amino acid substitution. The remainder of the conserved sites, including those at the dimer interface, are tolerant to most, if not all, of the substitutions tested. © 1992 Wiley-Liss, Inc.     

Common active site architecture and binding strategy of four phenylpropanoid P450s from Arabidopsis thaliana as revealed by molecular modeling

Despite extensive primary sequence diversity, crystal structures of several bacterial cytochrome P450 monooxygenases (P450s) and a single eukaryotic P450 indicate that these enzymes share a structural core of alpha-helices and beta-sheets and vary in the loop regions contacting individual substrates. To determine the extent to which individual structural features are conserved among divergent P450s existing in a single biosynthetic pathway, we have modeled the structures of four highly divergent P450s (CYP73A5, CYP84A1, CYP75B1, CYP98A3) in the Arabidopsis phenylpropanoid pathway synthesizing lignins, flavonoids and anthocyanins. Analysis of these models has indicated that, despite primary sequence identities as low as 13%, the structural cores and several loop regions of these P450s are highly conserved. Substrate docking indicated that all four enzymes employ a common strategy to identify their substrates in that their cinnamate-derived substrates align along helix I with their aromatic ring positioned towards the C-terminus of this helix and their aliphatic tails positioned towards the N-terminus. Further similarity was observed in the way the substrates contact the consensus P450 substrate recognition sites (SRS). Residues predicted to contact the aromatic ring region exist in SRS5, SRS6 and the C-terminal portion of SRS4 and residues contacting the distal end of each substrate exist in SRS1, SRS2 and the N-terminal portion of SRS4. Alignments of the regions contacting the aromatic ring region indicate that SRS4, SRS5 and SRS6 share higher degrees of sequence conservation than found in SRS1, SRS2 or the full-length protein.     

FlgM anti-sigma factors: identification of novel members of the family, evolutionary analysis, homology modeling, and analysis of sequence-structure-function relationships

FlgM proteins, also known as Anti-sigma-28 factor (sigma28), are negative regulators of flagellin synthesis. Recently, a three-dimensional structure of the Aquifex aeolicus sigma28/FlgM complex (PDB code: 1rp3) was determined by X-ray crystallography at 2.3 A resolution. Furthermore, experimental data on bacterial FlgM, including site-directed mutagenesis and structural characterization by NMR are also available. However, an interpretation of the sequence-structure-function relationships combining X-ray and NMR data with the evolutionary information extracted from the increasing number of FlgM-related sequences annotated in databases is not available. In the present study, we combined database sequence searches and sequence-analysis tools to update the multiple sequence alignment of a previously characterized cluster of orthologs (COG2747) and the PFAM classification of protein domains (PF04316) for the FlgM family. A phylogenetic analysis of 77 protein sequences revealed the presence of at least three major sequence clades within the FlgM family. Besides, we predicted functional residues using a SequenceSpace method. We also generated homology models for Bacillus subtilis and Salmonella typhimurium FlgM proteins, for which sequence-structure-function relationship data are available, and used the docking program ClusPro to hypothesize about the dimer association between FlgM proteins. In conclusion, the analysis presented in this work will be useful in designing new experiments to understand better protein-protein interactions between FglM, sigma factors, and putative molecules from the flagellar export apparatus. Electronic Supplementary Material is available in the online version of this article at http://link.springer.de/     

NMR structure of the Bordetella bronchiseptica protein NP_888769.1 establishes a new phage-related protein family PF13554

The solution structure of the hypothetical phage-related protein NP_888769.1 from the Gram-negative bacterium Bordetella bronchoseptica contains a well-structured core comprising a five-stranded, antiparallel β-sheet packed on one side against two α-helices and a short β-hairpin with three flexibly disordered loops extending from the central β-sheet. A homology search with the software DALI identified two Protein Data Bank deposits with Z-scores > 8, where both of these proteins have less than 8% sequence identity relative to NP_888769.1, and one has been functionally annotated as a lambda phage tail terminator protein. A sequence-homology analysis then confirmed that NP_888769.1 represents the first three-dimensional structural representative of a new protein family that was previously predicted by the Joint Center for Structural Genomics, which includes so far about 20 prophage proteins encoded in bacterial genomes.     

Triggering loops and enzyme function: identification of loops that trigger and modulate movements

Enzyme function often involves a conformational change. There is a general agreement that loops play a vital role in correctly positioning the catalytically important residues. Nevertheless, predicting the functional loops and most importantly their role in enzyme function remains a difficult task. A major reason for this difficulty is that loops that undergo conformational change are frequently not well conserved in their primary sequence. beta1,4-Galactosyltransferase is one such enzyme. There, the amino acid sequence of a long loop that undergoes a large conformational change upon substrate binding is not well conserved. Our molecular dynamics simulations show that the large conformational change in the long loop is brought about by a second, interacting loop. Interestingly, while the structural change of the second loop is much smaller than that of the long loop, its sequence (particularly glycine residues) is highly conserved. We further examine the generality of the proposition that there are loops that trigger movements but nevertheless show little or no structural changes in crystals. We focus on two other enzymes, enolase and lipase. We chose these enzymes, since they too undergo conformational change upon ligand binding, however, they have different folds and different functions. Through multiple sets of simulations we show that the conformational change of the functional loop(s) is brought about through communication of flexibility by triggering loops that have several glycine residues. We further propose that similar to the conservation of common favorable fold types and structural motifs, evolution has also conserved common "skillful" mechanisms. Mechanisms may be conserved across different folds, sequences and functions, with adaptation to specific enzymatic roles.     

Superfamily active site templates

We show that three-dimensional signatures consisting of only a few functionally important residues can be diagnostic of membership in superfamilies of enzymes. Using the enolase superfamily as a model system, we demonstrate that such a signature, or template, can identify superfamily members in structural databases with high sensitivity and specificity. This is remarkable because superfamilies can be highly diverse, with members catalyzing many different overall reactions; the unifying principle can be a conserved partial reaction or chemical capability. Our definition of a superfamily thus hinges on the disposition of residues involved in a conserved function, rather than on fold similarity alone. A clear advantage of basing structure searches on such active site templates rather than on fold similarity is the specificity with which superfamilies with distinct functional characteristics can be identified within a large set of proteins with the same fold, such as the (beta/alpha)8 barrels. Preliminary results are presented for an additional group of enzymes with a different fold, the haloacid dehalogenase superfamily, suggesting that this approach may be generally useful for assigning reading frames of unknown function to specific superfamilies and thereby allowing inference of some of their functional properties.     

A robust method to detect structural and functional remote homologues

With currently available sequence data, it is feasible to conduct extensive comparisons among large sets of protein sequences. It is still a much more challenging task to partition the protein space into structurally and functionally related families solely based on sequence comparisons. The ProtoNet system automatically generates a treelike classification of the whole protein space. It stands to reason that this classification reflects evolutionary relationships, both close and remote. In this article, we examine this hypothesis. We present a semiautomatic procedure that singles out certain inner nodes in the ProtoNet tree that should ideally correspond to structurally and functionally defined protein families. We compare the performance of this method against several expert systems. Some of the competing methods incorporate additional extraneous information on protein structure or on enzymatic activities. The ProtoNet-based method performs at least as well as any of the methods with which it was compared. This article illustrates the ProtoNet-based method on several evolutionarily diverse families. Using this new method, an evolutionary divergence scheme can be proposed for a large number of structural and functional related superfamilies.     

Unique Features of Human Protein Arginine Methyltransferase 9 (PRMT9) and Its Substrate RNA Splicing Factor SF3B2

Human protein arginine methyltransferase (PRMT) 9 symmetrically dimethylates arginine residues on splicing factor SF3B2 (SAP145) and has been functionally linked to the regulation of alternative splicing of pre-mRNA. Site-directed mutagenesis studies on this enzyme and its substrate had revealed essential unique residues in the double E loop and the importance of the C-terminal duplicated methyltransferase domain. In contrast to what had been observed with other PRMTs and their physiological substrates, a peptide containing the methylatable Arg-508 of SF3B2 was not recognized by PRMT9 in vitro. Although amino acid substitutions of residues surrounding Arg-508 had no great effect on PRMT9 recognition of SF3B2, moving the arginine residue within this sequence abolished methylation. PRMT9 and PRMT5 are the only known mammalian enzymes capable of forming symmetric dimethylarginine (SDMA) residues as type II PRMTs. We demonstrate here that the specificity of these enzymes for their substrates is distinct and not redundant. The loss of PRMT5 activity in mouse embryo fibroblasts results in almost complete loss of SDMA, suggesting that PRMT5 is the primary SDMA-forming enzyme in these cells. PRMT9, with its duplicated methyltransferase domain and conserved sequence in the double E loop, appears to have a unique structure and specificity among PRMTs for methylating SF3B2 and potentially other polypeptides.