Molecular dynamics analysis of the structural and dynamic properties of the functionally enhanced hepta-variant of mouse 5-aminolevulinate synthase

5-Aminolevulinate synthase catalysis: The catcher in heme biosynthesis

5-Aminolevulinate (ALA) synthase (ALAS), a homodimeric pyridoxal-5'-phosphate (PLP)-dependent enzyme, catalyzes the first step of heme biosynthesis in metazoa, fungi and α-proteobacteria. In this review, we focus on the advances made in unraveling the mechanism of the ALAS-catalyzed reaction during the past decade. The interplay between the PLP cofactor and the protein moiety determines and modulates the multi-intermediate reaction cycle of ALAS, which involves the decarboxylative condensation of two substrates, glycine and succinyl-CoA. Substrate binding and catalysis are rapid, and product (ALA) release dominates the overall ALAS kinetic mechanism. Interconversion between a catalytically incompetent, open conformation and a catalytically competent, closed conformation is linked to ALAS catalysis. Reversion to the open conformation, coincident with ALA dissociation, defines the slowest step of the reaction cycle. These findings were further substantiated by introducing seven mutations in the16-amino acid loop that gates the active site, yielding an ALAS variant with a greatly increased rate of catalytic turnover and heightened specificity constants for both substrates. Recently, molecular dynamics (MD) simulation analysis of various dimeric ALAS forms revealed that the seven active site loop mutations caused the proteins to adopt different conformations. In particular, the emergence of a β-strand in the mutated loop, which interacted with two preexisting β-strands to form an anti-parallel three-stranded β-sheet, conferred the murine heptavariant with a more stable open conformation and prompted faster product release than wild-type mALAS2. Moreover, the dynamics of the mALAS2 active site loop anti-correlated with that of the 35 amino acid C-terminal sequence. This led us to propose that this C-terminal extension, which is absent in prokaryotic ALASs, finely tunes mammalian ALAS activity. Based on the above results, we extend our previous proposal to include that discovery of a ligand inducing the mammalian C-terminal extension to fold offers a good prospect for the development of a new drug for X-linked protoporphyria and/or other porphyrias associated with enhanced ALAS activity and/or porphyrin accumulation.

Human p38α mitogen-activated protein kinase in the Asp168-Phe169-Gly170-in (DFG-in) state can bind allosteric inhibitor Doramapimod

Doramapimod (BIRB-796) is widely recognized as one of the most potent and selective type II inhibitors of human p38α mitogen-activated protein kinase (MAPK); however, the understanding of its binding mechanism remains incomplete. Previous studies indicated high affinity of the ligand to a so-called allosteric pocket revealed only in the 'out' state of the DFG motif (i.e. Asp168-Phe169-Gly170) when Phe169 becomes fully exposed to the solvent. The possibility of alternative binding in the DFG-in state was hypothesized, but the molecular mechanism was not known. Methods of bioinformatics, docking and long-time scale classical and accelerated molecular dynamics have been applied to study the interaction of Doramapimod with the human p38α MAPK. It was shown that Doramapimod can bind to the protein even when the Phe169 is fully buried inside the allosteric pocket and the kinase activation loop is in the DFG-in state. Orientation of the inhibitor in such a complex is significantly different from that in the known crystallographic complex formed by the kinase in the DFG-out state; however, the Doramapimod's binding is followed by the ligand-induced conformational changes, which finally improve accommodation of the inhibitor. Molecular modelling has confirmed that Doramapimod combines the features of type I and II inhibitors of p38α MAPK, i.e. can directly and indirectly compete with the ATP binding. It can be concluded that optimization of the initial binding in the DFG-in state and the final accommodation in the DFG-out state should be both considered at designing novel efficient type II inhibitors of MAPK and homologous proteins. Communicated by Ramaswamy H. Sarma.

Anti-Correlation between the Dynamics of the Active Site Loop and C-Terminal Tail in Relation to the Homodimer Asymmetry of the Mouse Erythroid 5-Aminolevulinate Synthase

Biosynthesis of heme represents a complex process that involves multiple stages controlled by different enzymes. The first of these proteins is a pyridoxal 5′-phosphate (PLP)-dependent homodimeric enzyme, 5-aminolevulinate synthase (ALAS), that catalyzes the rate-limiting step in heme biosynthesis, the condensation of glycine with succinyl-CoA. Genetic mutations in human erythroid-specific ALAS (ALAS2) are associated with two inherited blood disorders, X-linked sideroblastic anemia (XLSA) and X-linked protoporphyria (XLPP). XLSA is caused by diminished ALAS2 activity leading to decreased ALA and heme syntheses and ultimately ineffective erythropoiesis, whereas XLPP results from “gain-of-function” ALAS2 mutations and consequent overproduction of protoporphyrin IX and increase in Zn²⁺-protoporphyrin levels. All XLPP-linked mutations affect the intrinsically disordered C-terminal tail of ALAS2. Our earlier molecular dynamics (MD) simulation-based analysis showed that the activity of ALAS2 could be regulated by the conformational flexibility of the active site loop whose structural features and dynamics could be changed due to mutations. We also revealed that the dynamic behavior of the two protomers of the ALAS2 dimer differed. However, how the structural dynamics of ALAS2 active site loop and C-terminal tail dynamics are related to each other and contribute to the homodimer asymmetry remained unanswered questions. In this study, we used bioinformatics and computational biology tools to evaluate the role(s) of the C-terminal tail dynamics in the structure and conformational dynamics of the murine ALAS2 homodimer active site loop. To assess the structural correlation between these two regions, we analyzed their structural displacements and determined their degree of correlation. Here, we report that the dynamics of ALAS2 active site loop is anti-correlated with the dynamics of the C-terminal tail and that this anti-correlation can represent a molecular basis for the functional and dynamic asymmetry of the ALAS2 homodimer.

Investigating the bifunctionality of cyclizing and "classical" 5-aminolevulinate synthases

The precursor to all tetrapyrroles is 5-aminolevulinic acid, which is made either via the condensation of glycine and succinyl-CoA catalyzed by an ALA synthase (the C4 or Shemin pathway) or by a pathway that uses glutamyl-tRNA as a precursor and involves other enzymes (the C5 pathway). Certain ALA synthases also catalyze the cyclization of ALA-CoA to form 2-amino-3-hydroxycyclopent-2-en-1-one. Organisms with synthases that possess this second activity nevertheless rely upon the C5 pathway to supply ALA for tetrapyrrole biosynthesis. The C₅ N units are components of a variety of secondary metabolites. Here, we show that an ALA synthase used exclusively for tetrapyrrole biosynthesis is also capable of catalyzing the cyclization reaction, albeit at much lower efficiency than the dedicated cyclases. Two absolutely conserved serines present in all known ALA-CoA cyclases are threonines in all known ALA synthases, suggesting they could be important in distinguishing the functions of these enzymes. We found that purified mutant proteins having single and double substitutions of the conserved residues are not improved in their respective alternate activities; rather, they are worse. Protein structural modeling and amino acid sequence alignments were explored within the context of what is known about the reaction mechanisms of these two different types of enzymes to consider what other features are important for the two activities.