Characterization of the structure and iron uptake mechanisms of the Staphylococcus aureus ferric hydroxamate-binding lipoprotein FhuD2

The Reverse Vaccinology (RV) approach allows using genomic information for the delineation of new protein-based vaccines starting from an in silico analysis. The first powerful example of the application of the RV approach is given by the development of a protein-based vaccine against serogroup B Meningococcus. A similar approach was also used to identify new Staphylococcus aureus vaccine candidates, including the ferric hydroxamate-binding lipoprotein FhuD2. S. aureus is a widespread human pathogen, which employs various different strategies for iron uptake, including: (i) siderophore-mediated iron acquisition using the endogenous siderophores staphyloferrin A and B, (ii) siderophore-mediated iron acquisition using xeno-siderophores (the pathway exploited by FhuD2) and (iii) heme-mediated iron acquisition. In this work the high resolution crystal structure of FhuD2 in the iron (III)-siderophore-bound form was determined. FhuD2 belongs to the Periplasmic Binding Protein family (PBP ) class III, and is principally formed by two globular domains, at the N- and C-termini of the protein, that make up a cleft where ferrichrome-iron (III) is bound. The N- and C-terminal domains, connected by a single long α-helix, present Rossmann-like folds, showing a β-stranded core and an α-helical periphery, which do not undergo extensive structural rearrangement when they interact with the ligand, typical of class III PBP members. The structure shows that ferrichrome-bound iron does not come directly into contact with the protein; rather, the metal ion is fully coordinated by six oxygen donors of the hydroxamate groups of three ornithine residues, which, with the three glycine residues, make up the peptide backbone of ferrichrome. Furthermore, it was found that iron-free ferrichrome is able to subtract iron from transferrin. This study shows for the first time the structure of FhuD2, which was found to bind to siderophores ,and that the protein plays an important role in S. aureus colonization and infection phases.

Ablation of the tumor suppressor histidine triad nucleotide binding protein 1 is protective against hepatic ischemia/reperfusion injury

The identification of cellular pathways capable of limiting ischemia/reperfusion (I/R) injury remains a frontier in medicine, and its clinical relevance is urgent. Histidine triad nucleotide binding protein 1 (HINT1) is a tumor suppressor that influences apoptosis. Because apoptotic pathways are a feature of I/R injury, we asked whether Hint1 influences hepatic I/R injury. Hint1(-/-) and C57BL/6 mice were subjected to 70% liver ischemia followed by reperfusion for 3 or 24 hours or to a sham operation. The serum aminotransferase levels, histological lesions, apoptosis, reactive oxygen species, and expression of B cell lymphoma 2-associated X protein (Bax), heme oxygenase 1 (HO-1), interleukin-6 (IL-6), IL-10, tumor necrosis factor-a, Src, nuclear factor kappa B (p65/RelA), and c-Jun were quantified. The responses to toll-like receptor ligands and nicotinamide adenine dinucleotide phosphate oxidase activity in Kupffer cells were compared in Hint1(-/-) mice and C57BL/6 mice. After I/R, the levels of serum aminotransferases, parenchymal necrosis, and hepatocellular apoptosis were significantly lower in Hint1(-/-) mice versus control mice. Furthermore, Bax expression decreased more than 2-fold in Hint1(-/-) mice, and the increases in reactive oxygen species and HO-1 expression that were evident in wild-type mice after I/R were absent in Hint1(-/-) mice. The phosphorylation of Src and the nuclear translocation of p65 were increased in Hint1(-/-) mice, whereas the nuclear expression of phosphorylated c-Jun was decreased. The levels of the protective cytokines IL-6 and IL-10 were increased in Hint1(-/-) mice. These effects increased survival after I/R in mice lacking Hint1. Hint1(-/-) Kupffer cells were less activated than control cells after stimulation with lipopolysaccharides. CONCLUSION: The Hint1 protein influences the course of I/R injury, and its ablation in Kupffer cells may limit the extent of the injury.

Study Of The Binding Of Substrate By The Phe557Val Mutant Soybean Lipoxygenase-1

Lipoxygenases are a class of enzymes which consist of non-heme iron dioxygenases that are produced by fungi, plants, and mammals and catalyze the oxygenation of polyunsaturated fatty acid substrates to unsaturated fatty acid hydroperoxide products. The unsaturated fatty acid hydroperoxide products are stereo- and regiospecific. One such lipoxygenase, soybean lipoxygenase-1 (SBLO-1), catalyzes the conversion of linoleate to 13-hydroperoxy-9(Z),11(E)-octadecadienoate (13-HPOD) and a small amount of 9-hydroperoxy-10(E),12(Z)-octadecadienoate (9-HPOD). Although the structure of SBLO-1 is known and it is the most widely studied lipoxygenase, how it binds to substrate is still poorly understood. Two competing binding hypotheses that have been used to understand and explain the binding are the head first binding model and the tail first binding model. The head first binding model predicts linoleate binds with its polar carboxylate group in the binding pocket and the methyl terminus at the surface of the binding pocket. The tail first binding model predicts that linoleate binds with its methyl terminus end in the binding pocket and the polar carboxylate group at the surface of the binding pocket. Both binding models have been used in the explanation of previous work. In previous work the replacement of phenylalanine with valine has been performed to produce the phe557val mutant SBLO-1. The mutant SBLO-1 was then used in the enzymatic oxygenation of linoleate. With this mutant, the amount of 9-HPOD that is formed increases. This result has been interpreted using the head-first binding model in which the smaller valine residue allows linoleate to bind with the polar carboxylate group of linoleate interacting with arginine-707. The work presented in this thesis confirms the regiochemical results of the previous work and further tests the head-first binding model. If head-first binding occurs, the 9-HPOD is expected to have primarily S configuration. Utilizing chiral-phase HPLC, it was found that the 9-HPOD produced by the phe557val mutant SBLO-1 is primarily S, consistent with head-first binding. The head-first binding model was also tested using linoleyl dimethylamine (LDMA), which has been shown to be a good substrate for SBLO-1 at pH 7.0, where LDMA is thought to be positively charged. This model predicts that less of the 9-peroxide should be produced with this substrate. Through the use of gas chromatography/mass spectrometry, it was found that the conversion of LDMA by the phe557val mutant SBLO-1 resulted in the formation of a 46:54 mixture of the 13-peroxide:9-peroxide. The higher amount of 9-peroxide is the opposite of what is expected for the currently proposed model suggesting that the proposed model may not be entirely correct. The results thus far have been consistent with reverse binding but not with the proposed interaction of the polar end of the substrate with arginine-707.

Structural studies of cytochrome P4501A1: Identification of substrate and cumene hydroperoxide binding regions in the active site of P4501A1 and characterization of the P4501A1 interaction with cytochrome P450 reductase

Cytochromes P450 are a superfamily of heme-thiolate proteins that function in a concert with another protein, cytochrome P450 reductase, as terminal oxidases of an enzymatic system catalyzing the metabolism of a variety of foreign compounds and endogenous substrates. In order to better understand P450s catalytic mechanism and substrate specificity, information about the structure of the active site is necessary. Given the lack of a crystal structure of mammalian P450, other methods have been used to elucidate the substrate recognition and binding site structure in the active center. In this project I utilized the photoaffinity labeling technique and site-directed mutagenesis approach to gain further structural insight into the active site of mammalian cytochrome P4501AI and examine the role of surface residues in the interaction of P4501A1 with the reductase. ^ Four crosslinked peptides were identified by photoaffinity labeling using diazido benzphetamine as a substrate analog. Alignment of the primary structure of cytochrome P4501A1 with that of bacterial cytochrome P450102 (the crystal structure of which is known) revealed that two of the isolated crosslinked peptides can be placed in the vicinity of heme (in the L helix region and β10-β11 sheet region of cytochrome P450102) and could be involved in substrate binding. The other two peptides were located on the surface of the protein with the label bound specifically to Lys residues that were proposed to be involved in reductase-P450 interaction. ^ Alternatively, it has been shown that some of the organic hydroperoxides can support P450 catalyzed reactions in the absence of NADPH, O2 and reductase. By means of photoaffinity labeling the cumene hydroperoxide binding region was identified. Using azidocumene as the photoaffinity label, the tripeptide T501-L502-K503 was shown to be the site where azidocumene covalently binds to P4501A1. The sequence alignment of cytochrome P4501A1 with cytochrome P450102 predicts that this region might correspond to β-sheet structure localized on the distal side of the heme ring near the I helix and the oxygen binding pocket. The role of Thr501 in the cumene hydroperoxide binding was confirmed by mutations of this residue and kinetic analysis of the effects of the mutations. ^ In addition, the role of two lysine residues, Lys271 and Lys279, in the interaction with reductase was examined by means of site-directed mutagenesis. The lysine residues were substituted with isoleucine and enzymatic activity of the wild type and the mutants were compared in reductase- and cumene hydroperoxide-supported systems. The lysine 279 residue has been shown to play a critical role in the P4501A1-reductase interaction. ^

Structural factors controlling ligand binding to myoglobin: A kinetic hole-burning study

Using temperature-derivative spectroscopy in the temperature range below 100 K, we have studied the dependence of the Soret band on the recombination barrier in sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation at 12 K. The spectra were separated into contributions from the photodissociated species, Mb*CO, and CO-bound myoglobin. The line shapes of the Soret bands of both photolyzed and liganded myoglobin were analyzed with a model that takes into account the homogeneous bandwidth, coupling of the electronic transition to vibrational modes, and static conformational heterogeneity. The analysis yields correlations between the activation enthalpy for rebinding and the model parameters that characterize the homogeneous subensembles within the conformationally heterogeneous ensemble. Such couplings between spectral and functional parameters arise when they both originate from a common structural coordinate. This effect is frequently denoted as “kinetic hole burning.” The study of these correlations gives direct insights into the structure–function relationship in proteins. On the basis of earlier work that assigned spectral parameters to geometric properties of the heme, the connections with the heme geometry are discussed. We show that two separate structural coordinates influence the Soret line shape, but only one of the two is coupled to the enthalpy barrier for rebinding. We give evidence that this coordinate, contrary to widespread belief, is not the iron displacement from the mean heme plane.

Heme is an effector molecule for iron-dependent degradation of the bacterial iron response regulator (Irr) protein

The bacterial iron response regulator (Irr) protein mediates iron-dependent regulation of heme biosynthesis. Pulse–chase and immunoprecipitation experiments showed that Irr degraded in response to 6 μM iron with a half-life of ≈30 min and that this regulated stability was the principal determinant of control by iron. Irr contains a heme regulatory motif (HRM) near its amino terminus. A role for heme in regulation was implicated by the retention of Irr in heme synthesis mutants in the presence of iron. Addition of heme to low iron (0.3 μM) cultures was sufficient for the disappearance of Irr in cells of the wild-type and heme mutant strains. Spectral and binding analyses of purified recombinant Irr showed that the protein bound heme with high affinity and caused a blue shift in the absorption spectrum of heme to a shorter wavelength. A Cys29 → Ala substitution within the HRM of Irr (IrrC29A) abrogated both high affinity binding to heme and the spectral blue shift. In vivo turnover experiments showed that, unlike wild-type Irr, IrrC29A was stable in the presence of iron. We conclude that iron-dependent degradation of Irr involves direct binding of heme to the protein at the HRM. The findings implicate a regulatory role for heme in protein degradation and provide direct evidence for a functional HRM in a prokaryote.

Femtosecond resolution of ligand-heme interactions in the high-affinity quinol oxidase bd: A di-heme active site?

Interaction of the two high-spin hemes in the oxygen reduction site of the bd-type quinol oxidase from Escherichia coli has been studied by femtosecond multicolor transient absorption spectroscopy. The previously unidentified Soret band of ferrous heme b595 was determined to be centered around 440 nm by selective excitation of the fully reduced unliganded or CO-bound cytochrome bd in the α-band of heme b595. The redox state of the b-type hemes strongly affects both the line shape and the kinetics of the absorption changes induced by photodissociation of CO from heme d. In the reduced enzyme, CO photodissociation from heme d perturbs the spectrum of ferrous cytochrome b595 within a few ps, pointing to a direct interaction between hemes b595 and d. Whereas in the reduced enzyme no heme d-CO geminate recombination is observed, in the mixed-valence CO-liganded complex with heme b595 initially oxidized, a significant part of photodissociated CO does not leave the protein and recombines with heme d within a few hundred ps. This caging effect may indicate that ferrous heme b595 provides a transient binding site for carbon monoxide within one of the routes by which the dissociated ligand leaves the protein. Taken together, the data indicate physical proximity of the hemes d and b595 and corroborate the possibility of a functional cooperation between the two hemes in the dioxygen-reducing center of cytochrome bd.

Relative role of heme nitrosylation and β-cysteine 93 nitrosation in the transport and metabolism of nitric oxide by hemoglobin in the human circulation

To quantify the reactions of nitric oxide (NO) with hemoglobin under physiological conditions and to test models of NO transport on hemoglobin, we have developed an assay to measure NO–hemoglobin reaction products in normal volunteers, under basal conditions and during NO inhalation. NO inhalation markedly raised total nitrosylated hemoglobin levels, with a significant arterial–venous gradient, supporting a role for hemoglobin in the transport and delivery of NO. The predominant species accounting for this arterial–venous gradient is nitrosyl(heme)hemoglobin. NO breathing increases S-nitrosation of hemoglobin β-chain cysteine 93, however only to a fraction of the level of nitrosyl(heme)hemoglobin and without a detectable arterial–venous gradient. A strong correlation between methemoglobin and plasma nitrate formation was observed, suggesting that NO metabolism is a primary physiological cause of hemoglobin oxidation. Our results demonstrate that NO–heme reaction pathways predominate in vivo, NO binding to heme groups is a rapidly reversible process, and S-nitrosohemoglobin formation is probably not a primary transport mechanism for NO but may facilitate NO release from heme.

Inducible nitric oxide synthase: identification of amino acid residues essential for dimerization and binding of tetrahydrobiopterin.

Nitric oxide synthases (NOSs) require tetrahydrobiopterin (BH4) for dimerization and NO production. Mutation analysis of mouse inducible NOS (iNOS; NOS2) identified Gly-450 and Ala-453 as critical for NO production, dimer formation, and BH4 binding. Substitutions at five neighboring positions were tolerated, and normal binding of heme, calmodulin, and NADPH militated against major distortions affecting the NH2-terminal portion, midzone, or COOH terminus of the inactive mutants. Direct involvement of residues 450 and 453 in the binding of BH4 is supported by the striking homology of residues 448-480 to a region extensively shared by the three BH4-utilizing aromatic amino acid hydroxylases and is consistent with the conservation of these residues among all 10 reported NOS sequences, including mammalian NOSs 1, 2, and 3, as well as avian and insect NOSs. Altered binding of BH4 and/or L-arginine may explain how the addition of a single methyl group to the side chain of residue 450 or the addition of three methylenes to residue 453 can each abolish an enzymatic activity that reflects the concerted function of 1143 other residues.

Role of heme regulated eIF2α kinase in pancreatic beta cell function and viability

The eukaryotic translation initiation factor 2 alpha (eIF2α) is part of the initiation complex that drives the initiator amino acid methionine to the ribosome, a crucial step in protein translation. In stress conditions such as virus infection, endoplasmic reticulum (ER) stress, amino acid or heme deficiency eIF2α can be phosphorylated and thereby inhibit global protein synthesis. This adaptive mechanism prevents protein accumulation and consequent cytotoxic effects. Heme-regulated eIF2α kinase (HRI) is a member of the eIF2α kinase family that regulates protein translation in heme deficiency conditions. Although present in all tissues, HRI is predominantly expressed in erythroid cells where it remains inactive in the presence of normal heme concentrations. In response to heme deficiency, HRI is activated and phosphorylates eIF2α decreasing globin synthesis. This mechanism is important to prevent accumulation of heme-free globin chains which cause ER stress and apoptosis. RNA sequencing data from our group showed that in human islets and in primary rat beta cells HRI is the most expressed eIF2α kinase compared to the other family members. Despite its high expression levels, little is known about HRI function in beta cells. The aim of this project is to identify the role of HRI in pancreatic beta cells. This was investigated taking a loss-of-function approach. HRI knock down (KD) by RNA interference induced beta cell apoptosis in basal condition. HRI KD potentiated the apoptotic effects of palmitate or proinflammatory cytokines, two in vitro models for type 2 and type 1 diabetes, respectively. Increased cytokine-induced apoptosis was also observed in HRI-deficient primary rat beta cells. Unexpectedly, we observed a mild increase in eIF2α phosphorylation in HRI-deficient cells. The levels of mRNA or protein expression of C/EBP homologous protein (CHOP) and activating transcription factor 4 (ATF4) were not modified. HRI KD cells have decreased spliced X-box binding protein 1 (XBP1s), an important branch of the ER stress response. However, overexpression of XBP1s by adenovirus in HRI KD cells did not protect from HRI siRNA-induced apoptosis. HRI deficiency decreased phosphorylation of Akt and its downstream targets glycogen synthase kinase 3 (GSK3), forkhead box protein O1 (FOXO1) and Bcl-2-associated death promoter (BAD). Overexpression of a constitutively active form of Akt by adenovirus in HRI-deficient beta cells partially decreased HRI KD-mediated apoptosis. Interestingly, BAD silencing protected from apoptosis caused by HRI deficiency. HRI silencing in beta cells also induced JNK activation. These results suggest an important role of HRI in beta cell survival through modulation of the Akt/BAD pathway. Thus, HRI may be an interesting target to modulate beta cell fate in diabetic conditions.

Kidney dysfunction and hypertension: Role for cadmium, P450 and heme oxygenases?

Cadmium (Cd) is a metal toxin of continuing worldwide concern. Daily intake of Cd, albeit in small quantities, is associated with a number of adverse health effects which are attributable to distinct pathological changes in a variety of tissues and organs. In the present review, we focus on its renal tubular effects in people who have been exposed environmentally to Cd at levels below the provisional tolerable intake level set for the toxin. We highlight the data linking such low-level Cd intake with tubular injury, altered abundance of cytochromes P450 (CYPs) in the kidney and an expression of a hypertensive phenotype. We provide updated knowledge on renal and vascular effects of the eicosanoids 20-hydroxyeicosatetraenoic acid (20-HETE) and eicosatrienoic acids (EETs), which are biologically active metabolites from arachidonate metabolism mediated by certain CYPs in the kidney. We note the ability of Cd to elicit oxidative stress and to alter metal homeostasis notably of zinc which may lead to augmentation of the defense mechanisms involving induction of the antioxidant enzyme heme oxygenase-1 (HO-1) and the metal binding protein metallothionein (MT) in the kidney. We hypothesize that renal Cd accumulation triggers the host responses mediated by HO-I and MT in an attempt to protect the kidney against injurious oxidative stress and to resist a rise in blood pressure levels. This hypothesis predicts that individuals with less active HO-1 (caused by the HO-1 genetic polymorphisms) are more likely to have renal injury and express a hypertensive phenotype following chronic ingestion of low-level Cd, compared with those having more active HO-1. Future analytical and molecular epidemiologic research should pave the way to the utility of induction of heme oxygenases together with dietary antioxidants in reducing the risk of kidney injury and hypertension in susceptible people.

Proton nuclear magnetic resonance investigation of the native and modified active site structure of heme proteins

Hemoproteins are a very important class of enzymes in nature sharing the essentially same prosthetic group, heme, and are good models for exploring the relationship between protein structure and function. Three important hemoproteins, chloroperoxidase (CPO), horseradish peroxidase (HRP), and cytochrome P450cam (P450cam), have been extensively studied as archetypes for the relationship between structure and function. In this study, a series of 1D and 2D NMR experiments were successfully conducted to contribute to the structural studies of these hemoproteins. ^ During the epoxidation of allylbenzene, CPO is converted to an inactive green species with the prosthetic heme modified by addition of the alkene plus an oxygen atom forming a five-membered chelate ring. Complete assignment of the NMR resonances of the modified porphyrin extracted and demetallated from green CPO unambiguously established the structure of this porphyrin as an NIII-alkylated product. A novel substrate binding motif of CPO was proposed from this concluded regiospecific N-alkylation structure. ^ Soybean peroxidase (SBP) is considered as a more stable, more abundant and less expensive substitute of HRP for industrial applications. A NMR study of SBP using 1D and 2D NOE methods successfully established the active site structure of SBP and consequently fills in the blank of the SBP NMR study. All of the hyperfine shifts of the SBP-CN- complex are unambiguously assigned together with most of the prosthetic heme and all proximal His170 resonances identified. The active site structure of SBP revealed by this NMR study is in complete agreement with the recombinant SBP crystal structure and is highly similar to that of the HRP with minor differences. ^ The NMR study of paramagnetic P450cam had been greatly restricted for a long time. A combination of 2D NMR methods was used in this study for P450cam-CN - complexes with and without camphor bound. The results lead to the first unequivocal assignments of all heme hyperfine-shifted signals, together with certain correlated diamagnetic resonances. The observed alternation of the assigned novel proximal cysteine β-CH2 resonances induced by camphor binding indicated a conformational change near the proximal side.^

Conformational dynamics associated with ligand binding to vertebrate hexa-coordinate hemoglobins

Neuroglobin (Ngb) and cytoglobin (Cygb) are two new additions to the globin family, exhibiting heme iron hexa-coordination, a disulfide bond and large internal cavities. These proteins are implicated in cytoprotection under hypoxic-ischemic conditions, but the molecular basis of their cytoprotective function is unclear. Herein, a photothermal and spectroscopic study of the interactions of diatomic ligands with Ngb, Cygb, myoglobin and hemoglobin is presented. The impact of the disulfide bond in Ngb and Cygb and role of conserved residues in Ngb His64, Val68, Cys55, Cys120 and Tyr44 on conformational dynamics associated with ligand binding/dissociation were investigated. Transient absorption and photoacoustic calorimetry studies indicate that CO photo-dissociation from Ngb leads to a volume expansion (13.4±0.9 mL mol-1), whereas a smaller volume change was determined for Ngb with reduced Cys (ΔV=4.6±0.3 mL mol-1). Furthermore, Val68 side chain regulates ligand migration between the distal pocket and internal hydrophobic cavities since Val68Phe geminate quantum yield is ∼2.7 times larger than that of WT Ngb. His64Gln and Tyr44Phe mutations alter the thermodynamic parameters associated with CO photo-release indicating that electrostatic/hydrogen binding network that includes heme propionate groups, Lys 67, His64, and Tyr 44 in Ngb modulates the energetics of CO photo-dissociation. In Cygb, CO escape from the protein matrix is fast (< 40 ns) with a ΔH of 18±2 kcal mol-1 in Cygbred, whereas disulfide bridge formation promotes a biphasic ligand escape associated with an overall enthalpy change of 9±4 kcal mol-1. Therefore, the disulfide bond modulates conformational dynamics in Ngb and Cygb. I propose that in Cygb with reduced Cys the photo-dissociated ligand escapes through the hydrophobic tunnel as occurs in Ngb, whereas the CO preferentially migrates through the His64 gate in Cygbox. To characterize Cygb surface 1,8-ANS interactions with Cygb were investigated employing fluorescence spectroscopy, ITC and docking simulations. Two 1,8-ANS binding sites were identified. One binding site is located close to the extended N-terminus of Cygb and was also identified as a binding site for oleate. Furthermore, guanidinium hydrochloride-induced unfolding studies of Cygb reveal that the disulfide bond does not impact Cygb stability, whereas binding of cyanide slightly increases the protein stability.

Myoglobins: the link between discoloration and lipid oxidation in muscle and meat

Aerobic metabolism changes rapidly to glycolysis post-mortem resulting in a pH-decrease during the transformation of muscle in to meat affecting ligand binding and redox potential of the heme iron in myoglobin, the meat pigment. The inorganic chemistry of meat involves (i) redox-cycling between iron(II), iron(III), and iron(IV)/protein radicals; (ii) ligand exchange processes; and (iii) spin-equilibra with a change in coordination number for the heme iron. In addition to the function of myoglobin for oxygen storage, new physiological roles of myoglobin are currently being discovered, which notably find close parallels in the processes in fresh meat and nitrite-cured meat products. Myoglobin may be characterized as a bioreactor for small molecules like O2, NO, CO, CO2, H2O, and HNO with importance in bio-regulation and in protection against oxidative stress in vivo otherwise affecting lipids in membranes. Many of these processes may be recognised as colour changes in fresh meat and cured meat products under different atmospheric conditions, and could also be instructive for teaching purposes.

Estudo teórico da interação existente entre a artemisinina e o heme

The purpose of this work is to study theoretically stereoelectronic aspects of the interaction between heme and artemisinin in the transitional heme-artemisinin complex. Through semi-empirical calculations using the PM3 method, the potential energy barrier of artemisinin rotation relative to heme in the heme-artemisinin complex was studied in vacuum and in the partially solvated state. The minimum heat of formation obtained for the complex with and without water molecules is -702.39 and -100.86 kcal mol-1, respectively, which corresponds to the dihedral angle C-Fe-O1-O2 of 43.93º and 51.90º around the iron-oxygen O1 bond, respectively. The water molecules bind to heme via 13 hydrogen bonds and O-H O and 6 C-H O interactions, which accounts for -67.23 kcal mol-1. It is observed that the inclusion of water molecules does not affect significantly the stability of the heme-artemisinin complex.