ESAT-6 and Ag85A Synthetic Peptides as Candidates for an Immunodiagnostic Test in Children with a Clinical Suspicion of Tuberculosis.

BACKGROUND: Tuberculosis (TB) is being underdetected in children as most are smear-negative. This work was aimed at evaluating ESAT-6 and Ag85A synthetic peptides' serodiagnostic potential for diagnosing children having a clinical suspicion of TB. METHODS: The study involved 438 children: 77 Creole nonindigenous (13 suspected of having TB and 64 healthy ones) and 361 Warao indigenous children (39 suspected of TB and 322 healthy children). The approach's diagnostic information was compared using operational characteristics and receiver-operating characteristic (ROC) curves. RESULTS: Ag85A P-29879 had 94.6% sensitivity (AUC = 0.741: 0.651 to 0.819 95% CI) in indigenous children. ESAT-6 P-12036 and P-12037 had 100% and 92.3% of sensitivity (AUC = 0.929: 0.929: 0.846 to 0.975 95% CI and 0.791: 63.9 to 98.7 95% CI, respectively) in Creole children. ESAT-6 peptides also allowed a differentiation between children with TB and healthy ones. CONCLUSIONS: Further validation of this approach could lead to developing a complementary tool for rapid TB diagnosis in children.

Virtual screening to identify Leishmania braziliensis N-myristoyltransferase inhibitors: pharmacophore models, docking, and molecular dynamics.

Leishmaniasis is caused by several protozoa species belonging to genus Leishmania that are hosted by humans and other mammals. Millions of new cases are recorded every year and the drugs available on the market do not show satisfactory efficacy and safety. A hierarchical virtual screening approach based on the pharmacophore model, molecular docking, and molecular dynamics was conducted to identify possible Leishmania braziliensis N-misristoyltransferase (LbNMT) inhibitors. The adopted pharmacophore model had three main features: four hydrophobic centers, four hydrogen-bond acceptor atoms, and one positive nitrogen center. The molecules (n=15,000) were submitted to alignment with the pharmacophore model and only 27 molecules aligned to model. Six molecules were submitted to molecular docking, using receptor PDB ID 5A27. After docking, the ZINC35426134 was a top-ranked molecule (- 64.61 kcal/mol). The molecule ZINC35426134 shows hydrophobic interactions with Phe82, Tyr209, Val370, and Leu391 and hydrogen bonds with Asn159, Tyr318, and Val370. Molecular dynamics simulations were performed with the protein in its APO and HOLO forms for 37 ns in order to assess the stability of the protein-ligand complex. Results showed that the HOLO form was more stable than the APO one, and it suggests that the ZINC35426134 binding stabilizes the enzyme. Therefore, the selected molecule has the potential to meet the herein proposed target.

Phasins, Multifaceted Polyhydroxyalkanoate Granule-Associated Proteins.

Phasins are the major polyhydroxyalkanoate (PHA) granule-associated proteins. They promote bacterial growth and PHA synthesis and affect the number, size, and distribution of the granules. These proteins can be classified in 4 families with distinctive characteristics. Low-resolution structural studies and in silico predictions were performed in order to elucidate the structure of different phasins. Most of these proteins share some common structural features, such as a preponderant α-helix composition, the presence of disordered regions that provide flexibility to the protein, and coiled-coil interacting regions that form oligomerization domains. Due to their amphiphilic nature, these proteins play an important structural function, forming an interphase between the hydrophobic content of PHA granules and the hydrophilic cytoplasm content. Phasins have been observed to affect both PHA accumulation and utilization. Apart from their role as granule structural proteins, phasins have a remarkable variety of additional functions. Different phasins have been determined to (i) activate PHA depolymerization, (ii) increase the expression and activity of PHA synthases, (iii) participate in PHA granule segregation, and (iv) have both in vivo and in vitro chaperone activities. These properties suggest that phasins might play an active role in PHA-related stress protection and fitness enhancement. Due to their granule binding capacity and structural flexibility, several biotechnological applications have been developed using different phasins, increasing the interest in the study of these remarkable proteins.

Structural and Affinity Determinants in the Interaction between Alcohol Acyltransferase from F. x ananassa and Several Alcohol Substrates: A Computational Study.

Aroma and flavor are important factors of fruit quality and consumer preference. The specific pattern of aroma is generated during ripening by the accumulation of volatiles compounds, which are mainly esters. Alcohol acyltransferase (AAT) (EC 2.3.1.84) catalyzes the esterification reaction of aliphatic and aromatic alcohols and acyl-CoA into esters in fruits and flowers. In Fragaria x ananassa, there are different volatiles compounds that are obtained from different alcohol precursors, where octanol and hexanol are the most abundant during fruit ripening. At present, there is not structural evidence about the mechanism used by the AAT to synthesize esters. Experimental data attribute the kinetic role of this enzyme to 2 amino acidic residues in a highly conserved motif (HXXXD) that is located in the middle of the protein. With the aim to understand the molecular and energetic aspects of volatiles compound production from F. x ananassa, we first studied the binding modes of a series of alcohols, and also different acyl-CoA substrates, in a molecular model of alcohol acyltransferase from Fragaria x ananassa (SAAT) using molecular docking. Afterwards, the dynamical behavior of both substrates, docked within the SAAT binding site, was studied using routine molecular dynamics (MD) simulations. In addition, in order to correlate the experimental and theoretical data obtained in our laboratories, binding free energy calculations were performed; which previous results suggested that octanol, followed by hexanol, presented the best affinity for SAAT. Finally, and concerning the SAAT molecular reaction mechanism, it is suggested from molecular dynamics simulations that the reaction mechanism may proceed through the formation of a ternary complex, in where the Histidine residue at the HXXXD motif deprotonates the alcohol substrates. Then, a nucleophilic attack occurs from alcohol charged oxygen atom to the carbon atom at carbonyl group of the acyl CoA. This mechanism is in agreement with previous results, obtained in our group, in alcohol acyltransferase from Vasconcellea pubescens (VpAAT1).

The canonical DHHC motif is not absolutely required for the activity of the yeast S-acyltransferases Swf1 and Pfa4.

Protein S-acyltransferases, also known as palmitoyltransferases (PATs), are characterized by the presence of a 50-amino acid domain called the DHHC domain. Within this domain, these four amino acids constitute a highly conserved motif. It has been proposed that the palmitoylation reaction occurs through a palmitoyl-PAT covalent intermediate that involves the conserved cysteine in the DHHC motif. Mutation of this cysteine results in lack of function for several PATs, and DHHA or DHHS mutants are used regularly as catalytically inactive controls. In a genetic screen to isolate loss-of-function mutations in the yeast PAT Swf1, we isolated an allele encoding a Swf1 DHHR mutant. Overexpression of this mutant is able to partially complement a swf1Δ strain and to acylate the Swf1 substrates Tlg1, Syn8, and Snc1. Overexpression of the palmitoyltransferase Pfa4 DHHA or DHHR mutants also results in palmitoylation of its substrate Chs3. We also investigated the role of the first histidine of the DHHC motif. A Swf1 DQHC mutant is also partially active but a DQHR is not. Finally, we show that Swf1 substrates are differentially modified by both DHHR and DQHC Swf1 mutants. We propose that, in the absence of the canonical mechanism, alternative suboptimal mechanisms take place that are more dependent on the reactivity of the acceptor protein. These results also imply that caution must be exercised when proposing non-canonical roles for PATs on the basis of considering DHHC mutants as catalytically inactive and, more generally, contribute to an understanding of the mechanism of protein palmitoylation.

Identification of Giardia lamblia DHHC proteins and the role of protein S-palmitoylation in the encystation process.

Protein S-palmitoylation, a hydrophobic post-translational modification, is performed by protein acyltransferases that have a common DHHC Cys-rich domain (DHHC proteins), and provides a regulatory switch for protein membrane association. In this work, we analyzed the presence of DHHC proteins in the protozoa parasite Giardia lamblia and the function of the reversible S-palmitoylation of proteins during parasite differentiation into cyst. Two specific events were observed: encysting cells displayed a larger amount of palmitoylated proteins, and parasites treated with palmitoylation inhibitors produced a reduced number of mature cysts. With bioinformatics tools, we found nine DHHC proteins, potential protein acyltransferases, in the Giardia proteome. These proteins displayed a conserved structure when compared to different organisms and are distributed in different monophyletic clades. Although all Giardia DHHC proteins were found to be present in trophozoites and encysting cells, these proteins showed a different intracellular localization in trophozoites and seemed to be differently involved in the encystation process when they were overexpressed. dhhc transgenic parasites showed a different pattern of cyst wall protein expression and yielded different amounts of mature cysts when they were induced to encyst. Our findings disclosed some important issues regarding the role of DHHC proteins and palmitoylation during Giardia encystation.

Divergent metabolic phenotype between two sisters with congenital generalized lipodystrophy due to double AGPAT2 homozygous mutations. a clinical, genetic and in silico study.

Congenital generalized lipodystrophy (CGL) is a rare autosomal recessive disorder characterized by extreme reduction of white adipose tissue (WAT) mass. CGL type 1 is the most frequent form and is caused by mutations in AGPAT2. Genetic and clinical studies were performed in two affected sisters of a Chilean family. These patients have notoriously dissimilar metabolic abnormalities that correlate with differential levels of circulating leptin and soluble leptin receptor fraction. Sequencing of AGPAT2 exons and exon-intron boundaries revealed two homozygous mutations in both sisters. Missense mutation c.299G>A changes a conserved serine in the acyltransferase NHX4D motif of AGPAT2 (p.Ser100Asn). Intronic c.493-1G>C mutation destroy a conserved splicing site that likely leads to exon 4 skipping and deletion of whole AGPAT2 substrate binding domain. In silico protein modeling provided insights of the mechanisms of lack of catalytic activity owing to both mutations.

Zinc co-ordination by the DHHC cysteine-rich domain of the palmitoyltransferase Swf1.

S-acylation, commonly known as palmitoylation, is a widespread post-translational modification of proteins that consists of the thioesterification of one or more cysteine residues with fatty acids. This modification is catalysed by a family of PATs (palmitoyltransferases), characterized by the presence of a 50-residue long DHHC-CRD (Asp-His-His-Cys cysteine-rich domain). To gain knowledge on the structure-function relationships of these proteins, we carried out a random-mutagenesis assay designed to uncover essential amino acids in Swf1, the yeast PAT responsible for the palmitoylation of SNARE (soluble N-ethylmaleimide-sensitive fusion protein-attachment protein receptor) proteins. We identified 21 novel loss-of-function mutations, which are mostly localized within the DHHC-CRD. Modelling of the tertiary structure of the Swf1 DHHC domain suggests that it could fold as a zinc-finger domain, co-ordinating two zinc atoms in a CCHC arrangement. All residues predicted to be involved in the co-ordination of zinc were found to be essential for Swf1 function in the screen. Moreover, these mutations result in unstable proteins, in agreement with a structural role for these zinc fingers. The conservation of amino acids predicted to form each zinc-binding pocket suggests a shared function, as the selective pressure to maintain them is lost upon mutation of one of them. A Swf1 orthologue that lacks one of the zinc-binding pockets is able to complement a yeast swf1∆ strain, possibly because a similar fold can be stabilized by hydrogen bonds instead of zinc co-ordination. Finally, we show directly that recombinant Swf1 DHHC-CRD is able to bind zinc. Sequence analyses of DHHC domains allowed us to present models of the zinc-binding properties for all PATs.

Identification and subcellular localization of TcHIP, a putative Golgi zDHHC palmitoyl transferase of Trypanosoma cruzi.

Protein palmitoylation is a post-translational modification that contributes to determining protein localization and function. Palmitoylation has been described in trypanosomatid protozoa, but no zDHHC palmitoyl transferase has been identified in Trypanosoma cruzi, the etiological agent of Chagas disease in Latin America. In this study we identify and show the subcellular localization of TcHIP (Tc00.1047053508199.50), a putative T. cruzi zDHHC palmitoyl transferase. Analysis of the deduced protein sequence indicates that it contains ankyrin repeats (Ank and Ank2) and the zDHHC conserved domain, typical of zDHHC palmitoyl transferases. A TcHIP polyclonal antiserum obtained from mice immunized with the purified recombinant protein was used to study the presence and subcellular localization of the native enzyme. In western blots this antiserum recognized a protein of about 95 kDa, consistent with the predicted molecular mass of TcHIP (95.4 kDa), in whole extracts of T. cruzi epimastigotes, metacyclic trypomastigotes and intracellular amastigotes. Immunolocalization by confocal microscopy showed TcHIP labeling at the Golgi complex, co-localizing with the T. cruzi Golgi marker TcRab7-GFP. Transfectant T. cruzi epimastigotes containing a construct encoding TcHIP fused to proteins A and C (TcHIP/AC) were obtained. In western blotting experiments, the TcHIP polyclonal antiserum recognized both native and TcHIP/AC proteins in extracts of the transfectants. Confocal microscopy showed co-localization of native TcHIP with TcHIP/AC. These findings demonstrate the presence of a putative zDHHC palmitoyl transferase (TcHIP) containing ankyrin and zDHHC domains in different developmental forms of T. cruzi, and its association with the Golgi complex.

Striking natural diversity in glandular trichome acylsugar composition is shaped by variation at the Acyltransferase2 locus in the wild tomato Solanum habrochaites.

Acylsugars are polyesters of short- to medium-length acyl chains on sucrose or glucose backbones that are produced in secretory glandular trichomes of many solanaceous plants, including cultivated tomato (Solanum lycopersicum). Despite their roles in biotic stress adaptation and their wide taxonomic distribution, there is relatively little information about the diversity of these compounds and the genes responsible for their biosynthesis. In this study, acylsugar diversity was assessed for 80 accessions of the wild tomato species Solanum habrochaites from throughout the Andes Mountains. Trichome metabolites were analyzed by liquid chromatography-time of flight-mass spectrometry, revealing the presence of at least 34 structurally diverse acylsucroses and two acylglucoses. Distinct phenotypic classes were discovered that varied based on the presence of glucose or sucrose, the numbers and lengths of acyl chains, and the relative total amounts of acylsugars. The presence or absence of an acetyl chain on the acylsucrose hexose ring caused clustering of the accessions into two main groups. Analysis of the Acyltransferase2 gene (the apparent ortholog of Solyc01g105580) revealed differences in enzyme activity and gene expression correlated with polymorphism in S. habrochaites accessions that varied in acylsucrose acetylation. These results are consistent with the hypothesis that glandular trichome acylsugar acetylation is under selective pressure in some populations of S. habrochaites and that the gene mutates to inactivity in the absence of selection.

Homology modeling of the structure of acyl coA:isopenicillin N-acyltransferase (IAT) from Penicillium chrysogenum. IAT interaction studies with isopenicillin-N, combining molecular dynamics simulations and docking.

In the last step of penicillin biosynthesis, acyl-CoA:isopenicillin N acyltransferase (IAT) (E.C. 2.3.1.164) catalyzes the conversion of isopenicillin N (IPN) to penicillin G. IAT substitutes the α-aminoadipic acid side chain of IPN by a phenylacetic acid phenolate group (from phenylacetyl-CoA). Having a three-dimensional (3D) structure of IAT helps to determine the steps involved in side chain exchange by identifying the atomic details of substrate recognition. We predicted the IAT 3-D structure (α- and ß-subunits), as well as the manner of IPN and phenylacetyl-CoA bind to the mature enzyme (ß-subunit). The 3D IAT prediction was achieved by homology modeling and molecular docking in different snapshots, and refined by molecular dynamic simulations. Our model can reasonably interpret the results of a number of experiments, where key residues for IAT processing as well as strictly conserved residues most probably involved with enzymatic activity were mutated. Based on the results of docking studies, energies associated with the complexes, and binding constants calculated, we identified a site located in the region generated by ß1, ß2 and ß5 strands, which forms part of the central structure of ß-subunit, as the potential binding site of IPN. The site comprises the amino acid residues Cys103, Asp121, Phe122, Phe123, Ala168, Leu169, His170, Gln172, Phe212, Arg241, Leu262, Asp264, Arg302, Ser309, and Arg310. Through hydrogen bonds, the IPN binding site establishes interactions with Cys103, Leu169, Gln172, Asp264 and Arg310. Our model is also validated by a recently revealed crystal structure of the mature enzyme.

Specificity of transmembrane protein palmitoylation in yeast.

Many proteins are modified after their synthesis, by the addition of a lipid molecule to one or more cysteine residues, through a thioester bond. This modification is called S-acylation, and more commonly palmitoylation. This reaction is carried out by a family of enzymes, called palmitoyltransferases (PATs), characterized by the presence of a conserved 50- aminoacids domain called "Asp-His-His-Cys- Cysteine Rich Domain" (DHHC-CRD). There are 7 members of this family in the yeast Saccharomyces cerevisiae, and each of these proteins is thought to be responsible for the palmitoylation of a subset of substrates. Substrate specificity of PATs, however, is not yet fully understood. Several yeast PATs seem to have overlapping specificity, and it has been proposed that the machinery responsible for palmitoylating peripheral membrane proteins in mammalian cells, lacks specificity altogether.Here we investigate the specificity of transmembrane protein palmitoylation in S. cerevisiae, which is carried out predominantly by two PATs, Swf1 and Pfa4. We show that palmitoylation of transmembrane substrates requires dedicated PATs, since other yeast PATs are mostly unable to perform Swf1 or Pfa4 functions, even when overexpressed. Furthermore, we find that Swf1 is highly specific for its substrates, as it is unable to substitute for other PATs. To identify where Swf1 specificity lies, we carried out a bioinformatics survey to identify amino acids responsible for the determination of specificity or Specificity Determination Positions (SDPs) and showed experimentally, that mutation of the two best SDP candidates, A145 and K148, results in complete and partial loss of function, respectively. These residues are located within the conserved catalytic DHHC domain suggesting that it could also be involved in the determination of specificity. Finally, we show that modifying the position of the cysteines in Tlg1, a Swf1 substrate, results in lack of palmitoylation, as expected for a highly specific enzymatic reaction.

Structural characterization and substrate specificity of VpAAT1 protein related to ester biosynthesis in mountain papaya fruit.

The aroma in fruits is an important attribute of quality that influences consumer's acceptance. This attribute is a complex character determined by a set of low molecular weight volatile compounds. In mountain papaya fruit (Vasconcellea pubescens) the aroma is determined mainly by esters, which are produced through an esterification reaction catalyzed by the enzyme alcohol acyltransferase (AAT) that utilizes alcohols and acyl-CoAs as substrates. In order to understand the molecular mechanism involved in the production of esters in this fruit, an AAT gene which has been previously cloned and characterized from mountain papaya (VpAAT1) was expressed in yeasts, and the highest enzyme activity of the recombinant protein was obtained when the enzyme was tested for its ability to produce benzyl acetate. On the other hand, to gain insight the mechanism of action at the molecular level, a structural model for VpAAT1 protein was built by comparative modelling methodology, which was validated and refined by molecular dynamics simulation. The VpAAT1 structure consists of two domains connected by a large crossover loop, with a solvent channel in the center of the structure formed between the two domains. Residues H166 and D170, important for catalytic action, displayed their side chains towards the central cavity of the channel allowing their interaction with the substrates. The conformational interaction between the protein and several ligands was explored by molecular docking simulations, and the predictions obtained were tested through kinetic analysis. Kinetic results showed that the lowest K(M) values were obtained for acetyl-CoA and benzyl alcohol. In addition, the most favorable predicted substrate orientation was observed for benzyl alcohol and acetyl CoA, showing a perfect coincidence between kinetic studies and molecular docking analysis.

Bi- and multilinear PLS coupled to MIA-QSAR in the prediction of antifungal activities of some benzothiazole derivatives.

The activities of a series of benzothiazole derivatives, some Candida albicans N-myristoyltransferase (Nmt) inhibitors, were modeled through MIA-QSAR (multivariate image analysis applied to quantitative structure-activity relationship) by using two different regression methods: N-PLS, applied to the three-way array, and PLS, applied to the unfolded array. Both models demonstrated excellent predictive ability, with results comparable to those obtained through 3D approaches. In order to compare the results obtained through MIA descriptors with the predictions of a classical 2D QSAR, some representative physicochemical descriptors were calculated and regressed against the experimental pIC50 values through multiple linear regression, demonstrating that MIA-QSAR was superior for this series of compounds.

A novel motif at the C-terminus of palmitoyltransferases is essential for Swf1 and Pfa3 function in vivo.

S-acylation (commonly known as palmitoylation) is a widespread post-translational modification that consists of the addition of a lipid molecule to cysteine residues of a protein through a thioester bond. This modification is predominantly mediated by a family of proteins referred to as PATs (palmitoyltransferases). Most PATs are polytopic membrane proteins, with four to six transmembrane domains, a conserved DHHC motif and variable C-and N-terminal regions, that are probably responsible for conferring localization and substrate specificity. There is very little additional information on the structure-function relationship of PATs. Swf1 and Pfa3 are yeast members of the DHHC family of proteins. Swf1 is responsible for the S-acylation of several transmembrane SNAREs (soluble N-ethylmaleimide-sensitive fusion protein-attachment protein receptors) and other integral membrane proteins. Pfa3 is required for the palmitoylation of Vac8, a protein involved in vacuolar fusion. In the present study we describe a novel 16-amino-acid motif present at the cytosolic C-terminus of PATs, that is required for Swf1 and Pfa3 function in vivo. Within this motif, we have identified a single residue in Swf1, Tyr323, as essential for function, and this is correlated with lack of palmitoylation of Tlg1, a SNARE that is a substrate of Swf1. The equivalent mutation in Pfa3 also affects its function. These mutations are the first phenotype-affecting mutations uncovered that do not lie within the DHHC domain, for these or any other PATs. The motif is conserved in 70% of PATs from all eukaryotic organisms analysed, and may have once been present in all PATs. We have named this motif PaCCT ('Palmitoyltransferase Conserved C-Terminus').

The structurally constrained protein evolution model accounts for sequence patterns of the LbetaH superfamily.

BACKGROUND: Structure conservation constrains evolutionary sequence divergence, resulting in observable sequence patterns. Most current models of protein evolution do not take structure into account explicitly, being unsuitable for investigating the effects of structure conservation on sequence divergence. To this end, we recently developed the Structurally Constrained Protein Evolution (SCPE) model. The model starts with the coding sequence of a protein with known three-dimensional structure. At each evolutionary time-step of an SCPE simulation, a trial sequence is generated by introducing a random point mutation in the current coding DNA sequence. Then, a "score" for the trial sequence is calculated and the mutation is accepted only if its score is under a given cutoff, lambda. The SCPE score measures the distance between the trial sequence and a given reference sequence, given the structure. In our first brief report we used a "global score", in which the same reference sequence, the ancestral one, was used at each evolutionary step. Here, we introduce a new scoring function, the "local score", in which the sequence accepted at the previous evolutionary time-step is used as the reference. We assess the model on the UDP-N-acetylglucosamine acyltransferase (LPXA) family, as in our previous report, and we extend this study to all other members of the left-handed parallel beta helix fold (LbetaH) superfamily whose structure has been determined. RESULTS: We studied site-dependent entropies, amino acid probability distributions, and substitution matrices predicted by SCPE and compared with experimental data for several members of the LbetaH superfamily. We also evaluated structure conservation during simulations. Overall, SCPE outperforms JTT in the description of sequence patterns observed in structurally constrained sites. Maximum Likelihood calculations show that the local-score and global-score SCPE substitution matrices obtained for LPXA outperform the JTT model for the LPXA family and for the structurally constrained sites of class i of other members within the LbetaH superfamily. CONCLUSION: We extended the SCPE model by introducing a new scoring function, the local score. We performed a thorough assessment of the SCPE model on the LPXA family and extended it to all other members of known structure of the LbetaH superfamily.

Structural constraints and emergence of sequence patterns in protein evolution.

The aim of this work was to study the relationship between structure conservation and sequence divergence in protein evolution. To this end, we developed a model of structurally constrained protein evolution (SCPE) in which trial sequences, generated by random mutations at gene level, are selected against departure from a reference three-dimensional structure. Since at the mutational level SCPE is completely unbiased, any emergent sequence pattern will be due exclusively to structural constraints. In this first report, it is shown that SCPE correctly predicts the characteristic hexapeptide motif of the left-handed parallel beta helix (LbetaH) domain of UDP-N-acetylglucosamine acyltransferases (LpxA).

Evolutionary analysis of gamma-carbonic anhydrase and structurally related proteins.

We studied the evolutionary relationships between gamma-carbonic anhydrase (gamma-CA) and a very diverse group of proteins that share the sequence motif characteristic of the left-handed parallel beta-helix (LbetaH) fold. This sequence motif is characterized by the imperfect tandem repetition of short hexapeptide units, which makes it difficult to obtain a reliable alignment based on sequence information alone. To solve this problem, we used a structural alignment of three members of the group with known crystallographic structures as a seed to obtain a reliable sequence alignment. Then, we applied protein maximum-parsimony and maximum-likelihood phylogenetic inference methods to this alignment. We found that gamma-CA belongs to a diverse superfamily of proteins that share the LbetaH domain. This superfamily is composed mainly of acyltransferases. The most remarkable feature of the phylogenetic tree obtained is that its main branches group together functionally related proteins, so that the coarse topology can be rather easily explained in terms of functional diversification. Regarding the main branch of the tree containing gamma-CA, we found that, in addition to the group of its closest relatives that had already been studied, gamma-CA is closely related to the tetrahydrodipicolinate N-succinyltransferases.