Functional and structural impact of 10 ACADM missense mutations on human medium chain acyl-Coa dehydrogenase.

Medium chain acyl-CoA dehydrogenase (MCAD) deficiency (MCADD) is associated with ACADM gene mutations, leading to an impaired function and/or structure of MCAD. Importantly, after import into the mitochondria, MCAD must incorporate a molecule of flavin adenine dinucleotide (FAD) per subunit and assemble into tetramers. However, the effect of MCAD amino acid substitutions on FAD incorporation has not been investigated. Herein, the commonest MCAD variant (p.K304E) and 11 additional rare variants (p.Y48C, p.R55G, p.A88P, p.Y133C, p.A140T, p.D143V, p.G224R, p.L238F, p.V264I, p.Y372N, and p.G377V) were functionally and structurally characterized. Half of the studied variants presented a FAD content <65 % compared to the wild-type. Most of them were recovered as tetramers, except the p.Y372N (mainly as dimers). No correlation was found between the levels of tetramers and FAD content. However, a correlation between FAD content and the cofactor's affinity, proteolytic stability, thermostability, and thermal inactivation was established. We showed that the studied amino acid changes in MCAD may alter the substrate chain-length dependence and the interaction with electron-transferring-flavoprotein (ETF) necessary for a proper functioning electron transfer thus adding additional layers of complexity to the pathological effect of ACADM missense mutations. Although the majority of the variant MCADs presented an impaired capacity to retain FAD during their synthesis, some of them were structurally rescued by cofactor supplementation, suggesting that in the mitochondrial environment the levels and activity of those variants may be dependent of FAD's availability thus contributing for the heterogeneity of the MCADD phenotype found in patients presenting the same genotype.

Probing the Allosteric Modulation of P-Glycoprotein: A Medicinal Chemistry Approach Toward the Identification of Noncompetitive P-Gp Inhibitors.

A medicinal chemistry approach combining in silico and in vitro methodologies was performed aiming at identifying and characterizing putative allosteric drug-binding sites (aDBSs) at the interface of the transmembrane- and nucleotide-binding domains (TMD-NBD) of P-glycoprotein. Two aDBSs were identified, one in TMD1/NBD1 and another one in TMD2/NBD2, by means of in silico fragment-based molecular dynamics and characterized in terms of size, polarity, and lining residues. From a small library of thioxanthone and flavanone derivatives, experimentally described to bind at the TMD-NBD interfaces, several compounds were identified to be able to decrease the verapamil-stimulated ATPase activity. An IC50 of 81 ± 6.6 µM is reported for a flavanone derivative in the ATPase assays, providing evidence for an allosteric efflux modulation in P-glycoprotein. Molecular docking and molecular dynamics gave additional insights on the binding mode on how flavanone derivatives may act as allosteric inhibitors.

Long-range communication between transmembrane- and nucleotide-binding domains does not depend on drug binding to mutant P-glycoprotein.

In this study, the impact of four P-gp mutations (G185V, G830V, F978A and ΔF335) on drug-binding and efflux-related signal-transmission mechanism was comprehensively evaluated in the presence of ligands within the drug-binding pocket (DBP), experimentally related with changes in their drug efflux profiles. The severe repacking of the transmembrane helices (TMH), induced by mutations and exacerbated by the presence of ligands, indicates that P-gp is sensitive to perturbations in the transmembrane region. Alterations on drug-binding were also observed as a consequence of the TMH repacking, but were not always correlated with alterations on ligands binding mode and/or binding affinity. Finally, and although all P-gp variants holo systems showed considerable changes in the intracellular coupling helices/nucleotide-binding domain (ICH-NBD) interactions, they seem to be primarily induced by the mutation itself rather than by the presence of ligands within the DBP. The data further suggest that the changes in drug efflux experimentally reported are mostly related with changes on drug specificity rather than effects on signal-transmission mechanism. We also hypothesize that an increase in the drug-binding affinity may also be related with the decreased drug efflux, while minor changes in binding affinities are possibly related with the increased drug efflux observed in transfected cells.Communicated by Ramaswamy H. Sarma.

Theoretical insights on helix repacking as the origin of P-glycoprotein promiscuity.

P-glycoprotein (P-gp, ABCB1) overexpression is, currently, one of the most important multidrug resistance (MDR) mechanisms in tumor cells. Thus, modulating drug efflux by P-gp has become one of the most promising approaches to overcome MDR in cancer. Yet, more insights on the molecular basis of drug specificity and efflux-related signal transmission mechanism between the transmembrane domains (TMDs) and the nucleotide binding domains (NBDs) are needed to develop molecules with higher selectivity and efficacy. Starting from a murine P-gp crystallographic structure at the inward-facing conformation (PDB ID: 4Q9H), we evaluated the structural quality of the herein generated human P-gp homology model. This initial human P-gp model, in the presence of the "linker" and inserted in a suitable lipid bilayer, was refined through molecular dynamics simulations and thoroughly validated. The best human P-gp model was further used to study the effect of four single-point mutations located at the TMDs, experimentally related with changes in substrate specificity and drug-stimulated ATPase activity. Remarkably, each P-gp mutation is able to induce transmembrane α-helices (TMHs) repacking, affecting the drug-binding pocket volume and the drug-binding sites properties (e.g. volume, shape and polarity) finally compromising drug binding at the substrate binding sites. Furthermore, intracellular coupling helices (ICH) also play an important role since changes in the TMHs rearrangement are shown to have an impact in residue interactions at the ICH-NBD interfaces, suggesting that identified TMHs repacking affect TMD-NBD contacts and interfere with signal transmission from the TMDs to the NBDs.

Author Correction: Structure-function relationships in ABCG2: insights from molecular dynamics simulations and molecular docking studies.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Structure-function relationships in ABCG2: insights from molecular dynamics simulations and molecular docking studies.

Efflux pumps of the ATP-binding cassette transporters superfamily (ABC transporters) are frequently involved in the multidrug-resistance (MDR) phenomenon in cancer cells. Herein, we describe a new atomistic model for the MDR-related ABCG2 efflux pump, also named breast cancer resistance protein (BCRP), based on the recently published crystallographic structure of the ABCG5/G8 heterodimer sterol transporter, a member of the ABCG family involved in cholesterol homeostasis. By means of molecular dynamics simulations and molecular docking, a far-reaching characterization of the ABCG2 homodimer was obtained. The role of important residues and motifs in the structural stability of the transporter was comprehensively studied and was found to be in good agreement with the available experimental data published in literature. Moreover, structural motifs potentially involved in signal transmission were identified, along with two symmetrical drug-binding sites that are herein described for the first time, in a rational attempt to better understand how drug binding and recognition occurs in ABCG2 homodimeric transporters.

Unveiling the Pathogenic Molecular Mechanisms of the Most Common Variant (p.K329E) in Medium-Chain Acyl-CoA Dehydrogenase Deficiency by in Vitro and in Silico Approaches.

Medium-chain acyl-CoA dehydrogenase deficiency (MCADD) is the most common genetic disorder affecting the mitochondrial fatty acid ß-oxidation pathway. The mature and functional form of human MCAD (hMCAD) is a homotetramer assembled as a dimer of dimers (monomers A/B and C/D). Each monomer binds a FAD cofactor, necessary for the enzyme's activity. The most frequent mutation in MCADD results from the substitution of a lysine with a glutamate in position 304 of mature hMCAD (p.K329E in the precursor protein). Here, we combined in vitro and in silico approaches to assess the impact of the p.K329E mutation on the protein's structure and function. Our in silico results demonstrated for the first time that the p.K329E mutation, despite lying at the dimer-dimer interface and being deeply buried inside the tetrameric core, seems to affect the tetramer surface, especially the ß-domain that forms part of the catalytic pocket wall. Additionally, the molecular dynamics data indicate a stronger impact of the mutation on the protein's motions in dimer A/B, while dimer C/D remains similar to the wild type. For dimer A/B, severe disruptions in the architecture of the pockets and in the FAD and octanoyl-CoA binding affinities were also observed. The presence of unaffected pockets (C/D) in the in silico studies may explain the decreased enzymatic activity determined for the variant protein (46% residual activity). Moreover, the in silico structural changes observed for the p.K329E variant protein provide an explanation for the structural instability observed experimentally, namely, the disturbed oligomeric profile, thermal stability, and conformational flexibility, with respect to the wild-type.

Insights into Medium-chain Acyl-CoA Dehydrogenase Structure by Molecular Dynamics Simulations.

The medium-chain acyl-CoA dehydrogenase (MCAD) is a mitochondrial enzyme that catalyzes the first step of mitochondrial fatty acid ß-oxidation (mFAO) pathway. Its deficiency is the most common genetic disorder of mFAO. Many of the MCAD disease-causing variants, including the most common p.K304E variant, show loss of function due to protein misfolding. Herein, we used molecular dynamics simulations to provide insights into the structural stability and dynamic behavior of MCAD wild-type (MCADwt) and validate a structure that would allow reliable new studies on its variants. Our results revealed that in both proteins the flavin adenine dinucleotide (FAD) has an important structural role on the tetramer stability and also in maintaining the volume of the enzyme catalytic pockets. We confirmed that the presence of substrate changes the dynamics of the catalytic pockets and increases FAD affinity. A comparison between the porcine MCADwt (pMCADwt) and human MCADwt (hMCADwt) structures revealed that both proteins are essentially similar and that the reversion of the double mutant E376G/T255E of hMCAD enzyme does not affect the structure of the protein neither its behavior in simulation. Our validated hMCADwt structure is crucial for complementing and accelerating the experimental studies aiming for the discovery and development of potential stabilizers of MCAD variants as candidates for the treatment of MCAD deficiency (MCADD).