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.

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 inwardfacing 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. Multidrug resistance (MDR) to anticancer drugs is, at the moment, a major contributor to chemotherapy failure 1. In cancer, one of the most significant MDR mechanisms results from the overexpression of P-glycoprotein, a membrane efflux pump (P-gp, ABCB1) 2. Thus, a deeper understanding on P-gp substrate specificity and efflux-related signal transmission mechanism remains crucial for the development of more potent and selective compounds able to modulate drug efflux 3. P-glycoprotein exports a broad range of structurally unrelated compounds through an ATP-dependent mechanism 4. P-gp is organized in two homologous functional units (N-and C-terminal halves) with a pseudo-2-fold symmetry. Each halve comprises one transmembrane domain (TMD), formed by six transmembrane α-helices (TMHs), and one cytoplasmic nucleotide-binding domain (NBD). Both N-and C-terminal halves are connected by a small peptide sequence (the "linker"; residues 627-688) 5,6. The TMHs are directly linked to the respective NBD by the intracellular loops, through the functional TM helices 6 (NBD1) and 12 (NBD2) and non-covalently by short intracellular coupling helices (ICHs), located between the structural TMHs 2/3 (ICH1-NBD1), 4/5 (ICH2-NBD2), 8/9 (ICH3-NBD2) and 10/11 (ICH4-NBD1) 7,8 (Fig. 1). These ICHs were found to be important for the maturation and folding of the P-gp transpo...

P‐glycoprotein (P‐gp) has been considered an important molecular target in the reversal of multidrug resistance (MDR). As such, the development of P‐gp modulators able to restore drug sensitivity in resistant cells is still considered one of the most promising strategies for overcoming MDR. Since the identification of the P‐gp's role in MDR, several studies have been performed in order to develop effective P‐gp modulators and understand the efflux mechanism. However, no efflux modulator is still clinically available for treating multidrug‐resistant cancers. Nevertheless, recent experimental studies suggest that MDR can be surpassed by targeting a specific region within the ABC transporter structure rather than the polyspecific drug‐binding pocket. This article will focus on the information available about this new target region and on a brief overview of which scaffolds would be suitable for modulating P‐gp at this new location. WIREs Comput Mol Sci 2017, 7:e1316. doi: 10.1002/wcms.1316 This article is categorized under: Structure and Mechanism > Computational Biochemistry and Biophysics

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).

The occluded conformation suggested in a recent article that revealed a new inward-17 facing conformation for the human P-glycoprotein may not represent the closing of a gate region 18 but instead an artifact derived from lateral compression in a too small sized nanodisc, used to 19 stabilize the transmembrane domains of the transporter. 20 21

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.

No abstract

The modulation of drug efflux by P-glycoprotein (P-gp, ABCB1) represents one of the most promising approaches to overcome multidrug resistance (MDR) in cancer cells, however the mechanisms of drug specificity and signal-transmission are still poorly understood, hampering the development of more selective and efficient P-gp modulators. 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), which are 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.