Figure 1.
Mass fraction evolution across the membrane layer for different values of fKn for Kn = 0.1.
Citation: Eman Koshlaf, Andrew S Ball. Soil bioremediation approaches for petroleum hydrocarbon polluted environments[J]. AIMS Microbiology, 2017, 3(1): 25-49. doi: 10.3934/microbiol.2017.1.25
| [1] | Tianyuan Xu, Shanming Ji, Chunhua Jin, Ming Mei, Jingxue Yin . EARLY AND LATE STAGE PROFILES FOR A CHEMOTAXIS MODEL WITH DENSITY-DEPENDENT JUMP PROBABILITY. Mathematical Biosciences and Engineering, 2018, 15(6): 1345-1385. doi: 10.3934/mbe.2018062 |
| [2] | Wenjie Zhang, Lu Xu, Qiao Xin . Global boundedness of a higher-dimensional chemotaxis system on alopecia areata. Mathematical Biosciences and Engineering, 2023, 20(5): 7922-7942. doi: 10.3934/mbe.2023343 |
| [3] | Sunwoo Hwang, Seongwon Lee, Hyung Ju Hwang . Neural network approach to data-driven estimation of chemotactic sensitivity in the Keller-Segel model. Mathematical Biosciences and Engineering, 2021, 18(6): 8524-8534. doi: 10.3934/mbe.2021421 |
| [4] | Qianhong Zhang, Fubiao Lin, Xiaoying Zhong . On discrete time Beverton-Holt population model with fuzzy environment. Mathematical Biosciences and Engineering, 2019, 16(3): 1471-1488. doi: 10.3934/mbe.2019071 |
| [5] | Chichia Chiu, Jui-Ling Yu . An optimal adaptive time-stepping scheme for solving reaction-diffusion-chemotaxis systems. Mathematical Biosciences and Engineering, 2007, 4(2): 187-203. doi: 10.3934/mbe.2007.4.187 |
| [6] | Xu Song, Jingyu Li . Asymptotic stability of spiky steady states for a singular chemotaxis model with signal-suppressed motility. Mathematical Biosciences and Engineering, 2022, 19(12): 13988-14028. doi: 10.3934/mbe.2022652 |
| [7] | Tingting Yu, Sanling Yuan . Dynamics of a stochastic turbidostat model with sampled and delayed measurements. Mathematical Biosciences and Engineering, 2023, 20(4): 6215-6236. doi: 10.3934/mbe.2023268 |
| [8] | Lin Zhang, Yongbin Ge, Zhi Wang . Positivity-preserving high-order compact difference method for the Keller-Segel chemotaxis model. Mathematical Biosciences and Engineering, 2022, 19(7): 6764-6794. doi: 10.3934/mbe.2022319 |
| [9] | Changwook Yoon, Sewoong Kim, Hyung Ju Hwang . Global well-posedness and pattern formations of the immune system induced by chemotaxis. Mathematical Biosciences and Engineering, 2020, 17(4): 3426-3449. doi: 10.3934/mbe.2020194 |
| [10] | Marcin Choiński, Mariusz Bodzioch, Urszula Foryś . A non-standard discretized SIS model of epidemics. Mathematical Biosciences and Engineering, 2022, 19(1): 115-133. doi: 10.3934/mbe.2022006 |
Nanomedicine knows an increasing interest due to its estimated importance in efficient medication delivery [1]. Crucial medical topics deal with its application in several fields, such as overcoming the blood-brain barrier [2], anti-diabetic medication [3] and uptake in endothelial cells [4]. Nanomedicine deals with the novel chemical, physical or biological properties of medication at nanometric scale, as well as the application of nanotechnology aiming at making a medical diagnosis, treating or preventing the cause of illness. Therefore, nanomedicines have the ability, due to their particular size or surface functionality, to demonstrate physicochemical and/or biological properties that make them interesting for medication delivery, controlled medication release and/or improved medication transport through biological barriers, such as cell membranes, which would otherwise be much more difficult [5]. One of the key parameters for successful medication delivery is controlling or improving the permeation of medication through the destined place (e.g. vascular system, cell) in the human body [6,7,8]. When speaking about cell permeation and the body uptake efficiency of drugs into the body that it entails, many works take the apparent cell permeability, Papp, through monolayers of human epithelial cells as a widely accepted standard [9,10,11]. Values of Papp (often presented in medicine as a volumetric flux per unity of area) of such cells are also often used for the estimation of the blood-brain barrier [12,13]. It is useful to understand the apparent cell permeation as a series of parallel and serial resistances (from layer to another). This vision is based on the concept of solubility diffusion [14]. The permeability is seen here, therefore, not as an imposed force, but is rather passive in a way that only spontaneous diffusion governs the permeation process. This diffusion process is determined by a concentration gradient. The concentration gradient is classically described by Fick's law or, in case of sizes (within the membranes separating the monolayers) at nanoscale, by a more general evolution equation, describing the extended diffusion law, as will be discussed later. Active transport, as could be caused by transporter proteins, is not considered here. Although, in reality, chemicals/drugs could be transported by proteins (and thus have an effect on Papp), it is not taken into account here. It is well admitted that generally three pathways are possible through a cellular monolayer, i.e. in between the cells (the paracellular way), through the cell (considered here) and along the membrane solely (the lateral way) [15,16]. In this case, parallel resistances of the three pathways would apply. Here, we only consider the pathway through the cell. Following the pH partition hypothesis [17], it can be assumed that only the neutral form of the drug chemicals can pass the membrane resistances. Therefore, chemicals that are ionized should undergo an acid-base reaction in between the water layer and the membrane before they are able to pass the membrane as a neutral species. According to the aforementioned pH partition hypothesis, such a reaction can be considered to occur instantaneously, meaning that no kinetic hinderness (and therefore no additional resistance) needs to be considered.
It is often suggested that the diffusion through membranes are of non-Fickian nature. We propose that this is either because of non-local effects (the size of the membranes are of the same order of magnitude as the drug molecules or not much bigger) or because the relaxation time becomes important. As controlling and enhancing drug delivery is an important issue in the field of nanomedicine, next to in vitro experiments, modelling is crucial to understand the underlying mechanisms. In nanomedicine, the inertia of the diffusing particles is relevant in addition to the small size of the channels through which the medication passes, being in some cases of the same order of magnitude as the mean free path of the particles. Furthermore, purely microscopic models are often difficult to implement directly to real-life applications. Therefore, further progress in this field can be enhanced by models that are easily applicable, whilst keeping a qualitative efficient and understandable way to link to experiments in a multi-disciplinary way. Therefore, we propose here a model based on thermodynamic principles, where we extend them by incorporating effects, in a rigourous way, that account for the relaxation and size of the drug molecules with respect to the channels through which they diffuse. With the purpose to demonstrate the applicability of our new theory, we treat two cases: drug diffusion through a cell membrane and drug release out of a scaffold. For the first case, let us imagine the existence of two compartments of well-mixed water that are separated by a barrier, i.e. a membrane: the first containing the drug molecules as the donor and the second represents the inner cell or the receptor. The purpose is to follow the evolution of the mass fraction profile through the membrane at the beginning of the diffusion process, where we assume that the drug molecules (in that beginning process) concentration is that high that the mass fraction at the inlet of the membrane is constant. For the second case, we discuss drug release out of a scaffold and follow the concentration profile to assess the rapidity of the drug release. To account for the importance of the relaxation time of the drug molecules and the non-local effects, we propose the formulation of mass transport using Extended Non-Equilibrium Thermodynamics [18,19], whose main characteristic is to upgrade the flux and higher order fluxes to the rank of independent variables.
At micro and nanoscales, mass transport is mostly influenced by non-local effects and relaxation times. The classical Fick law
| J=−ρD∇c, | (1) |
relating the mass flux vector J to the mass fraction gradient ∇c, with ρ and D denoting the density and mass diffusion coefficient, respectively, is not applicable at short diffusional time (with respect to the relaxation time, i.e. large relaxation times) and small spatial scales. In order to account for large relaxation times, Fick's law can be generalized by a Cattaneo-like equation (in analogy to the one used for heat transport [20]) under the form
| τr∂tJ+J=−ρD∇c | (2) |
where τr designates the relaxation time of the mass flux and ∂t the partial time derivative. To account furthermore for non-local effects, Extended Non-Equilibrium Thermodynamics (see for underlying developments in previous work [18]) proposes a further generalization
| τr∂tJ+J=−ρD∇c+ℓ2∇2J | (3) |
It is convenient to reformulate (3) in dimensionless form. Therefore, the mass flux J is rescaled with ρDδ, the space coordinate with a characteristic length δ (standing also for the membrane thickness), the time t with the characteristic diffusion time τd=δ2D. Introducing the Knudsen number Kn≡ℓδ, a dimensional number M=τrτd and using the same notation for the non-dimensional quantities, Eq (3) becomes
| M∂tJ+J=−∇c+Kn2∇2J | (4) |
It should be noted here that by choosing τd as time scale, one implies that any change in M means a change in the relaxation time, which evokes non-local effects. This motivates proposing an adaptation for the expression for M. Indeed, kinetically speaking, one may define a relation D≡13ℓ2τr, in analogy to the thermal conductivity [1]. Using the definition of τd and Kn, one finds a relation M=13Kn2. However, in suspensions, these kinetic considerations are not always valid and differences could occur. Acknowledging that some partial proportionality is admitted between M and the Kn number (the relaxation time and the mean free path are related in that they both can cause non-local effects), this motivates to define rather M=fKn13Kn2. Note that using this definition does not change anything for the analysis to be performed (indeed, in principle, any number could be chosen for M in the framework of a dimensionless analysis conditioned by correct interpretation of the results). Moreover, fKn is a factor that accounts for the observation that in nanofluids the relaxation time is not necessarily fully proportional to the Kn number. This factor also allows investigating separately the effects of the relaxation time and the mean free path on drug diffusion.
We can see clearly differences between (3) and the conventional constitutional Fick's law for mass diffusion (1), but also some similarities. When we deal with systems where the characteristic size is found to be much larger than the mean free path of the molecules l, Eq (3) reduces to the Cattaneo-like equation (2). This would also correspond to neglect any higher order of mass fluxes. Now, when the changes in the mass flux occur in a time span that is much larger than the relaxation time, we obtain Fick's law. Let us stress that the latter issue depends on the physicochemical relation the medication molecule can have with the carrier fluid or the host environment. Other well-known methods that are used for solving models that deal with nanoscale systems are the phonon-Boltzmann equation, Monte Carlo simulations and Molecular Dynamics. The Monte Carlo simulations and Molecular Dynamics method are numerical methods that generally do not require higher order variables, but instead demand high computing times, whilst being less flexible for modelling various scales at a time, although some efforts have been performed concerning multi-scale molecular dynamics [21]. Nonetheless, the flexibility in scale remains less than the present approach. The phonon-Boltzmann equation also demands high computational time for solving a model. The aforementioned aspects are generally less of a problem when it comes to the model that is proposed in this work. It should be noted that the present development could be extended by introducing non-linear mass fluxes or higher-order (higher than the second order) fluxes. However, implementing such an extension, albeit theoretically not an impossible task (see [18]), would introduce numerical difficulties and higher computing time, which is an absent aspect of the present approach. Furthermore, it would also loose the ability to propose analytical solutions for some particular cases, offering a pedagogical presentation allowing easier interpretations. Non-linear contributions of the mass flux (at first and higher orders) could be interesting to consider when the length scale becomes much smaller than the mean free path of the molecules. For the purposes of nanomedicine applications, this would in many cases lead to unnecessary complications. Indeed, for many cases a linear mass flux limited to second-order mass fluxes (which leads to the term ∇2J in (3)) is sufficient. In comparison, it already appeared to be sufficient for heat transport applications [22] and for predicting nanoflow properties [23].
From typical clinical observations that lead to data in the literature, it appears that one of the most important issues to be studied for improving drug delivery are through enhancing the permeability. In porous media, the concentration gradient can be negligible with respect to the flux, which is often the case for nanoporous flow that is pressure-driven by some external pressure. There are some cases in the human body that are comparable from a conceptual point of view, such as oxygen and carbon dioxide transfer from the lungs to the veins within the alveoles. Nonetheless, porous flow within cell membrane pores presents often the difference that, due to the absence of externally imposed pressure, the concentration gradient can be quite important with respect to the flux. In case of passive permeability, one speaks often of diffusional processes. As the diffusion process through the membrane is often the limiting pathway, the flux across the membrane will undergo non-local effects. Flux variation in the perpendicular direction to the flux across the membrane can be neglected, assuming a well-dispersed medication solution at the entrance of the whole membrane. One would then use (3) to describe such a process, reduced to a one-dimensional equation, replacing the ∇ operator by the derivative ∂y. However, the concentration gradient is not simply a gradient of the concentrations from both sides of the membrane, since concentration jumps can exist at the entry or exit of those membranes due to lipid bilayer interaction or wettability effects. The gradient that counts is the one within the membrane. Therefore, it is important to consider the fact that flux by diffusion across a membrane (even if it were uniform) depends on the solubility. As there are different solubilities at both sides of the membrane and within the membrane as well, we can define partition coefficients ϑin and ϑout for the entry and exit of the membrane, respectively:
| ϑin=Nmb,inNin | (5) |
| ϑout=Nmb,outNout | (6) |
where Nmb, in and Nmb, out are the concentrations within the membrane from the entry and exit sides, respectively. Also, Nin and Nout stand for the same concentrations, but at the medium sides (outside the membrane), respectively. It is more convenient to write (5) and (6) in mass fractions, c, where the subscripts mean the same thing as for the concentrations. Assuming ϑin = ϑout≡ϑ, we have
| ϑ=cmb,in(1−cin)cin(1−cmb,in) | (7) |
| ϑ=cmb,out(1−cout)cout(1−cmb,out) | (8) |
where ϑ is to be determined experimentally as a material property. We can readily see that for a higher partition coefficient, the concentration gradient within the membrane is then larger than on the outer sides of the membrane, which results into a higher membrane solubility and thus a higher permeation. To model the whole process, next to the one-dimensional version of Eq (3), the mass fraction evolution is given by the classical species balance equation, which in dimensionless form gives
| ∂tc=−∇⋅J | (9) |
Since the thickness of a cell membrane is much smaller than its size in the other directions, we consider a one-dimensional model through the membrane with thickness δ, assuming that in all directions other than across the membrane the variables remain constant. The entrance of the membrane (at the drug donor side) is denoted by the dimensionless coordinate y = 0 and the exit (at the drug receptor side) by y = 1. The transport of the drug molecules through a cell membrane will then be modeled by the one-dimensional versions of Eqs. (4) and (9)
| M∂tJ+J=−∂yc+Kn2∂2yJ | (10) |
| ∂tc=−∂yJ | (11) |
Here, we omitted the subscript mb for the membrane and will do so in the following. The drug at the donor side of, but outside, the membrane is given by cin. The inner-side drug mass fraction of the membrane at y = 0 is then, via the partition theory, given by
| c|y=0=ϑcin1+(ϑ−−1)cin | (12) |
At the cell-side of the membrane, y = 1, we assume a continuity of fluxes so that
| −D∂c∂y=−Dce∂cce∂y | (13) |
where the subscript ce denotes the inner cell. Assuming that the inner cell is an unstirred well-mixed environment, the drug uptake can be considered to occur under a constant rising concentration. Therefore, the gradient is nearly zero, so that we can say that ∂cce∂y≈0. This allows to state approximately that at y = 1
| ∂c∂y|y=1=0 | (14) |
The initial mass fraction within the membrane is given by
| c|t=0=c0 | (15) |
A supplementary interest is the continuous delivery of drugs out of a scaffold. In this case, a scaffold, containing the drugs, is introduced into the area of interest (the manner in which this is done is not the point of discussion here). This scaffold is a porous structure, where we are interested to provide a tool for analysing the diffusional flow through this medium. One side of the scaffold contains the drugs, maintained at a certain concentration or mass fraction, cin. The other side is where the drugs enter the area of interest, where we assume that the uptake is immediate so that the mass fraction is zero there, i.e. c0 = 0. Let us approximate such a scaffold as cylindrical pores. Contrary to the previous subsection, the dimension of the flux direction (axial direction of the pores) is larger than the radius of such a cylinder (radial direction of the pores). In that case, we can say (and also motivated by previous work in porous systems [23]) that the flux in the axial direction of the pores (the direction of drug release) will be approximately constant. However, the flux will be variable over the radial direction of the flux. Indeed, due to the small pore size with respect to the axial dimension of the pores, the wall effects will influence the flux and a flux profile exists in the radial direction. As for the mass fraction gradient (being a slow variable), its radial gradient is negligible with respect to the axial one (influenced by the drug release). As such, Eq (4) gives in dimensionless form (using ρDR, R and L as scaling for the mass flux, the radial spatial coordinate and the axial spatial coordinate, respectively, where R is the pore radius and L the length of the pore)
| J=−α∂c∂x+Kn2r1r∂∂r(r∂J∂r), | (16) |
with Knr≡ℓR, r the radial spatial coordinate, x the axial spatial coordinate in the pores and α≡RL. Assuming maximum flux in the middle of the pore (r = 0), i.e. ∂J∂r|r=0=0, and first-order slip (defining Bs≡ℓsR, where ℓs stands for the slip length, related to the wall material) at the wall (r = 1), i.e.
| J=−Bs∂J∂r|r=1, | (17) |
the flux can be found as a function of r. Taking subsequently an average value over a cross-section of the cylinder gives the following mass flux
| J=−α∂c∂x(1−2Kn2rB[1,1Knr]KnrB[0,1Knr]+BsB[1,1Knr]) | (18) |
with B the Bessel-I function. It should be noted that Eq (18) depends on how slip is observed and how the mass flux is defined. Other effects can influence the mass flux, but here it is the purpose to show how physical models can be used in the field of nanomedicine. Such an effective flux (taking into account the material of the scaffold (via the slip factor) and non-local effects via the Knudsen number) is now filled in the dimensionless (time scale is taken as L2D) species balance in the flow direction
| ∂tc=−1α∂xJ | (19) |
The evolution of the mass fraction is then described as
| ∂tc=∂2c∂x2(1−2Kn2B[1,1Kn]KnB[0,1Kn]+BsB[1,1Kn]) | (20) |
The boundary condition on one side of the scaffold is given by
| ∂c∂x|x=0=0 | (21) |
at x = 0, which stands for assuming that the drug only goes on one direction and
| c|x=1=c0e−t | (22) |
at x = 1, which allows for a drug uptake relaxation time of the same order of magnitude as the characteristic diffusion time (which is reasonable, considering that transport in cells is mainly of diffusive character), describing the mass fraction going from c0 to 0 (complete drug uptake). This depends of course of the system in question and the area in the body. Nevertheless, studying the effect of non-locality and slip factor on drug uptake, the choice of a relaxation time is not of importance here and does not influence our analysis. Eq (22) also shows the initial condition by evaluating it at t = 0.
The model is solved using a semi-implicit difference method for the spatial discretization and the forward Euler method for the temporal discretization for each time step t. At the beginning of the simulation, the initial conditions are imposed. At a certain time t, the equations for the mass flux and mass fractions are written in matrix form A*V = B, where A is a matrix of dimension (n-2)×(n-2) (n being the number of nodes in the discretized spatial domain and two nodes being dedicated to the boundary conditions, the bulk equations are thus discretized on (n-2) nodes) containing the coefficients in the equations for the mass fractions and mass fluxes at time t+1, V a vector of dimension (n-2) containing the unknown mass fraction and mass flux values at the spatial nodes at time t+1 and B a vector of dimension (n-2) containing the known mass fraction and mass flulx values at time t. Due to the spatial gradients in the bulk equations, the values at nodes 2 and n-1 depend on the boundary values at nodes 1 and n, respectively. The vector V is calculated as V = A-1B. The boundary condition for at the membrane entrance is imposed and held at that value at each time step. The boundary condition at the exit of the membrane at time t+1 is then obtained through the bulk values at time t+1. This procedure is repeated until the relative errors of the bulk and boundary values of both the mass fraction and mass flux between the previous and present loop become smaller than 10-4. The obtained values at time t+1 are stocked in the output matrix and used as the known values for the next time step t+2.
We use our model, given by Eqs. (10)–(12), (14)–(15), and solve it by means of a semi-implicite finite difference algorithm, explained in section 4.1. We take cin = 0.8, c0 = 0 and ϑ = 1 (we neglect the effect of the partition coefficient for the moment). The results are presented in Figures 1-5. Figures 1-3 show the evolution of the mass fraction across the membrane for Kn numbers Kn = 0.1, Kn = 1 and Kn = 10, respectively. The different Kn numbers stand for different sizes of the membrane with respect to the molecular mean free path of the drug molecules in the fluid present in the membrane (which can be taken as a watery solution). Furthermore, for each Kn number different fKn values are used: fKn = 0.1, fKn = 1 and fKn = 10.
We can see that for a given Kn number a higher fKn number (higher relaxation time with respect to a fixed diffusion characteristic time), the evolution of the mass flux is significantly slower, with a large gradient at the membrane entrance. This implies that the mass flux is very high at the membrane entrance but very low at the membrane exit. This effect is precisely referred to as non-locality. It is important here that by choosing the diffusion characteristic time as the scale time, any change in fKn implies a change in the relaxation time. This means that one compares systems with different relaxation times, which illustrates the effect of medication diffusion as a function of different medication-carrier fluid systems. A higher Kn number has the same non-local effect, albeit on the geometrical level (for the same system). This latter effect is not so clearly visible in Figures 1–3, but can be more appreciated looking to the evolution of the mass fraction at the membrane exit y = 1 in Figure 4.
Effectively, Figure 4 shows that a higher Kn number, for all fKn values, results into a lower mass fraction, which is representative for a lower drug diffusion/solubility in the membrane. The non-local effect at high Kn numbers is responsible for a non-Fickian mass fraction profile, which is indicated by a mass fraction at the exit of the membrane that is less than the initial mass fraction at steady state (this in comparison with Fickian behaviour, where for a zero-flux boundary condition, the profile should become flat). This means that for drugs that have large relaxation times or mean free paths, the drug uptake could be troublesome. It is interesting to note that the values shown in Figure 4 are exactly at the membrane exit, but that a mass fraction jump exists just before that exit. Not being well discernible in Figures 1-3, a zoom of some of those curves has been plotted in Figure 5.
This mass fraction jump can be seen clearly in Figure 5 and stems from the non-local behaviour from both the mass flux time derivative (∂tJ) and mass flux diffusion term (∇2J). Now, if the partition coefficient was not at unity, we would have not only a jump due to a non-locality at the membrane side of the interface, but also at the cell side, due the partition coefficient. Figure 6 shows the influence of the partition coefficient on the mass fraction evolution.
It can be seen that a higher partition coefficient results into a higher mass fraction and a higher drug uptake (higher mass fraction at the membrane exit). Since a zero-flux boundary condition (standing for a quick drug diffusion in the cell and not for a wall-like boundary, as is discussed under Eq (13)) is maintained at the membrane exit, the partition coefficient does not affect the mass fraction here. However, it does effect the quantitative value of the mass fraction just outside the membrane (in the cell). The higher the partition coefficient, the lower the mass fraction will be here. Also, as Figure 6 seems to show, a higher partition coefficient does not exhibit a different behaviour of the mass fraction at the membrane exit. That is, a higher partition coefficient gives a higher mass fraction of the drug at the membrane entrance, which translates itself into a higher mass fraction at the membrane exit. Nonetheless, this increase is cancelled by the partition coefficient itself at the membrane exit. This motivates using drugs that have a higher partition coefficient at the membrane entrance, but a lower one at the membrane exit. This would, of course, not lead to a higher total uptake (mass conservation implies that the total quantity of the uptake is equal to the administered quantity), but may lead to a higher efficiency if a higher mass fraction is needed for a better healing in a shorter time period. Having discussed drug uptake by means of drug diffusion through a cell membrane, it is also interesting to provide tools to drug release by scaffolds, applicable in the field of nanomedicine. This is treated in the next subsection.
We use our model, given by Eqs (20)–(22), and solve it in the same way as for the non-local drug diffusion model. We take c0 = 0.8 and present the results in Figures 7 and 8 for different Kn numbers and slip factors Bs. When we take Bs = 0, this stands for a zero slip length of the mass flux. This means that the mass flux at the wall of the pores in the scaffold is equal to zero. However, when Bs→∞, this means that there is full slip at the pore walls and that the mass flux is everywhere equal in the radial direction. We take a value Bs = 10 to be quite large, standing for "large slip" and compare its results to "no-slip". The slip behaviour of the mass flux is not clearly defined in the literature, but is rather often used in the framework of fluid flow, related to its velocity. The very definition of a mass flux in a porous medium being [23]
| J=ϕρf(vf−v) | (23) |
where ϕ is the porosity of the porous medium, ρf the density of the fluid flowing through it, vf the fluid velocity and v the barycentric velocity, we see that the mass flux is directly related to the fluid velocity. If we note that in a two-component system the barycentric velocity is expressed as ρv = ϕρfvf+(1-ϕ)ρsvs, wherein the total mass density ρ is ρ = ϕρf+(1-ϕ)ρs with ρs and vs the density and velocity of the solid component, respectively, and that in case of a scaffold vs = 0, we can quickly see that in (23) the only variable the mass flux is related directly to is the fluid velocity vf. As such, the discussion on the slip factor becomes more physical and understandable.
Figures 7 and 8 show that the Kn number tends to delay the drug release. This suggests that for pores of sizes of the order of magnitude of the mean free path of the drug molecules the drug release is slowed down. On the other hand, in some applications, it is necessary to have pores of such a size. Figure 8 shows that the slower drug release can be compensated by a higher slip factor, which accelerates the release of medication, illustrated by mass fraction profiles that evolve faster with respect to Figure 7. Figure 9 shows the mass fraction evolution for different Kn numbers for "no-slip" (Bs = 0) and "high-slip" (Bs = 10) at the membrane exit (y = 1). It can be seen that indeed, when Bs is increased, the reduction of the mass fraction due to an increase in the Kn number (from Kn = 0.1 to Kn = 1) is compensated, resulting into a similar profile. This confirms that by increasing the slip factor, which amounts to say altering the material properties of the scaffold pores to properties that reduce the interaction with the drug molecule solutions, compensates the reduction of drug release in nanopores.
This work presents new developments in Extended Non-Equilibrium Thermodynamics, proposing a model for studying drug uptake through cell membranes and drug release out of scaffolds. The theory is based on elevating fluxes to the same independent status of the usual variables. In this particular case, the mass flux has the same status as the mass fraction or concentration of the medication solution. This means that the mass flux follows its own evolution next to that of the mass fraction, being connected on an equal level via a species balance. By taking into account the relative importance of the relaxation time with respect to the diffusion characteristic time, and of the mean free path with respect to the system characteristic size, non-local effects are included in this way.
As for the diffusion through a cell membrane, it appeared that for sizes and relaxation times of the order of magnitude of or smaller than the mean free path (Kn = O(1)) and diffusion characteristic time (M = O(1)), respectively, leads to a lower drug uptake. This is due to the non-local effect, mentioned before. Non-Fickian profiles result from this, where it can be seen that at high non-locality a large mass fraction gradient is present close to the membrane entrance, but a nearly horizontal one near the exit. The latter contributes to the reduction of drug uptake. Furthermore, it has also been shown that at the membrane exit, a mass fraction jump (within the membrane) occurs that is not due to the boundary condition, but also due the non-locality effect. This is an effect that has not only been observed for heat transport [25], where it is already an established phenomenon, but also for mass transport [26]. The effect of the partition coefficient on the mass fraction is also studied, since they are typically introduced to account for the typical mass fraction jump (outside the membrane) observed in membrane diffusion studies, though still presenting experimental uncertainty to its precise value [26,27,28]. No particular effect has been observed other than that a higher partition coefficient causes a higher mass fraction jump at the membrane boundaries. If the medication could be chosen in such a way that the partition coefficient is high at the entrance and low at the exit, the efficiency of the medication could be increased if a higher concentration would lead to a better healing. Of course, mass conservation implies that this would not influence the total quantity of the uptake, but rather the speed of uptake, which could influence the kinetics of healing.
As for the drug release from the scaffold, the non-locality effect also showed to have a reducing impact on the release of drugs from the scaffold. It appeared that such a reduction could be countered by choosing the material of the scaffold in such a way that (besides biocompatibility, which is not a complicated task on its own) the interaction between the scaffold pore walls and the medication solution is weakened. This results into a higher slip factor enhancing drug release.
Future developments of this work include experimental demonstration, using the proposed model in this work. As quantitative agreement between experiments and models are difficult to find in the literature, comparison with experimental data should first be performed by nanomaterials of which material properties are well known and documented. As example, we can mention the delivery of nucleic acids with DNA nanostructures as new means of drug delivery by nanocarriers [29,30]. In the case of the nanocarriers, which can also be seen as scaffolds, the diffusion of the carriers through the cell membranes could be described by the set of Eqs (10)–(12) and (14)–(15). Subsequently, the delivery of the drugs out of the carriers into the cells could be described as release from a scaffold by Eqs (20)–(22). The difficulty lies in determining the interaction the medication can have with the nanocarrier, such as the DNA nanostructures [29]. In other words, it is important to know the slip length or slipping coefficient in Eq (17). Enhancement of drug release does not depend only on the release of the drug out of the carrier, but also the solubility of the carrier itself through the cell membranes. It is reported that DNA nanostructures can overcome the usually large circulation times of carriers [29,30]. The aforementioned discussion implies a necessity of measuring some material properties in order to use the presented model in this paper. In some cases, information about material properties could be obtained by in vitro experiments with a controlled experimental setup in a similar way performed by [31] for high-precision regulation of cell behaviour. This brings the intermediate step between modelling and full experimental demonstration to be identified by bio-mimicking cell membranes or scaffolds. Finally, such information can complete an experimental demonstration of our model.
Financial support from BelSPo is acknowledged.
The author declares there is no conflict of interest.
| [1] | Petroleum B, BP statistical review of world energy, 2015. Available from: www. bp. com/statisticalreview. |
| [2] | Sihag S, Pathak H, Jaroli D (2014) Factors affecting the rate of biodegradation of polyaromatic hydrocarbons. Int J Pure App Biosci 2: 185-202. |
| [3] |
Stroud JL, Paton GI, Semple KT (2007) Microbe-aliphatic hydrocarbon interactions in soil: implications for biodegradation and bioremediation. J Appl Microbiol 102: 1239-1253. doi: 10.1111/j.1365-2672.2007.03401.x
|
| [4] | Das N, Chandran P (2011) Microbial degradation of petroleum hydrocarbon contaminants: an overview. Biotechnol Res Int 2011: 941810. |
| [5] |
Gallego JL, Loredo J, Llamas JF, et al. (2001) Bioremediation of diesel-contaminated soils: evaluation of potential in situ techniques by study of bacterial degradation. Biodegradation 12: 325-335. doi: 10.1023/A:1014397732435
|
| [6] | Gong R, Hochmuth M, Weatherburn DC (2003) Polycyclic aromatic hydrocarbons and diesel particulates. |
| [7] |
Latif I, Karim A, Zuki A, et al. (2010) Pulmonary modulation of benzo [a] pyrene-induced hemato-and hepatotoxicity in broilers. Poultry Sci 89: 1379-1388. doi: 10.3382/ps.2009-00622
|
| [8] |
Samanta SK, Singh OV, Jain RK (2002) Polycyclic aromatic hydrocarbons: environmental pollution and bioremediation. Trends Biotechnol 20: 243-248. doi: 10.1016/S0167-7799(02)01943-1
|
| [9] |
Huang XD, El-Alawi Y, Penrose DM, et al. (2004) A multi-process phytoremediation system for removal of polycyclic aromatic hydrocarbons from contaminated soils. Environ Pollut 130: 465-476. doi: 10.1016/j.envpol.2003.09.031
|
| [10] |
Arthur EL, Rice PJ, Rice PJ, et al. (2005) Phytoremediation—An Overview. CRC Crit Rev Plant Sci 24: 109-122. doi: 10.1080/07352680590952496
|
| [11] |
Guo J, Feng R, Ding Y, et al. (2014) Applying carbon dioxide, plant growth-promoting rhizobacterium and EDTA can enhance the phytoremediation efficiency of ryegrass in a soil polluted with zinc, arsenic, cadmium and lead. J Environ Manage 141: 1-8. doi: 10.1016/j.jenvman.2013.12.039
|
| [12] | Onwurah I, Ogugua V, Onyike N, et al. (2007) Crude oil spills in the environment, effects and some innovative clean-up biotechnologies. |
| [13] |
Yemashova NA, Murygina VP, Zhukov DV, et al. (2007) Biodeterioration of crude oil and oil derived products: a review. Rev Environ Sci Biotechnol 6: 315-337. doi: 10.1007/s11157-006-9118-8
|
| [14] |
Head IM, Jones DM, Röling WF (2006) Marine microorganisms make a meal of oil. Nat Rev Microbiol 4: 173-182. doi: 10.1038/nrmicro1348
|
| [15] | Yasin G, Bhanger MI, Ansari TM, et al. (2013) Quality and chemistry of crude oils. J Pet Technol Altern Fuels 4: 53-63. |
| [16] | Alvarez PJ, Illman WA (2005) Bioremediation and natural attenuation: process fundamentals and mathematical models. John Wiley & Sons. |
| [17] |
Martinez-Gomez C, Vethaak A, Hylland K, et al. (2010) A guide to toxicity assessment and monitoring effects at lower levels of biological organization following marine oil spills in European waters. ICES J Mar Sci 67: 1105-1118. doi: 10.1093/icesjms/fsq017
|
| [18] |
Haritash A, Kaushik C (2009) Biodegradation aspects of polycyclic aromatic hydrocarbons (PAHs): a review. J Hazard Mater 169: 1-15. doi: 10.1016/j.jhazmat.2009.03.137
|
| [19] |
Mishra S, Jyot J, Kuhad RC, et al. (2001) Evaluation of inoculum addition to stimulate in situ bioremediation of oily-sludge-contaminated soil. Appl Environ Microb 67: 1675-1681. doi: 10.1128/AEM.67.4.1675-1681.2001
|
| [20] |
Kim KH, Jahan SA, Kabir E, et al. (2013) A review of airborne polycyclic aromatic hydrocarbons (PAHs) and their human health effects. Environ Int 60: 71-80. doi: 10.1016/j.envint.2013.07.019
|
| [21] | Masih J, Singhvi R, Kumar K, et al. (2012) Seasonal variation and sources of polycyclic aromatic hydrocarbons (PAHs) in indoor and outdoor air in a semi arid tract of northern India. AAQR 12: 515-525. |
| [22] | Perelo LW (2010) Review: in situ and bioremediation of organic pollutants in aquatic sediments. J Hazard Mater 177: 81-89. |
| [23] |
Gaspar A, Zellermann E, Lababidi S, et al. (2012) Characterization of saturates, aromatics, resins, and asphaltenes heavy crude oil fractions by atmospheric pressure laser ionization Fourier transform ion cyclotron resonance mass spectrometry. Energ Fuel 26: 3481-3487. doi: 10.1021/ef3001407
|
| [24] |
Leyva C, Ancheyta J, Berrueco C, et al. (2013) Chemical characterization of asphaltenes from various crude oils. Fuel Process Technol 106: 734-738. doi: 10.1016/j.fuproc.2012.10.009
|
| [25] |
Vinas M, Grifoll M, Sabate J, et al. (2002) Biodegradation of a crude oil by three microbial consortia of different origins and metabolic capabilities. J Ind Microbiol Biot 28: 252-260. doi: 10.1038/sj.jim.7000236
|
| [26] |
He L, Lin F, Li X, et al. (2015) Interfacial sciences in unconventional petroleum production: from fundamentals to applications. Chem Soc Rev 44: 5446-5494. doi: 10.1039/C5CS00102A
|
| [27] |
Fuentes S, Méndez V, Aguila P, et al. (2014) Bioremediation of petroleum hydrocarbons: catabolic genes, microbial communities, and applications. Appl Microbiol Biot 98: 4781-4794. doi: 10.1007/s00253-014-5684-9
|
| [28] |
Dong TT, Lee BK (2009) Characteristics, toxicity, and source apportionment of polycylic aromatic hydrocarbons (PAHs) in road dust of Ulsan, Korea. Chemosphere 74: 1245-1253. doi: 10.1016/j.chemosphere.2008.11.035
|
| [29] |
Ou S, Zheng J, Zheng J, et al. (2004) Petroleum hydrocarbons and polycyclic aromatic hydrocarbons in the surficial sediments of Xiamen Harbour and Yuan Dan Lake, China. Chemosphere 56: 107-112. doi: 10.1016/j.chemosphere.2004.02.022
|
| [30] |
Chandra S, Sharma R, Singh K, et al. (2013) Application of bioremediation technology in the environment contaminated with petroleum hydrocarbon. Ann Microbiol 63: 417-431. doi: 10.1007/s13213-012-0543-3
|
| [31] |
Paruk JD, Adams EM, Uher-Koch H, et al. (2016) Polycyclic aromatic hydrocarbons in blood related to lower body mass in common loons. Sci Total Environ 565: 360-368. doi: 10.1016/j.scitotenv.2016.04.150
|
| [32] |
Militon C, Boucher D, Vachelard C, et al. (2010) Bacterial community changes during bioremediation of aliphatic hydrocarbon-contaminated soil. Fems Microbiol Ecol 74: 669-681. doi: 10.1111/j.1574-6941.2010.00982.x
|
| [33] |
Wasmund K, Burns KA, Kurtböke DI, et al. (2009) Novel alkane hydroxylase gene (alkB) diversity in sediments associated with hydrocarbon seeps in the Timor Sea, Australia. Appl Environ Microb 75: 7391-7398. doi: 10.1128/AEM.01370-09
|
| [34] |
Adgate JL, Goldstein BD, McKenzie LM (2014) Potential public health hazards, exposures and health effects from unconventional natural gas development. Environ Sci Technol 48: 8307-8320. doi: 10.1021/es404621d
|
| [35] | Zhu X, Venosa AD, Suidan MT, et al. (2001) Guidelines for the bioremediation of marine shorelines and freshwater wetlands. US EPA |
| [36] | Sadler R, Connell D (2003) Analytical methods for the determination of total petroleum hydrocarbons in soil. In: Proceedings of the fifth national workshop on the assessment of site contamination, 133-150. |
| [37] | Van Beilen JB, Funhoff EG (2007) Alkane hydroxylases involved in microbial alkane degradation. Appl Microbiol Biot 74: 13-21. |
| [38] |
Hatzinger PB, Alexander M (1995) Effect of aging of chemicals in soil on their biodegradability and extractability. Environ Sci Technol 29: 537-545. doi: 10.1021/es00002a033
|
| [39] |
Semple KT, Morriss A, Paton G (2003) Bioavailability of hydrophobic organic contaminants in soils: fundamental concepts and techniques for analysis. Eur J Soil Sci 54: 809-818. doi: 10.1046/j.1351-0754.2003.0564.x
|
| [40] |
Baboshin M, Golovleva L (2012) Aerobic bacterial degradation of polycyclic aromatic hydrocarbons (PAHs) and its kinetic aspects. Microbiology 81: 639-650. doi: 10.1134/S0026261712060021
|
| [41] |
Rojo F (2009) Degradation of alkanes by bacteria. Environ Microbiol 11: 2477-2490. doi: 10.1111/j.1462-2920.2009.01948.x
|
| [42] | Wentzel A, Ellingsen TE, Kotlar HK, et al. (2007) Bacterial metabolism of long-chain n-alkanes. Appl Microbiol Biot 76: 1209-1221. |
| [43] | Kotani T, Yurimoto H, Kato N, et al. (2007) Novel acetone metabolism in a propane-utilizing bacterium, Gordonia sp. strain TY-5. J Bacteriol 189: 886-893. |
| [44] |
Bamforth SM, Singleton I (2005) Bioremediation of polycyclic aromatic hydrocarbons: current knowledge and future directions. J Chem Technol Biot 80: 723-736. doi: 10.1002/jctb.1276
|
| [45] |
Kanaly RA, Harayama S (2010) Advances in the field of high-molecular-weight polycyclic aromatic hydrocarbon biodegradation by bacteria. Microb Biotechnol 3: 136-164. doi: 10.1111/j.1751-7915.2009.00130.x
|
| [46] |
Seo JS, Keum YS, Li QX (2009) Bacterial degradation of aromatic compounds. Int J Environ Res Public Health 6: 278-309. doi: 10.3390/ijerph6010278
|
| [47] | Maigari AU, Maigari MU (2015) Microbial metabolism of polycyclic aromatic hydrocarbons (PAHs): a review. Int J Sci Eng 6: 1449-1459. |
| [48] | Rojo F, (2010) Enzymes for aerobic degradation of alkanes, In: Handbook of hydrocarbon and lipid microbiology, Springer, 781-797. |
| [49] |
Tortella GR, Diez MC, Durán N (2005) Fungal diversity and use in decomposition of environmental pollutants. Crit Rev Microbiol 31: 197-212. doi: 10.1080/10408410500304066
|
| [50] |
Brooijmans RJ, Pastink MI, Siezen RJ (2009) Hydrocarbon-degrading bacteria: the oil-spill clean-up crew. Microb Biotechnol 2: 587-594. doi: 10.1111/j.1751-7915.2009.00151.x
|
| [51] |
Germida J, Frick C, Farrell R, et al. (2002) Phytoremediation of oil-contaminated soils. Dev Soil Sci 28: 169-186. doi: 10.1016/S0166-2481(02)80015-0
|
| [52] | Wang YN, Cai H, Chi CQ, et al. (2007) Halomonas shengliensis sp. nov., a moderately halophilic, denitrifying, crude-oil-utilizing bacterium. Int J Syst Evol Micr 57: 1222-1226. |
| [53] | Mnif S, Chamkha M, Sayadi S (2009) Isolation and characterization of Halomonas sp. strain C2SS100, a hydrocarbon-degrading bacterium under hypersaline conditions. J Appl Microbiol 107: 785-794. |
| [54] |
Borzenkov I, Milekhina E, Gotoeva M, et al. (2006) The properties of hydrocarbon-oxidizing bacteria isolated from the oilfields of Tatarstan, Western Siberia, and Vietnam. Microbiology 75: 66-72. doi: 10.1134/S0026261706010127
|
| [55] | Dastgheib SMM, Amoozegar MA, Khajeh K, et al. (2011) A halotolerant Alcanivorax sp. strain with potential application in saline soil remediation. Appl Microbiol Biot 90: 305-312. |
| [56] |
Chaillan F, Le Flèche A, Bury E, et al. (2004) Identification and biodegradation potential of tropical aerobic hydrocarbon-degrading microorganisms. Res Microbiol 155: 587-595. doi: 10.1016/j.resmic.2004.04.006
|
| [57] |
Moody JD, Freeman JP, Fu PP, et al. (2004) Degradation of benzo [a] pyrene by Mycobacterium vanbaalenii PYR-1. Appl Environ Microb 70: 340-345. doi: 10.1128/AEM.70.1.340-345.2004
|
| [58] |
Van Beilen JB, Li Z, Duetz WA, et al. (2003) Diversity of alkane hydroxylase systems in the environment. Oil Gas Sci Technol 58: 427-440. doi: 10.2516/ogst:2003026
|
| [59] |
McKew BA, Coulon F, Osborn AM, et al. (2007) Determining the identity and roles of oil-metabolizing marine bacteria from the Thames estuary, UK. Environ Microbiol 9: 165-176. doi: 10.1111/j.1462-2920.2006.01125.x
|
| [60] |
Pilon-Smits E (2005) Phytoremediation. Annu Rev Plant Biol 56: 15-39. doi: 10.1146/annurev.arplant.56.032604.144214
|
| [61] |
Holden P, LaMontagne M, Bruce A, et al. (2002) Assessing the role of Pseudomonas aeruginosa surface-active gene expression in hexadecane biodegradation in sand. Appl Environ Microb 68: 2509-2518. doi: 10.1128/AEM.68.5.2509-2518.2002
|
| [62] |
Makadia TH, Adetutu EM, Simons KL, et al. (2011) Re-use of remediated soils for the bioremediation of waste oil sludge. J Environ Manage 92: 866-871. doi: 10.1016/j.jenvman.2010.10.059
|
| [63] |
Sheppard PJ, Simons KL, Kadali KK, et al. (2012) The importance of weathered crude oil as a source of hydrocarbonoclastic microorganisms in contaminated seawater. J Microbiol Biotechnol 22: 1185-1192. doi: 10.4014/jmb.1201.01049
|
| [64] | Rodríguez-Blanco A, Antoine V, Pelletier E, et al. (2010) Effects of temperature and fertilization on total vs. active bacterial communities exposed to crude and diesel oil pollution in NW Mediterranean Sea. Environ Pollut 158: 663-673. |
| [65] |
Sarkar D, Ferguson M, Datta R, et al. (2005) Bioremediation of petroleum hydrocarbons in contaminated soils: comparison of biosolids addition, carbon supplementation, and monitored natural attenuation. Environ Pollut 136: 187-195. doi: 10.1016/j.envpol.2004.09.025
|
| [66] |
Tang J, Wang R, Niu X, et al. (2010) Enhancement of soil petroleum remediation by using a combination of ryegrass (Lolium perenne) and different microorganisms. Soil Till Res 110: 87-93. doi: 10.1016/j.still.2010.06.010
|
| [67] |
Li X, Wu Y, Lin X, et al. (2012) Dissipation of polycyclic aromatic hydrocarbons (PAHs) in soil microcosms amended with mushroom cultivation substrate. Soil Biol Biochem 47: 191-197. doi: 10.1016/j.soilbio.2012.01.001
|
| [68] | Agnello AC, Bagard M, Van Hullebusch ED, et al. (2016) Comparative bioremediation of heavy metals and petroleum hydrocarbons co-contaminated soil by natural attenuation, phytoremediation, bioaugmentation and bioaugmentation-assisted phytoremediation. Sci Total Environ 563: 693-703. |
| [69] |
Kadali KK, Simons KL, Sheppard PJ, et al. (2012) Mineralisation of weathered crude oil by a hydrocarbonoclastic consortia in marine mesocosms. Water Air Soil Poll 223: 4283-4295. doi: 10.1007/s11270-012-1191-8
|
| [70] |
Yu K, Wong A, Yau K, et al. (2005) Natural attenuation, biostimulation and bioaugmentation on biodegradation of polycyclic aromatic hydrocarbons (PAHs) in mangrove sediments. Mar Pollut Bull 51: 1071-1077. doi: 10.1016/j.marpolbul.2005.06.006
|
| [71] |
Margesin R, Hämmerle M, Tscherko D (2007) Microbial activity and community composition during bioremediation of diesel-oil-contaminated soil: effects of hydrocarbon concentration, fertilizers, and incubation time. Microbial Ecol 53: 259-269. doi: 10.1007/s00248-006-9136-7
|
| [72] |
Festa S, Coppotelli B, Morelli I (2016) Comparative bioaugmentation with a consortium and a single strain in a phenanthrene-contaminated soil: Impact on the bacterial community and biodegradation. Appl Soil Ecol 98: 8-19. doi: 10.1016/j.apsoil.2015.08.025
|
| [73] | Nikolopoulou M, Kalogerakis N, (2010) Biostimulation strategies for enhanced bioremediation of marine oil spills including chronic pollution. In: Handbook of hydrocarbon and lipid microbiology, Springer, 2521-2529. |
| [74] | Alexander M (1999) Biodegradation and bioremediation. Gulf Professional Publishing. |
| [75] |
Walworth J, Pond A, Snape I, et al. (2007) Nitrogen requirements for maximizing petroleum bioremediation in a sub-Antarctic soil. Cold Reg Sci Technol 48: 84-91. doi: 10.1016/j.coldregions.2006.07.001
|
| [76] |
Piehler MF, Swistak J, Pinckney J, et al. (1999) Stimulation of diesel fuel biodegradation by indigenous nitrogen fixing bacterial consortia. Microbial Ecol 38: 69-78. doi: 10.1007/s002489900157
|
| [77] |
Li H, Zhao Q, Boufadel MC, et al. (2007) A universal nutrient application strategy for the bioremediation of oil-polluted beaches. Mar Pollut Bull 54: 1146-1161. doi: 10.1016/j.marpolbul.2007.04.015
|
| [78] |
Chaillan F, Chaineau C, Point V, et al. (2006) Factors inhibiting bioremediation of soil contaminated with weathered oils and drill cuttings. Environ Pollut 144: 255-265. doi: 10.1016/j.envpol.2005.12.016
|
| [79] |
Sayara T, Borràs E, Caminal G, et al. (2011) Bioremediation of PAHs-contaminated soil through composting: Influence of bioaugmentation and biostimulation on contaminant biodegradation. Int Biodeter Biodegr 65: 859-865. doi: 10.1016/j.ibiod.2011.05.006
|
| [80] |
Ros M, Rodriguez I, Garcia C, et al. (2010) Microbial communities involved in the bioremediation of an aged recalcitrant hydrocarbon polluted soil by using organic amendments. Bioresour Technol 101: 6916-6923. doi: 10.1016/j.biortech.2010.03.126
|
| [81] |
Blyth W, Shahsavari E, Morrison PD, et al. (2015) Biosurfactant from red ash trees enhances the bioremediation of PAH contaminated soil at a former gasworks site. J Environ Manage 162: 30-36. doi: 10.1016/j.jenvman.2015.07.041
|
| [82] |
Mair J, Schinner F, Margesin R (2013) A feasibility study on the bioremediation of hydrocarbon-contaminated soil from an alpine former military site: effects of temperature and biostimulation. Cold Reg Sci Technol 96: 122-128. doi: 10.1016/j.coldregions.2013.07.006
|
| [83] |
Liu W, Luo Y, Teng Y, et al. (2010) Bioremediation of oily sludge-contaminated soil by stimulating indigenous microbes. Environ Gecochem Hlth 32: 23-29. doi: 10.1007/s10653-009-9262-5
|
| [84] |
Palmroth MR, Pichtel J, Puhakka JA (2002) Phytoremediation of subarctic soil contaminated with diesel fuel. Bioresour Technol 84: 221-228. doi: 10.1016/S0960-8524(02)00055-X
|
| [85] |
Schaefer M, Juliane F (2007) The influence of earthworms and organic additives on the biodegradation of oil contaminated soil. Appl Soil Ecol 36: 53-62. doi: 10.1016/j.apsoil.2006.11.002
|
| [86] |
Suja F, Rahim F, Taha MR, et al. (2014) Effects of local microbial bioaugmentation and biostimulation on the bioremediation of total petroleum hydrocarbons (TPH) in crude oil contaminated soil based on laboratory and field observations. Int Biodeter Biodegr 90: 115-122. doi: 10.1016/j.ibiod.2014.03.006
|
| [87] | Al-Sulaimani HS, Al-Wahaibi YM, Al-Bahry S, et al. (2010) Experimental investigation of biosurfactants produced by Bacillus species and their potential for MEOR in Omani oil field. Society of Petroleum Engineers. |
| [88] |
Reid BJ, Jones KC, Semple KT (2000) Bioavailability of persistent organic pollutants in soils and sediments—a perspective on mechanisms, consequences and assessment. Environ Pollut 108: 103-112. doi: 10.1016/S0269-7491(99)00206-7
|
| [89] |
Liu PWG, Chang TC, Chen CH, et al. (2014) Bioaugmentation efficiency investigation on soil organic matters and microbial community shift of diesel-contaminated soils. Int Biodeter Biodegr 95: 276-284. doi: 10.1016/j.ibiod.2014.05.004
|
| [90] |
Hyun S, Ahn MY, Zimmerman AR, et al. (2008) Implication of hydraulic properties of bioremediated diesel-contaminated soil. Chemosphere 71: 1646-1653. doi: 10.1016/j.chemosphere.2008.01.026
|
| [91] |
Nam K, Chung N, Alexander M (1998) Relationship between organic matter content of soil and the sequestration of phenanthrene. Environ Sci Technol 32: 3785-3788. doi: 10.1021/es980428m
|
| [92] | Alexander M (2000) Aging, bioavailability, and overestimation of risk from environmental pollutants. Environ Sci Technol 34: 4259-4265. |
| [93] |
Si-Zhong Y, Hui-Jun J, Zhi W, et al. (2009) Bioremediation of oil spills in cold environments: a review. Pedosphere 19: 371-381. doi: 10.1016/S1002-0160(09)60128-4
|
| [94] | Okoh AI (2006) Biodegradation alternative in the cleanup of petroleum hydrocarbon pollutants. Biotechnol Mol Biol Rev 1: 38-50. |
| [95] | Margesin R, Schinner F (2001) Biodegradation and bioremediation of hydrocarbons in extreme environments. Appl Microbiol Biot 56: 650-663. |
| [96] |
Perfumo A, Banat IM, Marchant R, et al. (2007) Thermally enhanced approaches for bioremediation of hydrocarbon-contaminated soils. Chemosphere 66: 179-184. doi: 10.1016/j.chemosphere.2006.05.006
|
| [97] |
Scherr K, Aichberger H, Braun R, et al. (2007) Influence of soil fractions on microbial degradation behavior of mineral hydrocarbons. Eur J Soil Biol 43: 341-350. doi: 10.1016/j.ejsobi.2007.03.009
|
| [98] |
Amellal N, Portal J-M, Berthelin J (2001) Effect of soil structure on the bioavailability of polycyclic aromatic hydrocarbons within aggregates of a contaminated soil. Appl Geochem 16: 1611-1619. doi: 10.1016/S0883-2927(01)00034-8
|
| [99] |
Chaerun S, Tazaki K (2005) How kaolinite plays an essential role in remediating oil-polluted seawater. Clay Miner 40: 481-491. doi: 10.1180/0009855054040185
|
| [100] | Atlas RM, Philp J (2005) Bioremediation. Applied microbial solutions for real-world environmental cleanup. ASM Press. |
| [101] |
Zucchi M, Angiolini L, Borin S, et al. (2003) Response of bacterial community during bioremediation of an oil-polluted soil. J Appl Microbiol 94: 248-257. doi: 10.1046/j.1365-2672.2003.01826.x
|
| [102] | Tang J, Wang M, Wang F, et al. (2011) Eco-toxicity of petroleum hydrocarbon contaminated soil. J Environ Sci 23: 845-851. |
| [103] |
Sheppard PJ, Adetutu EM, Makadia TH, et al. (2011) Microbial community and ecotoxicity analysis of bioremediated, weathered hydrocarbon-contaminated soil. Soil Research 49: 261-269. doi: 10.1071/SR10159
|
| 1. | Lu Xu, Chunlai Mu, Qiao Xin, Global boundedness and asymptotic behavior of solutions for a quasilinear chemotaxis model of multiple sclerosis with nonlinear signal secretion, 2023, 28, 1531-3492, 1215, 10.3934/dcdsb.2022118 |
Figures(6) / Tables(4)
Eman Koshlaf, Andrew S Ball. Soil bioremediation approaches for petroleum hydrocarbon polluted environments[J]. AIMS Microbiology, 2017, 3(1): 25-49. doi: 10.3934/microbiol.2017.1.25
DownLoad:
