Interaction between artemisinin and heme. A Density Functional Theory study of structures and interaction energies

Nanozyme-like single-atom catalyst combined with artesunate achieves photothermal-enhanced nanocatalytic therapy in the near-infrared biowindow

Selectively generating active free radical (AFR) in tumor microenvironment (TME) can promote irreversible oxidation of biomolecules and damage tumor cells, resulting in effective tumor inhibition. However, therapeutic efficacy of AFR-based tumor suppression approaches is often limited by insufficient amount of H₂O₂ or O₂ within TME. To overcome this obstacle, we design a pH/photothermal dual responsive nanosystem (PFeSA@AS) for combined photothermal and nanocatalytic therapy in the near-infrared biowindow. Here the Fe single-atom dispersed N, S-doped carbon nanosheets (FeSA) nanozyme is dispersed by phospholipid-polyethylene glycol-amine (DSPE-PEG-NH₂), and further loads artesunate (AS) via an amide reaction. Upon 808-nm laser irradiation in TME, the AS is released and further be catalyzed by the FeSA nanozyme to produce cytotoxic C-centered AFRs, and further be accelerated due to the photothermal conversion performance of FeSA (23.35%). The nanocatalytic process of FeSA nanozyme is realized by density functional theory (DFT). The tumor inhibition rates of a CT26 xenograft model is 92% through a photothermal-enhanced nanocatalytic synergistic therapy, and negligible systematic toxicity is observed. This work offers a potential paradigm of multifunctional single atomic catalysts (SACs) for enhancing tumor nanocatalytic therapy. STATEMENT OF SIGNIFICANCE: We designed a pH/photothermal dual responsive nanosystem (PFeSA@AS) for nanocatalytic therapy: (1) the nanosystem responsively releases AS under 808-nm laser irradiation in TME; (2) FeSA in the nanosystem can act as heme mimetic to convert AS into high cytotoxic C-centered free radicals for nanocatalytic therapy; (3) the photothermal conversion performance of FeSA further enhances the catalytic process to yield abundant AFR. Both in vitro and in vivo results demonstrate that this nanosystem can efficiently inhibit tumor growth through a photothermal-enhanced nanocatalytic synergistic therapy.

Engineered P450 Atom-Transfer Radical Cyclases are Bifunctional Biocatalysts: Reaction Mechanism and Origin of Enantioselectivity

New-to-nature radical biocatalysis has recently emerged as a powerful strategy to tame fleeting open-shell intermediates for stereoselective transformations. In 2021, we introduced a novel metalloredox biocatalysis strategy that leverages the innate redox properties of the heme cofactor of P450 enzymes, furnishing new-to-nature atom-transfer radical cyclases (ATRCases) with excellent activity and stereoselectivity. Herein, we report a combined computational and experimental study to shed light on the mechanism and origins of enantioselectivity for this system. Molecular dynamics and quantum mechanics/molecular mechanics (QM/MM) calculations revealed an unexpected role of the key beneficial mutation I263Q. The glutamine residue serves as an essential hydrogen bond donor that engages with the carbonyl moiety of the substrate to promote bromine atom abstraction and enhance the enantioselectivity of radical cyclization. Therefore, the evolved ATRCase is a bifunctional biocatalyst, wherein the heme cofactor enables atom-transfer radical biocatalysis, while the hydrogen bond donor residue further enhances the activity and enantioselectivity. Unlike many enzymatic stereocontrol rationales based on a rigid substrate binding model, our computations demonstrate a high degree of rotational flexibility of the allyl moiety in an enzyme-substrate complex and succeeding intermediates. Therefore, the enantioselectivity is controlled by the radical cyclization transition states rather than the substrate orientation in ground-state complexes in the preceding steps. During radical cyclization, anchoring effects of the Q263 residue and steric interactions with the heme cofactor concurrently control the π-facial selectivity, allowing for highly enantioselective C-C bond formation. Our computational findings are corroborated by experiments with ATRCase mutants generated from site-directed mutagenesis.

Chemometric methods in antimalarial drug design from 1,2,4,5-tetraoxanes analogues

A set of 23 steroidal 1,2,4,5-tetraoxane analogues were studied using quantum-chemical method (B3LYP/6-31 G*) and multivariate analyses (PCA, HCA, KNN and SIMCA) in order to calculate the properties and correlate them with antimalarial activity (log RA) against Plasmodium falciparum clone D-6 from Sierra Leone. PCA results indicated 99.94% of the total variance and it was possible to divide the compounds into two classes: less and more active. Descriptors responsible for separating were: highest occupied molecular orbital energy (HOMO), bond length (O1-O2), Mulliken electronegativity (χ) and Bond information content (BIC0). We use HCA, KNN and SIMCA to explain relationships between molecular properties and biological activity of a training set and to predict antimalarial activity (log RA) of 13 compounds (#24-36) with unknown biological activity. We apply molecular docking simulations to identify intermolecular interactions with a selected biological target. The results obtained in multivariate analysis aided in the understanding of the activity of the new compound's design (#24-36). Thus, through chemometric analyses and docking molecular study, we propose theoretical synthetic routes for the most promising compounds 28, 30, 32 and 36 that can proceed to synthesis steps and in vitro and in vivo assays.

Activation of artemisinin and heme degradation in Leishmania tarentolae promastigotes: A possible link

Endoperoxides (EPs) appear to be promising drug candidates against protozoal diseases, including malaria and leishmaniasis. Previous studies have shown that these drugs need an intracellular activation to exert their pharmacological potential. The efficiency of these drugs is linked to the extensive iron demand of these intracellular protozoal parasites. An essential step of the activation mechanism of these drugs is the formation of radicals in Leishmania. Iron is a known trigger for intracellular radical formation. However, the activation of EPs by low molecular iron or by heme iron may strongly depend on the structure of the EPs themselves. In this study, we focused on the activation of artemisinin (Art) in Leishmania tarentolae promastigotes (LtP) in comparison to reference compounds. Viability assays in different media in the presence of different iron sources (hemin/fetal calf serum) showed that IC50 values of Art in LtP were modulated by assay conditions, but overall were within the low micromolar range. Low temperature electron paramagnetic resonance (EPR) spectroscopy of LtP showed that Art shifted the redox state of the labile iron pool less than the EP ascaridole questioning its role as a major activator of Art in LtP. Based on the high reactivity of Art with hemin in previous biomimetic experiments, we focused on putative heme-metabolizing enzymes in Leishmania, which were so far not well described. Inhibitors of mammalian heme oxygenase (HO; tin and chromium mesoporphyrin) acted antagonistically to Art in LtP and boosted its IC50 value for several magnitudes. By inductively coupled plasma methods (ICP-OES, ICP-MS) we showed that these inhibitors do not block iron (heme) accumulation, but are taken up and act within LtP. These inhibitors blocked the conversion of hemin to bilirubin in LtP homogenates, suggesting that an HO-like enzyme activity in LtP exists. NADPH-dependent degradation of Art and hemin was highest in the small granule and microsomal fractions of LtP. Photometric measurements in the model Art/hemin demonstrated that hemin requires reduction to heme and that subsequently an Art/heme complex (λmax 474 nm) is formed. EPR spin-trapping in the system Art/hemin revealed that NADPH, ascorbate and cysteine are suitable reductants and finally activate Art to acyl-carbon centered radicals. These findings suggest that heme is a major activator of Art in LtP either via HO-like enzyme activities and/or chemical interaction of heme with Art.

Elimination of Schistosoma mansoni in infected mice by slow release of artemisone

The current treatment of schistosomiasis is based on the anti-helminthic drug praziquantel (PZQ). PZQ affects only the adult stages of schistosomes. In addition, resistance to PZQ is emerging. We suggest a drug, which could serve as a potential alternative or complement to PZQ, and as a means of treating infections at earlier, pre-granuloma stage. Derivatives of the peroxidic antimalarial drug artemisinin have been indicated as alternatives, because both plasmodia and schistosomes are blood-feeding parasites. The mechanism of action of artemisinins is related to oxidative effects of the artemisinins on intracellular reductants leading to formation of cytotoxic reactive oxygen species. We used artemisone, which has improved pharmacokinetics and anti-plasmodial activity, and reduced toxicity compared to other artemisinins in clinical use against malaria. We infected adult mice by subcutaneous injection of S. mansoni cercariae (about 200) and treated them at various times post infection by the following methods: i. artemisone suspension administered by gavage (400-450 mg/kg); ii. subcutaneous injection of a gel containing a known concentration of artemisone (115-120 mg/kg); iii. subcutaneous insertion of the drug incorporated in a solid polymer (56-60 mg/kg); iv. intraperitoneal injection of the drug solubilized in DMSO (115-120 mg/kg). Drug administration in polymers was performed to enable slow release of the artemisone that was verified in vivo and in vitro bioassays using drug-sensitive malaria parasites. We found superior strong anti-schistosome effects up to a total reduction of worm number, mainly following repetitive treatments with the drug absorbed in the polymers (73.1% and 95.9% reduction in mice treated with artemisone in gel 7 and 14, and 21, 28 and 35 days post infection, respectively). The results indicate that artemisone has a potent anti-schistosome activity. Its main importance in this context is its effectiveness in treating hosts harboring juvenile schistosomes, before egg-deposition and induction of deleterious immune responses.

The interaction of heme with plakortin and a synthetic endoperoxide analogue: new insights into the heme-activated antimalarial mechanism

In the present work we performed a combined experimental and computational study on the interaction of the natural antimalarial endoperoxide plakortin and its synthetic analogue 4a with heme. Obtained results indicate that the studied compounds produce reactive carbon radical species after being reductively activated by heme. In particular, similarly to artemisinin, the formation of radicals prone to inter-molecular reactions should represent the key event responsible for Plasmodium death. To our knowledge this is the first experimental investigation on the reductive activation of simple antimalarial endoperoxides (1,2-dioxanes) by heme and results were compared to the ones previously obtained from the reaction with FeCl₂. The obtained experimental data and the calculated molecular interaction models represent crucial tools for the rational optimization of our promising class of low-cost synthetic antimalarial endoperoxides.

A QSAR, pharmacokinetic and toxicological study of new artemisinin compounds with anticancer activity

The Density Functional Theory (DFT) method and the 6-31G** basis set were employed to calculate the molecular properties of artemisinin and 20 derivatives with different degrees of cytotoxicity against the human hepatocellular carcinoma HepG2 line. Principal component analysis (PCA) and hierarchical cluster analysis (HCA) were employed to select the most important descriptors related to anticancer activity. The significant molecular descriptors related to the compounds with anticancer activity were the ALOGPS_log, Mor29m, IC5 and GAP energy. The Pearson correlation between activity and most important descriptors were used for the regression partial least squares (PLS) and principal component regression (PCR) models built. The regression PLS and PCR were very close, with variation between PLS and PCR of R² = ±0.0106, R²_ajust = ±0.0125, s = ±0.0234, F_(4,11) = ±12.7802, Q² = ±0.0088, SEV = ±0.0132, PRESS = ±0.4808 and S_PRESS = ±0.0057. These models were used to predict the anticancer activity of eight new artemisinin compounds (test set) with unknown activity, and for these new compounds were predicted pharmacokinetic properties: human intestinal absorption (HIA), cellular permeability (P_CaCO2), cell permeability Maden Darby Canine Kidney (P_MDCK), skin permeability (P_Skin), plasma protein binding (PPB) and penetration of the blood-brain barrier (C_Brain/Blood), and toxicological: mutagenicity and carcinogenicity. The test set showed for two new artemisinin compounds satisfactory results for anticancer activity and pharmacokinetic and toxicological properties. Consequently, further studies need be done to evaluate the different proposals as well as their actions, toxicity, and potential use for treatment of cancers.

A SAR and QSAR study of new artemisinin compounds with antimalarial activity

The Hartree-Fock method and the 6-31G** basis set were employed to calculate the molecular properties of artemisinin and 20 derivatives with antimalarial activity. Maps of molecular electrostatic potential (MEPs) and molecular docking were used to investigate the interaction between ligands and the receptor (heme). Principal component analysis and hierarchical cluster analysis were employed to select the most important descriptors related to activity. The correlation between biological activity and molecular properties was obtained using the partial least squares and principal component regression methods. The regression PLS and PCR models built in this study were also used to predict the antimalarial activity of 30 new artemisinin compounds with unknown activity. The models obtained showed not only statistical significance but also predictive ability. The significant molecular descriptors related to the compounds with antimalarial activity were the hydration energy (HE), the charge on the O11 oxygen atom (QO11), the torsion angle O1-O2-Fe-N2 (D2) and the maximum rate of R/Sanderson Electronegativity (RTe⁺). These variables led to a physical and structural explanation of the molecular properties that should be selected for when designing new ligands to be used as antimalarial agents.

Synergy of ferrous ion on 5-aminolevulinic acid-mediated growth inhibition of Plasmodium falciparum

Haem biosynthesis appeared to be a target of malaria therapy because 5-aminolevulinic acid (ALA), a haem biosynthesis starting material, with light exposure or a high amount of ALA alone reduced Plasmodium falciparum growth to undetectable level. However, the administration of a high dose of ALA is unrealistic for clinical therapy. We found that Fe(2+) enhanced P. falciparum-killing potency of ALA and significantly inhibited the parasite growth. The intermediates of haem biosynthesis localized to the parasite organelles, and coproporphyrin III was the most accumulated intermediate. These novel findings may lead to development of a new anti-malarial drug using ALA and Fe(2+).

A plausible mechanism for the antimalarial activity of artemisinin: A computational approach

Artemisinin constitutes the frontline treatment to aid rapid clearance of parasitaemia and quick resolution of malarial symptoms. However, the widespread promiscuity about its mechanism of action is baffling. There is no consensus about the biochemical target of artemisinin but recent studies implicate haem and PfATP6 (a calcium pump). We investigated the role of iron and artemisinin on PfATP6, in search of a plausible mechanism of action, via density functional theory calculations, docking and molecular dynamics simulations. Results suggest that artemisinin gets activated by iron which in turn inhibits PfATP6 by closing the phosphorylation, nucleotide binding and actuator domains leading to loss of function of PfATP6 of the parasite and its death. The mechanism elucidated here should help in the design of novel antimalarials.

Molecular modeling and evaluation of binding mode and affinity of artemisinin-quinine hybrid and its congeners with Fe-protoporphyrin-IX as a putative receptor

A recent rational approach to anti-malarial drug design is characterized as "covalent biotherapy" involves linking of two molecules with individual intrinsic activity into a single agent, thus packaging dual activity into a single hybrid molecule. In view of this background and reported anti malaria synergism between artemisinin and quinine; we describe the computer-assisted docking to predict molecular interaction and binding affinity of Artemisinin-Quinine hybrid and its derivatives with the intraparasitic haeme group of human haemoglobin. Starting from a crystallographic structure of Fe-protoporphyrin-IX, binding modes, orientation of peroxide bridge (Fe-O distance), docking score and interaction energy are predicted using the docking molecular mechanics based on generalized Born/surface area (MM-GBSA) solvation model. Seven new ligands were identified with a favourable glide score (XP score) and binding free energy (ΔG) with reference to the experimental structure from a data set of thirty four hybrid derivatives. The result shows the conformational property of the drug-receptor interaction and may lead to rational design and synthesis of improved potent artemisinin based hybrid antimalarial that target haemozoin formation.

COMPUTATIONAL MODELING OF ANTIMALARIAL 10-SUBSTITUTED DEOXOARTEMISININS

Nineteen 10-substitued deoxoartemisinin derivatives and artemisinin with activity against D-6 strains of malarial falciparum designated as Sierra Leone are studied. We use molecular electrostatic potential maps in an attempt to identify key structural features of the artemisinins that are necessary for their activities and molecular docking to investigate the interaction with the molecular receptor (heme). Chemometric modeling: Principal Component Analysis (PCA), Hierarchical Cluster Analysis (HCA), K-Nearest Neighbor (KNN), Soft Independent Modeling of Class Analogy (SIMCA) and Stepwise Discriminant Analysis (SDA) are employed to reduce dimensionality and investigate which subset of descriptors are responsible for the classification between more active (MA) and less active (LA) artemisinins. The PCA, HCA, KNN, SIMCA and SDA studies showed that the descriptors LUMO (Lowest Unoccupied Molecular Orbital) energy, DFeO1 (Distance between the O 1 atom from ligand and iron atom from heme), X1A (Average Connectivity Index Chi-1) and Mor15u (Molecular Representation of Structure Based on Electron Diffraction) code of signal 15, unweighted, are responsible for separating the artemisinins according to their degree of antimalarial activity. The prediction study was done with a new set of eight artemisinins by using the chemometric methods and five of them were predicted as active against D-6 strains of falciparum malaria. In order to verify if the key structural features that are necessary for their antimalarial activities were investigated for the interaction with the heme, we also carried out calculations of the molecular electrostatic potential (MEP) and molecular docking. MEP maps and molecular docking were analyzed for more active compounds of the prediction set.

A topological study of the decomposition of 6,7,8-trioxabicyclo[3.2.2]nonane induced by Fe(II): modeling the artemisinin reaction with heme

We report a theoretical study on the electronic and topological aspects of the reaction of dihydrated Fe(OH)(2) with 6,7,8-trioxabicyclo[3.2.2]nonane, as a model for the reaction of heme with artemisinin. A comparison is made with the reaction of dihydrated ferrous hydroxide with O(2), as a model for the heme interaction with oxygen. We found that dihydrated Fe(OH)(2) reacts more efficiently with the artemisinin model than with O(2). This result suggests that artemisinin instead of molecular oxygen would interact with heme, disrupting its detoxification process by avoiding the initial heme to hemin oxidation, and killing in this way the malaria parasite. The ELF and AIM theories provide support for such a conclusion, which further clarifies our understanding on how artemisinin acts as an antimalarial agent.

A study of new antimalarial artemisinins through molecular modeling and multivariate analysis

Artemisinin and 18 derivatives with antimalarial activity against W-2 strains of Plasmodium falciparum were studied through quantum chemistry and multivariate analysis. The geometry optimization of the structures was realized with the Hartree-Fock (HF) theory and 3-21G basis set. Maps of molecular electrostatic potential (MEP) and molecular docking were used to investigate the interaction between the ligands and the receptor (heme). Principal Component Analysis (PCA) and Hierarchical Cluster Analysis (HCA) were employed to select the most important descriptors related to activity. A predictive model was generated by the Partial Least Square (PLS) method through 15 molecules and 4 used as an external validation set, which were selected in the training set, the validation parameters of which are Q2 = 0.85 and R2 = 0.86. The model included as molecular parameters, the radial distribution function, RDF060e, the hydration energy, HE, and the distance between the O1 atom from the ligand and the iron atom from heme, d(Fe-O1). Thus, the synthesis of new derivatives may follow the results of the MEP maps and the PLS analysis.