A catalytic mechanism for cysteine N-terminal nucleophile hydrolases, as revealed by free energy simulations

Discovery of Tetrahydro Tanshinone I as a Naturally Occurring Covalent Pan-Inhibitor Against Gut Microbial Bile Salt Hydrolases

Gut microbial bile salt hydrolases (gmBSHs), an important class of bacteria-produced cysteine hydrolases, play a crucial role in bile acid metabolism. Modulating the total gmBSH activity is a feasible way for ameliorating some metabolic diseases including colorectal cancer, type 2 diabetes, and obesity. This study reported the discovery and characterization of a botanical compound as a covalent pan-inhibitor of gmBSHs. Following the screening of more than 100 botanical compounds, tanshinones were found with strong time-dependent anti-EfBSH effects. After that, a total of 17 naturally occurring tanshinones were collected, and their anti-EfBSH potentials were tested. Among all tested tanshinones, tetrahydro tanshinone I (THTI) exhibited the most potent inhibitory effects against five gmBSHs (EfBSH, LsBSH, BtBSH, CpBSH, and BlBSH), showing the IC50 values ranging from 0.28 ± 0.05 μM to 1.62 ± 0.07 μM. Further investigations showed that THTI could covalently modify the conserved catalytic cysteine (Cys2) of all tested gmBSHs, while this agent could strongly inhibit the total gmBSHs activity in live microorganisms and murine gut luminal content. Collectively, THTI is identified as a naturally occurring covalent pan-inhibitor of gmBSHs, which offers a promising lead compound to develop more efficacious gmBSHs inhibitors for the management of bile acid metabolism and related metabolic disorders.

First Insight into the Degradome of Aspergillus ochraceus: Novel Secreted Peptidases and Their Inhibitors

Aspergillus fungi constitute a pivotal element within ecosystems, serving as both contributors of biologically active compounds and harboring the potential to cause various diseases across living organisms. The organism's proteolytic enzyme complex, termed the degradome, acts as an intermediary in its dynamic interaction with the surrounding environment. Using techniques such as genome and transcriptome sequencing, alongside protein prediction methodologies, we identified putative extracellular peptidases within Aspergillus ochraceus VKM-F4104D. Following manual annotation procedures, a total of 11 aspartic, 2 cysteine, 2 glutamic, 21 serine, 1 threonine, and 21 metallopeptidases were attributed to the extracellular degradome of A. ochraceus VKM-F4104D. Among them are enzymes with promising applications in biotechnology, potential targets and agents for antifungal therapy, and microbial antagonism factors. Thus, additional functionalities of the extracellular degradome, extending beyond mere protein substrate digestion for nutritional purposes, were demonstrated.

Bile salt hydrolase acyltransferase activity expands bile acid diversity

Bile acids (BAs) are steroid detergents in bile that contribute to the absorption of fats and fat-soluble vitamins while shaping the gut microbiome because of their antimicrobial properties1-4. Here we identify the enzyme responsible for a mechanism of BA metabolism by the gut microbiota involving amino acid conjugation to the acyl-site of BAs, thus producing a diverse suite of microbially conjugated bile acids (MCBAs). We show that this transformation is mediated by acyltransferase activity of bile salt hydrolase (bile salt hydrolase/transferase, BSH/T). Clostridium perfringens BSH/T rapidly performed acyl transfer when provided various amino acids and taurocholate, glycocholate or cholate, with an optimum at pH 5.3. Amino acid conjugation by C. perfringens BSH/T was diverse, including all proteinaceous amino acids except proline and aspartate. MCBA production was widespread among gut bacteria, with strain-specific amino acid use. Species with similar BSH/T amino acid sequences had similar conjugation profiles and several bsh/t alleles correlated with increased conjugation diversity. Tertiary structure mapping of BSH/T followed by mutagenesis experiments showed that active site structure affects amino acid selectivity. These MCBA products had antimicrobial properties, where greater amino acid hydrophobicity showed greater antimicrobial activity. Inhibitory concentrations of MCBAs reached those measured natively in the mammalian gut. MCBAs fed to mice entered enterohepatic circulation, in which liver and gallbladder concentrations varied depending on the conjugated amino acid. Quantifying MCBAs in human faecal samples showed that they reach concentrations equal to or greater than secondary and primary BAs and were reduced after bariatric surgery, thus supporting MCBAs as a significant component of the BA pool that can be altered by changes in gastrointestinal physiology. In conclusion, the inherent acyltransferase activity of BSH/T greatly diversifies BA chemistry, creating a set of previously underappreciated metabolites with the potential to affect the microbiome and human health.

Elucidation of the Active Form and Reaction Mechanism in Human Asparaginase Type III Using Multiscale Simulations

l-asparaginases catalyze the asparagine hydrolysis to aspartate. These enzymes play an important role in the treatment of acute lymphoblastic leukemia because these cells are unable to produce their own asparagine. Due to the immunogenic response and various side effects of enzymes of bacterial origin, many attempts have been made to replace these enzymes with mammalian enzymes such as human asparaginase type III (hASNaseIII). This study investigates the reaction mechanism of hASNaseIII through molecular dynamics simulations, quantum mechanics/molecular mechanics methods, and free energy calculations. Our simulations reveal that the dimeric form of the enzyme plays a vital role in stabilizing the substrate in the active site, despite the active site residues coming from a single protomer. Protomer-protomer interactions are essential to keep the enzyme in an active conformation. Our study of the reaction mechanism indicates that the self-cleavage process that generates an N-terminal residue (Thr168) is required to activate the enzyme. This residue acts as the nucleophile, attacking the electrophilic carbon of the substrate after a proton transfer from its hydroxyl group to the N-terminal amino group. The reaction mechanism proceeds with the formation of an acyl-enzyme complex and its hydrolysis, which turns out to be the rate-determining step. Our proposal of the enzymatic mechanism sheds light on the role of different active site residues and rationalizes the studies on mutations. The insights provided here about hASNaseIII activity could contribute to the comprehension of the disparities among different ASNases and might even guide the design of new variants with improved properties for acute lymphoblastic leukemia treatment.

Discovery of a spirocyclic 3-bromo-4,5-dihydroisoxazole covalent inhibitor of hGAPDH with antiproliferative activity against pancreatic cancer cells

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), a key glycolytic enzyme, plays a crucial role in the energy metabolism of cancer cells and has been proposed as a valuable target for the development of anticancer agents. Among a series of 5-substituted 3-bromo-4,5-dihydroisoxazole (BDHI) derivatives, we identified the spirocyclic compound 11, which is able to covalently inactivate recombinant human GAPDH (hGAPDH) with a faster reactivity than koningic acid, one of the most potent hGAPDH inhibitors known to date. Computational studies confirmed that conformational rigidification is crucial to stabilize the interaction of the inhibitor with the binding site, thus favoring the subsequent covalent bond formation. Investigation of intrinsic warhead reactivity at different pH disclosed the negligible reactivity of 11 with free thiols, highlighting its ability to selectively react with the activated cysteine of hGAPDH with respect to other sulfhydryl groups. Compound 11 strongly reduced cancer cell growth in four different pancreatic cancer cell lines and its antiproliferative activity correlated well with the intracellular inhibition of hGAPDH. Overall, our results qualify 11 as a potent hGAPDH covalent inhibitor with a moderate drug-like reactivity that could be further exploited to develop anticancer agents.

Role of Stereochemistry on the Biological Activity of Nature-Inspired 3-Br-Acivicin Isomers and Derivatives

Chiral natural compounds are often biosynthesized in an enantiomerically pure fashion, and stereochemistry plays a pivotal role in biological activity. Herein, we investigated the significance of chirality for nature-inspired 3-Br-acivicin (3-BA) and its derivatives. The three unnatural isomers of 3-BA and its ester and amide derivatives were prepared and characterized for their antimalarial activity. Only the (5S, αS) isomers displayed significant antiplasmodial activity, revealing that their uptake might be mediated by the L-amino acid transport system, which is known to mediate the acivicin membrane's permeability. In addition, we investigated the inhibitory activity towards Plasmodium falciparum glyceraldehyde 3-phosphate dehydrogenase (PfGAPDH) since it is involved in the multitarget mechanism of action of 3-BA. Molecular modeling has shed light on the structural and stereochemical requirements for an efficient interaction with PfGAPDH, leading to covalent irreversible binding and enzyme inactivation. While stereochemistry affects the target binding only for two subclasses (1a-d and 4a-d), it leads to significant differences in the antimalarial activity for all subclasses, suggesting that a stereoselective uptake might be responsible for the enhanced biological activity of the (5S, αS) isomers.

Characterization of the mechanism of bile salt hydrolase substrate specificity by experimental and computational analyses

Bile salt hydrolases (BSHs) are currently being investigated as target enzymes for metabolic regulators in humans and as growth promoters in farm animals. Understanding structural features underlying substrate specificity is necessary for inhibitor design. Here, we used a multidisciplinary workflow including mass spectrometry, mutagenesis, molecular dynamic simulations, machine learning, and crystallography to demonstrate substrate specificity in Lactobacillus salivarius BSH, the most abundant enzyme in human and farm animal intestines. We show the preference of substrates with a taurine head and a dehydroxylated sterol ring for hydrolysis. A regression model that correlates the relative rates of hydrolysis of various substrates in various enzyme mutants with the residue-substrate interaction energies guided the identification of structural determinants of substrate binding and specificity. In addition, we found T208 from another BSH protomer regulating the hydrolysis. The designed workflow can be used for fast and comprehensive characterization of enzymes with a broad range of substrates.

Mechanistic Modeling of Lys745 Sulfonylation in EGFR C797S Reveals Chemical Determinants for Inhibitor Activity and Discriminates Reversible from Irreversible Agents

Targeted covalent inhibitors hold promise for drug discovery, particularly for kinases. Targeting the catalytic lysine of epidermal growth factor receptor (EGFR) has attracted attention as a new strategy to overcome resistance due to the emergence of C797S mutation. Sulfonyl fluoride derivatives able to inhibit EGFR^{L858R/T790M/C797S} by sulfonylation of Lys745 have been reported. However, atomistic details of this process are still poorly understood. Here, we describe the mechanism of inhibition of an innovative class of compounds that covalently engage the catalytic lysine of EGFR, through a sulfur(VI) fluoride exchange (SuFEx) process, with the help of hybrid quantum mechanics/molecular mechanics (QM/MM) and path collective variables (PCVs) approaches. Our simulations identify the chemical determinants accounting for the irreversible activity of agents targeting Lys745 and provide hints for the further optimization of sulfonyl fluoride agents.

The Human Ntn-Hydrolase Superfamily: Structure, Functions and Perspectives

N-terminal nucleophile (Ntn)-hydrolases catalyze the cleavage of amide bonds in a variety of macromolecules, including the peptide bond in proteins, the amide bond in N-linked protein glycosylation, and the amide bond linking a fatty acid to sphingosine in complex sphingolipids. Ntn-hydrolases are all sharing two common hallmarks: Firstly, the enzymes are synthesized as inactive precursors that undergo auto-proteolytic self-activation, which, as a consequence, reveals the active site nucleophile at the newly formed N-terminus. Secondly, all Ntn-hydrolases share a structural consistent αββα-fold, notwithstanding the total lack of amino acid sequence homology. In humans, five subclasses of the Ntn-superfamily have been identified so far, comprising relevant members such as the catalytic active subunits of the proteasome or a number of lysosomal hydrolases, which are often associated with lysosomal storage diseases. This review gives an updated overview on the structural, functional, and (patho-)physiological characteristics of human Ntn-hydrolases, in particular.

Free-Energy Simulations Support a Lipophilic Binding Route for Melatonin Receptors

The effects of the neurohormone melatonin are mediated by the activation of the GPCRs MT₁ and MT₂ in a variety of tissues. Crystal structures suggest ligand access to the orthosteric binding site of MT₁ and MT₂ receptors through a lateral channel between transmembrane (TM) helices IV and V. We investigated the feasibility of this lipophilic entry route for 2-iodomelatonin, a nonselective agonist with a slower dissociation rate from the MT₂ receptor, applying enhanced sampling simulations and free-energy calculations. 2-Iodomelatonin unbinding was investigated with steered molecular dynamics simulations which revealed different trajectories passing through the gap between TM helices IV and V for both receptors. For one of these unbinding trajectories from the MT₁ receptor, an umbrella-sampling protocol with path-collective variables provided a calculated energy barrier consistent with the experimental dissociation rate. The side-chain flexibility of Tyr5.38 was significantly different in the two receptor subtypes, as assessed by metadynamics simulations, and during ligand unbinding it frequently assumes an open conformation in the MT₁ but not in the MT₂ receptor, favoring 2-iodomelatonin egress. Taken together, our simulations are consistent with the possibility that the gap between TM IV and V is a way of connecting the orthosteric binding site and the membrane core for lipophilic melatonin receptor ligands. Our simulations also suggest that the open state of Tyr5.38 generates a small pocket on the surface of MT₁ receptor, which could participate in the recognition of MT₁-selective ligands and may be exploited in the design of new selective compounds.

Functional and Phylogenetic Diversity of BSH and PVA Enzymes

Bile salt hydrolase (BSH) and penicillin V acylase (PVA) are related enzymes that are classified as choloylglycine hydrolases (CGH). BSH enzymes have attracted significant interest for their ability to modulate the composition of the bile acid pool, alter bile acid signaling events mediated by the host bile acid receptors FXR and TGR5 and influence cholesterol homeostasis in the host, while PVA enzymes have been widely utilised in an industrial capacity in the production of semi-synthetic antibiotics. The similarities between BSH and PVA enzymes suggest common evolution of these enzymes and shared mechanisms for substrate binding and catalysis. Here, we compare BSH and PVA through analysis of the distribution, phylogeny and biochemistry of these microbial enzymes. The development of new annotation approaches based upon functional enzyme analyses and the potential implications of BSH enzymes for host health are discussed.

Density-functional tight-binding: basic concepts and applications to molecules and clusters

The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.

Design and SAR Analysis of Covalent Inhibitors Driven by Hybrid QM/MM Simulations

Quantum mechanics/molecular mechanics (QM/MM) hybrid technique is emerging as a reliable computational method to investigate and characterize chemical reactions occurring in enzymes. From a drug discovery perspective, a thorough understanding of enzyme catalysis appears pivotal to assist the design of inhibitors able to covalently bind one of the residues belonging to the enzyme catalytic machinery. Thanks to the current advances in computer power, and the availability of more efficient algorithms for QM-based simulations, the use of QM/MM methodology is becoming a viable option in the field of covalent inhibitor design. In the present review, we summarized our experience in the field of QM/MM simulations applied to drug design problems which involved the optimization of agents working on two well-known drug targets, namely fatty acid amide hydrolase (FAAH) and epidermal growth factor receptor (EGFR). In this context, QM/MM simulations gave valuable information in terms of geometry (i.e., of transition states and metastable intermediates) and reaction energetics that allowed to correctly predict inhibitor binding orientation and substituent effect on enzyme inhibition. What is more, enzyme reaction modelling with QM/MM provided insights that were translated into the synthesis of new covalent inhibitor featured by a unique combination of intrinsic reactivity, on-target activity, and selectivity.

Design and synthesis of cyanamides as potent and selective N-acylethanolamine acid amidase inhibitors

N-acylethanolamine acid amidase (NAAA) inhibition represents an exciting novel approach to treat inflammation and pain. NAAA is a cysteine amidase which preferentially hydrolyzes the endogenous biolipids palmitoylethanolamide (PEA) and oleoylethanolamide (OEA). PEA is an endogenous agonist of the nuclear peroxisome proliferator-activated receptor-α (PPAR-α), which is a key regulator of inflammation and pain. Thus, blocking the degradation of PEA with NAAA inhibitors results in augmentation of the PEA/PPAR-α signaling pathway and regulation of inflammatory and pain processes. We have prepared a new series of NAAA inhibitors exploring the azetidine-nitrile (cyanamide) pharmacophore that led to the discovery of highly potent and selective compounds. Key analogs demonstrated single-digit nanomolar potency for hNAAA and showed >100-fold selectivity against serine hydrolases FAAH, MGL and ABHD6, and cysteine protease cathepsin K. Additionally, we have identified potent and selective dual NAAA-FAAH inhibitors to investigate a potential synergism between two distinct anti-inflammatory molecular pathways, the PEA/PPAR-α anti-inflammatory signaling pathway,^1-4 and the cannabinoid receptors CB1 and CB2 pathways which are known for their antiinflammatory and antinociceptive properties.^5-8 Our ligand design strategy followed a traditional structure-activity relationship (SAR) approach and was supported by molecular modeling studies of reported X-ray structures of hNAAA. Several inhibitors were evaluated in stability assays and demonstrated very good plasma stability (t_1/2 > 2 h; human and rodents). The disclosed cyanamides represent promising new pharmacological tools to investigate the potential role of NAAA inhibitors and dual NAAA-FAAH inhibitors as therapeutic agents for the treatment of inflammation and pain.

The complex structure of bile salt hydrolase from Lactobacillus salivarius reveals the structural basis of substrate specificity

The gut bacterial bile salt hydrolase (BSH) plays a critical role in host lipid metabolism and energy harvest. Therefore, BSH is a promising microbiome target to develop new therapies to regulate obesity in humans and novel non-antibiotic growth promoters for food animals. We previously reported the 1.90 Å apo crystal structure of BSH from Lactobacillus salivarius (lsBSH). In this study, we soaked the lsBSH crystal with glycocholic acid (GCA), a substrate, and obtained a 2.10 Å structure containing complex of lsBSH bound to GCA and cholic acid (CA), a product. The substrate/product sits in the water-exposed cavity molded by Loops 2 and 3. While the glycine moiety of GCA is exposed into a highly polar pocket, the sterane core of GCA is stabilized by aromatic and hydrophobic interactions. Comparison of product binding with BSH from Clostridium perfringenes reveals a distinct orientation of the sterane core in the binding site. The stability of the substrate-lsBSH complex and the putative catalytic mechanism were explored with molecular dynamics simulations. Site-directed mutagenesis of lsBSH demonstrated that Cys2 and Asn171 are critical for enzymatic activity, while Tyr24, Phe65 and Gln257 contribute to the substrate specificity. Together, this study provides structural insights into BSH-substrate interaction, the mechanism of catalysis and substrate specificity, which facilitate rational design of BSH inhibitors.

Covalent Inhibitors of Plasmodium falciparum Glyceraldehyde 3-Phosphate Dehydrogenase with Antimalarial Activity in Vitro

Covalent inhibitors of PfGAPDH characterized by a 3-bromoisoxazoline warhead were developed, and their mode of interaction with the target enzyme was interpreted by means of molecular modeling studies: some of them displayed a submicromolar antiplasmodial activity against both chloroquine sensitive and resistant strains of Plasmodium falciparum, with good selectivity indices.

Bile salt hydrolases: Structure and function, substrate preference, and inhibitor development

The worldwide trend of limiting the use of antibiotic growth promoters (AGPs) in animal production creates challenges for the animal feed industry, thus necessitating the development of effective non-antibiotic alternatives to improve animal performance. Increasing evidence has shown that the growth-promoting effect of AGPs is highly correlated with the reduced activity of bile salt hydrolase (BSH, EC 3.5.1.24), an intestinal bacteria-producing enzyme that has a negative impact on host fat digestion and energy harvest. Therefore, BSH inhibitors may become novel, attractive alternatives to AGPs. Detailed knowledge of BSH substrate preferences and the wealth of structural data on BSHs provide a solid foundation for rationally tailored BSH inhibitor design. This review focuses on the relationship between structure and function of BSHs based on the crystal structure, kinetic data, molecular docking and comparative structural analyses. The molecular basis for BSH substrate recognition is also discussed. Finally, recent advances and future prospectives in the development of potent, safe, and cost-effective BSH inhibitors are described.

A Comprehensive Genome Survey Provides Novel Insights into Bile Salt Hydrolase (BSH) in Lactobacillaceae

Bile salt hydrolase (BSH) is a well-known enzyme that has been commonly characterized in probiotic bacteria, as it has cholesterol-lowering effects. However, its molecular investigations are scarce. Here, we build a local database of BSH sequences from Lactobacillaceae (BSH⁻SDL), and phylogenetic analysis and homology searches were employed to elucidate their comparability and distinctiveness among species. Evolutionary study demonstrates that BSH sequences in BSH⁻SDL are divided into five groups, named BSH A, B, C, D and E here, which can be the genetic basis for BSH classification and nomenclature. Sequence analysis suggests the differences between BSH-active and BSH-inactive proteins clearly, especially on site 82. In addition, a total of 551 BSHs from 107 species are identified from 451 genomes of 158 Lactobacillaceae species. Interestingly, those bacteria carrying various copies of BSH A or B can be predicted to be potential cholesterol-lowering probiotics, based on the results of phylogenetic analysis and the subtypes that those previously reported BSH-active probiotics possess. In summary, this study elaborates the molecular basis of BSH in Lactobacillaceae systematically, and provides a novel methodology as well as a consistent standard for the identification of the BSH subtype. We believe that high-throughput screening can be efficiently applied to the selection of promising candidate BSH-active probiotics, which will advance the development of healthcare products in cholesterol metabolism.

L718Q mutant EGFR escapes covalent inhibition by stabilizing a non-reactive conformation of the lung cancer drug osimertinib

Impact of L718Q mutation on the inhibitory activity of osimertinib on EGFR revealed by free-energy simulations.

Osimertinib is a third-generation inhibitor approved for the treatment of non-small cell lung cancer. It overcomes resistance to first-generation inhibitors by incorporating an acrylamide group which alkylates Cys797 of EGFR T790M. The mutation of a residue in the P-loop (L718Q) was shown to cause resistance to osimertinib, but the molecular mechanism of this process is unknown. Here, we investigated the inhibitory process for EGFR T790M (susceptible to osimertinib) and EGFR T790M/L718Q (resistant to osimertinib), by modelling the chemical step (i.e., alkylation of Cys797) using QM/MM simulations and the recognition step by MD simulations coupled with free-energy calculations. The calculations indicate that L718Q has a negligible impact on both the activation energy for Cys797 alkylation and the free-energy of binding for the formation of the non-covalent complex. The results show that Gln718 affects the conformational space of the EGFR–osimertinib complex, stabilizing a conformation of acrylamide which prevents reaction with Cys797.

Structure and function of a highly active Bile Salt Hydrolase (BSH) from Enterococcus faecalis and post-translational processing of BSH enzymes

Bile Salt Hydrolase (BSH), a member of Cholylglycine hydrolase family, catalyzes the de-conjugation of bile acids and is evolutionarily related to penicillin V acylase (PVA) that hydrolyses a different substrate such as penicillin V. We report the three-dimensional structure of a BSH enzyme from the Gram-positive bacteria Enterococcus faecalis (EfBSH) which has manifold higher hydrolase activity compared to other known BSHs and displays unique allosteric catalytic property. The structural analysis revealed reduced secondary structure content compared to other known BSH structures, particularly devoid of an anti-parallel β-sheet in the assembly loop and part of a β-strand is converted to increase the length of a substrate binding loop 2. The analysis of the substrate binding pocket showed reduced volume owing to altered loop conformations and increased hydrophobicity contributed by a higher ratio of hydrophobic to hydrophilic groups present. The aromatic residues F18, Y20 and F65 participate in substrate binding. Thus, their mutation affects enzyme activity. Docking and Molecular Dynamics simulation studies showed effective polar complementarity present for the three hydroxyl (-OH) groups of GCA substrate in the binding site contributing to higher substrate specificity and efficient catalysis. These are unique features characteristics of this BSH enzyme and thought to contribute to its higher activity and specificity towards bile salts as well as allosteric effects. Further, mechanism of autocatalytic processing of Cholylglycine Hydrolases by the excision of an N-terminal Pre-peptide was examined by inserting different N-terminal pre-peptides in EfBSH sequence. The results suggest that two serine residues next to nucleophile cysteine are essential for autocalytic processing to remove precursor peptide. Since pre-peptide is absent in EfBSH the mutation of these serines is tolerated. This suggests that an evolution-mediated subordination of the pre-peptide excision site resulted in loss of pre-peptide in EfBSH and other related Cholylglycine hydrolases.

Competition between two cysteines in covalent binding of biliverdin to phytochrome domains

In this work, we disclose a mechanism of competing chemical reactions of protein assembly for a bacterial phytochrome using modern methods of molecular modeling.

Mechanistic Insights into the Reaction of Chlorination of Tryptophan Catalyzed by Tryptophan 7-Halogenase

Tryptophan 7-halogenase catalyzes chlorination of free tryptophan to 7-chlorotryptophan, which is the first step in the antibiotic pyrrolnitrin biosynthesis. Many biologically and pharmaceutically active natural products contain chlorine and thus, an understanding of the mechanism of its introduction into organic molecules is important. Whilst enzyme-catalyzed chlorination is accomplished with ease, it remains a difficult task for the chemists. Therefore, utilizing enzymes in the synthesis of chlorinated organic compounds is important, and providing atomistic mechanistic insights about the reaction mechanism of tryptophan 7-halogenase is vital and timely. In this work, we examined a mechanism for the reaction of tryptophan chlorination, performed by tryptophan 7-halogenase, by calculating potential energy and free energy surfaces using two different Combined Quantum Mechanical/Molecular Mechanical (QM/MM) methods both employing Density Functional Theory (DFT) for the QM region. Both computational strategies agree on the nature of the rate-limiting step and provided close results for the reaction barriers of the two reaction steps. The calculations for both the potential energy and the free energy profiles showed very similar geometric features and hydrogen bonding interactions for the characterized stationary points.

The Conformational Change in Elongation Factor Tu Involves Separation of Its Domains

Elongation factor Tu (EF-Tu) is a highly conserved GTPase that is responsible for supplying the aminoacylated tRNA to the ribosome. Upon binding to the ribosome, EF-Tu undergoes GTP hydrolysis, which drives a major conformational change, triggering the release of aminoacylated tRNA to the ribosome. Using a combination of molecular simulation techniques, we studied the transition between the pre- and post-hydrolysis structures through two distinct pathways. We show that the transition free energy is minimal along a non-intuitive pathway that involves "separation" of the GTP binding domain (domain 1) from the OB folds (domains 2 and 3), followed by domain 1 rotation, and, eventually, locking the EF-Tu conformation in the post-hydrolysis state. The domain separation also leads to a slight extension of the linker connecting domain 1 to domain 2. Using docking tools and correlation-based analysis, we identified and characterized the EF-Tu conformations that release the tRNA. These calculations suggest that EF-Tu can release the tRNA before the domains separate and after domain 1 rotates by 25°. We also examined the EF-Tu conformations in the context of the ribosome. Given the high degrees of sequence similarity with other translational GTPases, we predict a similar separation mechanism is followed.

Pharmacophore Identification and Scaffold Exploration to Discover Novel, Potent, and Chemically Stable Inhibitors of Acid Ceramidase in Melanoma Cells

Acid ceramidase (AC) hydrolyzes ceramides, which are central lipid messengers for metabolism and signaling of sphingolipids. A growing body of evidence links deregulation of sphingolipids to several diseases, including cancer. Indeed, AC expression is abnormally high in melanoma cells. AC inhibition may thus be key to treating malignant melanoma. Here, we have used a systematic scaffold exploration to design a general pharmacophore for AC inhibition. This pharmacophore comprises a 6 + 5 fused ring heterocycle linked to an aliphatic substituent via a urea moiety. We have thus identified the novel benzimidazole derivatives 10, 21, 27, and 30, which are highly potent AC inhibitors. Their chemical and metabolic stabilities are comparable or superior to those of previously reported AC inhibitors. Moreover, they are potent against endogenous AC in intact melanoma cells. These novel inhibitors merit further characterization and can serve as a promising starting point for the discovery of new antimelanoma therapeutics.

Lysosomal glycosphingolipid catabolism by acid ceramidase: formation of glycosphingoid bases during deficiency of glycosidases

Glycosphingoid bases are elevated in inherited lysosomal storage disorders with deficient activity of glycosphingolipid catabolizing glycosidases. We investigated the molecular basis of the formation of glucosylsphingosine and globotriaosylsphingosine during deficiency of glucocerebrosidase (Gaucher disease) and α-galactosidase A (Fabry disease). Independent genetic and pharmacological evidence is presented pointing to an active role of acid ceramidase in both processes through deacylation of lysosomal glycosphingolipids. The potential pathophysiological relevance of elevated glycosphingoid bases generated through this alternative metabolism in patients suffering from lysosomal glycosidase defects is discussed.

Structural analysis of a penicillin V acylase from Pectobacterium atrosepticum confirms the importance of two Trp residues for activity and specificity

Penicillin V acylases (PVA) catalyze the deacylation of the beta-lactam antibiotic phenoxymethylpenicillin (Pen V). They are members of the Ntn hydrolase family and possess an N-terminal cysteine as the main catalytic nucleophile residue. They form the evolutionarily related cholylglycine hydrolase (CGH) group which includes bile salt hydrolases (BSH) responsible for bile deconjugation. Even though a few PVA and BSH structures have been reported, no structure of a functional PVA from Gram-negative bacteria is available. Here, we report the crystal structure of a highly active PVA from Gram-negative Pectobacterium atrosepticum (PaPVA) at 2.5Å resolution. Structural comparison with PVAs from Gram-positive bacteria revealed that PaPVA had a distinctive tetrameric structure and active site organization. In addition, mutagenesis of key active site residues and biochemical characterization of the resultant variants elucidated the role of these residues in substrate binding and catalysis. The importance of residue Trp23 and Trp87 side chains in binding and correct positioning of Pen V by PVAs was confirmed using mutagenesis and substrate docking with a 15ns molecular dynamics simulation. These results establish the unique nature of Gram-negative CGHs and necessitate further research about their substrate spectrum.

Periplasmic Export of Bile Salt Hydrolase in Escherichia coli by the Twin-Arginine Signal Peptides

Bile salt hydrolase (BSH, EC 3.5.1.24) is considered as an ideal way with lower cost and less side effects to release the risk of coronary heart disease caused by hypercholesterolemia. As bile salt hydrolase from Lactobacillus plantarum BBE7 could not be efficiently exported by PelB signal peptide of the general secretory (Sec) pathway, three twin-arginine signal peptides from twin-arginine translocation (Tat) pathway were synthesized, fused with bsh gene, inserted into expression vectors pET-20b(+) and pET-22b(+), and transformed into four different Escherichia coli hosts, respectively. Among the 24 recombinant bacteria obtained, E. coli BL21 (DE3) pLysS (pET-20b(+)-dmsA-bsh) showed the highest BSH activity in periplasmic fraction, which was further increased to 1.21 ± 0.03 U/mL by orthogonal experimental design. And, signal peptide dimethyl sulfoxide reductase subunit DmsA (DMSA) had the best activity of exported BSH. More importantly, the presence of BSH in the periplasm had proven to be caused by the export rather than cell leakage. For the first time, we report the periplasmic expression of BSH by signal peptides from the Tat pathway. This will lay a solid foundation for the purification and biochemical characterization of BSH from the supernatant, and strategies adopted here could be used for the periplasmic expression of other proteins in E. coli.

Quantum mechanics/molecular mechanics modeling of covalent addition between EGFR-cysteine 797 and N-(4-anilinoquinazolin-6-yl) acrylamide

Irreversible epidermal growth factor receptor (EGFR) inhibitors can circumvent resistance to first-generation ATP-competitive inhibitors in the treatment of nonsmall-cell lung cancer. They covalently bind a noncatalytic cysteine (Cys797) at the surface of EGFR active site by an acrylamide warhead. Herein, we used a hybrid quantum mechanics/molecular mechanics (QM/MM) potential in combination with umbrella sampling in the path-collective variable space to investigate the mechanism of alkylation of Cys797 by the prototypical covalent inhibitor N-(4-anilinoquinazolin-6-yl) acrylamide. Calculations show that Cys797 reacts with the acrylamide group of the inhibitor through a direct addition mechanism, with Asp800 acting as a general base/general acid in distinct steps of the reaction. The obtained reaction free energy is negative (ΔA = -12 kcal/mol) consistent with the spontaneous and irreversible alkylation of Cys797 by N-(4-anilinoquinazolin-6-yl) acrylamide. Our calculations identify desolvation of Cys797 thiolate anion as a key step of the alkylation process, indicating that changes in the intrinsic reactivity of the acrylamide would have only a minor impact on the inhibitor potency.

Penicillin acylases revisited: importance beyond their industrial utility

It is of great importance to study the physiological roles of enzymes in nature; however, in some cases, it is not easily apparent. Penicillin acylases are pharmaceutically important enzymes that cleave the acyl side chains of penicillins, thus paving the way for production of newer semi-synthetic antibiotics. They are classified according to the type of penicillin (G or V) that they preferentially hydrolyze. Penicillin acylases are also used in the resolution of racemic mixtures and peptide synthesis. However, it is rather unfortunate that the focus on the use of penicillin acylases for industrial applications has stolen the spotlight from the study of the importance of these enzymes in natural metabolism. The penicillin acylases, so far characterized from different organisms, show differences in their structural nature and substrate spectrum. These enzymes are also closely related to the bacterial signalling phenomenon, quorum sensing, as detailed in this review. This review details studies on biochemical and structural characteristics of recently discovered penicillin acylases. We also attempt to organize the available insights into the possible in vivo role of penicillin acylases and related enzymes and emphasize the need to refocus research efforts in this direction.

A reaction path study of the catalysis and inhibition of the Bacillus anthracis CapD γ-glutamyl transpeptidase

The CapD enzyme of Bacillus anthracis is a γ-glutamyl transpeptidase from the N-terminal nucleophile hydrolase superfamily that covalently anchors the poly-γ-D-glutamic acid (pDGA) capsule to the peptidoglycan. The capsule hinders phagocytosis of B. anthracis by host cells and is essential for virulence. The role CapD plays in capsule anchoring and remodeling makes the enzyme a promising target for anthrax medical countermeasures. Although the structure of CapD is known, and a covalent inhibitor, capsidin, has been identified, the mechanisms of CapD catalysis and inhibition are poorly understood. Here, we used a computational approach to map out the reaction steps involved in CapD catalysis and inhibition. We found that the rate-limiting step of either CapD catalysis or inhibition was a concerted asynchronous formation of the tetrahedral intermediate with a barrier of 22-23 kcal/mol. However, the mechanisms of these reactions differed for the two amides. The formation of the tetrahedral intermediate with pDGA was substrate-assisted with two proton transfers. In contrast, capsidin formed the tetrahedral intermediate in a conventional way with one proton transfer. Interestingly, capsidin coupled a conformational change in the catalytic residue of the tetrahedral intermediate to stretching of the scissile amide bond. Furthermore, capsidin took advantage of iminol-amide tautomerism of its diacetamide moiety to convert the tetrahedral intermediate to the acetylated CapD. As evidence of the promiscuous nature of CapD, the enzyme cleaved the amide bond of capsidin by attacking it on the opposite side compared to pDGA.

Computational insights into function and inhibition of fatty acid amide hydrolase

The Fatty Acid Amide Hydrolase (FAAH) enzyme is a membrane-bound serine hydrolase responsible for the deactivating hydrolysis of a family of naturally occurring fatty acid amides. FAAH is a critical enzyme of the endocannabinoid system, being mainly responsible for regulating the level of its main cannabinoid substrate anandamide. For this reason, pharmacological inhibition of FAAH, which increases the level of endogenous anandamide, is a promising strategy to cure a variety of diseases including pain, inflammation, and cancer. Much structural, mutagenesis, and kinetic data on FAAH has been generated over the last couple of decades. This has prompted several informative computational investigations to elucidate, at the atomic-level, mechanistic details on catalysis and inhibition of this pharmaceutically relevant enzyme. Here, we review how these computational studies - based on classical molecular dynamics, full quantum mechanics, and hybrid QM/MM methods - have clarified the binding and reactivity of some relevant substrates and inhibitors of FAAH. We also discuss the experimental implications of these computational insights, which have provided a thoughtful elucidation of the complex physical and chemical steps of the enzymatic mechanism of FAAH. Finally, we discuss how computations have been helpful for building structure-activity relationships of potent FAAH inhibitors.

An improved method for specificity annotation shows a distinct evolutionary divergence among the microbial enzymes of the cholylglycine hydrolase family

Bile salt hydrolases (BSHs) are gut microbial enzymes that play a significant role in the bile acid modification pathway. Penicillin V acylases (PVAs) are enzymes produced by environmental microbes, having a possible role in pathogenesis or scavenging of phenolic compounds in their microbial habitats. The correct annotation of such physiologically and industrially important enzymes is thus vital. The current methods relying solely on sequence homology do not always provide accurate annotations for these two members of the cholylglycine hydrolase (CGH) family as BSH/PVA enzymes. Here, we present an improved method [binding site similarity (BSS)-based scoring system] for the correct annotation of the CGH family members as BSH/PVA enzymes, which along with the phylogenetic information incorporates the substrate specificity as well as the binding site information. The BSS scoring system was developed through the analysis of the binding sites and binding modes of the available BSH/PVA structures with substrates glycocholic acid and penicillin V. The 198 sequences in the dataset were then annotated accurately using BSS scores as BSH/PVA enzymes. The dataset presented contained sequences from Gram-positive bacteria, Gram-negative bacteria and archaea. The clustering obtained for the dataset using the method described above showed a clear distinction in annotation of Gram-positive bacteria and Gram-negative bacteria. Based on this clustering and a detailed analysis of the sequences of the CGH family in the dataset, we could infer that the CGH genes might have evolved in accordance with the hypothesis stating the evolution of diderms and archaea from the monoderms.

Predicting the reactivity of nitrile-carrying compounds with cysteine: a combined computational and experimental study

Here, we report on a mechanistic investigation based on DFT calculations and kinetic measures aimed at determining the energetics related to the cysteine nucleophilic attack on nitrile-carrying compounds. Activation energies were found to correlate well with experimental kinetic measures of reactivity with cysteine in phosphate buffer. The agreement between computations and experiments points to this DFT-based approach as a tool for predicting both nitrile reactivity toward cysteines and the toxicity of nitriles as electrophile agents.

Advances in the discovery of N-acylethanolamine acid amidase inhibitors

N-Acylethanolamine acid amidase (NAAA) is a cysteine amidase that hydrolyzes saturated or monounsaturated fatty acid ethanolamides, such as palmitoylethanolamide (PEA) and oleoylethanolamide (OEA). PEA has been shown to exert analgesic and anti-inflammatory effects by engaging peroxisome proliferator-activated receptor-α. Like other fatty acid ethanolamides, PEA is not stored in cells, but produced on demand from cell membrane precursors, and its actions are terminated by intracellular hydrolysis by either fatty acid amide hydrolase or NAAA. Endogenous levels of PEA and OEA have been shown to decrease during inflammation. Modulation of the tissue levels of PEA by inhibition of enzymes responsible for the breakdown of this lipid mediator may represent therefore a new therapeutic strategy for the treatment of pain and inflammation. While a large number of inhibitors of fatty acid amide hydrolase have been discovered, few compounds have been reported to inhibit NAAA activity. Here, we describe the most representative NAAA inhibitors and briefly highlight their pharmacological profile. A recent study has shown that a NAAA inhibitor attenuated heat hyperalgesia and mechanical allodynia caused by local inflammation or nerve damage in animal models of pain and inflammation. This finding encourages further exploration of the pharmacology of NAAA inhibitors.

Novel insights into eukaryotic γ-glutamyltranspeptidase 1 from the crystal structure of the glutamate-bound human enzyme

The enzyme γ-glutamyltranspeptidase 1 (GGT1) is a conserved member of the N-terminal nucleophile hydrolase family that cleaves the γ-glutamyl bond of glutathione and other γ-glutamyl compounds. In animals, GGT1 is expressed on the surface of the cell and has critical roles in maintaining cysteine levels in the body and regulating intracellular redox status. Expression of GGT1 has been implicated as a potentiator of asthma, cardiovascular disease, and cancer. The rational design of effective inhibitors of human GGT1 (hGGT1) has been delayed by the lack of a reliable structural model. The available crystal structures of several bacterial GGTs have been of limited use due to differences in the catalytic behavior of bacterial and mammalian GGTs. We report the high resolution (1.67 Å) crystal structure of glutamate-bound hGGT1, the first of any eukaryotic GGT. Comparisons of the active site architecture of hGGT1 with those of its bacterial orthologs highlight key differences in the residues responsible for substrate binding, including a bimodal switch in the orientation of the catalytic nucleophile (Thr-381) that is unique to the human enzyme. Compared with several bacterial counterparts, the lid loop in the crystal structure of hGGT1 adopts an open conformation that allows greater access to the active site. The hGGT1 structure also revealed tightly bound chlorides near the catalytic residue that may contribute to catalytic activity. These are absent in the bacterial GGTs. These differences between bacterial and mammalian GGTs and the new structural data will accelerate the development of new therapies for GGT1-dependent diseases.

String method for calculation of minimum free-energy paths in Cartesian space in freely-tumbling systems

The string method is a molecular-simulation technique that aims to calculate the minimum free-energy path of a chemical reaction or conformational transition, in the space of a pre-defined set of reaction coordinates that is typically highly dimensional. Any descriptor may be used as a reaction coordinate, but arguably the Cartesian coordinates of the atoms involved are the most unprejudiced and intuitive choice. Cartesian coordinates, however, present a non-trivial problem, in that they are not invariant to rigid-body molecular rotations and translations, which ideally ought to be unrestricted in the simulations. To overcome this difficulty, we reformulate the framework of the string method to integrate an on-the-fly structural-alignment algorithm. This approach, referred to as SOMA (String method with Optimal Molecular Alignment), enables the use of Cartesian reaction coordinates in freely tumbling molecular systems. In addition, this scheme permits the dissection of the free-energy change along the most probable path into individual atomic contributions, thus revealing the dominant mechanism of the simulated process. This detailed analysis also provides a physically-meaningful criterion to coarse-grain the representation of the path. To demonstrate the accuracy of the method we analyze the isomerization of the alanine dipeptide in vacuum and the chair-to-inverted-chair transition of β-D mannose in explicit water. Notwithstanding the simplicity of these systems, the SOMA approach reveals novel insights into the atomic mechanism of these isomerizations. In both cases, we find that the dynamics and the energetics of these processes are controlled by interactions involving only a handful of atoms in each molecule. Consistent with this result, we show that a coarse-grained SOMA calculation defined in terms of these subsets of atoms yields nearidentical minimum free-energy paths and committor distributions to those obtained via a highly-dimensional string.