Functional annotation of putative fadE9 of Mycobacterium tuberculosis as isobutyryl-CoA dehydrogenase involved in valine catabolism

Genome scale metabolic network modelling for metabolic profile predictions

Metabolic profiling (metabolomics) aims at measuring small molecules (metabolites) in complex samples like blood or urine for human health studies. While biomarker-based assessment often relies on a single molecule, metabolic profiling combines several metabolites to create a more complex and more specific fingerprint of the disease. However, in contrast to genomics, there is no unique metabolomics setup able to measure the entire metabolome. This challenge leads to tedious and resource consuming preliminary studies to be able to design the right metabolomics experiment. In that context, computer assisted metabolic profiling can be of strong added value to design metabolomics studies more quickly and efficiently. We propose a constraint-based modelling approach which predicts in silico profiles of metabolites that are more likely to be differentially abundant under a given metabolic perturbation (e.g. due to a genetic disease), using flux simulation. In genome-scale metabolic networks, the fluxes of exchange reactions, also known as the flow of metabolites through their external transport reactions, can be simulated and compared between control and disease conditions in order to calculate changes in metabolite import and export. These import/export flux differences would be expected to induce changes in circulating biofluid levels of those metabolites, which can then be interpreted as potential biomarkers or metabolites of interest. In this study, we present SAMBA (SAMpling Biomarker Analysis), an approach which simulates fluxes in exchange reactions following a metabolic perturbation using random sampling, compares the simulated flux distributions between the baseline and modulated conditions, and ranks predicted differentially exchanged metabolites as potential biomarkers for the perturbation. We show that there is a good fit between simulated metabolic exchange profiles and experimental differential metabolites detected in plasma, such as patient data from the disease database OMIM, and metabolic trait-SNP associations found in mGWAS studies. These biomarker recommendations can provide insight into the underlying mechanism or metabolic pathway perturbation lying behind observed metabolite differential abundances, and suggest new metabolites as potential avenues for further experimental analyses.

In Silico Identification of Novel Derivatives of Rifampicin Targeting Ribonuclease VapC2 of M. tuberculosis H37Rv: Rifampicin Derivatives Target VapC2 of Mtb H37Rv

The emergence of multi-drug-resistant Mycobacterium tuberculosis (Mtb) strains has rendered many of the currently available anti-TB drugs ineffective. Hence, there is a pressing need to discover new potential drug targets/candidates. In this study, attempts have been made to identify novel inhibitors of the ribonuclease VapC2 of Mtb H37Rv using various computational techniques. Ribonuclease VapC2 Mtb H37Rv's protein structure was retrieved from the PDB databank, 22 currently used anti-TB drugs were retrieved from the PubChem database, and protein-ligand interactions were analyzed by docking studies. Out of the 22 drugs, rifampicin (RIF), being a first-line drug, showed the best binding energy (-8.8 Kcal/mol) with Mtb H37Rv VapC2; hence, it was selected as a parent molecule for the design of its derivatives. Based on shape score and radial plot criteria, out of 500 derivatives designed through SPARK (Cresset^®, Royston, UK) program, the 10 best RIF derivatives were selected for further studies. All the selected derivatives followed the ADME criteria concerning drug-likeness. The docking of ribonuclease VapC2 with RIF derivatives revealed the best binding energy of -8.1 Kcal/mol with derivative 1 (i.e., RIF-155841). A quantitative structure-activity relationship study revealed that derivative 1's activity assists in the inhibition of ribonuclease VapC2. The stability of the VapC2-RIF155841 complex was evaluated using molecular dynamics simulations for 50 ns and the complex was found to be stable after 10 nsec. Further, a chemical synthesis scheme was designed for the newly identified RIF derivative (RIF-155841), which verified that its chemical synthesis is possible for future in vitro/in vivo experimental validation. Overall, this study evaluated the potential of the newly designed RIF derivatives with respect to the Mtb VapC2 protein, which is predicted to be involved in some indispensable processes of the related pathogen. Future experimental studies regarding RIF-155841, including the exploration of the remaining RIF derivatives, are warranted to verify our current findings.

PRSS37 deficiency leads to impaired energy metabolism in testis and sperm revealed by DIA-based quantitative proteomic analysis

Our previous studies have reported that a putative trypsin-like serine protease, PRSS37, is exclusively expressed in testicular germ cells during late spermatogenesis and essential for sperm migration from the uterus into the oviduct and sperm-egg recognition via mediating the interaction between PDILT and ADAM3. In the present study, the global proteome profiles of wild-type (wt) and Prss37^-/- mice in testis and sperm were compared employing data independent acquisition (DIA) technology. Overall, 2506 and 459 differentially expressed proteins (DEPs) were identified in Prss37-null testis and sperm, respectively, when compared to control groups. Bioinformatic analyses revealed that most of DEPs were related to energy metabolism. Of note, the DEPs associated with pathways for the catabolism such as glucose via glycolysis, fatty acids via β-oxidation, and amino acids via oxidative deamination were significantly down-regulated. Meanwhile, the DEPs involved in the tricarboxylic acid cycle (TCA cycle) and oxidative phosphorylation (OXPHOS) were remarkably decreased. The DIA data were further confirmed by a markedly reduction of intermediate metabolites (citrate and fumarate) in TCA cycle and terminal metabolite (ATP) in OXPHOS system after disruption of PRSS37. These outcomes not only provide a more comprehensive understanding of the male fertility of energy metabolism modulated by PRSS37 but also furnish a dynamic proteomic resource for further reproductive biology studies.

Multiple acyl-CoA dehydrogenase deficiency kills Mycobacterium tuberculosis in vitro and during infection

The human pathogen Mycobacterium tuberculosis depends on host fatty acids as a carbon source. However, fatty acid β-oxidation is mediated by redundant enzymes, which hampers the development of antitubercular drugs targeting this pathway. Here, we show that rv0338c, which we refer to as etfD, encodes a membrane oxidoreductase essential for β-oxidation in M. tuberculosis. An etfD deletion mutant is incapable of growing on fatty acids or cholesterol, with long-chain fatty acids being bactericidal, and fails to grow and survive in mice. Analysis of the mutant’s metabolome reveals a block in β-oxidation at the step catalyzed by acyl-CoA dehydrogenases (ACADs), which in other organisms are functionally dependent on an electron transfer flavoprotein (ETF) and its cognate oxidoreductase. We use immunoprecipitation to show that M. tuberculosis EtfD interacts with FixA (EtfB), a protein that is homologous to the human ETF subunit β and is encoded in an operon with fixB, encoding a homologue of human ETF subunit α. We thus refer to FixA and FixB as EtfB and EtfA, respectively. Our results indicate that EtfBA and EtfD (which is not homologous to human EtfD) function as the ETF and oxidoreductase for β-oxidation in M. tuberculosis and support this pathway as a potential target for tuberculosis drug development.

The pathogen Mycobacterium tuberculosis depends on host fatty acids and cholesterol as carbon sources. Here, Beites et al. identify a protein complex that is essential for fatty acid and cholesterol utilization and thus for survival of M. tuberculosis during infection, supporting this pathway as a potential target for tuberculosis drug development.

Role of a cysteine residue in substrate entry and catalysis in MtHIBADH: Analysis by chemical modifications and site-directed mutagenesis

Despite sharing conserved substrate-binding residues, members of 3-hydroxyisobutyrate dehydrogenase (HIBADH) superfamily show remarkable differences in substrate preference. Cysteine residues were identified within a radius of 6 Å surrounding both the active site and the substrate entry site of HIBADH enzyme from Mycobacterium tuberculosis (MtHIBADH). Chemical modification with thiol-modifying reagents, pCMB and DTNB, abrogated the dehydrogenase activity of the enzyme. The loss in activity followed pseudo-first-order kinetics as a function of the concentration of pCMB. S-HIBA (substrate) binding provided partial protection, while NAD (cofactor) binding provided ~70% protection from thiol-modifying reagent. Site-directed mutagenesis of cysteine residues present in the MtHIBADH enzyme identified the indispensable role of Cys-210 residue, located at C-terminal domain, for its dehydrogenase activity. Cys-210 mutation to serine reduced the dehydrogenase activity by ~2-fold while mutation to alanine strikingly reduced the activity by ~140-fold. C210A mutation did not perturb the state of oligomerization of the enzyme but perturbed the secondary structure content. Structural analysis revealed the involvement of Cys-210 residue in inter-chain interaction with Gln-178, which acts as hydrogen bond donor and coordinates with Cys-210 and Gly-208 of the adjacent subunit. The data demonstrate a critical role of Cys-210 residue in maintaining the conformation and rigidity of loop composed of substrate-interacting residues involved in the entry of S-HIBA substrate in MtHIBADH.

Structural basis for the broad substrate specificity of two acyl-CoA dehydrogenases FadE5 from mycobacteria

FadE, an acyl-CoA dehydrogenase, introduces unsaturation to carbon chains in lipid metabolism pathways. Here, we report that FadE5 from Mycobacterium tuberculosis (MtbFadE5) and Mycobacterium smegmatis (MsFadE5) play roles in drug resistance and exhibit broad specificity for linear acyl-CoA substrates but have a preference for those with long carbon chains. Here, the structures of MsFadE5 and MtbFadE5, in the presence and absence of substrates, have been determined. These reveal the molecular basis for the broad substrate specificity of these enzymes. FadE5 interacts with the CoA region of the substrate through a large number of hydrogen bonds and an unusual π-π stacking interaction, allowing these enzymes to accept both short- and long-chain substrates. Residues in the substrate binding cavity reorient their side chains to accommodate substrates of various lengths. Longer carbon-chain substrates make more numerous hydrophobic interactions with the enzyme compared with the shorter-chain substrates, resulting in a preference for this type of substrate.