Développement de modèles prédictifs de la toxicocinétique de substances organiques

Les modèles pharmacocinétiques à base physiologique (PBPK) permettent de simuler la dose interne de substances chimiques sur la base de paramètres spécifiques à l’espèce et à la substance. Les modèles de relation quantitative structure-propriété (QSPR) existants permettent d’estimer les paramètres spécifiques au produit (coefficients de partage (PC) et constantes de métabolisme) mais leur domaine d’application est limité par leur manque de considération de la variabilité de leurs paramètres d’entrée ainsi que par leur domaine d’application restreint (c. à d., substances contenant CH3, CH2, CH, C, C=C, H, Cl, F, Br, cycle benzénique et H sur le cycle benzénique). L’objectif de cette étude est de développer de nouvelles connaissances et des outils afin d’élargir le domaine d’application des modèles QSPR-PBPK pour prédire la toxicocinétique de substances organiques inhalées chez l’humain. D’abord, un algorithme mécaniste unifié a été développé à partir de modèles existants pour prédire les PC de 142 médicaments et polluants environnementaux aux niveaux macro (tissu et sang) et micro (cellule et fluides biologiques) à partir de la composition du tissu et du sang et de propriétés physicochimiques. L’algorithme résultant a été appliqué pour prédire les PC tissu:sang, tissu:plasma et tissu:air du muscle (n = 174), du foie (n = 139) et du tissu adipeux (n = 141) du rat pour des médicaments acides, basiques et neutres ainsi que pour des cétones, esters d’acétate, éthers, alcools, hydrocarbures aliphatiques et aromatiques. Un modèle de relation quantitative propriété-propriété (QPPR) a été développé pour la clairance intrinsèque (CLint) in vivo (calculée comme le ratio du Vmax (μmol/h/kg poids de rat) sur le Km (μM)), de substrats du CYP2E1 (n = 26) en fonction du PC n octanol:eau, du PC sang:eau et du potentiel d’ionisation). Les prédictions du QPPR, représentées par les limites inférieures et supérieures de l’intervalle de confiance à 95% à la moyenne, furent ensuite intégrées dans un modèle PBPK humain. Subséquemment, l’algorithme de PC et le QPPR pour la CLint furent intégrés avec des modèles QSPR pour les PC hémoglobine:eau et huile:air pour simuler la pharmacocinétique et la dosimétrie cellulaire d’inhalation de composés organiques volatiles (COV) (benzène, 1,2-dichloroéthane, dichlorométhane, m-xylène, toluène, styrène, 1,1,1 trichloroéthane et 1,2,4 trimethylbenzène) avec un modèle PBPK chez le rat. Finalement, la variabilité de paramètres de composition des tissus et du sang de l’algorithme pour les PC tissu:air chez le rat et sang:air chez l’humain a été caractérisée par des simulations Monte Carlo par chaîne de Markov (MCMC). Les distributions résultantes ont été utilisées pour conduire des simulations Monte Carlo pour prédire des PC tissu:sang et sang:air. Les distributions de PC, avec celles des paramètres physiologiques et du contenu en cytochrome P450 CYP2E1, ont été incorporées dans un modèle PBPK pour caractériser la variabilité de la toxicocinétique sanguine de quatre COV (benzène, chloroforme, styrène et trichloroéthylène) par simulation Monte Carlo. Globalement, les approches quantitatives mises en œuvre pour les PC et la CLint dans cette étude ont permis l’utilisation de descripteurs moléculaires génériques plutôt que de fragments moléculaires spécifiques pour prédire la pharmacocinétique de substances organiques chez l’humain. La présente étude a, pour la première fois, caractérisé la variabilité des paramètres biologiques des algorithmes de PC pour étendre l’aptitude des modèles PBPK à prédire les distributions, pour la population, de doses internes de substances organiques avant de faire des tests chez l’animal ou l’humain.

Quantitative structure-toxicity relationships for chlorophenols to bioluminescent lux-marked bacteria using atom-based semi-empirical molecular-orbital descriptors

Literature data on the toxicity of chlorophenols for three luminescent bacteria (Vibrio fischeri, and the lux-marked Pseudomonas fluorescens 10586s pUCD607 and Burkholderia spp. RASC c2 (Tn4431)) have been analyzed in relation to a set of computed molecular physico-chemical properties. The quantitative structure-toxicity relationships of the compounds in each species showed marked differences when based upon semi-empirical molecular-orbital molecular and atom based properties. For mono-, di- and tri-chlorophenols multiple linear regression analysis of V. fischeri toxicity showed a good correlation with the solvent accessible surface area and the charge on the oxygen atom. This correlation successfully predicted the toxicity of the heavily chlorinated phenols, suggesting in V. fischeri only one overall mechanism is present for all chlorophenols. Good correlations were also found for RASC c2 with molecular properties, such as the surface area and the nucleophilic super-delocalizability of the oxygen. In contrast the best QSTR for P. fluorescens contained the 2nd order connectivity index and ELUMO suggesting a different, more reactive mechanism. Cross-species correlations were examined, and between V. fischeri and RASC c2 the inclusion of the minimum value of the nucleophilic susceptibility on the ring carbons produced good results. Poorer correlations were found with P. fluorescens highlighting the relative similarity of V. fischeri and RASC c2, in contrast to that of P. fluorescens.

Quantitative structure-property relationships for octanol-air partition coefficients of polychlorinated naphthalenes, chlorobenzenes and p,p '-DDT

The octanol-air partition coefficient (K-OA) is a key descriptor of chemicals partitioning between the atmosphere and environmental organic phases. Quantitative structure-property relationships (QSPR) are necessary to model and predict KOA from molecular structures. Based on 12 quantum chemical descriptors computed by the PM3 Hamiltonian, using partial least squares (PLS) analysis, a QSPR model for logarithms of K-OA to base 10 (log K-OA) for polychlorinated naphthalenes (PCNs), chlorobenzenes and p,p'-DDT was obtained. The cross-validated Q(cum)(2) value of the model is 0.973, indicating a good predictive ability of the model. The main factors governing log K-OA of the PCNs, chlorobenzenes, and p,p'-DDT are, in order of decreasing importance, molecular size and molecular ability of donating/accepting electrons to participate in intermolecular interactions. The intermolecular dispersive interactions play a leading role in governing log K-OA. The more chlorines in PCN and chlorobenzene molecules, the greater the log K-OA values. Increasing E-LUMO (the energy of the lowest unoccupied molecular orbital) of the molecules leads to decreasing log K-OA values, implying possible intermolecular interactions between the molecules under study and octanol molecules. (C) 2002 Elsevier Science Ltd. All rights reserved.

Quantitative structure-property relationships for octanol-air partition coefficients of polychlorinated biphenyls

Based on nine quantum chemical descriptors computed by PM3 Hamiltonian, using partial least squares analysis, a significant quantitative structure-property relationship for the logarithm of octanol-air partition coefficients (log K-OA) of polychlorinated biphenyls (PCBs) was obtained. The cross-validated Q(cum)(2) value of the model is 0.962, indicating a good predictive ability. The intermolecular dispersive interactions and thus the size of the PCB molecules play a key role in governing log K-OA. The greater the size of PCB molecules, the greater the log K-OA values. Increasing E-LUMO (the energy of the lowest unoccupied molecular orbital) values of the PCBs leads to decreasing log K-OA values, indicating possible interactions between PCB and octanol molecules. Increasing Q(Cl)(+) (the most positive net atomic charges on a chlorine atom) and Q(C)(-) (the largest negative net atomic charge on a carbon atom) values of PCBs results in decreasing log K-OA values, implying possible intermolecular electrostatic interactions between octanol and PCB molecules. (C) 2002 Elsevier Science Ltd. All rights reserved.

Quantitative structure-property relationships (QSPRs) on direct photolysis of PCDDs

By the use of partial least squares (PLS) method and 27 quantum chemical descriptors computed by PM3 Hamiltonian, a statistically significant QSPR were obtained for direct photolysis quantum yields (Y) of selected Polychlorinated dibenzo-p-dioxins (PCDDs). The QSPR can be used for prediction. The direct photolysis quantum yields of the PCDDs are dependent on the number of chlorine atoms bonded with the parent structures, the character of the carbon-oxygen bonds, and molecular polarity. Increasing bulkness and polarity of PCDDs lead to decrease of log Y values. Increasing the frontier molecular orbital energies (E-lumo and E-homo) and heat of formation (HOF) values leads to increase of log Y values. (C) 2001 Elsevier Science Ltd. All rights reserved.

Quantitative structure-property relationship study on reductive dehalogenation of selected halogenated aliphatic hydrocarbons in sediment slurries

In this study, by the use of partial least squares (PLS) method and 26 quantum chemical descriptors computed by PM3 Hamiltonian, a quantitative structure-property relationship (QSPR) model was developed for reductive dehalogenation rate constants of 13 halogenated aliphatic compounds in sediment slurry under anaerobic conditions. The model can be used to explain the dehalogenation mechanism. Halogenated aliphatic compounds with great energy of the lowest unoccupied molecular orbital (E-lumo), total energy (TE), electronic energy (EE), the smallest bond order of the carbon-halogen bonds (BO) and the most positive net atomic charges on an atom of the molecule (q(+)) values tend to be reductively dehalogenated slow, whereas halogenated aliphatic compounds with high values of molecular weight (Mw), average molecular polarizability (a) and core-core repulsion energy (CCR) values tend to be reductively dehalogenated fastest. (C) 2001 Published by Elsevier Science Ltd.

Quantitative structure-property relationship studies on n-octanol/water partitioning coefficients of PCDD/Fs

Based on some fundamental quantum chemical descriptors computed by PM3 Hamiltonian, by the use of partial least-squares (PLS) analysis, a significant quantitative structure-property relationship (QSPR) model for logK(ow) of polychlorinated dibenzo-p-dioxins and dibenzo-p-furans (PCDD/Fs) was obtained. The QSPR can be used for prediction. The intermolecular dispersive interactions and thus the bulkness of the PCDD/Fs are the main factors affecting the logK(ow). The more chlorines in the PCDD/F molecule, the greater the logK(ow) values. (C) 2001 Elsevier Science Ltd. All rights reserved.

Peripheries of molecular projection in three-dimensional space and studies of quantitative structure-activity relationships

A new index, i.e., the periphery representation of the projection of a molecule from 3D space to a 2D plane is described. The results, correlation with toxicity of substituted nitrobenzenes, obtained by using periphery descriptors are much better than that obtained by using the areas (i.e., shadows) of projections of the compounds. Even better results were achieved by using the combination of periphery descriptors and the projections areas as well as the indicated variable K reflecting the action of group NO position on the benzene ring.

Applications of the method of three-dimensional projection of a molecule in quantitative structure-activity relationship of phenol derivatives

Five variables for phenol derivatives were calculated by molecular projection in three-dimensional space which were combined with eight quantum-chemical parameters and three Am indices. These variables were selected by using leaps-and-bounds regression analysis. Multiple linear regression analysis and artificial neural networks' were performed, and the results obtained by using. artificial neural networks are superior than that obtained by using multiple linear regression.

Applications of generalized torsion angles in studies on quantitative structure-activity relationship

In this article, generalized torsion angles of derivatives of 1-[(2-hydroxyethoxy)methy1]-6(phenylthio)thymine(HEPT) were calculated, which include abundant three dimensional information of molecules. Molecular similarity matrix was built based on the calculated generalized torsion angles. These similarities were taken as the new variables, and the new variables were selected by using Leaps-and-Bounds regression analysis. Multiple regression analysis and neural networks were performed, and the satisfactory results were achieved by using the neural networks.

Effective descriptions of molecular structures and the quantitative structure - Activity relationship studies

In this research. we found CoMFA alone could not obtain sufficiently a strong equation to allow confident prediction for aminobenzenes. When some other parameter. such as heat of molecular formation of the compounds, was introduced into the CoMFA model, the results Were improved greatly. It gives us a hint that a better description for molecular structures will yield a better prediction model, and this hint challenged us to look for another method-the projection areas of molecules in 3D space for 3D-QSAR. It is surprising that much better results than that obtained by using CoMFA Were achieved. Besides the CoMFA analysis. multiregression analysis and neural network methods for building the models were used in this paper.