On the hydrophobicity of peptides:comparing empirical predictions of peptide log P values

Peptides are of great therapeutic potential as vaccines and drugs. Knowledge of physicochemical descriptors, including the partition coefficient logP, is useful for the development of predictive Quantitative Structure-Activity Relationships (QSARs). We have investigated the accuracy of available programs for the prediction of logP values for peptides with known experimental values obtained from the literature. Eight prediction programs were tested, of which seven programs were fragment-based methods: XLogP, LogKow, PLogP, ACDLogP, AlogP, Interactive Analysis's LogP and MlogP; and one program used a whole molecule approach: QikProp. The predictive accuracy of the programs was assessed using r(2) values, with ALogP being the most effective (r( 2) = 0.822) and MLogP the least (r(2) = 0.090). We also examined three distinct types of peptide structure: blocked, unblocked, and cyclic. For each study (all peptides, blocked, unblocked and cyclic peptides) the performance of programs rated from best to worse is as follows: all peptides - ALogP, QikProp, PLogP, XLogP, IALogP, LogKow, ACDLogP, and MlogP; blocked peptides - PLogP, XLogP, ACDLogP, IALogP, LogKow, QikProp, ALogP, and MLogP; unblocked peptides - QikProp, IALogP, ALogP, ACDLogP, MLogP, XLogP, LogKow and PLogP; cyclic peptides - LogKow, ALogP, XLogP, MLogP, QikProp, ACDLogP, IALogP. In summary, all programs gave better predictions for blocked peptides, while, in general, logP values for cyclic peptides were under-predicted and those of unblocked peptides were over-predicted.

Development, validation and uncertainty analysis of quantitative structure and activity relationship models for Log P of disinfection by -products

Hydrophobicity as measured by Log P is an important molecular property related to toxicity and carcinogenicity. With increasing public health concerns for the effects of Disinfection By-Products (DBPs), there are considerable benefits in developing Quantitative Structure and Activity Relationship (QSAR) models capable of accurately predicting Log P. In this research, Log P values of 173 DBP compounds in 6 functional classes were used to develop QSAR models, by applying 3 molecular descriptors, namely, Energy of the Lowest Unoccupied Molecular Orbital (ELUMO), Number of Chlorine (NCl) and Number of Carbon (NC) by Multiple Linear Regression (MLR) analysis. The QSAR models developed were validated based on the Organization for Economic Co-operation and Development (OECD) principles. The model Applicability Domain (AD) and mechanistic interpretation were explored. Considering the very complex nature of DBPs, the established QSAR models performed very well with respect to goodness-of-fit, robustness and predictability. The predicted values of Log P of DBPs by the QSAR models were found to be significant with a correlation coefficient R2 from 81% to 98%. The Leverage Approach by Williams Plot was applied to detect and remove outliers, consequently increasing R 2 by approximately 2% to 13% for different DBP classes. The developed QSAR models were statistically validated for their predictive power by the Leave-One-Out (LOO) and Leave-Many-Out (LMO) cross validation methods. Finally, Monte Carlo simulation was used to assess the variations and inherent uncertainties in the QSAR models of Log P and determine the most influential parameters in connection with Log P prediction. The developed QSAR models in this dissertation will have a broad applicability domain because the research data set covered six out of eight common DBP classes, including halogenated alkane, halogenated alkene, halogenated aromatic, halogenated aldehyde, halogenated ketone, and halogenated carboxylic acid, which have been brought to the attention of regulatory agencies in recent years. Furthermore, the QSAR models are suitable to be used for prediction of similar DBP compounds within the same applicability domain. The selection and integration of various methodologies developed in this research may also benefit future research in similar fields.

A chemometric study on the analgesic activity of cannabinoid compounds using SDA, KNN and SIMCA methods

The supervised pattern recognition methods K-Nearest Neighbors (KNN), stepwise discriminant analysis (SDA), and soft independent modelling of class analogy (SIMCA) were employed in this work with the aim to investigate the relationship between the molecular structure of 27 cannabinoid compounds and their analgesic activity. Previous analyses using two unsupervised pattern recognition methods (PCA-principal component analysis and HCA-hierarchical cluster analysis) were performed and five descriptors were selected as the most relevants for the analgesic activity of the compounds studied: R (3) (charge density on substituent at position C(3)), Q (1) (charge on atom C(1)), A (surface area), log P (logarithm of the partition coefficient) and MR (molecular refractivity). The supervised pattern recognition methods (SDA, KNN, and SIMCA) were employed in order to construct a reliable model that can be able to predict the analgesic activity of new cannabinoid compounds and to validate our previous study. The results obtained using the SDA, KNN, and SIMCA methods agree perfectly with our previous model. Comparing the SDA, KNN, and SIMCA results with the PCA and HCA ones we could notice that all multivariate statistical methods classified the cannabinoid compounds studied in three groups exactly in the same way: active, moderately active, and inactive.

Ester prodrugs of a potent analgesic, morphine-6-sulfate: syntheses, spectroscopic and physicochemical properties

The aim of this work is to develop 3-acyl prodrugs of the potent analgesic morphine-6-sulfate (M6S). These are expected to have higher potency and/or exhibit longer duration of analgesic action than the parent compound. M6S and the prodrugs were synthesized, then purified either by recrystallization or by semi-preparative HPLC and the structures confirmed by mass spectrometry, IR spectrophotometry and by detailed 1- and 2-D NMR studies. The lipophilicities of the compounds were assessed by a combination of shake-flask, group contribution and HPLC retention methods. The octanol-buffer partition coefficient could only be obtained directly for 3-heptanoylmorphine-6-sulfate, using the shake-flask method. The partition coefficients (P) for the remaining prodrugs were estimated from known methylene group contributions. A good linear relationship between log P and the HPLC log capacity factors was demonstrated. Hydrolysis of the 3-acetyl prodrug, as a representative of the group, was found to occur relatively slowly in buffers (pH range 6.15-8.01), with a small buffer catalysis contribution. The rates of enzymatic hydrolysis of the 3-acyl group in 10% rat blood and in 10% rat brain homogenate were investigated. The prodrugs followed apparent first order hydrolysis kinetics, with a significantly faster hydrolysis rate found in 10% rat brain homogenate than in 10% rat blood for all compounds. (C) 1998 Elsevier Science B.V. All rights reserved.

Cationic drug pharmacokinetics in diseased livers determined by fibrosis index, hepatic protein content, microsomal activity, and nature of drug

The disposition kinetics of six cationic drugs in perfused diseased and normal rat livers were determined by multiple indicator dilution and related to the drug physicochemical properties and liver histopathology. A carbon tetrachloride (CCl4)induced acute hepatocellular injury model had a higher fibrosis index (FI), determined by computer-assisted image analysis, than did an alcohol-induced chronic hepatocellular injury model. The alcohol-treated group had the highest hepatic alpha(1)- acid glycoprotein, microsomal protein (MP), and cytochrome P450 (P450) concentrations. Various pharmacokinetic parameters could be related to the octanol-water partition coefficient (log P-app) of the drug as a surrogate for plasma membrane partition coefficient and affinity for MP or P450, the dependence being lower in the CCl4-treated group and higher in the alcohol-treated group relative to controls. Stepwise regression analysis showed that hepatic extraction ratio, permeability-surface area product, tissue-binding constant, intrinsic clearance, partition ratio of influx (k(in)) and efflux rate constant (k(out)), and k(in)/k(out) were related to physicochemical properties of drug (log P-app or pK(a)) and liver histopathology (FI, MP, or P450). In addition, hepatocyte organelle ion trapping of cationic drugs was evident in all groups. It is concluded that fibrosis-inducing hepatic disease effects on cationic drug disposition in the liver may be predicted from drug properties and liver histopathology.

Bovine-serum-albumin-containing receptor phase better predicts transdermal absorption parameters for lipophilic compounds

Based on the hypothesis that limited receptor solubility of lipophilic compounds may result in lower observed permeability parameters, the aim of this study was to determine the in vitro human epidermal permeability coefficients and membrane retention of a series of aliphatic alcohols (C1-C10, log p -0.72 to 4.06) using two different receptor solutions (water and 4% bovine serum albumin in phosphate-buffered saline). Aqueous solutions of radiolabeled alcohols were dosed into the stratum corneum side of membranes mounted in side-by-side glass diffusion cells. Appearance of alcohol in the receptor compartment filled with either of the two solutions was monitored over a 7 h period when both stratum corneum (assessed by tape stripping) and the remaining epidermis levels of radioactivity were determined. In a separate study the degree of binding of alcohols to 4% bovine serum albumin was determined. The data showed increased receptor phase solubility in the bovine serum albumin solution and higher permeability coefficients for the more lipophilic alcohols in the series. No changes were seen in the partitioning of the alcohols from the vehicle into either the stratum corneum or tape-stripped epidermis with the two receptor phases; however, a decrease in the amount of the more lipophilic alcohols partitioning into the water receptor phase from the tape-stripped epidermis was observed. We conclude that bovine serum albumin receptor phase allows better estimation of real permeability parameters for lipophilic compounds due to its increased solubility capacity and we question whether permeability parameters for lipophilic solutes from older data sets based on aqueous receptor phases are completely reliable.

Effect of solute lipophilicity on penetration through canine skin

Objective To investigate the effect of lipophilicity on the percutaneous penetration of a homologous series of alcohols through canine skin Design Skin harvested from Greyhound thorax was placed in Franz-type diffusion cells and the in vitro passage of radio-labelled (C-14) alcohols (ethanol, butanol, hexanol and octanol (Log P 0.19 - 3.0)) through separate skin sections was measured in replicates of five. Permeability coefficient (k(P), cm/h), maximum flux (J(max), mol/cm(2)/h) and residue remaining within the skin were determined. Results The k(P) increased with increasing lipophilicity (6.2 x 10(-4) +/- 1.6 x 10(-4) cm/h for ethanol to 1.8 x 10(-2) 3.6 x 10(-3) cm/h for octanol). Alcohol residues remaining within each skin sample followed a similar pattern. An exponential decrease in Jmax with increasing lipophilicity was observed. Conclusion Changes in canine skin permeability occur with increasing alcohol lipophilicity. This finding has practical consequences for the design of topical formulations and optimisation of drug delivery through animal skin.

Estudo comparativo de procedimentos experimentais e computacionais para cálculo da lipofilia molecular

Introdução – A lipofilia é uma das propriedades físico-químicas que mais influencia a capacidade de uma molécula se movimentar através de compartimentos biológicos. O coeficiente de partição octanol/água (log P) permite, assim, obter uma estimativa da absorção dos fármacos no organismo. A existência de métodos indirectos para um cálculo rápido do log P pode revelar-se de grande importância na análise de listas de compostos com potencial acção farmacológica, reduzindo-as àqueles que se prevêem ter um melhor comportamento biológico. Objectivos – O propósito deste estudo é dar a conhecer um método cromatográfico de RP-HPLC desenvolvido para a determinação indirecta da lipofilia molecular e avaliar a performance de vários programas de cálculo computacional desse mesmo parâmetro. Metodologias – Seleccionaram-se 25 compostos químicos, avaliou-se o log P de cada um deles por RP-HPLC e confrontaram-se os resultados obtidos com os de sete programas computacionais. Resultados – O método RP-HPLC testado demonstrou ser vantajoso em comparação com o convencional shake flask. O programa de cálculo indirecto que proporcionou resultados mais próximos dos experimentais foi o ALOGPS© 2.1. Conclusões – A escolha ideal para a determinação da lipofilia de compostos cujo log P estimado esteja entre 0 e 6 é, sobretudo no que diz respeito à rapidez e simplicidade do processo, o método experimental indirecto RP-HPLC. Quanto aos métodos computacionais concluiu-se que nenhum dos programas, incluindo o ALOGPS© 2.1, demonstrou ser eficaz na avaliação de isómeros pelo que, para estes compostos, será sempre necessário recorrer ao método shake flask ou RP-HPLC.

Towards the validation of "in silico" models and physicochemical filters to identify and characterize new chemical entities

Summary: Lipophilicity plays an important role in the determination and the comprehension of the pharmacokinetic behavior of drugs. It is usually expressed by the partition coefficient (log P) in the n-octanol/water system. The use of an additional solvent system (1,2-dichlorethane/water) is necessary to obtain complementary information, as the log Poct values alone are not sufficient to explain ail biological properties. The aim of this thesis is to develop tools allowing to predict lipophilicity of new drugs and to analyze the information yielded by those log P values. Part I presents the development of theoretical models used to predict lipophilicity. Chapter 2 shows the necessity to extend the existing solvatochromic analyses in order to predict correctly the lipophilicity of new and complex neutral compounds. In Chapter 3, solvatochromic analyses are used to develop a model for the prediction of the lipophilicity of ions. A global model was obtained allowing to estimate the lipophilicity of neutral, anionic and cationic solutes. Part II presents the detailed study of two physicochemical filters. Chapter 4 shows that the Discovery RP Amide C16 stationary phase allows to estimate lipophilicity of the neutral form of basic and acidic solutes, except of lipophilic acidic solutes. Those solutes present additional interactions with this particular stationary phase. In Chapter 5, 4 different IANI stationary phases are investigated. For neutral solutes, linear data are obtained whatever the IANI column used. For the ionized solutes, their retention is due to a balance of electrostatic and hydrophobie interactions. Thus no discrimination is observed between different series of solutes bearing the same charge, from one column to an other. Part III presents two examples illustrating the information obtained thanks to Structure-Properties Relationships (SPR). Comparing graphically lipophilicity values obtained in two different solvent systems allows to reveal the presence of intramolecular effects .such as internai H-bond (Chapter 6). SPR is used to study the partitioning of ionizable groups encountered in Medicinal Chemistry (Chapter7). Résumé La lipophilie joue un .rôle important dans la détermination et la compréhension du comportement pharmacocinétique des médicaments. Elle est généralement exprimée par le coefficient de partage (log P) d'un composé dans le système de solvants n-octanol/eau. L'utilisation d'un deuxième système de solvants (1,2-dichloroéthane/eau) s'est avérée nécessaire afin d'obtenir des informations complémentaires, les valeurs de log Poct seules n'étant pas suffisantes pour expliquer toutes les propriétés biologiques. Le but de cette thèse est de développer des outils permettant de prédire la lipophilie de nouveaux candidats médicaments et d'analyser l'information fournie par les valeurs de log P. La Partie I présente le développement de modèles théoriques utilisés pour prédire la lipophilie. Le chapitre 2 montre la nécessité de mettre à jour les analyses solvatochromiques existantes mais inadaptées à la prédiction de la lipophilie de nouveaux composés neutres. Dans le chapitre 3, la même méthodologie des analyses solvatochromiques est utilisée pour développer un modèle permettant de prédire la lipophilie des ions. Le modèle global obtenu permet la prédiction de la lipophilie de composés neutres, anioniques et cationiques. La Partie II présente l'étude approfondie de deux filtres physicochimiques. Le Chapitre 4 montre que la phase stationnaire Discovery RP Amide C16 permet la détermination de la lipophilie de la forme neutre de composés basiques et acides, à l'exception des acides très lipophiles. Ces derniers présentent des interactions supplémentaires avec cette phase stationnaire. Dans le Chapitre 5, 4 phases stationnaires IAM sont étudiées. Pour les composés neutres étudiés, des valeurs de rétention linéaires sont obtenues, quelque que soit la colonne IAM utilisée. Pour les composés ionisables, leur rétention est due à une balance entre des interactions électrostatiques et hydrophobes. Donc aucune discrimination n'est observée entre les différentes séries de composés portant la même charge d'une colonne à l'autre. La Partie III présente deux exemples illustrant les informations obtenues par l'utilisation des relations structures-propriétés. Comparer graphiquement la lipophilie mesurée dans deux différents systèmes de solvants permet de mettre en évidence la présence d'effets intramoléculaires tels que les liaisons hydrogène intramoléculaires (Chapitre 6). Cette approche des relations structures-propriétés est aussi appliquée à l'étude du partage de fonctions ionisables rencontrées en Chimie Thérapeutique (Chapitre 7) Résumé large public Pour exercer son effet thérapeutique, un médicament doit atteindre son site d'action en quantité suffisante. La quantité effective de médicament atteignant le site d'action dépend du nombre d'interactions entre le médicament et de nombreux constituants de l'organisme comme, par exemple, les enzymes du métabolisme ou les membranes biologiques. Le passage du médicament à travers ces membranes, appelé perméation, est un paramètre important à optimiser pour développer des médicaments plus puissants. La lipophilie joue un rôle clé dans la compréhension de la perméation passive des médicaments. La lipophilie est généralement exprimée par le coefficient de partage (log P) dans le système de solvants (non miscibles) n-octanol/eau. Les valeurs de log Poct seules se sont avérées insuffisantes pour expliquer la perméation à travers toutes les différentes membranes biologiques du corps humain. L'utilisation d'un système de solvants additionnel (le système 1,2-dichloroéthane/eau) a permis d'obtenir les informations complémentaires indispensables à une bonne compréhension du processus de perméation. Un grand nombre d'outils expérimentaux et théoriques sont à disposition pour étudier la lipophilie. Ce travail de thèse se focalise principalement sur le développement ou l'amélioration de certains de ces outils pour permettre leur application à un champ plus large de composés. Voici une brève description de deux de ces outils: 1)La factorisation de la lipophilie en fonction de certaines propriétés structurelles (telle que le volume) propres aux composés permet de développer des modèles théoriques utilisables pour la prédiction de la lipophilie de nouveaux composés ou médicaments. Cette approche est appliquée à l'analyse de la lipophilie de composés neutres ainsi qu'à la lipophilie de composés chargés. 2)La chromatographie liquide à haute pression sur phase inverse (RP-HPLC) est une méthode couramment utilisée pour la détermination expérimentale des valeurs de log Poct.

Descripció de propietats moleculars i activitats biològiques emprant l'energia de repulsió electró-electró

En aquest treball es descriu l'ús de les mesures de semblança molecular quàntica (MSMQ) per a caracteritzar propietats i activitats biològiques moleculars, i definir descriptors emprables per a construir models QSAR i QSPR. L'estudi que es presenta consisteix en la continuació d'un treball recent, on es descrivien relacions entre el paràmetre log P i MSMQ, donant així una alternativa a aquest parimetre hidrofòbic empíric. L'actual contribució presenta una nova mesura, capaç d'estendre l'ús de les MSMQ, que consisteix en l'energia de repulsió electró-electró (Vee). Aquest valor, disponible normalment a partir de programari de química quàntica, considera la molècula com una sola entitat, i no cal recórrer a l'ús decontribucions de fragments. La metodologia s'ha aplicat a cinc tipus diferents de compostos on diferents propietats moleculars i activitats biològiques s'han correlacionat amb Vee com a únic descriptor molecular. En tots els casos estudiats, s'han obtingut correlacions satisfactòries.

iLOGP: A Simple, Robust, and Efficient Description of n-Octanol/Water Partition Coefficient for Drug Design Using the GB/SA Approach.

The n-octanol/water partition coefficient (log Po/w) is a key physicochemical parameter for drug discovery, design, and development. Here, we present a physics-based approach that shows a strong linear correlation between the computed solvation free energy in implicit solvents and the experimental log Po/w on a cleansed data set of more than 17,500 molecules. After internal validation by five-fold cross-validation and data randomization, the predictive power of the most interesting multiple linear model, based on two GB/SA parameters solely, was tested on two different external sets of molecules. On the Martel druglike test set, the predictive power of the best model (N = 706, r = 0.64, MAE = 1.18, and RMSE = 1.40) is similar to six well-established empirical methods. On the 17-drug test set, our model outperformed all compared empirical methodologies (N = 17, r = 0.94, MAE = 0.38, and RMSE = 0.52). The physical basis of our original GB/SA approach together with its predictive capacity, computational efficiency (1 to 2 s per molecule), and tridimensional molecular graphics capability lay the foundations for a promising predictor, the implicit log P method (iLOGP), to complement the portfolio of drug design tools developed and provided by the SIB Swiss Institute of Bioinformatics.

Partition coefficient and molecular flexibility: the concept of lipophilicity space.

The rationale of this study was to investigate molecular flexibility and its influence on physicochemical properties with a view to uncovering additional information on the fuzzy concept of dynamic molecular structure. Indeed, it is now known that computed molecular interaction fields (MIFs) such as molecular electrostatic potentials (MEPs) and lipophilicity potentials (MLPs) are conformation-dependent, as are dipole moments. A database of 125 compounds was used whose conformational space was explored, while conformation-dependent parameters were computed for each non-redundant conformer found in the conformational space of the compounds. These parameters were the virtual log P (log P(MLP), calculated by a MLP approach), the apolar surface area (ASA), polar surface area (PSA), and solvent-accessible surface (SAS). For each compound, the range taken by each parameter (its property space) was divided by the number of rotors taken as an index of flexibility, yielding a parameter termed 'molecular sensitivity'. This parameter was poorly correlated with others (i.e., it contains novel information) and showed the compounds to fall into two broad classes. 'Sensitive' molecules are those whose computed property ranges are markedly sensitive to conformational effects, whereas 'insensitive' (in fact, less sensitive) molecules have property ranges which are comparatively less affected by conformational fluctuations. A pharmacokinetic application is presented.

Estudo QSPR sobre os coeficientes de partição: descritores mecânico-quânticos e análise multivariada

Quantum chemistry and multivariate analysis were used to estimate the partition coefficients between n-octanol and water for a serie of 188 compounds, with the values of the q 2 until 0.86 for crossvalidation test. The quantum-mechanical descriptors are obtained with ab initio calculation, using the solvation effects of the Polarizable Continuum Method. Two different Hartree-Fock bases were used, and two different ways for simulating solvent cavity formation. The results for each of the cases were analised, and each methodology proposed is indicated for particular case.