Team:Bielefeld-Germany/Modelling/Oxidation

From 2013.igem.org

Revision as of 02:49, 29 October 2013 by AnnaK (Talk | contribs)


Modelling - Mediator Oxidation

Mediator Oxidation

In a third electrochemical reaction the reduced mediator is regenerated at the electrode.This electrochemical oxidation at the anode surface occurs as shown in equation:

Bielefeld-germany-model-oxid-reaction-1.PNG


,where Mred is the reduced mediator and
Mox the oxidized mediator.
Then the current output can be calculated based on formula according to the Faraday's law:

Bielefeld-germany-model-oxid-reaction-2.PNG



,where I is the current density [A]
[Mred] is the concentration of reduced mediator in the chamber
n is the number of electrons taking part in the electrode reaction,
F is the Faradays constant (96 500 C) and
k3 is the reaction rate, mentioned above.


Simulation

The value k3 obtained as mentioned above was applied in the simulation of the current output in the final reaction of the three-reaction model. The start concentrations of the NAD+ and oxidized mediator were set as in the simulations preformed for the proceeding reactions at 100 µM. The simulation time span was set for 60 sec. The resulting plot for the mediator MB is shown in the Figure 1, below:

Figure 1: Curves for the current output, fluxes of concentration of the intermediates. Simulation performed for the start concentration of both NAD+ and oxidized mediator MB at 100 µM.



Further simulation has been performed for four different start concentrations of the oxidised mediator in order to investigate how varying start concentration influence the current output. The concentrations were set at 10 µM, 50 µM, 100 µM and 500 µM. The time span was set at 30 sec. The curves of the current outputs are presented in the Figure 2.

Figure 2: Curves for the current output and fluxes of concentration of the intermediates. Simulation performed for the varying start concentration of the oxidized mediator MB [Mox] .



The MATLAB source code for both simulations can be obtained here and here.


Analogous simulations have been performed for the mediator NMP. As the reaction rate k3 specific for this mediator could not be accessed, the value specific for the mediator MB was applied. The diagrams presenting the results of those simulations are shown in the Figure 3 and 4 respectively.

Figure 3: Curves for the current output, fluxes of concentration of the intermediates. Simulation performed for the start concentration of both NAD+ and oxidized mediator NMP at 100 µM.



Figure 4: Curves for the current output and fluxes of concentration of the intermediates. Simulation performed for the varying start concentration of the oxidized mediator NMP [Mox] .


The according .m files with MATLAB source code can be downloaded from here and here.

Discussion

The curve for the current output simulated for the mediator MB has been presented in the figure 1. It could be shown that the maximal output is reached rapidly after the start of the reaction system. Immediately after reaching the maximum, it decreases equally fast. At the 20 second mark it goes in the stable phase of about 20 per cent of the maximal current output.
The comparison of different start concentrations of the oxidized mediator can be obtained from the figure 2. The curve form is identical for all four concentrations and the reached current values differ only slightly for the different simulations. Also the maximum output is almost identical and is reached at the same time mark for all four approaches. Interestingly the curves obtained from the simulation for the mediator NMP are very similar to those for the mediator MB. This points that the reaction rate k3 for the oxidation of the mediator at the anode plays a predominant role in generating current.
Obviously there is necessity to constantly provide the system with the redox molecules in order to sustain the current at the high level.

Figure 5: .




References







Contents