Malaria parasites require a divergent heme oxygenase for apicoplast gene expression and biogenesis.

Malaria parasites have evolved unusual metabolic adaptations that specialize them for growth within heme-rich human erythrocytes. During blood-stage infection, Plasmodium falciparum parasites internalize and digest abundant host hemoglobin within the digestive vacuole. This massive catabolic process generates copious free heme, most of which is biomineralized into inert hemozoin. Parasites also express a divergent heme oxygenase (HO)-like protein (PfHO) that lacks key active-site residues and has lost canonical HO activity. The cellular role of this unusual protein that underpins its retention by parasites has been unknown. To unravel PfHO function, we first determined a 2.8 Å-resolution X-ray structure that revealed a highly α-helical fold indicative of distant HO homology. Localization studies unveiled PfHO targeting to the apicoplast organelle, where it is imported and undergoes N-terminal processing but retains most of the electropositive transit peptide. We observed that conditional knockdown of PfHO was lethal to parasites, which died from defective apicoplast biogenesis and impaired isoprenoid-precursor synthesis. Complementation and molecular-interaction studies revealed an essential role for the electropositive N-terminus of PfHO, which selectively associates with the apicoplast genome and enzymes involved in nucleic acid metabolism and gene expression. PfHO knockdown resulted in a specific deficiency in levels of apicoplast-encoded RNA but not DNA. These studies reveal an essential function for PfHO in apicoplast maintenance and suggest that Plasmodium repurposed the conserved HO scaffold from its canonical heme-degrading function in the ancestral chloroplast to fulfill a critical adaptive role in organelle gene expression.

Identification of a divalent metal transporter required for cellular iron metabolism in malaria parasites.

Plasmodium falciparum malaria parasites invade and multiply inside red blood cells (RBCs), the most iron-rich compartment in humans. Like all cells, P. falciparum requires nutritional iron to support essential metabolic pathways, but the critical mechanisms of iron acquisition and trafficking during RBC infection have remained obscure. Parasites internalize and liberate massive amounts of heme during large-scale digestion of RBC hemoglobin within an acidic food vacuole (FV) but lack a heme oxygenase to release porphyrin-bound iron. Although most FV heme is sequestered into inert hemozoin crystals, prior studies indicate that trace heme escapes biomineralization and is susceptible to non-enzymatic degradation within the oxidizing FV environment to release labile iron. Parasites retain a homolog of divalent metal transporter 1 (DMT1), a known mammalian iron transporter, but its role in P. falciparum iron acquisition has not been tested. Our phylogenetic studies indicate that P. falciparum DMT1 (PfDMT1) retains conserved molecular features critical for metal transport. We localized this protein to the FV membrane and defined its orientation in an export-competent topology. Conditional knockdown of PfDMT1 expression is lethal to parasites, which display broad cellular defects in iron-dependent functions, including impaired apicoplast biogenesis and mitochondrial polarization. Parasites are selectively rescued from partial PfDMT1 knockdown by supplementation with exogenous iron, but not other metals. These results support a cellular paradigm whereby PfDMT1 is the molecular gatekeeper to essential iron acquisition by blood-stage malaria parasites and suggest that therapeutic targeting of PfDMT1 may be a potent antimalarial strategy.

Biogenesis of cytochromes c and c 1 in the electron transport chain of malaria parasites.

Plasmodium malaria parasites retain an essential mitochondrional electron transport chain (ETC) that is critical for growth within humans and mosquitoes and a key antimalarial drug target. ETC function requires cytochromes c and c 1 that are unusual among heme proteins due to their covalent binding to heme via conserved CXXCH sequence motifs. Heme attachment to these proteins in most eukaryotes requires the mitochondrial enzyme holocytochrome c synthase (HCCS) that binds heme and the apo cytochrome to facilitate biogenesis of the mature cytochrome c or c 1. Although humans encode a single bifunctional HCCS that attaches heme to both proteins, Plasmodium parasites are like yeast and encode two separate HCCS homologs thought to be specific for heme attachment to cyt c (HCCS) or cyt c 1 (HCC1S). To test the function and specificity of P. falciparum HCCS and HCC1S, we used CRISPR/Cas9 to tag both genes for conditional expression. HCC1S knockdown selectively impaired cyt c 1 biogenesis and caused lethal ETC dysfunction that was not reversed by over-expression of HCCS. Knockdown of HCCS caused a more modest growth defect but strongly sensitized parasites to mitochondrial depolarization by proguanil, revealing key defects in ETC function. These results and prior heterologous studies in E. coli of cyt c hemylation by P. falciparum HCCS and HCC1S strongly suggest that both homologs are essential for mitochondrial ETC function and have distinct specificities for biogenesis of cyt c and c 1, respectively, in parasites. This study lays a foundation to develop novel strategies to selectively block ETC function in malaria parasites.

The interdependence of isoprenoid synthesis and apicoplast biogenesis in malaria parasites.

Isoprenoid precursor synthesis is an ancient and fundamental function of plastid organelles and a critical metabolic activity of the apicoplast in Plasmodium malaria parasites [1-3]. Over the past decade, our understanding of apicoplast properties and functions has increased enormously [4], due in large part to our ability to rescue blood-stage parasites from apicoplast-specific dysfunctions by supplementing cultures with isopentenyl pyrophosphate (IPP), a key output of this organelle [5,6]. In this Pearl, we explore the interdependence between isoprenoid metabolism and apicoplast biogenesis in P. falciparum and highlight critical future questions to answer.

Direct tests of cytochrome c and c1 functions in the electron transport chain of malaria parasites.

The mitochondrial electron transport chain (ETC) of Plasmodium malaria parasites is a major antimalarial drug target, but critical cytochrome (cyt) functions remain unstudied and enigmatic. Parasites express two distinct cyt c homologs (c and c-2) with unusually sparse sequence identity and uncertain fitness contributions. P. falciparum cyt c-2 is the most divergent eukaryotic cyt c homolog currently known and has sequence features predicted to be incompatible with canonical ETC function. We tagged both cyt c homologs and the related cyt c1 for inducible knockdown. Translational repression of cyt c and cyt c1 was lethal to parasites, which died from ETC dysfunction and impaired ubiquinone recycling. In contrast, cyt c-2 knockdown or knockout had little impact on blood-stage growth, indicating that parasites rely fully on the more conserved cyt c for ETC function. Biochemical and structural studies revealed that both cyt c and c-2 are hemylated by holocytochrome c synthase, but UV-vis absorbance and EPR spectra strongly suggest that cyt c-2 has an unusually open active site in which heme is stably coordinated by only a single axial amino acid ligand and can bind exogenous small molecules. These studies provide a direct dissection of cytochrome functions in the ETC of malaria parasites and identify a highly divergent Plasmodium cytochrome c with molecular adaptations that defy a conserved role in eukaryotic evolution.

Direct Tests of Cytochrome Function in the Electron Transport Chain of Malaria Parasites.

The mitochondrial electron transport chain (ETC) of Plasmodium malaria parasites is a major antimalarial drug target, but critical cytochrome functions remain unstudied and enigmatic. Parasites express two distinct cyt c homologs ( c and c -2) with unusually sparse sequence identity and uncertain fitness contributions. P. falciparum cyt c -2 is the most divergent eukaryotic cyt c homolog currently known and has sequence features predicted to be incompatible with canonical ETC function. We tagged both cyt c homologs and the related cyt c 1 for inducible knockdown. Translational repression of cyt c and cyt c 1 was lethal to parasites, which died from ETC dysfunction and impaired ubiquinone recycling. In contrast, cyt c -2 knockdown or knock-out had little impact on blood-stage growth, indicating that parasites rely fully on the more conserved cyt c for ETC function. Biochemical and structural studies revealed that both cyt c and c -2 are hemylated by holocytochrome c synthase, but UV-vis absorbance and EPR spectra strongly suggest that cyt c -2 has an unusually open active site in which heme is stably coordinated by only a single axial amino-acid ligand and can bind exogenous small molecules. These studies provide a direct dissection of cytochrome functions in the ETC of malaria parasites and identify a highly divergent Plasmodium cytochrome c with molecular adaptations that defy a conserved role in eukaryotic evolution. SIGNIFICANCE STATEMENT: Mitochondria are critical organelles in eukaryotic cells that drive oxidative metabolism. The mitochondrion of Plasmodium malaria parasites is a major drug target that has many differences from human cells and remains poorly studied. One key difference from humans is that malaria parasites express two cytochrome c proteins that differ significantly from each other and play untested and uncertain roles in the mitochondrial electron transport chain (ETC). Our study revealed that one cyt c is essential for ETC function and parasite viability while the second, more divergent protein has unusual structural and biochemical properties and is not required for growth of blood-stage parasites. This work elucidates key biochemical properties and evolutionary differences in the mitochondrial ETC of malaria parasites.

Multiple-Streams Focusing-Based Cell Separation in High Viscoelasticity Flow.

Viscoelastic flow has been widely used in microfluidic particle separation processes, in which particles get focused on the channel center in diluted viscoelastic flow. In this paper, the transition from single-stream focusing to multiple-streams focusing (MSF) in high viscoelastic flow is observed, which is applied for cell separation processes. Particle focusing stream bifurcation is caused by the balance between elastic force and viscoelastic secondary flow drag force. The influence of cell physical properties, such as cell dimension, shape, and deformability, on the formation of multiple-streams focusing is studied in detail. Particle separation is realized utilizing different separation criteria. The size-based separation of red (RBC) and white (WBC) blood cells is demonstrated in which cells get focused in different streams based on their dimension difference. Cells with different deformabilities get stretched in the viscoelastic flow, leading to the change of focusing streams, and this property is harnessed to separate red blood cells infected with the malaria parasite, Plasmodium falciparum. The achieved results promote our understanding of particle movement in the high viscoelastic flow and enable new particle manipulation and separation processes for sample treatment in biofluids.

Critical role for isoprenoids in apicoplast biogenesis by malaria parasites.

Isopentenyl pyrophosphate (IPP) is an essential metabolic output of the apicoplast organelle in Plasmodium falciparum malaria parasites and is required for prenylation-dependent vesicular trafficking and other cellular processes. We have elucidated a critical and previously uncharacterized role for IPP in apicoplast biogenesis. Inhibiting IPP synthesis blocks apicoplast elongation and inheritance by daughter merozoites, and apicoplast biogenesis is rescued by exogenous IPP and polyprenols. Knockout of the only known isoprenoid-dependent apicoplast pathway, tRNA prenylation by MiaA, has no effect on blood-stage parasites and thus cannot explain apicoplast reliance on IPP. However, we have localized an annotated polyprenyl synthase (PPS) to the apicoplast. PPS knockdown is lethal to parasites, rescued by IPP and long- (C50) but not short-chain (≤C20) prenyl alcohols, and blocks apicoplast biogenesis, thus explaining apicoplast dependence on isoprenoid synthesis. We hypothesize that PPS synthesizes long-chain polyprenols critical for apicoplast membrane fluidity and biogenesis. This work critically expands the paradigm for isoprenoid utilization in malaria parasites and identifies a novel essential branch of apicoplast metabolism suitable for therapeutic targeting.

Divergent acyl carrier protein decouples mitochondrial Fe-S cluster biogenesis from fatty acid synthesis in malaria parasites.

Most eukaryotic cells retain a mitochondrial fatty acid synthesis (FASII) pathway whose acyl carrier protein (mACP) and 4-phosphopantetheine (Ppant) prosthetic group provide a soluble scaffold for acyl chain synthesis and biochemically couple FASII activity to mitochondrial electron transport chain (ETC) assembly and Fe-S cluster biogenesis. In contrast, the mitochondrion of Plasmodium falciparum malaria parasites lacks FASII enzymes yet curiously retains a divergent mACP lacking a Ppant group. We report that ligand-dependent knockdown of mACP is lethal to parasites, indicating an essential FASII-independent function. Decyl-ubiquinone rescues parasites temporarily from death, suggesting a dominant dysfunction of the mitochondrial ETC. Biochemical studies reveal that Plasmodium mACP binds and stabilizes the Isd11-Nfs1 complex required for Fe-S cluster biosynthesis, despite lacking the Ppant group required for this association in other eukaryotes, and knockdown of parasite mACP causes loss of Nfs1 and the Rieske Fe-S protein in ETC complex III. This work reveals that Plasmodium parasites have evolved to decouple mitochondrial Fe-S cluster biogenesis from FASII activity, and this adaptation is a shared metabolic feature of other apicomplexan pathogens, including Toxoplasma and Babesia. This discovery unveils an evolutionary driving force to retain interaction of mitochondrial Fe-S cluster biogenesis with ACP independent of its eponymous function in FASII.

Doxycycline has distinct apicoplast-specific mechanisms of antimalarial activity.

Doxycycline (DOX) is a key antimalarial drug thought to kill Plasmodium parasites by blocking protein translation in the essential apicoplast organelle. Clinical use is primarily limited to prophylaxis due to delayed second-cycle parasite death at 1-3 µM serum concentrations. DOX concentrations > 5 µM kill parasites with first-cycle activity but are thought to involve off-target mechanisms outside the apicoplast. We report that 10 µM DOX blocks apicoplast biogenesis in the first cycle and is rescued by isopentenyl pyrophosphate, an essential apicoplast product, confirming an apicoplast-specific mechanism. Exogenous iron rescues parasites and apicoplast biogenesis from first- but not second-cycle effects of 10 µM DOX, revealing that first-cycle activity involves a metal-dependent mechanism distinct from the delayed-death mechanism. These results critically expand the paradigm for understanding the fundamental antiparasitic mechanisms of DOX and suggest repurposing DOX as a faster acting antimalarial at higher dosing whose multiple mechanisms would be expected to limit parasite resistance.

In-Cell Enzymology To Probe His-Heme Ligation in Heme Oxygenase Catalysis.

Heme oxygenase (HO) is a ubiquitous enzyme with key roles in inflammation, cell signaling, heme disposal, and iron acquisition. HO catalyzes the oxidative conversion of heme to biliverdin (BV) using a conserved histidine to coordinate the iron atom of bound heme. This His-heme interaction has been regarded as being essential for enzyme activity, because His-to-Ala mutants fail to convert heme to biliverdin in vitro. We probed a panel of proximal His mutants of cyanobacterial, human, and plant HO enzymes using a live-cell activity assay based on heterologous co-expression in Escherichia coli of each HO mutant and a fluorescent biliverdin biosensor. In contrast to in vitro studies with purified proteins, we observed that multiple HO mutants retained significant activity within the intracellular environment of bacteria. X-ray crystallographic structures of human HO1 H25R with bound heme and additional functional studies suggest that HO mutant activity inside these cells does not involve heme ligation by a proximal amino acid. Our study reveals unexpected plasticity in the active site binding interactions with heme that can support HO activity within cells, suggests important contributions by the surrounding active site environment to HO catalysis, and can guide efforts to understand the evolution and divergence of HO function.

Deconvoluting heme biosynthesis to target blood-stage malaria parasites.

Heme metabolism is central to blood-stage infection by the malaria parasite Plasmodium falciparum. Parasites retain a heme biosynthesis pathway but do not require its activity during infection of heme-rich erythrocytes, where they can scavenge host heme to meet metabolic needs. Nevertheless, heme biosynthesis in parasite-infected erythrocytes can be potently stimulated by exogenous 5-aminolevulinic acid (ALA), resulting in accumulation of the phototoxic intermediate protoporphyrin IX (PPIX). Here we use photodynamic imaging, mass spectrometry, parasite gene disruption, and chemical probes to reveal that vestigial host enzymes in the cytoplasm of Plasmodium-infected erythrocytes contribute to ALA-stimulated heme biosynthesis and that ALA uptake depends on parasite-established permeability pathways. We show that PPIX accumulation in infected erythrocytes can be harnessed for antimalarial chemotherapy using luminol-based chemiluminescence and combinatorial stimulation by low-dose artemisinin to photoactivate PPIX to produce cytotoxic reactive oxygen. This photodynamic strategy has the advantage of exploiting host enzymes refractory to resistance-conferring mutations.

Determination of Hydrogen Bond Structure in Water versus Aprotic Environments To Test the Relationship Between Length and Stability.

Hydrogen bonds profoundly influence the architecture and activity of biological macromolecules. Deep appreciation of hydrogen bond contributions to biomolecular function thus requires a detailed understanding of hydrogen bond structure and energetics and the relationship between these properties. Hydrogen bond formation energies (ΔGf) are enormously more favorable in aprotic solvents than in water, and two classes of contributing factors have been proposed to explain this energetic difference, focusing respectively on the isolated and hydrogen-bonded species: (I) water stabilizes the dissociated donor and acceptor groups much better than aprotic solvents, thereby reducing the driving force for hydrogen bond formation; and (II) water lengthens hydrogen bonds compared to aprotic environments, thereby decreasing the potential energy within the hydrogen bond. Each model has been proposed to provide a dominant contribution to ΔGf, but incisive tests that distinguish the importance of these contributions are lacking. Here we directly test the structural basis of model II. Neutron crystallography, NMR spectroscopy, and quantum mechanical calculations demonstrate that O-H···O hydrogen bonds in crystals, chloroform, acetone, and water have nearly identical lengths and very similar potential energy surfaces despite ΔGf differences >8 kcal/mol across these solvents. These results rule out a substantial contribution from solvent-dependent differences in hydrogen bond structure and potential energy after association (model II) and thus support the conclusion that differences in hydrogen bond ΔGf are predominantly determined by solvent interactions with the dissociated groups (model I). These findings advance our understanding of universal hydrogen-bonding interactions and have important implications for biology and engineering.

The heme biosynthesis pathway is essential for Plasmodium falciparum development in mosquito stage but not in blood stages.

Heme is an essential cofactor for aerobic organisms. Its redox chemistry is central to a variety of biological functions mediated by hemoproteins. In blood stages, malaria parasites consume most of the hemoglobin inside the infected erythrocytes, forming nontoxic hemozoin crystals from large quantities of heme released during digestion. At the same time, the parasites possess a heme de novo biosynthetic pathway. This pathway in the human malaria parasite Plasmodium falciparum has been considered essential and is proposed as a potential drug target. However, we successfully disrupted the first and last genes of the pathway, individually and in combination. These knock-out parasite lines, lacking 5-aminolevulinic acid synthase and/or ferrochelatase (FC), grew normally in blood-stage culture and exhibited no changes in sensitivity to heme-related antimalarial drugs. We developed a sensitive LC-MS/MS assay to monitor stable isotope incorporation into heme from its precursor 5-[(13)C4]aminolevulinic acid, and this assay confirmed that de novo heme synthesis was ablated in FC knock-out parasites. Disrupting the FC gene also caused no defects in gametocyte generation or maturation but resulted in a greater than 70% reduction in male gamete formation and completely prevented oocyst formation in female Anopheles stephensi mosquitoes. Our data demonstrate that the heme biosynthesis pathway is not essential for asexual blood-stage growth of P. falciparum parasites but is required for mosquito transmission. Drug inhibition of pathway activity is therefore unlikely to provide successful antimalarial therapy. These data also suggest the existence of a parasite mechanism for scavenging host heme to meet metabolic needs.

The peculiarities and paradoxes of Plasmodium heme metabolism.

For over a century, heme metabolism has been recognized to play a central role during intraerythrocytic infection by Plasmodium parasites, the causative agent of malaria. Parasites liberate vast quantities of potentially cytotoxic heme as a by-product of hemoglobin catabolism within the digestive vacuole, where heme is predominantly sequestered as inert crystalline hemozoin. Plasmodium spp. also utilize heme as a metabolic cofactor. Despite access to abundant host-derived heme, parasites paradoxically maintain a biosynthetic pathway. This pathway has been assumed to produce the heme incorporated into mitochondrial cytochromes that support electron transport. In this review, we assess our current understanding of the love-hate relationship between Plasmodium parasites and heme, we discuss recent studies that clarify several long-standing riddles about heme production and utilization by parasites, and we consider remaining challenges and opportunities for understanding and targeting heme metabolism within parasites.

Quantitative dissection of hydrogen bond-mediated proton transfer in the ketosteroid isomerase active site.

Hydrogen bond networks are key elements of protein structure and function but have been challenging to study within the complex protein environment. We have carried out in-depth interrogations of the proton transfer equilibrium within a hydrogen bond network formed to bound phenols in the active site of ketosteroid isomerase. We systematically varied the proton affinity of the phenol using differing electron-withdrawing substituents and incorporated site-specific NMR and IR probes to quantitatively map the proton and charge rearrangements within the network that accompany incremental increases in phenol proton affinity. The observed ionization changes were accurately described by a simple equilibrium proton transfer model that strongly suggests the intrinsic proton affinity of one of the Tyr residues in the network, Tyr16, does not remain constant but rather systematically increases due to weakening of the phenol-Tyr16 anion hydrogen bond with increasing phenol proton affinity. Using vibrational Stark spectroscopy, we quantified the electrostatic field changes within the surrounding active site that accompany these rearrangements within the network. We were able to model these changes accurately using continuum electrostatic calculations, suggesting a high degree of conformational restriction within the protein matrix. Our study affords direct insight into the physical and energetic properties of a hydrogen bond network within a protein interior and provides an example of a highly controlled system with minimal conformational rearrangements in which the observed physical changes can be accurately modeled by theoretical calculations.

Direct tests of enzymatic heme degradation by the malaria parasite Plasmodium falciparum.

Malaria parasites generate vast quantities of heme during blood stage infection via hemoglobin digestion and limited de novo biosynthesis, but it remains unclear if parasites metabolize heme for utilization or disposal. Recent in vitro experiments with a heme oxygenase (HO)-like protein from Plasmodium falciparum suggested that parasites may enzymatically degrade some heme to the canonical HO product, biliverdin (BV), or its downstream metabolite, bilirubin (BR). To directly test for BV and BR production by P. falciparum parasites, we DMSO-extracted equal numbers of infected and uninfected erythrocytes and developed a sensitive LC-MS/MS assay to quantify these tetrapyrroles. We found comparable low levels of BV and BR in both samples, suggesting the absence of HO activity in parasites. We further tested live parasites by targeted expression of a fluorescent BV-binding protein within the parasite cytosol, mitochondrion, and plant-like plastid. This probe could detect exogenously added BV but gave no signal indicative of endogenous BV production within parasites. Finally, we recombinantly expressed and tested the proposed heme degrading activity of the HO-like protein, PfHO. Although PfHO bound heme and protoporphyrin IX with modest affinity, it did not catalyze heme degradation in vivo within bacteria or in vitro in UV absorbance and HPLC assays. These observations are consistent with PfHO's lack of a heme-coordinating His residue and suggest an alternative function within parasites. We conclude that P. falciparum parasites lack a canonical HO pathway for heme degradation and thus rely fully on alternative mechanisms for heme detoxification and iron acquisition during blood stage infection.

Quantitative, directional measurement of electric field heterogeneity in the active site of ketosteroid isomerase.

Understanding the electrostatic forces and features within highly heterogeneous, anisotropic, and chemically complex enzyme active sites and their connection to biological catalysis remains a longstanding challenge, in part due to the paucity of incisive experimental probes of electrostatic properties within proteins. To quantitatively assess the landscape of electrostatic fields at discrete locations and orientations within an enzyme active site, we have incorporated site-specific thiocyanate vibrational probes into multiple positions within bacterial ketosteroid isomerase. A battery of X-ray crystallographic, vibrational Stark spectroscopy, and NMR studies revealed electrostatic field heterogeneity of 8 MV/cm between active site probe locations and widely differing sensitivities of discrete probes to common electrostatic perturbations from mutation, ligand binding, and pH changes. Electrostatic calculations based on active site ionization states assigned by literature precedent and computational pK(a) prediction were unable to quantitatively account for the observed vibrational band shifts. However, electrostatic models of the D40N mutant gave qualitative agreement with the observed vibrational effects when an unusual ionization of an active site tyrosine with a pK(a) near 7 was included. UV-absorbance and (13)C NMR experiments confirmed the presence of a tyrosinate in the active site, in agreement with electrostatic models. This work provides the most direct measure of the heterogeneous and anisotropic nature of the electrostatic environment within an enzyme active site, and these measurements provide incisive benchmarks for further developing accurate computational models and a foundation for future tests of electrostatics in enzymatic catalysis.

Hydrogen bonding in the active site of ketosteroid isomerase: electronic inductive effects and hydrogen bond coupling.

Computational studies are performed to analyze the physical properties of hydrogen bonds donated by Tyr16 and Asp103 to a series of substituted phenolate inhibitors bound in the active site of ketosteroid isomerase (KSI). As the solution pK(a) of the phenolate increases, these hydrogen bond distances decrease, the associated nuclear magnetic resonance (NMR) chemical shifts increase, and the fraction of protonated inhibitor increases, in agreement with prior experiments. The quantum mechanical/molecular mechanical calculations provide insight into the electronic inductive effects along the hydrogen bonding network that includes Tyr16, Tyr57, and Tyr32, as well as insight into hydrogen bond coupling in the active site. The calculations predict that the most-downfield NMR chemical shift observed experimentally corresponds to the Tyr16-phenolate hydrogen bond and that Tyr16 is the proton donor when a bound naphtholate inhibitor is observed to be protonated in electronic absorption experiments. According to these calculations, the electronic inductive effects along the hydrogen bonding network of tyrosines cause the Tyr16 hydroxyl to be more acidic than the Asp103 carboxylic acid moiety, which is immersed in a relatively nonpolar environment. When one of the distal tyrosine residues in the network is mutated to phenylalanine, thereby diminishing this inductive effect, the Tyr16-phenolate hydrogen bond becomes longer and the Asp103-phenolate hydrogen bond shorter, as observed in NMR experiments. Furthermore, the calculations suggest that the differences in the experimental NMR data and electronic absorption spectra for pKSI and tKSI, two homologous bacterial forms of the enzyme, are due predominantly to the third tyrosine that is present in the hydrogen bonding network of pKSI but not tKSI. These studies also provide experimentally testable predictions about the impact of mutating the distal tyrosine residues in this hydrogen bonding network on the NMR chemical shifts and electronic absorption spectra.