RESUMEN
Based on recently published initial experimental results on the intercalation of a class of broad spectrum antiparasitic compounds, we present a purely theoretical approach for determining if these compounds may preferentially intercalate with guanosine/cytosine (GC)-rich or adenosine/thymidine (TA)-rich regions of DNA. The predictive model presented herein is based upon utilization of density functional theory (DFT) to determine a priori how the best intercalator may energetically and sterically interact with each of the nucleoside base pairs. A potential new method using electrostatic potential maps (EPMs) to visually select the best poses is introduced and compared to the existing brute-force center of mass (COM) approach. The EPM and COM predictions are in agreement with each other, but the EPM method is potentially much more efficient. We report that 4-azatryptantrin, the best intercalator, is predicted to favor π-stacking with GC over that of TA by approximately 2-4 kcal/mol. This represents a significant difference if one takes into account the Boltzmann distribution at physiological temperature. This theoretical method will be utilized to guide future experimental studies on the elucidation of possible mechanism(s) for the action of these antiparasitic compounds at the molecular level.
Asunto(s)
ADN/química , ADN/metabolismo , Quinazolinas/química , Quinazolinas/metabolismo , Simulación por Computador , Modelos Químicos , Modelos Moleculares , Estructura MolecularRESUMEN
A new submolecular quantitative structure activity relationship (QSAR) descriptor was applied toward elucidating the anti-malarial pharmacophore of tryptanthrins, a class of compounds known for their anti-parasitic activities. The new descriptor is based on experimental and computational measurements of the tunneling barriers of individual lobes of the molecular orbitals. Lobe-by-lobe QSAR correlation plots revealed a single lobe of the LUMO to be strongly associated with tryptanthrin's anti-malarial activity. The correlation also showed a threshold behavior wherein barriers below a particular value show low IC50 values. Above the threshold, the correlation of IC50 vs the logarithm of the barrier is linear with R2â¯=â¯0.999. This barrier threshold may be applied as a design criterion for future tryptanthrin-based anti-malarial lead optimization. The new descriptor may be broadly applicable toward other molecular systems of interest, such as catalysts, pesticides, and herbicides. The authors have named the new descriptor: submolecular tunneling analysis of barriers (STAB).