Evaluation and prediction of drug permeation

A major challenge confronting the pharmaceutical scientist is to optimize the selective and efficient delivery of new active entities and drug candidates. Successful drug development requires not only optimization of specific and potent pharmacodynamic activity, but also efficient delivery to the target site. Following advances in rational drug design, combinatorial chemistry and high-throughput screening techniques, the number of newly discovered and promising active compounds has increased dramatically in recent years, often making delivery problems the rate-limiting step in drug research. To overcome these problems, a good knowledge of the pharmacokinetic barriers encountered by bioactive compounds is required. This review gives an overview of the properties of relevant physiological barriers and presents some important biological models for evaluation of drug permeation and transport. Physicochemical determinants in drug permeation and the relevance of quantitative and qualitative approaches to the prediction and evaluation of passive drug absorption are also discussed.     

Collateral sensitivity of resistant MRP1-overexpressing cells to flavonoids and derivatives through GSH efflux

The multidrug resistance protein 1 (MRP1) is involved in multidrug resistance of cancer cells by mediating drug efflux out of cells, often in co-transport with glutathione (GSH). GSH efflux mediated by MRP1 can be stimulated by verapamil. In cells overexpressing MRP1, we have previously shown that verapamil induced a huge intracellular GSH depletion which triggered apoptosis of the cells. That phenomenon takes place in the more global anticancer strategy called “collateral sensitivity” and could be exploited to eradicate some chemoresistant cancer cells. Seeking alternative compounds to verapamil, we screened a library of natural flavonoids and synthetic derivatives. A large number of these compounds stimulate MRP1-mediated GSH efflux and the most active ones have been evaluated for their cytotoxic effect on MRP1-overexpressing cells versus parental cells. Interestingly, some are highly and selectively cytotoxic for MRP1-cells, leading them to apoptosis. However, some others do not exhibit any cytotoxicity while promoting a strong GSH efflux, indicating that GSH efflux is necessary but not sufficient for MRP1-cells apoptosis. In support to this hypothesis, structure activity relationships show that the absence of a hydroxyl group at position 3 of the flavonoid C ring is an absolute requirement for induction of MRP1-cells death, but is not for GSH efflux stimulation. Chrysin (compound 8) and its derivatives, compounds 11 and 22, exhibit a high selectivity toward MRP1-cells with a IC₅₀ value of 4.1 μM for compound 11 and 4.9 μM for chrysin and compound 22, making them among the best described selective killer compounds of multidrug ABC transporter-overexpressing cells.     

Non-specific SIRT inhibition as a mechanism for the cytotoxicity of ginkgolic acids and urushiols

Ginkgolic acids and urushiols are natural alkylphenols known for their mutagenic, carcinogenic and genotoxic potential. However, the mechanism of toxicity of these compounds has not been thoroughly elucidated so far. Considering that the SIRT inhibitory potential of anacardic acids has been hypothesized by in silico techniques, we herein demonstrated through both in vitro and computational methods that structurally related compounds such as ginkgolic acids and urushiols are able to modulate SIRT activity. Moreover, their SIRT inhibitory profile and cytotoxicity were comparable to sirtinol, a non-specific SIRT inhibitor (SIRT1 and SIRT2), and different from EX-527, a SIRT1 specific inhibitor. This is the first report on the SIRT inhibition of ginkgolic acids and urushiols. The results reported here are in line with previously observed effects on the induction of apoptosis by this class of compounds, and the non-specific SIRT inhibition is suggested as a new mechanism for their in vitro cytotoxicity.     

Modeling the Met Form of Human Tyrosinase: A Refined and Hydrated Pocket for Antagonist Design

Tyrosinases are type 3 copper proteins involved in melanin biosynthesis, responsible for skin and hair color in mammals. To steer tyrosinase inhibitor discovery for therapeutic and cosmetic purposes, structural information about human tyrosinase is necessary. As this protein has never been crystallized so far, we derived a robust homology model built using structural information from Streptomyces castaneoglobisporus and Ipomea batata catecholoxidase enzymes. The active site containing two copper atoms in co‐ordination with six histidine residues was refined through an optimization protocol based on molecular mechanics parameters for copper co‐ordination and charges calculated by quantum mechanics methods. Five structural water molecules and a hydroxyl ion were found to be essential for optimization. The superimposition of the human homology model on crystallographic structures of tyrosinases from other species revealed similar overall backbone topologies, active site conformations, and conserved water molecules. Phenylthiourea (PTU), the tyrosinase inhibitor of reference, was then docked into the solvated human active pocket. A binding mode consistent with crystallographic information was obtained. Taken together, these findings demonstrated that the human tyrosinase model, deposited in the Protein Model Database, is a reliable structure for future rational inhibitor design projects.     

Chirality and drugs

The two enantiomers of a chiral drug may have vastly different pharmacodynamic and pharmacokinetic properties. As a result, the research and development of chiral drugs raises specific problems some of which are discussed here. Thus, various pharmacokinetic interactions may involve two enantiomers, as seen for example when one enantiomer inhibits the metabolism of the other and modifies its effects. A different situation occurs when a third compound stereoselectively inhibits the metabolism of one of the two enantiomers. Another problem examined here results from the lack of configurational stability of some chiral drugs, a little known phenomenon whose consequences can be of pharmacological or pharmaceutical significance depending on the rate of the reaction of racemization or epimerisation. In-depth investigations are needed before choosing between a eutomer or a racemate.     

The Hydrogen-Bond: computational approaches and applications to drug design

This mini-review begins with a presentation of the hydrogen-bond in the context of other intermolecular recognition forces of significance in molecular biology and molecular pharmacology. This is followed by a survey of the various computational methods available to calculate the hydrogen-bonding capacity of compounds. Such methods use quantum mechanics, molecular mechanics, or algorithms based on experimental fragmental values. A recent and promising advance in the computation of H-bonding capacity is the development of specific molecular interaction fields (MIFs) known as molecular hydrogen-bonding potentials (MHBPs). Their interest in relating molecular properties to pharmacokinetic behaviour is highlighted with two examples, namely oral drug absorption and blood-brain barrier permeation.     

Exploration of the pharmacophore of 3-alkyl-5-arylimidazolidinediones as new CB(1) cannabinoid receptor ligands and potential antagonists: synthesis, lipophilicity, affinity, and molecular modeling

A set of 29 3-alkyl 5-arylimidazolidinediones (hydantoins) with affinity for the human cannabinoid CB(1) receptor was studied for their lipophilicity and conformational properties in order to delineate a pharmacophore. These molecules constitute a new template for cannabinoid receptor recognition, since (a) their structure differs from that of classical cannabinoid ligands and (b) antagonism is the mechanism of action of at least three compounds (20, 21, and 23). Indeed, in the [(35)S]-GTP gamma S binding assay using rat cerebellum homogenates, they behave as antagonists without any inverse agonism component. Using a set of selected compounds, experimental lipophilicity was measured by RP-HPLC and calculated by a fragmental method (CLOGP) and a conformation-dependent method (CLIP based on the molecular lipophilicity potential). These approaches revealed two models which differentiate the binding mode of nonpolar and polar hydantoins and which could explain, at least for compounds 20, 21, and 23, the mechanism of action of this new family of cannabinoid ligands.     

Liposome/water lipophilicity: methods, information content, and pharmaceutical applications

This review discusses liposome/water lipophilicity in terms of the structure of liposomes, experimental methods, and information content. In a first part, the structural properties of the hydrophobic core and polar surface of liposomes are examined in the light of potential interactions with solute molecules. Particular emphasis is placed on the physicochemical properties of polar headgroups of lipids in liposomes. A second part is dedicated to three useful methods to study liposome/water partitioning, namely potentiometry, equilibrium dialysis, and (1)H-NMR relaxation rates. In each case, the principle and limitations of the method are discussed. The next part presents the structural information encoded in liposome/water lipophilicity, in other words the solutes' structural and physicochemical properties that determine their behavior and hence their partitioning in such systems. This presentation is based on a comparison between isotropic (i.e., solvent/water) and anisotropic (e.g., liposome/water) systems. An important factor to be considered is whether the anisotropic lipid phase is ionized or not. Three examples taken from the authors' laboratories are discussed to illustrate the factors or combinations thereof that govern liposome/water lipophilicity, namely (a) hydrophobic interactions alone, (b) hydrophobic and polar interactions, and (c) conformational effects plus hydrophobic and ionic interactions. The next part presents two studies taken from the field of QSAR to exemplify the use of liposome/water lipophilicity in structure-disposition and structure-activity relationships. In the conclusion, we summarize the interests and limitations of this technology and point to promising developments.     

Toxication of MPTP (1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine) and analogs by monoamine oxidase. A structure-reactivity relationship study

MPTP (1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine) elicits motor deficits similar to those observed in Parkinson's disease. Before exerting its neurotoxic action, MPTP must be activated by brain monoamine oxidase (MAO) to the neurotoxic metabolite MPP+ (1-methyl-4-phenylpyridinium). MPTP derivatives differ in their reactivity as MAO substrates and in their neurotoxicity. A structure-reactivity relationship study based on literature data was undertaken in order to determine the key features in the structure of MPTP and analogs that are responsible for the reactivity towards MAO. Thirty-three MPTP derivatives (including MPTP itself) were included in the study. To explain the reactivity towards MAO of the 33 MPTP analogs, different statistical methods (principal component analysis, multiple linear regression analysis) as well as the CoMFA (Comparative Molecular Field Analysis) approach, a new tool in structure-activity correlations, were used. Linear regression analysis failed to yield any predictive model, but suggested some trends. In contrast, the CoMFA approach was successful in correlating structural features and MAO reactivity. Coefficient contour maps showed where differences in the steric field (van der Waals' interactions) are most highly associated with differences in MAO reactivity. Several positive (in the ortho- and meta-position of the phenyl group) and negative (in the para-position of the phenyl group; beyond the N-methyl group) interaction regions were identified. Some structural features of the MAO active site could be postulated. First, the N-methyl group has the ideal size and elicits ideal interactions within the MAO active pocket, while smaller or larger groups are less favorable; second, para-substituent on the phenyl ring produce steric hindrances and are unfavorable to reactivity; third, ortho- and meta-substituents may have stabilizing interactions within the active pocket and are favorable to the reactivity. Moreover the model derived by CoMFA allowed us to make successful predictions of reactivity towards MAO for several additional tetrahydropyridines.