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.     

Ferryl-oxo heme intermediate in the antimalarial mode of action of artemisinin

Fourier transform infrared (FTIR) and resonance Raman (RR) spectroscopies have been employed to investigate the reductive cleavage of the O-O bond of the endoperoxide moiety of the antimalarial drug artemisinin and its analog trioxane alcohol by hemin dimer. We have recorded FTIR spectra in the nu(O-O) and nu(as)(Fe-O-Fe) regions of artemisinin and of the hemin dimer that show the cleavage of the endoperoxide and that of the hemin dimer, respectively. We observed similar results in the trioxane alcohol/hemin dimer reaction. The RR spectrum of the artemisinin/hemin dimer reaction displays a vibrational mode at 850 cm(-1) that shifts to 818 cm(-1) when the experiment is repeated with (18)O-O(18) endoperoxide enriched trioxane alcohol. The frequency of this vibration and the magnitude of the (18)O-O(18) isotopic shift led us to assign the 850 cm(-1) mode to the Fe(IV) = O stretching vibration of a ferryl-xoxo heme intermediate that occurs in the artemisinin/hemin dimer and trioxane alcohol/hemin reactions. These results provide the first direct characterization of the antimalarial mode of action of artemisinin and its trioxane analog, and suggest that artemisinin appears to react with heme molecules that have been incorporated into hemozoin and subsequently the heme performs cytochrome P450-type chemistry.     

The study on the interaction between seryl-histidine dipeptide and proteins by circular dichroism and molecular modeling

The selective cleavage of proteins is very important in key biological processes. Chemical (nonenzymatic) reagents such as cyanogen bromide and transition metal complexes are used extensively with great defects. In this paper, the binding of seryl-histidine dipeptide (abbreviated as SH) with bovine serum albumin (BSA) and lysozyme were investigated by the circular dichroism spectroscopy (CD) at 298K, molecular docking studies and quantum chemical calculations based on the previous results of polyacrylamide gel electrophoresis (PAGE). From the studies of CD, it showed that SH interacted strongly with BSA and lysozyme. The change percentages of the secondary structures of BSA and lysozyme were calculated. The contents of the beta-sheets decreased remarkably. It indicated that the interactions between SH and proteins could break the hydrogen bonds of beta-sheets selectively. The docking studies between SH and BSA showed that the position of the oxygen atom of the hydroxyl group of SH (O(12)) was in favor of a nucleophilic attack on carbon atom of the amide bond of a beta-sheet (C(34)) because the distance between O(12) and C(34) was 3.37A. Natural charges, natural atomic hybrid percentages and square sums of HOMO coefficients calculated by the NBO and population analysis at HF/6-31G* supported the suggested mechanism. And so SH may be an interesting agent for the therapeutic use.     

Artemisinin, an endoperoxide antimalarial, disrupts the hemoglobin catabolism and heme detoxification systems in malarial parasite.

Endoperoxide antimalarials based on the ancient Chinese drug Qinghaosu (artemisinin) are currently our major hope in the fight against drug-resistant malaria. Rational drug design based on artemisinin and its analogues is slow as the mechanism of action of these antimalarials is not clear. Here we report that these drugs, at least in part, exert their effect by interfering with the plasmodial hemoglobin catabolic pathway and inhibition of heme polymerization. In an in vitro experiment we observed inhibition of digestive vacuole proteolytic activity of malarial parasite by artemisinin. These observations were further confirmed by ex vivo experiments showing accumulation of hemoglobin in the parasites treated with artemisinin, suggesting inhibition of hemoglobin degradation. We found artemisinin to be a potent inhibitor of heme polymerization activity mediated by Plasmodium yoelii lysates as well as Plasmodium falciparum histidine-rich protein II. Interaction of artemisinin with the purified malarial hemozoin in vitro resulted in the concentration-dependent breakdown of the malaria pigment. Our results presented here may explain the selective and rapid toxicity of these drugs on mature, hemozoin-containing, stages of malarial parasite. Since artemisinin and its analogues appear to have similar molecular targets as chloroquine despite having different structures, they can potentially bypass the quinoline resistance machinery of the malarial parasite, which causes sublethal accumulation of these drugs in resistant strains.     

A model binding site for testing scoring functions in molecular docking

Prediction of interaction energies between ligands and their receptors remains a major challenge for structure-based inhibitor discovery. Much effort has been devoted to developing scoring schemes that can successfully rank the affinities of a diverse set of possible ligands to a binding site for which the structure is known. To test these scoring functions, well-characterized experimental systems can be very useful. Here, mutation-created binding sites in T4 lysozyme were used to investigate how the quality of atomic charges and solvation energies affects molecular docking. Atomic charges and solvation energies were calculated for 172,118 molecules in the Available Chemicals Directory using a semi-empirical quantum mechanical approach by the program AMSOL. The database was first screened against the apolar cavity site created by the mutation Leu99Ala (L99A). Compared to the electronegativity-based charges that are widely used, the new charges and desolvation energies improved ranking of known apolar ligands, and better distinguished them from more polar isosteres that are not observed to bind. To investigate whether the new charges had predictive value, the non-polar residue Met102, which forms part of the binding site, was changed to the polar residue glutamine. The structure of the resulting Leu99Ala and Met102Gln double mutant of T4 lysozyme (L99A/M102Q) was determined and the docking calculation was repeated for the new site. Seven representative polar molecules that preferentially docked to the polar versus the apolar binding site were tested experimentally. All seven bind to the polar cavity (L99A/M102Q) but do not detectably bind to the apolar cavity (L99A). Five ligand-bound structures of L99A/M102Q were determined by X-ray crystallography. Docking predictions corresponded to the crystallographic results to within 0.4A RMSD. Improved treatment of partial atomic charges and desolvation energies in database docking appears feasible and leads to better distinction of true ligands. Simple model binding sites, such as L99A and its more polar variants, may find broad use in the development and testing of docking algorithms.     

A computational comparison of the atomic models of the actomyosin interface

Several atomic models of the actomyosin interface have been proposed based on the docking together of their component structures using electron microscopy and resonance energy-transfer measurements. Although these models are in approximate agreement in the location of the binding interfaces when myosin is tightly bound to actin, their relationships to molecular docking simulations based on computational free-energy calculations are investigated here. Both rigid-docking and flexible-docking conformational search strategies were used to identify free-energy minima at the interfaces between atomic models of myosin and actin. These results suggest that the docking model produced by resonance energy-transfer data is closer to a free-energy minimum at the interface than are the available atomic models based on electron microscopy. The conformational searches were performed using both scallop and chicken skeletal muscle myosins and identified similarly oriented actin-binding interfaces that serve to validate that these models are at the global minimum. These results indicate that the existing docking models are close to but not precisely at the lowest-energy initial contact site for strong binding between myosin and actin that should represent an initial contact between the two proteins; therefore, conformational changes are likely to be important during the transition to a strongly bound complex.     

Surface Charges and Regulation of FMN to Heme Electron Transfer in Nitric-oxide Synthase

The nitric-oxide synthases (NOS, EC 1.14.13.39) are modular enzymes containing attached flavoprotein and heme (NOSoxy) domains. To generate nitric oxide (NO), the NOS FMN subdomain must interact with the NOSoxy domain to deliver electrons to the heme for O₂ activation during catalysis. The molecular basis and how the interaction is regulated is unclear. We explored the role of eight positively charged residues that create an electropositive patch on NOSoxy in enabling the electron transfer by incorporating mutations that neutralized or reversed their individual charges. Stopped-flow and steady-state experiments revealed that individual charges at Lys⁴²³, Lys⁶²⁰, and Lys⁶⁶⁰ were the most important in enabling heme reduction in nNOS. Charge reversal was more disruptive than neutralization in all cases, and the effects on heme reduction were not due to a weakening in the thermodynamic driving force for heme reduction. Mutant NO synthesis activities displayed a complex pattern that could be simulated by a global model for NOS catalysis. This analysis revealed that the mutations impact the NO synthesis activity only through their effects on heme reduction rates. We conclude that heme reduction and NO synthesis in nNOS is enabled by electrostatic interactions involving Lys⁴²³, Lys⁶²⁰, and Lys⁶⁶⁰, which form a triad of positive charges on the NOSoxy surface. A simulated docking study reveals how electrostatic interactions of this triad can enable an FMN-NOSoxy interaction that is productive for electron transfer.     

A crystallographic and theoretical study of an (E)-2-Hydroxyiminoethanone derivative: prediction of cyclooxygenase inhibition selectivity of stilbenoids by MM-PBSA and the role of atomic charge

We recently reported that the hydroxyiminoethanone derivative, (E)-OXM, behaves as a highly selective COX-1 inhibitor (COX-1 SI = 833), and also an interesting scaffold with unique characteristics. In the current study, a comprehensive crystallographic and computational study was performed to elucidate its conformational stability and pharmacological activity. Its conformational energy was studied at the B3LYP/6-311G** level of theory and compared to the single-crystal X-ray diffraction data. In addition, computational studies of three structurally different stilbenoid derivatives used as selective COX-1 or COX-2 inhibitors were undertaken to predict their COX selectivity potentials. Flexible docking was performed for all compounds at the active site of both COX-1 and COX-2 enzymes by considering some of the key residues as flexible during the docking operation. In the next step, molecular dynamic simulation and binding free energy calculations were performed by MM-PBSA. Final results were found to be highly dependent on the atomic charges of the inhibitors and the choice of force field used to calculate the atomic charges. The binding conformation of the hydroxyiminoethanone derivative is highly correlated with the type of COX isoform inhibited. Our predictive approach can truly predict the cyclooxygenase inhibition selectivity of stilbenoid inhibitors.     

青蒿素及其类似物抗疟构效关系的DFT研究

为了从原子水平上揭示青蒿素及其类似物的结构与抗疟活性之间的关系,运用密度泛函理论DFT方法,在B3LYP/6-31G*水平上对青蒿素及其类似物二氢青蒿素、蒿甲醚和青蒿琥酯的结构和性质进行了理论计算。从分子的平衡构型、Wiberg键级、溶剂化能、偶极矩和静电势等方面分析了青蒿素及其类似物的抗疟构效关系。结果表明,青蒿素及其类似物结构中七元环上的过氧桥键、醚氧键以及六元环上的内酯结构是其抗疟作用的关键活性位,过氧桥键处负的静电势越多,青蒿素与血红素的相互作用越强,分子的抗疟活性越强。理论预测四个药物分子的抗疟活性顺序为:青蒿素<二氢青蒿素<蒿甲醚<青蒿琥酯,与实验活性结果一致。     

Structure and electronic properties of heme complexes with various ligands

The electronic and atomic structures, and the molecular dynamics of the atomic structure at 310 K of a set of heme complexes with His and Gly amino acids in the 5th coordination position and some ligands (O2, NO) in the 6th position were studied by ab initio (3-21G basis set) and semiempirical (PM3) quantum chemistry methods and the method of molecular dynamics. It was shown that the type of coordination of the imidazole ring influences the constant of chemical bonding of molecular oxygen of the complexes. On the other hand, NO and O2 molecules have different transinfluence on the ligand in the 5th coordination position. It was shown that temperature affects profoundly the atomic and electronic structures of the complexes, the tightness of chemical bonding and their reactivity.     

Active site structure of SoxB-type cytochrome bo3 oxidase from thermophilic Bacillus

Two-subunit SoxB-type cytochrome c oxidase in Bacillus stearothermophilus was over-produced, purified, and examined for its active site structures by electron paramagnetic resonance (EPR) and resonance Raman (RR) spectroscopies. This is cytochrome bo3 oxidase containing heme B at the low-spin heme site and heme O at the high-spin heme site of the binuclear center. EPR spectra of the enzyme in the oxidized form indicated that structures of the high-spin heme O and the low-spin heme B were similar to those of SoxM-type oxidases based on the signals at g=6.1, and g=3.04. However, the EPR signals from the CuA center and the integer spin system at the binuclear center showed slight differences. RR spectra of the oxidized form showed that heme O was in a 6-coordinated high-spin (nu3 = 1472 cm(-1)), and heme B was in a 6-coordinated low-spin (nu3 = 1500 cm(-1)) state. The Fe2+-His stretching mode was observed at 211 cm(-1), indicating that the Fe2+-His bond strength is not so much different from those of SoxM-type oxidases. On the contrary, both the Fe2+-CO stretching and Fe2+-C-O bending modes differed distinctly from those of SoxM-type enzymes, suggesting some differences in the coordination geometry and the protein structure in the proximity of bound CO in cytochrome bo3 from those of SoxM-type enzymes.     

Molecular docking and 3-D-QSAR studies on the possible antimalarial mechanism of artemisinin analogues

Artemisinin (Qinghaosu) is a natural constituent found in Artemisia annua L, which is an effective drug against chloroquine-resistant Plasmodium falciparum strains and cerebral malaria. The antimalarial activities of artemisinin and its analogues appear to be mediated by the interactions of the drugs with hemin. In order to understand the antimalarial mechanism and the relationship between the physicochemical properties and the antimalarial activities of artemisinin analogues, we performed molecular docking simulations to probe the interactions of these analogues with hemin, and then performed three-dimensional quantitative structure-activity relationship (3-D-QSAR) studies on the basis of the docking models employing comparative molecular force fields analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Molecular docking simulations generated probable 'bioactive' conformations of artemisinin analogues and provided a new insight into the antimalarial mechanism. The subsequent partial least squares (PLS) analysis indicates that the calculate binding energies correlate well with the experimental activity values. The CoMFA and CoMSIA models based on the bioactive conformations proved to have good predictive ability and in turn match well with the docking result, which further testified the reliability of the docking model. Combining these results, that is molecular docking and 3-D-QSAR, together, the binding model and activity of new synthesized artemisinin derivatives were well explained.     

Partial atomic charges of amino acids in proteins

Using a semiempirical quantum mechanical procedure (FCPAC) we have calculated the partial atomic charges of amino acids from 494 high-resolution protein structures. To analyze the influence of the protein's environment, we considered each residue under two conditions: either as the center of a tripeptide with PDB structure geometry (free) or as the center of 13-16 amino acid clusters extracted from the PDB structure (buried). The partial atomic charges from residues in helices and in sheets were separated. The FCPAC partial atomic charges of the Cbeta and Calpha of most residues correlate with their helix propensity, positively for Cbeta and negatively for Calpha (r2 = 0.76 and 0.6, respectively). The main consequence of burying residues in proteins is the polarization of the backbone C=O bond, which is more pronounced in helices than in sheets. The average shift of the oxygen partial charges that results from burying is -0.120 in helix and -0.084 in sheet with the charge of the proton as unit. Linear correlations are found between the average NMR chemical shifts and the average FCPAC partial charges of Calpha (r2 = 0.8-0.85), N (r3 = 0.67-0.72), and Cbeta (r2 = 0.62) atoms. Correlations for helix and beta-sheet FCPAC partial charges show parallel regressions, suggesting that the charge variations due to burying in proteins differentiate between the dihedral angle effects and the polarization of backbone atoms.     

Structural dynamics of myoglobin: effect of internal cavities on ligand migration and binding

Using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures, we have studied CO binding to the heme and CO migration among cavities in the interior of sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation. Photoproduct intermediates, characterized by CO in different locations, were selectively enhanced by laser illumination at specific temperatures. Measurements were performed on the wild-type protein and a series of mutants (L104W, I107W, I28F, and I28W) in which bulky amino acid side chains were introduced to block passageways between cavities or to fill these sites. Binding of xenon was also employed as an alternative means of filling cavities. In all samples, photolyzed CO ligands were observed to initially bind at primary docking site B in the vicinity of the heme iron, from where they migrate to the secondary docking sites, the Xe4 and/or Xe1 cavities. To examine the relevance of these internal docking sites for physiological ligand binding, we have performed room-temperature flash photolysis on the entire set of proteins in the CO- and O(2)-bound form. Together with the cryospectroscopic results, these data provide a clear picture of the role of the internal sites for ligand escape from and binding to myoglobin.     

Heme Mediates Cytotoxicity from Artemisinin and Serves as a General Anti-Proliferation Target

Heme (Fe2+ protoporphyrin IX) is an essential molecule that has been implicated the potent antimalarial action of artemisinin and its derivatives, although the source and nature of the heme remain controversial. Artemisinins also exhibit selective cytotoxicity against cancer cells in vitro and in vivo. We demonstrate that intracellular heme is the physiologically relevant mediator of the cytotoxic effects of artemisinins. Increasing intracellular heme synthesis through the addition of aminolevulinic acid, protoporphyrin IX, or transferrin-bound iron increased the cytotoxicity of dihydroartemisinin, while decreasing heme synthesis through the addition of succinyl acetone decreased its cytotoxic activity. A simple and robust high throughput assay was developed to screen chemical compounds that were capable of interacting with heme. A natural products library was screened which identified the compound coralyne, in addition to artemisinin, as a heme interacting compound with heme synthesis dependent cytotoxic activity. These results indicate that cellular heme may serve a general target for the development of both anti-parasitic and anti-cancer therapeutics.     

Comparison of the heme structures of horseradish peroxidase compounds X and II by resonance Raman spectroscopy

Horseradish peroxidase will catalyze the chlorination of certain substrates by sodium chlorite through an intermediate known as compound X. A chlorite-derived chlorine atom is known to be retained by compound X and has been proposed to be located at the heme active site. Although several heme structures have been proposed for compound X, including an Fe(IV)-OCl group, preliminary data previously reported by our laboratory suggested that compound X contained a heme Fe(IV) = O group, based on the similarity of a compound X resonance Raman band at 788 cm-1 to resonance Raman Fe(IV) = O stretching vibrations recently identified for horseradish peroxidase compound II and ferryl myoglobin. Isotopic studies now confirm that the 788 cm-1 resonance Raman band of compound X is, in fact, due to a heme Fe(IV) = O group, with the oxygen atom derived from chlorite. The Fe(IV) = O frequency of compound X, of horseradish peroxidase isoenzymes B and C, undergoes a pH-induced frequency shift, with behavior which appears to be the same as that previously reported for compound II, formed from the same isoenzymes. These observations strongly suggest that compounds II and X have very similar, if not identical, heme structures. The chlorine atom thus appears not to be heme-bound and may rather be located on an amino acid residue. The studies on compound X reported here were done in a pH region above pH 8, where compound X is moderately stable. The present results do not necessarily apply to compound X below pH 8.     

Using situs for flexible and rigid-body fitting of multiresolution single-molecule data

We describe here a set of multiresolution visualization and docking procedures that we refer to as the Situs package. The package was developed to provide an efficient and robust method for the fitting of atomic structures into low-resolution data. The current release was optimized specifically for the visualization and docking of single molecules. A novel 3D graphics viewer, volslice3d, was developed for the package to provide an immersive virtual reality environment for measuring and rendering volumetric data sets. The precision of single-molecule, rigid-body docking was tested on simulated (noise-free) low-resolution density maps. For spatial resolutions near 20 A typically arising in electron microscopy image reconstructions, a docking precision on the order of 1 A can be achieved. The shape-matching score captured the correct solutions in all 10 trial cases and was sufficiently stringent to yield unique matches in 8 systems. Novel routines were developed for the flexible docking of atomic structures whose shape deviates from the corresponding low-resolution shape. Test calculations on isoforms of actin and lactoferrin demonstrate that the flexible docking faithfully reproduces conformational differences with a precision < 2 A if atomic structures are locally conserved.     

Theoretical study on electronic structures of FeOO, FeOOH, FeO(H2O), and FeO in hemes: as intermediate models of dioxygen reduction in cytochrome c oxidase

The electronic structures of heme-dioxygen complexes have been studied as intermediate models of dioxygen reduction mechanism catalyzed by the mixed valence (MV) and fully reduced (FR) cytochrome c oxidase (CcO). Dioxygen, protons and electrons were sequentially added to the heme along the proposed reaction path for the O(2) reduction mechanism. The electronic structures of [FeOO], [FeOO](-), [FeOOH](+), [FeOOH], [Fe=O, H(2)O](+), [Fe=O](+) and [Fe=O] were thoroughly investigated by using the unrestricted hybrid exchange-correlation functional B3LYP method. The additions of two protons and an electron to [FeOO] lead to the OO bond cleavage to produce a H(2)O molecule with the electron transfer from Fe(II) in heme to the OO moiety. It is shown that the intrinsic OO bond cleavage occurs by adding two protons and two electrons into the OO bond, indicating consistency with a H(2)O formation catalyzed by both MV and FR CcO. The study of the electronic structures of heme-dioxygen complexes gave the different proposals for the mechanisms of a H(2)O formation by both MV and FR CcO. For the mixed valence CcO, starting from the [FeOO] complex, the final products are single H(2)O molecule and compound I of the oxo heme. For the fully reduced CcO, the final products are single H(2)O molecule and compound II of the oxo heme which is a reduced state of the compound I.     

Alkylation of human hemoglobin A0 by the antimalarial drug artemisinin

In vitro, the heme cofactor of human iron(II) hemoglobin was efficiently and quickly alkylated at meso positions by the peroxide-based antimalarial drug artemisinin, leading to heme-artemisinin-derived covalent adducts. This reaction occurred in the absence of any added protease or in the presence of an excess of an extra non-heme protein, or even when artemisinin was added to hemolysed human blood. This activation of artemisinin by the heme moiety of non-digested hemoglobin clearly indicates the high affinity of this drug for heme, and its efficient alkylating ability under very mild conditions.