Mathematical modeling of aqueous humor flow and intraocular pressure under uncertainty: towards individualized glaucoma management

Supplementary Files

Figure 1.

How to Cite

Szopos M, Cassani S, Guidoboni G, Prud’homme C, Sacco R, Siesky B, Harris A. Mathematical modeling of aqueous humor flow and intraocular pressure under uncertainty: towards individualized glaucoma management. MAIO [Internet]. 2016 Dec. 15 [cited 2022 Jun. 25];1(2):29-3. Available from:

Copyright notice

Authors who publish with this journal agree to the following terms:

  1. Authors retain copyright and grant the journal right of first publication, with the work twelve (12) months after publication simultaneously licensed under a Creative Commons Attribution License that allows others to share the work with an acknowledgement of the work’s authorship and initial publication in this journal.

  2. After 12 months from the date of publication, authors are able to enter into separate, additional contractual arrangements for the non-exclusive distribution of the journal's published version of the work (e.g., post it to an institutional repository or publish it in a book), with an acknowledgement of its initial publication in this journal.


aqueous humor flow; glaucoma management; intraocular pressure; mathematical modeling; sensitivity analysis


Purpose: The aim of the proposed analysis is to provide both a qualitative description and a quantitative assessment of how variations in aqueous humor (AH) flow parameters influence intraocular pressure (IOP) and the outcome of IOP-lowering medications.

Methods: We developed a mathematical model that describes the steady-state value of IOP as the result of the balance between AH production and drainage. We performed stochastic simulations to assess the influence of different factors on the IOP distribution in ocular normotensive and ocular hypertensive subjects and on the IOP reduction following medications.

Results: The distribution of the relative frequency of a given IOP value for ocular normotensive subjects fits a right-skewed Gaussian curve with a frequency peak of 25% at 15.13 mmHg and a skewness of 0.2, in very good agreement with the results from a population-based study on approximately 12,000 individuals. The model also shows that the outcomes of IOP-lowering treatments depend on the levels of pre-treatment IOP and blood pressure. The model predicts mean IOP reductions of 2.55 mmHg and 4.31 mmHg when the pre-treatment IOP mean values are 15.13 mmHg and 20.12 mmHg, respectively; these predictions are in qualitative and quantitative agreement with clinical findings.

Conclusion: These findings may help identify patient-specific factors that influence the efficacy of IOP-lowering medications and aid the development of novel, effective, and individualized therapeutic approaches to glaucoma management.