Purification and biochemical characterization of the Mycobacterium tuberculosis beta-ketoacyl-acyl carrier protein synthases KasA and KasB

Mycolic acids are vital components of the Mycobacterium tuberculosis cell wall, and enzymes involved in their formation represent attractive targets for the discovery of novel anti-tuberculosis agents. Biosynthesis of the fatty acyl chains of mycolic acids involves two fatty acid synthetic systems, the multifunctional polypeptide fatty acid synthase I (FASI), which performs de novo fatty acid synthesis, and the dissociated FASII system, which consists of monofunctional enzymes, and acyl carrier protein (ACP) and elongates FASI products to long chain mycolic acid precursors. In this study, we present the initial characterization of purified KasA and KasB, two beta-ketoacyl-ACP synthase (KAS) enzymes of the M. tuberculosis FASII system. KasA and KasB were expressed in E. coli and purified by affinity chromatography. Both enzymes showed activity typical of bacterial KASs, condensing an acyl-ACP with malonyl-ACP. Consistent with the proposed role of FASII in mycolic acid synthesis, analysis of various acyl-ACP substrates indicated KasA and KasB had higher specificity for long chain acyl-ACPs containing at least 16 carbons. Activity of KasA and KasB increased with use of M. tuberculosis AcpM, suggesting that structural differences between AcpM and E. coli ACP may affect their recognition by the enzymes. Both enzymes were sensitive to KAS inhibitors cerulenin and thiolactomycin. These results represent important steps in characterizing KasA and KasB as targets for antimycobacterial drug discovery.     

A novel interaction linking the FAS-II and phthiocerol dimycocerosate (PDIM) biosynthetic pathways

The fatty acid biosynthesis (FAS-II) pathway in Mycobacterium tuberculosis generates long chain fatty acids that serve as the precursors to mycolic acids, essential components of the mycobacterial cell wall. Enzymes in the FAS-II pathway are thought to form one or more noncovalent multi-enzyme complexes within the cell, and a bacterial two-hybrid screen was used to search for missing components of the pathway and to furnish additional data on interactions involving these enzymes in vivo. Using the FAS-II beta-ketoacyl synthase, KasA, as bait, an extensive bacterial two-hybrid screen of a M. tuberculosis genome fragment library unexpectedly revealed a novel interaction between KasA and PpsB as well as PpsD, two polyketide modules involved in the biosynthesis of the virulence lipid phthiocerol dimycocerosate (PDIM). Sequence analysis revealed that KasA interacts with PpsB and PpsD in the region of the acyl carrier domain of each protein, raising the possibility that lipids could be transferred between the FAS-II and PDIM biosynthetic pathways. Subsequent studies utilizing purified proteins and radiolabeled lipids revealed that fatty acids loaded onto PpsB were transferred to KasA and also incorporated into long chain fatty acids synthesized using a Mycobacterium smegmatis lysate. These data suggest that in addition to producing PDIMs, the growing phthiocerol product can also be shuttled into the FAS-II pathway via KasA as an entry point for further elongation. Interactions between these biosynthetic pathways may exist as a simple means to increase mycobacterial lipid diversity, enhancing functionality and the overall complexity of the cell wall.     

Computational delineation of the catalytic step of a high‐fidelity DNA polymerase

The Bacillus fragment, belonging to a class of high‐fidelity polymerases, demonstrates high processivity (adding ～115 bases per DNA binding event) and exceptional accuracy (1 error in 10⁶ nucleotide incorporations) during DNA replication. We present analysis of structural rearrangements and energetics just before and during the chemical step (phosphodiester bond formation) using a combination of classical molecular dynamics, mixed quantum mechanics molecular mechanics simulations, and free energy computations. We find that the reaction is associative, proceeding via the two‐metal‐ion mechanism, and requiring the proton on the terminal primer O3′ to transfer to the pyrophosphate tail of the incoming nucleotide before the formation of the pentacovalent transition state. Different protonation states for key active site residues direct the system to alternative pathways of catalysis and we estimate a free energy barrier of ～12 kcal/mol for the chemical step. We propose that the protonation of a highly conserved catalytic aspartic acid residue is essential for the high processivity demonstrated by the enzyme and suggest that global motions could be part of the reaction free energy landscape.     

Substrate recognition by β-ketoacyl-ACP synthases

β-Ketoacyl-ACP synthase (KAS) enzymes catalyze Claisen condensation reactions in the fatty acid biosynthesis pathway. These reactions follow a ping-pong mechanism in which a donor substrate acylates the active site cysteine residue after which the acyl group is condensed with the malonyl-ACP acceptor substrate to form a β-ketoacyl-ACP. In the priming KASIII enzymes the donor substrate is an acyl-CoA while in the elongating KASI and KASII enzymes the donor is an acyl-ACP. Although the KASIII enzyme in Escherichia coli (ecFabH) is essential, the corresponding enzyme in Mycobacterium tuberculosis (mtFabH) is not, suggesting that the KASI or II enzyme in M. tuberculosis (KasA or KasB, respectively) must be able to accept a CoA donor substrate. Since KasA is essential, the substrate specificity of this KASI enzyme has been explored using substrates based on phosphopantetheine, CoA, ACP, and AcpM peptide mimics. This analysis has been extended to the KASI and KASII enzymes from E. coli (ecFabB and ecFabF) where we show that a 14-residue malonyl-phosphopantetheine peptide can efficiently replace malonyl-ecACP as the acceptor substrate in the ecFabF reaction. While ecFabF is able to catalyze the condensation reaction when CoA is the carrier for both substrates, the KASI enzymes ecFabB and KasA have an absolute requirement for an ACP substrate as the acyl donor. Provided that this requirement is met, variation in the acceptor carrier substrate has little impact on the k(cat)/K(m) for the KASI reaction. For the KASI enzymes we propose that the binding of ecACP (AcpM) results in a conformational change that leads to an open form of the enzyme to which the malonyl acceptor substrate binds. Finally, the substrate inhibition observed when palmitoyl-CoA is the donor substrate for the KasA reaction has implications for the importance of mtFabH in the mycobacterial FASII pathway.     

Protonation states of methionine aminopeptidase and their relevance for inhibitor binding and catalytic activity

We have performed a computational study of different protomeric states of the methionine aminopeptidase active site using a combined quantum-mechanical/molecular mechanical simulation approach. The aim of this study was to clarify the native protonation state of the enzyme, which is needed for the development of novel irreversible inhibitors that can possibly be used as antiangiogenic and antibiotic drugs by virtual screening and other drug design methods. The results of the simulations indicated that two protonation states are possible without disturbing the overall geometry of the active site. We then verified experimentally the presence of the two protonation states by studying the substrate hydrolysis and inhibitor binding reactions at different pH values and come to the conclusion that one of the protomeric states is relevant for inhibitor binding, whereas the other is relevant for substrate hydrolysis. This result has implications for the development of other inhibitors of this class of enzymes and adds a new perspective to the pharmacological properties of the antiangiogenic drug fumagillin, which is an irreversible inhibitor of the human methionine aminopeptidase type II.     

Isoniazid affects multiple components of the type II fatty acid synthase system of Mycobacterium tuberculosis

Genetic and biochemical evidence has implicated two different target enzymes for isoniazid (INH) within the unique type II fatty acid synthase (FAS) system involved in the production of mycolic acids. These two components are an enoyl acyl carrier protein (ACP) reductase, InhA, and a beta-ketoacyl-ACP synthase, KasA. We compared the consequences of INH treatment of Mycobacterium tuberculosis (MTB) with two inhibitors having well-defined targets: triclosan (TRC), which inhibits InhA; and thiolactomycin (TLM), which inhibits KasA. INH and TLM, but not TRC, upregulate the expression of an operon containing five FAS II components, including kasA and acpM. Although all three compounds inhibit mycolic acid synthesis, treatment with INH and TLM, but not with TRC, results in the accumulation of ACP-bound lipid precursors to mycolic acids that were 26 carbons long and fully saturated. TLM-resistant mutants of MTB were more cross-resistant to INH than TRC-resistant mutants. Overexpression of KasA conferred more resistance to TLM and INH than to TRC. Overexpression of InhA conferred more resistance to TRC than to INH and TLM. Co-overexpression of both InhA and KasA resulted in strongly enhanced levels of INH resistance, in addition to cross-resistance to both TLM and TRC. These results suggest that these components of the FAS II complex are not independently regulated and that alterations in the expression level of InhA affect expression levels of KasA. Nonetheless, INH appeared to resemble TLM more closely in overall mode of action, and KasA levels appeared to be tightly correlated with INH sensitivity.     

Chromophore protonation state controls photoswitching of the fluoroprotein asFP595

Fluorescent proteins have been widely used as genetically encodable fusion tags for biological imaging. Recently, a new class of fluorescent proteins was discovered that can be reversibly light-switched between a fluorescent and a non-fluorescent state. Such proteins can not only provide nanoscale resolution in far-field fluorescence optical microscopy much below the diffraction limit, but also hold promise for other nanotechnological applications, such as optical data storage. To systematically exploit the potential of such photoswitchable proteins and to enable rational improvements to their properties requires a detailed understanding of the molecular switching mechanism, which is currently unknown. Here, we have studied the photoswitching mechanism of the reversibly switchable fluoroprotein asFP595 at the atomic level by multiconfigurational ab initio (CASSCF) calculations and QM/MM excited state molecular dynamics simulations with explicit surface hopping. Our simulations explain measured quantum yields and excited state lifetimes, and also predict the structures of the hitherto unknown intermediates and of the irreversibly fluorescent state. Further, we find that the proton distribution in the active site of the asFP595 controls the photochemical conversion pathways of the chromophore in the protein matrix. Accordingly, changes in the protonation state of the chromophore and some proximal amino acids lead to different photochemical states, which all turn out to be essential for the photoswitching mechanism. These photochemical states are (i) a neutral chromophore, which can trans-cis photoisomerize, (ii) an anionic chromophore, which rapidly undergoes radiationless decay after excitation, and (iii) a putative fluorescent zwitterionic chromophore. The overall stability of the different protonation states is controlled by the isomeric state of the chromophore. We finally propose that radiation-induced decarboxylation of the glutamic acid Glu215 blocks the proton transfer pathways that enable the deactivation of the zwitterionic chromophore and thus leads to irreversible fluorescence. We have identified the tight coupling of trans-cis isomerization and proton transfers in photoswitchable proteins to be essential for their function and propose a detailed underlying mechanism, which provides a comprehensive picture that explains the available experimental data. The structural similarity between asFP595 and other fluoroproteins of interest for imaging suggests that this coupling is a quite general mechanism for photoswitchable proteins. These insights can guide the rational design and optimization of photoswitchable proteins.     

The protonation state of the catalytic aspartates in plasmepsin II

Assigning the correct protonation state to the catalytic residues is essential for a realistic modelling of an enzyme's active site. Plasmepsins are pharmaceutically relevant aspartic proteases involved in haemoglobin degradation by Plasmodium spp. In aspartic proteases, one of the two catalytic aspartates is protonated, while the other is negatively charged. Here, multiple explicit-water molecular dynamics simulations of plasmepsin II, uncomplexed and with a hydroxypropylamine peptidomimetic inhibitor, indicate that protonation of Asp214 favours a stable active site structure. Moreover, the protonation state of the catalytic aspartate has a strong influence on a linear chain of hydrogen bonds with the adjacent side chains.     

Inhibition of InhA activity,but not KasA activity,induces formation of a KasA-containing complex in mycobacteria

Isoniazid (INH) remains one of the key drugs used to control tuberculosis, with the enoyl-AcpM reductase InhA being the primary target. However, based on the observation that INH-treated Mycobacterium tuberculosis overproduces KasA, an enzyme involved in the biosynthesis of mycolic acids, and induces the formation of a covalent complex consisting of AcpM, KasA, and INH, it has been proposed that KasA represents the primary target of INH. However, the relevance of this complex to INH action remains obscure. This study was aimed at clarifying the role of InhA and KasA in relation to INH activity. By using anti-KasA antibodies we detected the KasA-containing complex in INH-treated Mycobacterium smegmatis. In addition, INH-treated cells also produced constant levels of KasA that were not sequestered in the complex and presumably were sufficient to ensure mycolic acid biosynthesis. Interestingly, a furA-lacking strain induced the complex at lower concentrations of INH compared with the control strain, whereas higher INH concentrations were necessary to induce the complex in a strain that lacks katG, suggesting that INH needs to be activated by KatG to induce the KasA-containing complex. The InhA inhibitors ethionamide and diazaborine also induced the complex; thus, its formation was not specifically relevant to INH action but was because of InhA inhibition. In addition, in vitro assays using purified InhA and KasA demonstrated that KatG-activated INH, triclosan, and diazaborine inhibited InhA but not KasA activity. Moreover, several thermosensitive InhA mutant strains of M. smegmatis constitutively expressed the KasA-containing complex. This study provides the biochemical and genetic evidence. 1) Only inhibition of InhA, but not KasA, induces the KasA-containing complex. 2) INH is not part of the complex. 3) INH does not target KasA, consistent with InhA being the primary target of INH.     

Conformations of higher gangliosides and their binding with cholera toxin - investigation by molecular modeling, molecular mechanics, and molecular dynamics

Molecular mechanics and molecular dynamics studies are performed to investigate the conformational preference of cell surface higher gangliosides (GT1A and GT1B) and their interaction with Cholera Toxin. The water mediated hydrogen bonding network exists between sugar residues in gangliosides. An integrated molecular modeling, molecular mechanics, and molecular dynamics calculation of cholera toxin complexed with GT1A and GT1B reveal that, the active site of cholera toxin can accommodate these higher gangliosides. Direct and water mediated hydrogen bonding interactions stabilize these binding modes and play an essential role in defining the order of specificity for different higher ganglioside towards cholera toxin. This study identifies that the binding site of cholera toxin is shallow and can accommodate a maximum of two NeuNAc residues. The NeuNAc binding site of cholera toxin may be crucial for the design of inhibitors that can prevent the infection of cholera.     

QM/MM investigation of structure and spectroscopic properties of a vanadium-containing peroxidase

We present a comprehensive analysis of the most likely ground state configuration of the resting state of vanadium dependent chloroperoxidase (VCPO) based on quantum mechanics/molecular mechanics (QM/MM) evaluations of ground state properties, UV-vis spectra and NMR chemical shifts. Within the QM/MM framework, density functional theory (DFT) calculations are used to characterize the resting state of VCPO via time-dependent density functional theory (TD-DFT) calculations of electronic excitation energies and NMR chemical shifts. Comparison with available experimental data allows us to determine the most likely protonation state of VCPO, a state which results in a doubly protonated axial oxygen, a site largely stabilized by hydrogen bonds. We found that the bulk of the protein that is beyond the immediate layer surrounding the cofactor, has an important electrostatic effect on the absorption maximum. Through examination of frontier orbitals, we analyze the nature of two bound water molecules and the extent to which relevant residues in the active site influence the spectroscopy calculations.     

The condensing activities of the Mycobacterium tuberculosis type II fatty acid synthase are differentially regulated by phosphorylation

Phosphorylation of proteins by Ser/Thr protein kinases (STPKs) has recently become of major physiological importance because of its possible involvement in virulence of bacterial pathogens. Although Mycobacterium tuberculosis has eleven STPKs, the nature and function of the substrates of these enzymes remain largely unknown. In this work, we have identified for the first time STPK substrates in M. tuberculosis forming part of the type II fatty acid synthase (FAS-II) system involved in mycolic acid biosynthesis: the malonyl-CoA::AcpM transacylase mtFabD, and the beta-ketoacyl AcpM synthases KasA and KasB. All three enzymes were phosphorylated in vitro by different kinases, suggesting a complex network of interactions between STPKs and these substrates. In addition, both KasA and KasB were efficiently phosphorylated in M. bovis BCG each at different sites and could be dephosphorylated by the M. tuberculosis Ser/Thr phosphatase PstP. Enzymatic studies revealed that, whereas phosphorylation decreases the activity of KasA in the elongation process of long chain fatty acids synthesis, this modification enhances that of KasB. Such a differential effect of phosphorylation may represent an unusual mechanism of FAS-II system regulation, allowing pathogenic mycobacteria to produce full-length mycolates, which are required for adaptation and intracellular survival in macrophages.     

Thiolactomycin and related analogues as novel anti-mycobacterial agents targeting KasA and KasB condensing enzymes in Mycobacterium tuberculosis

Prevention efforts and control of tuberculosis are seriously hampered by the appearance of multidrug-resistant strains of Mycobacterium tuberculosis, dictating new approaches to the treatment of the disease. Thiolactomycin (TLM) is a unique thiolactone that has been shown to exhibit anti-mycobacterial activity by specifically inhibiting fatty acid and mycolic acid biosynthesis. In this study, we present evidence that TLM targets two beta-ketoacyl-acyl-carrier protein synthases, KasA and KasB, consistent with the fact that both enzymes belong to the fatty-acid synthase type II system involved in fatty acid and mycolic acid biosynthesis. Overexpression of KasA, KasB, and KasAB in Mycobacterium bovis BCG increased in vivo and in vitro resistance against TLM. In addition, a multidrug-resistant clinical isolate was also found to be highly sensitive to TLM, indicating promise in counteracting multidrug-resistant strains of M. tuberculosis. The design and synthesis of several TLM derivatives have led to compounds more potent both in vitro against fatty acid and mycolic acid biosynthesis and in vivo against M. tuberculosis. Finally, a three-dimensional structural model of KasA has also been generated to improve understanding of the catalytic site of mycobacterial Kas proteins and to provide a more rational approach to the design of new drugs.     

Protein-protein interactions within the Fatty Acid Synthase-II system of Mycobacterium tuberculosis are essential for mycobacterial viability

Despite the existence of efficient chemotherapy, tuberculosis remains a leading cause of mortality worldwide. New drugs are urgently needed to reduce the potential impact of the emergence of multidrug-resistant strains of the causative agent Mycobacterium tuberculosis (Mtb). The front-line antibiotic isoniazid (INH), and several other drugs, target the biosynthesis of mycolic acids and especially the Fatty Acid Synthase-II (FAS-II) elongation system. This biosynthetic pathway is essential and specific for mycobacteria and still represents a valuable system for the search of new anti-tuberculous agents. Several data, in the literature, suggest the existence of protein-protein interactions within the FAS-II system. These interactions themselves might serve as targets for a new generation of drugs directed against Mtb. By using an extensive in vivo yeast two-hybrid approach and in vitro co-immunoprecipitation, we have demonstrated the existence of both homotypic and heterotypic interactions between the known components of FAS-II. The condensing enzymes KasA, KasB and mtFabH interact with each other and with the reductases MabA and InhA. Furthermore, we have designed and constructed point mutations of the FAS-II reductase MabA, able to disrupt its homotypic interactions and perturb the interaction pattern of this protein within FAS-II. Finally, we showed by a transdominant genetic approach that these mutants are dominant negative in both non-pathogenic and pathogenic mycobacteria. These data allowed us to draw a dynamic model of the organization of FAS-II. They also represent an important step towards the design of a new generation of anti-tuberculous agents, as being inhibitors of essential protein-protein interactions.     

Site-specific basicities regulate molecular recognition in receptor binding: in silico docking of thyroid hormones

Interactions between thyroid hormone α and β receptors and the eight protonation microspecies of each of the main thyroid hormones (thyroxine, liothyronine, and reverse liothyronine) were investigated and quantitated by molecular modeling. Flexible docking of the various protonation forms of thyroid hormones and high-affinity thyromimetics to the two thyroid receptors was carried out. In this method the role of the ionization state of each basic site could be studied in the composite process of molecular recognition. Our results quantitate at the molecular level how the ionization state and the charge distribution influence the protein binding. The anionic form of the carboxyl group (i.e., carboxylate site) is essential for protein binding, whereas the protonated form of amino group worsens the binding. The protonation state of the phenolate plays a less important role in the receptor affinity; its protonation, however, alters the electron density and the concomitant stacking propensity of the aromatic rings, resulting in a different binding score. The combined results of docking and microspeciation studies show that microspecies with the highest concentration at the pH of blood are not the strongest binding ones. The calculated binding free energy values can be well interpreted in terms of the interactions between the actual sites of the microspecies and the receptor amino acids. Our docking results were validated and compared with biological data from the literature. Since the thyroid hormone receptors influence several physiologic functions, such as metabolic rate, cholesterol and triglyceride levels, and heart frequency, our binding results provide a molecular basis for drug design and development in related therapeutic indications.     

Theoretical studies on the interaction of biphenyl inhibitors with Mycobacterium tuberculosis protein tyrosine phosphatase MptpB

MptpB is an essential secreted virulence factor for M. tuberculosis. Inhibition of MptpB impairs mycobacterial survival in host macrophages and thus helps reduce tuberculosis infections. However, the binding mode of the biphenyl inhibitors, which are known as some of the most potent MptpB inhibitors, remains unclear. In this study, to understand the interactions between biphenyl inhibitors and MptpB, docking and molecular dynamics simulations were carried out using AutoDock and GROMACS softwares. Calculation results show that all the biphenyl inhibitors can be docked to the binding site of MptpB, with the acid warheads forming a hydrogen bond network at the active site. But the binding modes of other terminals of these inhibitors are different. The cyclohexyl and trifluoromethyl substituents at R1 and R2 sites are necessary for the inhibitors to adopt their double-site binding mechanism. The estimated binding affinities are basically consistent with the experimental results. MD simulations show that these binding complexes display different stability.     

Experimental and 'in silico' analysis of the effect of pH on HIV-1 protease inhibitor affinity: implications for the charge state of the protein ionogenic groups

The pH dependence of the HIV-1 protease inhibitor affinity was studied by determining the interaction kinetics of a series of inhibitors at three pH values by surface plasmon resonance (SPR) biosensor analysis. The results were rationalized by molecular mechanics based protocols that have as a starting point the structures of the HIV-1 protease inhibitor complexes differing in the protonation states as predicted by our calculations. The SPR experiments indicate a variety of binding affinity pH dependencies which are rather well reproduced by our simulations. Moreover, our calculations are able to pinpoint the possible changes in the charged state of the protein binding site and of the inhibitor that underlie the observed effects of the pH on binding affinity. The combination of SPR and molecular mechanics calculations has afforded novel insights into the pH dependence of inhibitor interactions with their target. This work raises the possibility of designing inhibitors with different pH binding affinity profiles to the ones described here.     

A computational study of the protein-ligand interactions in CDK2 inhibitors: using quantum mechanics/molecular mechanics interaction energy as a predictor of the biological activity

We report a combined quantum mechanics/molecular mechanics (QM/MM) study to determine the protein-ligand interaction energy between CDK2 (cyclin-dependent kinase 2) and five inhibitors with the N(2)-substituted 6-cyclohexyl-methoxy-purine scaffold. The computational results in this work show that the QM/MM interaction energy is strongly correlated to the biological activity and can be used as a predictor, at least within a family of substrates. A detailed analysis of the protein-ligand structures obtained from molecular dynamics simulations shows specific interactions within the active site that, in some cases, have not been reported before to our knowledge. The computed interaction energy gauges the strength of protein-ligand interactions. Finally, energy decomposition and multiple regression analyses were performed to check the contribution of the electrostatic and van der Waals energies to the total interaction energy and to show the capabilities of the computational model to identify new potent inhibitors.     

The biosynthesis of mycolic acids in Mycobacterium tuberculosis relies on multiple specialized elongation complexes interconnected by specific protein-protein interactions

Tuberculosis kills about two million people every year and remains one of the leading causes of mortality worldwide. As a result of the increasing antibiotic resistance of Mycobacterium tuberculosis (Mtb) strains, there is an urgent need for new antitubercular drugs. Several efficient antibiotics, including isoniazid, specifically target the fatty acid synthase-II (FAS-II) complex of mycolic acid biosynthesis. We have previously shown that there are protein-protein interactions between the components of FAS-II that are essential for mycobacterial survival. We have now looked at the potential partners of FAS-II, mtFabD, the methyltransferases MmaAs, and Pks13. A combination of yeast two-hybrid and co-immunoprecipitation experiments showed that mtFabD interacts with each beta-ketoacyl-synthase (KasA, KasB and mtFabH) and with the core of FAS-II (InhA and MabA). The methyltransferases have a greater affinity for KasA and KasB than for mtFabH, suggesting that modifications on the meromycolic chains may occur during their elongation. Finally, Pks13, which catalyzes the final Claisen condensation of mycolic acids, interacts specifically with KasB. These data allowed us to determine the architecture of the multiple specialized FAS-II complexes, giving us insights into the organization of the complete mycolic acids biosynthesis. Our studies suggest a new and crucial interaction (KasB-Pks13) as a putative target for peptidomimetic antibiotics.