Stochastic Modelling of Solution Particle Movement: An Individual Case of Coupled Concentration Gradient Dependent and Independent Movements of Efavirenz

This work proposed a coupled model of diffusion. It adopted two forms of coupled movement, the interacting and non-interacting driven forms of movement of a solution particle of efavirenz concentration measured in blood plasma. Data from projected pharmacokinetics in a patient on efavirenz were used. A relationship between interacting and non-interacting diffusion was suggested through a stochastic differential equation. The solution particle with a small value of relative acceleration drift to its active neighbourhood was projected to have a corresponding high transport/interacting diffusion.


Introduction
The work attempted to combine two forms of diffusion starting with phenomenological approach motivated by Fick's laws and their mathematical solutions and the "random walk" of diffusing particles suggested by Robert Brown. This random walk of microscopic particles in suspension in a fluid has since been adopted as the "Brownian motion" in honour of Brown. The mathematical form of the Brownian motion was derived by Einstein [1]. The Fick's laws follow from Fourier's law of thermal conduction [1] [2] [3] [4].
The research encouraging modelling of random fluctuations found in PK/PD (Pharmacokinetic and Pharmacodynamic) relations has been suggested. It has been noted that there is an increasing need to extend the deterministic models which are currently favoured to models including a stochastic component in modelling pharmacological processes [5].
The transport diffusion is the concentration gradient dependent driven movement and describes the change in concentration of a solution particle. The transport diffusion can further be studied from Fick's laws by use of partial differential equations, however reseachers have used different variable spaces to study this flow [6] [7] [8]. Furthermore, in this work we tracked the transport diffusion by making use of the one to one function relationship between secondary saturation and concentration [9]. It is also noted that concentration is not a movement parameter as compared to secondary saturation. The self diffusion is a concentration gradient independent driven movement and describes the non-interacting movement of solution particle in the volume space. The observation of self diffusion component has been attributed to Brown [1].
The work has allowed for a proposition on the possibility to study the two movements in the 24 h dosing period that aid the process of diffusion by considering stochastic differential equations [5] [10]. It describes the movement of an interacting particle. Inferences were proposed on the possible relationship between the two movements (interacting and non-interacting). Numerical solutions on diffusional constituent movement characterisation were derived for a patient on efavirenz. This work presents a stochastic coupled diffusion model. This enables us to study the relationship between transport and self forms of diffusion. Stochastic coupled models have been applied on combined PK/PD processes [5]. In addition, in this case the coupling is in modelling the relationship of the two separate diffusional movements. A measure is developed to estimate the extent of gradient driven diffusion with the aid of a coupled stochastic differential equation. The Verhulst model is proposed to estimate the logistic relation obtained from modeling the self and transport forms of diffusion [11].

Methods
Simulated projected data on secondary saturation movement, time and concentration was taken from pharmacokinetic projections made on Patient P on 600 mg dose considered in Nemaura [9] [12]. Stochastic differential equations, Ordinary differential equations and linear and non-linear regression models are used to model diffusion in plasma. The following softwares, MATLAB and R were used for development of models. The solution of the Stochastic Differential Equations are approximated with the Euler-Maruyama scheme.

Deterministic Models for Concentration and Gradient Driven Movements
Initially there was modelling of the concentration profile of a solution particle of patient P in time. Letting ( ) x t model concentration of solution particle at time t and was given by, The process ( ) x t was sufficiently modelled by the following equation, We modelled the gradient driven movement ( ) The following was noted,

Stochastic Model
The secondary saturation movement followed a one to one relation with concentration and was thus used as compared to other forms of secondary movement that is convection and advection [9]. The process ( )( ) 1 y t h − modelled both the interaction and non-interacting diffusion movement of the solution particle in the volume space. The process ( ) y t modelled these two diffusion forms of the solution particle and was thus given by, where the process ( ) b t modelled the independent standard Brownian motion (Weiner process),

x t t t x t b t y uv
µ σ where, Furthermore, for and some C < ∞ the following conditions are satistified, Conditions 13 and 14 allowed for the existence and uniqueness of solutions for Equation (11) [10].

Modelling Concentration and the Concentration Gradient Driven Movement in the Absence of Self Diffusion
We estimated the parameters in Equation (1) for the concentration-time curve (Table 1).
Furthermore, we estimated the parameter values for the derived secondary saturation ( ) , F x t with respect to time and concentration in Patient P ( Table 2  and Table 3).   There was consideration of how the concentration ( )  (1) and (7) (Figure 1).
Informed by results herewith there was an adoption that transport diffusion was driven by concentration. Thus concentration could be used as a substitute parameter for transport diffusion.

Stochastic Model: Modelling the Solution Particle Movement with Variable Relative Acceleration (Drift)
The following case was considered where to have no transport diffusion. A solution particle with increasing relative acceleration drift to its neighbour would have less magnitude of extent with respect to transport diffusion. However, a solution particle with low relative drift to its neighbour had a high potential of transport diffusion. That is transport diffusion was more pronounced in neighbouring particles with relatively low acceleration drift. A relationship of the correlation between ( ) x t and ( ) y t , and the drift in self diffusion was proposed ( Figure 2).
An equation for the logistic decay relationship [11] was proposed, where k and l (decay) are constants (    It is important to note that the relationship between concentration change and solution particle movement was inferred to be predominantly affected by relative acceleration drift.

Conclusions
Other researchers have been able to correlate self and transport forms of diffusion from experimental observations [8]. Transport diffusion has been shown to thrive in neighbouring solution particles that have zero relative drift. The transport diffusion in a neighbourhood of a solution particle was projected to be affected by drift. Considering an increased level of self diffusion of a solution particle relative to its active neighbourhood, showed that it will give out less in terms of concentration material to be exchanged through gradient. A logistic relation was projected between this correlation and drift. An elaborate study of the dynamics of the logistic decay curve proposed was done in Gonze, and considering the "Malthusian parameter value" resulted in an unstable steady state and unstable cycles for the system derived for Patient P since it had a value greater than 1 [11]. Relatively slow movement of neighbouring solution particles was projected to enable optimum drug transfer through transport diffusion in the plasma.
Throughout this work, there was an inherent assumption that the particles were already in orientation that enabled transport diffusion. This was because for transport diffusion to occur the particles should be in such an orientation which enables progression of that process [8]. The extent of the transport diffu-sion of a solution particle was projected to decrease with increasing relative drift to an active neighbourhood.