Quantitative structure-activity relationships of mutagenic activity from quantum topological descriptors: triazenes and halogenated hydroxyfuranones (mutagen-X) derivatives

Effect of iontophoresis on dacarbazine cutaneous delivery for melanoma topical treatment

Dacarbazine (DTIC) is the drug of choice for melanoma treatment, but its systemic administration is related to several adverse effects. Here, DTIC topical delivery stimulated by iontophoresis is proposed to overcome such drawbacks. Hence, this work analyzed the impact of anodal iontophoresis on DTIC cutaneous delivery to provide an innovative topical alternative for melanoma treatment. The electrical stability of the drug was evaluated prior to the iontophoretic experiments, which demonstrated the need to add an antioxidant to the drug formulation. DTIC cutaneous permeation was evaluated in vitro for 6 h using three current densities (0.10, 0.25, and 0.50 mA/cm2). In addition, the effect of DTIC against skin cancer cells (MeWo and WM164) was investigated for 72 h of exposure to the drug. Iontophoresis stimulated skin drug permeation compared to the passive control. However, the antioxidant presence reduced DTIC permeation under the lower currents of 0.10 and 0.25 mA/cm2, which was compensated by increasing the current density to 0.50 mA/cm2. At 0.50 mA/cm2, iontophoresis enhanced topical cutaneous drug permeation 7-fold (p < 0.05) compared to the passive control. DTIC showed a concentration-dependent antiproliferative effect on melanoma cell lines. Thus, iontophoresis intensifies DTIC skin penetration in concentrations that can reduce cell viability and induce cell death. In conclusion, DTIC cutaneous delivery mediated by iontophoresis is a promising approach for treating melanomas and other skin tumors.

Chromatographic method for dacarbazine quantification in skin permeation experiments

Dacarbazine (DTIC) is a chemotherapeutic drug currently used for the systemic treatment of melanomas. Considering the easy access to these tumors, a topical route of drug administration could provide a more comfortable and less toxic treatment. However, DTIC quantification aiming at the design of topical formulations is challenging, pondering all the interferents present in the drug samples recovered from the skin. Hence, this work intended to validate a selective chromatographic method for DTIC determination in skin permeation studies. A reversed-phase C18 column was used as a stationary phase, and gradient elution of a mobile phase consisting of methanol and pH 6.5 sodium phosphate monohydrate buffer (0.01 mol/L) at a flow rate of 1.0 mL/min was implemented. DTIC was detected at 364 nm. The method was selective against skin interferents, linear (r = 0.9995) in a concentration range of 1.0-15.0 μg/mL, precise with an overall variation coefficient lower than 3.8%, accurate achieving recovery from the skin layers within 91-112%, and sensitive for the proposed application (detection limit = 0.10 μg/ mL, quantification limit = 0.30 μg/mL). Furthermore, the analytical method was successfully tested in in vitro skin permeation studies. In conclusion, the developed method is appropriate for DTIC analysis from the skin sample matrix.

Quantitative Structure-Activity/Property/Toxicity Relationships through Conceptual Density Functional Theory-Based Reactivity Descriptors

Developing effective structure-activity/property/toxicity relationships (QSAR/QSPR/QSTR) is very helpful in predicting biological activity, property, and toxicity of a given set of molecules. Regular change in these properties with the structural alteration is the main reason to obtain QSAR/QSPR/QSTR models. The advancement in making different QSAR/QSPR/QSTR models to describe activity, property, and toxicity of various groups of molecules is reviewed in this chapter. The successful implementation of Conceptual Density Functional Theory (CDFT)-based global as well as local reactivity descriptors in modeling effective QSAR/QSPR/QSTR is highlighted.

Design of novel leads: ligand based computational modeling studies on non-nucleoside reverse transcriptase inhibitors (NNRTIs) of HIV-1

Researchers are on the constant lookout for new antiviral agents for the treatment of AIDS. In the present work, ligand based modeling studies are performed on analogues of substituted phenyl-thio-thymines, which act as non-nucleoside reverse transcriptase inhibitors (NNRTIs) and novel leads are extracted. Using alignment-dependent descriptors, based on group center overlap (SALL, HDALL, HAALL and RALL), an alignment-independent descriptor (S log P), a topological descriptor (Balaban index (J)) and a 3D descriptor dipole moment (μ) and shape based descriptors (Kappa 2 index ((2)κ)), a correlation is derived with inhibitory activity. Linear and non-linear techniques have been used to achieve the goal. Support Vector Machine (SVM, R = 0.929, R(2) = 0.863) and Back Propagation Neural Network (BPNN, R = 0.928, R(2) = 0.861) methods yielded near similar results and outperformed Multiple Linear Regression (MLR, R = 0.915, R(2) = 0.837). The predictive ability of the models are cross-validated using a test dataset (SVM: R = 0.846, R(2) = 0.716, BPNN: R = 0.841, R(2) = 0.707 and MLR: R = 0.833, R(2) = 0.694). It is concluded that the hydrophobicity (S log P) and the polarity (μ) of a ligand and the presence of hydrogen donor (HDALL) moieties are the deciding factors in improving antiviral activity and pharmaco-therapeutic properties. Based on the above findings, a virtual dataset is created to extract probable leads with reasonable antiviral activity as well as better pharmacophoric properties.

Modeling biophysical and biological properties from the characteristics of the molecular electron density, electron localization and delocalization matrices, and the electrostatic potential

The electron density and the electrostatic potential are fundamentally related to the molecular hamiltonian, and hence are the ultimate source of all properties in the ground- and excited-states. The advantages of using molecular descriptors derived from these fundamental scalar fields, both accessible from theory and from experiment, in the formulation of quantitative structure-to-activity and structure-to-property relationships, collectively abbreviated as QSAR, are discussed. A few such descriptors encode for a wide variety of properties including, for example, electronic transition energies, pK(a)'s, rates of ester hydrolysis, NMR chemical shifts, DNA dimers binding energies, π-stacking energies, toxicological indices, cytotoxicities, hepatotoxicities, carcinogenicities, partial molar volumes, partition coefficients (log P), hydrogen bond donor capacities, enzyme-substrate complementarities, bioisosterism, and regularities in the genetic code. Electronic fingerprinting from the topological analysis of the electron density is shown to be comparable and possibly superior to Hammett constants and can be used in conjunction with traditional bulk and liposolubility descriptors to accurately predict biological activities. A new class of descriptors obtained from the quantum theory of atoms in molecules' (QTAIM) localization and delocalization indices and bond properties, cast in matrix format, is shown to quantify transferability and molecular similarity meaningfully. Properties such as "interacting quantum atoms (IQA)" energies which are expressible into an interaction matrix of two body terms (and diagonal one body "self" terms, as IQA energies) can be used in the same manner. The proposed QSAR-type studies based on similarity distances derived from such matrix representatives of molecular structure necessitate extensive investigation before their utility is unequivocally established.

Global and local chemical reactivities of mutagen X and simple derivatives

Registered by the World Health Organization (WHO), 3-chloro-4-(dichloromethyl)-5-hydroxy-2(5H)-furanone (MX) is one of the strongest bacterial mutagens ever tested, as highlighted by the Ames Salmonella typhimurium TA100 assay. We provide new insights concerning this mutagenic activity on the basis of global and local theoretically defined electrophilicity indices. Our results further support the idea that mutagenicity of MX and its analogues is related more closely to one-electron transfer processes from the electron-rich biological environment than to adduct formation processes. We also stress that, although the Z-open tautomers are intrinsically more electrophilic than furanone ring analogues, the observed mutagenic activity is significantly correlated only to the electrophilicity response of the ring forms. In that context, we also emphasize that it is electrophilicity at the C α in the α-β unsaturated carbonyl moiety that exhibits a strong correlation with the observed mutagenic activity.

Prediction of interaction energies of substituted hydrogen-bonded Watson-Crick cytosine:guanine(8X) base pairs

We investigated the variation in the interaction energy between the Watson-Crick hydrogen-bonded DNA base pairs guanine and cytosine (G(8X):C), where guanine is substituted in the C8 position by 37 different functional groups. Base pairs were optimized at the B3LYP/6-311+G(2d,p) level. A base pair complex containing a more strongly electron-withdrawing group remarkably forms a more stable base pair with C. Multivariate linear regression provided a quantitative relationship between the interaction energies and descriptors generated by the quantum chemical topology (QCT) approach. The descriptors were sampled from the monomers only, not the supermolecular base pair complexes. A model with r2 = 0.96 and a root-mean-square (rms) value of 0.6 kJ/mol was obtained for a training set of 28 base pair complexes. The model was tested by an external test set of 9 complexes, yielding r2 = 0.99 and an rms value of 0.2 kJ/mol. The results indicated that the bonds C6=O6 and N2-H2 at the hydrogen-bonded frontier of the guanine derivatives play an important role in transmitting the substituent effects. A linear correlation between substitution energies and Hammett constants (sigma(m)) was also obtained for all 37 substituents, yielding r2 = 0.82 and an rms value of 1.2 kJ/mol. The model based on QCT descriptors can therefore be used for the prediction of the interaction energy of the base pair G(8x):C, strictly based on data for the G(8x) monomers only.

Computation of relative bond dissociation enthalpies (DeltaBDE) of phenolic antioxidants from quantum topological molecular similarity (QTMS)

A recently proposed method called quantitative topological molecular similarity (QTMS) generated a model for the computation of the relative substituent effects on the bond dissociation enthalpies (DeltaBDEs) for a set of 39 phenols. The data set includes a diverse set of substituents with monosubstituted and poly-substituted derivatives that exhibit different electronic and steric effects. Many share common structural features with already well-established antioxidants. QTMS reveals the active region of the substituted phenols and identifies the electronic descriptors that best explain the range of DeltaBDEs observed. For substituents in the 4-X position (para) we find that our model requires a correction for radical stabilization enthalpy (RSE). Application of the QTMS methodology yields an unrivalled QSAR with r(2) = 0.98 and q(2) = 0.85 for the bond dissociation enthalpies of this phenolic antioxidant data set.

Validation of critical points in the electron density as descriptors by building quantitative structure-property relationships for the atomic polar tensor

A crucial component of research in the field of quantitative structure-activity/property relationships is the identification of molecular descriptors relevant to the activity or property of interest. Descriptors based on the topology of the electron density as formulated in Bader's theory of atoms in molecules are investigated in detail in this work. In a model study, the authors investigate their ability to predict the atomic polar tensor (the gradient of the molecular dipole moment), which contains information on the vibrational intensities in infrared spectroscopy and constitutes a scheme for partitioning the total charge distribution into atomic charges. The atomic polar tensor may therefore be used to investigate whether the descriptors give adequate information on the local electronic structure in the molecule. Both the trace of the atomic polar tensor and for planar molecules its out-of-plane component may be interpreted as definitions of atomic charges suitable for prediction. Hydrogen and carbon atoms in a set of 60 aromatic compounds with various substituents have been studied. Excellent results for prediction of hydrogen and carbon charges have been achieved with cross-validated squared correlation coefficients between predicted and theoretical values varying from 0.92 and 0.977 for the most complex set of substituents when the value, Laplacian, and ellipticity of the electron density in the bond critical points are used as descriptors. The carbon charges defined from the trace of the atomic polar tensor are correlated with its out-of-plane component whereas such relationship is not observed for the hydrogen charges studied in this work.

QSAR models based on quantum topological molecular similarity

A new method called quantum topological molecular similarity (QTMS) was fairly recently proposed [J. Chem. Inf. Comp. Sc., 41, 2001, 764] to construct a variety of medicinal, ecological and physical organic QSAR/QSPRs. QTMS method uses quantum chemical topology (QCT) to define electronic descriptors drawn from modern ab initio wave functions of geometry-optimised molecules. It was shown that the current abundance of computing power can be utilised to inject realistic descriptors into QSAR/QSPRs. In this article we study seven datasets of medicinal interest : the dissociation constants (pK(a)) for a set of substituted imidazolines , the pK(a) of imidazoles , the ability of a set of indole derivatives to displace [(3)H] flunitrazepam from binding to bovine cortical membranes , the influenza inhibition constants for a set of benzimidazoles , the interaction constants for a set of amides and the enzyme liver alcohol dehydrogenase , the natriuretic activity of sulphonamide carbonic anhydrase inhibitors and the toxicity of a series of benzyl alcohols. A partial least square analysis in conjunction with a genetic algorithm delivered excellent models. They are also able to highlight the active site, of the ligand or the molecule whose structure determines the activity. The advantages and limitations of QTMS are discussed.

The cytotoxicity of ortho alkyl substituted 4-X-phenols: A QSAR based on theoretical bond lengths and electron densities

A new method called quantum topological molecular similarity (QTMS) was recently proposed [O'Brien, S. E.; Popelier, P. L. A. J. Chem. Inf. Comp. Sci.2001, 41, 764] and has been shown to be successful in a variety of medicinal, ecological and physical organic QSAR/QSPRs. QTMS method uses electronic descriptors drawn from ab initio wavefunctions of geometry-optimized molecules. We investigated a remarkable and unusual set of ortho alkyl-substituted phenols [Selassie, C. D.; Verma, R. P.; Kapur, S.; Shusterman, A. J.; Hansch, C. J. Chem. Soc., Perkin2002, II, 1112], recently studied by the Hansch group. Our results do not support their proposal that a steric factor is important in the determination of the cytotoxicity of this set of substituted phenols. Thus, we conclude that the cytotoxicity of these sterically encumbered phenols is dependent primarily on electronic and radical effects, and that steric issues do not appear to be a critical distinguishing factor.

Quantum chemical topology (QCT) descriptors as substitutes for appropriate Hammett constants

A technique called quantum topological molecular similarity (QTMS) was recently proposed [J. Chem. Inf. Comput. Sci., 2001, 41, 764] in order to construct a variety of medicinal, ecological and physical organic QSAR/QSPRs, based on modern ab initio wave functions of geometry optimised molecules, in combination with quantum chemical topology (QCT). The current abundance of computing power can be utilised to inject realistic descriptors into QSAR/QSPRs. In previous work [J. Chem. Soc., Perkin Trans. 2, 2002, 1231] it was proven that a set of Hammett constants (sigma(p), sigma(m), sigma(I) and sigma(p)0) for a sizeable set of mono- and polysubstituted carboxylic acids can be replaced by QCT bond descriptors. Using QTMS and proper statistical validation we examined seven data sets in total. The first three sets (para-substituted phenols (sigma-), substituted toluenes (sigma+) and bromophenethylamines (sigma+)) corroborate that a wider class of Hammett constants can also be replaced by QCT descriptors. A fourth set (benzyl radicals) focuses on non-Hammett behaviour being superimposed on Hammett behaviour. QCT descriptors selectively correlate with Hammett behaviour. The QTMS analysis of the last three sets (toxicity of benzyl alcohols, chromatographic capacity factors of chalcones and herbicidal activity of 5-chloro-2,3-dicyanopyrazines) screens for false positives. This test is successfully passed in that QCT descriptors fail when lipophilicity/hydrophobicity is in charge. Hence, overall, the discriminatory capacity of QCT descriptors is established, in detecting Hammett behaviour and specifically replacing the Hammett constants by more modern and non-empirical descriptors.