Abstract

Malaria continues to be a major cause of children's morbidity and mortality worldwide, causing nearly one million deaths annually. The human malaria parasite, Plasmodium falciparum, synthesizes fatty acids employing the Type II fatty acid biosynthesis system (FAS II), unlike humans that rely on the Type I (FAS I) pathway. The FAS II system elongates acyl fatty acid precursors of the cell membrane in Plasmodium. Enoyl reductase (ENR) enzyme is a member of the FAS II system. Here we present steady-state kinetics, pre-steady-state kinetics, and equilibrium fluorescence spectroscopy data that allowed proposal of P. falciparum ENR (PfENR) enzyme mechanism. Moreover, building on previous results, the present study also evaluates the PfENR inhibition by the pentacyano(isoniazid)ferrateII compound. This inorganic complex represents a new class of lead compounds for the development of antimalarial agents focused on the inhibition of PfENR.

1. Introduction

Malaria, a disease caused by a protozoan pathogen of the genus Plasmodium, remains one of the most devastating human diseases, causing as many as 250 million new cases and nearly one million deaths annually worldwide. This disease is especially worrisome in Africa, where one out of every five infant deaths is related to its effects [1]. Although malaria can be caused by mainly four species of the genus, Plasmodium falciparum is responsible for the most severe and deadly form of the disease and is to blame for 90% of malaria-related deaths occurring in African children below the age of five [2]. At present, the available drugs used to fight the disease are (a) aryl aminoalcohol compounds, such as quinine [3], (b) antifolates-dihydrofolate reductase inhibitors, such as pyrimethamine [4, 5], or (c) artemisinin derivatives [6]. Therapeutic treatment has, however, encountered serious obstacles, as this parasite has developed at least partial resistance to nearly every antimalarial regimen introduced to date, including the first-line drugs chloroquine and sulfadoxine-pyrimethamine [7]. As single-drug treatments are no longer adequate, combinations of two or more drugs will probably offer improved efficacy and reduced risk of emergence of drug-resistant parasites. Current examples of drug combinations are artemisinin-amodiaquine and artemether-lumefantrine (Coartem) [8]. The lack of an effective prophylactic vaccine and the generation and dissemination of drug-resistant strains clearly indicate the need for the development of new therapeutic approaches and identification of novel targets.

P. falciparum, a malarial parasite of the phylum Apicomplexa, has been found to contain an apicoplast, an organelle that originally arose from a cyanobacterium through a secondary endosymbiotic process, and, thus, possesses four membranes. During these endosymbiotic processes, resultant of millions of years of adaptive selection, the endosymbiont and host were integrated in a relationship involving gene losses and gene transfers, resulting in a reduction of the apicoplast genome to a bare minimum [9]. The majority of the apicoplast proteins are encoded by genes in the nuclear genome of the host, and the proteins that reside in the apicoplast and are encoded by the nucleus are posttranslationally targeted to the apicoplast organelle via the secretory pathway, owing to the existence of a bipartite N-terminal leader sequence [10]. In fact, the apicoplast is indispensable for the malarial parasite, since several vital metabolic processes for the parasite occur in this organelle. Among these processes are isoprene, haem, and fatty acid biosynthesis [11]. Fatty acid biosynthesis is critical for the parasite development as fatty acids are the major components of cell membranes. Moreover, fatty acids are important source of energy, play a key role in signal transduction as well as in protein acylation, and are needed for the growth, differentiation, and homeostasis in P. falciparum. It is also known that lipid biosynthesis is highly increased during the erythrocytic phases, when the parasite grows and divides extremely fast [12]. As the parasite invades its host, it attempts to protect itself by creating a so-called parasitophorous vacuole, in part as a protection from the host’s immune system [13]. In this process, the parasite needs to produce its own fatty acids de novo so as to form and expand its membrane [14].

The enzymes involved in the biosynthesis of fatty acids are organized in two distinct ways within living systems. Fungi, mammals, and some mycobacteria accomplish fatty acid synthesis through the action of multifunctional proteins, each reaction being catalysed by a distinct region (domain) of these single-polypeptide proteins. These enzymes are classified as Type I fatty acid synthases (FAS I) and are also described as the “associative” type of fatty acid synthases, as successive steps in the fatty acid synthetic reaction occur at specific domains [15, 16]. In contrast, plants and many bacteria use the Type II or “dissociated” fatty acid synthases (FAS II), whose enzymes are best characterized in Escherichia coli [17]. Each step of the FAS II system is catalysed by a specific enzyme encoded by a separate gene. Since the FAS II is absent in humans, it represents a great opportunity for therapeutic applications. The fatty acid chain elongation module of FAS II consists of the cycling of the growing fatty acyl moiety through four key enzymes. Elongation is initiated by β-ketoacyl-ACP synthase I/II (FabB/FabF) enzyme, which catalyses the condensation of the growing fatty acyl ACP with malonyl-ACP to yield a β-ketoacyl-ACP with an additional two-carbon unit. The second step is the reduction of β-ketoacyl-ACP to β-hydroxyacyl-ACP catalysed by the β-ketoacyl-ACP reductase (FabG) enzyme. Water is then removed by β-hydroxydecanoyl-ACP dehydratase/isomerase (FabA) and/or β-hydroxyacyl-ACP dehydratase (FabZ) to produce 2-trans-enoyl-ACP. The final step in the elongation process is the NADH-dependent reduction of the double bond to produce fatty acyl ACP, which is catalysed by FabI, the only enoyl-ACP reductase (ENR) enzyme in Plasmodium (Figure 1). Successive elongation cycles typically lead to the formation of decanoic (C-10), lauric (C-12), and myristic (C-14) acids as the major fatty acyl ACPs in Plasmodium [13], whereas palmitic (C-16) acid is the most predominant fatty acid synthesized in humans [12].

ENR from Mycobacterium tuberculosis (MtInhA) was first identified as being the primary target for isoniazid (INH; Figure 2(a)) [18], the most prescribed drug to treat active tuberculosis. The enzyme is an NADH-dependent enoyl-ACP reductase specific for long-chain enoyl thioester substrates. MtInhA is a member of the mycobacterial FAS II, which elongates acyl fatty acid precursors of mycolic acids. INH is a prodrug that is activated by the de mycobacterial katG-encoded catalase-peroxidase enzyme in the presence of manganese ions, NAD (H) and oxygen [19, 20]. Approximately 50% of M. tuberculosis INH-resistant clinical isolates harbor deletions of, or missense mutations in, the katG gene [21, 22]. The KatG-produced acytylpyridine fragment of INH is covalently attached to the C4 position of NADH and forms a binary complex with the wild-type MtInhA [23]. This isonicotinoyl-NAD+ adduct has been characterized spectroscopically and shown to bind to the enzyme with a dissociation constant value lower than 0.4 nM [24]. The isonicotinoyl-NAD+ adduct has been shown to be a slow, tight-binding competitive inhibitor of wild-type MtInhA [25]. The initial rapidly reversible weak binding (𝐾𝑖=16 nM) is followed by a slow isomerization, leading to a tighter enzyme-inhibitor complex with an overall dissociation constant (𝐾𝑖) value of 0.75 nM.

The mechanism of action of INH requires its conversion, by the mycobacterial KatG, into a number of electrophilic intermediates [20]. The formation of the adduct has been proposed to be formed through the addition of either an isonicotinic acyl anion to NAD+ or an isonicotinic acyl radical to an NAD free radical [23]. The peroxidase activity of KatG catalyses the conversion of Mn2+ to Mn3+ [26]. Rapid freeze-quench electron paramagnetic resonance spectroscopy experiments have shown that hydrogen atom abstractions may be initiated by a KatG tyrosil radical [27]. Furthermore, it has been shown that the yield of isonicotinoyl-NAD+ adduct is about the same after oxidation of INH by KatG or Mn3+ [28]. Accordingly, oxidation by Mn3+-pyrophosphate has been proposed as an alternative method for nonenzymatic INH activation for simple chemical synthesis of various INH derivatives that mimic the isonicotinoyl-NAD+ adduct [29]. Based on the activation mechanism proposed for INH, via an electron transfer reaction, an alternative self-activation inner-sphere electron transfer reaction route was proposed for designing new drugs for the treatment of INH-resistant and INH-sensitive tuberculosis. These drugs would be activated by electron transfer reactions before interacting with their cellular targets, and a redox reversible metal complex coordinated to the prodrug was employed as a prototype. Accordingly, we have recently demonstrated that pentacyano(isoniazid)ferrateII or [FeII(CN)5(INH)]3− (Figure 2(b)) inhibits the activity of both wild-type and I21V mutant (INH-resistant) MtInhA enzymes [30]. The in vitro kinetics of inactivation indicates that this process requires no activation by KatG, no need for the presence of NADH, and is also effective against INH-resistant InhA mutants. An MIC value of 0.2 μg mL−1 for this inorganic complex was determined by the radiometric BACTEC AFB system for the M. tuberculosis H37Rv strain, and toxicity assays in HL60 leukemia and MCS-7 breast cancer cells yielded an IC50 value greater than 25 μg mL−1, thereby indicating a good selectivity index (SI = IC50/MIC > 125; the Tuberculosis Antimicrobial Acquisition and Coordinating Facility of USA suggests that in order to a compound move forward through screening programs, SI should be larger than 10). More recently, we have shown that the [FeII(CN)5(INH)]3− complex is a slow-onset inhibitor of MtInhA enzyme activity, with a true overall dissociation constant value of 70 nM [31]. In this mechanism of action an initial enzyme-inhibitor complex is rapidly formed, which then undergoes a slow isomerization reaction of an enzyme-inhibitor binary complex, where the inhibitor is more tightly bound to enzymes. The weakness in the use of classical enzyme inhibitors as drugs for clinical conditions is that inhibition results in the upstream accumulation of the substrate for the enzyme, which may overcome the inhibition. By contrast, the build-up of substrate cannot have any effect on the isomerization of the enzyme-inhibitor complex typical of the slow-onset mechanism and hence reversal of the inhibition [32].

The structural similarity between MtInhA and PfENR stimulated us to test the [FeII(CN)5(INH)]3− against the malarial enzyme. This work investigates the efficacy of the [FeII(CN)5(INH)]3− compound against PfENR. In addition, we propose an enzyme mechanism for PfENR based on steady-state kinetics, pre-steady-state kinetics, and fluorescence spectroscopy data.

2. Materials and Methods

2.1. Materials

All chemicals used were of analytical or reagent grade and required no further purification. Crotonoyl-CoA, NADH, and butyryl-CoA were from Sigma (Steinheim, Germany). PfENR (EC 1.3.1.9) was expressed and purified by Genscript Corp., New Jersey, USA. The [FeII(CN)5(INH)]3− compound was a kind gift of Dr. Luiz G.F Lopes and Dr. Eduardo H.S. Sousa from Universidade Federal do Ceará, Departamento de Química Orgânica e Inorgânica, Fortaleza, CE, Brazil. The steady-state activity assays were carried out in a UV-2550 UV-Visible Spectrophotometer, and fluorescence binding experiments were carried out in an RF-5301PC Spectrophotometer, both from Shimadzu (Kyoto, Japan).

All enzyme activity assays were performed under initial rate conditions at 25°C in 100 mM sodium phosphate buffer pH 7.5 in 500 μL total reaction volumes, and each individual initial rate datum was the average of duplicate or triplicate measurements. Crotonoyl-CoA, NADH, butyryl-CoA, [FeII(CN)5(INH)]3−, and the enzyme were dissolved in 100 mM sodium phosphate buffer pH 7.5. The reactions were monitored by the decrease in absorbance at 340 nm due to the conversion of NADH (ε340 = 6220 M−1 cm−1) to NAD+.

2.2. Initial Velocity

In order to determine the true steady-state kinetic parameters and initial velocity patterns, PfENR activity was measured in the presence of varying concentrations of crotonoyl-CoA (10–240 μM) and several fixed-varied NADH concentrations (20–200 μM), and reactions were initiated by the addition of 0.12 μM of PfENR into the reaction mixture.

2.3. Fluorescence Spectroscopy

Fluorescence titration was performed to assess the dissociation constant at equilibrium for binary complexes formed between PfENR and either substrate(s) or product(s) at 25°C. Crotonoyl-CoA, NADH, butyryl-CoA, and PfENR were dissolved in 100 mM sodium phosphate buffer pH 7.5. Fluorescence titration with crotonoyl-CoA was carried out by making microliter additions of 0.5 mM and 5 mM crotonoyl-CoA stock solutions (0.50–425.97 μM crotonoyl-CoA final concentrations) to 1 mL of 3 μM PfENR, keeping the dilution to a maximum of 10.1%. Measurements of changes in intrinsic protein fluorescence of PfENR upon crotonoyl-CoA binding employed excitation wavelength at 305 nm, and the emission wavelengths ranged from 325 to 370 nm (maximum PfENR 𝜆EM=348 nm); the slits for excitation and emission wavelengths were, respectively, 15 and 5 nm. Fluorescence titration with NADH was carried out by making microliter additions of 0.5 mM, 1 mM, and 5 mM NADH stock solutions (0.50–119.21 μM NADH final concentrations) to 1 mL of 3 μM PfENR, keeping the dilution to a maximum of 4.1%. Measurements of changes in NADH fluorescence upon binding to PfENR employed excitation wavelength at 360 nm, and the emission wavelengths ranged from 360 to 560 nm (maximum 𝜆EM=460 nm); the slits for excitation and emission wavelengths were, respectively, 1.5 and 5 nm. Fluorescence titration with butyryl-CoA was carried out by making microliter additions of 0.5 mM, 2.5 mM, 5 mM, and 10 mM butyryl-CoA stock solutions (0.50–550.85 μM butyryl-CoA final concentrations) to 1 mL of 3 μM PfENR, keeping the dilution to a maximum of 8.4%. Measurements of changes in intrinsic protein fluorescence upon butyryl-CoA binding to PfENR employed excitation wavelength at 300 nm, and the emission wavelengths ranged from 320 to 400 nm (maximum PfENR 𝜆EM=337 nm); the slits for excitation and emission wavelengths were, respectively, 5 and 10 nm. Control experiments were employed to both determine the maximum ligand concentrations to be used with no inner filter effect and determine dilution effects on protein fluorescence. No fluorescence titration with NAD+ could be performed due to a large inner filter on protein fluorescence.

2.4. Pre-Steady-State Kinetics

Pre-steady-state kinetic measurements of NADH→NAD+ conversion by PfENR were carried out to confirm the proposed kinetic mechanism and assess whether the apparent rate constant for this conversion contributes to the rate-limiting step. The decrease in absorbance upon NADH conversion to NAD+ was monitored at 340 nm (1 mm slit width = 4.65 nm spectral band), at 25°C, using a split time base (0.2–2 s; 200 data points for each time base). The experimental conditions were 3.14 μM PfENR and 200 μM NADH in 100 mM sodium phosphate buffer pH 7.5 (mixing chamber concentrations). The dead time of the stopped-flow equipment is 1.37 ms.

2.5. [FeII(CN)5(INH)]3− In Vitro Activity against PfENR

[FeII(CN)5(INH)]3− time-dependent inhibition of PfENR was performed since it had been observed in MtInhA [31]. The enzyme was preincubated with the inhibitor for 15 min (2.77 : 2455 μM enzyme : inhibitor solution), and the activity was monitored by adding 0.12 : 54 μM enzyme : inhibitor mixture aliquots, with different times of exposure, to the reaction mixtures which contained crotonoyl-CoA and NADH concentrations close to their 𝐾M values. The reaction was monitored by following enzyme activity during the time course.

2.6. Slow-Binding Inhibition Kinetics of PfENR by [FeII(CN)5·(INH)]3−

Typical slow-binding experiments were carried out in the presence of various [FeII(CN)5(INH)]3− concentrations as an attempt to elucidate the kinetic mechanism of the time-dependent inhibition observed for PfENR. In the case of slow-binding inhibitors, the binding, which is comprised by the establishment of the equilibrium between enzyme, inhibitor, and enzyme-inhibitor complexes, occurs slowly in a time scale of seconds to minutes [32]. These experiments were conducted under conditions in which [FeII(CN)5(INH)]3− concentrations were significantly larger than enzyme concentration, and measurements were made within a time range over which there was no significant substrate depletion. The assay mixtures contained 80 μM crotonoyl-CoA, 120 μM NADH, 2% (vol/vol) glycerol, 49.7 μg mL−1 BSA, and 0–35 μM [FeII(CN)5(INH)]3−. The reactions were initiated by adding enzyme (0.69 nM) to the assay mixture and were monitored by following the decrease in NADH absorbance at 340 nm for 21600 s (6 hours). Enzyme activity was stabilized by the addition of glycerol and BSA to the reaction mixture as described for MtInhA [25, 31].

2.7. Two-Step Mechanism of PfENR Inhibition by [FeII(CN)5·(INH)]3−

The two-step mechanism observed for MtInhA inhibition by [FeII(CN)5(INH)]3− [31] was investigated for PfENR. The latter was preincubated in the presence of various concentrations of [FeII(CN)5(INH)]3−, and the decrease of PfENR activity was assessed by withdrawing aliquots from the preincubation mixtures at appropriate intervals, followed by rapid measurement of initial velocity. Reactions of time-dependent inactivation of PfENR were carried out in mixtures containing 12 μM PfENR and various concentrations of inhibitor. Control experiments were performed with no inhibitor. At various times, the residual enzyme activity was determined by diluting an aliquot (100-fold) of the reaction mixture into the assay solution, which contained 40 μM crotonoyl-CoA and 50 μM NADH.

2.8. Data Analysis

The kinetic parameter values and their respective standard errors were obtained by fitting the data to the appropriate equations by using the nonlinear regression function of SigmaPlot 2004 (SPSS, Inc.).

Apparently parallel family of lines of initial velocity data were tentatively fitted to (1), which describes a double displacement mechanism [33]. For (1), 𝑣 is the measured reaction velocity, 𝑉 is the maximal velocity, 𝐴 and 𝐵 are the concentrations of the substrates (NADH and crotonoyl-CoA), and 𝐾𝑎 and 𝐾𝑏 are their respective Michaelis constants:𝑣=𝑉𝐴𝐵𝐾𝑎𝐵+𝐾𝑏𝐴+𝐴𝐵.(1)

Data from equilibrium fluorescence spectroscopy were fitted either to a hyperbolic function ((2); [34]) or to the Hill equation ((3); [35]). For (2) and (3), 𝐹 is the observed fluorescence signal, 𝐹max is the maximal fluorescence, 𝐹/𝐹max ratio represents the degree of saturation, 𝐴 is the concentration of the ligand (either NADH or butyryl-CoA), 𝑛 represents the total number of binding sites, 𝐾𝑑 is the overall dissociation constant, and 𝐾 represents a mean dissociation constant for enzyme: ligand binary complex formation, which is comprised of interaction factors and the intrinsic dissociation constant [36, 37]:𝐹𝐹max=𝐴𝐾𝑑𝐹+𝐴,(2)𝐹max=𝐴𝑛𝐾+𝐴𝑛.(3)

The equilibrium fluorescence data representing the fractional saturation (0.10–0.90) for the homotetrameric form [38] of PfENR were fitted to (4) (the Hill plot), in which 𝑌 is the fraction of substrate binding sites occupied by the substrate and is the Hill coefficient:𝑌log1𝑌=log𝐴log𝐾.(4)

The pre-steady-state time course of NADH→NAD+ conversion was fitted to (5) for a single exponential decay, in which 𝐴 is the absorbance at time 𝑡, 𝐴0 is the absorbance at time zero, and 𝑘 is the apparent first-order rate constant for product formation:𝐴=𝐴0𝑒𝑘𝑡.(5)

Apparent pseudofirst-order rate constants were calculated by fitting the progress curves for time-dependent inactivation of PfENR to integrated equation (5) for slow-binding inhibition [25, 32, 39]. For (6), 𝐴𝑡 and 𝐴0 are the 340-nm absorbances at time 𝑡 and time 0, 𝑣0, 𝑣𝑠, and 𝑘obs represent, respectively, the initial velocity, the final steady-state velocity, and the apparent first-order rate constant for the establishment of the equilibrium between EI and EI* complexes (the observed rate constant for formation of EI* or the apparent isomerization rate constant):𝐴𝑡=𝑣𝑠𝑣𝑡+0𝑣𝑠1𝑒𝑘obs𝑡𝑘obs.(6)

An estimate of the apparent dissociation constant (𝐾app𝑖) was obtained by fitting the variation of initial rate (𝑣0) with inhibitor concentration (𝐼) to (7) for competitive inhibition, in which 𝑣0max is the initial velocity in the absence of inhibitor [31]:𝑣0=𝑣0max[𝐼]1+/𝐾app𝑖.(7)

A “true” value for the apparent dissociation constant 𝐾𝑖 for the initial enzyme-inhibitor binary complex (EI) can be calculated by (8), in which 𝐾NADH and 𝐾Crotonoyl-CoA are the values for the Michaelis-Menten constants and [NADH] and [Crotonoyl-CoA] are the concentration values used in the slow-onset inhibition assays:𝐾𝑖=𝐾app𝑖[]1+NADH/𝐾NADH+Crotonoyl-CoA/𝐾Crotonoyl-CoA.(8)

The data from two-step mechanism studies of PfENR inhibition by [FeII(CN)5(INH)]3− were fitted to (9), in which 𝐾𝑖 is the dissociation constant for the initial enzyme-inhibitor binary complex (EI), 𝑘5 is the forward isomerization rate constant for conversion of EI to EI*, 𝑘6 is the reverse isomerization rate constant (conversion of EI* to EI), 𝐼 is the concentration of inhibitor ([FeII(CN)5(INH)]3−), and 𝑘inact is the apparent first-order constant for enzyme inhibition:𝑘inact=𝑘6+𝑘5𝐼𝐾𝑖+𝐼.(9)

3. Results and Discussion

3.1. Initial Velocity Patterns and Kinetic Constants

Initial velocity patterns were determined using either crotonoyl-CoA or NADH as the variable substrate to distinguish between a sequential and a double displacement (“ping-pong”) mechanism, providing information regarding the order of substrate binding. Double-reciprocal plot [40] analysis revealed a parallel line pattern for both substrates (Figure 3), which is consistent with a double displacement mechanism. In this mechanism, a catalytic event can occur when only one of the two substrates is bound to the enzyme. A catalytic event occurs after the binding of substrate A to the enzyme (E), but prior to the binding of substrate B. After the first catalytic event, a portion of substrate A remains in the active site of the enzyme, whereas the other portion denoted as product P dissociates. The first catalytic event produces a different stable enzyme form (F) to which substrate B can bind and proceed with catalysis to generate product Q and regenerate free enzyme (E). The data were tentatively fitted to (1), yielding the following values for the true steady-state kinetic constants: 𝑉max=5.8 (±0.3) U mg−1, 𝑘cat=6.2 (±0.4) s−1, 𝐾NADH=85 (±9) μM, 𝐾Crotonoyl-CoA=63 (±6) μM, 𝑘cat/𝐾NADH=7.3 (±0.9) ×104 M−1 s−1, and 𝑘cat/𝐾Crotonoyl-CoA=9.8 (±1.1) ×104 M−1 s−1.

There are, however, some conditions where lines that seem parallel really are not. For an Ordered Bi Bi system in which the dissociation constant for substrate A (𝐾𝑖𝑎) is very much smaller than the Michaelis-Menten constant for A (𝐾𝑎 or 𝐾MA), the family of reciprocal plots intersect far to the left of the 1/𝑣 axis and far below the 1/[A] axis [36]. In addition, seemingly parallel initial velocity patterns can also occur in rapid equilibrium random systems where the binding of one substrate strongly inhibits the binding of the other.

3.2. Equilibrium Binding of Ligands to PfENR

Equilibrium binary complex formation experiments were employed to confirm, or refute, the proposed enzyme mechanism. There was an enhancement in nucleotide fluorescence upon NADH binding to PfENR, whereas butyryl-CoA binding to PfENR resulted in intrinsic protein fluorescence quench. Titration of PfENR with NADH was hyperbolic (Figure 4(a)), and the data were fitted to (3) yielding a value of 398 (±53) μM for the overall dissociation constant (𝐾𝑑) of NADH. Enzyme titration with butyryl-CoA was sigmoidal (Figure 4(b)), and the data were best fitted to (4) yielding a value for 𝐾 for butyryl-CoA of 1151 (±638) μM. The value for 𝐾 represents a mean dissociation constant for PfENR : butyryl-CoA binary complex formation, which is comprised of interaction factors and the intrinsic dissociation constant [36, 37]. A Hill plot of the data (Figure 4(c)) fitted to (4) yielded a value larger than one for (1.6 ± 0.1). The latter indicates moderate positive homotropic cooperativity on butyryl-CoA binding to tetrameric PfENR, which is further supported by the upward-curved double-reciprocal plot (Figure 4(c)—inset). No intrinsic protein fluorescence change upon addition of crotonoyl-CoA to PfENR could be observed, suggesting that this compound cannot bind to free PfENR enzyme. Protein fluorescence titration with NAD+ could not be performed due to large inner filter effect.

3.3. Pre-Steady-State Kinetics

Steady-state kinetics data suggested that PfENR follows a double displacement (“ping-pong”) mechanism. Fluorescence spectroscopy data showed that NADH binds to free enzyme whereas it appears that crotonoyl-CoA cannot. Although butyryl-CoA binding to free enzyme resulted in protein fluorescence quench, it is not certain whether crotonoyl-CoA cannot bind to free enzyme or its binding results in no change in protein fluorescence. Hence, the Ordered Bi Bi and rapid equilibrium random mechanisms could not be ruled out. To try to address this uncertainty, pre-steady-state kinetics measurements were carried out to show whether or not there is hydride transfer in the absence of crotonoyl-CoA as expected for the double displacement mechanism with NADH binding first. Moreover, determination of a value for the apparent first-order rate constant for the course of NADH→NAD+ conversion could also indicate whether or not this step is rate limiting. It should be pointed out that NADH concentration was in excess as compared to PfENR concentration, so that any burst in product formation could be reliably assessed [41, 42]. Fitting the pre-steady-state trace showing a single exponential decay upon NADH→NAD+ conversion (Figure 5) to (5) yielded a value of 1100 (±29) s−1 for the apparent first-order rate constant for product formation (1 cm pathlength). The NADH→NAD+ conversion step has thus no contribution to rate-limiting step(s) as its value is much larger than the catalytic rate constant value. There was an apparent accumulation of observable product (NAD+) on the enzyme active sites when high concentration of PfENR (3.14 μM) was mixed with NADH (200 μM). Taken together, steady-state kinetics, pre-steady-state kinetics, and equilibrium binding data are consistent with a double displacement mechanism for PfENR (Figure 6). In this mechanism, a catalytic event occurs after the binding of NADH to free enzyme (E), generating and releasing NAD+ into solution prior to the binding of crotonoyl-CoA to a modified form of the enzyme (F), followed by formation of butyryl-CoA and regeneration of free enzyme (E) (Figure 6). The enzyme mechanism has been shown to be ordered sequential with enoyl-ACP substrate binding first for Staphylococcus aureus enoyl reductase [43] and ordered sequential with NADH binding first for fabV-encoded enoyl reductase from Burkholderia mallei [44]. Although uncommon for reductases/dehydrogenases, double displacement mechanisms have been described for Enterobacter cloacae nitroreductase [45], Salmonella typhimurium histidinol dehydrogenase [46], bovine liver mitochondrial dihydroorotate dehydrogenase [47], and mammalian pyruvate dehydrogenase complex [48]. Of course, no mechanism is ever proved; at best one can say that given the current data the likely mechanism for PfENR is double displacement. Further efforts are needed to provide support for either the double displacement mechanism or any other that may give a better description of the experimental results.

Based on substrate binding to immobilized PfENR by Surface Plasmon Resonance (SPR) results, it has been proposed that PfENR enzyme mechanism is random order of substrate addition [49]. Although there have been reports on determination of apparent steady-state kinetics parameters for PfENR [38, 50, 51], to the best of our knowledge, no data have been presented for determination of true steady-state kinetic parameters on which to base a proposal for PfENR enzyme mechanism in solution. To determine the true steady-state kinetic parameters and initial velocity patterns for distinguishing mechanisms, enzyme activity measurements have to be carried out in the presence of varying concentrations of one substrate and several fixed-varied concentrations of the other substrate. Interestingly, the 𝐾𝑑 value of 398 μM for NADH is approximately 8-fold larger than a previously reported value (51.6 μM) by SPR [49]. These results are somewhat puzzling as equilibrium constants determined by SPR usually match well the values obtained by solution techniques (fluorescence spectroscopy, isothermal titration calorimetry). In addition, a 𝐾𝑑 value of 62.5 μM for crotonoyl-CoA has been determined by SPR measurements [49], whereas no change in intrinsic PfENR protein fluorescence could be detected upon addition of crotonoyl-CoA to free enzyme. Incidentally, the true steady-state value for 𝑘cat (6.2 s−1) reported here is approximately 4-fold larger as compared to an apparent value (1.62 s−1) described elsewhere [51].

3.4. Slow-Onset Inhibition Kinetics of PfENR by [FeII(CN)5·(INH)]3−

We have previously shown that [FeII(CN)5·(INH)]3− compound (Figure 2(b)) is a slow-onset inhibitor of wild-type and drug-resistant MtInhAs [30, 31, 52]. We then anticipated that [FeII(CN)5(INH)]3− may also be an effective inhibitor of PfENR enzyme activity. Prior to embarking on a full analysis of the mode of inhibition of PfENR activity by [FeII(CN)5(INH)]3− compound, the enzyme was preincubated with inhibitor, and aliquots were withdrawn at different times to measure the remaining residual activity. PfENR inhibition by [FeII(CN)5(INH)]3− was time dependent (data not shown). Thereby, PfENR slow-binding experiments were carried out in the presence of various concentrations of [FeII(CN)5(INH)]3− to try to elucidate the mechanism of the observed time-dependent enzyme inhibition. The progress curves for enzyme activity indicate that [FeII(CN)5(INH)]3− is a slow-onset inhibitor of PfENR (Figure 7).

Three basic slow-onset inhibition mechanisms have been described by Morrison and Walsh (Figure 8; [32]). In mechanism A (Figure 8), it is assumed that the time-dependent enzyme inactivation observed is resultant of a slow bimolecular enzyme-inhibitor interaction, because the inhibitor concentration is low and/or due to barriers that the inhibitor encounters in its binding at the enzyme’s active site. In mechanism B (Figure 8), it is assumed that there is an initial rapid enzyme-inhibitor interaction to form EI, which then undergoes a slow isomerization (conformational change) to EI* [53]. For inhibitors that conform to this mechanism, the overall dissociation constant (𝐾𝑖) must be lower than the dissociation constant for initially formed complex EI (𝐾𝑖), a condition that can be satisfied only when 𝑘6<𝑘5(𝐾𝑖=𝐾𝑖𝑘6/(𝑘5+𝑘6) and when 𝑘5 and 𝑘6 are slower than all other steps [32]. According to mechanism B, in the presence of a fixed concentration of substrate (𝐴), the value of 𝑘obs increases hyperbolically as a function of inhibitor concentration with limiting values of 𝑘6 and 𝑘5+𝑘6 at, respectively, zero and infinite concentrations of inhibitor. As described in Section 2, 𝑘obs in (6) represents the apparent first-order rate constant for the establishment of the equilibrium between EI and EI* complexes (the observed rate constant for formation of EI* or the apparent isomerization rate constant). A third mechanism (mechanism C) has also been described to account for slow-binding inhibition of enzyme-catalysed reactions (Figure 8). This mechanism involves a slow initial isomerization of two forms of free enzyme in solution (E and E*), followed by a rapid binding of inhibitor to form E*I complex. The 𝑘f and 𝑘r stand, respectively, for the forward and reverse rate constants for the conversion of free enzyme (E) to a form of free enzyme (E*) to which the inhibitor binds. However, fitting the 𝑘obs values as a function of [FeII(CN)5(INH)]3− inhibitor concentration yielded a plot that did not allow to draw solid conclusions on the mechanism. On the other hand, fitting the initial rate values (𝑣0) as a function of [FeII(CN)5(INH)]3− inhibitor concentration to (7) yielded a value of 40 (±11) μM for 𝐾app𝑖. A “true” value of 11 (±3) μM for the apparent dissociation constant 𝐾𝑖 for the initial enzyme-inhibitor binary complex (EI) could thus be calculated by employing (8).

3.5. Two-Step Mechanism of PfENR Inhibition by [FeII(CN)5·(INH)]3−

The two-step mechanism of PfENR inhibition by [FeII(CN)5(INH)]3− was assessed by an experiment in which the enzyme was preincubated in the presence of various concentrations of inhibitor and the decrease of PfENR activity was assessed by withdrawing aliquots of enzyme from preincubation mixtures at appropriate intervals, followed by rapid measurement of initial velocity. These results were analyzed in terms of a two-step inhibition mechanism in which the initial rapid binding of the inhibitor to enzyme is followed by a second slow step that results in the final enzyme-inhibitor complex. Since the value of 𝐾𝑖 is much larger than 𝐾𝑖, the inhibitor concentration in the assay solution is lower than the 𝐾𝑖 value ([I]<𝐾𝑖), and there is a large excess of the competitive substrates NADH and crotonoyl-CoA over [FeII(CN)5(INH)]3−, the steady-state concentration of EI in the assay solution is kinetically insignificant [31]. The apparent first-order rate constants for PfENR inactivation (𝑘inact) could be calculated from a plot of residual enzyme activity against time for each inhibitor concentration. A plot of 𝑘inact values versus inhibitor concentration shows that the 𝑘inact values increase hyperbolically as a function of inhibitor concentration (Figure 9), which is consistent with mechanism B (Figure 8). These data were thus fitted to (9), which describes the slow-binding inhibition for mechanism B, yielding values of 58 (±12) μM for 𝐾𝑖, 53 (±2) s−1 for 𝑘5, and 4 (±2) s−1 for 𝑘6. A value of 4 μM can thus be calculated for the overall inhibition constant (𝐾𝑖=𝐾𝑖𝑘6/(𝑘5+𝑘6). It should be pointed out that an independent estimate for the 𝐾𝑖 value could be obtained by plotting the final steady-state velocity (𝑣𝑠) as a function of [FeII(CN)5(INH)]3− inhibitor concentration and fitting the data to (7) for competitive inhibition. However, the data were not consistent as the experimental data appear to depart somewhat from the values predicted by fitting the data to (6).

The two-step mechanism followed by [FeII(CN)5·(INH)]3− inactivation of PfENR involves an enzyme-inhibitor isomerization step. The enzyme isomerization could involve, for instance, a protein conformational change or alteration of water structure at the active site. For the inhibitor, the change could be hydration/dehydration of a carbonyl, change in ionization state, or a conformational change. These alterations could take place in the enzyme-inhibitor complex or in solution [54].

4. Concluding Remarks

The weakness in the use of classical enzyme inhibitors as drugs for clinical conditions is that inhibition results in the upstream accumulation of the substrate for the enzyme, which may overcome the inhibition. By contrast, the build-up of substrate cannot have any effect on the isomerization of enzyme-inhibitor complex typical of the slow-onset mechanism and hence reversal of the inhibition [32]. Moreover, the low half-time value for the limiting step for inhibitor dissociation from the binary complex is a desirable feature since it may be expected to enhance inhibitor’s effectiveness [54]. It has recently been proposed by scientists at GlaxoSmithKline (USA) that measurement of the dissociative half-life is a crucial metric of compound optimization and could be a key indicator of in vivo duration of efficacy and target selectivity [55]. It has recently been shown that [FeII(CN)5(INH)]3− compound has a promising profile in toxicological assays and that it is amenable to large-scale production [56]. These features indicate that [FeII(CN)5(INH)]3− compound is a viable candidate for further studies (including in vivo efficacy) and that it may represent a prototype for the development of chemotherapeutic agents to treat malaria.

Acknowledgments

This work was supported by Millennium Initiative Program (CNPq), Brazil, to L. H. P. Silva The authors also acknowledge support by National Institute of Science and Technology on Tuberculosis (Decit/SCTIE/MS-MCT-CNPq-FNDCT-CAPES). P. S. M. Medeiros and R. G. Ducati thank fellowships awarded by CNPq.