Molecular dynamics simulations and membrane protein structure quality

Despite a growing repertoire of membrane protein structures (currently ∼120 unique structures), considerations of low resolution and crystallization in the absence of a lipid bilayer require the development of techniques to assess the global quality of membrane protein folds. This is also the case for assessment of, e.g. homology models of human membrane proteins based on structures of (distant) bacterial homologues. Molecular dynamics (MD) simulations may be used to help evaluate the quality of a membrane protein structure or model. We have used a structure of the bacterial ABC transporter MsbA which has the correct transmembrane helices but an incorrect handedness and topology of their packing to test simulation methods of quality assessment. An MD simulation of the MsbA model in a lipid bilayer is compared to a simulation of another bacterial ABC transporter, BtuCD. The latter structure has demonstrated good conformational stability in the same bilayer environment and over the same timescale (20 ns) as for the MsbA model simulation. A number of comparative analyses of the two simulations were performed to assess changes in the structural integrity of each protein. The results show a significant difference between the two simulations, chiefly due to the dramatic structural deformations of MsbA. We therefore propose that MD could become a useful quality control tool for membrane protein structural biology. In particular, it provides a way in which to explore the global conformational stability of a model membrane protein fold.     

Characterization of the Walker A motif of MsbA using site-directed spin labeling electron paramagnetic resonance spectroscopy

MsbA is an ABC transporter that transports lipid A across the inner membrane of Gram-negative bacteria such as Escherichia coli. Without functional MsbA present, bacterial cells accumulate a toxic amount of lipid A within their inner membranes. A crystal structure of MsbA was recently obtained that provides an excellent starting point for functional dynamics studies in membranes [Chang and Roth (2001) Science 293, 1793-1800]. Although a structure of MsbA is now available, several functionally important motifs common to ABC transporters are unresolved in the crystal structure. The Walker A domain, one of the ABC transporter consensus motifs that is directly involved in ATP binding, is located within a large unresolved region of the MsbA ATPase domain. Site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy is a powerful technique for characterizing local areas within a large protein structure in addition to detecting and following changes in local structure due to dynamic interactions. MsbA reconstituted into lipid membranes has been evaluated by EPR spectroscopy, and it has been determined that the Walker A domain forms an alpha-helical structure, which is consistent with the structure of this motif observed in other crystallized ABC transporters. In addition, the interaction of the Walker A residues with ATP before, during, and after hydrolysis was followed using SDSL EPR spectroscopy in order to identify the residues directly involved in substrate binding and hydrolysis.     

Recognition of sulfonylurea receptor (ABCC8/9) ligands by the multidrug resistance transporter P-glycoprotein (ABCB1): functional similarities based on common structural features between two multispecific ABC proteins

ATP-sensitive K(+) (K(ATP)) channels are the target of a number of pharmacological agents, blockers like hypoglycemic sulfonylureas and openers like the hypotensive cromakalim and diazoxide. These agents act on the channel regulatory subunit, the sulfonylurea receptor (SUR), which is an ABC protein with homologies to P-glycoprotein (P-gp). P-gp is a multidrug transporter expressed in tumor cells and in some healthy tissues. Because these two ABC proteins both exhibit multispecific recognition properties, we have tested whether SUR ligands could be substrates of P-gp. Interaction with P-gp was assayed by monitoring ATPase activity of P-gp-enriched vesicles. The blockers glibenclamide, tolbutamide, and meglitinide increased ATPase activity, with a rank order of potencies that correlated with their capacity to block K(ATP) channels. P-gp ATPase activity was also increased by the openers SR47063 (a cromakalim analog), P1075 (a pinacidil analog), and diazoxide. Thus, these molecules bind to P-gp (although with lower affinities than for SUR) and are possibly transported by P-gp. Competition experiments among these molecules as well as with typical P-gp substrates revealed a structural similarity between drug binding domains in the two proteins. To rationalize the observed data, we addressed the molecular features of these proteins and compared structural models, computerized by homology from the recently solved structures of murine P-gp and bacterial ABC transporters MsbA and Sav1866. Considering the various residues experimentally assigned to be involved in drug binding, we uncovered several hot spots, which organized spatially in two main binding domains, selective for SR47063 and for glibenclamide, in matching regions of both P-gp and SUR.     

Coarse grain lipid-protein molecular interactions and diffusion with MsbA flippase

Coarse-grained (CG) modeling has proven effective for simulating lipid bilayer dynamics on scales of biological interest. Modeling the dynamics of flexible membrane proteins within the bilayer, on the other hand, poses a considerable challenge due to the complexity of the folding or conformational landscape. In the present work, the multiscale coarse-graining method is applied to atomistic peptide-lipid "soup" simulations to develop a general set of CG protein-lipid interaction potentials. The reduced model was constructed to be compatible with recent solvent-free CG models developed for protein-protein folding and lipid-lipid model bilayer interactions. The utility of the force field was demonstrated by molecular dynamics simulation of the MsbA ABC transporter in a mixed DOPC/DOPE bilayer. An elastic network was parameterized to restrain the MsbA dimer in its open, closed and hydrolysis intermediate conformations and its impact on domain flexibility was examined. Conformational stability enabled long-time dynamics simulation of MsbA freely diffusing in a 25 nm membrane patch. Three-dimensional density analysis revealed that a shell of weakly bound "annular lipids" solvate the membrane accessible surface of MsbA and its internal substrate-binding chamber. The annular lipid binding modes, along with local perturbations in head group structure, are a function of the orientation of grooves formed between transmembrane helices and may influence the alternating access mechanism of substrate entry and translocation.     

Extending the structure of an ABC transporter to atomic resolution: modeling and simulation studies of MsbA

Molecular modeling and simulation approaches have been use to generate a complete model of the prokaryotic ABC transporter MsbA from Escherichia coli, starting from the low-resolution structure-based Calpha trace (PDB code 1JSQ). MsbA is of some biomedical interest as it is homologous to mammalian transporters such as P-glycoprotein and TAP. The quality of the MsbA model is assessed using a combination of molecular dynamics simulations and static structural analysis. These results suggest that the approach adopted for MsbA may be of general utility for generating all atom models from low-resolution crystal structures of membrane proteins. Molecular dynamics simulations of the MsbA model inserted in a fully solvated octane slab (a membrane mimetic environment) reveal that while the monomer is relatively stable, the dimer is unstable and undergoes significant conformational drift on a nanosecond time scale. This suggests that the MsbA crystal dimer may not correspond to the MsbA dimer in vivo. An alternative model of the dimer is discussed in the context of available experimental data.     

Resting state conformation of the MsbA homodimer as studied by site-directed spin labeling

MsbA is the ABC transporter for lipid A and is found in the inner membranes of Gram-negative bacteria such as Escherichia coli. Without MsbA present, bacterial cells accumulate a toxic amount of lipid A within their inner membranes. A crystal structure of MsbA was recently obtained that provides an excellent starting point for functional dynamics studies in membranes [Chang, and Roth (2001) Science 293, 1793-1800]. Although a structure of MsbA is now available, many questions remain concerning its mechanism of transport. Site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy is a powerful approach for characterizing local areas within a large protein structure in addition to detecting and following changes in local structure due to dynamic interactions within a protein. The quaternary structure of the resting state of the MsbA homodimer reconstituted into lipid membranes has been evaluated by SDSL EPR spectroscopy and chemical cross-linking techniques. SDSL and cross-linking results are consistent with the controversial resting state conformation of the MsbA homodimer found in the crystal structure, with the tips of the transmembrane helices forming a dimer interface. The position of MsbA in the membrane bilayer along with the relative orientation of the transmembrane helical bundles with respect to one another has been determined. Characterization of the resting state of the MsbA homodimer is essential for future studies on the functional dynamics of this membrane transporter.     

Characterization of the LSGGQ and H motifs from the Escherichia coli lipid A transporter MsbA

ATP-binding cassette (ABC) transporters make up one of the largest superfamilies of proteins known and have been shown to transport substrates ranging from lipids and antibiotics to sugars and amino acids. The dysfunction of ABC transporters has been linked to human pathologies such as cystic fibrosis, hyperinsulinemia, and macular dystrophy. Several bacterial ABC transporters are also necessary for bacterial survival and transport of virulence factors in an infected host. MsbA is a 65 kDa protein that forms a functional homodimer consisting of two six-helix transmembrane domains and two approximately 250 amino acid nucleotide-binding domains (NBD). The NBDs contain several conserved regions such as the Walker A, LSGGQ, and H motif that bind directly to ATP and align it for hydrolysis. MsbA transports lipid A, its native substrate, across the inner membrane of Gram-negative bacteria. The loss or dysfunction of MsbA results in a toxic accumulation of lipid A inside the cell, leading to cell-membrane instability and cell death. Using site-directed spin labeling electron paramagnetic resonance spectroscopy, conserved motifs within the MsbA NBD have been evaluated for structure and dynamics upon substrate binding. It has been determined that the LSGGQ NBD consensus sequence is consistent with an alpha-helical conformation and that these residues maintain extensive tertiary contacts throughout hydrolysis. The dynamics of the LSGGQ and the H-motif region have been studied in the presence of ATP, ADP, and ATP plus vanadate to identify the residues that are directly affected by interactions with the substrate before, after, and during hydrolysis, respectively.     

ATP hydrolysis at one of the two sites in ABC transporters initiates transport related conformational transitions

ABC transporters play important roles in all types of organisms by participating in physiological and pathological processes. In order to modulate the function of ABC transporters, detailed knowledge regarding their structure and dynamics is necessary. Available structures of ABC proteins indicate three major conformations, a nucleotide-bound "bottom-closed" state with the two nucleotide binding domains (NBDs) tightly closed, and two nucleotide-free conformations, the "bottom-closed" and the "bottom-open", which differ in the extent of separation of the NBDs. However, it remains a question how the widely open conformation should be interpreted, and whether hydrolysis at one of the sites can drive conformational transitions while the NBDs remain in contact. To extend our knowledge, we have investigated the dynamic properties of the Sav1866 transporter using molecular dynamics (MD) simulations. We demonstrate that the replacement of one ATP by ADP alters the correlated motion patterns of the NBDs and the transmembrane domains (TMD). The results suggest that the hydrolysis of a single nucleotide could lead to extracellular closure, driving the transport cycle. Essential dynamics analysis of simulations suggests that single nucleotide hydrolysis can drive the system toward a "bottom-closed" apo conformation similar to that observed in the structure of the MsbA transporter. We also found significant structural instability of the "bottom-open" form of the transporters in simulations. Our results suggest that ATP hydrolysis at one of the sites promotes transport related conformational changes leading to the "bottom-closed" apo conformation, which could thus be physiologically more relevant for describing the structure of the apo state.     

The structures of MsbA: Insight into ABC transporter-mediated multidrug efflux

ATP-binding cassette (ABC) transporters are integral membrane proteins that couple ATP hydrolysis to the transport of various molecules across cellular membranes. Found in both prokaryotes and eukaryotes, a sub-group of these transporters are involved in the efflux of hydrophobic drugs and lipids, causing anti-microbial and chemotherapeutic multidrug resistance. In this review, we examine recent structural and functional analysis of the ABC transporter MsbA and implications on the mechanism of multidrug efflux.     

Time-resolved Fourier transform infrared spectroscopy of the nucleotide-binding domain from the ATP-binding Cassette transporter MsbA: ATP hydrolysis is the rate-limiting step in the catalytic cycle

MsbA is an essential Escherichia coli ATP-binding cassette (ABC) transporter involved in the flipping of lipid A across the cytoplasmic membrane. It is a close homologue of human P-glycoprotein involved in multidrug resistance, and it similarly accepts a variety of small hydrophobic xenobiotics as transport substrates. X-ray structures of three full-length ABC multidrug exporters (including MsbA) have been published recently and reveal large conformational changes during the transport cycle. However, how ATP hydrolysis couples to these conformational changes and finally the transport is still an open question. We employed time-resolved FTIR spectroscopy, a powerful method to elucidate molecular reaction mechanisms of soluble and membrane proteins, to address this question with high spatiotemporal resolution. Here, we monitored the hydrolysis reaction in the nucleotide-binding domain of MsbA at the atomic level. The isolated MsbA nucleotide-binding domain hydrolyzed ATP with V(max) = 45 nmol mg(-1) min(-1), similar to the full-length transporter. A Hill coefficient of 1.49 demonstrates positive cooperativity between the two catalytic sites formed upon dimerization. Global fit analysis of time-resolved FTIR data revealed two apparent rate constants of ~1 and 0.01 s(-1), which were assigned to formation of the catalytic site and hydrolysis, respectively. Using isotopically labeled ATP, we identified specific marker bands for protein-bound ATP (1245 cm(-1)), ADP (1101 and 1205 cm(-1)), and free phosphate (1078 cm(-1)). Cleavage of the β-phosphate-γ-phosphate bond was found to be the rate-limiting step; no protein-bound phosphate intermediate was resolved.     

The ATP binding cassette multidrug transporter LmrA and lipid transporter MsbA have overlapping substrate specificities

LmrA is an ATP binding cassette (ABC) multidrug transporter in Lactococcus lactis that is a structural and functional homologue of the human multidrug resistance P-glycoprotein MDR1 (ABCB1). LmrA is also homologous to MsbA, an essential ABC transporter in Escherichia coli involved in the trafficking of lipids, including Lipid A. We have compared the substrate specificities of LmrA and MsbA in detail. Surprisingly, LmrA was able to functionally substitute for a temperature-sensitive mutant MsbA in E. coli WD2 at non-permissive temperatures, suggesting that LmrA could transport Lipid A. LmrA also exhibited a Lipid A-stimulated, vanadate-sensitive ATPase activity. Reciprocally, the expression of MsbA conferred multidrug resistance on E. coli. Similar to LmrA, MsbA interacted with photoactivatable substrate [3H]azidopine, displayed a daunomycin, vinblastine, and Hoechst 33342-stimulated vanadate-sensitive ATPase activity, and mediated the transport of ethidium from cells and Hoechst 33342 in proteoliposomes containing purified and functionally reconstituted protein. Taken together, these data demonstrate that MsbA and LmrA have overlapping substrate specificities. Our observations imply the presence of structural elements in the recently published crystal structures of MsbA in E. coli and Vibrio cholera (Chang, G., and Roth, C. B. (2001) Science 293, 1793-1800; Chang, G. (2003) J. Mol. Biol. 330, 419-430) that support drug-protein interactions and suggest a possible role for LmrA in lipid trafficking in L. lactis.     

Characterization of the E506Q and H537A dysfunctional mutants in the E. coli ABC transporter MsbA

MsbA is a member of the ABC transporter superfamily that is specifically found in Gram-negative bacteria and is homologous to proteins involved in both bacterial and human drug resistance. The E506Q and H537A mutations have been introduced and used for crystallization of other members of the ABC transporter protein family, including BmrA and the ATPase domains MalK, HlyB-NBD, and MJ0796, but have not been previously studied in detail or investigated in the MsbA lipid A exporter. We utilized an array of biochemical and EPR spectroscopy approaches to characterize the local and global effects of these nucleotide binding domain mutations on the E. coli MsbA homodimer. The lack of cell viability in an in vivo growth assay confirms that the presence of the E506Q or H537A mutations within MsbA creates a dysfunctional protein. To further investigate the mode of dysfunction, a fluorescent ATP binding assay was used and showed that both mutant proteins maintain their ability to bind ATP, but ATPase assays indicate hydrolysis is severely inhibited by each mutation. EPR spectroscopy data using previously identified and characterized reporter sites within the nucleotide binding domain along with ATP detection assays show that hydrolysis does occur over time in both mutants, though more readily in the H537A protein. DEER spectroscopy demonstrates that both proteins studied are purified in a closed dimer conformation, indicating that events within the cell can induce a stable, closed conformation of the MsbA homodimer that does not reopen even in the absence of nucleotide.     

An Escherichia coli mutant defective in lipid export

Escherichia coli phospholipids and lipopolysaccharide, made on the inner surface of the inner membrane, are rapidly transported to the outer membrane by mechanisms that are not well characterized. We now report a temperature-sensitive mutant (WD2) with an A270T substitution in a trans-membrane region of the ABC transporter MsbA. As shown by (32)P(i) and (14)C-acetate labeling, export of all major lipids to the outer membrane is inhibited by approximately 90% in WD2 after 30 min at 44 degrees C. Transport of newly synthesized proteins is not impaired. Electron microscopy shows reduplicated inner membranes in WD2 at 44 degrees C, consistent with a key role for MsbA in lipid trafficking.     

The flexibility of P-glycoprotein for its poly-specific drug binding from molecular dynamics simulations

The multidrug efflux pump P-glycoprotein (P-gp) contributes to multidrug resistance in about half of human cancers. Recently, high resolution X-ray crystal structures of mouse P-gp (inward-facing) were reported, which significantly facilitates the understanding of the function of P-gp and the structure-based design of inhibitors for P-gp. Here we perform 20?ns molecular dynamics simulations of inward-facing P-gp with/without ligand in explicit lipid and water to investigate the flexibility of P-gp for its poly-specific drug binding. By analyzing the interactions between P-gp and QZ59-RRR or QZ59-SSS, we summarize the important residues and the flexibility of different parts of P-gp. Particularly, the flexibility of the side chains of aromatic residues (Phe and Tyr) allows them to form rotamers with different orientations in the binding pocket, which plays a critical role for the poly-specificity of the drug-binding cavity of P-gp. MD simulations reveal that trans-membrane (TM) TM12 and TM6 are flexible and contribute to the poly-specific drug binding, while TM4 and TM5 are rigid and stabilize the whole structure. We also construct outward-facing P-gp based on the MsbA structure and perform 20?ns MD simulations. The comparison between the MD results for outward-facing P-gp and those for inward-facing P-gp shows that the TM parts in outward-facing P-gp undergo significant conformational change to facilitate the export of small molecules.     

A structural model for the open conformation of the mdr1 P-glycoprotein based on the MsbA crystal structure

The validity of the structure of the Escherichia coli MsbA lipid transporter as a model from the mdr1 P-glycoprotein has been evaluated. Comparative sequence analyses, motif search and secondary structure prediction indicated that each of the two P-glycoprotein halves is structurally similar to the MsbA monomer and also suggested that the open dimer structure is valid for P-glycoprotein. Homology modeling was used to predict the structure of P-glycoprotein using MsbA as a template. The resulting modeled structure allowed a detailed study of the interactions between the intracellular domain and the nucleotide binding domain and suggested that these contacts are involved in mediating the coupling between nucleotide binding domain conformational changes and transmembrane helices reorientation during transport. In P-glycoprotein, the internal chamber open to the inner leaflet and the inner medium is significantly different in size and charge than in MsbA. These differences can be related to those of the transported substrates. Moreover an ensemble of 20 conserved aromatic residues appears to border the periphery of each side of the chamber in P-glycoprotein. These may be important for size selection and proper positioning of drugs for transport. The relevance of the modeled conformation to P-gp function is discussed.     

Substrate Binding Stabilizes a Pre-translocation Intermediate in the ATP-binding Cassette Transport Protein MsbA

ATP-binding cassette (ABC) transporters belong to one of the largest protein superfamilies that expands from prokaryotes to man. Recent x-ray crystal structures of bacterial and mammalian ABC exporters suggest a common alternating access mechanism of substrate transport, which has also been biochemically substantiated. However, the current model does not yet explain the coupling between substrate binding and ATP hydrolysis that underlies ATP-dependent substrate transport. In our studies on the homodimeric multidrug/lipid A ABC exporter MsbA from Escherichia coli, we performed cysteine cross-linking, fluorescence energy transfer, and cysteine accessibility studies on two reporter positions, near the nucleotide-binding domains and in the membrane domains, for transporter embedded in a biological membrane. Our results suggest for the first time that substrate binding by MsbA stimulates the maximum rate of ATP hydrolysis by facilitating the dimerization of nucleotide-binding domains in a state, which is markedly distinct from the previously described nucleotide-free, inward-facing and nucleotide-bound, outward-facing conformations of ABC exporters and which binds ATP.     

Structure of MsbA from Vibrio cholera: a multidrug resistance ABC transporter homolog in a closed conformation

The spread of multidrug resistance (MDR) is a world health crisis that presents a significant challenge to the treatment of cancer and infection. MDR can be caused by a group of ABC (MDR-ABC) transporters that move hydrophobic drug molecules and lipids across the cell membrane. To gain insight into the conformational changes these transporters undergo when flipping hydrophobic substrates across the lipid bilayer, we have determined the structure of the lipid flippase MsbA from Vibrio cholera (VC-MsbA) to 3.8A. Structural comparison of VC-MsbA to MsbA from Escherichia coli reveals that the transporters share a structurally conserved core of transmembrane alpha-helices, but differ in the relative orientations of their nucleotide-binding domains (NBD). The transmembrane domain of VC-MsbA is captured in a closed conformation and the structure supports a "power stroke" model of transporter dynamics where opposing NBDs associate upon ATP binding. The separation of the alpha and beta domains of the NBD suggests the possibility that their association could make them competent to bind ATP and gives further insight into the structural basis for catalytic regulation.     

Effects of the L511P and D512G mutations on the Escherichia coli ABC transporter MsbA

MsbA is a member of the ABC transporter superfamily and is homologous to ABC transporters linked to multidrug resistance. The nucleotide binding domains (NBDs) of these proteins include conserved motifs that are involved in ATP binding, including conserved SALD residues (D-loop) that are diagnostic in identifying ABC transporters but whose roles have not been identified. Within the D-loop, single point mutations L511P and D512G were discovered by random mutational analysis of MsbA to disrupt protein function in the cell [Polissi, A., and Georgopoulos, C. (1996) Mol. Microbiol. 20, 1221-1233] but have not been further studied in MsbA or in detail in any other ABC transporter. In these studies, we show that both L511P and D512G mutants of MsbA are able to bind ATP at near-wild-type levels but are unable to maintain cell viability in an in vivo growth assay, verifying the theory that they are dysfunctional at some point after ATP binding. An ATPase assay further suggests that the L511P mutation prevents effective ATP hydrolysis, and an ATP detection assay reveals that only small amounts of ATP are hydrolyzed; D512G is able to hydrolyze ATP at a rate 3-fold faster than that of the wild type. EPR spectroscopy studies using reporter sites within the NBDs also indicate that at least some hydrolysis occurs in L511P or D512G MsbA but show fewer spectral changes than observed for the same reporters in the wild-type background. These studies indicate that L511 is necessary for efficient ATP hydrolysis and D512 is essential for conformational rearrangements required for flipping lipid A.     

Molecular dynamics simulation and molecular docking studies of 1,4-Dihydropyridines as P-glycoprotein’s allosteric inhibitors

P-glycoprotein (P-gp) is a main factor contributing to multidrug resistance. The effect of this transporter protein on limiting the effectiveness of chemotherapy has been shown by various studies. In a previous report, we synthesized some 14-dihydropyridine (DHP) derivatives as inhibitors of human P-gp. In the present study, a computational approach has been exploited to reveal the main interactions between DHPs and P-gp. In order to do this, homology modeling was performed to obtain a model of the protein. Then, molecular dynamics simulation was used to refine the constructed model of P-gp in the presence of the lipids bilayer. Model validation was performed with several tools. Finally, molecular docking followed by MD simulation of ligand–protein complex was employed to elucidate the binding mode and the dynamical changes of protein with/without DHPs bound. The results emphasized that interaction of the residues Gln912, Ser909, Arg905, Ser474, Val472 with DHPs play a crucial role in the inhibitory of these ligands and this was in a relatively good accordance with the results reported in the experimental studies.     

PfMDR1: mechanisms of transport modulation by functional polymorphisms

ATP-Binding Cassette (ABC) transporters are efflux pumps frequently associated with multidrug resistance in many biological systems, including malaria. Antimalarial drug-resistance involves an ABC transporter, PfMDR1, a homologue of P-glycoprotein in humans. Twenty years of research have shown that several single nucleotide polymorphisms in pfmdr1 modulate in vivo and/or in vitro drug susceptibility. The underlying physiological mechanism of the effect of these mutations remains unclear. Here we develop structural models for PfMDR1 in different predicted conformations, enabling the study of transporter motion. Such analysis of functional polymorphisms allows determination of their potential role in transport and resistance. The bacterial MsbA ABC pump is a PfMDR1 homologue. MsbA crystals in different conformations were used to create PfMDR1 models with Modeller software. Sequences were aligned with ClustalW and analysed by Ali2D revealing a high level of secondary structure conservation. To validate a potential drug binding pocket we performed antimalarial docking simulations. Using aminoquinoline as probe drugs in PfMDR1 mutated parasites we evaluated the physiology underlying the mechanisms of resistance mediated by PfMDR1 polymorphisms. We focused on the analysis of well known functional polymorphisms in PfMDR1 amino acid residues 86, 184, 1034, 1042 and 1246. Our structural analysis suggested the existence of two different biophysical mechanisms of PfMDR1 drug resistance modulation. Polymorphisms in residues 86/184/1246 act by internal allosteric modulation and residues 1034 and 1042 interact directly in a drug pocket. Parasites containing mutated PfMDR1 variants had a significant altered aminoquinoline susceptibility that appears to be dependent on the aminoquinoline lipophobicity characteristics as well as vacuolar efflux by PfCRT. We previously described the in vivo selection of PfMDR1 polymorphisms under antimalarial drug pressure. Now, together with recent PfMDR1 functional reports, we contribute to the understanding of the specific structural role of these polymorphisms in parasite antimalarial drug response.