A method for efficient calculation of diffusion and reactions of lipophilic compounds in complex cell geometry

A general description of effects of toxic compounds in mammalian cells is facing several problems. Firstly, most toxic compounds are hydrophobic and partition phenomena strongly influence their behaviour. Secondly, cells display considerable heterogeneity regarding the presence, activity and distribution of enzymes participating in the metabolism of foreign compounds i.e. bioactivation/biotransformation. Thirdly, cellular architecture varies greatly. Taken together, complexity at several levels has to be addressed to arrive at efficient in silico modelling based on physicochemical properties, metabolic preferences and cell characteristics. In order to understand the cellular behaviour of toxic foreign compounds we have developed a mathematical model that addresses these issues. In order to make the system numerically treatable, methods motivated by homogenization techniques have been applied. These tools reduce the complexity of mathematical models of cell dynamics considerably thus allowing to solve efficiently the partial differential equations in the model numerically on a personal computer. Compared to a compartment model with well-stirred compartments, our model affords a more realistic representation. Numerical results concerning metabolism and chemical solvolysis of a polycyclic aromatic hydrocarbon carcinogen show good agreement with results from measurements in V79 cell culture. The model can easily be extended and refined to include more reactants, and/or more complex reaction chains, enzyme distribution etc, and is therefore suitable for modelling cellular metabolism involving membrane partitioning also at higher levels of complexity.