Algoritmos evolutivos e modelo HP para predição de estruturas de proteínas

Predição de estruturas de proteínas (PSP) é um problema computacionalmente complexo. Modelos simplificados da molécula proteica (como o Modelo HP) e o uso de Algoritmos Evolutivos (AEs) estão entre as principais técnicas investigadas para PSP. Entretanto, a avaliação de uma estrutura representada pelo Modelo HP considera apenas o número de contatos hidrofóbicos, não possibilitando distinguir entre estruturas com o mesmo número de contatos hidrofóbicos. Neste trabalho, é apresentada uma nova formulação multiobjetivo para PSP em Modelo HP. Duas métricas são avaliadas: o número de contatos hidrofóbicos e a distância entre os aminoácidos hidrofóbicos, as quais são tratados pelo AE Multiobjetivo em Tabelas (AEMT). O algoritmo mostrou-se rápido e robusto.

Relationship between global structural parameters and Enzyme Commission hierarchy: Implications for function prediction

In protein databases there is a substantial number of proteins structurally determined but without function annotation. Understanding the relationship between function and structure can be useful to predict function on a large scale. We have analyzed the similarities in global physicochemical parameters for a set of enzymes which were classified according to the four Enzyme Commission (EC) hierarchical levels. Using relevance theory we introduced a distance between proteins in the space of physicochemical characteristics. This was done by minimizing a cost function of the metric tensor built to reflect the EC classification system. Using an unsupervised clustering method on a set of 1025 enzymes, we obtained no relevant clustering formation compatible with EC classification. The distance distributions between enzymes from the same EC group and from different EC groups were compared by histograms. Such analysis was also performed using sequence alignment similarity as a distance. Our results suggest that global structure parameters are not sufficient to segregate enzymes according to EC hierarchy. This indicates that features essential for function are rather local than global. Consequently, methods for predicting function based on global attributes should not obtain high accuracy in main EC classes prediction without relying on similarities between enzymes from training and validation datasets. Furthermore, these results are consistent with a substantial number of studies suggesting that function evolves fundamentally by recruitment, i.e., a same protein motif or fold can be used to perform different enzymatic functions and a few specific amino acids (AAs) are actually responsible for enzyme activity. These essential amino acids should belong to active sites and an effective method for predicting function should be able to recognize them. (C) 2012 Elsevier Ltd. All rights reserved.

Extensive structural change of the envelope protein of dengue virus induced by a tuned ionic strength: conformational and energetic analyses

The Dengue has become a global public health threat, with over 100 million infections annually; to date there is no specific vaccine or any antiviral drug. The structures of the envelope (E) proteins of the four known serotype of the dengue virus (DENV) are already known, but there are insufficient molecular details of their structural behavior in solution in the distinct environmental conditions in which the DENVs are submitted, from the digestive tract of the mosquito up to its replication inside the host cell. Such detailed knowledge becomes important because of the multifunctional character of the E protein: it mediates the early events in cell entry, via receptor endocytosis and, as a class II protein, participates determinately in the process of membrane fusion. The proposed infection mechanism asserts that once in the endosome, at low pH, the E homodimers dissociate and insert into the endosomal lipid membrane, after an extensive conformational change, mainly on the relative arrangement of its three domains. In this work we employ all-atom explicit solvent Molecular Dynamics simulations to specify the thermodynamic conditions in that the E proteins are induced to experience extensive structural changes, such as during the process of reducing pH. We study the structural behavior of the E protein monomer at acid pH solution of distinct ionic strength. Extensive simulations are carried out with all the histidine residues in its full protonated form at four distinct ionic strengths. The results are analyzed in detail from structural and energetic perspectives, and the virtual protein movements are described by means of the principal component analyses. As the main result, we found that at acid pH and physiological ionic strength, the E protein suffers a major structural change; for lower or higher ionic strengths, the crystal structure is essentially maintained along of all extensive simulations. On the other hand, at basic pH, when all histidine residues are in the unprotonated form, the protein structure is very stable for ionic strengths ranging from 0 to 225 mM. Therefore, our findings support the hypothesis that the histidines constitute the hot points that induce configurational changes of E protein in acid pH, and give extra motivation to the development of new ideas for antivirus compound design.

On the impact of disproportional samples in credit scoring models: An application to a Brazilian bank data

Statistical methods have been widely employed to assess the capabilities of credit scoring classification models in order to reduce the risk of wrong decisions when granting credit facilities to clients. The predictive quality of a classification model can be evaluated based on measures such as sensitivity, specificity, predictive values, accuracy, correlation coefficients and information theoretical measures, such as relative entropy and mutual information. In this paper we analyze the performance of a naive logistic regression model (Hosmer & Lemeshow, 1989) and a logistic regression with state-dependent sample selection model (Cramer, 2004) applied to simulated data. Also, as a case study, the methodology is illustrated on a data set extracted from a Brazilian bank portfolio. Our simulation results so far revealed that there is no statistically significant difference in terms of predictive capacity between the naive logistic regression models and the logistic regression with state-dependent sample selection models. However, there is strong difference between the distributions of the estimated default probabilities from these two statistical modeling techniques, with the naive logistic regression models always underestimating such probabilities, particularly in the presence of balanced samples. (C) 2012 Elsevier Ltd. All rights reserved.

Peptide fingerprinting of the neurotoxic fractions isolated from the secretions of sea anemones Stichodactyla helianthus and Bunodosoma granulifera. New members of the APETx-like family identified by a 454 pyrosequencing approach

Sea anemones are known to contain a wide diversity of biologically active peptides, mostly unexplored according to recent peptidomic and transcriptomic studies. In the present work, the neurotoxic fractions from the exudates of Stichodactyla helianthus and Bunodosoma granulifera were analyzed by reversed-phase chromatography and mass spectrometry. The first peptide fingerprints of these sea anemones were assessed, revealing the largest number of peptide components (156) so far found in sea anemone species, as well as the richer peptide diversity of B. granulifera in relation to S. helianthus. The transcriptomic analysis of B. granulifera, performed by massive cDNA sequencing with 454 pyrosequencing approach allowed the discovery of five new APETx-like peptides (U-AITX-Bg1a-e - including the full sequences of their precursors for four of them), which together with type 1 sea anemone sodium channel toxins constitute a very distinguishable feature of studied sea anemone species belonging to genus Bunodosoma. The molecular modeling of these new APETx-like peptides showed a distribution of positively charged and aromatic residues in putative contact surfaces as observed in other animal toxins. On the other hand, they also showed variable electrostatic potentials, thus suggesting a docking onto their targeted channels in different spatial orientations. Moreover several crab paralyzing toxins (other than U-AITX-Bg1a-e), which induce a variety of symptoms in crabs, were isolated. Some of them presumably belong to new classes of crab-paralyzing peptide toxins, especially those with molecular masses below 2 kDa, which represent the smallest peptide toxins found in sea anemones. (C) 2011 Elsevier Inc. All rights reserved.

The Diversity of Polyketide Synthase Genes from Sugarcane-Derived Fungi

The chemical ecology and biotechnological potential of metabolites from endophytic and rhizosphere fungi are receiving much attention. A collection of 17 sugarcane-derived fungi were identified and assessed by PCR for the presence of polyketide synthase (PKS) genes. The fungi were all various genera of ascomycetes, the genomes of which encoded 36 putative PKS sequences, 26 shared sequence homology with beta-ketoacyl synthase domains, while 10 sequences showed homology to known fungal C-methyltransferase domains. A neighbour-joining phylogenetic analysis of the translated sequences could group the domains into previously established chemistry-based clades that represented non-reducing, partially reducing and highly reducing fungal PKSs. We observed that, in many cases, the membership of each clade also reflected the taxonomy of the fungal isolates. The functional assignment of the domains was further confirmed by in silico secondary and tertiary protein structure predictions. This genome mining study reveals, for the first time, the genetic potential of specific taxonomic groups of sugarcane-derived fungi to produce specific types of polyketides. Future work will focus on isolating these compounds with a view to understanding their chemical ecology and likely biotechnological potential.

The leishmanicidal flavonols quercetin and quercitrin target Leishmania (Leishmania) amazonensis arginase

Polyamine biosynthesis enzymes are promising drug targets for the treatment of leishmaniasis, Chagas' disease and African sleeping sickness. Arginase, which is a metallohydrolase, is the first enzyme involved in polyamine biosynthesis and converts arginine into ornithine and urea. Ornithine is used in the polyamine pathway that is essential for cell proliferation and ROS detoxification by trypanothione. The flavonols quercetin and quercitrin have been described as antitrypanosomal and antileishmanial compounds, and their ability to inhibit arginase was tested in this work. We characterized the inhibition of recombinant arginase from Leishmania (Leishmania) amazonensis by quercetin, quercitrin and isoquercitrin. The IC50 values for quercetin, quercitrin and isoquercitrin were estimated to be 3.8, 10 and 4.3 mu M, respectively. Quercetin is a mixed inhibitor, whereas quercitrin and isoquercitrin are uncompetitive inhibitors of L. (L.) amazonensis arginase. Quercetin interacts with the substrate L-arginine and the cofactor Mn2+ at pH 9.6, whereas quercitrin and isoquercitrin do not interact with the enzyme's cofactor or substrate. Docking analysis of these flavonols suggests that the cathecol group of the three compounds interact with Asp129, which is involved in metal bridge formation for the cofactors Mn-A(2+) and Mn-B(2+) in the active site of arginase. These results help to elucidate the mechanism of action of leishmanicidal flavonols and offer new perspectives for drug design against Leishmania infection based on interactions between arginase and flavones. (C) 2012 Elsevier Inc. All rights reserved.

Calcium Binding to Leptospira Outer Membrane Antigen LipL32 Is Not Necessary for Its Interaction with Plasma Fibronectin, Collagen Type IV, and Plasminogen

LipL32 is the most abundant outer membrane protein from pathogenic Leptospira and has been shown to bind extracellular matrix (ECM) proteins as well as Ca2+. Recent crystal structures have been obtained for the protein in the apo-and Ca2+-bound forms. In this work, we produced three LipL32 mutants (D163-168A, Q67A, and S247A) and evaluated their ability to interact with Ca2+ and with ECM glycoproteins and human plasminogen. The D163-168A mutant modifies aspartate residues involved in Ca2+ binding, whereas the other two modify residues in a cavity on the other side of the protein structure. Loss of calcium binding in the D163-D168A mutant was confirmed using intrinsic tryptophan fluorescence, circular dichroism, and thermal denaturation whereas the Q67A and S247A mutants presented the same Ca2+ affinity as the wild-type protein. We then evaluated if Ca2+ binding to LipL32 would be crucial for its interaction with collagen type IV and plasma proteins fibronectin and plasminogen. Surprisingly, the wild-type protein and all three mutants, including the D163-168A variant, bound to these ECM proteins with very similar affinities, both in the presence and absence of Ca2+ ions. In conclusion, calcium binding to LipL32 may be important to stabilize the protein, but is not necessary to mediate interaction with host extracellular matrix proteins.

Mutation Spectrum in the Large GTPase Dynamin 2, and Genotype-Phenotype Correlation in Autosomal Dominant Centronuclear Myopathy

Centronuclear myopathy (CNM) is a genetically heterogeneous disorder associated with general skeletal muscle weakness, type I fiber predominance and atrophy, and abnormally centralized nuclei. Autosomal dominant CNM is due to mutations in the large GTPase dynamin 2 (DNM2), a mechanochemical enzyme regulating cytoskeleton and membrane trafficking in cells. To date, 40 families with CNM-related DNM2 mutations have been described, and here we report 60 additional families encompassing a broad genotypic and phenotypic spectrum. In total, 18 different mutations are reported in 100 families and our cohort harbors nine known and four new mutations, including the first splice-site mutation. Genotype-phenotype correlation hypotheses are drawn from the published and new data, and allow an efficient screening strategy for molecular diagnosis. In addition to CNM, dissimilar DNM2 mutations are associated with Charcot-Marie-Tooth (CMT) peripheral neuropathy (CMTD1B and CMT2M), suggesting a tissue-specific impact of the mutations. In this study, we discuss the possible clinical overlap of CNM and CMT, and the biological significance of the respective mutations based on the known functions of dynamin 2 and its protein structure. Defects in membrane trafficking due to DNM2 mutations potentially represent a common pathological mechanism in CNM and CMT. Hum Mutat 33: 949-959, 2012. (C) 2012 Wiley Periodicals, Inc.

Mechanical Alterations in airway smooth muscle after interaction with indoor pollutant – A mathematical model.

The viscoelasticity of mammalian lung is determined by the mechanical properties and structural regulation of the airway smooth muscle (ASM). The exposure to polluted air may deteriorate these properties with harmful consequences to individual health. Formaldehyde (FA) is an important indoor pollutant found among volatile organic compounds. This pollutant permeates through the smooth muscle tissue forming covalent bonds between proteins in the extracellular matrix and intracellular protein structure changing mechanical properties of ASM and inducing asthma symptoms, such as airway hyperresponsiveness, even at low concentrations. In the experimental scenario, the mechanical effect of FA is the stiffening of the tissue, but the mechanism behind this effect is not fully understood. Thus, the aim of this study is to reproduce the mechanical behavior of the ASM, such as contraction and stretching, under FA action or not. For this, it was created a two-dimensional viscoelastic network model based on Voronoi tessellation solved using Runge-Kutta method of fourth order. The equilibrium configuration was reached when the forces in different parts of the network were equal. This model simulates the mechanical behavior of ASM through of a network of dashpots and springs. This dashpot-spring mechanical coupling mimics the composition of the actomyosin machinery of ASM through the contraction of springs to a minimum length. We hypothesized that formation of covalent bonds, due to the FA action, can be represented in the model by a simple change in the elastic constant of the springs, while the action of methacholine (MCh) reduce the equilibrium length of the spring. A sigmoid curve of tension as a function of MCh doses was obtained, showing increased tension when the muscle strip was exposed to FA. Our simulations suggest that FA, at a concentration of 0.1 ppm, can affect the elastic properties of the smooth muscle ¯bers by a factor of 120%. We also analyze the dynamic mechanical properties, observing the viscous and elastic behavior of the network. Finally, the proposed model, although simple, incorporates the phenomenology of both MCh and FA and reproduces experimental results observed with in vitro exposure of smooth muscle to FA. Thus, this new mechanical approach incorporates several well know features of the contractile system of the cells in a tissue level model. The model can also be used in different biological scales.

Mechanical alterations in airway smooth muscle after interaction with indoor pollutant - A mathematical model

The viscoelasticity of mammalian lung is determined by the mechanical properties and structural regulation of the airway smooth muscle (ASM). The exposure to polluted air may deteriorate these properties with harmful consequences to individual health. Formaldehyde (FA) is an important indoor pollutant found among volatile organic compounds. This pollutant permeates through the smooth muscle tissue forming covalent bonds between proteins in the extracellular matrix and intracellular protein structure changing mechanical properties of ASM and inducing asthma symptoms, such as airway hyperresponsiveness, even at low concentrations. In the experimental scenario, the mechanical effect of FA is the stiffening of the tissue, but the mechanism behind this effect is not fully w1derstood. Thus, the aim of this study is to reproduce the mechanical behavior of the ASM, such as contraction and stretching, under FA action or not. For this, it was created a two-dimensional viscoelastic network model based on Voronoi tessellation solved using Runge-Kutta method of fourth order. The equilibrium configuration was reached when the forces in different parts of the network were equal. This model simulates the mechanical behavior of ASM through of a network of dashpots and springs. This dashpot-spring mechanical coupling mimics the composition of the actomyosin machinery of ASM through the contraction of springs to a minimum length. We hypothesized that formation of covalent bonds, due to the FA action, can be represented in the model by a simple change in the elastic constant of the springs, while the action of methacholinc (MCh) reduce the equilibrium length of the spring. A sigmoid curve of tension as a function of MCh doses was obtained, showing increased tension when the muscle strip was exposed to FA. Our simulations suggest that FA, at a concentration of 0.1 ppm, can affect the elastic properties of the smooth muscle fibers by a factor of 120%. We also analyze the dynamic mechanical properties, observing the viscous and elastic behavior of the network. Finally, the proposed model, although simple, ir1corporates the phenomenology of both MCh and FA and reproduces experirnental results observed with ir1 vitro exposure of smooth muscle to .FA. Thus, this new mechanical approach incorporates several well know features of the contractile system of the cells ir1 a tissue level model. The model can also be used in different biological scales.

Prediction of apparent protein digestibility of ingredients and diets by in vitro pH-stat degree of protein hydrolysis with species-specific enzymes for juvenile Pacific white shrimp Litopenaeus vannamei

Aquafeed production faces global issues related to availability of feed ingredients. Feed manufacturers require greater flexibility in order to develop nutritional and cost-effective formulations that take into account nutrient content and availability of ingredients. The search for appropriate ingredients requires detailed screening of their potential nutritional value and variability at the industrial level. In vitro digestion of feedstuffs by enzymes extracted from the target species has been correlated with apparent protein digestibility (APD) in fish and shrimp species. The present study verified the relationship between APD and in vitro degree of protein hydrolysis (DH) with Litopenaeus vannamei hepatopancreas enzymes in several different ingredients (n = 26): blood meals, casein, corn gluten meal, crab meal, distiller`s dried grains with solubles, feather meal, fish meals, gelatin, krill meals, poultry by-product meal, soybean meals, squid meals and wheat gluten. The relationship between APD and DH was further verified in diets formulated with these ingredients at 30% inclusion into a reference diet. APD was determined in vivo (30.1 +/- 0.5 degrees C, 32.2 +/- 0.4%.) with juvenile L vannamei (9 to 12 g) after placement of test ingredients into a reference diet (35 g kg(-1) CP: 8.03 g kg(-1) lipid; 2.01 kcal g(-1)) with chromic oxide as the inert marker. In vitro DH was assessed in ingredients and diets with standardized hepatopancreas enzymes extracted from pond-reared shrimp. The DH of ingredients was determined under different assay conditions to check for the most suitable in vitro protocol for APD prediction: different batches of enzyme extracts (HPf5 or HPf6), temperatures (25 or 30 degrees C) and enzyme activity (azocasein): crude protein ratios (4 U: 80 mg CP or 4 U: 40 mg CP). DH was not affected by ingredient proximate composition. APD was significantly correlated to DH in regressions considering either ingredients or diets. The relationships between APD and DH of the ingredients could be suitably adjusted to a Rational Function (y = (a + bx)/(1 + cx + dx2), n = 26. Best in vitro APD predictions were obtained at 25 degrees C, 4 U: 80 mg CP both for ingredients (R(2) = 0.86: P = 0.001) and test diets (R(2) = 0.96; P = 0.007). The regression model including all 26 ingredients generated higher prediction residuals (i.e., predicted APD - determined APD) for corn gluten meal, feather meal. poultry by-product meal and krill flour. The remaining test ingredients presented mean prediction residuals of 3.5 points. A model including only ingredients with APD>80% showed higher prediction precision (R(2) = 0.98: P = 0.000004; n = 20) with average residual of 1.8 points. Predictive models including only ingredients from the same origin (e.g., marine-based, R(2) = 0.98; P = 0.033) also displayed low residuals. Since in vitro techniques have been usually validated through regressions against in vivo APD, the DH predictive capacity may depend on the consistency of the in vivo methodology. Regressions between APD and DH suggested a close relationship between peptide bond breakage by hepatopancreas digestive proteases and the apparent nitrogen assimilation in shrimp, and this may be a useful tool to provide rapid nutritional information. (C) 2009 Elsevier B.V. All rights reserved.

Prediction of culture performance of juvenile Litopenaeus vannamei by in vitro (pH-stat) degree of feed protein hydrolysis with species-specific enzymes

Rapid in vitro methods for measuring digestibility may be useful in analysing aqua feeds if the extent and limits of their application are clearly defined. The pH-stat protein digestibility routine with shrimp hepatopancreas enzymes was previously related to apparent protein digestibility with juvenile Litopenaeus vannamei fed diets containing different protein ingredients. The potential of the method to predict culture performance of shrimp fed six commercial feeds (T3, T4, T5, T6, T7 and T8) with 350 g kg(-1) declared crude-protein content was assessed. The consistency of results obtained using hepatopancreas enzyme extracts from either pond or clear water-raised shrimp was further verified in terms of reproducibility and possible diet history effects upon in vitro outputs. Shrimps were previously acclimated and then maintained over 56 days (initial mean weight 3.28 g) on each diet in 500-L tanks at 114 ind m(-2), clear water closed system with continuous renewal and mechanical filtering (50 mu m), with four replicates per treatment. Feeds were offered four times daily (six days a week) delivered in trays at feeding rates ranging from 4.0% to 7.0% of stocked shrimp biomass. Feed was accessible to shrimp 4 h daily for 1-h feeding period after which uneaten feed was recovered. Growth and survival were determined every 14 days from a sample of 16 individuals per tank. Water quality was monitored daily (pH, temperature and salinity) and managed by water back flushing filter cleaning every 7-10 days. Feeds were analysed for crude protein, gross energy, amino acids and pepsin digestibility. In vitro pH-stat degree of protein hydrolysis (DH%) was determined for each feed using hepatopancreas enzyme extracts from experimental (clear water) or pond-raised shrimp. Feeds resulted in significant differences in shrimp performance (P < 0.05) as seen by the differences in growth rates (0.56-0.98 g week(-1)), final weight and feed conversion ratio (FCR). Shrimp performance and in vitro DH% with pond-raised shrimp enzymes showed significant correlation (P < 0.05) for yield (R-2 = 0.72), growth rates (R-2 = 0.72-0.80) and FCR (R-2 = -0.67). Other feed attributes (protein : energy ratio, amino acids, true protein, non-protein nitrogen contents and in vitro pepsin digestibility) showed none or limited correlation with shrimp culture performance. Additional correlations were found between growth rates and methionine (R-2 = 0.73), FCR and histidine (R-2 = -0.60), and DH% and methionine or methionine+cystine feed contents (R-2 = 0.67-0.92). pH-stat assays with shrimp enzymes generated reproducible DH% results with either pond (CV <= 6.5%) or clear water (CV <= 8.5%) hepatopancreas enzyme sources. Moreover, correlations between shrimp growth rates and feed DH% were significant regardless of the enzyme origin (pond or clear water-raised shrimp) and showed consistent R-2 values. Results suggest the feasibility of using standardized hepatopancreas enzyme extracts for in vitro protein digestibility.

Water-in-Oil Micro-Emulsion Enhances the Secondary Structure of a Protein by Confinement

A scheme is presented in which an organic solvent environment in combination with surfactants is used to confine a natively unfolded protein inside an inverse microemulsion droplet. This type of confinement allows a study that provides unique insight into the dynamic structure of an unfolded, flexible protein which is still solvated and thus under near-physiological conditions. In a model system, the protein osteopontin (OPN) is used. It is a highly phosphorylated glycoprotein that is expressed in a wide range of cells and tissues for which limited structural analysis exists due to the high degree of flexibility and large number of post-translational modifications. OPN is implicated in tissue functions, such as inflammation and mineralisation. It also has a key function in tumour metastasis and progression. Circular dichroism measurements show that confinement enhances the secondary structural features of the protein. Small-angle X-ray scattering and dynamic light scattering show that OPN changes from being a flexible protein in aqueous solution to adopting a less flexible and more compact structure inside the microemulsion droplets. This novel approach for confining proteins while they are still hydrated may aid in studying the structure of a wide range of natively unfolded proteins.

Identification of novel components of NAD-utilizing metabolic pathways and prediction of their biochemical functions

Nicotinamide adenine dinucleotide (NAD) is a ubiquitous cofactor participating in numerous redox reactions. It is also a substrate for regulatory modifications of proteins and nucleic acids via the addition of ADP-ribose moieties or removal of acyl groups by transfer to ADP-ribose. In this study, we use in-depth sequence, structure and genomic context analysis to uncover new enzymes and substrate-binding proteins in NAD-utilizing metabolic and macromolecular modification systems. We predict that Escherichia coli YbiA and related families of domains from diverse bacteria, eukaryotes, large DNA viruses and single strand RNA viruses are previously unrecognized components of NAD-utilizing pathways that probably operate on ADP-ribose derivatives. Using contextual analysis we show that some of these proteins potentially act in RNA repair, where NAD is used to remove 2'-3' cyclic phosphodiester linkages. Likewise, we predict that another family of YbiA-related enzymes is likely to comprise a novel NAD-dependent ADP-ribosylation system for proteins, in conjunction with a previously unrecognized ADP-ribosyltransferase. A similar ADP-ribosyltransferase is also coupled with MACRO or ADP-ribosylglycohydrolase domain proteins in other related systems, suggesting that all these novel systems are likely to comprise pairs of ADP-ribosylation and ribosylglycohydrolase enzymes analogous to the DraG-DraT system, and a novel group of bacterial polymorphic toxins. We present evidence that some of these coupled ADP-ribosyltransferases/ribosylglycohydrolases are likely to regulate certain restriction modification enzymes in bacteria. The ADP-ribosyltransferases found in these, the bacterial polymorphic toxin and host-directed toxin systems of bacteria such as Waddlia also throw light on the evolution of this fold and the origin of eukaryotic polyADP-ribosyltransferases and NEURL4-like ARTs, which might be involved in centrosomal assembly. We also infer a novel biosynthetic pathway that might be involved in the synthesis of a nicotinate-derived compound in conjunction with an asparagine synthetase and AMPylating peptide ligase. We use the data derived from this analysis to understand the origin and early evolutionary trajectories of key NAD-utilizing enzymes and present targets for future biochemical investigations.