Improving the Prediction of Local Drug Distribution Profiles in the Brain with a New 2D Mathematical Model

The development of drugs that target the brain is very challenging. A quantitative understanding is needed of the complex processes that govern the concentration–time profile of a drug (pharmacokinetics) within the brain. So far, there are no studies on predicting the drug concentration within the brain that focus not only on the transport of drugs to the brain through the blood–brain barrier (BBB), but also on drug transport and binding within the brain. Here, we develop a new model for a 2D square brain tissue unit, consisting of brain extracellular fluid (ECF) that is surrounded by the brain capillaries. We describe the change in free drug concentration within the brain ECF, by a partial differential equation (PDE). To include drug binding, we couple this PDE to two ordinary differential equations that describe the concentration–time profile of drug bound to specific as well as non-specific binding sites that we assume to be evenly distributed over the brain ECF. The model boundary conditions reflect how free drug enters and leaves the brain ECF by passing the BBB, located at the level of the brain capillaries. We study the influence of parameter values for BBB permeability, brain ECF bulk flow, drug diffusion through the brain ECF and drug binding kinetics, on the concentration–time profiles of free and bound drug.