Prediction of high- and low-affinity quinol-analogue-binding sites in the aa3 and bo3 terminal oxidases from Bacillus subtilis and Escherichia coli1

The Dynamics of OXA-23 β-Lactamase from Acinetobacter baumannii

Antibiotic resistance is a pressing topic, which also affects β-lactam antibiotic molecules. Until a few years ago, it was considered no more than an interesting species from an academic point of view, Acinetobacter baumanii is today one of the most serious threats to public health, so much so that it has been declared one of the species for which the search for new antibiotics, or new ways to avoid its resistance, is an absolute priority according to WHO. Although there are several molecular mechanisms that are responsible for the extreme resistance of A. baumanii to antibiotics, a class D β-lactamase is the main cause for the clinical concern of this bacterial species. In this work, we analyzed the A. baumanii OXA-23 protein via molecular dynamics. The results obtained show that this protein is able to assume different conformations, especially in some regions around the active site. Part of the OXA-23 protein has considerable conformational motility, while the rest is less mobile. The importance of these observations for understanding the functioning mechanism of the enzyme as well as for designing new effective molecules for the treatment of A. baumanii is discussed.

Pyruvate dehydrogenase operates as an intramolecular nitroxyl generator during macrophage metabolic reprogramming

M1 macrophages enter a glycolytic state when endogenous nitric oxide (NO) reprograms mitochondrial metabolism by limiting aconitase 2 and pyruvate dehydrogenase (PDH) activity. Here, we provide evidence that NO targets the PDH complex by using lipoate to generate nitroxyl (HNO). PDH E2-associated lipoate is modified in NO-rich macrophages while the PDH E3 enzyme, also known as dihydrolipoamide dehydrogenase (DLD), is irreversibly inhibited. Mechanistically, we show that lipoate facilitates NO-mediated production of HNO, which interacts with thiols forming irreversible modifications including sulfinamide. In addition, we reveal a macrophage signature of proteins with reduction-resistant modifications, including in DLD, and identify potential HNO targets. Consistently, DLD enzyme is modified in an HNO-dependent manner at Cys477 and Cys484, and molecular modeling and mutagenesis show these modifications impair the formation of DLD homodimers. In conclusion, our work demonstrates that HNO is produced physiologically. Moreover, the production of HNO is dependent on the lipoate-rich PDH complex facilitating irreversible modifications that are critical to NO-dependent metabolic rewiring.

AI-Aided Search for New HIV-1 Protease Ligands

The availability of drugs capable of blocking the replication of microorganisms has been one of the greatest triumphs in the history of medicine, but the emergence of an ever-increasing number of resistant strains poses a serious problem for the treatment of infectious diseases. The search for new potential ligands for proteins involved in the life cycle of pathogens is, therefore, an extremely important research field today. In this work, we have considered the HIV-1 protease, one of the main targets for AIDS therapy. Several drugs are used today in clinical practice whose mechanism of action is based on the inhibition of this enzyme, but after years of use, even these molecules are beginning to be interested by resistance phenomena. We used a simple artificial intelligence system for the initial screening of a data set of potential ligands. These results were validated by docking and molecular dynamics, leading to the identification of a potential new ligand of the enzyme which does not belong to any known class of HIV-1 protease inhibitors. The computational protocol used in this work is simple and does not require large computational power. Furthermore, the availability of a large number of structural information on viral proteins and the presence of numerous experimental data on their ligands, with which it is possible to compare the results obtained with computational methods, make this research field the ideal terrain for the application of these new computational techniques.

Targeting mitochondrial impairment for the treatment of cardiovascular diseases: From hypertension to ischemia-reperfusion injury, searching for new pharmacological targets

Mitochondria and mitochondrial proteins represent a group of promising pharmacological target candidates in the search of new molecular targets and drugs to counteract the onset of hypertension and more in general cardiovascular diseases (CVDs). Indeed, several mitochondrial pathways result impaired in CVDs, showing ATP depletion and ROS production as common traits of cardiac tissue degeneration. Thus, targeting mitochondrial dysfunction in cardiomyocytes can represent a successful strategy to prevent heart failure. In this context, the identification of new pharmacological targets among mitochondrial proteins paves the way for the design of new selective drugs. Thanks to the advances in omics approaches, to a greater availability of mitochondrial crystallized protein structures and to the development of new computational approaches for protein 3D-modelling and drug design, it is now possible to investigate in detail impaired mitochondrial pathways in CVDs. Furthermore, it is possible to design new powerful drugs able to hit the selected pharmacological targets in a highly selective way to rescue mitochondrial dysfunction and prevent cardiac tissue degeneration. The role of mitochondrial dysfunction in the onset of CVDs appears increasingly evident, as reflected by the impairment of proteins involved in lipid peroxidation, mitochondrial dynamics, respiratory chain complexes, and membrane polarization maintenance in CVD patients. Conversely, little is known about proteins responsible for the cross-talk between mitochondria and cytoplasm in cardiomyocytes. Mitochondrial transporters of the SLC25A family, in particular, are responsible for the translocation of nucleotides (e.g., ATP), amino acids (e.g., aspartate, glutamate, ornithine), organic acids (e.g. malate and 2-oxoglutarate), and other cofactors (e.g., inorganic phosphate, NAD⁺, FAD, carnitine, CoA derivatives) between the mitochondrial and cytosolic compartments. Thus, mitochondrial transporters play a key role in the mitochondria-cytosol cross-talk by leading metabolic pathways such as the malate/aspartate shuttle, the carnitine shuttle, the ATP export from mitochondria, and the regulation of permeability transition pore opening. Since all these pathways are crucial for maintaining healthy cardiomyocytes, mitochondrial carriers emerge as an interesting class of new possible pharmacological targets for CVD treatments.

Structural evolution of Delta lineage of SARS-CoV-2

One of the main obstacles in prevention and treatment of COVID-19 is the rapid evolution of the SARS-CoV-2 Spike protein. Given that Spike is the main target of common treatments of COVID-19, mutations occurring at this virulent factor can affect the effectiveness of treatments. The B.1.617.2 lineage of SARS-CoV-2, being characterized by many Spike mutations inside and outside of its receptor-binding domain (RBD), shows high infectivity and relative resistance to existing cures. Here, utilizing a wide range of computational biology approaches, such as immunoinformatics, molecular dynamics (MD), analysis of intrinsically disordered regions (IDRs), protein-protein interaction analyses, residue scanning, and free energy calculations, we examine the structural and biological attributes of the B.1.617.2 Spike protein. Furthermore, the antibody design protocol of Rosetta was implemented for evaluation the stability and affinity improvement of the Bamlanivimab (LY-CoV55) antibody, which is not capable of interactions with the B.1.617.2 Spike. We observed that the detected mutations in the Spike of the B1.617.2 variant of concern can cause extensive structural changes compatible with the described variation in immunogenicity, secondary and tertiary structure, oligomerization potency, Furin cleavability, and drug targetability. Compared to the Spike of Wuhan lineage, the B.1.617.2 Spike is more stable and binds to the Angiotensin-converting enzyme 2 (ACE2) with higher affinity.

New Insights Regarding Hemin Inhibition of the Purified Rat Brain 2-Oxoglutarate Carrier and Relationships with Mitochondrial Dysfunction

A kinetic analysis of the transport assays on the purified rat brain 2-oxoglutarate/malate carrier (OGC) was performed starting from our recent results reporting about a competitive inhibitory behavior of hemin, a physiological porphyrin derivative, on the OGC reconstituted in an active form into proteoliposomes. The newly provided transport data and the elaboration of the kinetic equations show evidence that hemin exerts a mechanism of partially competitive inhibition, coupled with the formation of a ternary complex hemin-carrier substrate, when hemin targets the OGC from the matrix face. A possible interpretation of the provided kinetic analysis, which is supported by computational studies, could indicate the existence of a binding region responsible for the inhibition of the OGC and supposedly involved in the regulation of OGC activity. The proposed regulatory binding site is located on OGC mitochondrial matrix loops, where hemin could establish specific interactions with residues involved in the substrate recognition and/or conformational changes responsible for the translocation of mitochondrial carrier substrates. The regulatory binding site would be placed about 6 Å below the substrate binding site of the OGC, facing the mitochondrial matrix, and would allow the simultaneous binding of hemin and 2-oxoglutarate or malate to different regions of the carrier. Overall, the presented experimental and computational analyses help to shed light on the possible existence of the hemin-carrier substrate ternary complex, confirming the ability of the OGC to bind porphyrin derivatives, and in particular hemin, with possible consequences for the mitochondrial redox state mediated by the malate/aspartate shuttle led by the mitochondrial carriers OGC and AGC.

Proline/Glycine residues of the PG-levels guide conformational changes along the transport cycle in the mitochondrial carnitine/acylcarnitine carrier (SLC25A20)

Mitochondrial carnitine/acylcarnitine carrier (CAC) is a member of the mitochondrial carrier (MC) family and imports acylcarnitine into the mitochondrial matrix in exchange for carnitine, playing a pivotal role in carnitine shuttle, crucial for fatty acid oxidation. The crystallized structure of CAC has not been solved yet, however, the availability of several in vitro/in silico studies, also based on the crystallized structures of the ADP/ATP carrier in the cytosolic-conformation and in the matrix-conformation, has made possible to confirm the hypothesis of the single-binding centered-gated pore mechanism for all the members of the MC family. In addition, our recent bioinformatics analyses allowed quantifying in silico the importance of protein residues of MC substrate binding region, of those involved in the formation of the matrix and cytosolic gates, and of those belonging to the Pro/Gly (PG) levels, proposed to be crucial for the tilting/kinking/bending of the six MC transmembrane helices, funneling the substrate translocation pathway. Here we present a combined in silico/in vitro analysis employed for investigating the role played by a group of 6 proline residues and 6 glycine residues, highly conserved in CAC, belonging to MC PG-levels. Residues of the PG-levels surround the similarly located MC common substrate binding region, and were proposed to lead conformational changes and substrate translocation, following substrate binding. For our analysis, we employed 3D molecular modeling approaches, alanine scanning site-directed mutagenesis and in vitro transport assays. Our analysis reveals that P130 (H3), G268 (H6) and G220 (H5), mutated in alanine, affect severely CAC transport activity (mutant catalytic efficiency lower than 5 % compared to the wild type CAC), most likely due to their major role in triggering CAC conformational changes, following carnitine binding. Notably, P30A (H1) and G121A (H3) CAC mutants, increase the carnitine uptake up to 217 % and 112 %, respectively, compared to the wild type CAC.

BRAFV600E;K601Q metastatic melanoma patient-derived organoids and docking analysis to predict the response to targeted therapy

The V600E mutation in BRAF is associated with increased phosphorylation of Erk1/2 and high sensitivity to BRAFi/MEKi combination in metastatic melanoma. In very few patients, a tandem mutation in BRAF, V600 and K601, causes a different response to BRAFi/MEKi combination. BRAF^V600E;K601Q patient-derived organoids (PDOs) were generated to investigate targeted therapy efficacy and docking analysis was used to assess BRAF^V600E;K601Q interactions with Vemurafenib. PDOs were not sensitive to Vemurafenib and Cobimetinib given alone and sensitive to their combination, although not as responsive as BRAF^V600E PDOs. The docking analysis justified such a result showing that the tandem mutation in BRAF reduced the affinity for Vemurafenib. Tumor analysis showed that BRAF^V600E;K601Q displayed both increased phosphorylation of Erk1/2 at cytoplasmic level and activation of Notch resistance signaling. This prompted us to inhibit Notch signaling with Nirogacestat, achieving a greater antitumor response and providing PDOs-based evaluation of treatment efficacy in such rare metastatic melanoma.

Personalized Medicine in Mitochondrial Health and Disease: Molecular Basis of Therapeutic Approaches Based on Nutritional Supplements and Their Analogs

Mitochondrial diseases (MDs) may result from mutations affecting nuclear or mitochondrial genes, encoding mitochondrial proteins, or non-protein-coding mitochondrial RNA. Despite the great variability of affected genes, in the most severe cases, a neuromuscular and neurodegenerative phenotype is observed, and no specific therapy exists for a complete recovery from the disease. The most used treatments are symptomatic and based on the administration of antioxidant cocktails combined with antiepileptic/antipsychotic drugs and supportive therapy for multiorgan involvement. Nevertheless, the real utility of antioxidant cocktail treatments for patients affected by MDs still needs to be scientifically demonstrated. Unfortunately, clinical trials for antioxidant therapies using α-tocopherol, ascorbate, glutathione, riboflavin, niacin, acetyl-carnitine and coenzyme Q have met a limited success. Indeed, it would be expected that the employed antioxidants can only be effective if they are able to target the specific mechanism, i.e., involving the central and peripheral nervous system, responsible for the clinical manifestations of the disease. Noteworthily, very often the phenotypes characterizing MD patients are associated with mutations in proteins whose function does not depend on specific cofactors. Conversely, the administration of the antioxidant cocktails might determine the suppression of endogenous oxidants resulting in deleterious effects on cell viability and/or toxicity for patients. In order to avoid toxicity effects and before administering the antioxidant therapy, it might be useful to ascertain the blood serum levels of antioxidants and cofactors to be administered in MD patients. It would be also worthwhile to check the localization of mutations affecting proteins whose function should depend (less or more directly) on the cofactors to be administered, for estimating the real need and predicting the success of the proposed cofactor/antioxidant-based therapy.

The Non-Gastric H+/K+ ATPase (ATP12A) Is Expressed in Mammalian Spermatozoa

H⁺/K⁺ ATPase Type 2 is an heteromeric membrane protein involved in cation transmembrane transport and consists of two subunits: a specific α subunit (ATP12A) and a non-specific β subunit. The aim of this study was to demonstrate the presence and establish the localization of ATP12A in spermatozoa from Bubalus bubalis, Bos taurus and Ovis aries. Immunoblotting revealed, in all three species, a major band (100 kDa) corresponding to the expected molecular mass. The ATP12A immunolocalization pattern showed, consistently in the three species, a strong signal at the acrosome. These results, described here for the first time in spermatozoa, are consistent with those observed for the β₁ subunit of Na⁺/K⁺ ATPase, suggesting that the latter may assemble with the α subunit to produce a functional ATP12A dimer in sperm cells. The above scenario appeared to be nicely supported by 3D comparative modeling and interaction energy calculations. The expression of ATP12A during different stages of bovine sperm maturation progressively increased, moving from epididymis to deferent ducts. Based on overall results, we hypothesize that ATP12A may play a role in acrosome reactions. Further studies will be required in order to address the functional role of this target protein in sperm physiology.

Modeling SARS-CoV-2 spike/ACE2 protein-protein interactions for predicting the binding affinity of new spike variants for ACE2, and novel ACE2 structurally related human protein targets, for COVID-19 handling in the 3PM context

Aims

The rapid spread of new SARS-CoV-2 variants has highlighted the crucial role played in the infection by mutations occurring at the SARS-CoV-2 spike receptor binding domain (RBD) in the interactions with the human ACE2 receptor. In this context, it urgently needs to develop new rapid tools for quickly predicting the affinity of ACE2 for the SARS-CoV-2 spike RBD protein variants to be used with the ongoing SARS-CoV-2 genomic sequencing activities in the clinics, aiming to gain clues about the transmissibility and virulence of new variants, to prevent new outbreaks and to quickly estimate the severity of the disease in the context of the 3PM.

Methods

In our study, we used a computational pipeline for calculating the interaction energies at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface for a selected group of characterized infectious variants of concern/interest (VoC/VoI). By using our pipeline, we built 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for the VoC B.1.1.7-United Kingdom (carrying the mutations of concern/interest N501Y, S494P, E484K at the RBD), P.1-Japan/Brazil (RBD mutations: K417T, E484K, N501Y), B.1.351-South Africa (RBD mutations: K417N, E484K, N501Y), B.1.427/B.1.429-California (RBD mutations: L452R), the B.1.141 (RBD mutations: N439K), and the recent B.1.617.1-India (RBD mutations: L452R; E484Q) and the B.1.620 (RBD mutations: S477N; E484K). Then, we used the obtained 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for predicting the interaction energies at the protein–protein interface.

Results

Along SARS-CoV-2 mutation database screening and mutation localization analysis, it was ascertained that the most dangerous mutations at VoC/VoI spike proteins are located mainly at three regions of the SARS-CoV-2 spike “boat-shaped” receptor binding motif, on the RBD domain. Notably, the P.1 Japan/Brazil variant present three mutations, K417T, E484K, N501Y, located along the entire receptor binding motif, which apparently determines the highest interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface, among those calculated. Conversely, it was also observed that the replacement of a single acidic/hydrophilic residue with a basic residue (E484K or N439K) at the “stern” or “bow” regions, of the boat-shaped receptor binding motif on the RBD, appears to determine an interaction energy with ACE2 receptor higher than that observed with single mutations occurring at the “hull” region or with other multiple mutants. In addition, our pipeline allowed searching for ACE2 structurally related proteins, i.e., THOP1 and NLN, which deserve to be investigated for their possible involvement in interactions with the SARS-CoV-2 spike protein, in those tissues showing a low expression of ACE2, or as a novel receptor for future spike variants. A freely available web-tool for the in silico calculation of the interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface, starting from the sequences of the investigated spike and/or ACE2 variants, was made available for the scientific community at: https://www.mitoairm.it/covid19affinities.

Conclusion

In the context of the PPPM/3PM, the employment of the described pipeline through the provided webservice, together with the ongoing SARS-CoV-2 genomic sequencing, would help to predict the transmissibility of new variants sequenced from future patients, depending on SARS-CoV-2 genomic sequencing activities and on the specific amino acid replacement and/or on its location on the SARS-CoV-2 spike RBD, to put in play all the possible counteractions for preventing the most deleterious scenarios of new outbreaks, taking into consideration that a greater transmissibility has not to be necessarily related to a more severe manifestation of the disease.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-021-00267-w.

N-acetylaspartate release by glutaminolytic ovarian cancer cells sustains protumoral macrophages

Glutaminolysis is known to correlate with ovarian cancer aggressiveness and invasion. However, how this affects the tumor microenvironment is elusive. Here, we show that ovarian cancer cells become addicted to extracellular glutamine when silenced for glutamine synthetase (GS), similar to naturally occurring GS-low, glutaminolysis-high ovarian cancer cells. Glutamine addiction elicits a crosstalk mechanism whereby cancer cells release N-acetylaspartate (NAA) which, through the inhibition of the NMDA receptor, and synergistically with IL-10, enforces GS expression in macrophages. In turn, GS-high macrophages acquire M2-like, tumorigenic features. Supporting this in␣vitro model, in silico data and the analysis of ascitic fluid isolated from ovarian cancer patients prove that an M2-like macrophage phenotype, IL-10 release, and NAA levels positively correlate with disease stage. Our study uncovers the unprecedented role of glutamine metabolism in modulating macrophage polarization in highly invasive ovarian cancer and highlights the anti-inflammatory, protumoral function of NAA.

Zoledronic Acid as a Novel Dual Blocker of KIR6.1/2-SUR2 Subunits of ATP-Sensitive K+ Channels: Role in the Adverse Drug Reactions

Zoledronic acid (ZOL) is used as a bone-specific antiresorptive drug with antimyeloma effects. Adverse drug reactions (A.D.R.) are associated with ZOL-therapy, whose mechanics are unknown. ZOL is a nitrogen-containing molecule whose structure shows similarities with nucleotides, ligands of ATP-sensitive K⁺ (KATP) channels. We investigated the action of ZOL by performing in vitro patch-clamp experiments on native KATP channels in murine skeletal muscle fibers, bone cells, and recombinant subunits in cell lines, and by in silico docking the nucleotide site on KIR and SUR, as well as the glibenclamide site. ZOL fully inhibited the KATP currents recorded in excised macro-patches from Extensor digitorum longus (EDL) and Soleus (SOL) muscle fibers with an IC₅₀ of 1.2 ± 1.4 × 10⁻⁶ and 2.1 ± 3.7 × 10⁻¹⁰ M, respectively, and the KATP currents recorded in cell-attached patches from primary long bone cells with an IC₅₀ of 1.6 ± 2.8 × 10⁻¹⁰ M. ZOL fully inhibited a whole-cell KATP channel current of recombinant KIR6.1-SUR2B and KIR6.2-SUR2A subunits expressed in HEK293 cells with an IC₅₀ of 3.9 ± 2.7 × 10⁻¹⁰ M and 7.1 ± 3.1 × 10⁻⁶ M, respectively. The rank order of potency in inhibiting the KATP currents was: KIR6.1-SUR2B/SOL-KATP/osteoblast-KATP > KIR6.2-SUR2A/EDL-KATP >>> KIR6.2-SUR1 and KIR6.1-SUR1. Docking investigation revealed that the drug binds to the ADP/ATP sites on KIR6.1/2 and SUR2A/B and on the sulfonylureas site showing low binding energy <6 Kcal/mol for the KIR6.1/2-SUR2 subunits vs. the <4 Kcal/mol for the KIR6.2-SUR1. The IC₅₀ of ZOL to inhibit the KIR6.1/2-SUR2A/B channels were correlated with its musculoskeletal and cardiovascular risks. We first showed that ZOL blocks at subnanomolar concentration musculoskeletal KATP channels and cardiac and vascular KIR6.2/1-SUR2 channels.

Manganese impairs the QoxABCD terminal oxidase leading to respiration-associated toxicity

Cell physiology relies on metalloenzymes and can be easily disrupted by imbalances in metal ion pools. Bacillus subtilis requires manganese for growth and has highly regulated mechanisms for import and efflux that help maintain homeostasis. Cells defective for manganese (Mn) efflux are highly sensitive to intoxication, but the processes impaired by Mn excess are often unknown. Here, we employed a forward genetics approach to identify pathways affected by manganese intoxication. Our results highlight a central role for the membrane-localized electron transport chain in metal intoxication during aerobic growth. In the presence of elevated manganese, there is an increased generation of reactive radical species associated with dysfunction of the major terminal oxidase, the cytochrome aa₃ heme-copper menaquinol oxidase (QoxABCD). Intoxication is suppressed by diversion of menaquinol to alternative oxidases or by a mutation affecting heme A synthesis that is known to convert QoxABCD from an aa₃ to a bo₃ -type oxidase. Manganese sensitivity is also reduced by derepression of the MhqR regulon, which protects cells against reactive quinones. These results suggest that dysfunction of the cytochrome aa₃ -type quinol oxidase contributes to metal-induced intoxication.

SARS-CoV-2 Main Protease Active Site Ligands in the Human Metabolome

In late 2019, a global pandemic occurred. The causative agent was identified as a member of the Coronaviridae family, called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In this study, we present an analysis on the substances identified in the human metabolome capable of binding the active site of the SARS-CoV-2 main protease (Mpro). The substances present in the human metabolome have both endogenous and exogenous origins. The aim of this research was to find molecules whose biochemical and toxicological profile was known that could be the starting point for the development of antiviral therapies. Our analysis revealed numerous metabolites-including xenobiotics-that bind this protease, which are essential to the lifecycle of the virus. Among these substances, silybin, a flavolignan compound and the main active component of silymarin, is particularly noteworthy. Silymarin is a standardized extract of milk thistle, Silybum marianum, and has been shown to exhibit antioxidant, hepatoprotective, antineoplastic, and antiviral activities. Our results-obtained in silico and in vitro-prove that silybin and silymarin, respectively, are able to inhibit Mpro, representing a possible food-derived natural compound that is useful as a therapeutic strategy against COVID-19.

Targeting mitochondrial metabolite transporters in Penicillium expansum for reducing patulin production

There is an increasing need of alternative treatments to control fungal infection and consequent mycotoxin accumulation in harvested fruits and vegetables. Indeed, only few biological targets of antifungal agents have been characterized and can be used for limiting fungal spread from decayed fruits/vegetables to surrounding healthy ones during storage. On this concern, a promising target of new antifungal treatments may be represented by mitochondrial proteins due to some species-specific functions played by mitochondria in fungal morphogenesis, drug resistance and virulence. One of the most studied mycotoxins is patulin produced by several species of Penicillium and Aspergillus genera. Patulin is toxic to many biological systems including bacteria, higher plants and animalia. Although precise biochemical mechanisms of patulin toxicity in humans are not completely clarified, its high presence in fresh and processed apple fruits and other apple-based products makes necessary developing a strategy for limiting its presence/accumulation. Patulin biosynthetic pathway consists of an enzymatic cascade, whose first step is represented by the synthesis of 6-methylsalicylic acid, obtained from the condensation of one acetyl-CoA molecule with three malonyl-CoA molecules. The most abundant acetyl-CoA precursor is represented by citrate produced by mitochondria. In the present investigation we report about the possibility to control patulin production through the inhibition of mitochondrial/peroxisome transporters involved in the export of acetyl-CoA precursors from mitochondria and/or peroxisomes, with specific reference to the predicted P. expansum mitochondrial Ctp1p, DTC, Sfc1p, Oac1p and peroxisomal PXN carriers.

Targeting Penicillium expansum GMC Oxidoreductase with High Affinity Small Molecules for Reducing Patulin Production

Simple Summary

With the urgent necessity of potential treatments for limiting mycotoxin production and postharvest fungal rots, we propose a combined in silico/in vitro/in vivo strategy for the rapid and effective identification of bioactive small molecules, chosen among a chemical library hosting approved drugs and phytochemicals, to be used after harvest. The molecular target of our analysis was the GMC oxidoreductase from Penicillium expansum involved in the biosynthesis of patulin, a mycotoxin that can contaminate many foods, especially fruits and fruit-based products. The employed in silico/in vitro/in vivo assays described in our study proved the effectiveness of our strategy and in particular of two small molecules, 6-hydroxycoumarin (structurally related to umbelliferon, an already characterized patulin synthase inhibitor) and meticrane (an already approved drug) in reducing patulin accumulation. Our findings highly recommend the mentioned ligands to be subjected to further analysis for being used in the next future in place of other more toxic compounds, in postharvest treatments based on dipping or drenching methods.

Abstract

Flavine adenine dinucleotide (FAD) dependent glucose methanol choline oxidoreductase (GMC oxidoreductase) is the terminal key enzyme of the patulin biosynthetic pathway. GMC oxidoreductase catalyzes the oxidative ring closure of (E)-ascladiol to patulin. Currently, no protein involved in the patulin biosynthesis in Penicillium expansum has been experimentally characterized or solved by X-ray diffraction. Consequently, nothing is known about P. expansum GMC oxidoreductase substrate-binding site and mode of action. In the present investigation, a 3D comparative model for P. expansum GMC oxidoreductase has been described. Furthermore, a multistep computational approach was used to identify P. expansum GMC oxidoreductase residues involved in the FAD binding and in substrate recognition. Notably, the obtained 3D comparative model of P. expansum GMC oxidoreductase was used for performing a virtual screening of a chemical/drug library, which allowed to predict new GMC oxidoreductase high affinity ligands to be tested in in vitro/in vivo assays. In vitro assays performed in presence of 6-hydroxycoumarin and meticrane, among the highly affinity predicted binders, confirmed a dose-dependent inhibition (17–81%) of patulin production by 6-hydroxycoumarin (10 µM–1 mM concentration range), whereas the approved drug meticrane inhibited patulin production by 43% already at 10 µM. Furthermore, 6-hydroxycoumarin and meticrane caused a 60 and 41% reduction of patulin production, respectively, in vivo on apples at 100 µg/wound.

Glufosinate constrains synchronous and metachronous metastasis by promoting anti-tumor macrophages

Glutamine synthetase (GS) generates glutamine from glutamate and controls the release of inflammatory mediators. In macrophages, GS activity, driven by IL10, associates to the acquisition of M2-like functions. Conditional deletion of GS in macrophages inhibits metastasis by boosting the formation of anti-tumor, M1-like, tumor-associated macrophages (TAMs). From this basis, we evaluated the pharmacological potential of GS inhibitors in targeting metastasis, identifying glufosinate as a specific human GS inhibitor. Glufosinate was tested in both cultured macrophages and on mice bearing metastatic lung, skin and breast cancer. We found that glufosinate rewires macrophages toward an M1-like phenotype both at the primary tumor and metastatic site, countering immunosuppression and promoting vessel sprouting. This was also accompanied to a reduction in cancer cell intravasation and extravasation, leading to synchronous and metachronous metastasis growth inhibition, but no effects on primary tumor growth. Glufosinate treatment was well-tolerated, without liver and brain toxicity, nor hematopoietic defects. These results identify GS as a druggable enzyme to rewire macrophage functions and highlight the potential of targeting metabolic checkpoints in macrophages to treat cancer metastasis.

The oxygen-oxygen distance of water in crystallographic data sets

Water is a key component of cellular biochemistry and numerous water molecules are visible in crystallographic structures. Here we report a series of data sets of crystallographic water: a high resolution data set, a cytochrome c oxidase (subunit I) data set and a carbonic anhydrase data set. These data support the evidence that short distance water molecule pairs are present both at the surface and inside the cavities of proteins. These data are related to article entitled “Oxygen-oxygen distances in protein-bound crystallographic water suggest the presence of protonated clusters” (Palese, 2020) [1].

FAD/NADH Dependent Oxidoreductases: From Different Amino Acid Sequences to Similar Protein Shapes for Playing an Ancient Function

Flavoprotein oxidoreductases are members of a large protein family of specialized dehydrogenases, which include type II NADH dehydrogenase, pyridine nucleotide-disulphide oxidoreductases, ferredoxin-NAD+ reductases, NADH oxidases, and NADH peroxidases, playing a crucial role in the metabolism of several prokaryotes and eukaryotes. Although several studies have been performed on single members or protein subgroups of flavoprotein oxidoreductases, a comprehensive analysis on structure-function relationships among the different members and subgroups of this great dehydrogenase family is still missing. Here, we present a structural comparative analysis showing that the investigated flavoprotein oxidoreductases have a highly similar overall structure, although the investigated dehydrogenases are quite different in functional annotations and global amino acid composition. The different functional annotation is ascribed to their participation in species-specific metabolic pathways based on the same biochemical reaction, i.e., the oxidation of specific cofactors, like NADH and FADH₂. Notably, the performed comparative analysis sheds light on conserved sequence features that reflect very similar oxidation mechanisms, conserved among flavoprotein oxidoreductases belonging to phylogenetically distant species, as the bacterial type II NADH dehydrogenases and the mammalian apoptosis-inducing factor protein, until now retained as unique protein entities in Bacteria/Fungi or Animals, respectively. Furthermore, the presented computational analyses will allow consideration of FAD/NADH oxidoreductases as a possible target of new small molecules to be used as modulators of mitochondrial respiration for patients affected by rare diseases or cancer showing mitochondrial dysfunction, or antibiotics for treating bacterial/fungal/protista infections.

Explaining leak states in the proton pump of heme-copper oxidases observed in single-molecule experiments

Heme-copper oxidases couple the exergonic oxygen reduction with the endergonic proton translocation. Redox-linked structural changes have been localized in deeply buried regions of the protein, near the low-potential heme. How these movements can modulate distant gating events along the intramolecular proton path, where the entry (exit) of pumped proton occurs, is a major concern for the proton pump models. Generally, these models associate, more or less directly, all translocation events with redox transitions. Although they can account for many phenomenological aspects of the pump, evidences from single-molecules experiments about leak states of the pump represent a formidable challenge. Disconnecting the redox-linked pK_a shifts of the proton loading site from the external barriers, we obtain a simple stochastic mechanism which behaves similarly to the real enzyme, able to reverse the flow of the proton transfer.

The Ubiquinol Binding Site of Cytochrome bo3 from Escherichia coli Accommodates Menaquinone and Stabilizes a Functional Menasemiquinone

Cytochrome bo₃, one of three terminal oxygen reductases in the aerobic respiratory chain of Escherichia coli, has been well characterized as a ubiquinol oxidase. The ability of cytochrome bo₃ to catalyze the two-electron oxidation of ubiquinol-8 requires the enzyme to stabilize the one-electron oxidized ubisemiquinone species that is a transient intermediate in the reaction. Cytochrome bo₃ has been shown recently to also utilize demethylmenaquinol-8 as a substrate that, along with menaquinol-8, replaces ubiquinol-8 when E. coli is grown under microaerobic or anaerobic conditions. In this work, we show that its steady-state turnover with 2,3-dimethyl-1,4-naphthoquinol, a water-soluble menaquinol analogue, is just as efficient as with ubiquinol-1. Using pulsed electron paramagnetic resonance spectroscopy, we demonstrate that the same residues in cytochrome bo₃ that stabilize the semiquinone state of ubiquinone also stabilize the semiquinone state of menaquinone, with the hydrogen bond strengths and the distribution of unpaired spin density accommodated for the different substrate. Catalytic function with menaquinol is more tolerant of mutations at the active site than with ubiquinol. A mutation of one of the stabilizing residues (R71H in subunit I) that eliminates the ubiquinol oxidase activity of cytochrome bo₃ does not abolish activity with soluble menaquinol analogues.

Cytochrome c oxidase structures suggest a four-state stochastic pump mechanism

A simple stochastic model for a cytochrome c oxidase proton pump.

Conformations of the HIV-1 protease: A crystal structure data set analysis

The HIV protease is an important drug target for HIV/AIDS therapy, and its structure and function have been extensively investigated. This enzyme performs an essential role in viral maturation by processing specific cleavage sites in the Gag and Gag-Pol precursor polyproteins so as to release their mature forms. This 99 amino acid aspartic protease works as a homodimer, with the active site localized in a central cavity capped by two flexible flap regions. The dimer presents closed or open conformations, which are involved in the substrate binding and release. Here the results of the analysis of a HIV-1 protease data set containing 552 dimer structures are reported. Different dimensionality reduction methods have been used in order to get information from this multidimensional database. Most of the structures in the data set belong to two conformational clusters. An interesting observation that comes from the analysis of these data is that some protease sequences are localized preferentially in specific areas of the conformational landscape of this protein.

Location of the Substrate Binding Site of the Cytochrome bo3 Ubiquinol Oxidase from Escherichia coli

Cytochrome bo₃ is a respiratory proton-pumping oxygen reductase that is a member of the heme-copper superfamily that utilizes ubiquinol-8 (Q₈H₂) as a substrate. The current consensus model has Q₈H₂ oxidized at a low affinity site (Q_L), passing electrons to a tightly bound quinone cofactor at a high affinity site (Q_H site) that stabilizes the one-electron reduced ubisemiquinone, facilitating the transfer of electrons to the redox active metal centers where O₂ is reduced to water. The current work shows that the Q₈ bound to the Q_H site is more dynamic than previously thought. In addition, mutations of residues at the Q_H site that do not abolish activity have been re-examined and shown to have properties expected of mutations at the substrate binding site (Q_L): an increase in the K_M of the substrate ubiquinol-1 (up to 4-fold) and an increase in the apparent K_i of the inhibitor HQNO (up to 8-fold). The data suggest that there is only one binding site for ubiquinol in cyt bo₃ and that site corresponds to the Q_H site.

Searching for the low affinity ubiquinone binding site in cytochrome bo3 from Escherichia coli

The cytochrome bo₃ ubiquinol oxidase is one of three respiratory oxygen reductases in the aerobic respiratory chain of Escherichia coli. The generally accepted model of catalysis assumes that cyt bo₃ contains two distinct ubiquinol binding sites: (i) a low affinity (Q_L) site which is the traditional substrate binding site; and (ii) a high affinity (Q_H) site where a "permanently" bound quinone acts as a cofactor, taking two electrons from the substrate quinol and passing them one-by-one to the heme b component of the enzyme which, in turn, transfers them to the heme o₃/Cu_B active site. Whereas the residues at the Q_H site are well defined, the location of the Q_L site remains unknown. The published X-ray structure does not contain quinone, and substantial amounts of the protein are missing as well. A recent bioinformatics study by Bossis et al. [Biochem J. (2014) 461, 305-314] identified a sequence motif G¹⁶³EFX₃GWX₂Y¹⁷³ as the likely Q_L site in the family of related quinol oxidases. In the current work, this was tested by site-directed mutagenesis. The results show that these residues are not important for catalytic function and do not define the Q_L substrate binding site.

Oxygen as Acceptor

Like most bacteria, Escherichia coli has a flexible and branched respiratory chain that enables the prokaryote to live under a variety of environmental conditions, from highly aerobic to completely anaerobic. In general, the bacterial respiratory chain is composed of dehydrogenases, a quinone pool, and reductases. Substrate-specific dehydrogenases transfer reducing equivalents from various donor substrates (NADH, succinate, glycerophosphate, formate, hydrogen, pyruvate, and lactate) to a quinone pool (menaquinone, ubiquinone, and dimethylmenoquinone). Then electrons from reduced quinones (quinols) are transferred by terminal reductases to different electron acceptors. Under aerobic growth conditions, the terminal electron acceptor is molecular oxygen. A transfer of electrons from quinol to O₂ is served by two major oxidoreductases (oxidases), cytochrome bo₃ encoded by cyoABCDE and cytochrome bd encoded by cydABX. Terminal oxidases of aerobic respiratory chains of bacteria, which use O₂ as the final electron acceptor, can oxidize one of two alternative electron donors, either cytochrome c or quinol. This review compares the effects of different inhibitors on the respiratory activities of cytochrome bo₃ and cytochrome bd in E. coli. It also presents a discussion on the genetics and the prosthetic groups of cytochrome bo₃ and cytochrome bd. The E. coli membrane contains three types of quinones that all have an octaprenyl side chain (C₄₀). It has been proposed that the bo₃ oxidase can have two ubiquinone-binding sites with different affinities. "WHAT'S NEW" IN THE REVISED ARTICLE: The revised article comprises additional information about subunit composition of cytochrome bd and its role in bacterial resistance to nitrosative and oxidative stresses. Also, we present the novel data on the electrogenic function of appBCX-encoded cytochrome bd-II, a second bd-type oxidase that had been thought not to contribute to generation of a proton motive force in E. coli, although its spectral properties closely resemble those of cydABX-encoded cytochrome bd.

Molecular modeling of antibodies for the treatment of TNFα-related immunological diseases

Therapeutic monoclonal antibodies (mAbs) have high efficacy in treating TNF α-related immunological diseases. Other than neutralizing TNF α, these IgG1 antibodies exert Fc receptor-mediated effector functions such as the complement-dependent cytotoxicity (CDC) and antibody-dependent cell cytotoxicity (ADCC). The crystallizable fragment (Fc) of these IgG1 contains a single glycosylation site at Asn 297/300 that is essential for the CDC and ADCC. Glycosylated antibodies lacking core fucosylation showed an improved ADCC. However, no structural data are available concerning the ligand-binding interaction of these mAbs used in TNF α-related diseases and the role of the fucosylation. We therefore used comparative modeling for generating complete 3D mAb models that include the antigen-binding fragment (Fab) portions of infliximab, complexed with TNF α (4G3Y.pdb), the Fc region of the human IGHG1 fucosylated (3SGJ) and afucosylated (3SGK) complexed with the Fc receptor subtype Fcγ RIIIA, and the Fc region of a murine immunoglobulin (1IGT). After few thousand steps of energy minimization on the resulting 3D mAb models, minimized final models were used to quantify interactions occurring between Fcγ RIIIA and the fucosylated/afucosylated Fc fragments. While fucosylation does not affect Fab-TNF α interactions, we found that in the absence of fucosylation the Fc-mAb domain and Fcγ RIIIA are closer and new strong interactions are established between G129 of the receptor and S301 of the Chimera 2 Fc mAb; new polar interactions are also established between the Chimera 2 Fc residues Y299, N300, and S301 and the Fcγ RIIIA residues K128, G129, R130, and R155. These data help to explain the reduced ADCC observed in the fucosylated mAbs suggesting the specific AA residues involved in binding interactions.