The aim of this model is to analyze iron uptake by the introduced Pft ferritin within the statolith. The analysis would demonstrate (i) the expected rate of uptake and (ii) the expected amount of loaded iron at equilibrium (ie. once steady state is reached).
Using chemical kinetics equations for Langmuir adsorption adapted for single-pore diffusion,
Total intra-pore sites = Vacant sites + Sum of sites occupied by iron
[Iron]external + Intra-pore site → [Iron]internal
The rate laws are given as first order
Forward rate = k[Iron][Vacant Sites]
Reverse rate = k[Internalized iron]
Using these equations, we can develop steady-state profiles for internalized iron concentration as a concentration of bulk iron concentrations.
We see from the following plot that as external iron concentration increases, internalized iron concentration increases. The maximum iron load is approximately 80 M Fe2+ (4000 atoms in an ~8nm diameter spherical cage). Although the iron loaded within the ferritin at any given time has a higher concentration than the bulk, this is somewhat expected since biologically, ferritin serves as an iron storage molecule. The large disparity is due to a limitation in our current preliminary model, where bulk ferritin concentration is considered constant.
We classify this as chemisoprtion. Technically, the overall process includes physisorption for the limiting diffusion step and chemisorption is the negligible redox step that follows.
The parameters could be improved through experimental determinations of iron loading rate . Iron uptake and clearance kinetics profiles could indicate whether our system follows the kinetic trend described by the model or not.
The parameters could also be improved through analysis of plants grown in various concentrations of iron: this would yield the equilibrium iron stored within the plant as a function of iron media concentration
The aim of this model is to analyze the diffusion of essential nutrients and iron into arabidopsis roots from various plant growth mediums. The goal of this model is to compare the transport efficacy of various mediums. Specifically, our goal was to compare diffusion rates from the medium to arabidopsis roots in recyclable mediums, which we learned from our discussion with NASA scientist Dr. Gioia Massa are preferred for plant growth in space.
We can calculate concentration profiles using Fick’s Second Law. Ignoring all terms aside from radial diffusion, the concentration profile in each medium is given by:
At the medium-root interface, we know that the fluxes are equal. By Fick’s First Law:
As in the prior model, experimental work to determine diffusion coefficients would improve our existing model significantly. We could also seek to consider a more diverse set of media and nutrients.
The aim of this model is to take the ideal parameter for current determined by solving an ODE and find the combination of passives and NPN devices that will work for generation of a constant current in that range. Ideally, the output current of the device should be tunable via a 5 KOhm rheostat.
Unfortunately, CircuitLab is closed-source software but does state that it uses “SPICE-like component models” to give “accurate results for nonlinear circuit effects”. This is most important for modeling of the transistors, where nonlinear behavior is responsible for changes in gain current. SPICE uses a Gummel-Poon model to more accurately model changes in transistor gain over different base currents.
Initially, an adjustable Zener voltage reference was modeled as a zener diode with a set reverse voltage and a resistor pair in parallel as described in Texas Instruments’ typical application schematic for the LM4041 adjustable Zener voltage reference. A single NPN power transistor (D44VH) was used to control the output current.
Simulations were performed with a 200 mOhm load as an overestimate in order to supply the requisite 10A.
The circuit utilizing the Zener voltage reference was a valid result, but was very susceptible to changes in the load resistance and had rigid requirements for passive component values to prevent overcurrenting the Zener diode. Furthermore, due to the low-resistance series resistors lots of power was lost as heat in the constant current generation side of the circuit.
A Zener diode was chosen to replace the Zener voltage reference, and the resistor pair in parallel with it was removed. Less power was lost as heat when an NPN darlington was used to drive the D44VH. This means that the overall gain is very significantly amplified, and even with small currents of 10 mA through the MCP4642 we can modulate the output current between:
The circuits for the primary and secondary coil are identical except that the desired output current is set by a chassis mount resistor.
The resistor values are as follows:
The result of these changes is that the circuit is much more efficient (it takes advantage of the additional gain provided by a darlington pair) and more adjustable. The only device that has to be changed to grossly alter the coil current is the chassis mount resistor, and lots of fine adjustment is achievable by altering the tap fraction of the MCP4642. We went on to perform DC sweep simulations for this tap fraction (between 0 and 1) and plotted the resulting load current for the primary and secondary coils.
Circuit for constant current driver with 220 mOhm chassis resistor (R4) installed.
Dependency of load current (Amperes) on tap fraction of the MCP4642 (R10.K) for a 220 mOhm chassis resistor. The relationship is approximately linear (y = 3.7646x + 6.8596, R2 = 0.9848) and the current can vary between 7.117 to 10.928 A. This provides evidence that our revised circuit could meet the calculated theoretical current (9.87 A) and can still be adjusted to a wider range (+ 1.058A, -2.753A). This would allow us to compensate for real-world inaccuracies when building the constant current driver and still achieve the desired 9.87A.
Circuit for constant current driver with 3.3 Ohm chassis resistor (R4) installed.
Dependency of load current (Amperes) on tap fraction of the MCP4642 (R10.K) for a 220 mOhm chassis resistor. The relationship is still approximately linear (y = 0.0146x + 1.0347, R2 = 0.9723) but a second order polynomial fit the data better (y = 0.009x2 + 0.0056x + 1.0362, R² = 0.9974). This displays that for lower currents about 1A, the constant current driver’s rheostat-based current adjustment may not be able to adjust the current quite effectively. For our application, the accuracy and modulability of the current used to drive the secondary coil is less consequential than the primary. However, a separate selection of transistors may be desired for more accurate low-current operation in other projects.
The goal of this model is to get a better understanding of the mechanisms of iron loading in ferritin, as well as to estimate values such as diffusion constants. We intend to use molecular simulation packages such as Pyrosetta along with tools such as AlphaFold2. This model is still in very early stages as these simulations are very resource and time intensive.
We intend to continue developing this model with the following assumptions:
There are no explicit equations we are directly solving. We intend to use the pre-existing packages of pyRosetta, which are primarily Monte Carlo in nature, and AlphaFold2, which is a machine-learning algorithm
For the purpose of calculating diffusion through ferritin's ion channels, we intend on simplifying the system to a single ion-channel and it's immediate surroundings. Simulating the entire ferritin is unfeasible due to it's molecular size.
As this model is still in early development, we do not have significant results beyond evidence that different ferritin subunits are able to dock with one another
The following images are acquired using the OVITO software