Active Site Loop Research Articles

Conformational Selection of a Tryptophan Side Chain Drives the Generalized Increase in Activity of PET Hydrolases through a Ser/Ile Double Mutation.

Poly(ethylene terephthalate) (PET) is the most common polyester plastic in the packaging industry and a major source of environmental pollution due to its single use. Several enzymes, termed PET hydrolases, have been found to hydrolyze this polymer at different temperatures, with the enzyme from Ideonella sakaiensis (IsPETase) having optimal catalytic activity at 30-35 °C. Crystal structures of IsPETase have revealed that the side chain of a conserved tryptophan residue within an active site loop (W185) shifts between three conformations to enable substrate binding and product release. This is facilitated by two residues unique to IsPETase, S214 and I218. When these residues are inserted into other PET hydrolases in place of the otherwise strictly conserved histidine and phenylalanine residues found at their respective positions, they enhance activity and decrease T opt. Herein, we combine molecular dynamics and well-tempered metadynamics simulations to investigate dynamic changes of the S214/I218 and H214/F218 variants of IsPETase, as well as three other mesophilic and thermophilic PET hydrolases, at their respective temperature and pH optima. Our simulations show that the S214/I218 insertion both increases the flexibility of active site loop regions harboring key catalytic residues and the conserved tryptophan and expands the conformational plasticity of this tryptophan side chain, enabling the conformational transitions that allow for substrate binding and product release in IsPETase. The observed catalytic enhancement caused by this substitution in other PET hydrolases appears to be due to conformational selection, by capturing the conformational ensemble observed in IsPETase.

Room-temperature serial synchrotron crystallography of the human phosphatase PTP1B.

Room-temperature X-ray crystallography provides unique insights into protein conformational heterogeneity, but obtaining sufficiently large protein crystals is a common hurdle. Serial synchrotron crystallography (SSX) helps to address this hurdle by allowing the use of many medium- to small-sized crystals. Here, a recently introduced serial sample-support chip system has been used to obtain the first SSX structure of a human phosphatase, specifically protein tyrosine phosphatase 1B (PTP1B) in the unliganded (apo) state. In previous apo room-temperature structures, the active site and allosteric sites adopted alternate conformations, including open and closed conformations of the active-site WPD loop and of a distal allosteric site. By contrast, in our SSX structure the active site is best fitted with a single conformation, but the distal allosteric site is best fitted with alternate conformations. This observation argues for additional nuance in interpreting the nature of allosteric coupling in this protein. Overall, our results illustrate the promise of serial methods for room-temperature crystallography, as well as future avant-garde crystallography experiments, for PTP1B and other proteins.

An active site loop toggles between conformations to control antibiotic hydrolysis and inhibition potency for CTX-M β-lactamase drug-resistance enzymes

β-lactamases inactivate β-lactam antibiotics leading to drug resistance. Consequently, inhibitors of β-lactamases can combat this resistance, and the β-lactamase inhibitory protein (BLIP) is a naturally occurring inhibitor. The widespread CTX-M-14 and CTX-M-15 β-lactamases have an 83% sequence identity. In this study, we show that BLIP weakly inhibits CTX-M-14 but potently inhibits CTX-M-15. The structure of the BLIP/CTX-M-15 complex reveals that binding is associated with a conformational change of an active site loop of β-lactamase. Surprisingly, the loop structure in the complex is similar to that in a drug-resistant variant (N106S) of CTX-M-14. We hypothesized that the pre-established favorable loop conformation of the N106S mutant would facilitate binding. The N106S substitution results in a ~100- and 10-fold increase in BLIP inhibition potency for CTX-M-14 and CTX-M-15, respectively. Thus, this indicates that an active site loop in β-lactamase toggles between conformations that control antibiotic hydrolysis and inhibitor susceptibility. These findings highlight the role of accessible active site conformations in controlling enzyme activity and inhibitor susceptibility as well as the influence of mutations in selectively stabilizing discrete conformations.

Thermostabilizing mechanisms of canonical single amino acid substitutions at a GH1 β-glucosidase probed by multiple MD and computational approaches.

β-glucosidases play a pivotal role in second-generation biofuel (2G-biofuel) production. For this application, thermostable enzymes are essential due to the denaturing conditions on the bioreactors. Random amino acid substitutions have originated new thermostable β-glucosidases, but without a clear understanding of their molecular mechanisms. Here, we probe by different molecular dynamics simulation approaches with distinct force fields and submitting the results to various computational analyses, the molecular bases of the thermostabilization of the Paenibacillus polymyxa GH1 β-glucosidase by two-point mutations E96K (TR1) and M416I (TR2). Equilibrium molecular dynamic simulations (eMD) at different temperatures, principal component analysis (PCA), virtual docking, metadynamics (MetaDy), accelerated molecular dynamics (aMD), Poisson-Boltzmann surface analysis, grid inhomogeneous solvation theory and colony method estimation of conformational entropy allow to converge to the idea that the stabilization carried by both substitutions depend on different contributions of three classic mechanisms: (i) electrostatic surface stabilization; (ii) efficient isolation of the hydrophobic core from the solvent, with energetic advantages at the solvation cap; (iii) higher distribution of the protein dynamics at the mobile active site loops than at the protein core, with functional and entropic advantages. Mechanisms i and ii predominate for TR1, while in TR2, mechanism iii is dominant. Loop A integrity and loops A, C, D, and E dynamics play critical roles in such mechanisms. Comparison of the dynamic and topological changes observed between the thermostable mutants and the wildtype protein with amino acid co-evolutive networks and thermostabilizing hotspots from the literature allow inferring that the mechanisms here recovered can be related to the thermostability obtained by different substitutions along the whole family GH1. We hope the results and insights discussed here can be helpful for future rational approaches to the engineering of optimized β-glucosidases for 2G-biofuel production for industry, biotechnology, and science.

Role of Half-of-Sites Reactivity and Inter-Subunit Communications in DAHP Synthase Catalysis and Regulation.

α-Carboxyketose synthases, including 3-deoxy-d-arabinoheptulosonate 7-phosphate synthase (DAHPS), are long-standing targets for inhibition. They are challenging targets to create tight-binding inhibitors against, and inhibitors often display half-of-sites binding and partial inhibition. Half-of-sites inhibition demonstrates the existence of inter-subunit communication in DAHPS. We used X-ray crystallography and spatially resolved hydrogen-deuterium exchange (HDX) to reveal the structural and dynamic bases for inter-subunit communication in Escherichia coli DAHPS(Phe), the isozyme that is feedback-inhibited by phenylalanine. Crystal structures of this homotetrameric (dimer-of-dimers) enzyme are invariant over 91% of its sequence. Three variable loops make up 8% of the sequence and are all involved in inter-subunit contacts across the tight-dimer interface. The structures have pseudo-twofold symmetry indicative of inter-subunit communication across the loose-dimer interface, with the diagonal subunits B and C always having the same conformation as each other, while subunits A and D are variable. Spatially resolved HDX reveals contrasting responses to ligand binding, which, in turn, affect binding of the second substrate, erythrose-4-phosphate (E4P). The N-terminal peptide, M1-E12, and the active site loop that binds E4P, F95-K105, are key parts of the communication network. Inter-subunit communication appears to have a catalytic role in all α-carboxyketose synthase families and a regulatory role in some members.

Structure of lactate oxidase from Enterococcus hirae revealed new aspects of active site loop function: Product-inhibition mechanism and oxygen gatekeeper.

l-Lactate oxidase (LOx) is a flavin mononucleotide (FMN)-dependent triose phosphate isomerase (TIM) barrel fold enzyme that catalyzes the oxidation of l-lactate using oxygen as a primary electron acceptor. Although reductive half-reaction mechanism of LOx has been studied by structure-based kinetic studies, oxidative half-reaction and substrate/product-inhibition mechanisms were yet to be elucidated. In this study, the structure and enzymatic properties of wild-type and mutant LOxs from Enterococcus hirae (EhLOx) were investigated. EhLOx structure showed the common TIM-barrel fold with flexible loop region. Noteworthy observations were that the EhLOx crystal structures prepared by co-crystallization with product, pyruvate, revealed the complex structures with "d-lactate form ligand," which was covalently bonded with a Tyr211 side chain. This observation provided direct evidence to suggest the product-inhibition mode of EhLOx. Moreover, this structure also revealed a flip motion of Met207 side chain, which is located on the flexible loop region as well as Tyr211. Through a saturation mutagenesis study of Met207, one of the mutants Met207Leu showed the drastically decreased oxidase activity but maintained dye-mediated dehydrogenase activity. The structure analysis of EhLOx Met207Leu revealed the absence of flipping in the vicinity of FMN, unlike the wild-type Met207 side chain. Together with the simulation of the oxygen-accessible channel prediction, Met207 may play as an oxygen gatekeeper residue, which contributes oxygen uptake from external enzyme to FMN. Three clades of LOxs are proposed based on the difference of the Met207 position and they have different oxygen migration pathway from external enzyme to active center FMN.

HIV-1 protease with 10 lopinavir and darunavir resistance mutations exhibits altered inhibition, structural rearrangements and extreme dynamics

Antiretroviral drug resistance is a therapeutic obstacle for people with HIV. HIV protease inhibitors darunavir and lopinavir are recommended for resistant infections. We characterized a protease mutant (PR10x) derived from a highly resistant clinical isolate including 10 mutations associated with resistance to lopinavir and darunavir. Compared to the wild-type protease, PR10x exhibits ∼3-fold decrease in catalytic efficiency and Ki values of 2–3 orders of magnitude worse for darunavir, lopinavir, and potent investigational inhibitor GRL-519. Crystal structures of the mutant were solved in a ligand-free form and in complex with GRL-519. The structures show altered interactions in the active site, flap-core interface, hydrophobic core, hinge region, and 80s loop compared to the corresponding wild-type protease structures. The ligand-free crystal structure exhibits a highly curled flap conformation which may amplify drug resistance. Molecular dynamics simulations performed for 1 μs on ligand-free dimers showed extremely large fluctuations in the flaps for PR10x compared to equivalent simulations on PR with a single L76V mutation or wild-type protease. This analysis offers insight about the synergistic effects of mutations in highly resistant variants.

An Extended C-Terminus, the Possible Culprit for Differential Regulation of 5-Aminolevulinate Synthase Isoforms.

5-Aminolevulinate synthase (ALAS; E.C. 2.3.1.37) is a pyridoxal 5′-phosphate (PLP)-dependent enzyme that catalyzes the key regulatory step of porphyrin biosynthesis in metazoa, fungi, and α-proteobacteria. ALAS is evolutionarily related to transaminases and is therefore classified as a fold type I PLP-dependent enzyme. As an enzyme controlling the key committed and rate-determining step of a crucial biochemical pathway ALAS is ideally positioned to be subject to allosteric feedback inhibition. Extensive kinetic and mutational studies demonstrated that the overall enzyme reaction is limited by subtle conformational changes of a hairpin loop gating the active site. These findings, coupled with structural information, facilitated early prediction of allosteric regulation of activity via an extended C-terminal tail unique to eukaryotic forms of the enzyme. This prediction was subsequently supported by the discoveries that mutations in the extended C-terminus of the erythroid ALAS isoform (ALAS2) cause a metabolic disorder known as X-linked protoporphyria not by diminishing activity, but by enhancing it. Furthermore, kinetic, structural, and molecular modeling studies demonstrated that the extended C-terminal tail controls the catalytic rate by modulating conformational flexibility of the active site loop. However, the precise identity of any such molecule remains to be defined. Here we discuss the most plausible allosteric regulators of ALAS activity based on divergences in AlphaFold-predicted ALAS structures and suggest how the mystery of the mechanism whereby the extended C-terminus of mammalian ALASs allosterically controls the rate of porphyrin biosynthesis might be unraveled.

Interplay between the Enamine and Imine Forms of the Hydrolyzed Imipenem in the Active Sites of Metallo-β-lactamases and in Water Solution.

Deactivation of the β-lactam antibiotics in the active sites of the β-lactamases is among the main mechanisms of bacterial antibiotic resistance. As drugs of last resort, carbapenems are efficiently hydrolyzed by metallo-β-lactamases, presenting a serious threat to human health. Our study reveals mechanistic aspects of the imipenem hydrolysis by bizinc metallo-β-lactamases, NDM-1 and L1, belonging to the B1 and the B3 subclasses, respectively. The results of QM(PBE0-D3/6-31G**)/MM simulations show that the enamine product with the protonated nitrogen atom is formed as the major product in NDM-1 and as the only product in the L1 active site. In NDM-1, there is also another reaction pathway that leads to the formation of the (S)-enantiomer of the imine form of the hydrolyzed imipenem; this process occurs with the higher energy barriers. The absence of the second pathway in L1 is due to the different amino acid composition of the active site loop. In L1, the hydrophobic Pro226 residue is located above the pyrroline ring of imipenem that blocks protonation of the carbon atom. Electron density analysis is performed at the stationary points to compare reaction pathways in L1 and NDM-1. Tautomerization from the enamine to the imine form likely happens in solution after the dissociation of the hydrolyzed imipenem from the active site of the enzyme. Classical molecular dynamics simulations of the hydrolyzed imipenem in solution, both with the neutral enamine and the negatively charged N-C2-C3 fragment, demonstrate a huge diversity of conformations. The vast majority of conformations blocks the C3-atom from the side required for the (S)-imine formation upon tautomerization. Thus, according to our calculations, formation of the (R)-imine is more likely. QM(PBE0-D3/6-31G**)/MM molecular dynamics simulations of the hydrolyzed imipenem with the negatively charged N-C2-C3 fragment followed by the Laplacian bond order analysis demonstrate that the N═C2-C3- resonance structure is the most pronounced that facilitates formation of the imine form. The proposed mechanism of the enzymatic enamine formation and its subsequent tautomerization to the imine form in solution is in agreement with the recent spectroscopic and NMR studies.

Biochemical, Enzymatic, and Computational Characterization of Recurrent Somatic Mutations of the Human Protein Tyrosine Phosphatase PTP1B in Primary Mediastinal B Cell Lymphoma.

Human protein tyrosine phosphatase 1B (PTP1B) is a ubiquitous non-receptor tyrosine phosphatase that serves as a major negative regulator of tyrosine phosphorylation cascades of metabolic and oncogenic importance such as the insulin, epidermal growth factor receptor (EGFR), and JAK/STAT pathways. Increasing evidence point to a key role of PTP1B-dependent signaling in cancer. Interestingly, genetic defects in PTP1B have been found in different human malignancies. Notably, recurrent somatic mutations and splice variants of PTP1B were identified in human B cell and Hodgkin lymphomas. In this work, we analyzed the molecular and functional levels of three PTP1B mutations identified in primary mediastinal B cell lymphoma (PMBCL) patients and located in the WPD-loop (V184D), P-loop (R221G), and Q-loop (G259V). Using biochemical, enzymatic, and molecular dynamics approaches, we show that these mutations lead to PTP1B mutants with extremely low intrinsic tyrosine phosphatase activity that display alterations in overall protein stability and in the flexibility of the active site loops of the enzyme. This is in agreement with the key role of the active site loop regions, which are preorganized to interact with the substrate and to enable catalysis. Our study provides molecular and enzymatic evidence for the loss of protein tyrosine phosphatase activity of PTP1B active-site loop mutants identified in human lymphoma.

Evolution of homo-oligomerization of methionine S-adenosyltransferases is replete with structure-function constrains.

Homomers are prevalent in bacterial proteomes, particularly among core metabolic enzymes. Homomerization is often key to function and regulation, and interfaces that facilitate the formation of homomeric enzymes are subject to intense evolutionary change. However, our understanding of the molecular mechanisms that drive evolutionary variation in homomeric complexes is still lacking. How is the diversification of protein interfaces linked to variation in functional regulation and structural integrity of homomeric complexes? To address this question, we studied quaternary structure evolution of bacterial methionine S‐adenosyltransferases (MATs)—dihedral homotetramers formed along a large and conserved dimeric interface harboring two active sites, and a small, recently evolved, interdimeric interface. Here, we show that diversity in the physicochemical properties of small interfaces is directly linked to variability in the kinetic stability of MAT quaternary complexes and in modes of their functional regulation. Specifically, hydrophobic interactions within the small interface of Escherichia coli MAT render the functional homotetramer kinetically stable yet impose severe aggregation constraints on complex assembly. These constraints are alleviated by electrostatic interactions that accelerate dimer‐dimer assembly. In contrast, Neisseria gonorrhoeae MAT adopts a nonfunctional dimeric state due to the low hydrophobicity of its small interface and the high flexibility of its active site loops, which perturbs small interface integrity. Remarkably, in the presence of methionine and ATP, N. gonorrhoeae MAT undergoes substrate‐induced assembly into a functional tetrameric state. We suggest that evolution acts on the interdimeric interfaces of MATs to tailor the regulation of their activity and stability to unique organismal needs.

Abstract 3265: Discovery of a first-in-class reversible DNMT1-selective inhibitor

Abstract DNA methylation, a key epigenetic driver of transcriptional silencing, is universally dysregulated in cancer. Reversal of DNA methylation by hypomethylating agents, such as the cytidine analogs decitabine or azacytidine, has demonstrated clinical benefit in hematologic malignancies. These nucleoside analogs are incorporated into replicating DNA where they inhibit DNA cytosine methyltransferases DNMT1, DNMT3A, and DNMT3B through irreversible covalent interactions. These agents induce significant toxicity to normal blood cells thus limiting their clinical doses. We report the discovery of a series of potent first-in-class DNMT1-selective inhibitors that compete with the active-site loop of DNMT1 for penetration into hemi-methylated DNA. These compounds induce robust loss of DNA methylation, transcriptional activation, and cancer cell growth inhibition in vitro leading to superior tumor regression and survival in mouse models of acute myeloid leukemia due to improved in vivo tolerability compared with decitabine. Citation Format: Melissa B. Pappalardi. Discovery of a first-in-class reversible DNMT1-selective inhibitor [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3265.

Evidence of Conformational Changes in the Structure of Y‐Family Human Polymerase Kappa

Y‐family DNA polymerases, unlike replicative DNA polymerases, are enzymes characterized by their ability to bypass DNA lesions and continue extension of damaged DNA via Translesion Synthesis (TLS). One such polymerase, human polymerase kappa (Pol κ), can efficiently bypass minor groove adducts, most notably N2‐furfuryl‐deoxy‐guanosine (N2‐ffdG), but it is inhibited by major groove adducts such as O6‐methylated‐deoxy‐guanosine (O6‐Met‐dG). Pol k is composed of five domains including its novel N‐clasp domain, as well as the thumb, palm, fingers, and little finger (polymerase associated domain or PAD) domains which are found across Y‐family polymerases. The specificity that Pol κ demonstrates is thought to be related to Y‐family polymerases having smaller finger domains when compared to replicative DNA polymerases thus limiting major groove contact, while Pol κ’s PAD and catalytic core move to open the minor groove side. It is believed that Pol κ and its Escherichia colianalogue, DinB (DNA polymerase IV), undergo one of two conformational changes, open or closed, in their active site loops when in contact with undamaged DNA or minor groove adducts. Using molecular dynamics (MD) simulations, we examined conformational changes and utilized RMSD analysis to observe activity within Pol κ’s domains across three different Pol κ ternary structures in the presence of different templating lesions and their correct incoming nucleotides. These structures include, one with a minor groove N2‐furfuryl‐deoxy‐guanosine adduct (FDG for short), one with a major groove O6‐methylated‐deoxy‐guanosine adduct, and finally one with a mutation of a residue in Pol κ’s PAD (R507K) with an FDG lesion on the DNA templating strand. All systems were solvated using CHARMMgui, minimized, equilibrated, and run for 110 ns with NAMD, and observed in VMD. Our findings support that Pol κ exists in two conformations, which depend on the kind of lesion on the DNA and/or the incoming nucleotide. In particular, the movement is concentrated in the N‐clasp and fingers domains. Some residues which seem important for this transition are N52, M135, N464, and R507. Interestingly, kinetics studies have already identified some of these same residues as important for extension. Our MD studies give an atomistic explanation of why such residues are crucial for the correct functioning of this protein.

The Multiple Roles of the Active Site Loop of Watermelon Glyoxysomal Malate Dehydrogenase

All malate dehydrogenases have a flexible loop that starts to close over the active site region of malate dehydrogenase when cofactor binds and locks down completing the active site when the dicarboxylic acid substrate binds. As such the active site loop appears to play a role in both substrate specificity and correct alignment of the substrate to allow hydride transfer and proton abstraction/donation depending on the direction of the reaction. Combined with available crystal structures, detailed analysis of the 1smk.pdb structure shows the loop in a series of microstates from fully open to a closed state induced by the inhibitor citrate binding (figure 1), sequence conservation suggests that the loop region contains a 19 amino acid stretch from P119 to N137 using the watermelon glyoxysomal enzyme as reference. Clustal Omega Analysis of malate dehydrogenases from prokaryotes to mammals shows 5 families of conservation of the loop region (figure 2) allowing a series of hypotheses about structure function relationships of specific residues within the sequence to be made which have been tested by a combination of computational and site directed mutagenesis wetlab approaches. In addition to the two canonical arginines involved in substrate binding, we have explored every other residue in the loop region. Overall, mutations in the first half of the loop had small impact on Specific Activity while mutations in the second half of the loop had larger impacts. R130 mutants to E or S had the most impact while mutation to A or Q had lesser effects. Three other residues in particular had striking effects on Oxaloacetate saturation, G121A, L133S and N137. L133 and N137, at the C terminal base of the loop, were hypothesized to play a role in the repulsive interactions between substrate and cofactor in the active site. Consistent with this hypothesis, creation of an L133S or D137A mutant showed that concomitant with a reduction in specific activity to 2% and 5% wildtype (respectively), Km for Oxaloacetate was dramatically increased while Km for NADH was decreased somewhat. Finally, a series of mutations of M128 (Q,I and A) were explored with varying effects on specific activity (Q>I>A) and significant effects on Km for Oxaloacetate but lesser or no effects on Km for NADH. The impact of M128 mutations in particular were further explored computationally using DeepDDG analysis and effects on stability correlated with impact on kinetic properties. Overall these studies provide a detailed picture of the complex interactions the flexible loop makes that contribute to both substrate binding and catalysis.

Model for Calcium‐Dependent Activation of a Type I Metacaspase From the Fungus Schizophyllum commune

Apoptosis or programmed cell death (PCD) is one of the most important cellular processes which functions to avoid the accumulation of mutations or damage within a cell. PCD has been extensively studied in eukaryotic cells, demonstrating that caspase proteins have a significant role in this process. However, unlike mammalian cells, caspase genes are absent in plants, fungi, and protozoa. Instead, these organisms contain homologous proteins known as metacaspases which have demonstrated similar proteolytic functions in aiding PCD. Determining an accurate model for how these metacaspases are activated could provide insight into the methods behind this protease‐dependent cell death which could lead to further strategies for the regulation of fungal infections. This research focuses on the model of activation for the Type I metacaspases from the fungus Schizophyllum commune. The S. commune Type I metacaspases are dependent on calcium ions for activation; however, the molecular mechanism behind how calcium activates this protease remains unknown.It has been proposed that when calcium binds to the metacapsape, these ions induce a structural rearrangement of the active‐site region leading to activation. To further understand how calcium and other divalent metallic cations influence the activation of this protease, we used kinetic activity assays to quantify metacaspase activity and differential scanning fluorimetry to measure divalent metallic cation‐metacaspase binding.In addition to strong activation by calcium, manganese ions have demonstrated a moderate effect on metaspace activity, while magnesium, copper, samarium, and cobalt don’t affect activity. We will present evidence that calcium binds to two sites on the protein, while manganese binds to only a single site. Furthermore, calcium binding to the Type I metacapses occurs at two sites with affinities that differ by several orders of magnitude. Initial results suggest that it is necessary for calcium to saturate the high affinity site prior to binding to the low affinity site. Thus, it has been predicted that calcium’s initial saturation of the high affinity site induces slight conformation changes to the active‐site; however, once the low affinity site has been saturated, an active‐site loop undergoes further rearrangements resulting in metacaspase activation. Further analysis of binding interactions between the Type I metacapse and calcium (along with other divalent cations) will provide insight into how these structural rearrangements lead to increased metacaspase activity.

Effects of the T337M and G391V disease-related variants on human phosphoglucomutase 1: structural disruptions large and small.

Phosphoglucomutase 1 (PGM1) plays a central role in glucose homeostasis in human cells. Missense variants of this enzyme cause an inborn error of metabolism, which is categorized as a congenital disorder of glycosylation. Here, two disease-related variants of PGM1, T337M and G391V, which are both located in domain 3 of the four-domain protein, were characterized via X-ray crystallography and biochemical assays. The studies show multiple impacts resulting from these dysfunctional variants, including both short- and long-range structural perturbations. In the T337M variant these are limited to a small shift in an active-site loop, consistent with reduced enzyme activity. In contrast, theG391V variant produces a cascade of structural perturbations, including displacement of both the catalytic phosphoserine and metal-binding loops. This work reinforces several themes that were found in prior studies of dysfunctional PGM1 variants, including increased structural flexibility and the outsized impacts of mutations affecting interdomain interfaces. The molecular mechanisms of PGM1 variants have implications for newly described inherited disorders of related enzymes.

Structural characterization of dicyanopyridine containing DNMT1-selective, non-nucleoside inhibitors

DNMT1 maintains the parental DNA methylation pattern on newly replicated hemimethylated DNA. The failure of this maintenance process causes aberrant DNA methylation that affects transcription and contributes to the development and progression of cancers such as acute myeloid leukemia. Here, we structurally characterized a set of newly discovered DNMT1-selective, reversible, non-nucleoside inhibitors that bear a core 3,5-dicyanopyridine moiety, as exemplified by GSK3735967, to better understand their mechanism of inhibition. All of the dicyanopydridine-containing inhibitors examined intercalate into the hemimethylated DNA between two CpG base pairs through the DNA minor groove, resulting in conformational movement of the DNMT1 active-site loop. In addition, GSK3735967 introduces two new binding sites, where it interacts with and stabilizes the displaced DNMT1 active-site loop and it occupies an open aromatic cage in which trimethylated histone H4 lysine 20 is expected to bind. Our work represents a substantial step in generating potent, selective, and non-nucleoside inhibitors of DNMT1.

Reducing substrate inhibition of malate dehydrogenase from Geobacillus stearothermophilus by C-terminal truncation.

Malate dehydrogenase (MDH) catalyzes the reduction of oxaloacetate to L-malate. Geobacillus stearothermophilus MDH (gs-MDH) is used as a diagnostic reagent; however, gs-MDH is robustly inhibited at high substrate concentrations, which limits its reaction rate. Here, we reduced substrate inhibition of gs-MDH by deleting its C-terminal residues. Computational analysis showed that C-terminal residues regulate the position of the active site loop. C-terminal deletions of gs-MDH successfully increased Ki values by 5- to 8-fold with maintained thermal stability (>90% of the wild-type enzyme), although kcat/Km values were decreased by <2-fold. The structure of the mutant showed a shift in the location of the active site loop and a decrease in its volume, suggesting that substrate inhibition was reduced by eliminating the putative substrate binding site causing inhibition. Our results provide an effective method to reduce substrate inhibition of the enzyme without loss of other parameters, including binding and stability constants.

Identification of new putative inhibitors of Mycobacterium tuberculosis 3-dehydroshikimate dehydratase from a combination of ligand- and structure-based and deep learning in silico approaches

The development of new drugs against Mycobacterium tuberculosis is an essential strategy for fighting drug resistance. Although 3-dehydroquinate dehydratase (MtDHQ) is known to be a highly relevant target for M. tuberculosis, current research shows new putative inhibitors of MtDHQ selected by a large-scale ensemble-docking strategy combining ligand- and target-based chemoinformatic methods to deep learning. Initial chemical library was reduced from 216 million to approximately 460 thousand after pharmacophore, toxicity and molecular weight filters. Final library was subjected to an ensemble-docking protocol in GOLD which selected the top 300 molecules (GHITS). GHITS displayed different structures and characteristics when compared to known inhibitors (KINH). GHITS were further screened by post-docking analysis in AMMOS2 and deep learning virtual screening in DeepPurpose. DeepPurpose predicted that a number of GHITS had comparable or better affinity for the target than KINH. The best molecule was selected by consensus ranking using GOLD, AMMOS2 and DeepPurpose scores. Molecular dynamics revealed that the top hit displayed consistent and stable binding to MtDHQ, making strong interactions with active-site loop residues. Results forward new putative inhibitors of MtDHQ and reinforce the potential application of artificial intelligence methods for drug design. This work represents the first step in the validation of these molecules as inhibitors of MtDHQ.

The Structural Basis of Babesia orientalis Lactate Dehydrogenase.

Glycolytic enzymes play a crucial role in the anaerobic glycolysis of apicomplexan parasites for energy generation. Consequently, they are considered as potential targets for new drug development. Previous studies revealed that lactate dehydrogenase (LDH), a glycolytic enzyme, is a potential drug target in different parasites, such as Plasmodium, Toxoplasma, Cryptosporidium, and Piroplasma. Herein, in order to investigate the structural basis of LDH in Babesia spp., we determined the crystal structure of apo Babesia orientalis (Bo) LDH at 2.67-Å resolution in the space group P1. A five-peptide insertion appears in the active pocket loop of BoLDH to create a larger catalytic pocket, like other protozoa (except for Babesia microti LDH) and unlike its mammalian counterparts, and the absence of this extra insertion inactivates BoLDH. Without ligands, the apo BoLDH takes R-state (relaxed) with the active-site loop open. This feature is obviously different from that of allosteric LDHs in T-state (tense) with the active-site loop open. Compared with allosteric LDHs, the extra salt bridges and hydrogen bonds make the subunit interfaces of BoLDH more stable, and that results in the absence of T-state. Interestingly, BoLDH differs significantly from BmLDH, as it exhibits the ability to adapt quickly to the synthetic co-factor APAD+. In addition, the enzymatic activity of BoLDH was inhibited non-competitively by polyphenolic gossypol with a K i value of 4.25 μM, indicating that BoLDH is sensitive to the inhibition of gossypol and possibly to its new derivative compounds. The current work provides the structural basis of BoLDH for the first time and suggests further investigation on the LDH structure of other Babesia spp. That knowledge would indeed facilitate the screening and designing of new LDH inhibitors to control the intracellular proliferation of Babesia spp.