Identification of Possible Binding Sites for Morphine and Nicardipine on the Multidrug Transporter P-Glycoprotein Using Umbrella Sampling Techniques

Membrane-assisted tariquidar access and binding mechanisms of human ATP-binding cassette transporter P-glycoprotein

The human multidrug transporter P-glycoprotein (P-gp) is physiologically essential and of key relevance to biomedicine. Recent structural studies have shed light on the mode of inhibition of the third-generation inhibitors for human P-gp, but the molecular mechanism by which these inhibitors enter the transmembrane sites remains poorly understood. In this study, we utilized all-atom molecular dynamics (MD) simulations to characterize human P-gp dynamics under a potent inhibitor, tariquidar, bound condition, as well as the atomic-level binding pathways in an explicit membrane/water environment. Extensive unbiased simulations show that human P-gp remains relatively stable in tariquidar-free and bound states, while exhibiting a high dynamic binding mode at either the drug-binding pocket or the regulatory site. Free energy estimations by partial nudged elastic band (PNEB) simulations and Molecular Mechanics Generalized Born Surface Area (MM/GBSA) method identify two energetically favorable binding pathways originating from the cytoplasmic gate with an extended tariquidar conformation. Interestingly, free tariquidar in the lipid membrane predominantly adopts extended conformations similar to those observed at the regulatory site. These results suggest that membrane lipids may preconfigure tariquidar into an active ligand conformation for efficient binding to the regulatory site. However, due to its conformational plasticity, tariquidar ultimately moves toward the drug-binding pocket in both pathways, explaining how it acts as a substrate at low concentrations. Our molecular findings propose a membrane-assisted mechanism for the access and binding of the third-generation inhibitors to the binding sites of human P-gp, and offer deeper insights into the molecule design of more potent inhibitors against P-gp-mediated drug resistance.

Interaction of a Homologous Series of Amphiphiles with P-glycoprotein in a Membrane Environment-Contributions of Polar and Non-Polar Interactions

The transport of drugs by efflux transporters in biomembranes limits their bioavailability and is a major determinant of drug resistance development by cancer cells and pathogens. A large number of chemically dissimilar drugs are transported, and despite extensive studies, the molecular determinants of substrate specificity are still not well understood. In this work, we explore the role of polar and non-polar interactions on the interaction of a homologous series of fluorescent amphiphiles with the efflux transporter P-glycoprotein. The interaction of the amphiphiles with P-glycoprotein is evaluated through effects on ATPase activity, efficiency in inhibition of [¹²⁵I]-IAAP binding, and partition to the whole native membranes containing the transporter. The results were complemented with partition to model membranes with a representative lipid composition, and details on the interactions established were obtained from MD simulations. We show that when the total concentration of amphiphile is considered, the binding parameters obtained are apparent and do not reflect the affinity for P-gp. A new formalism is proposed that includes sequestration of the amphiphiles in the lipid bilayer and the possible binding of several molecules in P-gp's substrate-binding pocket. The intrinsic binding affinity thus obtained is essentially independent of amphiphile hydrophobicity, highlighting the importance of polar interactions. An increase in the lipophilicity and amphiphilicity led to a more efficient association with the lipid bilayer, which maintains the non-polar groups of the amphiphiles in the bilayer, while the polar groups interact with P-gp's binding pocket. The presence of several amphiphiles in this orientation is proposed as a mechanism for inhibition of P-pg function.

Vinca alkaloid binding to P-glycoprotein occurs in a processive manner

A mechanistic understanding of how P-glycoprotein (Pgp) is able to bind and transport its astonishing range of substrates remains elusive. Pharmacological data demonstrated the presence of at least four distinct binding sites, but their locations have not been fully elucidated. The combination of biochemical and structural data suggests that initial binding may occur in the central cavity or at the lipid-protein interface. Our objective was to define the binding sites for two transported substrates of Pgp; the anticancer drug vinblastine and the fluorescent probe rhodamine 123. A series of mutations was generated in positions proximal to previously defined drug-interacting residues on Pgp. The protein was purified and reconstituted into styrene-maleic acid lipid particles (SMALPs) to measure the apparent drug binding constant or into liposomes for assessment of drug-stimulated ATP hydrolysis. The biochemical data were reconciled with structural models of Pgp using molecular docking. The data indicated that the binding of rhodamine 123 occurred predominantly within the central cavity of Pgp. In contrast, the significantly more hydrophobic vinblastine bound to both the lipid-protein interface and within the central cavity. The data suggest that the initial interaction of vinca alkaloids with Pgp occurs at the lipid interface followed by internalisation into the central cavity, which also provides the transport conduit. This model is supported by recent structural observations with Pgp and early biophysical and cross-linking approaches. Moreover, the proposed model illustrates that the broad substrate profile for Pgp is underpinned by a combination of multiple initial interaction sites and an accommodating transport conduit.

Efflux dynamics of the antiseizure drug, levetiracetam, through the P-glycoprotein channel revealed by advanced comparative molecular simulations

Understanding the precise mechanistic details of the possible binding and transport of antiseizure medications (ASMs) through the P-glycoprotein (P-gp) efflux pump is essential to find strategies for the treatment of patients with epilepsy resistant to ASMs. In the present work, conventional molecular dynamics, binding free energy calculations, steered molecular dynamics and umbrella sampling were applied to study the interactions of levetiracetam and brivaracetam with P-gp and their possible egress path from the binding site. Comparative results for the control drugs, zosuquidar and verapamil, confirmed their established P-gp inhibitory activity. Brivaracetam, a non-substrate of P-gp, demonstrated stronger static and dynamic interactions with the exporter protein, than levetiracetam. The potential of mean force calculations indicated that the energy barriers through the ligand export were the lowest for levetiracetam, suggesting the drug as a P-gp substrate with facile passage through the transporter channel. Our findings also stressed the contribution of nonpolar interactions with P-gp channel lining as well as with membrane lipid molecules to hamper the ASM efflux by the transmembrane exporter. Appropriate structural engineering of the ASMs is thus recommended to address drug-resistant epilepsy.

Investigation of the molecular and mechanistic basis for the regioselective metabolism of midazolam by cytochrome P450 3A4

Cytochrome P450 3A4 (CYP3A4) is the most important P450 enzyme for drug metabolism and drug-drug interaction, due to it being responsible for the biotransformation of approximately 50% of clinically used drugs. Advance knowledge of the molecular and mechanistic basis of CYP3A4 regioselective metabolism is beneficial for understanding the production of metabolites, and may allow personalized metabolic pathways or designing pathway-specific therapeutics. In this work, we focus on investigating the ligand-receptor interactions, substrate conformational transition, and key factors regulating the specificity of metabolic pathways using midazolam (MDZ) as a probe. Here, three types of substrate-binding conformations related to the diversity of MDZ metabolites are identified. The results also suggest that an allosteric site for MDZ is located near the F'-helix, A-anchor, and C-terminal loop of CYP3A4. The presence of an effector in the allosteric site can accelerate the conformational transition of the substrate via modulating a "sandwich" structure, and may affect the proportion of metabolites at high substrate concentration. We hope that the results can improve the understanding of the CYP3A4 structure and function, and provide a new perspective for drug development.

Extended-ensemble docking to probe dynamic variation of ligand binding sites during large-scale structural changes of proteins

Proteins can sample a broad landscape as they undergo conformational transition between different functional states. At the same time, as key players in almost all cellular processes, proteins are important drug targets. Considering the different conformational states of a protein is therefore central for a successful drug-design strategy. Here we introduce a novel docking protocol, termed extended-ensemble docking, pertaining to proteins that undergo large-scale (global) conformational changes during their function. In its application to multidrug ABC-transporter P-glycoprotein (Pgp), extensive non-equilibrium molecular dynamics simulations employing system-specific collective variables are first used to describe the transition cycle of the transporter. An extended set of conformations (extended ensemble) representing the full transition cycle between the inward- and the outward-facing states is then used to seed high-throughput docking calculations of known substrates, non-substrates, and modulators of the transporter. Large differences are predicted in the binding affinities to different conformations, with compounds showing stronger binding affinities to intermediate conformations compared to the starting crystal structure. Hierarchical clustering of the binding modes shows all ligands preferably bind to the large central cavity of the protein, formed at the apex of the transmembrane domain (TMD), whereas only small binding populations are observed in the previously described R and H sites present within the individual TMD leaflets. Based on the results, the central cavity is further divided into two major subsites, first preferably binding smaller substrates and high-affinity inhibitors, whereas the second one shows preference for larger substrates and low-affinity modulators. These central subsites along with the low-affinity interaction sites present within the individual TMD leaflets may respectively correspond to the proposed high- and low-affinity binding sites in Pgp. We propose further an optimization strategy for developing more potent inhibitors of Pgp, based on increasing its specificity to the extended ensemble of the protein, instead of using a single protein structure, as well as its selectivity for the high-affinity binding site. In contrast to earlier in silico studies using single static structures of Pgp, our results show better agreement with experimental studies, pointing to the importance of incorporating the global conformational flexibility of proteins in future drug-discovery endeavors.

Effect of the Force Field on Molecular Dynamics Simulations of the Multidrug Efflux Protein P-Glycoprotein

Molecular dynamics (MD) simulations have been used extensively to study P-glycoprotein (P-gp), a flexible multidrug transporter that is a key player in the development of multidrug resistance to chemotherapeutics. A substantial body of literature has grown from simulation studies that have employed various simulation conditions and parameters, including AMBER, CHARMM, OPLS, GROMOS, and coarse-grained force fields, drawing conclusions from simulations spanning hundreds of nanoseconds. Each force field is typically parametrized and validated on different data and observables, usually of small molecules and peptides; there have been few comparisons of force field performance on large protein-membrane systems. Here we compare the conformational ensembles of P-gp embedded in a POPC/cholesterol bilayer generated over 500 ns of replicate simulation with five force fields from popular biomolecular families: AMBER 99SB-ILDN, CHARMM 36, OPLS-AA/L, GROMOS 54A7, and MARTINI. We find considerable differences among the ensembles with little conformational overlap, although they correspond to similar extents to structural data obtained from electron paramagnetic resonance and cross-linking studies. Moreover, each trajectory was still sampling new conformations at a high rate after 500 ns of simulation, suggesting the need for more sampling. This work highlights the need to consider known limitations of the force field used (e.g., biases toward certain secondary structures) and the simulation itself (e.g., whether sufficient sampling has been achieved) when interpreting accumulated results of simulation studies of P-gp and other transport proteins.

Development of the computational antibiotic screening platform (CLASP) to aid in the discovery of new antibiotics

Bacterial colonization of biotic and abiotic surfaces and antibiotic resistance are grand challenges with paramount societal impacts. However, in the face of increasing bacterial resistance to all known antibiotics, efforts to discover new classes of antibiotics have languished, creating an urgent need to accelerate the antibiotic discovery pipeline. A major deterrent in the discovering of new antibiotics is the limited permeability of molecules across the bacterial envelope. Notably, the Gram-negative bacteria have nutrient specific protein channels (or porins) that restrict the permeability of non-essential molecules, including antibiotics. Here, we have developed the Computational Antibiotic Screening Platform (CLASP) for screening of potential drug molecules through the porins. The CLASP takes advantage of coarse grain (CG) resolution, advanced sampling techniques, and a parallel computing environment to maximize its performance. The CLASP yields comprehensive thermodynamic and kinetic output data of a potential drug molecule within a few hours of wall-clock time. Its output includes the potential of mean force profile, energy barrier, the rate constant, and contact analysis of the molecule with the pore-lining residues, and the orientational analysis of the molecule in the porin channel. In our first CLASP application, we report the transport properties of six carbapenem antibiotics-biapenem, doripenem, ertapenem, imipenem, meropenem, and panipenem-through OccD3, a major channel for carbapenem uptake in Pseudomonas aeruginosa. The CLASP is designed to screen small molecule libraries with a fast turnaround time to yield structure-property relationships to discover antibiotics with high permeability. The CLASP will be freely distributed to enable accelerated antibiotic drug discovery.

Efflux mechanism and pathway of verapamil pumping by human P-glycoprotein

Multidrug resistance (MDR) caused by overexpressed permeability-glycoprotein (P-gp) in cancer cells is the main barrier for the cure of cancers. P-gp can pump many chemotherapeutic drugs, which is a viable target to overcome P-gp-mediated MDR by efficient inhibitors of P-gp. However, limited understanding of the efflux mechanism by human P-gp hinders the development of efficient inhibitors. Herein, the transport of a P-gp inhibitor, verapamil, by human P-gp has been investigated using targeted molecular dynamics simulations and energetics analysis based on our previous research on the transport of a drug (doxorubicin). The energetics analysis identifies that the driving forces for the transport of verapamil are electrostatic repulsions contributed by the positively charged residues in the initial stage and then hydrophobic interactions contributed by the important residues in the later stage. This scenario is generally consistent with that in the transport of doxorubicin. However, the positively charged residues and the important residues for the transport of verapamil are incompletely consistent with the relative residues for the transport of doxorubicin. Moreover, the binding free energy contributions of the positively charged residues for the transport of verapamil are generally higher than them for the transport of doxorubicin, while the important residues constitute significantly different binding free energy compositions in the transports of the two substrates. Consequently, the pathway for the transport of verapamil is identified, which shares only two residues (F336 and M986) with the pathway of doxorubicin. This may imply the weak competitiveness of verapamil with doxorubicin in the substrate efflux. Taken together, this work provided new insights into the efflux mechanisms by human P-gp and would be beneficial in the design of potent P-gp inhibitors.

On interactions of P-glycoprotein with various anti-tumor drugs by binding free energy calculations

P-glycoprotein (P-gp, MDR1), one of ATP-binding cassette transporters, may confer tumor cells cross-resistance to chemotherapeutics. A large amount of P-gp inhibitors were designed to inhibit the multidrug resistance (MDR) feature of P-gp. However, no sufficient researches were reported to explore the correlation between binding capacity and drug property by experiment. Without particular drug property found to inhibit the MDR feature of P-gp, the orientation of drug design is indefinite. In this work, 10 representative cancer drugs with various properties are used to bind with P-gp by molecular dynamics simulation. Binding free energy between P-gp and 10 drugs ranges -139 to -253 kJ/mol. It reveals that the promiscuity nature of P-gp is in light of the similar binding free energy in separate P-gp-ligand binding systems. The binding effect of P-gp and drugs correlates well with the size of drugs and has no apparent correlation with the polarity of each drug. The key reason is that van der Waal's interaction occupies most of the total binding free energy, and it is led by the number of atoms in the drugs. Two transmembrane segments (TM6 and TM12) and three types of amino acids (PHE, MET, and GLN) are vital in binding drugs with van der Waal's energy, which evident the influence between binding stability and size of drugs. This work provides the cause and theoretical basis for the promiscuity nature of P-gp.Communicated by Ramaswamy H. Sarma.

Endogenous Opiates and Behavior: 2016

This paper is the thirty-ninth consecutive installment of the annual review of research concerning the endogenous opioid system. It summarizes papers published during 2016 that studied the behavioral effects of molecular, pharmacological and genetic manipulation of opioid peptides, opioid receptors, opioid agonists and opioid antagonists. The particular topics that continue to be covered include the molecular-biochemical effects and neurochemical localization studies of endogenous opioids and their receptors related to behavior, and the roles of these opioid peptides and receptors in pain and analgesia, stress and social status, tolerance and dependence, learning and memory, eating and drinking, drug abuse and alcohol, sexual activity and hormones, pregnancy, development and endocrinology, mental illness and mood, seizures and neurologic disorders, electrical-related activity and neurophysiology, general activity and locomotion, gastrointestinal, renal and hepatic functions, cardiovascular responses, respiration and thermoregulation, and immunological responses.

Coarse-grained molecular dynamics simulations reveal lipid access pathways in P-glycoprotein

P-glycoprotein (P-gp) exports a broad range of dissimilar compounds, including drugs, lipids, and lipid-like molecules. Because of its substrate promiscuity, P-gp is a key player in the development of cancer multidrug resistance. Although P-gp is one of the most studied ABC transporters, the mechanism by which its substrates access the cavity remains unclear. In this study, we perform coarse-grained molecular dynamics simulations to explore possible lipid access pathways in the inward-facing conformation of P-gp embedded in bilayers of different lipid compositions. In the inward-facing orientation, only lipids from the lower leaflet access the cavity of the transporter. We identify positively charged residues at the portals of P-gp that favor lipid entrance to the cavity, as well as lipid-binding sites at the portals and within the cavity, which is in good agreement with previous experimental studies. This work includes several examples of lipid pathways for phosphatidylcholine and phosphatidylethanolamine lipids that help elucidate the molecular mechanism of lipid binding in P-gp.

The reliability of molecular dynamics simulations of the multidrug transporter P-glycoprotein in a membrane environment

Despite decades of research, the mechanism of action of the ABC multidrug transporter P-glycoprotein (P-gp) remains elusive. Due to experimental limitations, many researchers have turned to molecular dynamics simulation studies in order to investigate different aspects of P-gp function. However, such studies are challenging and caution is required when interpreting the results. P-gp is highly flexible and the time scale on which it can be simulated is limited. There is also uncertainty regarding the accuracy of the various crystal structures available, let alone the structure of the protein in a physiologically relevant environment. In this study, three alternative structural models of mouse P-gp (3G5U, 4KSB, 4M1M), all resolved to 3.8 Å, were used to initiate sets of simulations of P-gp in a membrane environment in order to determine: a) the sensitivity of the results to differences in the starting configuration; and b) the extent to which converged results could be expected on the times scales commonly simulated for this system. The simulations suggest that the arrangement of the nucleotide binding domains (NBDs) observed in the crystal structures is not stable in a membrane environment. In all simulations, the NBDs rapidly associated (within 10 ns) and changes within the transmembrane helices were observed. The secondary structure within the transmembrane domain was best preserved in the 4M1M model under the simulation conditions used. However, the extent to which replicate simulations diverged on a 100 to 200 ns timescale meant that it was not possible to draw definitive conclusions as to which structure overall was most stable, or to obtain converged and reliable results for any of the properties examined. The work brings into question the reliability of conclusions made in regard to the nature of specific interactions inferred from previous simulation studies on this system involving similar sampling times. It also highlights the need to demonstrate the statistical significance of any results obtained in simulations of large flexible proteins, especially where the initial structure is uncertain.

Comparison of mechanistic transport cycle models of ABC exporters

ABC (ATP binding cassette) transporters, ubiquitous in all kingdoms of life, carry out essential substrate transport reactions across cell membranes. Their transmembrane domains bind and translocate substrates and are connected to a pair of nucleotide binding domains, which bind and hydrolyze ATP to energize import or export of substrates. Over four decades of investigations into ABC transporters have revealed numerous details from atomic-level structural insights to their functional and physiological roles. Despite all these advances, a comprehensive understanding of the mechanistic principles of ABC transporter function remains elusive. The human multidrug resistance transporter ABCB1, also referred to as P-glycoprotein (P-gp), is one of the most intensively studied ABC exporters. Using ABCB1 as the reference point, we aim to compare the dominating mechanistic models of substrate transport and ATP hydrolysis for ABC exporters and to highlight the experimental and computational evidence in their support. In particular, we point out in silico studies that enhance and complement available biochemical data. “This article is part of a Special Issue entitled: Beyond the Structure Function Horizon of Membrane Proteins edited by Ute Hellmich, Rupak Doshi and Benjamin McIlwain.”

Homology modelling of human P-glycoprotein

P-glycoprotein (P-gp) is an ATP-binding cassette transporter that exports a huge range of compounds out of cells and is thus one of the key proteins in conferring multi-drug resistance in cancer. Understanding how it achieves such a broad specificity and the series of conformational changes that allow export to occur form major, on-going, research objectives around the world. Much of our knowledge to date has been derived from mutagenesis and assay data. However, in recent years, there has also been great progress in structural biology and although the structure of human P-gp has not yet been solved, there are now a handful of related structures on which homology models can be built to aid in the interpretation of the vast amount of experimental data that currently exists. Many models for P-gp have been built with this aim, but the situation is complicated by the apparent flexibility of the system and by the fact that although many potential templates exist, there is large variation in the conformational state in which they have been crystallized. In this review, we summarize how homology modelling has been used in the past, how models are typically selected and finally illustrate how MD simulations can be used as a means to give more confidence about models that have been generated via this approach.

Structural and dynamic perspectives on the promiscuous transport activity of P-glycoprotein

The multidrug transporter P-glycoprotein (P-gp) is expressed in the blood-brain barrier endothelium where it effluxes a range of drug substrates, preventing their accumulation within the brain. P-gp has been studied extensively for 40 years because of its crucial role in the absorption, distribution, metabolism and elimination of a range of pharmaceutical compounds. Despite this, many aspects of the structure-function mechanism of P-gp are unresolved. Here we review the emerging role of molecular dynamics simulation techniques in our understanding of the membrane-embedded conformation of P-gp. We discuss its conformational plasticity in the presence and absence of ATP, and recent efforts to characterize the drug binding sites and uptake pathways.

Do adsorbed drugs onto P-glycoprotein influence its efflux capability?

The membrane biophysical aspects by which multidrug resistance (MDR) relate to the ABC transporter function still remain largely unknown. Notwithstanding the central role that efflux pumps like P-glycoprotein have in MDR onset, experimental studies classified additionally the lipid micro-environment where P-gp is inserted as a determinant for the increased efflux capability demonstrated in MDR cell lines. Recently, a nonlinear model for drug-membrane interactions showed that, upon drug adsorption, long-range mechanical alterations are predicted to affect the P-gp ATPase function at external drug concentrations of ∼10-100 μM. However, our results also show that drug adsorption may also occur at P-gp nucleotide-binding domains where conformational changes drive the efflux cycle. Thus, we assessed the effect of drug adsorption to both protein-water and lipid-water interfaces by means of molecular dynamics simulations. The results show that free energies of adsorption are lower for modulators in both lipid/water and protein/water interfaces. Important differences in drug-protein interactions, protein dynamics and membrane biophysical characteristics were observed between the different classes. Therefore, we hypothesize that drug adsorption to the protein and lipid-water interface accounts for a complex network of events that affect the ability of transporters to efflux drugs.