A Mathematical Model of Aqueous Humor Production and Composition

Purpose: We develop a mathematical model that predicts aqueous humor (AH) production rate by the ciliary processes and aqueous composition in the posterior chamber (PC), with the aim of estimating how the aqueous production rate depends on the controlling parameters and how it can be manipulated. Methods: We propose a compartmental mathematical model that considers the stromal region, ciliary epithelium, and PC. All domains contain an aqueous solution with different chemical species. We impose the concentration of all species on the stromal side and exploit the various ion channels present on the cell membrane to compute the water flux produced by osmosis, the solute concentrations in the AH and the transepithelial potential difference. Results: With a feasible set of parameters, the model predictions of water flux from the stroma to the PC and of the solute concentrations in the AH are in good agreement with measurements. Key parameters which impact the aqueous production rate are identified. A relevant role is predicted to be played by cell membrane permeability to \(\text{K}^+\) and \(\text{Cl}^-\), by the level of transport due to the Na+-H+ exchanger and to the co-transporter of Na+/K+/2Cl-; and by carbonic anhydrase. Conclusions: The mathematical model predicts the formation and composition of AH, based on the structure of the ciliary epithelium. The model provides insight into the physical processes underlying the functioning of drugs that are adopted to regulate the aqueous production. It also suggests ion channels and cell membrane properties that may be targeted to manipulate the aqueous production rate.