Development of a physiologically based toxicokinetic model for butadiene and four major metabolites in humans: global sensitivity analysis for experimental design issues

1,3-Butadiene (BD) is metabolized in humans and rodents to mutagenic and carcinogenic species. Our previous work has focused on developing a physiologically based toxicokinetic (PBTK) model for BD to estimate its metabolic rate to 1,2-epoxy-3-butene (EB), using exhaled breath BD concentrations in human volunteers exposed by inhalation. In this paper, we extend our BD model to describe the kinetics of its four major metabolites EB, 1,2:3,4-diepoxybutane (DEB), 3-butene-1,2-diol (BDD), and 3,4-epoxy-1,2-butanediol (EBD), and to test whether the extended model and experimental data (to be collected for BD and metabolites in humans) are together adequate to estimate the metabolic rate constants of each of the above chemicals. Global sensitivity analyses (GSA) were conducted to evaluate the relative importance of the model parameters on model outputs during the 20min of exposure and the 40min after exposure ended. All model parameters were studied together with various potentially measurable model outputs: concentrations of BD and EB in exhaled air, concentrations of BD and all metabolites in venous blood, and cumulated amounts of urinary metabolites excreted within 24h. Our results show that pulmonary absorption of BD and subsequent distribution and metabolism in the well-perfused tissues compartment are the critical processes in the toxicokinetics of BD and metabolites. In particular, three parameters influence numerous outputs: the blood:air partition coefficient for BD, the metabolic rate of BD to EB, and the volume of the well-perfused tissues. Other influential parameters include other metabolic rates, some partition coefficients, and parameters driving the gas exchanges (in particular, for BD outputs). GSA shows that the impact of the metabolic rate of BD to EB on the BD concentrations in exhaled air is greatly increased if a few of the model's important parameters (such as the blood:air partition coefficient for BD) are measured experimentally. GSA also shows that all the transformation pathways described in the PBTK model may not be estimable if only data on the studied outputs are collected, and that data on a specific output for a chemical may not inform all the transformations involving that chemical.     

Lipophilic phenolic antioxidants: Correlation between antioxidant profile,partition coefficients and redox properties

Lipophilic compounds structurally based on caffeic, hydrocaffeic, ferulic and hydroferulic acids were synthesized. Subsequently, their antioxidant activity was evaluated as well as their partition coefficients and redox potentials. The structure–property–activity relationship (SPAR) results revealed the existence of a clear correlation between the redox potentials and the antioxidant activity. In addition, some compounds showed a proper lipophilicity to cross the blood–brain barrier. Their predicted ADME properties are also in accordance with the general requirements for potential CNS drugs. Accordingly, one can propose these phenolic compounds as potential antioxidants for tackling the oxidative status linked to the neurodegenerative processes.     

Erythrocyte ghost--buffer partition coefficients of phenobarbital, pentobarbital, and thiopental support the pH-partition hypothesis.

The partition coefficient (lambda) between red cell ghosts and buffer has been determined for three barbiturates over a range of pH. Experimental partition coefficients were linearly proportional to the calculated degree of association of the barbiturates. Lambda was 9.5 +/- 0.52 for phenobarbital, 12.7 +/- 0.91 for pentobarbital, and 27 +/- 4.9 for thiopental in their acid forms. Lambda for all three barbiturates in their anionic forms was zero. Our data support the assumption of the pH-partition hypothesis that the dependence of lambda on pH in biological membranes behaves essentially like that in organic solvents. However, the relative magnitudes of the erythrocyte partition coefficients correlate much more closely with the physiological permeability constants than do those of organic solvents, which tend to overestimate the differences between these compounds.     

Partition coefficients of beta-blockers in bile salt/lecithin micelles as a tool to assess the role of mixed micelles in gastrointestinal absorption

The objective of this study was to develop non-invasive spectroscopic methods to quantify the partition coefficients of two beta-blockers, atenolol and nadolol, in aqueous solutions of bile salt micelles and to assess the effect of lecithin on the partition coefficients of amphiphilic drugs in mixed bile salt/lecithin micelles, which were used as a simple model for the naturally occurring mixed micelles in the gastrointestinal tract. The partition coefficients (Kp) at 25.0 +/- 0.1degreesC and at 0.1 M NaCl ionic strength were determined by spectrofluorimetry and by derivative spectrophotometry, by fitting equations that relate molar extinction coefficients and relative fluorescence intensities to the partition constant Kp. Drug partition was controlled by the: (i) drug properties, with the more soluble drug in water (atenolol) exhibiting smaller values of Kp, and with both drugs interacting more extensively in the protonated form; and by (ii) the bile salt monomers, with the dihydroxylic salts producing larger values of Kp for the beta-blockers, and with glycine conjugation of the bile acid increasing the values of Kp for the beta-blockers. Addition of lecithin to bile salt micelles decreases the values of Kp of the beta-blockers. Mixed micelles incorporate hydrophobic compounds due to their large size and the fluidity of their core, but amphiphilic drugs, for which the interactions are predominantly polar/electrostatic, are poorly incorporated in mixed micelles of bile salts/lecithin.     

Partition coefficients of the iron (III) complexes of pyridoxal isonicotinoyl hydrazone and its analogs and the correlation to iron chelation efficacy. Correction of some reported partition coefficients

Pyridoxal isonicotinoyl hydrazone (PIH), salicylaldehydebenzoyl hydrazone (SBH), and their analogschelate iron(III) and show promise asorally effective drugs for treating diseases of iron overload. Theirbiological activity isrelated to their lipophilicity, as measured by their partition coefficients P betweenn-octanoland water. However, the method of calculating log P described in an article in this journal(Edwardet al. 1995; BioMetals, 8, 209-217) is faulty for compounds such as PIH, SBH andtheir analogs whichcontain adjacent hydrophilic groups. Consequently, the calculations reportedin the article, based on erro-neouslog P values of the chelating molecules, giveerroneous log P values of the iron(III) complexes. Thechelators most effective inmobilizing 59 Fe from reticulocytes have log P < 2.8, not log P < 0 and theiron(III)complexes of the most effective chelators have log P < 3.1, not log P < 0.     

The Metabolites of the Herbicide L-Phosphinothricin (Glufosinate) (Identification,Stability, and Mobility in Transgenic,Herbicide-Resistant,and Untransformed Plants)

The metabolism of the herbicide L-phosphinothricin (L-Pt) was analyzed in tobacco (Nicotiana tabacum), alfalfa (Medicago sativa), and carrot (Daucus carota). In transgenic, Pt-resistant plants expressing the Pt-N-acetyltransferase gene (pat), L-Pt was acetylated, resulting in two forms of N-acetyl-Pt (ac-Pt). In transgenic plants expressing only low pat-encoded acetylating activity as well as in genetically unmodified plants, three metabolic compounds 4-methylphosphinico-2-oxo-butanoic acid, 3-methylphosphinico-propanoic acid (MPP), and 4-methylphosphinico-2-hydroxy-butanoic acid (MHB) were identified. Hence, the transgene-encoded acetylation of L-Pt competes with a plant-specific degradation. The compounds MPP, MHB, and ac-Pt were found to be the final, stable products of the plant's metabolic pathways. The mobility of these stable compounds in the plant was investigated: L-Pt as well as the derived metabolites were found to be preferentially transported to the upper regions of the plant.     

Evaluation of lipophilicity and antitumour activity of parallel carboxamide libraries

Searching for molecules possessing antitumour activity, a parallel molecule library of aromatic carboxamides has been designed and synthesised. This work resulted in a "thiophene" sub-library containing a thiophene core and of a "furoyl" sub-library with a furoyl core, respectively. In both sub-libraries substitutions were carried out with six different groups resulting in six pairs of compounds differing in only the heteroatom of aromatic ring of the cores. To study the importance of the type of cores and the specific substitutions in relation to their lipophilicity and antitumour activity, lipophilicity of carboxamides was determined by chromatographical data (log k') and by software calculated parameters (CLOGP). Pairs of compounds were tested for their ability to inhibit the proliferation of the A431 cells by MTT assay. The isosteric molecule pairs were successfully separated. Our results showed that the experimentally determined (log k') and the calculated (CLOGP) lipophilicity parameters correlated well with each other. Furthermore, lipophilicity values of the thiophene sub-library were always higher than those in the furoyl sub-library. Moreover, compounds of the thiophene sub-library were more active than their respective furoyl pairs in our MTT antiproliferative assay. From these observations we can conclude that the higher the lipophilicity values the higher the antitumour activity of the carboxamides synthesised. Therefore, determination of lipophilicity by measuring the log k' or by calculating the CLOGP values of the carboxamide sub-libraries may help to predict their biological activities.     

GEPASI: a software package for modelling the dynamics, steady states and control of biochemical and other systems

GEPASI is a software system for modelling chemical and biochemicalreaction networks on computers running Microsoft Windows. Forany system of up to 45 metabolites and 45 reactions, each withany user-defined or one of 35 predefined rate equations, onecan produce trajectories of the metabolite concentrations andobtain a steady state (if it does exist). When steady-statesolutions are produced, elasticity and control coefficients,as defined in metabolic control analysis, are calculated. GEPASIalso allows the automatic generation of a sequence of simulationswith different combinations of parameter values, effectivelyscanning a hyper-solid in parameter space. Together with theability to produce user-defined columnar data files, these featuresallow for both very quick and systematic study of biochemicalpathway models. The source code (in C) is available on requestfrom the author, and while the user interface is dependent onhaving MS-Windows as the operating system, the numerical partis portable to other operating systems. GEPASI is suitable bothfor research and educational purposes. Although GEPASI was writtenwith biochemical pathways in mind, it can equally be used tosimulate other dynamical systems.     

Impact of the carbon chain length of novel platinum complexes derived from N-alkyl-propanediamines on their cytotoxic activity and cellular uptake

This work describes the synthesis and characterization of four new ligands derived from 1,3-propanediamine in addition to the preparation and characterization of their respective platinum(II) complexes by reaction with K₂PtCl₄. These ligands were obtained by the reaction of the corresponding alkyl mesylate with 1,3-propanediamine. We have prepared compounds having different carbon chains lengths in an attempt to correlate this factor, which influences the lipophilicity of the compounds, with cytotoxic activity. Octanol/water partition coefficients, the effect of the four complexes on the growth of two tumoral cell lines, and their cellular uptake were investigated. Increasing lipophilicity enhances the rate of cellular uptake and, consequently, the cytotoxic activity.     

Seed Dormancy in Red Rice : VII. Structure-Activity Studies of Germination Stimulants

Many chemically dissimilar substances break dormancy of seeds, but the relationship between chemical structure and physiological activity is unknown. In this study, the concentrations of organic acids, esters, aldehydes, alcohols, and inorganic weak acids required to elicit 50% germination of initially dormant, dehulled red rice seeds (Oryza sativa) were determined. The activity of most substances was very highly and inversely correlated to lipophilicity as measured by octanol/water partition coefficients; chemicals with the highest partition coefficients required the lowest concentrations to elicit the germination response. Relative efficacy was also dependent upon the functional group; generally, monocarboxylic acids were more effective than aldehydes, esters, hydroxyacids, and alcohols. Relative hydrophobicity plots supported a modulating role of the functional group. Dormancy-breaking activity of methyl formate, formic acid, nitrite, azide, and cyanide was higher than predicted based on lipophilicity and apparently was related to molecular size; compounds with smaller molecular widths were required at lower concentrations to achieve the 50% germination response.     

Potential drug targets on insomnia and intervention effects of Jujuboside A through metabolic pathway analysis as revealed by UPLC/ESI-SYNAPT-HDMS coupled with pattern recognition approach

Potential metabolites from the metabolic pathways could be therapeutic targets and useful for the discovery of broad spectrum drugs. UPLC/ESI-SYNAPT-HDMS coupled with pattern recognition methods including PCA, PLS-DA, OPLS-DA and Heatmap were integrated to examine the global metabolic signature of insomnia and intervention effects of Jujuboside A (JuA). Six unique pathways of the insomnia were identified using Ingenuity Pathway Analysis (IPA) software. The VIP-value threshold cutoff of the metabolites was set to 10, above this threshold, were filtered out as potential target biomarkers. Sixteen distinct metabolites were identified from these pathways, and 6 of them can be considered for rational drug design. It was further experimental validation that the changes in metabolic profiling were restored to their baseline values after JuA treatment according to the multivariate data analysis. Potential metabolite network of the insomnia was preliminarily predicted JuA-target interaction networks, and could be further explored for in silico docking studies with suitable drugs. Thus, our method is an efficient procedure for drug target identification through metabolic analysis. It can guide testable predictions, provide insights into drug action mechanisms and enable us to increase research productivity toward metabolomic drug discovery.     

Bile acid structure-activity relationship: evaluation of bile acid lipophilicity using 1-octanol/water partition coefficient and reverse phase HPLC

Two independent methods have been developed and compared to determine the lipophilicity of a representative series of naturally occurring bile acids (BA) in relation to their structure. The BA included cholic acid (CA), chenodeoxycholic acid (CDCA), ursodeoxycholic acid (UDCA), deoxycholic acid (DCA), hyodeoxycholic acid (HDCA), ursocholic acid (UCA), hyocholic acid (HCA), as well as their glycine and taurine amidates. Lipophilicity was determined using a 1-octanol/water shake-flask procedure and the experiments were performed at different pH and ionic strengths and at initial BA concentrations below their critical micellar concentrations (CMC) and the water solubility of the protonated form. The experimental data show that both the protonated (HA) and ionized (A-) forms of BA can distribute in 1-octanol, and consequently a partition coefficient for HA (logP' HA) and for A- (logP' A-) must be defined. An equation to predict a weighted apparent distribution coefficient (D) value as a function of pH and pKa has been developed and fits well with the experimental data. Differences between logP for protonated and ionized species for unconjugated BA were in the order of 1 log unit, which increased to 2 for glycine-amidate BA. The partition coefficient of the A- form increased with Na+ concentration and total ionic strength, suggesting an ion-pair mechanism for its partition into 1-octanol. Lipophilicity was also assessed using reverse phase chromatography (C-18-HPLC), and a capacity factor (K') for ionized species was determined. Despite a broad correlation with the logP data, some BA behaved differently. The logP values showed that the order of lipophilicity was DCA greater than CDCA greater than UDCA greater than HDCA greater than HCA greater than CA greater than UCA for both the protonated and ionized unconjugated and glycine-amidate BA, while the K' data showed an inversion for some BA, i.e., DCA greater than CDCA greater than CA greater than HCA greater than UDCA greater than HDCA greater than UCA. The logP data fitted well with other indirect measurements of BA monomeric lipophilicity such as albumin binding or accessible total hydrophobic surface area data calculated by energy minimization and molecular computer graphics. Differences between unconjugated and amidated BA are consistent with the presence of an amide bond and a lower pKa when pH dependence was studied. Capacity factors, on the other hand, were related to properties of BA micelles such as cholesterol-solubilizing capacity and membrane disruption, reflecting the BA detergency.(ABSTRACT TRUNCATED AT 400 WORDS)     

Relationship between lipophilicity and antitumor activity of molecule library of Mannich ketones determined by high-performance liquid chromatography, clogP calculation and cytotoxicity test

A series of Mannich ketones were synthesized in order to study the relative importance of structure and specific substitutions in relation to their lipophilicity and antitumor activity. Substitutions were carried out with morpholinyl, pirrolidinyl, piperidyl and tetrahydro-isoquinolyl groups in various positions on three different skeletons. Lipophilicity of Mannich ketones was characterised by chromatography data (log k') and by software calculated parameters (clogP). Compounds were tested on their ability to inhibit the proliferation of the A431 human adenocarcinoma cell line evaluated by MTT and apoptosis assays. The results suggest that the higher the lipophilicity values (log k' and clogP), the higher the antitumor and apoptotic activity of Mannich ketones. Determination of lipophilicity by measuring the log k' or by calculating the clogP values of the compounds may help to predict their biological activities.     

Mesenteric Organ Production, Hepatic Metabolism, and Renal Elimination of Norepinephrine and Its Metabolites in Humans

Abstract: This study used regional differences in plasma concentrations of norepinephrine and its metabolites to examine how production of the transmitter by sympathetic nerves, in particular, those innervating mesenteric organs, is integrated with metabolism by the liver and elimination by the kidneys. Higher concentrations of norepinephrine, its glycol metabolites 3,4-dihydroxyphenylglycol and 3-methoxy-4-hydroxyphenylglycol and their sulfate conjugates in portal venous than arterial plasma indicate substantial production of norepinephrine by mesenteric organs (15.5 nmol/min). Much lower concentrations of norepinephrine and its glycol metabolites in plasma leaving than entering the liver indicate their efficient hepatic removal (20 nmol/min). Higher concentrations of vanillylmandelic acid in the hepatic outflow than inflow indicate that this metabolic end product is produced largely from the norepinephrine and glycol metabolites removed by the liver. Renal elimination of vanillylmandelic acid (18–20 nmol/min), produced mainly by the liver (17 nmol/min), and of 3-methoxy-4-hydroxyphenylglycol sulfate (7–9 nmol/min), produced largely by mesenteric organs (7 nmol/min), comprised 86–91% of the total renal elimination of norepinephrine metabolites. The results show that mesenteric organs produce about one-half of the norepinephrine formed in the body. The liver removes substantial amounts of circulating norepinephrine and its glycol metabolites and converts these compounds to vanillylmandelic acid, which is then eliminated from the body by the kidneys. The sulfate conjugates are also metabolic end products eliminated by the kidneys. However, these metabolites are produced by extrahepatic tissues, in particular, mesenteric organs, which represent a significant source of sulfate-conjugated norepinephrine and 3,4-dihydroxyphenylglycol, and the main source of sulfate-conjugated 3-methoxy-4-hydroxyphenylglycol.     

Control in channelled pathways. A matrix method calculating the enzyme control coefficients

The usual equations expressing the enzyme control coefficients (quantitative indicators of 'global' control properties of a pathway) via the elasticity coefficients (reflecting local kinetic properties of an enzyme reaction), cannot be applied to a variety of 'non-ideal' pathways, in particular to pathways with metabolic channelling. Here we show that the relationship between the control and elasticity coefficients can be obtained by considering such a metabolic pathway as a network of elemental chemical conversions (steps). To calculate the control coefficients of enzymes one should first determine the elasticity coefficients of such elemental steps and then take their appropriate combinations. Although the method is illustrated for a channelled pathway it can be used for any non-ideal pathway including those with high enzyme concentrations where the sequestration of metabolites by enzymes cannot be neglected.     

Metabolism studies of the Kratom alkaloids mitraciliatine and isopaynantheine, diastereomers of the main alkaloids mitragynine and paynantheine, in rat and human urine using liquid chromatography-linear ion trap-mass spectrometry

Mitragyna speciosa (Kratom in Thai), native in Southeast Asia, is increasingly misused as a herbal drug of abuse. During metabolism studies on the Kratom alkaloids mitragynine, its diastereomers speciogynine and speciociliatine as well as paynantheine in rats and humans, further isomeric compounds were detected in Kratom users' urine. The question arose whether these compounds were formed from the low abundant, isomeric alkaloids mitraciliatine (MC) and isopaynantheine (ISO-PAY). Therefore, the aim of the presented study was to identify using liquid chromatography-linear ion trap-mass spectrometry their phase I and II metabolites in rat urine after administration of pure MC or ISO-PAY, to confirm their formation in humans, and finally to confirm whether the above-mentioned isomeric compounds in human urine represent MC and ISO-PAY and/or their metabolites. The metabolic pathways of both alkaloids in rats were found to be comparable to those of their corresponding diastereomers. In the human urines tested, not all metabolites found in rats could be detected because of the much lower amounts of MC and ISO-PAY in Kratom. However, all the above-mentioned so far unknown isomeric compounds could be identified in the human urine samples as MC, ISO-PAY and/or their metabolites. The used LC separation was also suitable for the differentiation of all other Kratom alkaloids and their metabolites in human urine.     

The role of lipophilicity in determining binding affinity and functional activity for 5-HT2A receptor ligands

The lipophilicity of a set of 5-HT(2A) ligands was determined using immobilized-artificial-membrane chromatography, a method that generates values well correlated with octanol-water partition coefficients. For agonists, a highly significant linear correlation was observed between binding affinity and lipophilicity. For ligands exhibiting partial agonist or antagonist properties, the lipophilicity was consistently higher than would be expected for an agonist of comparable affinity. The results suggest a possible method for distinguishing agonists from antagonists in high-throughput screening when a direct assay for functional activity is either unavailable or impractical.