Conformational changes in the nucleotide-binding domains of P-glycoprotein induced by ATP hydrolysis

Investigating how HIV-1 antiretrovirals differentially behave as substrates and inhibitors of P-glycoprotein via molecular dynamics simulations

HIV-1 can rapidly infect the brain upon initial infection, establishing latent reservoirs that induce neuronal damage and/or death, resulting in HIV-Associated Neurocognitive Disorder. Though anti-HIV-1 antiretrovirals (ARVs) suppress viral load, the blood-brain barrier limits drug access to the brain, largely because of highly expressed efflux proteins like P-glycoprotein (P-gp). While no FDA-approved P-gp inhibitor currently exists, HIV-1 protease inhibitors show promise as partial P-gp inhibitors, potentially enhancing drug delivery to the brain. Herein, we employed docking and molecular dynamics simulations to elucidate key differences in P-gp's interactions with several antiretrovirals, including protease inhibitors, with known inhibitory or substrate-like behaviors towards P-gp. Our results led us to hypothesize new mechanistic details of small-molecule efflux by and inhibition of P-gp, where the "Lower Pocket" in P-gp's transmembrane domain serves as the primary initial site for small-molecule binding. Subsequently, this pocket merges with the more traditionally studied drug binding site-the "Upper Pocket"-thus funneling small-molecule drugs, such as ARVs, towards the Upper Pocket for efflux. Furthermore, our results reinforce the understanding that both binding energetics and changes in protein dynamics are crucial in discerning small molecules as non-substrates, substrates, or inhibitors of P-gp. Our findings indicate that interactions between P-gp and inhibitory ARVs induce bridging of transmembrane domain helices, impeding P-gp conformational changes and contributing to the inhibitory behavior of these ARVs. Overall, insights gained in this study could serve to guide the design of future P-gp-targeting therapeutics for a wide range of pathological conditions and diseases, including HIV-1.

Cadmium transport by mammalian ATP-binding cassette transporters

Cellular responses to toxic metals depend on metal accessibility to intracellular targets, reaching interaction sites, and the intracellular metal concentration, which is mainly determined by uptake pathways, binding/sequestration and efflux pathways. ATP-binding cassette (ABC) transporters are ubiquitous in the human body-usually in epithelia-and are responsible for the transfer of indispensable physiological substrates (e.g. lipids and heme), protection against potentially toxic substances, maintenance of fluid composition, and excretion of metabolic waste products. Derailed regulation and gene variants of ABC transporters culminate in a wide array of pathophysiological disease states, such as oncogenic multidrug resistance or cystic fibrosis. Cadmium (Cd) has no known physiological role in mammalians and poses a health risk due to its release into the environment as a result of industrial activities, and eventually passes into the food chain. Epithelial cells, especially within the liver, lungs, gastrointestinal tract and kidneys, are particularly susceptible to the multifaceted effects of Cd because of the plethora of uptake pathways available. Pertinent to their broad substrate spectra, ABC transporters represent a major cellular efflux pathway for Cd and Cd complexes. In this review, we summarize current knowledge concerning transport of Cd and its complexes (mainly Cd bound to glutathione) by the ABC transporters ABCB1 (P-glycoprotein, MDR1), ABCB6, ABCC1 (multidrug resistance related protein 1, MRP1), ABCC7 (cystic fibrosis transmembrane regulator, CFTR), and ABCG2 (breast cancer related protein, BCRP). Potential detoxification strategies underlying ABC transporter-mediated efflux of Cd and Cd complexes are discussed.

Exploring the dynamics of the ABCB1 membrane transporter P-glycoprotein in the presence of ATP and active/non-active compounds through molecular dynamics simulations

P-glycoprotein (Pgp) is a member of the ATP-binding cassette family of transporters that confers multidrug resistance to cancer cells and is actively involved in the pharmacokinetics and toxicokinetics of a big variety of drugs. Extensive studies have provided insights into the binding of many compounds, but the precise mechanism of translocation across the membrane remains unknown; in this context, the major challenge has been to understand the basis for its polyspecificity. In this study, molecular dynamics (MD) simulations of human P-gp (hP-gp) in an explicit membrane-and-water environment were performed to investigate the dynamic behavior of the transporter in the presence of different compounds (active and inactive) in the binding pocket and ATP molecules within the nucleotide binding domains (NBDs). The complexes studied involve four compounds: cyclosporin A (CSA), amiodarone (AMI), pamidronate (APD), and valproic acid (VPA). While CSA and AMI are known to interact with P-gp, APD and VPA do not. The results highlighted how the presence of ATP notably contributed to increased flexibility of key residues in NBD1 of active systems, indicating potential conformational changes activating the translocation mechanism. MD simulations reveal how these domains adapt and respond to the presence of different substrates, as well as the influence of ATP binding on their flexibility. Furthermore, distinctive behavior was observed in the presence of active and inactive compounds, particularly in the arrangement of ATP between NBDs, supporting the proposed nucleotide sandwich dimer mechanism for ATP binding. This study provides comprehensive insights into P-gp behavior with various ligands and ATP, offering implications for drug development, toxicity assessment and demonstrating the validity of the results derived from the MD simulations.

Study Models of Drug-Drug Interactions Involving P-Glycoprotein: The Potential Benefit of P-Glycoprotein Modulation at the Kidney and Intestinal Levels

P-glycoprotein (P-gp) is a crucial membrane transporter situated on the cell's apical surface, being responsible for eliminating xenobiotics and endobiotics. P-gp modulators are compounds that can directly or indirectly affect this protein, leading to changes in its expression and function. These modulators can act as inhibitors, inducers, or activators, potentially causing drug-drug interactions (DDIs). This comprehensive review explores diverse models and techniques used to assess drug-induced P-gp modulation. We cover several approaches, including in silico, in vitro, ex vivo, and in vivo methods, with their respective strengths and limitations. Additionally, we explore the therapeutic implications of DDIs involving P-gp, with a special focus on the renal and intestinal elimination of P-gp substrates. This involves enhancing the removal of toxic substances from proximal tubular epithelial cells into the urine or increasing the transport of compounds from enterocytes into the intestinal lumen, thereby facilitating their excretion in the feces. A better understanding of these interactions, and of the distinct techniques applied for their study, will be of utmost importance for optimizing drug therapy, consequently minimizing drug-induced adverse and toxic effects.

Lipid-based nanocarriers challenging the ocular biological barriers: Current paradigm and future perspectives

Eye is the most specialized and sensory body organ and treating eye diseases efficiently is necessary. Despite various attempts, the design of a consummate ophthalmic drug delivery system remains unsolved because of anatomical and physiological barriers that hinder drug transport into the desired ocular tissues. It is important to advance new platforms to manage ocular disorders, whether they exist in the anterior or posterior cavities. Nanotechnology has piqued the interest of formulation scientists because of its capability to augment ocular bioavailability, control drug release, and minimize inefficacious drug absorption, with special attention to lipid-based nanocarriers (LBNs) because of their cellular safety profiles. LBNs have greatly improved medication availability at the targeted ocular site in the required concentration while causing minimal adverse effects on the eye tissues. Nevertheless, the exact mechanisms by which lipid-based nanocarriers can bypass different ocular barriers are still unclear and have not been discussed. Thus, to bridge this gap, the current work aims to highlight the applications of LBNs in the ocular drug delivery exploring the different ocular barriers and the mechanisms viz. adhesion, fusion, endocytosis, and lipid exchange, through which these platforms can overcome the barrier characteristics challenges.

The battle of lipid-based nanocarriers against blood-brain barrier: a critical review

Central nervous system integrity is the state of brain functioning across sensory, cognitive, emotional-social behaviors, and motor domains, allowing a person to realise his full potential. Thus, brain disorders seriously affect patients' quality of life. Efficient drug delivery to treat brain disorders remains a crucial challenge due to numerous brain barriers, particularly the blood-brain barrier (BBB), which greatly impacts the ultimate drug therapeutic efficacy. Lately, nanocarrier technology has made huge progress in overcoming these barriers by improving drug solubility, ameliorating its retention, reducing its toxicity, and targeting the encapsulated agents to different brain tissues. The current review primarily offers an overview of the different components of BBB and the progress, strategies, and contemporary applications of the nanocarriers, specifically lipid-based nanocarriers (LBNs), in treating various brain disorders.

Differential dynamics and direct interaction of bound ligands with lipids in multidrug transporter ABCG2

ABCG2 is an ATP-binding cassette (ABC) transporter that extrudes a wide range of xenobiotics and drugs from the cell and contributes to multidrug resistance in cancer cells. Following our recent structural characterization of topotecan-bound ABCG2, here, we present cryo-EM structures of ABCG2 under turnover conditions in complex with a special modulator and slow substrate, tariquidar, in nanodiscs. The structures reveal that similar to topotecan, tariquidar induces two distinct ABCG2 conformations under turnover conditions (turnover-1 and turnover-2). μs-scale molecular dynamics simulations of drug-bound and apo ABCG2 in native-like lipid bilayers, in both topotecan- and tariquidar-bound states, characterize the ligand size as a major determinant of its binding stability. The simulations highlight direct lipid-drug interactions for the smaller topotecan, which exhibits a highly dynamic binding mode. In contrast, the larger tariquidar occupies most of the available volume in the binding pocket, thus leaving little space for lipids to enter the cavity and interact with it. Similarly, when simulating ABCG2 in the apo inward-open state, we also observe spontaneous penetration of phospholipids into the binding cavity. The captured phospholipid diffusion pathway into ABCG2 offers a putative general path to recruit any hydrophobic/amphiphilic substrates directly from the membrane. Our simulations also reveal that ABCG2 rejects cholesterol as a substrate, which is omnipresent in plasma membranes that contain ABCG2. At the same time, cholesterol is found to prohibit the penetration of phospholipids into ABCG2. These molecular findings have direct functional ramifications on ABCG2's function as a transporter.

Lipid-mediated prestin organization in outer hair cell membranes and its implications in sound amplification

Prestin is a high-density motor protein in the outer hair cells (OHCs), whose conformational response to acoustic signals alters the shape of the cell, thereby playing a major role in sound amplification by the cochlea. Despite recent structures, prestin's intimate interactions with the membrane, which are central to its function remained unresolved. Here, employing a large set (collectively, more than 0.5 ms) of coarse-grained molecular dynamics simulations, we demonstrate the impact of prestin's lipid-protein interactions on its organization at densities relevant to the OHCs and its effectiveness in reshaping OHCs. Prestin causes anisotropic membrane deformation, which mediates a preferential membrane organization of prestin where deformation patterns by neighboring copies are aligned constructively. The resulting reduced membrane rigidity is hypothesized to maximize the impact of prestin on OHC reshaping. These results demonstrate a clear case of protein-protein cooperative communication in membrane, purely mediated by interactions with lipids.

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.

Membrane Mixer: A Toolkit for Efficient Shuffling of Lipids in Heterogeneous Biological Membranes

Molecular dynamics (MD) simulations of biological membranes have achieved such levels of sophistication that are commonly used to predict unresolved structures and various properties of lipids and to substantiate experimental data. While achieving sufficient sampling of lipid dynamics remains a major challenge, a commonly used method to improve lipid sampling, e.g., in terms of specific interactions with membrane-associated proteins, is to randomize the initial arrangement of lipid constituents in multiple replicas of simulations, without changing the overall lipid composition of the membrane of interest. Here, we introduce a method that can rapidly generate multiple replicas of lipid bilayers with different spatial and conformational configurations for any given lipid composition. The underlying algorithm, which allows one to shuffle lipids at any desired level, relies on the application of an external potential, here referred to as the "carving potential", that removes clashes/entanglements before lipid positions are exchanged (shuffled), thereby minimizing the energy penalty due to abrupt lipid repositioning. The method is implemented as "Membrane Mixer Plugin (MMP) 1.0" in VMD, with a convenient graphical user interface that guides the user in setting various options and parameters. The plugin is fully automated and generates new membrane replicas more rapidly and conveniently than other analogous tools. The plugin and its capabilities introduced here can be extended to include additional features in future versions.

Active participation of membrane lipids in inhibition of multidrug transporter P-glycoprotein

P-glycoprotein (Pgp) is a major efflux pump in humans, overexpressed in a variety of cancers and associated with the development of multi-drug resistance. Allosteric modulation by various ligands (e.g., transport substrates, inhibitors, and ATP) has been biochemically shown to directly influence structural dynamics, and thereby, the function of Pgp. However, the molecular details of such effects, particularly with respect to the role and involvement of the surrounding lipids, are not well established. Here, we employ all-atom molecular dynamics (MD) simulations to study the conformational landscape of Pgp in the presence of a high-affinity, third-generation inhibitor, tariquidar, in comparison to the nucleotide-free (APO) and the ATP-bound states, in order to characterize the mechanical effects of the inhibitor that might be of relevance to its blocking mechanism of Pgp. Simulations in a multi-component lipid bilayer show a dynamic equilibrium between open(er) and more closed inward-facing (IF) conformations in the APO state, with binding of ATP shifting the equilibrium towards conformations more prone to ATP hydrolysis and subsequent events in the transport cycle. In the presence of the inhibitor bound to the drug-binding pocket within the transmembrane domain (TMD), Pgp samples more open IF conformations, and the nucleotide binding domains (NBDs) become highly dynamic. Interestingly, and reproduced in multiple independent simulations, the inhibitor is observed to facilitate recruitment of lipid molecules into the Pgp lumen through the two proposed drug-entry portals, where the lipid head groups from the cytoplasmic leaflet penetrate into and, in some cases, translocate inside the TMD, while the lipid tails remain extended into the bulk lipid environment. These "wedge" lipids likely enhance the inhibitor-induced conformational restriction of the TMD leading to the differential modulation of coupling pathways observed with the NBDs downstream. We suggest a novel inhibitory mechanism for tariquidar, and potentially for related third-generation Pgp inhibitors, where lipids are seen to enhance the inhibitory role in the catalytic cycle of membrane transporters.