MNDO/d calculations on the interaction between artemisinin and heme

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.

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

Malaria is an infectious disease caused by the unicellular parasite Plasmodium sp. Currently, the malaria parasite is becoming resistant to the traditional pharmacological alternatives, which are ineffective. Artemisinin is the most recent advance in the chemotherapy of malaria. Since it has been proven that artemisinin may act on intracellular heme, we have undertaken a systematic study of several interactions and arrangements between artemisinin and heme. Density Functional Theory calculations were employed to calculate interaction energies, electronic states, and geometrical arrangements for the complex between the heme group and artemisinin. The results show that the interaction between the heme group and artemisinin at long distances occurs through a complex where the iron atom of the heme group retains its electronic features, leading to a quintet state as the most stable one. However, for interaction at short distances, due to artemisinin reduction by the heme group, the most stable complex has a septet spin state. These results suggest that a thermodynamically favorable interaction between artemisinin and heme may happen.

DFT study of the reductive decomposition of artemisinin

Artemisinin is a sesquiterpene lactone with an endoperoxide function that is essential for its antimalarial activity. The DFT B3LYP method, together with the 6-31G(d) and 6-31+G(d,p) basis set, is employed to calculate a set of radical anions and neutral species supposed to be formed during the rearrangement of artemisinin from the two radicals (C-centered and O-centered) that are supposed to play a relevant role in the mechanism of action. The B3LYP results show that the primary and the secondary radicals centered on C(4), generated by homolytic break of the C(3)-C(4) bond and by 1,5 hydrogen shift, respectively, are more stable than radicals centered on oxygen. The calculations show that the activation barriers for rearrangements are low, leading to a thermodynamically favorable process. These results reinforce our previous conclusions based on semi-empirical calculations but also give additional information on the reductive decomposition of artemisinin.

On the O2 binding of Fe–porphyrin, Fe–porphycene, and Fe–corrphycene complexes

Based on our previous study for the O2 binding of the Fe-Por complex, this study investigates the O2 binding mechanism in the Fe-porphyrin isomers, Fe-porphycene (FePc), and Fe-corrphycene (FeCor) complexes. By calculating the potential energy surface of the O2 binding, the present study explains the reason for the dramatic increase of O2 affinities observed in the FePc complex. In the case of FeCor-O2, the O2 binding process includes the intersystem crossing from a triplet to singlet state, as in the FePor-O2 complex. However, FePc-O2 uses only a singlet surface. This is because the ground state of the FePc complex in the deoxy state is a triplet state, while those of FePor and FeCor are a quintet state. Such difference originates from character of the SOMO. We estimated an equilibrium constant for the O2 binding that reasonably reproduced the trend observed in the experiments.

On the reversible O2 binding of the Fe–porphyrin complex

Electronic mechanism of the reversible O(2) binding by heme was studied by using Density Functional Theory calculations. The ground state of oxyheme was calculated to be open singlet state [Fe(S =1/2) + O(2)(S = 1/2)]. The potential energy surface for singlet state is associative, while that for triplet state is dissociative. Because the ground state of the O(2)+ deoxyheme system is triplet in the dissociation limit [Fe(S = 2) + O(2)(S = 1)], the O(2) binding process requires relativistic spin-orbit interaction to accomplish the intersystem crossing from triplet to singlet states. Owing to the singlet-triplet crossing, the activation energies for both O(2) binding and dissociation become moderate, and hence reversible. We also found that the deviation of the Fe atom from the porphyrin plane is also important reaction coordinate for O(2) binding. The potential surface is associative/dissociative when the Fe atom locates in-plane/out-of-plane.

Quantitative structure-activity relationship in aziridinyl-1,4-naphthoquinone antimalarials: study of theoretical correlations by the PM3 method

Several molecular parameters for 2,3,5-substituted 1,4-naphthoquinones including 2-aziridinyl and 2,3-aziridinyl-1,4-naphthoquinones with antimalarial activities were obtained with the semi-empirical PM3 method. The descriptor related to the Gibbs free energy of an isodesmic equation defining the reduction of the naphthoquinones was found to have high correlation with activity. The quantitative models reported clearly show a dependence of activity on the redox potential for reduction of the naphthoquinones. Compounds with lower values of DeltaG for reduction are more active than those with higher values of DeltaG.