A density functional study of flavonoid compounds with anti-HIV activity

Deep eutectic solvent assisted synthesis of dihydropyrimidinones/thiones via Biginelli reaction: theoretical investigations on their electronic and global reactivity descriptors

DES 2 (ZrOCl2·8H2O-ethylene glycol at 1 : 2 ratio) was used as a catalyst for the synthesis of dihydropyrimidinones/thiones via a Biginelli reaction.

Metal-free synthesis of 3-chalcogenyl chromones from alkynyl aryl ketones and diorganyl diselenides/disulfides mediated by PIFA

3-Selenyl/sulfenyl chromones/thiochromones were conveniently synthesized from the PIFA-mediated reactions between alkynyl aryl ketones bearing an ortho-methoxy/methylthio group and diorganyl diselenides/disulfides.

Structure characteristics of flavonoids for cyclooxygenase-2 mRNA inhibition in lipopolysaccharide-induced inflammatory macrophages

Flavonoids are natural active components. They distribute widely in edible plants with high activity of anti-inflammation. Inhibition of cyclooxygenase-2(COX-2) was determined by real-time fluorescent quantitative olymerase chain reaction (RTFQ-PCR). And structure characteristics of flavonoids for COX-2 inhibition were mainly analyzed by a quantitative structure activity relationship (QSAR) model. Descriptors such as SMR_VSA5, vsurf_DD12, reactive were the top three important independent variables to COX-2 mRNA inhibiton in RAW264.7. Low SMR_VSA5 value meant a lower molecular refractivity resulting in a lower COX-2 mRNA inhibition. High vsurf_DD12 value related to poor molecular balance and showed profound adverse to COX-2 mRNA inhibition. Reactive group in this paper referred to C2-C3 double bond contributed negatively to COX-2 mRNA inhibition. Glycosidic and C3-OH substitutions may lower SMR_VSA5 value. It indicated that flavanones such as hesperetin, naringenin, liquiritigenin were efficient to repress COX-2 mRNA and they were potential anti-inflammatory natural products. Further, substitution with a glucopyranosyl at C-6 resulted in a poorer molecule balance than that at C-8 and a lower COX-2 mRNA inhibiton accordingly. This may expain why orientin and vitexin exhibited better anti-inflammatory activity than their isomers homoorientin and isovitexin. Also, methoxyl groups at C-4' may also be a favorable flavonoid structural characteristic for COX-2 mRNA inhibiton. These results provide valuable information on understanding the high anti-inflammatory activity of flavonoids.

Estrogenic Active Stilbene Derivatives as Anti-Cancer Agents: A DFT and QSAR Study

Exploring different quantum chemical quantities for lead compounds is an ongoing approach in identifying crucial structural activity related features that are contributing into their biological activities. Herein, activity-related quantum chemical calculations were performed for the selected estrogenic stilbene derivatives using density functional theory (DFT) with B3LYP functional and 6-311++G** basis set. In addition, specific activity-related geometry-independent drug-like properties are discussed for these derivatives. To obtain the mathematical model that correlates the chemical descriptors with their measured estrogenic activities, the quantitative structure activity relationship (QSAR) is established using multiple linear regression (MLR) and support vector regression (SVR) methods. Satisfactory fit with a reasonable regression correlation coefficient (${\rm{R}}^{2}= 0.78$R2=0.78) between predicted and experimental $pEC_{50}$pEC50 values is observed using MLR method. The present study identifies the essential physicochemical descriptors that effectively contribute in the estrogenic activity. The applied approach provides helpful insight into the designing novel estrogenic agents with improved anticancer activities.

Theoretical insights into the competitive metal bioaffinity of lactoferrin as a metal ion carrier: a DFT study

Metal–protein complexes, specifically lactoferrin (Lf), an iron-binding glycoprotein found naturally in milk and several other body fluids play a pivotal role in all living organisms.

Study of Structure and Permeability Relationship of Flavonoids in Caco-2 Cells

Flavonoids exhibit a broad range of biological activities. However, poor absorption of some flavonoids is a major limitation for use of flavonoids as nutraceuticals. To investigate the structure requirements for flavonoids intestinal absorption, transepithelial transport and cellular accumulation (CA) of 30 flavonoids were determined using the Caco-2 cell monolayer. The bilateral permeation of five types of flavonoids followed the order: flavanones ≥ isoflavones > flavones ≥ chalcones > flavonols. The concentration of flavonoids accumulated in cells did not correlate with cell penetration since the correlation coefficient between the apparent permeability coefficient (P_app) and their corresponding CA was poor (R² < 0.3). Most flavonoids exhibited a ratio of 0.8–1.5 for P_{app A to B}/P_{app B to A}, suggesting passive diffusion pathways. However, luteolin, morin and taxifolin may involve the efflux mechanisms. The quantitative structure-permeability relationship (QSPR) study demonstrated that the intestinal absorption of flavonoids can be related to atomic charges on carbon 3′ (Q_C3′), molecule surface area (SlogP_V3), balance between the center of mass and position of hydrophobic region (vsurf_ID1) and solvation energy of flavonoids (E_sol). These results provide useful information for initially screening of flavonoids with high intestinal absorption.

Density Functional Theory Study on the Complexation of NOTA as a Bifunctional Chelator with Radiometal Ions

1,4,7-Triazacyclononane-1,4,7-triacetic acid (NOTA) is a key bifunctional chelator utilized for the complexation of metal ions in radiopharmaceutical applications; the ability of these chelators depends on the strength of their binding with ions. The focus of the present work is to evaluate the complexation of Cu²⁺, Ga³⁺, Sc³⁺, and In³⁺ radiometal ions with NOTA using density functional theory (B3LYP functional) and 6-311+G(2d,2p)/DGDZVP basis sets. The significant role of ion-water interactions in the chelation interaction energies in solution reflects the competition between ion-water and NOTA-ion interaction in the chelation process. There is reasonable agreement between experimental and theoretical binding constants, geometries, and ¹H NMR chemical shifts. Chelation interaction energies, Gibbs free energies, and entropies in solution show that the NOTA-Ga³⁺ and NOTA-Cu²⁺ are the most and least stable complexes, respectively. The natural atomic charges and second order perturbation analysis reveal charge transfer between NOTA and radiometal ions. The theoretical ¹H NMR chemical shifts of NOTA are in good agreement with experiment; these values are influenced by the presence of the ions, which have a deshielding effect on the protons of NOTA. Global scalar properties such as E_HOMO/E_LUMO, ΔE_LUMO-HOMO, and chemical hardness/softness confirm that the NOTA-Cu²⁺ complex, which has a singly occupied molecular orbital, has the lowest ΔE_LUMO-HOMO value, the least chemical hardness, and the highest chemical softness. The significant variation of the hardness and ΔE_LUMO-HOMO values of the complexes can be attributed to the different positions of the metal ions on the periodic table. This study affirms that, among the radiometal ions, Ga³⁺ can be used to effectively radiolabel NOTA chelator for radiopharmaceutical usage as it binds most stably with NOTA.

Synthesis, structural characterization, QSAR and docking studies of a new binuclear nickel (II) complex based on the flexible tetradentate N-donor ligand as a potent antibacterial and anticancer agent

A new nickel (II)complex namely [Ni₂(Lt)Cl₄] derived from the NiCl₂.6H₂O and 1,1,3,3-tetrakis(3,5-dimethyl-1-pyrazolyl)propane (Lt) has been synthesized and fully characterized by the single crystal X-ray diffraction, elemental analysis, FT-IR, UV-vis, density functional theory (DFT) calculations, antibacterial and anticancer activities. In the title complex, each of the Ni(II) atoms is tetrahedrally coordinated by two N atoms from one of the chelating bidentate bis(3,5-dimethylpyrazolyl)methane units of the Lt ligand and two Cl as terminal ligands. The neighboring [Ni₂(Lt)Cl₄] molecules are linked together by the intermolecular CH⋯Cl hydrogen bonds to generate a 1D chain structure. The chains are further stabilized by the intermolecular CH⋯π interactions to form a two-dimensional non-covalent bonded structure. The antibacterial activity of the free Lt ligand and its Ni (II) complex shows that the ability of these compounds to inhibit growth of the tested bacteria increase from the Lt to binuclear nickel (II) complex. Scanning probe microscopy (SPM) study of the treated B. subtilis and E. coli bacteria was implemented to understand the structural changes caused by the interactions between the nickel (II) complex and the target bacteria. The cytotoxicity test of the Lt ligand and its complex was evaluated against the human carcinoma cell line (Caco-2) using the MTT assay. The results indicate that the Lt ligand and its complex display cytotoxicity against Caco-2 with the IC50 values of 36.29μM and 12.97μM, respectively. Therefore, the complex can be nominated as a potential anticancer agent. Molecular docking investigations on the five standard antibiotic, five standard anticancer drugs, free Lt ligand, title complex and twelve receptors were performed by Autodock vina function. The results of docking and DFT calculations are in line with the in vitro data obtained via the antibacterial and anticancer activity of Lt ligand and its made-complex.

Theoretical Reactivity Study of Indol-4-Ones and Their Correlation with Antifungal Activity

Chemical reactivity descriptors of indol-4-ones obtained via density functional theory (DFT) and hard–soft acid–base (HSAB) principle were calculated to prove their contribution in antifungal activity. Simple linear regression was made for global and local reactivity indexes. Results with global descriptors showed a strong relationship between antifungal activity vs. softness (S) (r = 0.98) for series I (6, 7a–g), and for series II (8a–g) vs. chemical potential (µ), electronegativity (χ) and electrophilicity (ω) (r = 0.86), p < 0.05. Condensed reactivity descriptors s_k⁺, ω_k⁻ for different fragments had strong relationships for series I and II (r = 0.98 and r = 0.92). Multiple linear regression was statistically significant for S (r = 0.98), η (r = 0.91), and µ/ω (r = 0.91) in series I. Molecular electrostatic potential maps (MEP) showed negative charge accumulation around oxygen of carbonyl group and positive accumulation around nitrogen. Fukui function isosurfaces showed that carbons around nitrogen are susceptible to electrophilic attack, whereas the carbon atoms of the carbonyl and phenyl groups are susceptible to nucleophilic attack for both series. The above suggest that global softness in conjunction with softness and electrophilicity of molecular fragments in enaminone systems and pyrrole rings contribute to antifungal activity of indol-4-ones.

The interaction of NOTA as a bifunctional chelator with competitive alkali metal ions: a DFT study

This work investigates NOTA–alkali metal (Li⁺, Na⁺ and K⁺ and Rb⁺) complexation using density functional theory.

The effect of N-methylation of amino acids (Ac-X-OMe) on solubility and conformation: a DFT study

N-Methylation of amino acid derivatives (Ac-X-OMe, X = Gly, Val, Leu, Ile, Phe, Met, Cys, Ser, Asp and His) leads to an increase in aqueous solubility, lipophilicity and lowering of the cis/trans amide conformational energy barrier (E_A).

A structure-activity relationship study of flavonoids as inhibitors of E. coli by membrane interaction effect

Flavonoids exhibit a broad range of biological activities including antibacterial activity. However, the mechanism of their antibacterial activity has not been fully investigated. The antibacterial activity and membrane interaction of 11 flavonoids (including 2 polymethoxyflavones and 4 isoflavonoids) against Escherichia coli were examined in this study. The antibacterial capacity was determined as flavonoids>polymethoxyflavones>isoflavonoids. Using fluorescence, it was observed that the 5 flavonoids rigidified the liposomal membrane, while the polymethoxyflavones and isoflavonoids increased membrane fluidity. There was a significant positive correlation between antibacterial capacity and membrane rigidification effect of the flavonoids. A quantitative structure-activity relationship (QSAR) study demonstrated that the activity of the flavonoid compounds can be related to molecular hydrophobicity (CLogP) and charges on C atom at position3 (C3). The QSAR model could be used to predict the antibacterial activity of flavonoids which could lead to natural compounds having important use in food and medical industry.

QSAR AND ACTION MECHANISM OF TROXACITABINE PRODRUGS WITH ANTITUMOR ACTIVITY

The quantitative structure–activity relationship (QSAR) of troxacitabine prodrugs with antitumor activity has been studied by using the density functional theory (DFT), molecular mechanics (MM2), and statistical methods. The established QSAR model shows not only significant statistical quality, but also predictive ability, with the square of adjusted correlation coefficient [Formula: see text] and the square of the cross-validation coefficient (q2 = 0.807). The antitumor activity is expressed as p IC50, which is defined as the negative value of the logarithm of necessary molar concentration of a compound to cause 50% growth inhibition against the human non-small-cell lung cancer cell line SW1573. It appears to be mainly governed by two factors (or original variables), i.e. the calculated hydrophobic coefficient (C log P) of whole molecule and the net charges of the first atom of substituent R (Q FR ), although three descriptors, i.e. C log P, (C log P)2, and Q FR , were selected in our multiple linear regression model. The factor C log P shows parabolic relation to p IC50 and its suitable range is around 5.6, and the other factor Q FR shows a significant negative correlation with p IC50. In this paper, a detailed discussion on these two factors was carried out, and their close correlation with the action mechanism of these prodrugs was reasonably revealed. Such results can offer some useful theoretical references for understanding the action mechanism and directing the molecular design of this kind of compound with antitumor activity.

Comparative study on antioxidant capacity of flavonoids and their inhibitory effects on oleic acid-induced hepatic steatosis in vitro

Nonalcoholic fatty liver disease (NAFLD) is the most common liver disease and its incidence is rising worldwide. We compared the antioxidant capacity of seventeen flavonoids with their inhibitory effects on oleic acid-induced triglyceride (TG) over-accumulation in HepG2 cells. The results showed significant correlations (P < 0.01) between the inhibition of intracellular TG levels and the suppression effects on reactive oxygen species. Nevertheless, the radical-reducing activities of flavonoids assessed by chemical assays (cyclic voltammetry and Folin-Ciocalteu reagent assay) were poorly correlated with their intracellular TG inhibitory effects. The relationships between structural properties of flavonoids and their inhibitory effects on TG over-accumulation were discussed.

A systematic computational study on flavonoids

17 selected flavones derivatives, flavonoids, were analyzed through a systematic B3LYP/6-311++G** computational study with the aim of understanding the molecular factors that determine their structural and energetic properties in gas phase. Flavonoids were selected in a systematic way to infer the effect of the number and relative positions of hydroxyl groups on molecular properties. Different conformers for each flavonoid were analyzed and the strength and topology of the intramolecular hydrogen bonds studied through the computation of the corresponding torsional profiles. Atoms in a Molecule, and Natural Bond Orbital methodology was applied to the analysis of charge distribution along the studied molecules, and the intramolecular hydrogen bonds. Molecular shapes were studied through full geometry optimization, and the position of the catechol ring analyzed through dihedral scans.

QSAR modeling for lipid peroxidation inhibition potential of flavonoids using topological and structural parameters

In the present study, Quantitative Structure-Activity Relationship (QSAR) modeling has been carried out for lipid peroxidation (LPO)-inhibition potential of a set of 27 flavonoids, using structural and topological parameters. For the development of models, three methods were used: (1) stepwise regression, (2) factor analysis followed by multiple linear regressions (FA-MLR) and (3) partial least squares (PLS) analysis. The best equation was obtained from stepwise regression analysis (Q2 = 0.626) considering the leave-oneout prediction statistics.

QSAR modeling of antiradical and antioxidant activities of flavonoids using electrotopological state (E-State) atom parameters

In the present paper QSAR modeling using electrotopological state atom (E-state) parameters has been attempted to determine the antiradical and the antioxidant activities of flavonoids in two model systems reported by Burda et al. (2001). The antiradical property of a methanolic solution of 1, 1-diphenyl-2-picrylhydrazyl (DPPH) and the antioxidant activity of flavonoids in a β-carotenelinoleic acid were the two model systems studied. Different statistical tools used in this communication are stepwise regression analysis, multiple linear regressions with factor analysis as the preprocessing step for variable selection (FA-MLR) and partial least squares analysis (PLS). In both the activities the best equation is obtained from stepwise regression analysis, considering, both equation statistics and predictive ability (antiradical activity: R 2 = 0.927, Q2 = 0.871 and antioxidant activity: R 2 = 0.901, Q2 = 0.841).

QSAR, action mechanism and molecular design of flavone and isoflavone derivatives with cytotoxicity against HeLa

The quantitative structure-activity relationship (QSAR) of 32 flavone and isoflavone derivatives with cytotoxicity expressed as pGC50, which is defined as the negative value of the logarithm of necessary molar concentration of this series of compounds to cause 50% growth inhibition against the human cervical epithelioid carcinoma cell line (HeLa), has been studied by using the density functional theory (DFT), molecular mechanics (MM2) and statistical methods. In order to obtain QSAR model with high predictive ability, the original dataset was randomly divided into a training set comprising 26 compounds and a test set comprising the rest 6 compounds. An optimal model for the training set with significant statistical quality (RA2=0.852) and predictive ability (q2=0.818) was established. The same model was further applied to predict pGC50 values of the 6 compounds in the test set, and the resulting predictive correlation coefficient Rpred2 reaches 0.738, further showing that this QSAR model has high predictive ability. It is very interesting to find that the cytotoxicities of these compounds against HeLa appear to be mainly governed by two quantum-chemical factors, i.e., the energy (ELUMO) of the lowest unoccupied molecular orbital (LUMO) and the net charges of C atom at site 6 on aromatic rings (QC6). Here the possible action mechanism of these compounds was analyzed and discussed in detail, in particular, the fact why the flavone derivatives have considerably higher cytotoxicity than isoflavone derivatives was reasonably explained. Based on this QSAR equation, 5 new compounds with higher cytotoxicity have been theoretically designed. Such results can offer useful theoretical references for experimental works.

A quantum mechanic/molecular mechanic study of the wild-type and N155S mutant HIV-1 integrase complexed with diketo acid

Integrase (IN) is one of the three human immunodeficiency virus type 1 (HIV-1) enzymes essential for effective viral replication. Recently, mutation studies have been reported that have shown that a certain degree of viral resistance to diketo acids (DKAs) appears when some amino acid residues of the IN active site are mutated. Mutations represent a fascinating experimental challenge, and we invite theoretical simulations for the disclosure of still unexplored features of enzyme reactions. The aim of this work is to understand the molecular mechanisms of HIV-1 IN drug resistance, which will be useful for designing anti-HIV inhibitors with unique resistance profiles. In this study, we use molecular dynamics simulations, within the hybrid quantum mechanics/molecular mechanics (QM/MM) approach, to determine the protein-ligand interaction energy for wild-type and N155S mutant HIV-1 IN, both complexed with a DKA. This hybrid methodology has the advantage of the inclusion of quantum effects such as ligand polarization upon binding, which can be very important when highly polarizable groups are embedded in anisotropic environments, for example in metal-containing active sites. Furthermore, an energy terms decomposition analysis was performed to determine contributions of individual residues to the enzyme-inhibitor interactions. The results reveal that there is a strong interaction between the Lys-159, Lys-156, and Asn-155 residues and Mg(2+) cation and the DKA inhibitor. Our calculations show that the binding energy is higher in wild-type than in the N155S mutant, in accordance with the experimental results. The role of the mutated residue has thus been checked as maintaining the structure of the ternary complex formed by the protein, the Mg(2+) cation, and the inhibitor. These results might be useful to design compounds with more interesting anti-HIV-1 IN activity on the basis of its three-dimensional structure.

A theoretical study of phenolic compounds with antioxidant properties

Quantum chemical calculations at the DFT/B3LYP, HF, and AM1 and PM3 semiempirical levels were employed to calculate a set of molecular properties for 41 phenol compounds with antioxidant activity. The significant molecular descriptors related to the compounds were the vertical ionization potentials (IPvs) and the charge on oxygen atom 7. The IPv has been calculated using Koopman's theorem IPv=-epsilon(HOMO(DFT)), IPv=-epsilon(HOMO(HF)) and as the difference of energy calculations for the corresponding cation and for the neutral form IPv(E(cation)-E(neutral)) obtained with the DFT/B3LYP method. The best model obtained showed not only the statistical significance but also predictive ability.

Structure–activity relationship study of flavone compounds with anti-HIV-1 integrase activity: A density functional theory study

Human immunodeficiency virus type-1 integrase (HIV-1 IN) is an essential enzyme for effective viral replication. Flavone compounds have been very much studied due to their activity during the inhibition process of HIV-1 IN. In this study, we employed density functional theory (DFT) using the B3LYP hybrid functional to calculate a set of molecular properties for 32 flavonoid compounds with anti-HIV-1 IN activity. The stepwise discriminant analysis (SDA), principal component analysis (PCA) and hierarchical cluster analysis (HCA) methods were employed to reduce dimensionality and investigate possible relationship between the calculated properties and the anti-HIV-1 IN activity. These analyses showed that the molecular hydrophobicity (ClogP), charge on atom 11 and electrophilic index (omega) are responsible for the separation between anti-HIV-1 IN active and inactive compounds.