Lipid-mediated activation of plasma membrane-localized deubiquitylating enzymes modulate endosomal trafficking.

The abundance of plasma membrane-resident receptors and transporters has to be tightly regulated by ubiquitin-mediated endosomal degradation for the proper coordination of environmental stimuli and intracellular signaling. Arabidopsis OVARIAN TUMOR PROTEASE (OTU) 11 and OTU12 are plasma membrane-localized deubiquitylating enzymes (DUBs) that bind to phospholipids through a polybasic motif in the OTU domain. Here we show that the DUB activity of OTU11 and OTU12 towards K63-linked ubiquitin is stimulated by binding to lipid membranes containing anionic lipids. In addition, we show that the DUB activity of OTU11 against K6- and K11-linkages is also stimulated by anionic lipids, and that OTU11 and OTU12 can modulate the endosomal degradation of a model cargo and the auxin efflux transporter PIN2-GFP in vivo. Our results suggest that the catalytic activity of OTU11 and OTU12 is tightly connected to their ability to bind membranes and that OTU11 and OTU12 are involved in the fine-tuning of plasma membrane proteins in Arabidopsis.

Knowledge-driven design and optimization of potent symmetric anticancer molecules: A case study on PKM2 activators.

BACKGROUND: Pyruvate kinase M2 (PKM2) is preferentially expressed as a low-activity dimer over the active tetramer in proliferating tumor cells, resulting in metabolic reprogramming to achieve high energy requirements and nutrient uptake. This leads to a shift from the normal glycolytic pathway causing tumor cells to proliferate uncontrollably. This study utilizes knowledge-based drug discovery to determine the critical features from experimentally known PKM2 activators and design compounds that would significantly confer a stable structural and functional edge over the known compounds which are still at the preclinical stage. METHODS: Conscientious molecular modeling studies were carried out and critical structural features were identified and validated from the knowledge of experimentally known PKM2 activators to confer high-binding affinities. A virtual library of 200 palindromic and non-palindromic activators was designed based on these identified critical features to target a distinct activator binding-site. This binding would favor specific dimer-dimer association and subsequent protein tetramerization. The resultant compounds strongly correlated with identified structural features and binding affinities which further strengthened our findings. The designed activators were then subjected to pharmacokinetic profiling and toxicity prediction, followed by free-binding energy calculations and MD simulations. RESULTS: All the virtually designed activators comprising the identified critical features were observed to confer high-binding affinities ranging from -9.1 to -15.0 kcal/mol to the receptor protein. The designed activators also demonstrated optimum pharmacokinetic and toxicity profiles. CONCLUSION: The best activators selected for MD simulations studies were conclusively observed to stabilize the required tetrameric conformation suggesting that these activators could potentially target PKM2 tetramerization that might restore the normal glycolytic pathway and suppress tumor progression.

Electrostatic and steric effects underlie acetylation-induced changes in ubiquitin structure and function.

Covalent attachment of ubiquitin (Ub) to proteins is a highly versatile posttranslational modification. Moreover, Ub is not only a modifier but itself is modified by phosphorylation and lysine acetylation. However, the functional consequences of Ub acetylation are poorly understood. By generation and comprehensive characterization of all seven possible mono-acetylated Ub variants, we show that each acetylation site has a particular impact on Ub structure. This is reflected in selective usage of the acetylated variants by different E3 ligases and overlapping but distinct interactomes, linking different acetylated variants to different cellular pathways. Notably, not only electrostatic but also steric effects contribute to acetylation-induced changes in Ub structure and, thus, function. Finally, we provide evidence that p300 acts as a position-specific Ub acetyltransferase and HDAC6 as a general Ub deacetylase. Our findings provide intimate insights into the structural and functional consequences of Ub acetylation and highlight the general importance of Ub acetylation.

Structural basis of Omicron immune evasion: A comparative computational study.

BACKGROUND: The vaccines used against SARS-CoV-2 by now have been able to develop some neutralising antibodies in the vaccinated population and their effectiveness has been challenged by the emergence of the new strains with numerous mutations in the spike protein of SARS-CoV-2. Since S protein is the major immunogenic protein of the virus which contains Receptor Binding Domain (RBD) that interacts with the human Angiotensin-Converting Enzyme 2 (ACE2) receptors, any mutations in this region should affect the neutralisation potential of the antibodies leading to the immune evasion. Several variants of concern of the virus have emerged so far, amongst which the most critical are Delta and recently reported Omicron. In this study, we have mapped and reported mutations on the modelled RBD and evaluated binding affinities of various human antibodies with it. METHOD: Docking and molecular dynamics simulation studies have been used to explore the effect of mutations on the structure of RBD and RBD-antibody interaction. RESULTS: These analyses show that the mutations mostly at the interface of a nearby region lower the binding affinity of the antibody by ten to forty percent, with a downfall in the number of interactions formed as a whole. It implies the generation of immune escape variants. CONCLUSIONS: Notable mutations and their effect was characterised that explain the structural basis of antibody efficacy in Delta and a compromised neutralisation effect for the Omicron variant. Thus, our results pave the way for robust vaccine design that can be effective for many variants.

Novel imidazopyrimidines-based molecules induce tetramerization of tumor pyruvate kinase M2 and exhibit potent antiproliferative profile.

Discovery of novel and potent lead molecules for the specific therapeutic targets by de novo drug design is still in infancy. Here, we disclose the unprecedented development of imidazopyri(mi)dine-based tumor pyruvate kinase M2 (PKM2) modulators by subsequent link and grow strategy. The most potent modulator 15n acts as a PKM2 activator with an AC50 of 90 nM, with considerable cancer cell-selectivity and membrane-permeability. NMR metabolomics studies also revealed that treatment with 15n results in diminution in lactate concentrations in MCF-7 cells. 15n binds to a previously reported site at PKM2 adjacent to the interface of two monomers. In molecular dynamics (MD) simulation studies, it was observed that 15n stabilizes the PKM2 at the dimeric interface, assisting in the formation of a biologically active tetramer conformation. 15n was also screened on MCF-7 breast cancer cell lines grown on 3-D scaffolds, and the results exhibited better anticancer potential compared to control, paving the way for future clinical studies.

Discovery of boronic acid-based potent activators of tumor pyruvate kinase M2 and development of gastroretentive nanoformulation for oral dosing.

Several studies have established that cancer cells explicitly over-express the less active isoform of pyruvate kinase M2 (PKM2) is critical for tumorigenesis. The activation of PKM2 towards tetramer formation may increase affinity towards phosphoenolpyruvate (PEP) and avoidance of the Warburg effect. Herein, we describe the design, synthesis, and development of boronic acid-based molecules as activators of PKM2. The designed molecules were inspired by existing anticancer scaffolds and several fragments were assembled in the derivatives. 6a-6d were synthesized using a multi-step synthetic strategy in 55-70% yields, starting from cheap and readily available materials. The compounds were selectively cytotoxic to kill the cancerous cells at 80 nM, while they were non-toxic to the normal cells. The kinetic studies established the compounds as novel activators of PKM2 and (E/Z)-(4-(3-(2-((4-chlorophenyl)amino)-4-(dimethylamino)thiazol-5-yl)-2-(ethoxycarbonyl)-3-oxoprop-1-en-1-yl) phenyl)boronic acid (6c) emerged as the most potent derivative. 6c was further evaluated using various in silico tools to understand the molecular mechanism of tetramer formation. Docking studies revealed that 6c binds to the PKM2 dimer at the dimeric interface. Further to ascertain the binding site and mechanism of action, rigorous MD (molecular dynamics) simulations were undertaken, which led to the conclusion that 6c stabilizes the center of the dimeric interface that possibly promotes tetramer formation. We further planned to make a tablet of the developed molecule for oral delivery, but it was seriously impeded owing to poor aqueous solubility of 6c. To improve aqueous solubility and retain 6c at the lower gastrointestinal tract, thiolated chitosan-based nanoparticles (TCNPs) were prepared and further developed as tablet dosage form to retain anticancer potency in the excised goat colon. Our findings may provide a valuable pharmacological mechanism for understanding metabolic underpinnings that may aid in the clinical development of new anticancer agents targeting PKM2.

Deciphering molecular details of the RAC-ribosome interaction by EPR spectroscopy.

The eukaryotic ribosome-associated complex (RAC) plays a significant role in de novo protein folding. Its unique interaction with the ribosome, comprising contacts to both ribosomal subunits, suggests a RAC-mediated coordination between translation elongation and co-translational protein folding. Here, we apply electron paramagnetic resonance (EPR) spectroscopy combined with site-directed spin labeling (SDSL) to gain deeper insights into a RAC-ribosome contact affecting translational accuracy. We identified a local contact point of RAC to the ribosome. The data provide the first experimental evidence for the existence of a four-helix bundle as well as a long α-helix in full-length RAC, in solution as well as on the ribosome. Additionally, we complemented the structural picture of the region mediating this functionally important contact on the 40S ribosomal subunit. In sum, this study constitutes the first application of SDSL-EPR spectroscopy to elucidate the molecular details of the interaction between the 3.3 MDa translation machinery and a chaperone complex.

Interactions Determining the Structural Integrity of the Trimer of Plant Light Harvesting Complex in Lipid Membranes.

The structural basis for the stability of the trimeric form of the light harvesting complex (LHCII), a pigmented protein from green plants pivotal for photosynthesis, remains elusive till date. The protein embedded in a dipalmitoylphosphatidylcholine (DPPC) lipid membrane is investigated using all-atom molecular dynamics simulations to find out the interactions responsible for the structural integrity of the trimer and its relation to antenna function. Central association of chlorophyll a (CLA) molecules near the LHCII chains is attributed to a conserved coordination between the Mg of CLA and the oxygen of a specific residue of the first helix of a chain. The residue forms a salt-bridge with the fourth helix of the same chain of the trimer, not of the monomer. In an earlier experiment, three residues (WYR) at each chain of the trimer have been found indispensable for the trimerization and referred to as trimerization motif. We find that the residues of the trimerization motif are connected to the lipids or pigments by a chain of interactions rather than a direct contact. Synergistic effects of sequentially located hydrogen bonds and salt-bridges within monomers of the trimer keep the trimer conformation stable in association with the pigments or the lipids. These interactions are exclusively present in the pigmented trimer and not present in the monomer or in the unpigmented trimer. Thus, our results provide a molecular basis for the inherent stability of the LHCII trimer in a lipid membrane and explain many pre-existing experimental data.

A Ligand Selection Strategy Identifies Chemical Probes Targeting the Proteases of SARS-CoV-2.

Activity-based probes are valuable tools for chemical biology. However, finding probes that specifically target the active site of an enzyme remains a challenging task. Herein, we present a ligand selection strategy that allows to rapidly tailor electrophilic probes to a target of choice and showcase its application for the two cysteine proteases of SARS-CoV-2 as proof of concept. The resulting probes were specific for the active site labeling of 3CLpro and PLpro with sufficient selectivity in a live cell model as well as in the background of a native human proteome. Exploiting the probes as tools for competitive profiling of a natural product library identified salvianolic acid derivatives as promising 3CLpro inhibitors. We anticipate that our ligand selection strategy will be useful to rapidly develop customized probes and discover inhibitors for a wide range of target proteins also beyond corona virus proteases.

Coarse-Grained Simulations of Peptide Nanoparticle Formation: Role of Local Structure and Nonbonded Interactions.

Biocompatible nanostructures play an important role in drug delivery and tissue engineering applications. Controlled growth of peptide-based nanoparticles with specific morphology needs an understanding of the role of the sequence and solvation properties. In a previous combined experimental-computational study, we identified factors that govern the formation of well-defined aggregates by self-assembled pentapeptides using single amino acid substitution ( Mishra , N. K. ; Jain , A. ; Peter , C. ; Verma , S. J. Phys. Chem. B 2017 , 121 , 8155 - 8161 ). The atomistic simulation study suggested a subtle interplay between various peptide properties like rigidity/flexibility, hydrogen bonding, partitioning of aromatic residues, and dimerization of peptides that determine the different morphologies, while the overall aggregation propensity was mostly determined by the composition of the methanol/water solvent mixture. The size of the simulated aggregates and the time scales were rather restricted due to the atomistic character of the study. Here, we present an extension to a coarse-grained representation that allows for much larger system sizes and longer time scales. To this end, we have optimized a MARTINI model so that it can deal with a system that relies on local structure formation. We combine information on local behavior from atomistic studies and apply supportive dihedral angles together with local adjustment of the bead types to find the right interplay of solvent and peptides. Finally, to mimic the dimers, an introduction of additional bonds between the monomers was necessary. By adding the modifications stepwise, we were able to disentangle the influences of the various contributions, like the rigidity/flexibility of the peptides, dimer formation, or nonbonded properties of the beads, on the overall aggregation propensity and morphology of the nanoparticles. The obtained models resemble the experimental and atomistic behavior and are able to provide mechanistic insight into peptide nanoparticle formation.

Tipping the Scale from Disorder to Alpha-helix: Folding of Amphiphilic Peptides in the Presence of Macroscopic and Molecular Interfaces.

Secondary amphiphilicity is inherent to the secondary structural elements of proteins. By forming energetically favorable contacts with each other these amphiphilic building blocks give rise to the formation of a tertiary structure. Small proteins and peptides, on the other hand, are usually too short to form multiple structural elements and cannot stabilize them internally. Therefore, these molecules are often found to be structurally ambiguous up to the point of a large degree of intrinsic disorder in solution. Consequently, their conformational preference is particularly susceptible to environmental conditions such as pH, salts, or presence of interfaces. In this study we use molecular dynamics simulations to analyze the conformational behavior of two synthetic peptides, LKKLLKLLKKLLKL (LK) and EAALAEALAEALAE (EALA), with built-in secondary amphiphilicity upon forming an alpha-helix. We use these model peptides to systematically study their aggregation and the influence of macroscopic and molecular interfaces on their conformational preferences. We show that the peptides are neither random coils in bulk water nor fully formed alpha helices, but adopt multiple conformations and secondary structure elements with short lifetimes. These provide a basis for conformation-selection and population-shift upon environmental changes. Differences in these peptides' response to macroscopic and molecular interfaces (presented by an aggregation partner) can be linked to their inherent alpha-helical tendencies in bulk water. We find that the peptides' aggregation behavior is also strongly affected by presence or absence of an interface, and rather subtly depends on their surface charge and hydrophobicity.

Sticky water surfaces: helix-coil transitions suppressed in a cell-penetrating peptide at the air-water interface.

GALA is a 30 amino acid synthetic peptide consisting of a Glu-Ala-Leu-Ala repeat and is known to undergo a reversible structural transition from a disordered to an α-helical structure when changing the pH from basic to acidic values. In its helical state GALA can insert into and disintegrate lipid membranes. This effect has generated much interest in GALA as a candidate for pH triggered, targeted drug delivery. GALA also serves as a well-defined model system to understand cell penetration mechanisms and protein folding triggered by external stimuli. Structural transitions of GALA in solution have been studied extensively. However, cell penetration is an interfacial effect and potential biomedical applications of GALA would involve a variety of surfaces, e.g., nanoparticles, lipid membranes, tubing, and liquid-gas interfaces. Despite the apparent importance of interfaces in the functioning of GALA, the effect of surfaces on the reversible folding of GALA has not yet been studied. Here, we use sum frequency generation vibrational spectroscopy (SFG) to probe the structural response of GALA at the air-water interface and IR spectroscopy to follow GALA folding in bulk solution. We combine the SFG data with molecular dynamics simulations to obtain a molecular-level picture of the interaction of GALA with the air-water interface. Surprisingly, while the fully reversible structural transition was observed in solution, at the water-air interface, a large fraction of the GALA population remained helical at high pH. This "stickiness" of the air-water interface can be explained by the stabilizing interactions of hydrophobic leucine and alanine side chains with the water surface.

Optimization of an elastic network augmented coarse grained model to study CCMV capsid deformation.

The major protective coat of most viruses is a highly symmetric protein capsid that forms spontaneously from many copies of identical proteins. Structural and mechanical properties of such capsids, as well as their self-assembly process, have been studied experimentally and theoretically, including modeling efforts by computer simulations on various scales. Atomistic models include specific details of local protein binding but are limited in system size and accessible time, while coarse grained (CG) models do get access to longer time and length scales but often lack the specific local interactions. Multi-scale models aim at bridging this gap by systematically connecting different levels of resolution. Here, a CG model for CCMV (Cowpea Chlorotic Mottle Virus), a virus with an icosahedral shell of 180 identical protein monomers, is developed, where parameters are derived from atomistic simulations of capsid protein dimers in aqueous solution. In particular, a new method is introduced to combine the MARTINI CG model with a supportive elastic network based on structural fluctuations of individual monomers. In the parametrization process, both network connectivity and strength are optimized. This elastic-network optimized CG model, which solely relies on atomistic data of small units (dimers), is able to correctly predict inter-protein conformational flexibility and properties of larger capsid fragments of 20 and more subunits. Furthermore, it is shown that this CG model reproduces experimental (Atomic Force Microscopy) indentation measurements of the entire viral capsid. Thus it is shown that one obvious goal for hierarchical modeling, namely predicting mechanical properties of larger protein complexes from models that are carefully parametrized on elastic properties of smaller units, is achievable.

Coarse-Grained and Atomistic Simulations of the Salt-Stable Cowpea Chlorotic Mottle Virus (SS-CCMV) Subunit 26-49: ß-Barrel Stability of the Hexamer and Pentamer Geometries.

A combination of coarse-grained (CG) and atomistic simulations provides a suitable computational framework to study unstructured regions of proteins, for which experimental data are often lacking or limited. In this work, we combine CG and atomistic simulations with clustering algorithms and free energy reweighting methods to explore the conformational equilibrium of certain regions of the salt-stable cowpea chlorotic mottle virus (SS-CCMV). In particular, we focus on the geometry of converging strands (residues 26-49) from contacting subunits at the 3-fold (hexamer) and 5-fold (pentamer) symmetry points of the capsid. We show the following: (i) The simulations reproduce the experimentally observed ß-barrel for the hexamer. (ii) The pentamer geometry is unable to stabilize a ß-barrel conformation; it assumes various states instead, again in accordance with the experimental results which do not indicate a well-defined structure for the pentameric interface. (iii) Atomistic simulations of the backmapped CG structures remain relatively stable, indicative of plausible CG conformations and slow kinetics on the atomistic level.

Comparison of the inward- and outward-open homology models and ligand binding of human P-glycoprotein.

An homology model of human P-glycoprotein, based on the X-ray structure of the recently resolved mouse P-glycoprotein, is presented. The model corresponds to the inward-facing conformation competent for drug binding. From the model, the residues involved in the protein-binding cavity are identified and compared with those in the outward-facing conformation of human P-glycoprotein developed previously based on the Sav1866 structure. A detailed analysis of the interactions of the cyclic peptides QZ59-RRR and QZ59-SSS is presented in both the X-ray structures of mouse P-glycoprotein and the human P-glycoprotein model generated by ligand docking. The results confirm the functional role of transmembrane domains TM4, TM6, TM10 and TM12 as entrance gates to the protein cavity, and also imply differences in their functions. The analysis of the cavities in both models suggests that the ligands remain bound to the same residues during the transition from the inward- to the outward-facing conformations. The analysis of the ligand-protein interactions in the X-ray complexes shows differences in the residues involved, as well as in the specific interactions performed by the same ligand within the same protein. This observation is supported by docking of the QZ59 ligands into human P-glycoprotein, thus aiding in the understanding of the complex behavior of P-glycoprotein substrates and inhibitors. The results confirm the possibility for multispecific drug interactions of the protein, and are important for elucidating the P-glycoprotein function and ligand interactions.

Combined pharmacophore modeling, docking, and 3D QSAR studies of ABCB1 and ABCC1 transporter inhibitors.

Quinazolinones, indolo- and pyrrolopyrimidines with inhibitory effects toward ABCB1 (P-gp) and ABCC1 (MRP1) transporters were studied by pharmacophore modeling, docking, and 3D QSAR to describe the binding preferences of the proteins. The pharmacophore overlays between dual and/or highly selective inhibitors point to binding sites of different topology and physiochemical properties for MRP1 and P-gp. Docking of selective inhibitors into the P-gp binding cavity by the use of a structural model based on the recently resolved P-gp structure confirms the P-gp pharmacophore features identified, and reveals the interactions of some functional groups and atoms in the structures with particular protein residues. The 3D QSAR analysis of the dual-effect inhibitors allows satisfactory prediction of the selectivity index of the compounds and outlines electrostatics as most important for selectivity. The results from the combined modeling approach complement each other and could improve our understanding of the protein-ligand interactions involved, and could aid in the development of highly selective and potent inhibitors of P-gp and MRP1.

Synthesis and biological evaluation of a small molecule library of 3rd generation multidrug resistance modulators.

The development of new modulators possessing high efficacy, low toxicity and high selectivity is a pivotal approach to overcoming P-glycoprotein (P-gp) mediated multidrug resistance (MDR) in tumour cells. In this study 39 compounds are presented which have been synthesized and pharmacologically investigated in our laboratory. Similarly to the potent 3rd generation MDR modulator tariquidar (XR9576) the compounds contain a tetrahydroisoquinoline-ethyl-phenylamine substructure that, in contrast to XR9576, is connected to a smaller hydrophobic part, thus leading to molecules of lower molecular weight. The connection between the tetrahydroisoquinoline-ethyl-phenylamine substructure and the hydrophobic part was achieved through four different types of linkers: amide, urea, amide-ether and amide-styryl. A number of structural modifications in the hydrophobic part were created. The calcein AM assay served as test system to determine the P-gp transport inhibitory potencies of the compounds. For the amide linker derivatives a structure-activity relationship analysis was performed outlining which structural modifications contributed to the inhibitory potency. The compounds containing a bicyclic hydrophobic part with a particular substituent in a specific orientation were identified as the most potent amide derivatives. Among the urea derivatives the compounds with highest inhibitory potency possessed an ortho-nitro substituent. The conformational analysis revealed that this position enables the formation of a hydrogen bond to the urea linker thus stabilizing the conformation. Regarding the amide-styryl derivatives the elongation of the amide linker seemed to be most decisive for the observed increase in activity. The most promising candidate in the whole library possess an amide-ether linker and an ortho-nitro substituent in the hydrophobic part. This compound inhibites P-gp slightly less than tariquidar and can serve as a lead structure for new potent P-gp modulators.

Methoxylation of 3',4'-aromatic side chains improves P-glycoprotein inhibitory and multidrug resistance reversal activities of 7,8-pyranocoumarin against cancer cells.

The overexpression of P-glycoprotein (Pgp), an ATP-driven membrane exporter of hydrophobic xenobiotics, is one of the major causes of multidrug resistance (MDR) in cancer cells. Through extensive screening we have found that the extracts of Peucedanum praeruptorum Dunn. and one of the major components (+/-)-praeruptorin A (PA) may reverse Pgp-mediated multidrug resistance. Studies on novel PA derivatives have shown that (+/-)-3'-O,4'-O-dicinnamoyl-cis-khellactone (DCK) is more active than PA or verapamil and is a non-competitive inhibitor of Pgp. Here, we report that methoxylation of the cinnamoyl groups on DCK may further enhance its bioactivity. The structure-activity relationship is demonstrated by comparing two new pyranocoumarins (+/-)-3'-O,4'-O-bis(3,4-dimethoxycinnamoyl)-cis-khellactone (DMDCK) and (+/-)-3'-O,4'-O-bis(4-methoxycinnamoyl)-cis-khellactone (MMDCK). While the co-existence of 3- and 4-methoxy groups on cinnamoyl remarkably enhanced the Pgp-inhibitory activity, the lone existence of the 4-methoxy group on cinnamoyl reduced the activity. Contrary to DCK, DMDCK promoted the binding of UIC2 antibody to Pgp which signifies a conformational change of Pgp similar to that induced by transport substrates. While DCK moderately stimulated the basal Pgp-ATPase activity, DMDCK inhibited the activity. A pharmacophore search with verapamil-based template revealed that four functional groups of DMDCK could be simultaneously involved in interaction with Pgp whereas for DCK or MMDCK only three groups were involved. It is speculated that the additional 3-methoxy group on cinnamoyl allows DMDCK to interact more efficiently with Pgp substrate site(s). If DMDCK was tightly bind to Pgp substrate site(s) the complexes could be inactive with regard to transportation and ATP hydrolysis could also be inhibited.

Identification of putative binding sites of P-glycoprotein based on its homology model.

A homology model of P-glycoprotein based on the crystal structure of the multidrug transporter Sav1866 is developed, incorporated into a membrane environment, and optimized. The resulting model is analyzed in relation to the functional state and potential binding sites. The comparison of modeled distances to distances reported in experimental studies between particular residues suggests that the model corresponds most closely to the first ATP hydrolysis step of the protein transport cycle. Comparison to the protein 3D structure confirms this suggestion. Using SiteID and Site Finder programs three membrane related binding regions are identified: a region at the interface between the membrane and cytosol and two regions located in the transmembrane domains. The regions contain binding pockets of different size, orientation, and amino acids. A binding pocket located inside the membrane cavity is also identified. The pockets are analyzed in relation to amino acids shown experimentally to influence the protein function. The results suggest that the protein has multiple binding sites and may bind and/or release substrates in multiple pathways.

Functional assay and structure-activity relationships of new third-generation P-glycoprotein inhibitors.

Twenty-eight compounds, including 24 structurally related derivatives of tariquidar synthesized in our laboratory, and four XR compounds, reported by Xenova group Ltd, were investigated by the Hoechst 33342 and Calcein AM functional assays for estimation of their inhibitory effects on the transport activity of P-glycoprotein (P-gp). A high correlation between the effects obtained in both assays was observed at the substrate concentrations used. The analyses of kinetics data from experiments at different substrate concentrations revealed non-competitive inhibition in the Calcein AM assay and competitive inhibition in the Hoechst 33342 assay. The 3D structures of the compounds were further aligned on Hoechst 33342 using flexible and pharmacophore alignments. The results suggested that inhibitors could interact with the H-binding site of P-gp and this could potentially be achieved by different ways of binding. The best 3D-QSAR models, generated by CoMFA and CoMSIA, yielded an internal predictive squared correlation coefficient higher than 0.8 and included electrostatic, steric, hydrogen bond acceptor, and hydrophobic fields. Validation of the models on an external test set of 30 XR compounds gave predictive squared correlation coefficients of up to 0.66. An excellent correspondence between the experimental and modeled activities of the test compounds was observed. The models can be used for prediction and rational design of new P-gp inhibitors.