Regulation-Structured Dynamic Metabolic Model Provides a Potential Mechanism for Delayed Enzyme Response in Denitrification Process.

In a recent study of denitrification dynamics in hyporheic zone sediments, we observed a significant time lag (up to several days) in enzymatic response to the changes in substrate concentration. To explore an underlying mechanism and understand the interactive dynamics between enzymes and nutrients, we developed a trait-based model that associates a community's traits with functional enzymes, instead of typically used species guilds (or functional guilds). This enzyme-based formulation allows to collectively describe biogeochemical functions of microbial communities without directly parameterizing the dynamics of species guilds, therefore being scalable to complex communities. As a key component of modeling, we accounted for microbial regulation occurring through transcriptional and translational processes, the dynamics of which was parameterized based on the temporal profiles of enzyme concentrations measured using a new signature peptide-based method. The simulation results using the resulting model showed several days of a time lag in enzymatic responses as observed in experiments. Further, the model showed that the delayed enzymatic reactions could be primarily controlled by transcriptional responses and that the dynamics of transcripts and enzymes are closely correlated. The developed model can serve as a useful tool for predicting biogeochemical processes in natural environments, either independently or through integration with hydrologic flow simulators.

Targeted quantification of functional enzyme dynamics in environmental samples for microbially mediated biogeochemical processes.

Microbial enzymes catalytically drive biogeochemical processes in environments. The dynamic linkage between functional enzymes and biogeochemical species transformation has, however, rarely been investigated for decades because of the challenges to directly quantify enzymes in environmental samples. The diversity of microorganisms, the low amount of available biomass and the complexity of chemical composition in environmental samples represent the main challenges. To address the diversity challenge, we first identify several signature peptides that are conserved in the targeted enzymes with the same functionality across many phylogenetically diverse microorganisms using metagenome-based protein sequence data. Quantification of the signature peptides then allows estimation of the targeted enzyme abundance. To achieve analyses of the requisite sensitivity for complex environmental samples with low available biomass, we adapted a recently developed ultrasensitive targeted quantification technology, termed high-pressure high-resolution separations with intelligent selection and multiplexing (PRISM) by improving peptide separation efficiency and method detection sensitivity. Nitrate reduction dynamics catalyzed by dissimilatory and assimilatory enzymes in a hyporheic zone sediment was used as an example to demonstrate the application of the enzyme quantification approach. Together with the measurements of biogeochemical species, the approach enables investigating the dynamic linkage between functional enzymes and biogeochemical processes.

Effect of Water Chemistry and Hydrodynamics on Nitrogen Transformation Activity and Microbial Community Functional Potential in Hyporheic Zone Sediment Columns.

Hyporheic zones (HZ) are active biogeochemical regions where groundwater and surface water mix. N transformations in HZ sediments were investigated in columns with a focus on understanding how the dynamic changes in groundwater and surface water mixing affect microbial community and its biogeochemical function with respect to N transformations. The results indicated that denitrification, DNRA, and nitrification rates and products changed quickly in response to changes in water and sediment chemistry, fluid residence time, and groundwater-surface water exchange. These changes were accompanied by the zonation of denitrification functional genes along a 30 cm advective flow path after a total of 6 days' elution of synthetic groundwater with fluid residence time >9.8 h. The shift of microbial functional potential toward denitrification was correlated with rapid NO3- reduction collectively affected by NO3- concentration and fluid residence time, and was resistant to short-term groundwater-surface water exchange on a daily basis. The results implied that variations in microbial functional potential and associated biogeochemical reactions in the HZ may occur at space scales where steep concentration gradients present along the flow path and the variations would respond to dynamic HZ water exchange over different time periods common to natural and managed riverine systems.

Direct involvement of ombB, omaB, and omcB genes in extracellular reduction of Fe(III) by Geobacter sulfurreducens PCA.

The tandem gene clusters orfR-ombB-omaB-omcB and orfS-ombC-omaC-omcC of the metal-reducing bacterium Geobacter sulfurreducens PCA are responsible for trans-outer membrane electron transfer during extracellular reduction of Fe(III)-citrate and ferrihydrite [a poorly crystalline Fe(III) oxide]. Each gene cluster encodes a putative transcriptional factor (OrfR/OrfS), a porin-like outer-membrane protein (OmbB/OmbC), a periplasmic c-type cytochrome (c-Cyt, OmaB/OmaC) and an outer-membrane c-Cyt (OmcB/OmcC). The individual roles of OmbB, OmaB and OmcB in extracellular reduction of Fe(III), however, have remained either uninvestigated or controversial. Here, we showed that replacements of ombB, omaB, omcB, and ombB-omaB with an antibiotic gene in the presence of ombC-omaC-omcC had no impact on reduction of Fe(III)-citrate by G. sulfurreducens PCA. Disruption of ombB, omaB, omcB, and ombB-omaB in the absence of ombC-omaC-omcC, however, severely impaired the bacterial ability to reduce Fe(III)-citrate as well as ferrihydrite. These results unequivocally demonstrate an overlapping role of ombB-omaB-omcB and ombC-omaC-omcC in extracellular Fe(III) reduction by G. sulfurreducens PCA. Involvement of both ombB-omaB-omcB and ombC-omaC-omcC in extracellular Fe(III) reduction reflects the importance of these trans-outer membrane protein complexes in the physiology of this bacterium. Moreover, the kinetics of Fe(III)-citrate and ferrihydrite reduction by these mutants in the absence of ombC-omaC-omcC were nearly identical, which suggests that absence of any protein subunit eliminates function of OmaB/OmbB/OmcB protein complex. Finally, orfS was found to have a negative impact on the extracellular reduction of Fe(III)-citrate and ferrihydrite in G. sulfurreducens PCA probably by serving as a transcriptional repressor.

Redox Linked Flavin Sites in Extracellular Decaheme Proteins Involved in Microbe-Mineral Electron Transfer.

Extracellular microbe-mineral electron transfer is a major driving force for the oxidation of organic carbon in many subsurface environments. Extracellular multi-heme cytochromes of the Shewenella genus play a major role in this process but the mechanism of electron exchange at the interface between cytochrome and acceptor is widely debated. The 1.8 Å x-ray crystal structure of the decaheme MtrC revealed a highly conserved CX8C disulfide that, when substituted for AX8A, severely compromised the ability of S. oneidensis to grow under aerobic conditions. Reductive cleavage of the disulfide in the presence of flavin mononucleotide (FMN) resulted in the reversible formation of a stable flavocytochrome. Similar results were also observed with other decaheme cytochromes, OmcA, MtrF and UndA. The data suggest that these decaheme cytochromes can transition between highly reactive flavocytochromes or less reactive cytochromes, and that this transition is controlled by a redox active disulfide that responds to the presence of oxygen.

Genomic analyses of bacterial porin-cytochrome gene clusters.

The porin-cytochrome (Pcc) protein complex is responsible for trans-outer membrane electron transfer during extracellular reduction of Fe(III) by the dissimilatory metal-reducing bacterium Geobacter sulfurreducens PCA. The identified and characterized Pcc complex of G. sulfurreducens PCA consists of a porin-like outer-membrane protein, a periplasmic 8-heme c-type cytochrome (c-Cyt) and an outer-membrane 12-heme c-Cyt, and the genes encoding the Pcc proteins are clustered in the same regions of genome (i.e., the pcc gene clusters) of G. sulfurreducens PCA. A survey of additionally microbial genomes has identified the pcc gene clusters in all sequenced Geobacter spp. and other bacteria from six different phyla, including Anaeromyxobacter dehalogenans 2CP-1, A. dehalogenans 2CP-C, Anaeromyxobacter sp. K, Candidatus Kuenenia stuttgartiensis, Denitrovibrio acetiphilus DSM 12809, Desulfurispirillum indicum S5, Desulfurivibrio alkaliphilus AHT2, Desulfurobacterium thermolithotrophum DSM 11699, Desulfuromonas acetoxidans DSM 684, Ignavibacterium album JCM 16511, and Thermovibrio ammonificans HB-1. The numbers of genes in the pcc gene clusters vary, ranging from two to nine. Similar to the metal-reducing (Mtr) gene clusters of other Fe(III)-reducing bacteria, such as Shewanella spp., additional genes that encode putative c-Cyts with predicted cellular localizations at the cytoplasmic membrane, periplasm and outer membrane often associate with the pcc gene clusters. This suggests that the Pcc-associated c-Cyts may be part of the pathways for extracellular electron transfer reactions. The presence of pcc gene clusters in the microorganisms that do not reduce solid-phase Fe(III) and Mn(IV) oxides, such as D. alkaliphilus AHT2 and I. album JCM 16511, also suggests that some of the pcc gene clusters may be involved in extracellular electron transfer reactions with the substrates other than Fe(III) and Mn(IV) oxides.

A trans-outer membrane porin-cytochrome protein complex for extracellular electron transfer by Geobacter sulfurreducens†PCA.

The multi-heme, outer membrane c-type cytochrome (c-Cyt) OmcB of Geobacter sulfurreducens was previously proposed to mediate electron transfer across the outer membrane. However, the underlying mechanism has remained uncharacterized. In G. sulfurreducens, the omcB gene is part of two tandem four-gene clusters, each is predicted to encode a transcriptional factor (OrfR/OrfS), a porin-like outer membrane protein (OmbB/OmbC), a periplasmic c-type cytochrome (OmaB/OmaC) and an outer membrane c-Cyt (OmcB/OmcC) respectively. Here, we showed that OmbB/OmbC, OmaB/OmaC and OmcB/OmcC of G. sulfurreducens PCA formed the porin-cytochrome (Pcc) protein complexes, which were involved in transferring electrons across the outer membrane. The isolated Pcc protein complexes reconstituted in proteoliposomes transferred electrons from reduced methyl viologen across the lipid bilayer of liposomes to Fe(III)-citrate and ferrihydrite. The pcc clusters were found in all eight sequenced Geobacter and 11 other bacterial genomes from six different phyla, demonstrating a widespread distribution of Pcc protein complexes in phylogenetically diverse bacteria. Deletion of ombB-omaB-omcB-orfS-ombC-omaC-omcC gene clusters had no impact on the growth of G. sulfurreducens PCA with fumarate but diminished the ability of G. sulfurreducens PCA to reduce Fe(III)-citrate and ferrihydrite. Complementation with the ombB-omaB-omcB gene cluster restored the ability of G. sulfurreducens PCA to reduce Fe(III)-citrate and ferrihydrite.

Investigation of U(VI) adsorption in quartz-chlorite mineral mixtures.

A batch and cryogenic laser-induced time-resolved luminescence spectroscopy investigation of U(VI) adsorbed on quartz-chlorite mixtures with variable mass ratios have been performed under field-relevant uranium concentrations (5×10(-7) M and 5×10(-6) M) in pH 8.1 synthetic groundwater. The U(VI) adsorption Kd values steadily increased as the mass fraction of chlorite increased, indicating preferential sorption to chlorite. For all mineral mixtures, U(VI) adsorption Kd values were lower than that calculated from the assumption of component additivity possibly caused by surface modifications stemming from chlorite dissolution; The largest deviation occurred when the mass fractions of the two minerals were equal. U(VI) adsorbed on quartz and chlorite displayed characteristic individual luminescence spectra that were not affected by mineral mixing. The spectra of U(VI) adsorbed within the mixtures could be simulated by one surface U(VI) species on quartz and two on chlorite. The luminescence intensity decreased in a nonlinear manner as the adsorbed U(VI) concentration increased with increasing chlorite mass fraction-likely due to ill-defined luminescence quenching by both structural Fe/Cr in chlorite, and trace amounts of solubilized and reprecipitated Fe/Cr in the aqueous phase. However, the fractional spectral intensities of U(VI) adsorbed on quartz and chlorite followed the same trend of fractional adsorbed U(VI) concentration in each mineral phase with approximate linear correlations, offering a method to estimate of U(VI) concentration distribution between the mineral components with luminescence spectroscopy.

The X-ray crystal structure of Shewanella oneidensis OmcA reveals new insight at the microbe-mineral interface.

The X-ray crystal structure of Shewanella oneidensis OmcA, an extracellular decaheme cytochrome involved in mineral reduction, was solved to a resolution of 2.7 Å. The four OmcA molecules in the asymmetric unit are arranged so the minimum distance between heme 5 on adjacent OmcA monomers is 9 Å, indicative of a transient OmcA dimer capable of intermolecular electron transfer. A previously identified hematite binding motif was identified near heme 10, forming a hydroxylated surface that would bring a heme 10 electron egress site to ∼10 Å of a mineral surface.

Effect of subgrid heterogeneity on scaling geochemical and biogeochemical reactions: a case of U(VI) desorption.

The effect of subgrid heterogeneity in sediment properties on the rate of uranyl[U(VI)] desorption was investigated using a sediment collected from the U.S. Department of Energy Hanford site. The sediment was sieved into 7 grain size fractions that each exhibited different U(VI) desorption properties. Six columns were assembled using the sediment with its grain size fractions arranged in different spatial configurations to mimic subgrid heterogeneity in reactive transport properties. The apparent rate of U(VI) desorption varied significantly in the columns. Those columns with sediment structures leading to preferential transport had much lower rates of U(VI) desorption than those with relatively homogeneous transport. Modeling analysis indicated that the U(VI) desorption model and parameters characterized from well-mixed reactors significantly overpredicted the measured U(VI) desorption in the columns with preferential transport. A dual domain model, which operationally separates reactive transport properties into two subgrid domains, improved the predictions significantly. A similar effect of subgrid heterogeneity, albeit to a lesser degree, was observed for denitrification, which also occurred in the columns. The results imply that subgrid heterogeneity is an important consideration in extrapolating reaction rates from the laboratory to field.

Establishing a geochemical heterogeneity model for a contaminated vadose zone--aquifer system.

A large set of sediment samples from a 1600 m² experimental plot within a 2.2 km² vadose zone and groundwater uranium (VI) plume was subject to physical, chemical, and mineralogic characterization. The plot is being used for field experimentation on U(VI) recharge and transport processes within a persistent groundwater plume that exists in the groundwater-river interaction zone of the Columbia River at the U.S. DOE Hanford site. The samples were obtained during the installation of 35 tightly spaced (10 m separation) groundwater monitoring wells. The characterization measurements for each sample included total contaminant concentrations (U and Cu primarily), bicarbonate extractable U(VI), sequential ²³8U(VI) contaminant desorption Kd, ²³³U(VI) adsorption K(d), grain size distribution, surface area, extractable poorly crystalline Fe(III) oxides, and mineralogy. The characterization objective was to inform a conceptual model of coupled processes controlling the anomalous longevity of the plume, and to quantify the spatial heterogeneity of the contaminant inventory and the primary properties effecting reactive transport. Correlations were drawn between chemical, physical, and reaction properties, and Gaussian simulation was used to compute multiple 3-D realizations of extractable U(VI), the ²³³U(VI) adsorption K(d), and the distribution of the reactive <2 mm fraction. Adsorbed contaminant U(VI) was highest in the vadose zone and the zone of seasonal water table fluctuation lying at its base. Adsorbed U(VI) was measureable, but low, in the groundwater plume region where very high hydraulic conductivities existed. The distribution of adsorbed U(VI) displayed no apparent correlation with sediment physical or chemical properties. Desorption [²³8U(IV)] and adsorption [²³³U(VI)] K(d) values showed appreciable differences due to mass transfer controlled surface complexation and the effects of long subsurface residence times. The ²³³U(VI) adsorption K(d), a combined measure of surface complexation strength and site concentration, was relatively uniform throughout the domain, displaying correlation with fines distribution and surface area. The characterization results revealed U(VI) supplied to the groundwater plume through spatially heterogeneous recharge from residual contamination in the zone of seasonal water table fluctuation, and transport of U(VI) controlled by weak, kinetically-controlled surface complexation in the coarse-textured saturated zone. Geostatistical relationships for the adsorbed contaminant U distribution in the characterization domain allow an extrapolation to inventory at the plume scale, a critical unknown for remedial action.

Persistence of uranium groundwater plumes: contrasting mechanisms at two DOE sites in the groundwater-river interaction zone.

We examine subsurface uranium (U) plumes at two U.S. Department of Energy sites that are located near large river systems and are influenced by groundwater-river hydrologic interaction. Following surface excavation of contaminated materials, both sites were projected to naturally flush remnant uranium contamination to levels below regulatory limits (e.g., 30 µg/L or 0.126 µmol/L; U.S. EPA drinking water standard), with 10 years projected for the Hanford 300 Area (Columbia River) and 12 years for the Rifle site (Colorado River). The rate of observed uranium decrease was much lower than expected at both sites. While uncertainty remains, a comparison of current understanding suggests that the two sites have common, but also different mechanisms controlling plume persistence. At the Hanford 300 A, the persistent source is adsorbed U(VI) in the vadose zone that is released to the aquifer during spring water table excursions. The release of U(VI) from the vadose zone and its transport within the oxic, coarse-textured aquifer sediments is dominated by kinetically-limited surface complexation. Modeling implies that annual plume discharge volumes to the Columbia River are small (

Rapid electron exchange between surface-exposed bacterial cytochromes and Fe(III) minerals.

The mineral-respiring bacterium Shewanella oneidensis uses a protein complex, MtrCAB, composed of two decaheme cytochromes, MtrC and MtrA, brought together inside a transmembrane porin, MtrB, to transport electrons across the outer membrane to a variety of mineral-based electron acceptors. A proteoliposome system containing a pool of internalized electron carriers was used to investigate how the topology of the MtrCAB complex relates to its ability to transport electrons across a lipid bilayer to externally located Fe(III) oxides. With MtrA facing the interior and MtrC exposed on the outer surface of the phospholipid bilayer, the established in vivo orientation, electron transfer from the interior electron carrier pool through MtrCAB to solid-phase Fe(III) oxides was demonstrated. The rates were 10(3) times higher than those reported for reduction of goethite, hematite, and lepidocrocite by S. oneidensis, and the order of the reaction rates was consistent with those observed in S. oneidensis cultures. In contrast, established rates for single turnover reactions between purified MtrC and Fe(III) oxides were 10(3) times lower. By providing a continuous flow of electrons, the proteoliposome experiments demonstrate that conduction through MtrCAB directly to Fe(III) oxides is sufficient to support in vivo, anaerobic, solid-phase iron respiration.

Electron transport at the microbe-mineral interface: a synthesis of current research challenges.

Many bacterial and archaeal species can couple growth to the respiratory reduction or oxidation of insoluble mineral oxides of transition metals. These solid substrates are abundant electron sinks and sources for life on Earth, but, since they are insoluble in water, they cannot enter the bacterial cells. So, to exploit these electron sinks and sources, specific respiratory electron-transfer mechanisms must overcome the physical limitations associated with electron transfer between a microbe and extracellular metal oxides. Recent microbiological, geochemical, biochemical, spectroscopic and structural work is beginning to shed light on the molecular mechanism and impacts of electron transfer at the microbe-mineral interface from a nanometre to kilometre scale. The research field is attracting attention in applied quarters from those with interests in nanowires, microbial fuel cells, bioremediation and microbial cell factories.

Analysis of structural MtrC models based on homology with the crystal structure of MtrF.

The outer-membrane decahaem cytochrome MtrC is part of the transmembrane MtrCAB complex required for mineral respiration by Shewanella oneidensis. MtrC has significant sequence similarity to the paralogous decahaem cytochrome MtrF, which has been structurally solved through X-ray crystallography. This now allows for homology-based models of MtrC to be generated. The structure of these MtrC homology models contain ten bis-histidine-co-ordinated c-type haems arranged in a staggered cross through a four-domain structure. This model is consistent with current spectroscopic data and shows that the areas around haem 5 and haem 10, at the termini of an octahaem chain, are likely to have functions similar to those of the corresponding haems in MtrF. The electrostatic surfaces around haem 7, close to the ß-barrels, are different in MtrF and MtrC, indicating that these haems may have different potentials and interact with substrates differently.

Molecular structure and free energy landscape for electron transport in the decahaem cytochrome MtrF.

The free energy profile for electron flow through the bacterial decahaem cytochrome MtrF has been computed using thermodynamic integration and classical molecular dynamics. The extensive calculations on two versions of the structure help to validate the method and results, because differences in the profiles can be related to differences in the charged amino acids local to specific haem groups. First estimates of reorganization free energies λ yield a range consistent with expectations for partially solvent-exposed cofactors, and reveal an activation energy range surmountable for electron flow. Future work will aim at increasing the accuracy of λ with polarizable forcefield dynamics and quantum chemical energy gap calculations, as well as quantum chemical computation of electronic coupling matrix elements.

The roles of CymA in support of the respiratory flexibility of Shewanella oneidensis MR-1.

Shewanella species are isolated from the oxic/anoxic regions of seawater and aquatic sediments where redox conditions fluctuate in time and space. Colonization of these environments is by virtue of flexible respiratory chains, many of which are notable for the ability to reduce extracellular substrates including the Fe(III) and Mn(IV) contained in oxide and phyllosilicate minerals. Shewanella oneidensis MR-1 serves as a model organism to consider the biochemical basis of this flexibility. In the present paper, we summarize the various systems that serve to branch the respiratory chain of S. oneidensis MR-1 in order that electrons from quinol oxidation can be delivered the various terminal electron acceptors able to support aerobic and anaerobic growth. This serves to highlight several unanswered questions relating to the regulation of respiratory electron transport in Shewanella and the central role(s) of the tetrahaem-containing quinol dehydrogenase CymA in that process.

Development of a proteoliposome model to probe transmembrane electron-transfer reactions.

The mineral-respiring bacterium Shewanella oneidensis uses a protein complex, MtrCAB, composed of two decahaem cytochromes brought together inside a transmembrane porin to transport electrons across the outer membrane to a variety of mineral-based electron acceptors. A proteoliposome system has been developed that contains Methyl Viologen as an internalized electron carrier and valinomycin as a membrane-associated cation exchanger. These proteoliposomes can be used as a model system to investigate MtrCAB function.

Mtr extracellular electron-transfer pathways in Fe(III)-reducing or Fe(II)-oxidizing bacteria: a genomic perspective.

Originally discovered in the dissimilatory metal-reducing bacterium Shewanella oneidensis MR-1 (MR-1), key components of the Mtr (i.e. metal-reducing) pathway exist in all strains of metal-reducing Shewanella characterized. The protein components identified to date for the Mtr pathway of MR-1 include four multihaem c-Cyts (c-type cytochromes), CymA, MtrA, MtrC and OmcA, and a porin-like outer membrane protein MtrB. They are strategically positioned along the width of the MR-1 cell envelope to mediate electron transfer from the quinone/quinol pool in the inner membrane to Fe(III)-containing minerals external to the bacterial cells. A survey of microbial genomes has identified homologues of the Mtr pathway in other dissimilatory Fe(III)-reducing bacteria, including Aeromonas hydrophila, Ferrimonas balearica and Rhodoferax ferrireducens, and in the Fe(II)-oxidizing bacteria Dechloromonas aromatica RCB, Gallionella capsiferriformans ES-2 and Sideroxydans lithotrophicus ES-1. The apparent widespread distribution of Mtr pathways in both Fe(III)-reducing and Fe(II)-oxidizing bacteria suggests a bidirectional electron transfer role, and emphasizes the importance of this type of extracellular electron-transfer pathway in microbial redox transformation of iron. The organizational and electron-transfer characteristics of the Mtr pathways may be shared by other pathways used by micro-organisms for exchanging electrons with their extracellular environments.

Redox reactions of reduced flavin mononucleotide (FMN), riboflavin (RBF), and anthraquinone-2,6-disulfonate (AQDS) with ferrihydrite and lepidocrocite.

Flavins are secreted by the dissimilatory iron-reducing bacterium Shewanella and can function as endogenous electron transfer mediators. To assess the potential importance of flavins in Fe(III) bioreduction, we investigated the redox reaction kinetics of reduced flavin mononucleotide (i.e., FMNH(2)) and reduced riboflavin (i.e., RBFH(2)) with ferrihydrite and lepidocrocite. The organic reductants rapidly reduced and dissolved ferrihydrite and lepidocrocite in the pH range 4-8. The rate constant k for 2-line ferrihydrite reductive dissolution by FMNH(2) was 87.5 ± 3.5 M(-1)·s(-1) at pH 7.0 in batch reactors, and k was similar for RBFH(2). For lepidocrocite, k was 500 ± 61 M(-1)·s(-1) for FMNH(2) and 236 ± 22 M(-1)·s(-1) for RBFH(2). The surface area normalized initial reaction rates (r(a)) were between 0.08 and 77 µmol·m(-2)·s(-1) for various conditions in stopped-flow experiments. Initial rates (r(o)) were first-order with respect to iron(III) oxide concentration, and r(a) increased with decreasing pH. Poorly crystalline 2-line ferrihydrite yielded the highest r(a), followed by more crystalline 6-line ferrihydrite and crystalline lepidocrocite. Compared to a previous whole-cell study with Shewanella oneidensis strain MR-1, our findings suggest that the reduction of electron transfer mediators by the Mtr (i.e., metal-reducing) pathway coupled to lactate oxidation is rate limiting, rather than heterogeneous electron transfer to the iron(III) oxide.