Analysis of the application of the generalized monod kinetics model to describe the human corneal oxygen-consumption rate during soft contact lens wear

A three-dimensional model to describe complete human corneal oxygenation during contact lens wear

We perform a novel 3D study to quantify the corneal oxygen consumption and diffusion in each part of the cornea with different contact lens materials. The oxygen profile is calculated as a function of oxygen tension at the cornea-tear interface and the oxygen transmissibility of the lens, with values used in previous studies. We aim to determine the influence of a detailed geometry of the cornea in their modeling compared to previous low dimensional models used in the literature. To this end, a 3-D study based on an axisymmetric volume element analysis model was applied to different contact lenses currently on the market. We have obtained that the model provides a valuable tool for understanding the flux and cornea oxygen profiles through the epithelium, stroma, and endothelium. The most important results are related to the dependence of the oxygen flux through the cornea-lens system on the contact lens thickness and geometry. Both parameters play an important role in the corneal flux and oxygen tension distribution. The decline in oxygen consumption experienced by the cornea takes place just inside the epithelium, where the oxygen tension falls to between 95 and 16 mmHg under open eye conditions, and 30 to 0.3 mmHg under closed eye conditions, depending on the contact lens worn. This helps to understand the physiological response of the corneal tissue under conditions of daily and overnight contact lens wear, and the importance of detailed geometry of the cornea in the modeling of diffusion for oxygen and other species.

A refined model on flow and oxygen consumption in the human cornea depending on the oxygen tension at the interface cornea/post lens tear film during contact lens wear

The study of oxygen consumption rate under" in vivo" human cornea during contact lens wear has been technically a challenge and several attempts have been made in the last 20 years to model the physiology of the human cornea during contact lens wear. Unfortunately, some of these models, based on a constant corneal oxygen consumption rate, produce areas on the cornea where the oxygen tension is negative, which has no physical sense. In order to avoid such inconsistency, different researchers have developed alternative models of oxygen consumption, which predict the likely oxygen metrics available at the interface cornea/post lens tear film by determination of oxygen flux, oxygen consumption, and oxygen tension through the different layers (endothelium, stroma, and epithelium). Although oxygen deficiency produces corneal edema, corneal swelling, hypoxia, acidosis, and other abnormalities, the estimation of the oxygen distribution below the impact of a contact lens wear is interesting to know which lens transmissibility was adequate to maintain the cornea and avoid epithelial and stromal anoxia. The estimation of minimum transmissibility for a lens for extended wear applications will be very useful for both clinicians and manufacturers. The aim of this work is to present a complete discussion based on Monod kinetics model that permits give an estimation of oxygen partial pressure distribution, the profile distribution of corneal flux and oxygen consumption rate, and finally the estimation of the relaxation mechanism of the cornea depending on the oxygen tension at the interface cornea/post lens tear film. Relaxation time in this context can quantify the capability of the corneal tissue to adapt to increasing concentrations of oxygen. It is proposed this parameter as a biological meaningful indicator of the interaction between contact lens polymers and living tissues such as the corneal cellular layer.

Simulation of time-fractional oxygen diffusion in cornea coated by contact-lens

In this paper, the time-fractional oxygen diffusion has been simulated in a one-dimensional (1D) corneal-contact lens (CL) system. Different CLs have been employed as Balafilcon, thin- and thick-Polymacon. It is assumed that homogeneous and isotropic porous mediums of cornea and CL is saturated with compressible oxygen. The computations of the time-fractional derivations are done based on the Caputo method. The obtained results show that the fractional derivative order (FDO) severely affects pressure distribution in cornea and CL. Consequently, the magnitudes of post-lens-tear-film (PoLTF) pressure change due to diverse FDOs. Particularly, maximum changes have been observed in the results gained from the CLs with thicknesses more than 100 μm. The agreement of the results obtained from the time-fractional modeling with the experimental data compared to the standard diffusion modeling has been improved by more than 36%. Finally, it has been demonstrated that high-thickness CLs can cause exist anomalous diffusion process in cornea tissue.

Corneal relaxation time estimation as a function of tear oxygen tension in human cornea during contact lens wear

The purpose is to estimate the oxygen diffusion coefficient and the relaxation time of the cornea with respect to the oxygen tension at the cornea-tears interface. Both findings are discussed. From the experimental data provided by Bonanno et al., the oxygen tension measurements in vivo for human cornea-tears-contact lens (CL), the relaxation time of the cornea, and their oxygen diffusion coefficient were obtained by numerical calculation using the Monod-kinetic model. Our results, considering the relaxation time of the cornea, observe a different behavior. At the time less than 8 s, the oxygen diffusivity process is upper-diffusive, and for the relaxation time greater than 8 s, the oxygen diffusivity process is lower-diffusive. Both cases depend on the partial pressure of oxygen at the entrance of the cornea. The oxygen tension distribution in the cornea-tears interface is separated into two different zones: one for conventional hydrogels, which is located between 6 and 75 mmHg, with a relaxation time included between 8 and 19 s, and the other zone for silicone hydrogel CLs, which is located at high oxygen tension, between 95 and 140 mmHg, with a relaxation time in the interval of 1.5-8 s. It is found that in each zone, the diffusion coefficient varies linearly with the oxygen concentration, presenting a discontinuity in the transition of 8 s. This could be interpreted as an aerobic-to-anaerobic transition. We attribute this behavior to the coupling formalism between oxygen diffusion and biochemical reactions to produce adenosine triphosphate. © 2019 Wiley Periodicals, Inc. J Biomed Mater Res Part B: Appl Biomater 108B:14-21, 2020.