# Diagnostic biomarker kinetics: how brain-derived biomarkers distribute through the human body, and how this affects their diagnostic significance: the case of S100B

Robert Murcko, Nicola Marchi, Damian Bailey, Damir Janigro

► To cite this version:

Robert Murcko, Nicola Marchi, Damian Bailey, Damir Janigro. Diagnostic biomarker kinetics: how brain-derived biomarkers distribute through the human body, and how this affects their diagnostic significance: the case of S100B. Fluids and Barriers of the CNS, 2022, 19 (1), pp.32. 10.1186/s12987-022-00329-9. http://hal-03672124

HAL Id: hal-03672124

https://hal.science/hal-03672124

Submitted on 27 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

Check for
updates

# Diagnostic biomarker kinetics: how brain-derived biomarkers distribute through the human body, and how this affects their diagnostic significance: the case of S100B

Robert Murcko<fnref n="1" />, Nicola Marchi³, Damian Bailey⁴ and Damir Janigro<fnref n="1" />,²*

## Abstract

Blood biomarkers of neurological diseases are often employed to rule out or confirm the presence of significant intracranial or cerebrovascular pathology or for the differential diagnosis of conditions with similar presentations (e.g., hemorrhagic vs. embolic stroke). More widespread utilization of biomarkers related to brain health is hampered by our incomplete understanding of the kinetic properties, release patterns, and excretion of molecules derived from the brain. This is, in particular, true for S100B, an astrocyte-derived protein released across the blood-brain barrier (BBB). We developed an open-source pharmacokinetic computer model that allows investigations of biomarker's movement across the body, the sources of biomarker's release, and its elimination. This model was derived from a general in silico model of drug pharmacokinetics adapted for protein biomarkers. We improved the model's predictive value by adding realistic blood flow values, organ levels of S100B, lymphatic and glymphatic circulation, and glomerular filtration for excretion in urine. Three key variables control biomarker levels in blood or saliva: blood-brain barrier permeability, the S100B partition into peripheral organs, and the cellular levels of S100B in astrocytes. A small contribution to steady-state levels of glymphatic drainage was also observed; this mechanism also contributed to the uptake of organs of circulating S100B. This open-source model can also mimic the kinetic behavior of other markers, such as GFAP or NF-L. Our results show that S100B, after uptake by various organs from the systemic circulation, can be released back into systemic fluids at levels that do not significantly affect the clinical significance of venous blood or salivary levels after an episode of BBB disruption.

**Keywords:** Computer model, MATLAB, Simbiology, Astrocytes, Physiologically-based pharmacokinetic model, Glymphatics, Extracranial sources, Traumatic brain injury, Brain barriers, Saliva

## Background

Pharmacokinetic in silico models are essential for pharmacological studies and drug development. During the drug discovery and development process, potential clinical candidates are screened for their absorption, distribution, metabolism, and excretion (ADME) properties

to avoid clinic failures related to inappropriate ADME properties. Until recently, most pharmacokinetic models were aimed at predicting the properties of small (< 1 kD) molecules after oral or intravenous (i.v.) administration. Recently, biologics have become a significant portion of therapeutic agents, and the old small molecule software strategies had to be reformulated to adapt to large (> 10kD) molecular weight proteic therapeutics. While ad hoc software has been developed by Industry, academic

<footnote>*Correspondence: djanigro@flocel.com</footnote>

<footnote n="1">¹ FloTBI Inc., Cleveland, OH, USA</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
efforts have used available platforms (e.g., MATLAB) to model how drugs distribute in the body.

A perhaps unexpected utilization of pharmacokinetic modeling of large, proteic agents is the development of modified strategies to study the movement of diagnostic molecules in the human body. Several of these biomarkers are proteins with varying molecular properties and sizes. These are, most commonly, not administered conventionally but are instead released or synthesized ex novo by a specific organ, neoplasm, or cell type. For example, troponins are proteins found in skeletal and cardiac muscle fibers that regulate muscular contraction. Troponin tests measure the level of cardiac-specific troponin in the blood to help detect heart injury [1]. When there is damage to heart muscle cells, troponin is released into the blood, thus becoming detectable by a simple blood test. The necessity of cellular death for biomarker release is not universal since many other biomarkers are released by healthy cells (see below for S100B). Our previous work has focused primarily on brain-derived diagnostic molecules used to diagnose CNS or neurological diseases. These include GFAP, S100B, UCHL-1, and other less-studied reporters of brain disease or health [2]. An example of how pharmacokinetic models can be applied to brain diagnostic markers was published [3].

We adapted and refined MATLAB-based models [3–5] for the present study using the published data obtained by real-life experiments (direct measurements of S100B from human tissues; see [3–5]). We specifically wished to explore the pharmacokinetics of S100B, a reporter of blood–brain barrier (BBB) dysfunction (BBBD) and brain health [6–9]. While several studies have promoted its use in neurology and psychiatry [3, 7, 9, 10], others expressed doubts about its reliability for human diagnostics. These concerns primarily derive from the pitfalls listed below.

It was suggested that S100B not only derives from the brain but also has extracranial sources [11–16]. Thus, when both brain and peripheral trauma are involved, it is impossible to dissect a central vs. peripheral origin of the biomarker. This is a pitfall in studies where S100B is elevated in individuals with multi-trauma of orthopedic nature [17]. The same issues were reported for other brain-derived biomarkers [18]. Several counterarguments have been made, showing, for example, that extracranial sources where S100B is synthesized from mRNA are few (e.g., testis, descending tubules in kidney [3]) and that S100B content in other organs derives from uptake from circulation [4]. It was recently shown that time-dependent internalization of circulating S100B by mesenchymal stem cells occurs via the pathways of clathrin- and lipid raft-mediated endocytosis [19]. Others

have demonstrated that S100B in fat tissue does not contribute to peripherally detected levels [20], but the opposite was also suggested [21]. Therefore, controversy exists on the extent and relevance of extracranial sources of biomarkers used for CNS diagnostics.

In the field of sports medicine, it was shown that blood S100B increases after sub-concussive head hits [22, 23]: this was ascribed to increased BBB permeability as also documented by MRI [24]. Other studies have shown that S100B is increased by exercise alone [12, 25, 26], while others found no effect of strenuous exercise on S100B levels [27–30]. An explanation of these contrasting findings points to BBB damage induced by extreme exercise [2]. According to this hypothesis, strenuous exercise or performance in extreme sports results in BBB “opening,” possibly due to a mechanism involving free radical formation, as suggested by ref. [31]. In any case, it is not known how different sources of S100B contribute to the peripheral signal in blood (or saliva) [2, 26, 32, 33]. Lastly, a common motif in S100B diagnostics is that S100B is not specific for any neurological disease [3]. This is due to the fact that BBB leakage allowing S100B appearance in peripheral body fluids is a common feature of many neurological conditions [3].

Another point of contention relates to how the brain releases S100B during an insult. The leading hypothesis calls for a breach of the BBB as described above and in [34, 35]. An alternative hypothesis calls for the recently described glymphatic circulation as a means of brain release of biomarkers in blood [36]. The contribution of glymphatics in human subjects is unknown.

The scope of the present work is to answer, when possible, these questions by using two advanced full-body models of cerebrovascular and peripheral circulation after the release of S100B by the brain or other organs. A lymphatic compartment was also added to the model together with realistic urinary excretion pathways. The initial parameters of the model were derived from experimental observations and available human data [3].

## Methods

We used MATLAB 2019–2021b (MathWorks, Natick, MA) to design, test and simulate the model. The toolbox used was MATLAB’s Simbiology app (versions 5.0–6.2), aided by packages for partial differential equations, statistics, and parallel computing. Data were plotted in MATLAB and exported to CorelDraw (Corel Co.) as extended metafiles.

We developed two separate and independent models to mimic the behavior of circulating brain-derived small molecular weight proteic biomarkers. Model 1 was used
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Schematic representation of the structure of Model 1. The continuous red lines depict flow through arteries (except for the lung), while the continuous blue lines refer to venous flow. The dotted blue lines show the lymphatic vessels connecting the organs directly to the venous compartment. The right panel inset shows a graphical representation of the mechanism of BBBD presented herein, underscoring that venous levels are greatly influenced by leakage of biomarkers from the brain into the circulation</figcaption></figure>

primarily to assess the relevance of various peripheral organs to the signal measured in blood (Fig. 1). Model 2 was developed after ruling out the contributions by heart, bone, and skin (Figs. 4 and 5): these organs are not included in Model 2, which uses a different set of equations to focus on the contributions of adipose, muscle, and gut tissues to measured levels of the biomarker in blood. Model 2 also adopts a more complex brain modeling, as detailed below. The following sections highlight the shared and specific modeling strategies used.

Most of the simulations shown were run to steady-state with or without an accompanying BBB disruption event. This allows to follow up the kinetics of S100B in each organ or compartments. For Fig. 10, the BBBD was triggered after steady-state was achieved.

## Model 1

A human full-body physiologically-based pharmacokinetic (PBPK) model was adapted from [37]. This model contains lung, brain, skin, bone, adipose tissue, heart, kidney, muscle, and gut (Fig. 1). The volume of these organs is specified in Table 1. The organs are connected
<table><caption>**Table 1** Parameter values used for Model 1 (Fig. 2)</caption><thead><tr><th>Quantity name</th><th>Initial value</th><th>Units</th></tr></thead><tbody><tr><td>Arterial blood</td><td>2500.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Brain</td><td>1340.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>9.00</td><td>Nanogram/milliliter</td></tr><tr><td>Kidney</td><td>325.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.25</td><td>Nanogram/milliliter</td></tr><tr><td>Urine</td><td>600.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Venous blood</td><td>2500.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Lung</td><td>547.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.37</td><td>Nanogram/milliliter</td></tr><tr><td>Heart</td><td>359.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.28</td><td>Nanogram/milliliter</td></tr><tr><td>Bone</td><td>3950.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Adipose tissue</td><td>22,700.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.51</td><td>Nanogram/milliliter</td></tr><tr><td>Skin</td><td>3150.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.28</td><td>Nanogram/milliliter</td></tr><tr><td>Muscle</td><td>31,300.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.30</td><td>Nanogram/milliliter</td></tr><tr><td>Gut</td><td>1220.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.40</td><td>Nanogram/milliliter</td></tr><tr><td>Lymph</td><td>12.00</td><td>Liter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Muscle tissue vascular fraction</td><td>845.10</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Gut tissue vascular fraction</td><td>61.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Lung tissue vascular fraction</td><td>101.20</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Heart tissue vascular fraction</td><td>15.08</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.28</td><td>Nanogram/milliliter</td></tr><tr><td>Bone tissue vascular fraction</td><td>197.50</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Kidney tissue vascular fraction</td><td>22.75</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Skin tissue vascular fraction</td><td>157.50</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Adipose tissue vascular fraction</td><td>703.70</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Brain tissue vascular fraction</td><td>67.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>BloodFlowLungToArtery</td><td>313,980.00</td><td>milliliter/hour</td></tr><tr><td>BloodFlowVenousToLung</td><td>313,980.00</td><td>milliliter/hour</td></tr><tr><td>BloodFlowArteryToKidney</td><td>66,000.00</td><td>milliliter/hour</td></tr><tr><td>BBB_Index</td><td>0.00</td><td>Dimensionless</td></tr><tr><td>GFR</td><td>125.00</td><td>Milliliter/minute</td></tr><tr><td>BloodFlowArteryToSaliva</td><td>0.00</td><td>1/hour</td></tr></tbody></table>

<table><caption>**Table 1** (continued)</caption><thead><tr><th>Quantity name</th><th>Initial value</th><th>Units</th></tr></thead><tbody><tr><td>BloodFlowArteryToVenous</td><td>3139.00</td><td>Milliliter/hour</td></tr><tr><td>BloodFlowArteryToMuscle</td><td>45,000.00</td><td>Milliliter/hour</td></tr><tr><td>BloodFlowArteryToHeart</td><td>9000.00</td><td>Milliliter/hour</td></tr><tr><td>BloodFlowHeartToVenous</td><td>9000.00</td><td>Milliliter/hour</td></tr><tr><td>BloodFlowArterialToBone</td><td>15,000.00</td><td>Milliliter/hour</td></tr><tr><td>BloodFlowBoneToVenous</td><td>15,000.00</td><td>Milliliter/hour</td></tr><tr><td>BloodFlowKidneyToVenous</td><td>66,000.00</td><td>Milliliter/hour</td></tr><tr><td>BloodFlowAdiposeToVenous</td><td>15,600.00</td><td>Milliliter/hour</td></tr><tr><td>BloodFlowArterialToAdipose</td><td>15,600.00</td><td>Milliliter/hour</td></tr><tr><td>BloodFlowSkinToVenous</td><td>18,000.00</td><td>Milliliter/hour</td></tr><tr><td>BloodFlowArterialToSkin</td><td>18,000.00</td><td>Milliliter/hour</td></tr><tr><td>TissueFactor</td><td>0.18</td><td>Dimensionless</td></tr><tr><td>BloodFlowMuscleToVenous</td><td>45,000.00</td><td>Milliliter/hour</td></tr><tr><td>BloodFlowArteryToGut</td><td>66,000.00</td><td>Milliliter/hour</td></tr><tr><td>BloodFlowGutToVenous</td><td>66,000.00</td><td>Milliliter/hour</td></tr><tr><td>LymphFlowMuscle</td><td>90.00</td><td>Milliliter/hour</td></tr><tr><td>LymphFlowGut</td><td>132.00</td><td>Milliliter/hour</td></tr><tr><td>LymphFlowHeart</td><td>132.00</td><td>Milliliter/hour</td></tr><tr><td>LymphFlowAdipose</td><td>31.20</td><td>Milliliter/hour</td></tr><tr><td>LymphFlowSkin</td><td>36.00</td><td>Milliliter/hour</td></tr><tr><td>LymphFlowKidney</td><td>132.00</td><td>Milliliter/hour</td></tr><tr><td>LymphFlowToVein</td><td>3100.00</td><td>Milliliter/hour</td></tr><tr><td>kf_1</td><td>1.00</td><td>Dimensionless</td></tr><tr><td>BloodFlowLungToArterial</td><td>313,980.00</td><td>Milliliter/hour</td></tr><tr><td>Excretion</td><td>0.00</td><td>1/minute</td></tr><tr><td>RenalEliminationFactor</td><td>100.00</td><td>Dimensionless</td></tr><tr><td>kf_brain</td><td>0.02</td><td>Dimensionless</td></tr><tr><td>AdiposeTissueFactor</td><td>1.00</td><td>Dimensionless</td></tr><tr><td>PoreRatioMuscle</td><td>2000.00</td><td>Dimensionless</td></tr><tr><td>PoreRatioAdipose</td><td>500.00</td><td>Dimensionless</td></tr><tr><td>PoreRatioGut</td><td>500.00</td><td>Dimensionless</td></tr><tr><td>PoreRatioSkin</td><td>500.00</td><td>Dimensionless</td></tr><tr><td>PoreRatioLung</td><td>45.00</td><td>Dimensionless</td></tr><tr><td>PoreRatioKidney</td><td>200.00</td><td>Dimensionless</td></tr><tr><td>kr_brain</td><td>0.00</td><td>Dimensionless</td></tr><tr><td>kf_muscle</td><td>24.44</td><td>Dimensionless</td></tr><tr><td>kr_muscle</td><td>0.00</td><td>Dimensionless</td></tr><tr><td>kf_gut</td><td>104.17</td><td>Dimensionless</td></tr><tr><td>kr_gut</td><td>0.01</td><td>Dimensionless</td></tr><tr><td>kf_lung</td><td>617.28</td><td>Dimensionless</td></tr><tr><td>kr_lung</td><td>0.06</td><td>Dimensionless</td></tr><tr><td>kf_heart</td><td>130.21</td><td>Dimensionless</td></tr><tr><td>kr_heart</td><td>0.01</td><td>Dimensionless</td></tr><tr><td>kf_bone</td><td>797.10</td><td>Dimensionless</td></tr><tr><td>kr_bone</td><td>0.08</td><td>Dimensionless</td></tr><tr><td>kf_kidney</td><td>135.14</td><td>Dimensionless</td></tr><tr><td>kr_kidney</td><td>0.01</td><td>Dimensionless</td></tr><tr><td>kf_skin</td><td>66.67</td><td>Dimensionless</td></tr><tr><td>kr_skin</td><td>0.01</td><td>Dimensionless</td></tr></tbody></table>
To create a more pronounced differential within each tissue but keep the ratio of k values between each tissue standardized, the kf value was the kr value multiplied by a factor of 10,000, thus kf_organ = kr_organ * 10,000. This value was empirically chosen to match the rank order results for levels of S100B measured in various organs [3]. The order of measured values was brain > adipose > kidney > heart > muscle > lung > gut. The multiplier was derived by running the model with appropriate values to match the rank order of measured values. These values were then used as initial conditions.

The general equation for the organ's uptake or release of protein biomarkers was:

$$ \frac{d(\text{Organ.S100B})}{dt} = \frac{1}{\text{Organ}} * \left( - \left( \left( \begin{array}{l} \text{kf\_Organ} * \text{Organ.S100B} - \text{kr\_Organ} \\ * \{ \text{Organ Tissue Vascular Fraction} \} * \text{S100B} \end{array} \right) * \text{interstitial Flow} \right) \right) \quad (2) $$

where the value *Interstitial flow* represents the flow rate of the protein within the organ, *Organ* and *Tissue Vascular Fraction* volumes were derived from ref. [37]. Organ. S100B refers to the concentration of S100B within the specified organ. In previous and subsequent equations all of these variables (Organ, Organ.S100B, Organ Tissue Vascular Fraction, etc.) are labeled with tags to the specific organ that they are referring to.

Excretion of the biomarker protein was modeled by kidney filtration:

$$ \frac{d(\text{Organ.S100B})}{dt} = \frac{1}{\text{Urine}} * \left( \left( \text{Renal Elimination Factor} * \frac{\text{GFR}}{\text{Kidney}} * \text{Kidney.S100B} \right) * \text{Kidney} \right) \quad (3) $$

where *GFR* is the glomerular filtration rate, and the *Renal elimination factor* is an additional dimensionless tuning parameter ranging from 0 to 1.

Sensitivity analysis is the study of how the uncertainty in the output of a mathematical model can be divided and allocated to different sources of uncertainty in its inputs. In Simbiology, the routine of sensitivity analysis allows determining which rate constants and concentrations in a model significantly influence the overall behavior of the model (https://www.mathworks.com/help/simbio/ug/global-local-sensitivity-analysis-gsa-lsa-simbiology.html). SimBiology supports two types of sensitivity analyses: local and global sensitivity analysis (GSA). GSA uses Monte Carlo simulations, where a representative (global) set of parameter sample values are used to explore the effects of variations in model parameters of interest on the model response. In this approach, SimBiology performs a decomposition of the model output

(response) variance by calculating the first- and total-order Sobol indices. The first-order Sobol indices give the fractions of the overall response variance that can be attributed to variations in an input parameter alone. The total-order Sobol index gives the fraction of the overall response variance that can be attributed to joint parameter variations (see [40]). We used global sensitivity analysis to interpret the impact of S100B in various organs on venous biomarker levels (Figs. 3 and 6; See also Additional file 1: Fig. S1, Additional file 2: Fig. S2, Additional file 3: Fig. S3). In Model 1, sensitivity analyses were run with BBB_Index set to either 0 or 1 (Fig. 3C and D, respectively). For Model 2, we explored the effect

of changing S100B in the brain interstitium (1 and 10 ng/ml). The data in Supplemental figures were obtained by a Sobol sampling interpolation method, with 1000 samples; the simulation was run to steady state. The data are shown as time course (Additional file 1: Fig. S1, Additional file 2: Fig. S2) or bar graphs (Additional file 3: Fig. S3).

## Model 2

Model 2 follows the general structural backbone of Model 1 (Fig. 4). However, organs (except for the kidney, see below) are subdivided into vascular and interstitial compartments (Table 2). To describe the passage of protein from the interstitial (parenchymal) space into vascular space, we used the coefficient of vascular reflection (Sigma, or σ) as per reference [41]. The size-dependent restriction of large pores and small pores can, in fact, be represented as the vascular reflection coefficient, an indirect measure of the density of exchange pores. The model used to mimic the brain (Fig. 5) used BBB_Index and Trauma_Index to describe the passage of S100B across the interstitial, vascular, and cellular compartments. Note that unlike the dimensionless rate constants in Model 1, kinetic variables have dimensions of quantity/time in this model. The equation governing changes of biomarker's levels in the vascular compartment was:
<figure><img src="image_4.png" /><figcaption>**Fig. 3** Results from simulations and sensitivity analysis (Model 1). The left panel shows the actual levels of S100B in organs before (A: BBBD = 0) and after BBB disruption (B) obtained by setting the BBB\_Index value to 1. Note the pronounced increase of S100B in venous blood with comparably smaller changes in other compartments. Figures C and D show the results of sensitivity analysis queries under the same conditions. Note that before BBBD (A) venous levels at steady state never reached the 0.1 ng/ml thresholds, the upper ceiling for control values of S100B. The x axis in the bar plots specifies the sensitivity inputs and y axes the sensitivity outputs integrated over time</figcaption></figure>

$$ \frac{d\{\text{Vascular Compartment Brain}\}.S100B}{dt} = \frac{1}{\text{Vascular Compartment Brain}} \quad (4) \\ * ((\text{Arterial To Brain Blood Flow}) * \{\text{Arterial Blood}\}.S100B) \\ + ((\text{BBB\_Index} * \{\text{Interstitium Brain}\}.S100B) * \{\text{Interstitium Brain}\}) $$
where BBB_Index can change between 0 and 1 to mimic increased permeability of the cerebral vasculature. In addition to having a three compartment structure, brain modeling also included glymphatic drainage into central lymph and venous blood. The equation for brain interstitium S100B was thus:

$$ \frac{d(\text{Interstitium Brain}.S100B)}{\frac{dt}{1}} = \frac{\text{Interstitium Brain}}{\text{Interstitium Brain} \cdot \text{Glia}} \cdot \left( \begin{array}{l} ((\text{Trauma Index} \cdot \text{Glia.S100B}) \cdot \text{Glia}) \\ - (\text{Glymphatics} \cdot \text{Interstitium Brain}.S100B) \\ - ((\text{BBB Index} \cdot \text{Interstitium Brain}.S100B) \cdot \text{Interstitium Brain}) \end{array} \right) \quad (5) $$

where the term Trauma_Index refers to the passage of S100B from astrocytes in the cellular compartment (Glia) released directly into the brain interstitium. Glymphaticds is the rate of interstitial flow to Central lymph. Changes of S100B in the cellular compartment were described by:

$$ \text{Interstitium Brain}.S100B = \text{Glia.S100B} \cdot \text{TraumaIndex} \quad (6) $$

Glia.S100B was set constant at 50 ng/ml as per in vitro measurements by others [42].

The kidney was modeled by a single compartment with input from arterial blood and an output to urine. The process was described by:

<figure><img src="image_5.png" /><figcaption>Fig. 4 Schematic representation of Model 2. Note the addition of a salivary and glymphatic component. Peripheral organs are subdivided into interstitial and vascular compartments, while the brain is represented by three compartments, see Fig. 5. At the beginning of our simulations, S100B values in venous and lymphatic compartments were set to 0</figcaption></figure>

## Results

The structure of Model 1 is shown in Fig. 1, together with the graphic rendition of the process of BBB disruption. Figure 2 shows the formalism used to describe each organ in Model 1. The simulation of Model 1 led to the results shown in Fig. 3, which represents the changes in organs’

$$ \frac{d(\text{Kidney.S100B})}{dt} = \frac{1}{\text{Kidney}} \cdot \left( \begin{array}{l} - \left( \left( \text{Elimination Factor} \cdot \frac{\text{GFR}}{\text{Kidney}} \cdot \frac{\text{Kidney.S100B}}{\text{Tissue Partition Kidney}} \right) \cdot \text{Kidney} \right) \\ + (\text{Arterial To Kidney Blood Flow} \cdot \{\text{Arterial Blood}\} \cdot \text{S100B}) \\ - (\text{Kidney To Venous Blood Flow} \cdot \text{Kidney.S100B}) \end{array} \right) \quad (7) $$

and by:

$$ \frac{d(\text{Urine.S100B})}{dt} = \frac{1}{\text{Urine}} \cdot \left( \left( \left( \text{Elimination Factor} \cdot \frac{\text{GFR}}{\text{Kidney}} \cdot \frac{\text{Kidney.S100B}}{\text{Tissue Partition Kidney}} \right) \cdot \text{Kidney} \right) - (\text{kf} \cdot \text{Urine.S100B}) \right) \quad (8) $$

Tissue partition kidney was set at 0.1/min. GFR was set at 10 ml/h; this non-constant value was explored during simulations (e.g., Fig. 8).

Note that in all figures, except Fig. 10, the simulation started before steady-state conditions were reached, thus allowing the variables to express the kinetic significance of the underlying code. See, for example, Fig. 3A, B, where the time-dependent changes in S100B are shown.

S100B levels under normal conditions (A, BBBD = 0) or after BBB disruption (B; BBB_Index = 1). Without BBB disruption, individual organs displayed a change in parenchymal S100B content to eventually reach steady state. Also, note that venous levels, at steady state, were <0.1 ng/ml, which is consistent with clinical studies in normal adults when using the Roche Diagnostics test [43, 44]. Panels C and D show the results of the simulation in
<table><caption>**Table 2** Parameter values used for Model 2 (Fig. 4)</caption><thead><tr><th>Quantity name</th><th>Initial value</th><th>Units</th></tr></thead><tbody><tr><td>Arterial blood</td><td>1.16</td><td>Liter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Venous blood</td><td>2.33</td><td>Liter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Central lymph</td><td>12.00</td><td>Liter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Brain</td><td>1.34</td><td>Liter</td></tr><tr><td>Vascular compartment brain</td><td>67.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Interstitium brain</td><td>130.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>10.00</td><td>Nanogram/milliliter</td></tr><tr><td>Glia</td><td>0.20</td><td>Liter</td></tr><tr><td>S100B</td><td>50.00</td><td>Nanogram/milliliter</td></tr><tr><td>Urine</td><td>200.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Kidney</td><td>280.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.25</td><td>Nanogram/milliliter</td></tr><tr><td>Muscle</td><td>31.30</td><td>Liter</td></tr><tr><td>Interstitium muscle</td><td>2.80</td><td>Liter</td></tr><tr><td>S100B</td><td>0.30</td><td>Nanogram/milliliter</td></tr><tr><td>Vascular space muscle</td><td>0.85</td><td>Liter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Adipose</td><td>22.70</td><td>Liter</td></tr><tr><td>Interstitium adipose</td><td>3.20</td><td>Liter</td></tr><tr><td>S100B</td><td>0.50</td><td>Nanogram/milliliter</td></tr><tr><td>Vascular space adipose</td><td>703.10</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Gut</td><td>1220.00</td><td>Milliliter</td></tr><tr><td>Vascular space gut</td><td>61.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>Interstitium gut</td><td>325.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.18</td><td>Nanogram/milliliter</td></tr><tr><td>Saliva</td><td>10.00</td><td>Milliliter</td></tr><tr><td>S100B</td><td>0.00</td><td>Nanogram/milliliter</td></tr><tr><td>ArterialToBrainBloodFlow</td><td>42,000.00</td><td>Milliliter/hour</td></tr><tr><td>BrainToVenousBloodFlow</td><td>42,000.00</td><td>Milliliter/hour</td></tr><tr><td>BBB_index</td><td>0.00</td><td>1/hour</td></tr><tr><td>Trauma_index</td><td>0.00</td><td>1/hour</td></tr><tr><td>Glymphatics</td><td>8.40</td><td>Milliliter/hour</td></tr><tr><td>VenousToArterialBloodFlow</td><td>201,874.29</td><td>Milliliter/hour</td></tr><tr><td>LymphaticToVenous</td><td>100.00</td><td>Milliliter/hour</td></tr><tr><td>GFR</td><td>10.00</td><td>Milliliter/hour</td></tr><tr><td>ArterialToKidneyBloodFlow</td><td>66,000.00</td><td>Milliliter/hour</td></tr><tr><td>EliminationFactor</td><td>100.00</td><td>1/minute</td></tr><tr><td>KidneyToVenousBloodFlow</td><td>66,000.00</td><td>Milliliter/hour</td></tr><tr><td>ReverseSigmaGut</td><td>0.05</td><td>1/hour</td></tr><tr><td>ReverseSigmaAdipose</td><td>0.05</td><td>1/hour</td></tr><tr><td>ReverseSigmaMuscle</td><td>0.05</td><td>1/hour</td></tr><tr><td>BBB</td><td>100.00</td><td>1/minute</td></tr><tr><td>SalivaryBloodFlow</td><td>1.00</td><td>Milliliter/minute</td></tr></tbody></table>

<table><caption>**Table 2** (continued)</caption><thead><tr><th>Quantity name</th><th>Initial value</th><th>Units</th></tr></thead><tbody><tr><td>EmptyingSaliva</td><td>1.00</td><td>Dimensionless</td></tr><tr><td>kf</td><td>1.00</td><td>Dimensionless</td></tr><tr><td>kf_1</td><td>1.00</td><td>Dimensionless</td></tr><tr><td>SigmaBrainLymphatic</td><td>0.10</td><td>Dimensionless</td></tr><tr><td>TissuePartitionKidney</td><td>0.10</td><td>1/minute</td></tr><tr><td>ArterialToMuscleBloodFlow</td><td>45,000.00</td><td>Milliliter/hour</td></tr><tr><td>MuscleToVenousBloodFlow</td><td>45,000.00</td><td>Milliliter/hour</td></tr><tr><td>SigmaMuscleInterstitiumVascular</td><td>0.05</td><td>1/hour</td></tr><tr><td>MuscleLymphFlow</td><td>90.00</td><td>Milliliter/hour</td></tr><tr><td>SigmaMuscleLymphatic</td><td>0.80</td><td>Dimensionless</td></tr><tr><td>SigmaAdiposeInterstitiumVascular</td><td>0.05</td><td>1/hour</td></tr><tr><td>ArterialToAdiposeBloodFlow</td><td>15,600.00</td><td>Milliliter/hour</td></tr><tr><td>AdiposeToVenousBloodFlow</td><td>15,600.00</td><td>Milliliter/hour</td></tr><tr><td>SigmaAdiposeLymphatic</td><td>0.80</td><td>Dimensionless</td></tr><tr><td>AdiposeLymphFlow</td><td>31.20</td><td>Milliliter/hour</td></tr><tr><td>GutToVenousBloodFlow</td><td>66,000.00</td><td>Milliliter/hour</td></tr><tr><td>ArterialToGutBloodFlow</td><td>66,000.00</td><td>Milliliter/hour</td></tr><tr><td>SigmaGutInterstitiumVascular</td><td>0.05</td><td>1/hour</td></tr><tr><td>SigmaGutLymphatic</td><td>0.80</td><td>Dimensionless</td></tr><tr><td>GutLymphFlow</td><td>132.00</td><td>Milliliter/hour</td></tr></tbody></table>

A and B processed for sensitivity analysis (see “Methods” section). The main contributors to venous blood levels were muscle and gut tissues, with minor contributions by adipose, lung, and skin. After BBB disruption, sensitivity analysis pointed to brain sources as primary contributors to venous levels.

Since only a few organs contributed to the overall venous signal, we developed Model 2 based on three organs (muscle, adipose, and gut) plus the kidney and a “virtual” urine container mimicking the bladder (Fig. 4). Salivary production was also added to the model. The main difference between the two models is the description of brain S100B movements within and outside the brain parenchyma. For the brain, three compartments were used: vascular (i.e., the cerebrovascular circulation), interstitial (the brain extracellular space), and glia, referring to astrocytes, the primary cell type expressing S100B in the body (Fig. 5). The correspondence of the model with brain physiology is shown in Fig. 5B. In addition to the arterial influx and venous efflux, a glymphatic distribution process draining into Central lymph was added to the model. Another difference in Model 2 is that the structure of the organs and S100B movements within was based on the reflection coefficient (Sigma) rather than two-pore theory calculations (see “Methods” section).

We ran a sensitivity analysis for steady-state values of Central lymph, Arterial blood S100B, and Venous blood
shown in Fig. 8A. Note the small but measurable contri-
bution of glymphatic drainage to the venous signal.

For the sensitivity analysis shown in Fig. 6, we used an
interstitial concentration of S100B of 10 ng/ml. This value
is of course central to the model since it governs the lev-
els of S100B in peripheral organs and blood under nor-
mal conditions or after BBBD or trauma. We rerun the
simulation and sensitivity analysis with a low value of
interstitial S100B and compared the results with what
obtained with 10 ng/ml. The results are shown in Fig. 7;
A) refers to 10 and B) to 1 ng/ml S100B in the brain inter-
stitium. Note that no significant differences were seen in
overall sensitivity analysis. Similarly, we run a simulation
of venous values under various conditions using these
two values of interstitial S100B (compare Fig. 8A to C1)
to demonstrate an overall reduction of signal in venous
blood at low concentrations of S100B, as expected. Addi-
tional results for sensitivity analysis are shown in Addi-
tional file 1: Fig. S1, Additional file 2: Fig. S2, Additional
file 3: Fig. S3.

We previously measured S100B in several periph-
eral organs (see [3, 4]) and assigned these values as ini-
tial conditions for the simulations presented herein
(see Tables 1 and 2). We tested the hypothesis that the
levels measured in peripheral organs lacking mRNA for
S100B were due to diffusion of S100B from the blood.
Figure 8B shows these changes with 10 ng/ml interstitial
S100B, while Fig. 8C2 refers to 1 ng/ml. We started the
initial conditions with all organ levels set arbitrarily at 0
to test the extent of organs’ uptake of circulating S100B.
Note (Fig. 8B1) the increase in S100B due to the vascu-
lar uptake over a long period of control conditions (BBB_
Index=0). Figure 8B2 shows the contribution to the
venous levels of glymphatics and when the BBB_Index is
set to 0. The data in CD1-D3 show the changes in vari-
ables when BBBD and trauma were modeled. We then
studied the changes in several compartments (Fig. 8D)
under the same conditions. Note the effects of BBB dis-
ruption (D2) and trauma (D3) to organs and blood. Thus,
peripheral organs take up S100B from the circulation to,
in turn, contribute to blood levels. The amount of organs’

uptake of S100B depends on the assumed interstitial con-
centration of S100B.

The impact of glomerular filtration rate (GFR) and
urine formation on blood S100B levels was evaluated
(Fig. 9). Under intact or BBBD conditions, GFR greatly
influenced the levels of S100B measured in blood, lym-
phatics, urine, and organs. When setting GFR to zero
(Fig. 9A), we found a profound effect of kidney excre-
tion on both organ (left panel) and biological fluids (right
panel). In Fig. 9B, GFR was set at 10 or 100 while also
varying BBB_Index from 0 to 1.

Recent reports used salivary S100B and compared
its values to venous blood levels [33]. We simulated the
passive extravasation of arterial blood to form crevicu-
lar fluid [45], see Fig. 10. The levels of S100B in saliva, at
steady state, were larger than those in blood when blood
flow to saliva was adjusted to 4 ml (upper end of physi-
ologic levels [46]).

We formulated the hypothesis that after BBB dis-
ruption the half-life of S100B in blood is determined
in part by the availability of S100B in the brain inter-
stitium. This was tested as follows (Fig. 11). We simu-
lated a BBBD after reaching steady state at two time
points (arrows in Figure). Note that a broad range of
BBB_Index was explored (indicated in Figure). Also,
note that when the *Trauma_Index* was 0, the second
BBBD episode had little effect on S100B, unless the first
BBBD was minimal (0.1). We then repeated the simula-
tion with *Trauma_Index* set to 1. The secondary BBBD
response was restored parallel to a decreased depletion
of interstitial S100B in the brain (not shown). This sug-
gests that levels of S100B in the interstitium of the brain
are in part responsible for the time-dependent changes
in S100B in blood.

Discussion

The simulation efforts presented herein revealed sev-
eral surprising findings requiring human trials or animal
experiments to be confirmed or refuted. These unex-
pected results are listed in separate paragraphs below.

(See figure on next page.)

**Fig. 6** Sensitivity analysis for model 2. Steady-state conditions refer to the sensitivity measured at 10 h of simulation. The x axis in the heatmap plots specifies the sensitivity inputs. Y axes show the normalized sensitivity of venous, arterial, and lymphatic fluid obtained by varying the values of S100B in the compartments indicated by the numbers. The values of the outputs were integrated over time. We used global sensitivity analysis to interpret the impact of S100B levels in various organs on venous, lymphatic, and arterial biomarker levels. The key for the numbers at the bottom of each panel is 1: Venous blood; 2: Vascular compartment brain; 3: Vascular compartment adipose; 4: Vascular compartment gut; 5: Vascular compartment muscle; 6: Interstitium adipose; 7: Interstitium brain; 8: Interstitium gut; 9: Interstitium muscle; 10: Arterial blood; 11: Central lymph. When the brain is isolated from the periphery, and the only source of S100B available is the content of peripheral organs (top panel), the gut is the chief controller of body fluids S100B. However, when a communication brain to periphery is established via glymphatic drainage, the brain becomes the most influential organ for circulating S100B. This remains true after BBBD and the opening of the communication between the astrocyte content of S100B and the interstitium in the brain. For the brain interstitium S100B, in this simulation we used 10 ng/ml a concentration between CSF values (~ 3 ng/ml) and the measured interstitial value reported in [55]. See also Additional file 1: Fig. S1, Additional file 2: Fig. S2, Additional file 3: Fig. S3
<figure><img src="image_6.png" /><figcaption>**Fig. 6** (See legend on previous page.)</figcaption></figure>
<figure><img src="image_7.png" /></figure>

## Extracranial sources

It is well known that the distribution of S100B protein is not restricted to the brain. Several extracranial sources have been hypothesized to contribute to the blood levels used clinically [4, 5, 13, 15, 18]. Our simulation revealed that the contribution of “usual suspects” skin, and adipose tissue is relatively minor compared to the impact of gut and muscle release of S100B (Figs. 3 and 6). The possible explanation of these levels of S100B may depend on either local synthesis or uptake from blood. Since mRNA for S100B is lacking in gut tissue [4], the second explanation was tested (Fig. 7B1) by running a simulation where the initial values for organs’ S100B were arbitrarily set to zero while allowing for glymphatic-mediated contribution to blood (see also below). At steady state, peripheral tissues were loaded with levels of S100B comparable to those measured in vivo (Fig. 8A, B1) [3]. While the effects of glymphatics on venous levels was small (Fig. 8A), a prolonged stimulation (100 h) allowed to unveil a powerful effect on organs’ levels of S100B (Fig. 8B1), suggesting that glymphatic connection between brain and periphery was sufficient to load previously depleted organs with the biomarker. Note (B2) that when organs’ levels, BBBD, and glymphatics were set to 0, no venous signal was seen, suggesting that these three parameters are the exclusive contributors to steady-state organs’ levels of S100B. Thus, the most parsimonious explanation for the peripheral presence of S100B is the uptake of circulating protein, as

also shown in an animal model [4]. Conversely, these levels remain relatively stable once achieved until an event, such as BBB disruption occurs (Fig. 8). We also run simulations with low levels of interstitial S100B in brain (1 in lieu of 10 ng/ml). No major qualitative differences were seen at lower levels, albeit the responses to BBBD and trauma were reduced.

The question of whether peripheral levels of S100B contribute to the venous signal was answered by simulating control conditions or by adding BBB disruption events (Fig. 8A). At pre-BBBD time points, the brain influenced the blood signal via a mechanism involving glymphatic circulation (see above). BBBD (Fig. 8D2) and the trauma index in D3 dwarfed the control changes in S100B caused by other sources (organs and glymphatics). Sensitivity analysis showed that gut and muscle, but not kidney or adipose tissue, influence venous levels pre-BBBD, but brain contribution dominates after BBB disruption (Fig. 7 and Additional file 1: Fig. S1, Additional file 2: Fig. S2, Additional file 3: Fig. S3). Our results have thus shown that circulating S100B released by glymphatics in lymph and venous fluids is a likely scenario explaining tissue levels in the absence of transcription in peripheral organs. Additional BBB disruption did not increase S100B in organs (Fig. 8D1–D3). We have also shown that serum levels of S100B are only marginally affected by the release of S100B from organs, since the increased venous S100B never approached the cut-off
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /><figcaption>**Fig. 8** A Effects of glymphatics, BBB disruption, and the trauma index on venous blood S100B. Note that activating release from astrocytes (*Trauma_Index*) prevents return of S100B to pre-BBBD values, suggesting that parenchymal S100B is involved in the half-life of S100B. **B1** and **B2** "Filling" of empty organs (S100B set at 0 in adipose, muscle, and gut tissues) after a 100-h simulation with glymphatic communication between brain interstitial S100B and periphery. **B2** Changes in venous S100B in empty or full organs in the presence or absence of glymphatic contribution. When no S100B is available (*blue line*) venous levels are clamped at 0. When S100B in peripheral organs and glymphatic communication are present, a slight transient increase in venous S100B is seen. This is amplified by "opening" the BBB and establishment of a communication between glial cells reservoir and brain interstitial S100B. The data show that peripheral levels in organs can derive from brain reservoirs. **C1** Simulation identical to **A** but at a lower concentration of brain interstitial S100B. **C2** Simulation as in **B1** but with low levels of brain interstitial S100B. **D** Time course of S100B changes in various compartments under normal (**D1**), BBB disruption (**D2**), and trauma (**D3**).</figcaption></figure>
value of 0.1 or 0.15 nanogram/milliliter, which is the clinical ceiling for control subjects [43, 47].

### Do glymphatics contribute to biomarker blood levels?

Please note that as a semantic and scientific explanation of how brain effluxes solutes is still in progress [48–51], the term glymphatics is used here simply as a conveyor of the concept of brain clearance and not as an endorsement of a particular hypothesis. The “glymphatic flow” (in ml/h) may be paravascular or not; the only assumption in the model is that a flux from brain to lymphatic system exists. It was suggested that the primary source of S100B after traumatic brain injury is the brain’s communication with blood via glymphatic drainage [36]. We found no significant contribution of glymphatics to the overall signal in blood after BBBD (Figs. 6 and 8D1–D3). However, a small steady-state contribution of glymphatics to the pre-BBBD signal was observed (Fig. 8A). This contribution was however sharply decreased by lowering the interstitial brain S100B levels to 1 ng/ml (Fig. 8C1). At 10 ng/ml interstitial brain values of S100B, this finding suggests a continuous “trickle” of brain protein from the brain extracellular space into the blood via lymphatic drainage under physiological conditions. If this were the case, one expects that levels in blood will continuously increase, which is not what has been shown in human subjects. A fraction of what is being released from the brain is taken up by peripheral tissues as discussed above, but the factor that fully counteracts this constant source of S100B is kidney excretion of S100B (Fig. 10). In fact, when GFR was set to zero (in A), a constant increase was observed in peripheral fluids and organs. Thus, an equilibrium exists between glomerular filtration of small molecular weight protein [3] and S100B release from the brain interstitium via the glymphatic system. This finding predicts that patients with reduced glomerular filtration rate may have elevated levels of S100B (and other biomarkers) in the absence of a BBB contribution. A recent study [31] has shown that a constant source of S100B from brain to blood exists. This is an indirect validation of our modeling effort.

### Effects of parenchymal trauma on biomarker's levels

Brain damage and BBB disruption contribute to the overall levels of S100B in blood [34, 35]. However, in clinical practice is impossible to dissect out the contributions of these two factors independently. We have developed a subroutine in our software model that allows us to quantify and describe these two sources of blood S100B (Fig. 11). We

ran a simulation where levels of blood S100B in response to two BBB disruption events were monitored. When the *Trauma_index* was set to zero (no contribution of cellular release of S100B on peripheral or brain interstitial levels), we noted that the second BBB disruption episode did not cause an increase of S100B in venous blood unless a minimal level of disruption (*BBB_index* = 0.1) was used for the first event. We monitored the reserve of S100B sources in the interstitium to show that depletion of interstitial S100B occurred after the first, more significant, episodes. Therefore, the subsequent BBB “opening” was consequential only if a minimal depletion of S100B occurred during the first episode. When the *Trauma_index* was set to 1, replenishment of S100B in the extracellular space of the brain was repristinated, allowing for S100B release after the second BBB disruption event. This is a potentially important finding since it suggests that astrocytic sources of S100B are crucial in controlling the extent and duration of S100B during BBB disruptions.

### Comparison with existing models of blood biomarkers

In addition to our own prototype model [3, 32, 52], an effort to mimic biomarkers’ fate after TBI has been published [53]. The Authors use a much simplified, one-compartment model derived from oral absorption of therapeutic drugs. The limitations we found to be most relevant compared to the present study are: (1) Lack of distribution variables. Because only one compartment is used, the marker undergoes only blood distribution and thus disallows understanding of the impact of extracranial sources or the uptake of the marker by organs; (2) The model has only one path for the biomarker to leave the brain, ignoring glymphatic drainage; (3) Being a single-compartment model, there is no effort to reproduce organ size (including the brain) or cerebral and organ blood flow; and (4) The excretion data are presented only as a means to balance brain release.

### Brain levels of S100B

We initially used a middle-of-the-road concentration of S100B in brain (10 ng/ml). This value is supported by a recent paper where interstitial S100B levels were measured in brain slices [54]. Much higher levels have been measured after stroke and TBI [55, 56]. The pathological levels of S100B in brain tissue were modeled by the *Trauma_index*, which provides a replenishment of brain S100B by release from a reservoir with 50 ng/ml S100B. We already presented and discussed the outcome of trauma on S100B

(See figure on next page.)

**Fig. 9** **A** Effects of removing kidney filtration from the model. Note the increase in muscle S100B and venous levels, showing that glomerular filtration rate controls peripheral levels and kinetic behavior of S100B. **B** Effect of varying glomerular filtration on S100B. Note the drop of venous S100B with increased glomerular filtration rate and the lack of effect of GFR on brain interstitial levels and gut S100B. Also note organ-dependent changes in S100B with low or high GFR
A

<figure><img src="image_16.png" /></figure>

B

<figure><img src="image_17.png" /><figcaption>**Fig. 9** (See legend on previous page.)</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>**Fig. 10** Salivary levels in control or after BBB disruption. Note the delayed progression of salivary S100B towards steady state. Swallowing of saliva was not modeled for clarity and because nothing is known about the half-life of S100B in saliva. See refs. [2, 32, 33]. Saliva production was modeled at 4 and 2 ml/min</figcaption></figure>

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /><figcaption>**Fig. 11** Effect of the previous extent of BBB disruption on secondary BBB insults. Note that a delayed opening of the BBB after an earlier disruption translates into venous S100B levels that are inversely proportional to the extent of the earlier episode. In other words, a supramaximal increase in BBB release of S100B appears to deplete the brain sources, allowing only a minimal release of S100B by a subsequent episode. See sections "Results" and "Discussion"</figcaption></figure>

<figure><img src="image_21.png" /></figure>

levels after BBBD. We, however, also explored the possibility that under normal conditions, S100B levels in brain interstitium are equal to those typically reported for cerebrospinal fluid [57, 58]. When results with 10 vs 1 ng/ml were analyzed, no qualitative differences were found in terms of sensitivity or overall S100B dynamics (Fig. 7, Additional file 1: Fig. S1, Additional file 2: Fig. S2, Additional file 3: Fig. S3). However, the responses to BBBD and trauma were greatly reduced (compare 8A to 8C1). Thus,

the results from our simulation were independent from the levels of S100B in brain used as initial conditions.

## Limitations

The main limitation of our study is that we did not attempt to adapt the model to existing data on S100B, except when using realistic quantities of S100B in peripheral organs and a comparison of data with a previously established control value ceiling. All the data
sets available to us report S100B values in venous blood in individuals affected by a certain pathology or control subjects. In a previous work [3, 32] we focused on these pathophysiologic conditions. In the current study, the independent variable is time. To our knowledge, only a few studies reported the time course of S100B in blood; usually, only 2–3 time points were published. This makes it impossible to directly validate our model with existing data. We also used blood flow and volumetric data from the literature and accepted equilibrium values derived from the simulation of protein distribution after local injection [37]. Therefore, our modeling effort was not geared toward reproducing existing data but rather to allow for a discovery process of mechanisms that human subject-derived data make impossible to study.

## Future directions

The open-source format of the software developed herein (available at https://www.mathworks.com/matlab/bcentral/fileexchange/106145-diagnostic-pbpk-model-for-s100b?tid=srchtitle_Damir_2) will enable other researchers to adapt the core model to other situations and answer questions perhaps related to other biomarkers. Future developments will add the effect of molecular size (molecular weight and radius; see) on the movement across different compartments. This was already done in an older version of this model [3]. The main changes due to molecular size are likely to affect kidney filtration, passage across the BBB, and the overall kinetic properties of the marker. As new markers of CNS function are unveiled, we will focus on the physical and chemical properties of these proteins related to the voyage across organs and biofluids. In addition, we will develop a model that considers other biological variables, such as sex and age [3]. Finally, the model based on human subjects can be allometrically manipulated to include laboratory animals which are often used as surrogate experimental targets.

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s12987-022-00329-9.

**Additional file 1: Fig. S1.** Comparison of sensitivity analysis at two levels of brain interstitial S100B. The plot refers to time-dependent sensitivity. Note that unexplained variance was 0 in both simulations, demonstrating that the changes shown explain the variance of the model simulation. First and total order Sobol indices for model responses are shown.

**Additional file 2: Fig. S2.** Comparison of sensitivity analysis at two levels of brain interstitial S100B. The plot refers to time-dependent sensitivity. Note that unexplained variance was 0 in both simulations, demonstrating that the changes shown explain the variance of the model simulation. First and total order Sobol indices for model responses are shown.

**Additional file 3: Fig. S3.** Bar graph of the results shown in 1 and 2. See [40].

## Acknowledgements

We would like to thank Mr. Aaron Dadas and Mr. Jolewis Washington for helping with earlier versions of this model.

## Author contributions

DJ and RM constructed the model and tested it performance. NM and DB provided data for validation steps. All authors read and approved the final manuscript.

## Funding

ANR/Era-Net NeuVasc, ANR EpiCatcher (NM).

## Availability of data and materials

The model is currently available online.

## Declarations

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

All authors approved this submission.

### Competing interests

DJ has shares and IP in FloTBI a company centered around S100B as a marker of BBB disruption. The other authors report no conflict.

### Author details

¹FloTBI Inc., Cleveland, OH, USA. ²Department of Physiology and Biophysics, Case Western Reserve University, Cleveland, OH, USA. ³Laboratory of Cerebrovascular and Glia Research, Department of Neuroscience, Institute of Functional Genomics (UMR 5203 CNRS - U 1191 INSERM), University of Montpellier, Montpellier, France. ⁴Neurovascular Research Laboratory, Faculty of Life Sciences and Education, University of South Wales, Newport, UK.

Received: 6 April 2022 Accepted: 19 April 2022

Published online: 11 May 2022

## References

1. Segraves JM, Frishman WH. Highly sensitive cardiac troponin assays: a comprehensive review of their clinical utility. Cardiol Rev. 2015;23(6):282-9.

2. Janigro D, Bailey DM, Lehmann S, Badaut J, O'Flynn R, Hirtz C, et al. Peripheral blood and salivary biomarkers of blood-brain barrier permeability and neuronal damage: clinical and applied concepts. Front Neurol. 2020;11:577312.

3. Dadas A, Washington J, Marchi N, Janigro D. Improving the clinical management of traumatic brain injury through the pharmacokinetic modeling of peripheral blood biomarkers. Fluids Barriers CNS. 2016;13(1):21.

4. Bargerstock E, Puvenna V, Iffland P, Falcone T, Hossain M, Vetter S, et al. Is peripheral immunity regulated by blood-brain barrier permeability changes? PLoS ONE. 2014;9(7): e101477.

5. Pham N, Fazio V, Cucullo L, Teng Q, Biberthaler P, Bazarian JJ, et al. Extracranial sources of S100B do not affect serum levels. PLoS ONE. 2010;5(9): e12691.

6. Preston E, Webster J, Small D. Characteristics of sustained blood-brain barrier opening and tissue injury in a model for focal trauma in the rat. J Neurotrauma. 2001;18(1):83-92.

7. Kanner AA, Marchi N, Fazio V, Mayberg MR, Koltz MT, Siomin V, et al. Serum S100beta: a noninvasive marker of blood-brain barrier function and brain lesions. Cancer. 2003;97(11):2806-13.

8. Marchi N, Rasmussen P, Kapural M, Fazio V, Kight K, Mayberg MR, et al. Peripheral markers of brain damage and blood-brain barrier dysfunction. Restor Neurol Neurosci. 2003;21(3-4):109-21.

9. Marchi N, Cavaglia M, Bhudia S, Hallene K, Janigro D. Peripheral markers of blood-brain barrier damage. Clin Chim Acta. 2004;342(1-2):1-12.
10. Falcone T, Fazio V, Lee C, Simon B, Franco K, Marchi N, et al. Serum S100B: a potential biomarker for suicidality in adolescents? PLoS ONE. 2010;5(6): e11089.

11. Bouzat P, Francony G, Declety P, Brun J, Kaddour A, Renversez JC, et al. Can serum protein S100 beta predict neurological deterioration after moderate or minor traumatic brain injury? Annales Francaises D Anesthesia et de Reanimation. 2009;28(2):135-9.

12. Hasselblatt M, Mooren FC, von Ahsen N, Keyvani K, Fromme A, Schwarze-Eicker K, et al. Serum S100beta increases in marathon runners reflect extracranial release rather than glial damage. Neurology. 2004;62(9):1634-6.

13. Korfias S, Stranjalis G, Papadimitriou A, Psachoulia C, Daskalakis G, Antsaklis A, et al. Serum S-100B protein as a biochemical marker of brain injury: a review of current concepts. Curr Med Chem. 2006;13(30):3719-31.

14. Papa L, Silvestri S, Brophy GM, Giordano P, Falk JL, Braga CF, et al. GFAP out-performs S100beta in detecting traumatic intracranial lesions on computed tomography in trauma patients with mild traumatic brain injury and those with extracranial lesions. J Neurotrauma. 2014;31(22):1815-22.

15. Savola O, Pyhtinen J, Leino TK, Siitonen S, Niemela O, Hillbom M. Effects of head and extracranial injuries on serum protein S100B levels in trauma patients. J Trauma. 2004;56(6):1229-34.

16. Thelin EP, Jeppsson E, Frostell A, Svensson M, Mondello S, Bellander BM, et al. Utility of neuron-specific enolase in traumatic brain injury; relations to S100B levels, outcome, and extracranial injury severity. Crit Care. 2016;20:285.

17. Anderson RE, Hansson LO, Nilsson O, Dijlai-Merzoug R, Settergren G. High serum S100B levels for trauma patients without head injuries. Neurosurgery. 2001;48(6):1255-8.

18. Posti JP, Hossain I, Takala RS, Liedes H, Newcombe V, Outtrim J, et al. Glial fibrillary acidic protein and ubiquitin C-terminal hydrolase-L1 are not specific biomarkers for mild CT-negative traumatic brain injury. J Neurotrauma. 2017. https://doi.org/10.1089/neu.2016.4442.

19. Zhang Y, Zhu J, Xu H, Yi Q, Yan L, Ye L, et al. Time-dependent internalization of S100B by mesenchymal stem cells via the pathways of clathrin- and lipid raft-mediated endocytosis. Front Cell Dev Biol. 2021;9: 674995.

20. Pham N, Fazio V, Cucullo L, Teng Q, Bibberthaler P, Bazarian JJ, et al. Extracranial sources of S100B do not affect serum levels. PLoS ONE. 2010. https://doi.org/10.1371/journal.pone.0012691.

21. Steiner J, Schiltz K, Walter M, Wunderlich MT, Keilhoff G, Brisch R, et al. S100B serum levels are closely correlated with body mass index: an important caveat in neuropsychiatric research. Psychoneuroendocrinology. 2010;35(2):321-4.

22. Kawata K, Rubin LH, Takahagi M, Lee JH, Sim T, Szwanki V, et al. Sub-concussive impact-dependent increase in plasma S100beta levels in collegiate football players. J Neurotrauma. 2017;34(14):2254-60.

23. Zonner SW, Ejima K, Bevilacqua ZW, Huibregtse ME, Charleston C, Fulgar C, et al. Association of increased serum S100B levels with high school football subconcussive head impacts. Front Neurol. 2019;10:327.

24. Tagge CA, Fisher AM, Minaeva OV, Gaudreau-Balderrama A, Moncaster JA, Zhang XL, et al. Concussion, microvascular injury, and early tauopathy in young athletes after impact head injury and an impact concussion mouse model. Brain. 2018;141(2):422-58.

25. Schulte S, Podlog LW, Hamson-Utley JJ, Strathmann FG, Struder HK. A systematic review of the biomarker S100B: implications for sport-related concussion management. J Athl Train. 2014;49(6):830-50.

26. Michetti F, Bruschettini M, Frigiola A, Abella R, Giamberti A, Marchese N, et al. Saliva S100B in professional sportsmen: high levels at resting conditions and increased after vigorous physical activity. Clin Biochem. 2011;44(2-3):245-7.

27. Watson P, Shirreffs SM, Maughan RJ. Blood-brain barrier integrity may be threatened by exercise in a warm environment. Am J Physiol Regul Integr Comp Physiol. 2005;288(6):R1689-94.

28. Watson P, Black KE, Clark SC, Maughan RJ. Exercise in the heat: effect of fluid ingestion on blood-brain barrier permeability. Med Sci Sports Exerc. 2006;38(12):2118-24.

29. Schulte S, Schiffer T, Sperlich B, Kleinoder H, Holmberg HC. Serum concentrations of S100B are not affected by cycling to exhaustion with or without vibration. J Hum Kinet. 2011;30:59-63.

30. Koh SX, Lee JK. S100B as a marker for brain damage and blood-brain barrier disruption following exercise. Sports Med. 2014;44(3):369-85.

31. Bailey DM, et al. Hypoxemia promotes blood-brain barrier destabilization of the neurovascular unit during extreme apnea in humans. J Cereb Blood Flow Metab. 2022. https://doi.org/10.1177/0271678X221075967.

32. Dadas A, Janigro D. The role and diagnostic significance of cellular barriers after concussive head trauma. Concussion. 2018;3(1):Cnc53.

33. Janigro D, Kawata K, Silverman E, Marchi N, Diaz-Arrastia R. Is Salivary S100B a biomarker of traumatic brain injury? A pilot study. Front Neurol. 2020;11:528.

34. Marchi N, Cavaglia M, Fazio V, Bhudia S, Hallene K, Janigro D. Peripheral markers of blood-brain barrier damage. Clin Chim Acta. 2004;342(1-2):1-12.

35. Marchi N, Rasmussen PA, Kapural M, Fazio V, Cavaglia M, Janigro D. Peripheral markers of brain damage and blood-brain barrier dysfunction. Restor Neurol Neurosci. 2003;21(3-4):109-21.

36. Plog BA, Dashnaw ML, Hitomi E, Peng WG, Liao YH, Lou NH, et al. Biomarkers of traumatic injury are transported from brain to blood via the glymphatic system. J Neurosci. 2015;35(2):518-26.

37. Gill KL, Gardner I, Li L, Jamei M. A bottom-up whole-body physiologically based pharmacokinetic model to mechanistically predict tissue distribution and the rate of subcutaneous absorption of therapeutic proteins. AAPS J. 2016;18(1):156-70.

38. Heizmann CW, Fritz G, Schafer BW. S100 proteins: structure, functions and pathology. Front Biosci. 2002;7:D1356-68.

39. Sepp A, Meno-Tetang G, Weber A, Sanderson A, Schon O, Berges A. Computer-assembled cross-species/cross-modalities two-pore physiologically based pharmacokinetic model for biologics in mice and rats. J Pharmacokinet Pharmacodyn. 2019;46(4):339-59.

40. Tiemann CA, Vanlier J, Oosterveer MH, Groen AK, Hilbers PA, van Riel NA. Parameter trajectory analysis to identify treatment effects of pharmacological interventions. PLoS Comput Biol. 2013;9(8):e1003166.

41. Shah DK, Betts AM. Towards a platform PBPK model to characterize the plasma and tissue disposition of monoclonal antibodies in preclinical species and human. J Pharmacokinet Pharmacodyn. 2012;39(1):67-86.

42. Morquette P, Verdier D, Kadala A, Fethiere J, Philippe AG, Robitaille R, et al. An astrocyte-dependent mechanism for neuronal rhythmogenesis. Nat Neurosci. 2015;18(6):844-54.

43. Iverson GL, Posti JP, Ohman J, Blennow K, Zetterberg H, Luoto TM. Reliability of serum S100B measurement following mild traumatic brain injury: a comparison of assay measurements from two laboratories. Brain Inj. 2020;34(9):1237-44.

44. Bouvier D, Duret T, Rouzaire P, Jabaudon M, Rouzaire M, Nourrisson C, et al. Preanalytical, analytical, gestational and pediatric aspects of the S100B immuno-assays. Clin Chem Lab Med. 2016;54(5):833-42.

45. Rahim MA, Rahim ZH, Ahmad WA, Hashim OH. Can saliva proteins be used to predict the onset of acute myocardial infarction among high-risk patients? Int J Med Sci. 2015;12(4):329-35.

46. Jasim H, Olausson P, Hedenberg-Magnusson B, Ernberg M, Ghafouri B. The proteomic profile of whole and glandular saliva in healthy pain-free subjects. Sci Rep. 2016;6:39073.

47. Okonkwo DO, Puffer RC, Puccio AM, Yuh EL, Yue JK, Diaz-Arrastia R, et al. Point-of-care platform blood biomarker testing of glial fibrillary acidic protein versus S100 calcium-binding protein B for prediction of traumatic brain injuries: a transforming research and clinical knowledge in traumatic brain injury study. J Neurotrauma. 2020. https://doi.org/10.1089/neu.2020.7140.

48. Hladky SB, Barrand MA. The glymphatic hypothesis: the theory and the evidence. Fluids Barriers CNS. 2022;19(1):9.

49. Abbott NJ, Pizzo ME, Preston JE, Janigro D, Thorne RG. The role of brain barriers in fluid movement in the CNS: is there a "glymphatic" system? Acta Neuropathol. 2018;135(3):387-407.

50. Dadas A, Washington J, Janigro D. Cerebral waste accumulation and glymphatic clearance as mechanisms of human neurological diseases. J Neurol Neuromed. 2016;1(7):15-9.

51. Plog BA, Nedergaard M. The glymphatic system in central nervous system health and disease: past, present, and future. Annu Rev Pathol. 2017. https://doi.org/10.1146/annurev-pathol-051217-111018.

52. Dadas A, Washington J, Diaz-Arrastia R, Janigro D. Biomarkers in traumatic brain injury (TBI): a review. Neuropsychiatr Dis Treat. 2018;14:2989-3000.

53. Azizi S, Hier DB, Allen B, Obafemi-Ajayi T, Olbricht GR, Thimgan MS, et al. A kinetic model for blood biomarker levels after mild traumatic brain injury. Front Neurol. 2021;12: 668606.
54. Sakatani S, Seto-Ohshima A, Shinohara Y, Yamamoto Y, Yamamoto H, Itohara S, et al. Neural-activity-dependent release of S100B from astrocytes enhances kainate-induced gamma oscillations in vivo. J Neurosci. 2008;28(43):10928–36.

55. Winter CD, Clough GF, Pringle AK, Church MK. Outcome following severe traumatic brain injury TBI correlates with serum S100B but not brain extracellular fluid S100B: an intracerebral microdialysis study. World J Neurosci. 2013;3:93–9.

56. Sen J, Belli A, Petzold A, Russo S, Keir G, Thompson EJ, et al. Extracellular fluid S100B in the injured brain: a future surrogate marker of acute brain injury? Acta Neurochir. 2005;147(8):897–900.

57. Hajdukova L, Sobek O, Prchalova D, Bilkova Z, Koudelkova M, Lukaskova J, et al. Biomarkers of brain damage: S100B and NSE concentrations in cerebrospinal fluid—a normative study. Biomed Res Int. 2015;2015: 379071.

58. Reiber H. Dynamics of brain-derived proteins in cerebrospinal fluid. Clin Chim Acta. 2001;310(2):173–86.

**Publisher’s Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_22.png" /></figure>