Modeling the iron storage protein ferritin reveals how residual ferrihydrite iron determines initial ferritin iron sequestration kinetics

Computational models can be created more efficiently by composing them from smaller, well-defined sub-models that represent specific cellular structures that appear often in different contexts. Cellular iron metabolism is a prime example of this as multiple cell types tend to rely on a similar set of components (proteins and regulatory mechanisms) to ensure iron balance. One recurrent component, ferritin, is the primary iron storage protein in mammalian cells and is necessary for cellular iron homeostasis. Its ability to sequester iron protects cells from rising concentrations of ferrous iron limiting oxidative cell damage. The focus of the present work is establishing a model that tractably represents the ferritin iron sequestration kinetics such that it can be incorporated into larger cell models, in addition to contributing to the understanding of general ferritin iron sequestration dynamics within cells. The model’s parameter values were determined from published kinetic and binding experiments and the model was validated against independent data not used in its construction. Simulation results indicate that FT concentration is the most impactful on overall sequestration dynamics, while the FT iron saturation (number of iron atoms sequestered per FT cage) fine tunes the initial rates. Finally, because this model has a small number of reactions and species, was built to represent important details of FT kinetics, and has flexibility to include subtle changes in subunit composition, we propose it to be used as a building block in a variety of specific cell type models of iron metabolism.