The switching mechanism of the mitochondrial ADP/ATP carrier explored by free-energy landscapes

Nucleotide substitutions at the p.Gly117 and p.Thr180 mutational hot-spots of SKI alter molecular dynamics and may affect cell cycle

Heterozygous deleterious variants in SKI cause Shprintzen-Goldberg Syndrome, which is mainly characterized by craniofacial features, neurodevelopmental disorder and thoracic aorta dilatations/aneurysms. The encoded protein is a member of the transforming growth factor beta signaling. Paucity of reported studies exploring the SGS molecular pathogenesis hampers disease recognition and clinical interpretation of private variants. Here, the unpublished c.349G>A, p.[Gly117Ser] and the recurrent c.539C>T, p.[Thr180Met] SKI variants were studied combining in silico and in vitro approach. 3D comparative modeling and calculation of the interaction energy predicted that both variants alter the SKI tertiary protein structure and its interactions. Computational data were functionally corroborated by the demonstration of an increase of MAPK phosphorylation levels and alteration of cell cycle in cells expressing the mutant SKI. Our findings confirmed the effects of SKI variants on MAPK and opened the path to study the role of perturbations of the cell cycle in SGS.

Layered mechanisms regulating the human mitochondrial NAD+ transporter SLC25A51

SLC25A51 is the primary mitochondrial NAD+ transporter in humans and controls many local reactions by mediating the influx of oxidized NAD+. Intriguingly, SLC25A51 lacks several key features compared with other members in the mitochondrial carrier family, thus its molecular mechanism has been unclear. A deeper understanding would shed light on the control of cellular respiration, the citric acid cycle, and free NAD+ concentrations in mammalian mitochondria. This review discusses recent insights into the transport mechanism of SLC25A51, and in the process highlights a multitiered regulation that governs NAD+ transport. The aspects regulating SLC25A51 import activity can be categorized as contributions from (1) structural characteristics of the transporter itself, (2) its microenvironment, and (3) distinctive properties of the transported ligand. These unique mechanisms further evoke compelling new ideas for modulating the activity of this transporter, as well as new mechanistic models for the mitochondrial carrier family.

Adenine nucleotide carrier protein dysfunction in human disease

Adenine nucleotide translocase (ANT) is the prototypical member of the mitochondrial carrier protein family, primarily involved in ADP/ATP exchange across the inner mitochondrial membrane. Several carrier proteins evolutionarily related to ANT, including SLC25A24 and SLC25A25, are believed to promote the exchange of cytosolic ATP-Mg2+ with phosphate in the mitochondrial matrix. They allow a net accumulation of adenine nucleotides inside mitochondria, which is essential for mitochondrial biogenesis and cell growth. In the last two decades, mutations in the heart/muscle isoform 1 of ANT (ANT1) and the ATP-Mg2+ transporters have been found to cause a wide spectrum of human diseases by a recessive or dominant mechanism. Although loss-of-function recessive mutations cause a defect in oxidative phosphorylation and an increase in oxidative stress which drives the pathology, it is unclear how the dominant missense mutations in these proteins cause human diseases. In this review, we focus on how yeast was productively used as a model system for the understanding of these dominant diseases. We also describe the relationship between the structure and function of ANT and how this may relate to various pathologies. Particularly, mutations in Aac2, the yeast homolog of ANT, were recently found to clog the mitochondrial protein import pathway. This leads to mitochondrial precursor overaccumulation stress (mPOS), characterized by the toxic accumulation of unimported mitochondrial proteins in the cytosol. We anticipate that in coming years, yeast will continue to serve as a useful model system for the mechanistic understanding of mitochondrial protein import clogging and related pathologies in humans.

FA Sliding as the Mechanism for the ANT1-Mediated Fatty Acid Anion Transport in Lipid Bilayers

Mitochondrial adenine nucleotide translocase (ANT) exchanges ADP for ATP to maintain energy production in the cell. Its protonophoric function in the presence of long-chain fatty acids (FA) is also recognized. Our previous results imply that proton/FA transport can be best described with the FA cycling model, in which protonated FA transports the proton to the mitochondrial matrix. The mechanism by which ANT1 transports FA anions back to the intermembrane space remains unclear. Using a combined approach involving measurements of the current through the planar lipid bilayers reconstituted with ANT1, site-directed mutagenesis and molecular dynamics simulations, we show that the FA anion is first attracted by positively charged arginines or lysines on the matrix side of ANT1 before moving along the positively charged protein-lipid interface and binding to R79, where it is protonated. We show that R79 is also critical for the competitive binding of ANT1 substrates (ADP and ATP) and inhibitors (carboxyatractyloside and bongkrekic acid). The binding sites are well conserved in mitochondrial SLC25 members, suggesting a general mechanism for transporting FA anions across the inner mitochondrial membrane.

In Silico Analysis of the Structural Dynamics and Substrate Recognition Determinants of the Human Mitochondrial Carnitine/Acylcarnitine SLC25A20 Transporter

The Carnitine-Acylcarnitine Carrier is a member of the mitochondrial Solute Carrier Family 25 (SLC25), known as SLC25A20, involved in the electroneutral exchange of acylcarnitine and carnitine across the inner mitochondrial membrane. It acts as a master regulator of fatty acids β-oxidation and is known to be involved in neonatal pathologies and cancer. The transport mechanism, also known as "alternating access", involves a conformational transition in which the binding site is accessible from one side of the membrane or the other. In this study, through a combination of state-of-the-art modelling techniques, molecular dynamics, and molecular docking, the structural dynamics of SLC25A20 and the early substrates recognition step have been analyzed. The results obtained demonstrated a significant asymmetry in the conformational changes leading to the transition from the c- to the m-state, confirming previous observations on other homologous transporters. Moreover, analysis of the MD simulations' trajectories of the apo-protein in the two conformational states allowed for a better understanding of the role of SLC25A20 Asp231His and Ala281Val pathogenic mutations, which are at the basis of Carnitine-Acylcarnitine Translocase Deficiency. Finally, molecular docking coupled to molecular dynamics simulations lend support to the multi-step substrates recognition and translocation mechanism already hypothesized for the ADP/ATP carrier.

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.

Function-Related Asymmetry of the Interactions between Matrix Loops and Conserved Sequence Motifs in the Mitochondrial ADP/ATP Carrier

The ADP/ATP carrier (AAC) plays a central role in oxidative metabolism by exchanging ATP and ADP across the inner mitochondrial membrane. Previous experiments have shown the involvement of the matrix loops of AAC in its function, yet potential mechanisms remain largely elusive. One obstacle is the limited information on the structural dynamics of the matrix loops. In the current work, unbiased all-atom molecular dynamics (MD) simulations were carried out on c-state wild-type AAC and mutants. Our results reveal that: (1) two ends of a matrix loop are tethered through interactions between the residue of triplet 38 (Q38, D143 and Q240) located at the C-end of the odd-numbered helix and residues of the [YF]xG motif located before the N-end of the short matrix helix in the same domain; (2) the initial progression direction of a matrix loop is determined by interactions between the negatively charged residue of the [DE]G motif located at the C-end of the short matrix helix and the capping arginine (R30, R139 and R236) in the previous domain; (3) the two chemically similar residues D and E in the highly conserved [DE]G motif are actually quite different; (4) the N-end of the M3 loop is clamped by the [DE]G motif and the capping arginine of domain 2 from the two sides, which strengthens interactions between domain 2 and domain 3; and (5) a highly asymmetric stable core exists within domains 2 and 3 at the m-gate level. Moreover, our results help explain almost all extremely conserved residues within the matrix loops of the ADP/ATP carriers from a structural point of view. Taken together, the current work highlights asymmetry in the three matrix loops and implies a close relationship between asymmetry and ADP/ATP transport.

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.

Metabolic regulation of follicular helper T cell differentiation in a mouse model of lupus

Follicular helper T (T_FH) cells are expanded in systemic lupus erythematosus (SLE), where they are required for production of high affinity autoantibodies. A better understanding of the mechanisms that regulate the differentiation of T_FH cells is critical. Naïve T cells from lupus-prone B6.NZM2410.Sle1.Sle2.Sle3 (TC) mice showed an intrinsic higher capacity to differentiate into T_FH cells. Metabolic reprogramming is a vital regulatory mechanism for T cell differentiation, but how metabolic pathways contribute to T_FH cell expansion in SLE remains elusive. Here we show that glycolysis, mTOR signaling, FAO, and the activity of complex V of the electron transport chain support T_FH lineage commitment. Blocking complex I uniquely decreased the expansion of T_FH cells from lupus-prone mice, and inhibition of some pathways had a greater effect in lupus-prone than control T_FH cells. However, blocking glutaminolysis, complex III and ADP/ATP translocase did not affect T_FH cell expansion. Together, our results identified novel intrinsic metabolic requirements for T_FH cell differentiation, and further defined the differential metabolic pathways that support the expansion of T_FH cells in lupus-prone mice. Together, our data indicates the crucial but distinct roles for metabolic pathways in T_FH cell differentiation and provide a comprehensive experimental basis for fully understanding the precise roles of distant metabolic signaling in regulating the T_FH cell differentiation.

Evidence for Non-Essential Salt Bridges in the M-Gates of Mitochondrial Carrier Proteins

Mitochondrial carriers, which transport metabolites, nucleotides, and cofactors across the mitochondrial inner membrane, have six transmembrane α-helices enclosing a translocation pore with a central substrate binding site whose access is controlled by a cytoplasmic and a matrix gate (M-gate). The salt bridges formed by the three PX[DE]XX[RK] motifs located on the odd-numbered transmembrane α-helices greatly contribute to closing the M-gate. We have measured the transport rates of cysteine mutants of the charged residue positions in the PX[DE]XX[RK] motifs of the bovine oxoglutarate carrier, the yeast GTP/GDP carrier, and the yeast NAD⁺ transporter, which all lack one of these charged residues. Most single substitutions, including those of the non-charged and unpaired charged residues, completely inactivated transport. Double mutations of charged pairs showed that all three carriers contain salt bridges non-essential for activity. Two double substitutions of these non-essential charge pairs exhibited higher transport rates than their corresponding single mutants, whereas swapping the charged residues in these positions did not increase activity. The results demonstrate that some of the residues in the charged residue positions of the PX[DE]XX[KR] motifs are important for reasons other than forming salt bridges, probably for playing specific roles related to the substrate interaction-mediated conformational changes leading to the M-gate opening/closing.

Mitochondrial transport and metabolism of the vitamin B-derived cofactors thiamine pyrophosphate, coenzyme A, FAD and NAD+ , and related diseases: A review

Multiple mitochondrial matrix enzymes playing key roles in metabolism require cofactors for their action. Due to the high impermeability of the mitochondrial inner membrane, these cofactors need to be synthesized within the mitochondria or be imported, themselves or one of their precursors, into the organelles. Transporters belonging to the protein family of mitochondrial carriers have been identified to transport the coenzymes: thiamine pyrophosphate, coenzyme A, FAD and NAD⁺ , which are all structurally similar to nucleotides and derived from different B-vitamins. These mitochondrial cofactors bind more or less tightly to their enzymes and, after having been involved in a specific reaction step, are regenerated, spontaneously or by other enzymes, to return to their active form, ready for the next catalysis round. Disease-causing mutations in the mitochondrial cofactor carrier genes compromise not only the transport reaction but also the activity of all mitochondrial enzymes using that particular cofactor and the metabolic pathways in which the cofactor-dependent enzymes are involved. The mitochondrial transport, metabolism and diseases of the cofactors thiamine pyrophosphate, coenzyme A, FAD and NAD⁺ are the focus of this review.

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.

Structural determinants of ligands recognition by the human mitochondrial basic amino acids transporter SLC25A29. Insights from molecular dynamics simulations of the c-state

In mitochondria, metabolic processes require the trafficking of solutes and organic molecules, such as amino acids. This task is accomplished by the Mitochondrial Carrier Family members (also known as SLC25), among which the SLC25A29 is responsible for the translocation of basic amino acids. In this regard, nitric oxide levels originated by the arginine mitochondrial catabolism have been shown to strongly affect cancer cells' metabolic status. Furthermore, the metabolic disease saccharopinuria has been linked to a mitochondrial dysregulation caused by a toxic intermediate of the lysine catabolism. In both cases, a reduction of the activity of SLC25A29 has been shown to ameliorate these pathological conditions. However, no detailed structural data are available on SLC25A29. In the present work, molecular modelling, docking and dynamics simulations have been employed to analyse the structural determinants of ligands recognition by SLC25A29 in the c-state. Results confirm and reinforce earlier predictions that Asn73, Arg160 and Glu161, and Arg257 represent the ligand contact points I, II, and III, respectively, and that Arg160, Trp204 and Arg257 form a stable interaction, likely critical for ligand binding and translocation. These results are discussed in view of the experimental data available for SLC25A29 and other homologous carriers of the same family.

The Interaction of Hemin, a Porphyrin Derivative, with the Purified Rat Brain 2-Oxoglutarate Carrier

The mitochondrial 2-oxoglutarate carrier (OGC), isolated and purified from rat brain mitochondria, was reconstituted into proteoliposomes to study the interaction with hemin, a porphyrin derivative, which may result from the breakdown of heme-containing proteins and plays a key role in several metabolic pathways. By kinetic approaches, on the basis of the single binding centre gated pore mechanism, we analyzed the effect of hemin on the transport rate of OGC in uptake and efflux experiments in proteoliposomes reconstituted in the presence of the substrate 2-oxoglutarate. Overall, our experimental data fit the hypothesis that hemin operates a competitive inhibition in the 0.5-10 µM concentration range. As a consequence of the OGC inhibition, the malate/aspartate shuttle might be impaired, causing an alteration of mitochondrial function. Hence, considering that the metabolism of porphyrins implies both cytoplasmic and mitochondrial processes, OGC may participate in the regulation of porphyrin derivatives availability and the related metabolic pathways that depend on them (such as oxidative phosphorylation and apoptosis). For the sake of clarity, a simplified model based on induced-fit molecular docking supported the in vitro transport assays findings that hemin was as good as 2-oxoglutarate to bind the carrier by engaging specific ionic hydrogen bond interactions with a number of key residues known for participating in the similarly located mitochondrial carrier substrate binding site.

Computational studies of the mitochondrial carrier family SLC25. Present status and future perspectives

The members of the mitochondrial carrier family, also known as solute carrier family 25 (SLC25), are transmembrane proteins involved in the translocation of a plethora of small molecules between the mitochondrial intermembrane space and the matrix. These transporters are characterized by three homologous domains structure and a transport mechanism that involves the transition between different conformations. Mutations in regions critical for these transporters’ function often cause several diseases, given the crucial role of these proteins in the mitochondrial homeostasis. Experimental studies can be problematic in the case of membrane proteins, in particular concerning the characterization of the structure–function relationships. For this reason, computational methods are often applied in order to develop new hypotheses or to support/explain experimental evidence. Here the computational analyses carried out on the SLC25 members are reviewed, describing the main techniques used and the outcome in terms of improved knowledge of the transport mechanism. Potential future applications on this protein family of more recent and advanced in silico methods are also suggested.

The Mitochondrial Carnitine Acyl-carnitine Carrier (SLC25A20): Molecular Mechanisms of Transport, Role in Redox Sensing and Interaction with Drugs

The SLC25A20 transporter, also known as carnitine acyl-carnitine carrier (CAC), catalyzes the transport of short, medium and long carbon chain acyl-carnitines across the mitochondrial inner membrane in exchange for carnitine. The 30-year story of the protein responsible for this function started with its purification from rat liver mitochondria. Even though its 3D structure is not yet available, CAC is one of the most deeply characterized transport proteins of the inner mitochondrial membrane. Other than functional, kinetic and mechanistic data, post-translational modifications regulating the transport activity of CAC have been revealed. CAC interactions with drugs or xenobiotics relevant to human health and toxicology and the response of the carrier function to dietary compounds have been discovered. Exploiting combined approaches of site-directed mutagenesis with chemical targeting and bioinformatics, a large set of data on structure/function relationships have been obtained, giving novel information on the molecular mechanism of the transport catalyzed by this protein.

ANT1 Activation and Inhibition Patterns Support the Fatty Acid Cycling Mechanism for Proton Transport

Adenine nucleotide translocase (ANT) is a well-known mitochondrial exchanger of ATP against ADP. In contrast, few studies have shown that ANT also mediates proton transport across the inner mitochondrial membrane. The results of these studies are controversial and lead to different hypotheses about molecular transport mechanisms. We hypothesized that the H⁺-transport mediated by ANT and uncoupling proteins (UCP) has a similar regulation pattern and can be explained by the fatty acid cycling concept. The reconstitution of purified recombinant ANT1 in the planar lipid bilayers allowed us to measure the membrane current after the direct application of transmembrane potential ΔΨ, which would correspond to the mitochondrial states III and IV. Experimental results reveal that ANT1 does not contribute to a basal proton leak. Instead, it mediates H⁺ transport only in the presence of long-chain fatty acids (FA), as already known for UCPs. It depends on FA chain length and saturation, implying that FA’s transport is confined to the lipid-protein interface. Purine nucleotides with the preference for ATP and ADP inhibited H⁺ transport. Specific inhibitors of ATP/ADP transport, carboxyatractyloside or bongkrekic acid, also decreased proton transport. The H⁺ turnover number was calculated based on ANT1 concentration determined by fluorescence correlation spectroscopy and is equal to 14.6 ± 2.5 s⁻¹. Molecular dynamic simulations revealed a large positively charged area at the protein/lipid interface that might facilitate FA anion’s transport across the membrane. ANT’s dual function—ADP/ATP and H⁺ transport in the presence of FA—may be important for the regulation of mitochondrial membrane potential and thus for potential-dependent processes in mitochondria. Moreover, the expansion of proton-transport modulating drug targets to ANT1 may improve the therapy of obesity, cancer, steatosis, cardiovascular and neurodegenerative diseases.

Molecular Dynamics Simulations of Mitochondrial Uncoupling Protein 2

Molecular dynamics (MD) simulations of uncoupling proteins (UCP), a class of transmembrane proteins relevant for proton transport across inner mitochondrial membranes, represent a complicated task due to the lack of available structural data. In this work, we use a combination of homology modelling and subsequent microsecond molecular dynamics simulations of UCP2 in the DOPC phospholipid bilayer, starting from the structure of the mitochondrial ATP/ADP carrier (ANT) as a template. We show that this protocol leads to a structure that is impermeable to water, in contrast to MD simulations of UCP2 structures based on the experimental NMR structure. We also show that ATP binding in the UCP2 cavity is tight in the homology modelled structure of UCP2 in agreement with experimental observations. Finally, we corroborate our results with conductance measurements in model membranes, which further suggest that the UCP2 structure modeled from ANT protein possesses additional key functional elements, such as a fatty acid-binding site at the R60 region of the protein, directly related to the proton transport mechanism across inner mitochondrial membranes.

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.

Antimicrobial Peptide TP4 Targets Mitochondrial Adenine Nucleotide Translocator 2

Tilapia piscidin (TP) 4 is an antimicrobial peptide derived from Nile tilapia (Oreochromis niloticus), which shows broad-spectrum antibacterial activity and excellent cancer-killing ability in vitro and in vivo. Like many other antimicrobial peptides, TP4 treatment causes mitochondrial toxicity in cancer cells. However, the molecular mechanisms underlying TP4 targeting of mitochondria remain unclear. In this study, we used a pull-down assay on A549 cell lysates combined with LC-MS/MS to discover that TP4 targets adenine nucleotide translocator (ANT) 2, a protein essential for adenine nucleotide exchange across the inner membrane. We further showed that TP4 accumulates in mitochondria and colocalizes with ANT2. Moreover, molecular docking studies showed that the interaction requires Phe1, Ile2, His3, His4, Ser11, Lys14, His17, Arg21, Arg24 and Arg25 residues in TP4 and key residues within the cavity of ANT2. These findings suggest a mechanism by which TP4 may induce mitochondrial dysfunction to disrupt cellular energy metabolism.

On the Detection and Functional Significance of the Protein-Protein Interactions of Mitochondrial Transport Proteins

Protein-protein assemblies are highly prevalent in all living cells. Considerable evidence has recently accumulated suggesting that particularly transient association/dissociation of proteins represent an important means of regulation of metabolism. This is true not only in the cytosol and organelle matrices, but also at membrane surfaces where, for example, receptor complexes, as well as those of key metabolic pathways, are common. Transporters also frequently come up in lists of interacting proteins, for example, binding proteins that catalyze the production of their substrates or that act as relays within signal transduction cascades. In this review, we provide an update of technologies that are used in the study of such interactions with mitochondrial transport proteins, highlighting the difficulties that arise in their use for membrane proteins and discussing our current understanding of the biological function of such interactions.

Fernie A. Plant Mitochondrial Carriers: Molecular Gatekeepers That Help to Regulate Plant Central Carbon Metabolism

The evolution of membrane-bound organelles among eukaryotes led to a highly compartmentalized metabolism. As a compartment of the central carbon metabolism, mitochondria must be connected to the cytosol by molecular gates that facilitate a myriad of cellular processes. Members of the mitochondrial carrier family function to mediate the transport of metabolites across the impermeable inner mitochondrial membrane and, thus, are potentially crucial for metabolic control and regulation. Here, we focus on members of this family that might impact intracellular central plant carbon metabolism. We summarize and review what is currently known about these transporters from in vitro transport assays and in planta physiological functions, whenever available. From the biochemical and molecular data, we hypothesize how these relevant transporters might play a role in the shuttling of organic acids in the various flux modes of the TCA cycle. Furthermore, we also review relevant mitochondrial carriers that may be vital in mitochondrial oxidative phosphorylation. Lastly, we survey novel experimental approaches that could possibly extend and/or complement the widely accepted proteoliposome reconstitution approach.

Metadynamics to Enhance Sampling in Biomolecular Simulations

Molecular dynamics is a powerful simulation method to provide detailed atomic-scale insight into a range of biological processes including protein folding, biochemical reactions, ligand binding, and many others. Over the last several decades, enhanced sampling methods have been developed to address the large separation in time scales between a molecular dynamics simulation (usually microseconds or shorter) and the time scales of biological processes (often orders of magnitude longer). This chapter specifically focuses on the metadynamics family of methods, which achieves enhanced sampling through the introduction of a history-dependent bias potential that is based on one or more slow degrees of freedom, called collective variables. We introduce the method and its recent variants related to biomolecular studies and then discuss frontier areas of the method. A large part of this chapter is devoted to helping new users of the method understand how to choose metadynamics parameters properly and apply the method to their system of interest.

Molecular Basis of the Glucose Transport Mechanism in Plants

The SWEET family belongs to a class of transporters in plants that undergoes large conformational changes to facilitate transport of sugar molecules across the cell membrane (SWEET, Sugars Will Eventually Be Exported Transporter). However, the structures of their functionally relevant conformational states in the transport cycle have not been reported. In this study, we have characterized the conformational dynamics and complete transport cycle of glucose in the OsSWEET2b transporter using extensive molecular dynamics simulations. Using Markov state models, we estimated the free energy barrier associated with different states as well as for the glucose transport mechanism. SWEETs undergo a structural transition to outward-facing (OF), occluded (OC), and inward-facing (IF) and strongly support an alternate access transport mechanism. The glucose diffuses freely from outside to inside the cell without causing major conformational changes which means that the conformations of glucose unbound and bound snapshots are exactly the same for OF, OC, and IF states. We identified a network of hydrophobic core residues at the center of the transporter that restricts the glucose entry to the cytoplasmic side and acts as an intracellular hydrophobic gate. The mechanistic predictions from molecular dynamics simulations are validated using site-directed mutagenesis experiments. Our simulation also revealed hourglass-like intermediate states making the pore radius narrower at the center. This work provides new fundamental insights into how substrate-transporter interactions actively change the free energy landscape of the transport cycle to facilitate enhanced transport activity.

Mutagenesis Reveals That the OsPPa6 Gene Is Required for Enhancing the Alkaline Tolerance in Rice

Alkaline stress (AS) is one of the abiotic stressful factors limiting plant's growth and development. Inorganic pyrophosphatase is usually involved in a variety of biological processes in plant in response to the abiotic stresses. Here, to clarify the responsive regulation of inorganic pyrophosphatase in rice under AS, the mutagenesis of the OsPPa6 gene encoding an inorganic pyrophosphatase in rice cv. Kitaake (Oryza sativa L. ssp. japonica) was performed by the CRISPR/Cas9 system. Two homozygous independent mutants with cas9-free were obtained by continuously screening. qPCR reveals that the OsPPa6 gene was significantly induced by AS, and the mutagenesis of the OsPPa6 gene apparently delayed rice's growth and development, especially under AS. Measurements demonstrate that the contents of pyrophosphate in the mutants were higher than those in the wild type under AS, however, the accumulation of inorganic phosphate, ATP, chlorophyll, sucrose, and starch in the mutants were decreased significantly, and the mutagenesis of the OsPPa6 gene remarkably lowered the net photosynthetic rate of rice mutants, thus reducing the contents of soluble sugar and proline, but remarkably increasing MDA, osmotic potentials and Na⁺/K⁺ ratio in the mutants under AS. Metabonomics measurement shows that the mutants obviously down-regulated the accumulation of phosphorylcholine, choline, anthranilic acid, apigenin, coniferol and dodecanoic acid, but up-regulated the accumulation of L-valine, alpha-ketoglutarate, phenylpyruvate and L-phenylalanine under AS. This study suggests that the OsPPa6 gene is an important osmotic regulatory factor in rice, and the gene-editing of CRISPR/Cas9-guided is an effective method evaluating the responsive regulation of the stress-induced gene, and simultaneously provides a scientific support for the application of the gene encoding a soluble inorganic pyrophosphatase in molecular breeding.

Molecular dynamics simulations on apo ADP/ATP carrier shed new lights on the featured motif of the mitochondrial carriers

The ADP/ATP carrier (AAC) is a transporter responsible for the equal molar exchange of cytosolic ADP and ATP synthesized within mitochondrial matrix across the mitochondrial membrane. Its primary structure consists of three homologous repeats, and each repeat contains a conserved motif that is shared by all members of the mitochondrial carrier family (MCF). Although these MCF motif residues cluster together in the crystal structure of AAC, detailed analyses on the interactions among the motif residues are still limited. In the present study, all-atom molecular dynamics (MD) simulations of up to 10 μs have been carried out on AAC, and interactions and structural dynamics of the MCF motif residues have been specifically investigated. Our simulations have revealed: i) a very asymmetrical electrostatic network at the bottom of the pocket of apo AAC, ii) the asymmetrical interactions between the Pro kink region and the [YWF][KR] G motif in three repeats, iii) the role of the conserved Arg residues in stabilizing the C-ends of the odd-numbered helices, iv) the structural change of the [YWF][KR] G motif and its potential involvement in substrate translocation process. Our results highlight the asymmetry of the MCF residues in the three repeats, which might contribute to the ability of the carriers to transport the asymmetrical substrates. Our observations provide microscopic basis for further research on the translocation mechanism of mitochondrial carriers.

Perturbations of Native Membrane Protein Structure in Alkyl Phosphocholine Detergents: A Critical Assessment of NMR and Biophysical Studies

Membrane proteins perform a host of vital cellular functions. Deciphering the molecular mechanisms whereby they fulfill these functions requires detailed biophysical and structural investigations. Detergents have proven pivotal to extract the protein from its native surroundings. Yet, they provide a milieu that departs significantly from that of the biological membrane, to the extent that the structure, the dynamics, and the interactions of membrane proteins in detergents may considerably vary, as compared to the native environment. Understanding the impact of detergents on membrane proteins is, therefore, crucial to assess the biological relevance of results obtained in detergents. Here, we review the strengths and weaknesses of alkyl phosphocholines (or foscholines), the most widely used detergent in solution-NMR studies of membrane proteins. While this class of detergents is often successful for membrane protein solubilization, a growing list of examples points to destabilizing and denaturing properties, in particular for α-helical membrane proteins. Our comprehensive analysis stresses the importance of stringent controls when working with this class of detergents and when analyzing the structure and dynamics of membrane proteins in alkyl phosphocholine detergents.

Atomistic modeling of alternating access of a mitochondrial ADP/ATP membrane transporter with molecular simulations

The mitochondrial ADP/ATP carrier (AAC) is a membrane transporter that exchanges a cytosolic ADP for a matrix ATP. Atomic structures in an outward-facing (OF) form which binds an ADP from the intermembrane space have been solved by X-ray crystallography, and revealed their unique pseudo three-fold symmetry fold which is qualitatively different from pseudo two-fold symmetry of most transporters of which atomic structures have been solved. However, any atomic-level information on an inward-facing (IF) form, which binds an ATP from the matrix side and is fixed by binding of an inhibitor, bongkrekic acid (BA), is not available, and thus its alternating access mechanism for the transport process is unknown. Here, we report an atomic structure of the IF form predicted by atomic-level molecular dynamics (MD) simulations of the alternating access transition with a recently developed accelerating technique. We successfully obtained a significantly stable IF structure characterized by newly formed well-packed and -organized inter-domain interactions through the accelerated simulations of unprecedentedly large conformational changes of the alternating access without a prior knowledge of the target protein structure. The simulation also shed light on an atomistic mechanism of the strict transport selectivity of adenosine nucleotides over guanosine and inosine ones. Furthermore, the IF structure was shown to bind ATP and BA, and thus revealed their binding mechanisms. The present study proposes a qualitatively novel view of the alternating access of transporters having the unique three-fold symmetry in atomic details and opens the way for rational drug design targeting the transporter in the dynamic functional cycle.

Modelling the free energy profile of the mitochondrial ADP/ATP carrier

The mitochondrial ADP/ATP carrier catalyses the equimolar exchange of adenosine di- and tri-phosphates. It operates by an alternating access mechanism in which a single substrate-binding site is made available either to the mitochondrial matrix or the intermembrane space through conformational changes. These changes are prevented in the absence of substrate by a large energy barrier due to the need for sequential disruption and formation of a matrix and cytoplasmic salt bridge network that are located on either side of the central cavity. In analogy to enzyme catalysis, substrate lowers the energy barrier by binding tighter in the intermediate state. Here we provide an in-silico kinetic model that captures the free energy profile of these conformational changes and treats the carrier as a nanomachine moving stochastically from the matrix to cytoplasmic conformation under the influence of thermal energy. The model reproduces the dependency of experimentally determined k_cat and K_M values on the cytoplasmic network strength with good quantitative accuracy, implying that it captures the transport mechanism and can provide a framework to understand the structure-function relationships of this class of transporter. The results show that maximum transport occurs when the interaction energies of the cytoplasmic network, matrix network and substrate binding are approximately equal such that the energy barrier is minimized. Consequently, the model predicts that there will be other interactions in addition to those of the cytoplasmic network that stabilise the matrix conformation of the ADP/ATP carrier.

Lipid-Loving ANTs: Molecular Simulations of Cardiolipin Interactions and the Organization of the Adenine Nucleotide Translocase in Model Mitochondrial Membranes

The exchange of ADP and ATP across the inner mitochondrial membrane is a fundamental cellular process. This exchange is facilitated by the adenine nucleotide translocase, the structure and function of which are critically dependent on the signature phospholipid of mitochondria, cardiolipin (CL). Here we employ multiscale molecular dynamics simulations to investigate CL interactions within a membrane environment. Using simulations at both coarse-grained and atomistic resolutions, we identify three CL binding sites on the translocase, in agreement with those seen in crystal structures and inferred from nuclear magnetic resonance measurements. Characterization of the free energy landscape for lateral lipid interaction via potential of mean force calculations demonstrates the strength of interaction compared to those of binding sites on other mitochondrial membrane proteins, as well as their selectivity for CL over other phospholipids. Extending the analysis to other members of the family, yeast Aac2p and mouse uncoupling protein 2, suggests a degree of conservation. Simulation of large patches of a model mitochondrial membrane containing multiple copies of the translocase shows that CL interactions persist in the presence of protein–protein interactions and suggests CL may mediate interactions between translocases. This study provides a key example of how computational microscopy may be used to shed light on regulatory lipid–protein interactions.

UCP1 - A sophisticated energy valve

This review focuses on the biochemical work of UCP1 starting from the early observation by Ricquier and Kader in 1976. We entered this field in 1980 with the isolation of native UCP1 and then reported the amino acid sequence structure discovering a strong homology to the ADP/ATP carrier. With the isolated native UCP1 we studied structural and functional features, in particular the complex characteristics of nucleotide binding. A strong pH dependence of binding and herein the differences between diphopho- and triphopho-nucleotides were observed, resulting in the identification of residues which control binding site access by their H⁺ dissociation. Newly synthesized fluorescent nucleotide derivatives provided tools to determine a two state nucleotide binding in line with loose and tight UCP1 conformations and H⁺ transport inhibition. The slow transition between these states were a notable feature. The reconstitution of isolated UCP1 in vesicles demonstrated that UCP1 protein is in fact the uncoupling factor and not only a nucleotide controlled regulator. The H⁺ transport was shown to be electrophoretic with a linear relation to the membrane potential. The dependence of H⁺ transport on fatty acids (FA) was characterized and is elaborated here with a view of the experimental conditions of other research groups which had different views of the role of FA in H⁺ transport. Furthermore, to explain the contrast of the FA - nucleotide competition between mitochondria and reconstituted system, indirect paths for FA to relieve the inhibition in mitochondria are here proposed, such as a FA induced upward pH shift and a FA induced increase of cardiolipin level around UCP1 since cardiolipin has been found by us to relieve nucleotide binding on isolated UCP1. Recently reported patch clamp results on mitoplasts led to a reformulation of the H⁺ transport mechanism of FA in UCP1 in which bound FA shuttles with the carboxyl group between the two membrane sides along the translocation channel outward as FA^- and inward as FA^-H⁺. We propose here a modified version, where FA forms an immobile prosthetic group surrounded by the inner and outer gate of the H⁺ translocation channel. By alternating opening of the gates FA takes up H⁺ from the cytosol side and releases H⁺ to the matrix.

Dataset of the AAC2 conformations in the c-, intermediate- and m-states obtained from free-energy simulations

The data reported herein are related to the article entitled: "The switching mechanism of the mitochondrial ADP/ATP carrier explored by free-energy landscapes" (Pietropaolo et al., 2016) [1]. We report the coordinates of the ADP/ATP carrier (AAC2) in the presence and absence of adenine and guanine nucleotides in the c-, intermediate- and m-states obtained from the free-energy simulations and corresponding to the free-energy minima.