Development of a 3D finite element model of lens microcirculation
© Vaghefi et al.; licensee BioMed Central Ltd. 2012
Received: 26 July 2012
Accepted: 10 August 2012
Published: 19 September 2012
It has been proposed that in the absence of a blood supply, the ocular lens operates an internal microcirculation system. This system delivers nutrients, removes waste products and maintains ionic homeostasis in the lens. The microcirculation is generated by spatial differences in membrane transport properties; and previously has been modelled by an equivalent electrical circuit and solved analytically. While effective, this approach did not fully account for all the anatomical and functional complexities of the lens. To encapsulate these complexities we have created a 3D finite element computer model of the lens.
Initially, we created an anatomically-correct representative mesh of the lens. We then implemented the Stokes and advective Nernst-Plank equations, in order to model the water and ion fluxes respectively. Next we complemented the model with experimentally-measured surface ionic concentrations as boundary conditions and solved it.
Our model calculated the standing ionic concentrations and electrical potential gradients in the lens. Furthermore, it generated vector maps of intra- and extracellular space ion and water fluxes that are proposed to circulate throughout the lens. These fields have only been measured on the surface of the lens and our calculations are the first 3D representation of their direction and magnitude in the lens.
Values for steady state standing fields for concentration and electrical potential plus ionic and fluid fluxes calculated by our model exhibited broad agreement with observed experimental values. Our model of lens function represents a platform to integrate new experimental data as they emerge and assist us to understand how the integrated structure and function of the lens contributes to the maintenance of its transparency.
The microcirculation model was initially based on a combination of electrical impedance measurements and theoretical modelling [8–11]. The observation that the measured ionic currents were directed inward at the poles and outward around the equator [Figure 1A, led to the suggestion that these currents represent the external portion of a circulating ionic current that drives the internal microcirculatory system within the lens. Briefly, the working model is that the current, carried primarily by Na+, enters at all locations around the lens via the extracellular space between fibre cells. Na+ eventually crosses the fibre cell membranes, and then flows from cell-to-cell towards the surface via an intracellular pathway mediated by gap junction channels [Figure 1B. The gap junction coupling conductance in the outer shell of differentiating fibres is concentrated at the equator [12, 13]. Hence, the intracellular current is directed to the equatorial epithelial cells. Here, the highest densities of Na+/K+ pumps are located to actively transport Na+ out of the lens . Thus, the intracellular current effluxes are highly concentrated at the equator causing the net current to be outward. At the poles, there is very little intracellular current so the net current is predominantly inward, along the extracellular spaces [Figure 1B. The driving force for these fluxes is hypothesized to be the difference in the electromotive potential of surface cells that contain Na+/K+ pumps and K+-channels, and inner fiber cells that lack functional Na+/K+ pumps and K+-channels and whose permeability is dominated by non-selective cation and Cl- conductances . This electrical connection together with the different membrane properties of the surface and inner cells causes the standing current to flow. In this model, the circulating current creates a net flux of solute that in turn generates fluid flow. The extracellular flow of water convects nutrients towards the deeper lying fiber cells, while the intracellular flow removes wastes and creates a well-stirred intracellular compartment.
Thus in this model, it is the circulating Na+ current that generates a circulation of fluid inside the lens. It is important to note that while the existence of the circulating ionic currents are firmly based on existing experimental data, circulating fluid flows in the lens have proven more difficult to measure. At present, these fluxes are only predicted to occur from indirect measurements and models of the measured electrical properties of the lens. The current model of the lens fluid dynamics is based on an equivalent circuit analysis of the microcirculation system [3, 16]. This analytical model inherently relied on approximate solutions and simplification of the underlying physics. To improve our understanding of the circulation system, here we have adopted a finite element modelling (FEM) approach. Using FEM, we have developed a 3D computer model of the microcirculation system that encapsulates the complex interplay of its important features and can also be solved numerically. This computer-based modelling approach allowed structural features such as fiber cell orientation, extracellular space dimensions and gap junction distribution to be included in the model. It also included functional information on the spatial differences in membrane permeability between surface and inner cells, thought to drive the circulating currents.
In this paper we describe our expansion of the equations that govern ion and fluid dynamics in the lens [2, 3, 10, 15–18] to 3D, and their subsequent implementation into a new FEM that encapsulates known structural and functional parameters of the mouse lens. This model is based on our continuous imaging and modelling iterative investigation into the fluid dynamics of the lens [19–25]. To test the ability of this computer model to reproduce the functional properties of the lens, we have used a series of experimentally derived boundary conditions [26–28] to allow it to be solved. We show that our model is capable of predicting the experimentally measured steady state lens properties and to generate circulating ion and water fluxes. Hence we believe that our computer model is a useful tool to study how lens structure and function influence its optical properties.
In this section, we first present the derivation of the fundamental mathematical equations, originally formulated by Mathias et al. to describe the microcirculation system [2, 3, 10, 15–18]. We then develop a computer mesh to represent the structure of the mouse lens to enable these equations to be implemented in 3D. Next we implement the model using the C++ programming language and solve the model using experimentally derived boundary conditions [26–28]. The resultant 3D model calculates ion concentration gradients, membrane potential gradients, and intra- and extracellular ion and water fluxes in different regions of the lens.
Derivations of general equations
Glossary of symbols used in this manuscript
extracellular cleft width
maximum Na/K ATPase pump current density
Na/K ATPase pump current density
conductivity per membrane area
surface hydraulic permeability
intercellular hydraulic permeability
membrane solute permeability
equatorial radius of the model
posterior radius of the model
anterior radius of the model
fluid mass source
permittivity of vacuum
zeta potential (cell membrane potential)
Fluid flow parameters
8.314 × 103
4.0 × 10−8
4.0 × 10−8
2.7 × 10−7
3.05 × 10−14
1.45 × 10−8
Solute transport properties
1.380 × 10−11
1.6 × 10−10
9.648 × 104
Free solution/cytoplasm Na+ diffusion1
1.39 × 10−6
Free solution/cytoplasm K+ diffusion1
2.04 × 10−6
Free solution/cytoplasm Cl- diffusion1
2.12 × 10−6
Intracellular Na+ diffusion
1.39 × 10−8
Intracellular K+ diffusion
2.04 × 10−8
Intracellular Cl- diffusion
2.12 × 10−8
Extracellular Na+ diffusion
1.39 × 10−6
Extracellular K+ diffusion
2.04 × 10−6
Extracellular Cl- diffusion
2.12 × 10−6
Na+ fibre cell membrane conductivity2
Cl- fibre cell membrane conductivity2
K+ surface membrane conductivity2
Na+/K+ pump max. pump rate
Na+/K+ pump 1/2 max Na+ concentration
Na+/K+ pump 1/2 max K+ concentration
We solved these equations [Eq. 1 to Eq. 3] for the movement of ions and water in the extracellular and intracellular spaces of the lens. These two domains were linked by cell membranes which could be crossed by ions and water from one space to another [Figure 1B]. Since the water and ions enter the lens from the extracellular space, we started with these equations.
Extracellular flux equations
The fluid flow in the extracellular clefts can be partially described by the Stokes flow equations [Eq. 1 & Eq. 2. Another part of the extracellular fluid fluxes has been shown to be due to electro-osmosis . Electro-osmosis is due to the osmotic gradient created by uneven charge distribution in a fluid affected by an electric field. It has been shown that this osmosis is essential to the modelling of extracellular fluxes in the ocular lens .
The parameters and their units are listed in [Table 1. We didn’t consider electro-osmosis in the intracellular space portion of our model . This was since the electric field gradient across cell cytoplasm was considered negligible. We then coupled the extracellular water fluxes with the ionic fluxes in this domain.
Extracellular ion fluxes
We computationally solved the advective Nernst-Plank equation [Eq. 3] for the lenticular extracellular ion fluxes. We modelled the coupled water and ion extracellular fluxes throughout the 3D model. These fluxes could then become trans-membranous flows at any point of the model, crossing the modelled cell membranes.
Trans-membrane flux equations
The parameters and their units are listed in [Table 1]. In the above equations, E is the Nernst potential. At the Nernst potential, there is no net flow of ions across the cell membrane through the channels. Hence, in our model we linked the trans-membrane water and ion fluxes by ionic concentrations (i.e. osmosis) and membrane potential (i.e. Nernst potential). After crossing the membrane, we treated the water and ion fluxes as intracellular fluxes.
Intracellular flux equations
We modelled the fluxes in the intracellular space to pass through a mixture of cells cytoplasm and membranes . This was due to the large size of each element in our model compared to the lens cell volumes. It was impractical to discretely model the intracellular flow through each cell cytoplasm and membrane in a given element. Instead, we homogenized the flow equations to obtain one formula, describing the net cytoplasmic and membranous fluxes through every element.
Here ϕ is the angular degree from the equator (i.e. ϕ = 0° at equator, ϕ = 90° at anterior pole and ϕ = −90° at posterior pole). It has been determined that G max = 0.6mS/cm in the radial orientation, and in the other directions G = 0.23mS/cm . We modelled the intracellular ionic fluxes to be directed by the above equations and the regional distribution of gap junctions [Eq. 15 - Eq. 17. According to the microcirculation theory, these ionic fluxes then moved towards the periphery of the lens accompanied by water flows until they reach the boundary of the model.
Surface flux equations
The parameters and their units are listed in [Table 1] and the subscript “o” points to the outside (i.e. boundary) conditions. We considered the surface fluxes as the output of the system, while the extracellular fluxes were its input. For an incompressible stable system, like the current lens model, the input and output levels should equate at all times and this conservation of mass was controlled for as part of the Stokes equations [Eq. 1].
We assumed the initial condition of the system to be electrically neutral. Hence, the equation above ensures the electro-neutrality of the model at any point of time. We applied this equation over all the ion species modelled here. This was based on the assumption that the ion species not modelled (e.g. calcium) had no significant influence on electro-neutrality of the model.
Solving the model
Finite Element Mesh Creation
Finite element mesh specifications for the current model
Number of nodes
Number of elements
Initial conditions at outer lens boundary for the present model
Extracellular Na+ concentration
Extracellular K+ concentration
Extracellular Cl- concentration
Intracellular Na+ concentration
Intracellular K+ concentration
Intracellular Cl- concentration
We defined the maximum change of concentrations between two consecutive iterations (C n+1 - C n ) in all the elements as the convergence parameter. We observed that using the above boundary and initial conditions, this error was decreased with each iteration cycle. We set the model’s convergence criterion to be less than 5 mM. We believed that level of iteration error was adequate, since most of the solution fluctuations were caused by [K + ] i and [Na + ] e fields and the enforced iteration error threshold was less than 5% of the initial field value.
We solved the 3D model on our high performance computer (HPC) at the Auckland Bioengineering Institute (ABI). Our mainframe was comprised of an IBM server with 64 processors of 1.9 GHz calculation speed, peak performance of 306.21Gfps and 256 GB of physical memory.
Displaying the data
Text formatted file of the model’s computed fields were linked via JAVA programming language format to CMGUI ( http://www.cmgui.org) [48, 49], an advanced 3D visualisation software package with modelling capabilities. CMGUI was used for model visualisation and manipulation and allowed for automated the scaling and false colouring of the data. The following presentations were all created using this method.
Using the boundary conditions listed in [Table 5], we solved the model and generated 3D maps of standing fields of intracellular and extracellular ion concentrations, electrical potentials and pressure, plus circulating ionic and water fluxes. In this section we first use quarter section views of the 3D model [Figure 3C] to represent regional differences in standing electrochemical and pressure fields, allowing these predicted properties to be compared to experimentally derived values. Then we use the full model view to generate 3D vector maps that visualize the predicted ion and water fluxes in the lens for the first time.
Standing electrochemical and pressure fields
Comparison of modelled standing electrochemical gradients outcomes and the existing empirical data
3D model radial range
Experimental radial range
3D model averaged
6.9 to 16
5 to 16 1
13.8 3, 21.5 4
91 to 100
10 to 12
16.8 5, 12.3 5
−57 to −64
−50 to −70 2
0.0 to −27
0 to −36 2
0.0 to 10
0.0 to 48 7
Comparison of modelled circulating currents and the existing experimental mouse data
Net current density (μA/cm2)
22 1, 20 ± 2.6 2, outward
26 1, 32 2, inward
7 1, 42 2, inward
Fluid velocity (nm/s)
Vequator – intracellular
Vanterior pole – extracellular
Vposterior pole – extracellular
It has been proposed that the circulating currents observed experimentally in lenses from a variety of species [5, 14, 41, 52] drive an internal microcirculation system. In the absence of a blood supply the microcirculation delivers nutrients and removes metabolic waste products from inner fibre cells, while maintaining ionic homeostasis . The circulating fluxes are thought to be the net result of spatial differences in the location of ion channels and transporters that determine local membrane permeability, and the density of gap junction channels that direct intercellular fluxes within the lens . We first created a 3D mesh of the lens structure. We then implemented a series of equations that describe the local transport properties of cells in different regions in the lens. Finally we solved these equations using a FEM approach at each location of the model. In summary, we have produced a computer model of the lens that not only accurately predicts standing ionic concentrations [Figure 4 and electrical [Figure 5 gradients and the existence of a pressure gradient [Figure 6, but also produces the first 3D vector maps of predicted current [Figure 8 and fluid [Figure 10 fluxes inside the lens. The observed agreement between experimentally measured values and those calculated by our model [Table 6 & Table 7, suggest that our computer model is mimicking lens physiology and generates a microcirculation. However, the underestimation of the magnitude of the intracellular pressure gradient highlights the fact that the model is only as robust as its underlying assumptions and will require additional refinement and revision as new experimental data on lens structure and function becomes available. In the following sections we discuss the assumptions and limitations of our current model with the view to highlight areas where further refinements of the model are required.
The solute flux in a fluid is governed by diffusion, electro-diffusion (if the solute is charged) and advection. Diffusion is the random walk of particles due to Brownian motion. Electro-diffusion is the flux of a charged particle due to the force applied by an electric field. Advection is the transport of a solute by a fluid that is moving. The physical origins of these transport processes and the derivations of the equations have been discussed previously in the literature [31, 32]. The implementation of these equations in a computational platform is very well explained elsewhere. In this study, we implemented the driven equations and solved for a set of converged-upon 3D fields. This method produced the calculated standing ionic concentration gradients in the lens predicted by our model [Figure 4.These gradients in turn gave rise to a trans-membrane electric-potential gradient field [Figure 5, based on the Hodgkin–Huxley model. Consider a cell membrane that is not equally permeable to all the present ionic species on either side of it. Here, the permeable ions will tend to move down their concentration gradient taking their electrical charge with them as they go. Therefore, an electrical potential will be generated, which will drive them in the reverse direction. An electrochemical equilibrium will be reached when the diffusive force equals the electromotive force.
Our model has calculated a hydrostatic pressure gradient [Figure 6, the existence of which has recently been confirmed experimentally . While the orientations of the calculated and measured pressure gradients are similar, they differ in magnitude (48kPa compared to 10kPa). It has been proposed that this pressure gradient is generated by the restricted flow of water from the centre to the periphery of the lens by a pathway mediated by gap junction channels; since genetic manipulations to increase or decrease gap junction numbers produces inverse changes in the pressure gradient. This illustrates that structural components of the lens can influence the magnitude of the pressure gradient. This suggests that the difference between calculated and measured pressure fields may reflect the absence of a structural feature not currently captured in our model. In this regard we have recently identified a zone in the inner cortex of the lens that exhibits a reduction in the penetration of solutes and water [19, 22] that could influence the magnitude of the calculated pressure gradient. Future updates of the model will include such newly discovered structural elements, allowing their effect on the calculated pressure fields to be assessed.
The electrochemical fields, combined with the hydrostatic pressure gradients [Figure 6 in our model generate the circulating ionic currents throughout the lens [Figure 8 & Figure 9. The ionic fluxes are accompanied by water flows in the microcirculation system. The water fluxes are generally described by the Navier–Stokes equations, which are derived from the conservation of mass, momentum, and energy principals. Also, the fluid flow in the lens can be described as slow or low-Reynolds number flow. It is also reasonable to assume that the fluid flow in a normal lens is near or at steady-state at all times . These simplifications reduce the Navier–Stokes equations to the Stokes equations. The derivation and implementation of these equations into a computational framework is discussed elsewhere [16, 18]. In our model, solving this set of equations results in the calculation of the water flow velocity fields [Figure 10. Although the current model appears to accurately mimic the physiological homeostasis of the lens [Table 6 & Table 7, there are still some aspects that need future improvement.
Our model presently has been implemented with a line suture structure which runs from the anterior to the posterior of the lens, through its core [Figure 3A&B]. However the mouse has a Y-shaped suture that rotates 180 °C from the anterior to the posterior pole . Hence, we are planning to improve the current model with an anatomically accurate asymmetrical 3D structure of the sutures. In our current model, the 3D extracellular solute diffusion coefficients are constant throughout. However, we [21, 53] and others [54–56] have recently identified a barrier to the movement of solutes in the lens, using variety of techniques. We strongly believe that the existence of this barrier is very important in shaping the fluid dynamics of the lens and in general its microcirculation. Hence, we will implement the imaged barrier as a part of continuous improvement of our computational model.
Another important structural feature of the ocular lens is the presence of a Gradient of Refractive Index (GRIN), which acts to correct for inherent spherical aberration to improve the optical properties of the lens . This GRIN profile of the lens is directly dependant on the local water/protein concentration makeup of the lens [58, 59]. Recently, we have experimentally showed that the water/protein gradient of the lens is actively upheld by the microcirculation system . Hence, another step in the improvement of the current model is to add equations to estimate the concentration gradient of water in the lens. Using that calculated field, we will be able to estimate the 3D GRIN maps of the lens. We would then use these 3D GRIN maps and optical ray-tracing software to produce a model that links lens physiology to the optical properties of the lens.
Our first generation 3D finite element model of lens structure and function describes ion and fluid dynamics in the mouse lens. We chose to model the mouse lens as ion and fluid dynamics have been extensively studied in this species [3, 4, 15, 16]. We also believe the model is an essential first step towards creating a comprehensive model of the human lens. Any model of the human lens model would need to include its more complex structural features and would need to be created so that its dimensions could be altered to study the effects of lens growth and ageing on the circulation system . This future model would enable us to study the changes in lens physiology thought to underlie the initiation of age related nuclear cataract.
During this project, a 3D model of the flux movements inside the mouse ocular lens was designed and executed using our high performance computer (HPC). Reviewing the results of the current model, it appears that solute fluxes, accompanied by water, enter the lens via the extracellular space all around it but with larger magnitudes around the polar regions. Among these inwardly extracellular solute fluxes, the Na+ fluxes were seemed to be dominant followed closely by Cl- fluxes. Conversely, the solute effluxes appear to be via the intracellular space and seemed to be more pronounced around the equatorial region of the lens. The K+ fluxes were found to be the primary intracellular fluxes, caused mainly by exterior Na+/K+ ATPase pumps. The net effect of these influx and effluxes were thought to be best explained by the calculated net current densities. The pattern of these net current densities at the surface of the lens was similar to previous experimental findings [5, 6, 52]. Same fields modelled inside the lens found to be in agreement with the microcirculation theory of the lens [3, 9, 16, 61].
This study brings together all the available experimental and theoretical data on the fluid dynamics of the ocular lens in order to create a comprehensive 3D model of this tissue. Previous studies have investigated the links between the steady state fluxes in the lens and its physiological homeostasis [62–65]. Using our computational model, we would be able to study the connections between the biodynamic natures of these perturbations and their functional consequences.
We wish to thank Prof Richard T Mathias for his insightful inputs in the development of the current model.
- Bassnett S: Lens organelle degradation. Exp Eye Res 2002, 74: 1–6. 10.1006/exer.2001.1111View ArticleGoogle Scholar
- Mathias RT: Steady-state voltages, ion fluxes, and volume regulation in syncytial tissues. Biophys J 1985, 48: 435–448. 10.1016/S0006-3495(85)83799-1View ArticleGoogle Scholar
- Mathias RT, Rae JL, Baldo GJ: Physiological properties of the normal lens. Physiol Rev 1997, 77: 21–50.Google Scholar
- Donaldson P, Kistler J, Mathias RT: Molecular solutions to mammalian lens transparency. News Physiol Sci 2001, 16: 118–123.Google Scholar
- Robinson KR, Patterson JW: Localization of steady currents in the lens. Curr Eye Res 1982, 2: 843–7. 10.3109/02713688209020020View ArticleGoogle Scholar
- Patterson JW: Characterization of the equatorial current of the lens. Ophthalmic Res 1988, 20: 139–42. 10.1159/000266570View ArticleGoogle Scholar
- Tamiya S, Dean WL, Paterson CA, Delamere NA: Regional distribution of Na, K-ATPase activity in porcine lens epithelium. Invest Ophthalmol Vis Sci 2003, 44: 4395–9. 10.1167/iovs.03-0287View ArticleGoogle Scholar
- Mathias RT, Rae JL, Eisenberg RS: Electrical properties of structural components of the crystalline lens. Biophys J 1979, 25: 181–201. 10.1016/S0006-3495(79)85284-4View ArticleGoogle Scholar
- Mathias RT, Rae JL: Steady state voltages in the frog lens. Curr Eye Res 1985, 4: 421–430. 10.3109/02713688509025156View ArticleGoogle Scholar
- Mathias RT, Rae JL, Ebihara L, McCarthy RT: The localization of transport properties in the frog lens. Biophys J 1985, 48: 423–434. 10.1016/S0006-3495(85)83798-XView ArticleGoogle Scholar
- Mathias RT, Rae JL, Eisenberg RS: The lens as a nonuniform spherical syncytium. Biophys J 1981, 34: 61–83. 10.1016/S0006-3495(81)84837-0View ArticleGoogle Scholar
- Rae JL, Kuszak JR: The electrical coupling of epithelium and fibers in the frog lens. Exp Eye Res 1983, 36: 317–326. 10.1016/0014-4835(83)90114-8View ArticleGoogle Scholar
- Zampighi GA: Distribution of connexin50 channels and hemichannels in lens fibers: a structural approach. Cell Commun Adhes 2003, 10: 265–270.View ArticleGoogle Scholar
- Candia OA, Zamudio AC: Regional distribution of the Na+ and K+ currents around the crystalline lens of rabbit. Am J Physiol Cell Physiol 2002, 282: C252.View ArticleGoogle Scholar
- Mathias RT, White TW, Gong X: Lens gap junctions in growth, differentiation, and homeostasis. Physiol Rev 2010, 90: 179–206. 10.1152/physrev.00034.2009View ArticleGoogle Scholar
- Mathias RT, Rae JL: Transport properties of the lens. Am J Physiol Cell Physiol 1985, 249: 181–190.Google Scholar
- Mathias RT, Wang H: Local osmosis and isotonic transport. J Membr Biol 2005, 208: 39–53. 10.1007/s00232-005-0817-9View ArticleGoogle Scholar
- Malcolm DTK: A Computational Model of the Ocular Lens [Internet]. ResearchSpace@ Auckland, Auckland, New Zealand; 2006. [cited 2012 Sep 11]. Available from: https://researchspace.auckland.ac.nz/handle/2292/754 Google Scholar
- Vaghefi E, Pontre BP, Jacobs MD, Donaldson PJ: Visualizing ocular lens fluid dynamics using MRI: manipulation of steady state water content and water fluxes. Am J Physiol Regul Integr Comp Physiol 2011, 301: R335-R342. 10.1152/ajpregu.00173.2011View ArticleGoogle Scholar
- Vaghefi E, Jacobs MD: Uptake and distribution of gadolinium in the ocular lens. Conf Proc IEEE Eng Med Biol Soc 2008, 2008: 843–846.Google Scholar
- Vaghefi E, Walker K, Pontre BP, Jacobs MD, Donaldson PJ: Magnetic resonance and confocal imaging of solute penetration into the lens reveals a zone of restricted extracellular space diffusion. Am J Physiol Regul Integr Comp Physiol 2012, 302(11):R1250–1259. 10.1152/ajpregu.00611.2011View ArticleGoogle Scholar
- Vaghefi E, Pontre B, Donaldson PJ, Hunter PJ, Jacobs MD: Visualization of transverse diffusion paths across fiber cells of the ocular lens by small animal MRI. Physiol Meas 2009, 30: 1061. 10.1088/0967-3334/30/10/007View ArticleGoogle Scholar
- Vaghefi E: Computational modeling and magnetic resonance imaging of microcirculation in the ocular lens. (Bioengineering)--University of Auckland, Auckland, New Zealand; 2010.Google Scholar
- Vaghefi E, Osman NAA, Abas WABW, Wahab AKA, Ting H-N: Ocular Lens Microcirculation Model, A Web-Based Bioengineering Educational Tool. In 5th Kuala Lumpur International Conference on Biomedical Engineering 2011. Springer Berlin Heidelberg, Volume 35. Berlin, Heidelberg; 2011:25–28.View ArticleGoogle Scholar
- Vaghefi E, Hunter PJ, Jacobs MD: 3D Finite Element Modeling of Avascular Circulation in the Ocular Lens. In 4th Kuala Lumpur International Conference on Biomedical Engineering 2008. Edited by: Abu Osman NA, Ibrahim F, Wan Abas WAB, Abdul Rahman HS, Ting H-N. Springer Berlin Heidelberg, Berlin, Heidelberg; 2008:469–472.View ArticleGoogle Scholar
- Duncan G: Permeability of amphibian lens membranes to water. Exp Eye Res 1970, 9: 188–197. 10.1016/S0014-4835(70)80075-6View ArticleGoogle Scholar
- Delamere NA, Duncan G: A comparison of ion concentrations, potentials and conductances of amphibian, bovine and cephalopod lenses. J Physiol 1977, 272: 167–186.View ArticleGoogle Scholar
- Guerschanik SN, Reinach PS, Candia OA: Chloride compartments of the frog lens and chloride permeabilities of its isolated surfaces. Invest Ophthalmol Vis Sci 1977, 16: 512–520.Google Scholar
- Gao J, Sun X, Moore LC, White TW, Brink PR, Mathias RT: Lens intracellular hydrostatic pressure is generated by the circulation of sodium and modulated by gap junction coupling. J Gen Physiol 2011, 137: 507–520. 10.1085/jgp.201010538View ArticleGoogle Scholar
- Benedek GB, Villars FMH 2nd ed. Springer. Physics with Illustrative Examples from Medicine and Biology: Statistical Physics 2000.Google Scholar
- Benedek GB, Villars FMH 2nd ed. Springer. Physics With Illustrative Examples from Medicine and Biology: Electricity and Magnetism 2000.Google Scholar
- Fischbarg J, Diecke FP, Kuang K, Yu B, Kang F, Iserovich P, Li Y, Rosskothen H, Koniarek JP: Transport of fluid by lens epithelium. Am J Physiol 1999, 276: C548–557.Google Scholar
- Candia OA, Gerometta R: Fluid movement across the surface of the isolated bovine lens. Invest Ophthalmol Vis Sci 2003, 44: 3455.Google Scholar
- Candia OA: Electrolyte and fluid transport across corneal, conjunctival and lens epithelia. Exp Eye Res 2004, 78: 527–535. 10.1016/j.exer.2003.08.015View ArticleGoogle Scholar
- Varadaraj K, Kumari S, Shiels A, Mathias RT: Regulation of aquaporin water permeability in the lens. Invest Ophthalmol Vis Sci 2005, 46: 1393–1402. 10.1167/iovs.04-1217View ArticleGoogle Scholar
- McLaughlin S, Mathias RT: Electro-osmosis and the reabsorption of fluid in renal proximal tubules. J Gen Physiol 1985, 85: 699–728. 10.1085/jgp.85.5.699View ArticleGoogle Scholar
- Benedek GB, Villars F: Physics, with Illustrative Examples from Medicine and Biology: Mechanics. 2000.View ArticleGoogle Scholar
- Cannell MB, Jacobs MD, Donaldson PJ, Soeller C: Probing microscopic diffusion by 2-photon flash photolysis: measurement of isotropic and anisotropic diffusion in lens fiber cells. Microsc Res Tech 2004, 63: 50–57. 10.1002/jemt.10422View ArticleGoogle Scholar
- Gao J, Sun X, Yatsula V, Wymore RS, Mathias RT: Isoform-specific function and distribution of Na/K pumps in the frog lens epithelium. J Membr Biol 2000, 178: 89–101. 10.1007/s002320010017View ArticleGoogle Scholar
- Jacobs MD, Soeller C, Sisley AMG, Cannell MB, Donaldson PJ: Gap junction processing and redistribution revealed by quantitative optical measurements of connexin46 epitopes in the lens. Invest Ophthalmol Vis Sci 2004, 45: 191–9. 10.1167/iovs.03-0148View ArticleGoogle Scholar
- Parmelee JT: Measurement of steady currents around the frog lens. Exp Eye Res 1986, 42: 433–41. 10.1016/0014-4835(86)90003-5View ArticleGoogle Scholar
- Kuszak JR, Mazurkiewicz M, Zoltoski R: Computer modeling of secondary fiber development and growth: I. Nonprimate lenses. Molecular Vision 2006, 12: 271–282.Google Scholar
- Fischer G, Tilg B, Modre R, Huiskamp G, Fetzer J, Rucker W, Wach P: A bidomain model based BEM-FEM coupling formulation for anisotropic cardiac tissue. Ann Biomed Eng 2000, 28: 1229–1243.View ArticleGoogle Scholar
- Duncan G, Jacob TJ: Influence of external calcium and glucose on internal total and ionized calcium in the rat lens. J Physiol 1984, 357: 485–493.View ArticleGoogle Scholar
- Zampighi GA, Eskandari S, Kreman M: Epithelial organization of the mammalian lens. Exp Eye Res 2000, 71: 415–435. 10.1006/exer.2000.0895View ArticleGoogle Scholar
- Blankenship T, Bradshaw L, Shibata B, FitzGerald P: Structural specializations emerging late in mouse lens fiber cell differentiation. Invest Ophthalmol Vis Sci 2007, 48: 3269. 10.1167/iovs.07-0109View ArticleGoogle Scholar
- Duncan G: Movement of sodium and chloride across amphibian lens membranes. Exp Eye Res 1970, 10: 117. 10.1016/S0014-4835(70)80017-3View ArticleGoogle Scholar
- Hunter PJ, Borg TK: Integration from proteins to organs: the Physiome Project. Nat Rev Mol Cell Biol 2003, 4: 237–243. 10.1038/nrm1054View ArticleGoogle Scholar
- Rubin J, Wünsche BC, Cameron L, Stevens C: Animation and modelling of cardiac performance for patient monitoring. In Proceedings of IVCNZ 2005, 5: 476–481.Google Scholar
- Wang H, Gao J, Sun X, Martinez-Wittinghan FJ, Li L, Varadaraj K, Farrell M, Reddy VN, White TW, Mathias RT: The effects of GPX-1 knockout on membrane transport and intracellular homeostasis in the lens. J Membr Biol 2009, 227: 25–37. 10.1007/s00232-008-9141-5View ArticleGoogle Scholar
- Delamere NA, Paterson CA: The influence of calcium-free solutions upon permeability characteristics of the rabbit lens. Exp Eye Res 1979, 28: 45–53. 10.1016/0014-4835(79)90104-0View ArticleGoogle Scholar
- Parmelee JT, Robinson KR, Patterson JW: Effects of calcium on the steady outward currents at the equator of the rat lens. Invest Ophthalmol Vis Sci 1985, 26: 1343–8.Google Scholar
- Grey AC, Jacobs MD, Gonen T, Kistler J, Donaldson PJ: Insertion of MP20 into lens fibre cell plasma membranes correlates with the formation of an extracellular diffusion barrier. Exp Eye Res 2003, 77: 567–574. 10.1016/S0014-4835(03)00192-1View ArticleGoogle Scholar
- Moffat BA, Pope JM: Anisotropic water transport in the human eye lens studied by diffusion tensor NMR micro-imaging. Exp Eye Res 2002, 74: 677–87. 10.1006/exer.2001.1164View ArticleGoogle Scholar
- Truscott RJW: Age-related nuclear cataract-oxidation is the key. Exp Eye Res 2005, 80: 709–25. 10.1016/j.exer.2004.12.007View ArticleGoogle Scholar
- Moffat BA, Landman KA, Truscott RJ, Sweeney MH, Pope JM: Age-related changes in the kinetics of water transport in normal human lenses. Exp Eye Res 1999, 69: 663–9. 10.1006/exer.1999.0747View ArticleGoogle Scholar
- Smith G, Cox MJ, Calver R, Garner LF: The spherical aberration of the crystalline lens of the human eye. Vision Res 2001, 41: 235–243. 10.1016/S0042-6989(00)00206-6View ArticleGoogle Scholar
- Pierscionek BK, Chan DY: Refractive index gradient of human lenses. Optom Vis Sci 1989, 66: 822–829. 10.1097/00006324-198912000-00004View ArticleGoogle Scholar
- Ioseliani OR, PIERSCIONEK B illustrated. Nova Publishers. Species variations in the refractive index of the eye lens. Focus on Eye Research 2006, 238.Google Scholar
- Dubbelman M, Van der Heijde G: The shape of the aging human lens: curvature, equivalent refractive index and the lens paradox. Vision Res 2001, 41: 1867–1877. 10.1016/S0042-6989(01)00057-8View ArticleGoogle Scholar
- Mathias RT, Kistler J, Donaldson P: The lens circulation. J Membr Biol 2007, 216: 1–16. 10.1007/s00232-007-9019-yView ArticleGoogle Scholar
- Kobatashi S, Roy D, Spector A: Sodium/potassium ATPase in normal and cataractous human lenses. Curr Eye Res 1982, 2: 327–334. 10.3109/02713688209000778View ArticleGoogle Scholar
- Spector A: Oxidative stress-induced cataract: mechanism of action. FASEB J 1995, 9: 1173.Google Scholar
- White TW, Goodenough DA, Paul DL: Targeted ablation of connexin50 in mice results in microphthalmia and zonular pulverulent cataracts. J Cell Biol 1998, 143: 815–25. 10.1083/jcb.143.3.815View ArticleGoogle Scholar
- Buzhynskyy N, Girmens JF, Faigle W, Scheuring S: Human cataract lens membrane at subnanometer resolution. J Mol Biol 2007, 374: 162–169. 10.1016/j.jmb.2007.09.022View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.