Next Article in Journal
Ultrasonically Fabricated Beta-Carotene Nanoemulsion: Optimization, Characterization and Evaluation of Combinatorial Effect with Quercetin on Streptozotocin-Induced Diabetic Rat Model
Next Article in Special Issue
Tailoring the Lamellarity of Liposomes Prepared by Dual Centrifugation
Previous Article in Journal
Delivery Systems for Mitochondrial Gene Therapy: A Review
Previous Article in Special Issue
Preliminary Assessment of Intramuscular Depot of Lipid-Based Decoquinate Formulation for Long-Term Chemoprophylaxis of Malaria
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interplay of Hydropathy and Heterogeneous Diffusion in the Molecular Transport within Lamellar Lipid Mesophases

by
Antonio M. Bosch
1 and
Salvatore Assenza
1,2,3,*
1
Departamento de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, 28049 Madrid, Spain
2
Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, 28049 Madrid, Spain
3
Instituto Nicolás Cabrera, Universidad Autónoma de Madrid, 28049 Madrid, Spain
*
Author to whom correspondence should be addressed.
Pharmaceutics 2023, 15(2), 573; https://doi.org/10.3390/pharmaceutics15020573
Submission received: 14 December 2022 / Revised: 18 January 2023 / Accepted: 3 February 2023 / Published: 8 February 2023

Abstract

:
Lipid mesophases are being intensively studied as potential candidates for drug-delivery purposes. Extensive experimental characterization has unveiled a wide palette of release features depending on the nature of the host lipids and of the guest molecule, as well as on the environmental conditions. However, only a few simulation works have addressed the matter, which hampers a solid rationalization of the richness of outcomes observed in experiments. Particularly, to date, there are no theoretical works addressing the impact of hydropathy on the transport of a molecule within lipid mesophases, despite the significant fraction of hydrophobic molecules among currently-available drugs. Similarly, the high heterogeneity of water mobility in the nanoscopic channels within lipid mesophases has also been neglected. To fill this gap, we introduce here a minimal model to account for these features in a lamellar geometry, and systematically study the role played by hydropathy and water–mobility heterogeneity by Brownian-dynamics simulations. We unveil a fine interplay between the presence of free-energy barriers, the affinity of the drug for the lipids, and the reduced mobility of water in determining the net molecular transport. More in general, our work is an instance of how multiscale simulations can be fruitfully employed to assist experiments in release systems based on lipid mesophases.

1. Introduction

Lipid mesophases are self-organized structures where nanoscopic solvent channels emerge from the self-arrangement of lipids in the host solvent. Despite often maintaining the fluidity of the membrane, these aggregates show spatial periodicity following standard crystallographic space groups. In previous studies, the geometry and topology of lipid mesophases have been shown to depend on various parameters, such as pH, temperature, and lipid–solvent (usually water) ratio [1,2,3,4,5,6,7,8]. This has recently made these systems of significant interest for biotechnological applications in material design and drug delivery, as well as fundamental research on ion pumps, membrane protein crystallization, and cryo-enzymatic reactions [9,10,11,12,13,14,15].
In particular, heterogeneity, reproducibility, and high biocompatibility of lipid mesophases have made them a potential tool for drug and nutraceutical delivery [14]. Moreover, the rich structural landscape of lipid mesophases has recently been shown to be naturally explored during digestion of tryglycerides, which likely impacts the delivery of drugs embedded in such hosts [16]. As a consequence, a large amount of research has been devoted to understanding how the chemical and structural features of lipid mesophases influence their transport properties [5,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33]. Experimental efforts in characterizing the molecular transport through lipid mesophases have unveiled an extensive palette of delivery performance depending on factors such as the geometry and symmetry of the lipid mesophase, the size of the diffusion domain, the hydropathy of the diffusing molecule, and lipid–solvent composition. While opening many possibilities for devising tailored applications, this rich heterogeneity of relevant features makes it lengthy and costly to fully characterize the role played by each of them via experimental assays. Theoretical considerations and computational simulations can provide a complementary insight into these systems, due to their faster and cheaper implementation, but also because they allow to address the effect of the various factors one at a time, which is not always possible in experiments. Although there is a vast literature focused on the modelling of diffusion in confined environments [34,35] and on the exotic water mobility at the interface with hydrophilic objects [36,37,38,39], little work has studied these topics within the context of lipid mesophases [8,28,29,40].
Strikingly, there is a lack of theoretical studies addressing the influence of molecular affinity for the lipids on the transport of the guest particle. This is a particularly important gap in view of the potential applications of lipid mesophases for drug delivery, since about one third of current drugs show low solubility in water [41,42]. The importance of this aspect in drug delivery is being increasingly acknowledged, for instance by investigations focused on solubility aspects in biorelevant media [43]. In this work, we take a step in understanding the impact of hydropathy by studying the diffusion kinetics of a particle spending a finite amount of time in both the water channels and the lipid bilayers. We analyze the role played by the complex interaction free energy between the guest particle and the lipid molecules, as well as the impact of the reduced mobility of water in the vicinity of lipid heads. The latter is expected to strongly affect the quantitative determination of transport properties, as water with lower mobility extends beyond roughly 1 nm starting from the lipid heads [8,40,44], which is comparable to the overall size of the water nanochannels inside lipid mesophases. We focus on the case of a lamellar arrangement with geometrical parameters obtained from reported experimental data. We find that diffusion in the direction perpendicular to the lipid/water interface is strongly regulated by the free-energy barriers obstaculating particle exchange between the lipid and water phases, while parallel diffusion is determined by the hydropathy of the guest molecule, as quantified by the partition coefficient. Heterogeneity of water mobility enters the picture in quantitatively regulating the large-scale diffusion coefficient.

2. Materials and Methods

In a lamellar geometry, the system can be conveniently expressed by considering a reference frame with the z-axis oriented along the normal to the lipid/water interface. The diffusion of a guest molecule can be assessed by considering a Brownian dynamics in a periodic free-energy landscape U ( z ) and with a space-dependent diffusion coefficient D ( z ) . Within this framework, the physico-chemical properties of the particle, as well as its interaction with the solvent and the lipids, are concealed within the profiles U ( z ) and D ( z ) . Ideally, these functions are obtained by running dedicated atomistic simulations of the system; an example is provided by the amino acids study presented later in the manuscript, where U ( z ) is implemented according to previously-published results of atomistic simulations. Here, we also consider a physically-sound toy model, which accounts, in a minimalistic fashion, for the main features of a typical system, as detailed below. Although simplistic, this framework enables accounting for the features of the system which are the focus of the present study; for example, in the context of strongly-hydrophilic particles in cubic phases, this approach has enabled estimation of the amount of bound water by comparison of theoretical and experimental diffusion coefficients, obtaining values in agreement with direct experimental observations [28]. However, one always has to keep in mind that various microscopical features are not accounted for, including for instance the impact of electrostatics in the diffusion of charged molecules, or the presence of perturbations in the bilayer such as Helfrich undulations [45].

2.1. Model

The hydropathy of the guest molecule is described by considering a periodic potential of mean force U ( z ) such as the one depicted in Figure 1b (we show only one periodicity), which corresponds to the equation
U ( z ) = Δ U if z l h Δ U + 8 Δ U b Δ U h 2 l h z 2 if   l h < z l 3 4 h Δ U b 8 Δ U b Δ U h 2 l 1 2 h z 2 if   l 3 4 h < z l 1 2 h Δ U b 8 Δ U b h 2 l 1 2 h z 2 if   l 1 2 h < z l 1 4 h 8 Δ U b h 2 l z 2 if   l 1 4 h < z l 0 if   l < z a 2 .
In the previous formula, a is the lattice parameter, l is the total length of a lipid within the bilayer, and h is the size of the lipid head (see Figure 1). In order to set reasonable values for these parameters, we considered a = 6.65 nm, which corresponds to fully-hydrated L α lamellar mesophases obtained by water/dipalmitoylphosphatidylcholine mixtures at 43 °C [46]. Moreover, we also set l = 2.365 nm and h = 1 nm based on the electron-density profile computed in Ref. [47] for the same mixture. As for the energy parameters, Δ U is the free-energy difference between the plateaus corresponding to the lipid tails and the water region. Therefore, Δ U > 0 for hydrophilic molecules (such as in Figure 1b), while Δ U < 0 in the case of hydrophobicity. The parameter Δ U b introduces a barrier in correspondence with the lipid heads, which can mimic the kinetic barriers associated with the permeability of the membrane; moreover, Δ U b can be employed to introduce depletion of molecules from the lipid heads ( Δ U b > 0 ) or the tendency to sit at the water/lipid interface typical of amphiphilic molecules, which is the case for many proteins ( Δ U b < 0 ) [48]. The use of parabolic fragments in Equation (1) allows for tuning the potential between 0 ,   Δ U and Δ U b , while ensuring that both U ( z ) and its derivative are continuous throughout space (see Figure 1b), which avoids undesirable numerical instabilities in the simulations.
A similar approach was used to account for heterogeneity in molecular transport. To this aim, we introduced a space-dependent diffusion coefficient D ( z ) (Figure 1c):
D ( z ) = D l i p if z l D l i p + 2 D w a t D l i p w 2 l z 2 if   l < z l + 1 2 w D w a t 2 D w a t D l i p w 2 l + w z 2 if   l + 1 2 w < z l + w D w a t if   l + w < z 1 2 a .
In the previous formula, D w a t and D l i p correspond to the diffusion coefficients of the guest molecule when considered in pure water and in the lipid bilayer, respectively, while w is the thickness of the water layer in which the continuous change between D l i p and D w a t takes place; therefore, w accounts for the reduced mobility of water molecules in the vicinity of the lipid heads [44].
Typical values of D w a t for nanoscopic objects are found in the range 10 10 10 9 m2/s. For instance, at 25 °C, one has D w a t = 0.7 0.9 × 10 9 m2/s for amino acids [49,50,51,52,53], and D w a t = 0.5 , 0.7 , 0.8 × 10 9 m2/s for ibuprofen [54], aspirin [55], and paracetamol [56], respectively. The value of D w a t is expected to be dependent on temperature, T. When small temperature differences are considered (such as estimation of D w a t at physiological temperature starting from room-temperature measurements), a simple yet effective approach to estimate the effect of T is to assume a Stokes–Einstein relation D w a t = k B T / ( 6 π η ( T ) R ) , where k B is Boltzmann’s constant, R is the size of the particle, and η ( T ) is the temperature-dependent viscosity of water. This approach has enabled accurate predictions of transport of glucose molecules in monolinolein-based cubic phases [57]. Unless stated otherwise, in our simulations, we consider D w a t = 0.7 × 10 9 m2/s.
As for the diffusion coefficient in the lipid phase, D l i p , one expects its value to be significantly smaller than D w a t due to the lower fluidity of the lipid membrane as compared to water. For instance, the three-dimensional self-diffusion of lipids for various monoacylglycerols with cubic symmetry has been reported to be 1.1 1.3 × 10 11 m2/s [26], which gives values in the range 1.7 2 × 10 11 m2/s for the lateral diffusion coefficient when accounting for the geometric constraint imposed by the minimal surface at the mid-plane of the lipid bilayer [58]. Amino acids and drugs such as the ones mentioned above are smaller than lipid molecules, so that D l i p is expected to be somewhat larger for them. Here, we fix D l i p = 0.09 D w a t , based on molecular dynamics simulations of paracetamol in DPPC [47].
Finally, the parameter w was set in accordance with experimental evidence and molecular dynamics simulations, which point to the existence of 3–4 layers of water with reduced mobility in proximity of the lipid heads [8,40,44]. The specific value of this thickness was selected to be w = 0.96 nm (Figure 1c), in order to ensure that D w a t is reached exactly at z = a / 2 , thus avoiding a discontinuity in the derivative of D ( z ) , which would have occurred for larger values of w.

2.2. Brownian Dynamics

Based on the model introduced above, we run Brownian Dynamics simulations to study the diffusion kinetics. Despite the one-dimensional nature of the potential U ( z ) and the diffusion coefficient D ( z ) , we integrated the motion in three dimensions, thus also including the movement along planes parallel to the lipid/water interface. This was motivated by the expected effect of a space-dependent diffusion coefficient also on lateral diffusion. Following standard Euler integration, and including a drift term to correctly implement the spatial dependence of the diffusion coefficient [59,60], the updating rule for the position r x , y , z of a particle is:
x ( t + d t ) = x ( t ) + 2 D d t ξ x , y ( t + d t ) = y ( t ) + 2 D d t ξ y , z ( t + d t ) = z ( t ) D k B T d U d z d t + d D d z d t + 2 D d t ξ z .
In the previous formula, t is time and d t is the integration timestep, while ξ x , ξ y and ξ z are random variables distributed according to a Gaussian function with zero average and unit variance. Although molecular dynamics simulations indicate that diffusion in parallel and perpendicular directions to the lipid/water interface are distinct [8,39], for simplicity in Equation (3), we consider, for a given position, the same diffusion coefficient for the random movement in any direction. More accurate quantitative estimations will require a proper account of this feature for the system under study.
The choice of length units ( σ = 1 nm) together with the value chosen for D w a t determines a “natural” simulation timescale τ = σ 2 / D w a t for the time t. For instance, if D w a t = 0.7 × 10 9 m2/s = 0.7 nm2/ns, one has τ 1.4 ns. The value of τ is important in the determination of the integration of the timestep d t . Our rationale in its choice was to consider the largest possible value of d t which correctly recovers the equilibrium distribution of a collection of particles. In this regard, for the system corresponding to U ( z ) and D ( z ) as in Figure 1, we performed simulations of 10 4 particles initially located at z = 0 for a total time 2 × 10 4 τ , by considering several values of d t ; the optimal choice turned out to be d t = 3 × 10 4 τ , which at the end of the simulation led to the correct sampling of the implemented U ( z ) (orange points Figure 1b).
For each study reported in Section 3, we performed simulations for ensembles of 10 3 particles up to times ranging between 5 × 10 3 τ and 5 × 10 5 τ . Particles were pre-equilibrated by randomly extracting their initial position from their Boltzmann distribution by using the Ziggurat algorithm [61]. The value of the total simulation time was adapted according to the system under study, to ensure that the mean-square displacement in the perpendicular direction reached at least 100 nm2. We found that this threshold value was sufficient to ensure a good estimation of the long-time diffusion coefficient.
The diffusion kinetics of the system was monitored by computing the mean-square displacement:
M S D ( t ) = x ( t ) x 0 2 + y ( t ) y 0 2 M S D ( t ) = z ( t ) z 0 2 .
In the previous formula, r 0 { x 0 , y 0 , z 0 } is the initial position of each particle, while denotes averaging over the whole ensemble. M S D ( t ) and M S D ( t ) describe the diffusion kinetics in the directions parallel and perpendicular to the lipid/water interface, respectively. At long times, in both cases, the MSD is expected to be purely diffusive [28]: M S D ( t ) = 4 D t and M S D ( t ) = 2 D t , where D and D are the effective diffusion coefficients. The different prefactors account for the different dimensionalities of the two diffusion kinetics (two and one dimensions, respectively).

2.3. Relationship between log P and D

The logarithm of the partition coefficient log P is defined as log P = log 10 ( c l i p / c w a t ) , where c l i p and c w a t are the concentrations of particles in the lipid and water phases [62]. Denoting as l the length of the lipids and as a the lattice parameter, one has in general
10 log P = c l i p c w a t = 0 l e U ( z ) k B T d z / l l a 2 e U ( z ) k B T d z / a 2 l
where we consider only half repeating units due to the symmetry of the system. Note that this formula is not restricted to our choice of U ( z ) operated above; instead, it holds, in general, under the approximation that one can sharply distinguish the lipid and water phases. As a further approximation, we assume that within the water phase U ( z ) = 0 . This neglects desolvation effects on the guest molecule within the water layer in contact with the lipid heads, as well as electrostatic interactions. Nevertheless, actual computation of U ( z ) from atomistic simulations suggests that this assumption is pretty reasonable (see, e.g., Figure 7a). Within this approximation, the denominator in the previous formula is 1, thus yielding
0 l e U ( z ) k B T d z = l 10 log P .
The average diffusion coefficient D is obtained by thermal averaging of D ( z ) :
D = 0 a 2 D ( z ) e U ( z ) k B T d z 0 a 2 e U ( z ) k B T d z .
The denominator in the previous formula can be rearranged by means of Equation (6) as
0 a 2 e U ( z ) k B T d z = 0 l e U ( z ) k B T d z + l a 2 e U ( z ) k B T d z = l 10 log P + a 2 l .
Similarly, the numerator in Equation (7) can be rearranged as
0 a 2 D ( z ) e U ( z ) k B T d z = 0 l D ( z ) e U ( z ) k B T d z + l l + w D ( z ) e U ( z ) k B T d z + l + w a 2 D ( z ) e U ( z ) k B T d z .
In the first and third term on the right-hand side of the previous formula, one can assume constant values for the diffusion coefficient equal to D l i p and D w a t , respectively. The second term can be rewritten as w D ¯ , where D ¯ is the average value of D ( z ) in the layer of water molecules with non-trivial mobility (compare Figure 1c). We assume that D ¯ = ( D w a t + D l i p ) / 2 , which considers a symmetric profile of D ( z ) within this region, as is the case for the toy model introduced in Equation (2). This allows for rewriting the numerator as
0 a 2 D ( z ) e U ( z ) k B T d z = D l i p 0 l e U ( z ) k B T d z + w D w a t + D l i p 2 + a 2 l w D w a t
that is, by rearranging and making use of Equation (6),
0 a 2 D ( z ) e U ( z ) k B T d z = D l i p l 10 log P + w 2 + D w a t a 2 l w 2 .
Plugging Equations (8) and (11) in Equation (7) finally yields the formula reported in Equation (14) in the main text.

2.4. Amino Acids Simulations

We simulated the diffusion of 16 amino acids by considering for U ( z ) the potential of mean force derived in Ref. [63], where these residues were considered in the presence of DOPC. For the four charged amino acids (Arg, Lys, Glu, Asp), we considered U ( z ) = m i n { U c ( z ) , U 0 ( z ) } , where U c ( z ) and U 0 ( z ) are the free-energy profiles obtained for the charged and neutral variants of the amino acids, shifted according to the free-energy cost of their neutralization in pure water [63]. We deem this approximation to be reasonable, as seen explicitly in constant-pH simulations performed for similar systems [64]. The U ( z ) profiles are reported in Figure 7a. As for the diffusion profile, we considered Equation (2) with parameters adapted to the present system. We kept the value a = 6.65 nm for the lattice parameter. Following Ref. [63], we assigned to the lipid length the value l = 2.5 nm. Accordingly, we set w = 0.825 nm in order to attain D w a t exactly at z = a / 2 . The diffusion coefficients in water, D w a t , were obtained by considering for each residue the values obtained experimentally [49,50,51,52,53,65,66] and subsequently estimating the corresponding values at 37 °C according to the Stokes–Einstein equation, as discussed above [57]. The obtained values range between 0.9 nm2/ns (obtained for Trp) and 1.3 nm2/ns (obtained for Cys), with an average equal to 1.04 nm2/ns. As in the toy model, the value of D l i p was set to D l i p = 0.09 D w a t for each residue.
Brownian dynamics simulations were performed with the same parameters as discussed above. Due to the large values of the free-energy barriers, the diffusion coefficients for Arg, Asn, Lys, Asp, Gln, and Ile could not be assessed by brute-force simulations. Instead, for each of these systems, we considered various sets of simulations in which U ( z ) was renormalized by a factor α < 1 . The value of D ( α ) was computed for each simulation set and its dependence on α was fitted via an exponential decay. The sought value for the original system was then obtained by extrapolating to α = 1 . For Asn, Lys, Asp, Gln, and Ile, we considered α = 0.3 , 0.4 , 0.5 , 0.6 . In the case of Arg, the strong barrier imposed lower values α = 0.15 , 0.20 , 0.25 , 0.30 . As for D , for each value of α , we checked the quantitative agreement between the value obtained from the simulations and the result found by applying Equation (13); then, we considered the predicted value computed for α = 1 .

3. Results

Periodic lipidic mesophases exist in a wide variety of arrangements, including lamellar, hexagonal, and cubic symmetries [14]. Moreover, for each symmetry, an extended range of geometrical features can be obtained, e.g., channel swelling by addition of co-surfactants [5]. For the lamellar symmetry, the lipids in the bilayer can be arranged in different ways, such as crystals ( L c phase), gels ( L β ), or fluid membranes ( L α ) [1,46]. As the main goal of this work is to understand the combined effect of heterogeneous diffusion and hydropathy on molecular transport, we focused on the simple case of a lamellar symmetry with fixed geometrical features (Figure 1), leaving the important topic of the impact of topology and geometry to future work. We fixed the lattice parameter a = 6.65 nm, the lipid length l = 2.365 nm, and the size of lipid heads h = 1 nm by following data reported in the literature [46,47] (see Methods for further details). It is expected that quantitative results depend on the choice of these parameters, but the qualitative impact of hydrophobicity and spatial dependence of diffusion on molecular transport does not change once the values of a , l and h have been fixed. We note that, for any system of interest, these parameters can be obtained directly from experiments such as small-angle X-ray or neutron scattering; particularly, a can be obtained by analyzing the peak positions of the scattering profiles [17], while l and l h correspond to half the thickness of the bilayer and of its hydrophobic part, respectively, which can be computed by suitably fitting the scattering profiles [67] or by employing geometric arguments based on the known composition of the system [1].
In order to establish a connection between the microscopic insights gathered by molecular dynamics simulations and the macroscopic transport properties relevant for experiments and practical applications, it is necessary to introduce a mesoscopic description of lipid mesophases, which takes as input the microscopic details of the system and provides predictions on the macroscopic diffusion. We can describe the heterogeneous environment offered by lipid mesophases at the nanoscopic scale by introducing a periodic potential of mean force U ( z ) and a space-dependent diffusion coefficient D ( z ) , both depending on the position z of the particle with respect to the the lipid–water interface (Figure 1). These features can be extracted from molecular dynamics simulations or, although only partially, from experimental data. To assess the impact of hydropathy, we first consider a minimalistic model of U ( z ) , for which we report a representative plot for a single periodicity in Figure 1b. The parameter Δ U gives the overall change in free energy upon inclusion of the guest molecule in the lipid membrane as compared to water, while Δ U b represents the free-energy barrier (or well, for negative values) which regulates the timescale of particle exchange between the lipid and water regions. Heterogeneous transport is captured by a position-dependent diffusion coefficient D ( z ) (Figure 1c), which continuously changes from D l i p in the lipid membrane to D w a t in water, with the change happening within a distance range of about w = 1 nm, as reported in the literature [8,40,44] (see Methods for further discussion on the matter). U ( z ) and D ( z ) were then employed to run Brownian Dynamics simulations aimed at assessing the diffusion kinetics of the system at large scales. Full details of the model and of the simulation setup are described in the Methods. In a later section, we also consider a practical case study focused on amino acids, for which more realistic potentials were extracted from molecular dynamics simulations reported in the literature.

3.1. Assessing the Importance of Each Physical Ingredient

3.1.1. Impact of Δ U b

The simulation setup enables devising systems which, albeit unrealistic, provide clear insights on the impact of each feature taken separately. In this section, we focus on how the barrier Δ U b affects large-scale diffusion. To this aim, we thus fix Δ U = 0 and D l i p = D w a t = 0.7 nm2/ns, so as to isolate the effect of Δ U b alone. Examples of potentials U ( z ) are reported in the insets of Figure 2. In the figure, the main plots show the mean-square displacement (MSD) for selected values of Δ U b , in order to highlight the role played by the magnitude and sign of this parameter. Particularly, in Figure 2a, we show the MSD for Δ U b = 4 k B T (green squares) and Δ U b = 4 k B T (orange circles). The empty symbols correspond to the MSD computed along planes parallel to the lipid/water interface, MSD. Since in the considered systems the diffusion coefficient is constant throughout space, lateral diffusion is unaffected by the value of U ( z ) . Therefore, MSD shows a standard diffusive behavior characterized by a diffusion coefficient D = D w a t , i.e., MSD =   4 D w a t t (grey dashed line in Figure 2a). In a log-log plot such as the ones considered in the figure, this corresponds to a linear function with slope one shifted according to the value of D .
In contrast, the kinetics along z (filled symbols in Figure 2a) are characterized by richer dynamics, where three regimes can be identified. At short times (e.g., t 0.1 ns for the orange circles in Figure 2a), the particles diffuse with a diffusion coefficient D = D w a t . In contrast, a subdiffusive behavior is detected at intermediate times (0.1 ns t 10 ns), as evidenced by a local slope lower than one. At longer times ( t 10 ns), standard diffusion is retrieved, but with an effective diffusion coefficient D < D w a t . Notably, the two sets of data collapse onto the same curve in this diffusive regime, suggesting that the long-time behavior of the system is independent of the sign of Δ U b , although the onset of this regime is shifted towards larger values of t for positive values of the barrier. In Figure 2b, we show that similar considerations apply for Δ U b = ± 8 k B T , but the larger magnitude of the barriers results in lower values of the effective diffusion coefficient D along the direction perpendicular to the lipid/water interface.
In Figure 3a, we report the long-time diffusion coefficients obtained by varying the magnitude of Δ U b in the range 1 10 k B T , for both positive (purple hexagons) and negative (golden stars) values of the barrier. Coherently with Figure 2, the lateral diffusion coefficients D (empty symbols) are independent of the height of the barrier, and correspond to the value of D w a t (grey dashed line). On the other hand, the perpendicular diffusion coefficients (filled symbols) are strongly affected by Δ U b , displaying an exponential decay starting from 5 k B T , for which a two-parameters best fit yields D ( n m 2 / n s ) = 8.0 e 0.93 Δ U b / k B T (red continuous line). This can be rationalized by observing that, for large values of Δ U b , jumping events from one side of the barrier to the other are rare. The rate k at which these events take place is approximately captured by the Arrhenius equation, k A e Δ U b / k B T . At time t (assumed to be large enough for the system to be found in the long-time diffusing regime), the cumulate number of expected events is k t . Each jump corresponds to a certain length λ within the same order of magnitude of the lattice parameter a. The total distance L travelled by the random set of jumps thus satisfies the relation L 2 = k t λ 2 = A e Δ U b / k B T λ 2 t . Since one has also L 2 = 2 D t , one thus obtains D = ( A / 2 ) e Δ U b / k B T λ 2 e Δ U b / k B T . Note that the Arrhenius formula predicts a prefactor equal to 1 in the exponent, which is in good agreement with the best-fitting value. Finally, the results obtained for Δ U b > 0 and Δ U b < 0 fall on the same master curve, confirming the trend observed in the long-time regime of M S D in Figure 2.

3.1.2. Impact of Δ U

We performed a similar study focused on the hydropathy parameter Δ U , for which a representative profile of U ( z ) is reported as an inset in Figure 3b. Qualitatively, the trend is very similar to the previous study. This is expected, since the system is being described by a set of periodically-placed barriers with equal heights, which is similar to the case reported in Figure 3a. The same reasoning as in the previous section applies, although the quantitative details of the exponential fit change. In this case, the fitting formula for the perpendicular diffusion coefficient is D ( n m 2 / n s ) = 2.6 e 0.97 Δ U / k B T (red continuous line in Figure 3b).

3.1.3. Impact of D l i p / D w a t

Finally, we also considered systems with a position dependent diffusion coefficient D ( z ) as in Figure 1c, but in which no potential was present, i.e., U ( z ) = 0 , which corresponds to Δ U b = Δ U = 0 . We report the long-time diffusion coefficients D (empty symbols) and D (filled symbols) in Figure 4. In contrast with the previous cases, here the z dependence of the local diffusion coefficient affects D , which is found to increase linearly with the ratio D l i p / D w a t . Considering that in the Brownian motion (Equation (3)) the drift term does not involve directly the parallel direction, it is expected that the long-term value of D is the average value of D ( z ) . From Equation (2), one thus finds
D = 1 a a 2 a 2 D ( z ) d z = D w a t ( D w a t D l i p ) 2 l + w a
which is reported in Figure 4 as a red continuous line, and quantitatively reproduces the simulation data.
As for the perpendicular direction, from Figure 4, one finds that D is systematically lower than D . The simulation results can be rationalized by considering the two extreme cases. When D l i p D w a t , one retrieves the diffusion coefficient obtained in the trivial case D ( z ) = D w a t , i.e., D = D w a t due to the absence of an external potential. In contrast, for D l i p D w a t , one expects that the particles are effectively confined along the perpendicular direction, i.e., that they do not diffuse transversally to the lipid/water interface. Indeed, in order to increase M S D , they need to traverse the full periodic repeat of the system; however, within the lipid region, there is practically no diffusion, so that the particle will need a large amount of time to cross it. In other terms, one expects in this limit that the transition to the long-time regime (see, e.g., Figure 2) moves towards larger and larger times, thus resulting in a vanishing value of D . Note that this reasoning does not apply to D , since there is no need for the system to cross the lipid phase in order to increase M S D . Quantitatively, we found that a good description of the data is obtained by means of a power law D = D w a t D l i p / D w a t 0.92 , which provides the correct limiting behavior and is reported as a blue dashed line in Figure 4.

3.2. Putting the Physical Ingredients Together

Having assessed the role played by each feature of the model, we now consider their interplay in more complex scenarios. In this regard, we performed simulations with the full potential U ( z ) (Figure 1a), by varying simultaneously the values of Δ U and Δ U b . We run two sets of simulations: in Set 1, the diffusion coefficient was considered to be constant ( D l i p = D w a t ), while, in Set 2, we implemented a more realistic profile for D ( z ) , with D l i p = 0.09 D w a t (see Methods for further details). The results of these simulations are reported in Figure 5 and Figure 6.
In Figure 5, we report the results obtained for D in the simulations of Set 1 (Figure 5a) and Set 2 (Figure 5b). In both cases, the largest values of D are obtained for lower values of Δ U b and Δ U , i.e., close to the center of each figure. This is intuitively understandable, as this region corresponds to lower barriers to be overcome. Similarly to what is observed in Figure 3, increasing the magnitude of Δ U b and Δ U has a dramatic effect on D , which rapidly decreases to low values (note that, in Figure 5, the scale in the bar is logarithmic). When comparing the results of the two sets of simulations for given values of Δ U b and Δ U , it is evident that, for Set 2, the value of D is systematically lower than for Set 1. This is also expected, since the average diffusion coefficient D for the variable case is lower. However, the overall change of D cannot trivially be ascribed to a normalization of the results by D , i.e., in general, D , S e t 2 D , S e t 1 · D .
As for the lateral diffusion coefficient D , in the simulations of Set 1, one finds trivially that D = D w a t , analogous to what was reported in Figure 3. The results obtained for Set 2 are instead reported in Figure 6a. In this case, the largest values are obtained for Δ U b > 0 , Δ U > 0 with a large magnitude, i.e., by maximizing the depletion from the lipid phase. This can be understood by observing that, to achieve an efficient lateral diffusion, one does not need the particles to cross a full periodicity, but rather to maximize the time spent in the water phase, characterized by a larger mobility (Figure 1c). This is best achieved by increasing the energy penalty for particle localization in the lipid phase, i.e., by large, positive values of Δ U b and Δ U (Figure 1b). Quantitatively, as discussed above, the absence of a direct drift term in the parallel direction suggests that D can be identified with the average D . Hence, in the presence of a potential U ( z ) , one can generalize Equation (12) as
D = 1 a a 2 a 2 e U ( z ) k B T D ( z ) d z .
To enable a direct link with experimentally-measurable quantities, we quantify the relative amount of time spent in the lipid phase by means of the logarithm of the partition coefficient, log P = log 10 ( c l i p / c w a t ) [62], where c l i p and c w a t are the concentrations in the lipid and water phase. Rewriting Equation (13) by means of the partition coefficient yields (see Methods),
D = D l i p l 10 log P + w 2 + D w a t a 2 l w 2 l 10 log P + a 2 l .
Note that, for strongly hydrophobic molecules, log P is positive and has a large magnitude, so that, in the previous formula, the terms containing the factor 10 log P are much larger than the rest. In this case, one thus obtains D D l i p , which is expected since the molecule spends virtually all the time within the lipid bilayer. In contrast, for strongly hydrophilic molecules log P < 0 and large in magnitude, so that 10 log P 0 . The value obtained for D in this case is the average diffusion coefficient within the water phase, which does not trivially correspond to D w a t due to the position dependence of D ( z ) (Figure 1c). In Figure 6b, we compare the numerical results from the simulations (golden stars) with Equation (14) (red continuous line). The excellent agreement confirms the quantitative correspondence between D and D . To fully appreciate the dependence of D on log P , in Figure 6b, we also considered an extended set of systems (purple triangles), with Δ U and Δ U b going beyond the maximum magnitude 5 k B T considered in the simulations. For these systems, D was computed according to Equation (13).

3.3. Application: Large-Scale Transport of Amino Acids

As a practical example of the usage of the present approach, we dedicate this section to the study of diffusion of amino acids through lamellar phases. The potential of mean force U ( z ) for amino acids within phospholipidic bilayers has been the focus of previous investigation work [63,68]. Rather than using the toy model proposed in Figure 1b, we thus consider here U ( z ) for 16 amino acids as computed in Ref. [63], where the authors performed enhanced-sampling molecular dynamics of various residues embedded in DOPC bilayers. The resulting potentials are reported in Figure 7a. Although being more complex, they resemble the toy model introduced in Figure 1b, showing a plateau in correspondence of the lipid tails (z close to zero) and barriers or wells in proximity of the lipid heads ( | z | around 2–2.5 nm). Albeit being accessible to atomistic simulations [69], the position dependence of the diffusion coefficient was unfortunately not addressed in Ref. [63], so that we consider the toy model for D ( z ) presented in Figure 1c. Further details of the simulations are provided in the Methods.
Based on our analysis in the previous sections, we expect that the presence of barriers in U ( z ) has a dramatic effect on the perpendicular diffusion coefficient D . To quantify the extent of such barriers, for each amino acid, we computed the minimum and maximum values of the potential of mean force, which we denote as U m i n and U m a x , respectively. For instance, for arginine, we obtain U m i n 8.5 k B T and U m a x 23.5 k B T (second panel from the top-left corner in Figure 7a). In Figure 7b, we report D as a function of the difference U m a x U m i n , which we took as an indicator of the strength of the barriers. The heterogeneity of profiles for U ( z ) results in extremely different values for D , which span as many as ten orders of magnitude (compare main plot and inset in Figure 7b). Moreover, an Arrhenius-like exponential dependence captures the data over all the different scales, with the best-fitting formula being D 0.5 · e 0.77 ( U m a x U m i n ) / k B T (grey dashed line). The data cluster according to the physico-chemical properties of the amino acids. Particularly, the ones with charged or polar side chains are characterized by the lowest values of D (golden stars and red circles in Figure 7b, respectively). This is not expected a priori, since free-energy barriers are expected also for apolar residues. However, as stressed by the different values of U m a x U m i n obtained for the various sets of amino acids (see also full profiles of U ( z ) in Figure 7a), the depletion of charged and polar residues from the lipid bilayer is significantly stronger than the free energy gained by embedding apolar amino acids. A closer look to the free-energy profiles indicates that, for apolar residues, localization within the whole lipid region is either energetically favourable or it comes at a negligible cost. The barrier for perpendicular diffusion is thus provided by the depth of the free-energy wells. In contrast, for charged and polar residues, localization in the region corresponding to the lipid tails is highly costly, while the lipid heads are favoured or have low associated cost. Overall, this provides a significantly steeper barrier to perpendicular diffusion.
As for D , analogous to Figure 6b, in Figure 7c, we plot the results obtained from the simulations as a function of log P . The grey dashed line is the prediction according to Equation (14), which is again found to quantitatively capture the simulation data. Importantly, D is always found within the same order of magnitude ( 10 1 nm2/ns), which indicates that D D for large enough barriers, hence suggesting parallel diffusion to be dominant in such scenario. Similar to D , also in this case, we find a clustering of the points according to the physico-chemical properties of the residues. Particularly, apolar residues (green empty squares) have a stronger affinity for the lipid phase, as denoted by the larger values of log P (taken with sign). This implies that a larger fraction of time is spent in the slowly-diffusing region corresponding to the lipids, thus yielding low values of D . Similarly, polar residues (red empty circles) have lower affinity for lipids, thus resulting in faster parallel diffusion. Unintuitive results are obtained instead for charged residues (golden empty stars), for which one would imagine a strong depletion from the lipid phase. While Glu (E) and Asp (D) abide by the expected behavior (low affinity for lipids, a large value of D ), quite surprisingly, Lys (K) and Arg (R) show instead the opposite behavior. An inspection of the corresponding U ( z ) profiles for these residues (Figure 7a) reveals the presence of a deep well ( 10 k B T ) in correspondence with the lipid heads. Thus, despite being strongly depleted from the lipid tails, Lys and Arg spend most of the time sitting at the lipid/water interface, which is characterized by slow diffusion. In terms of the toy model, one can locate these residues in Figure 6a in correspondence with Δ U > 0 and Δ U b < 0 , both with large magnitudes. However, it is important to stress that the quantitative value of D for Lys and Arg (and more in general for all the amino acids with positive values of log P ) strongly depends on the profile chosen for D ( z ) (Figure 1c) and on the value of D l i p , so that a more precise quantitative estimation requires a proper determination of D ( z ) from ad hoc atomistic simulations. The cases of Lys and Arg are also an instructive example of how the size of barriers ( U m a x U m i n ) and the hydropathy of the molecule ( log P ) are not necessarily correlated with each other.

4. Discussion

Wrapping up, the systematic study of the toy model and the simulation of the diffusion of amino acids lead us to the following key points:
  • The key determinant for perpendicular diffusion is the overall height of the free-energy barriers. For barrier heights larger than roughly 5 k B T , one finds that D decreases exponentially with the size of the barrier (Figure 3, Figure 5 and Figure 7a).
  • Parallel diffusion is determined by the relative time spent in the lipid phase as compared to water, which provides a direct relation between D and the partition coefficient log P (Equation (14), Figure 6 and Figure 7b);
  • The lower boundary of D is equal to D l i p , obtained for highly-hydrophobic guest molecules. Together with point 1, this indicates that, for large enough barriers, parallel diffusion is dominant.
Experimental assays, such as pulsed-field gradient NMR [26] or macroscopic release setups [17], can access a macroscopic, three-dimensional diffusion coefficient D eff , probing scales equal to or larger than microns. How is D eff related to the diffusion coefficients considered in the simulations? It is key to observe that, although being ordered at the nanoscale, lipid mesophases are formed by micrometer-sized domains separated by grain boundaries, which usually lack orientational order at larger scales [70]. Hence, assuming neighboring domains to be randomly oriented with respect to each other, there is no net distinction between parallel and perpendicular diffusion at experimentally-relevant scales. The assumed random orientation enables accounting for large-scale diffusion by averaging over the various domains, hence the particles are expected to experience three-dimensional diffusion with an effective diffusion coefficient D eff = 2 D / 3 + D / 3 . The weights 2 / 3 and 1 / 3 associated with parallel and perpendicular diffusion are chosen so as to account properly for the dimensionality of the corresponding process. For large barriers, one can thus approximate D eff = 2 D / 3 ; by means of Equation (14), one obtains
D eff = 2 3 D l i p l 10 log P + w 2 + D w a t a 2 l w 2 l 10 log P + a 2 l .
Remarkably, Equation (15) enables estimating the macroscopic diffusion coefficient D eff from knowledge of local geometrical ( a , w , l ) and transport ( D l i p , D w a t ) properties, as well as from the overall thermodynamic equilibrium distribution ( log P ). In Figure 8, we test the accuracy of Equation (15) by comparing its prediction with the value obtained directly from the simulations, both for the toy model (panel a) and for the amino-acids’ simulations (panel b). For barriers larger than 5 k B T , the prediction gives values within 10% of the numerical ones. Nevertheless, it should be kept in mind that the threshold value of the barrier for which the approximation works reasonably is expected to depend on the numerical values chosen for the parameters, particularly for the ratio D l i p / D w a t .
As seen from Figure 8b, in the case of amino acids, the approximation works quite well in virtually all cases. From the literature, typical values of the barriers for various drugs are typically well beyond the threshold value of 5 k B T . For instance, molecular dynamics simulations have accessed the potential of mean force for embedding paracetamol in DPPC, obtaining a barrier 50 k B T [47]. Similarly, for aspirin and ibuprofen in DPPC, the free-energy barriers have been estimated to be 20 30 k B T [71]. These values suggest that, in practical applications, parallel diffusion is often dominant, thus enabling the employment of the approximation given by Equation (15). Based on this, as a further example, we collected experimental values of log P and D w a t for various drugs from the literature, and plugged them into Equation (15) to estimate the effective diffusion coefficient D eff characterizing the release from a L α lamellar mesophase with a 6.7 nm, l 2.4 m, and w 1 nm. These values are the geometrical parameters fixed in the toy model, and correspond to experimentally-observed structures for water/DPPC mixtures at 43 °C [46]. The estimations of D eff are listed in Table 1, and can provide a reference for researchers interested in studying the release of these drugs from lipidic mesophases. Following our treatment, we assume that, for each case, D l i p = 0.09 D w a t .
Although providing a useful starting point, the results reported in Table 1 are, however, only estimations based on the limited knowledge provided by the experimental partition coefficient. As mentioned above, the ideal strategy is to combine the present approach with atomistic simulations of the system of interest, so as to compute the free-energy profile U ( z ) and the diffusion profile D ( z ) to be used as input [69]. As a word of caution, some care has to be taken when dealing with molecules with pKa values close to the pH of the solution. Indeed, in this case, the possibility of dynamic protonation has to be taken into account. For instance, for the charged amino acids in Figure 7a, we considered for U ( z ) the value corresponding to the charged or neutral state according to their relative thermodynamic stability, which changes with z [63]. The best practice is to run constant-pH atomistic simulations [81], which yield directly the correct free-energy profile [64]. However, we believe that, in practice, this will result in small changes in the predicted transport values, as a molecule with preference for charged states (in pure water) tends to be depleted by the hydropobic center of the lipid bilayer, thus resulting in a high free-energy barrier for bilayer permeation [82]. In turn, this yields negligible values for D (Figure 7b), making parallel diffusion the dominant transport mechanism. While charge neutralization lowers the barrier at the center of the bilayer, it is expected that such barrier is still present and large in magnitude [82]; hence, while promoting perpendicular diffusion, the value of D is still expected to be too small to significantly affect overall diffusion.
To summarize, we have presented a multiscale approach to predict the macroscopic diffusion coefficient by a combination of atomistic simulations, from which the profiles for U ( z ) and D ( z ) (Figure 1b,c) can be extracted, and Brownian-dynamics simulations, which enable access to the effective large-scale diffusion emerging from the interplay of the nanoscopic features. Based on a minimalistic toy model and a case study focused on amino acids, we have discerned the impact of the main dynamic and thermodynamic features of the system on molecular transport at macroscopic scales. A further possibility is to use our results the other way round: from experimental knowledge of D eff , D w a t and log P , and under the assumption of large barriers and a diffusion profile with the shape considered here (Figure 1c), one can access the value of D l i p , characterizing diffusion of the inspected molecule within the lipid bilayer. Future work will consider ad hoc studies to obtain the detailed shape of D ( z ) for selected systems based on established procedures from the literature [69]. Moreover, we will also adapt the present framework to more complex topologies of lipid mesophases of direct relevance for release studies, including, for instance, hexagonal and cubic phases. From a wider perspective, the great potential of multiscale simulations is already being exploited in affine fields, such as the study of biological membranes [83] or the pathways of antibiotic intake from bacteria [84]; it is our hope that the results presented here will spark a similar interest for multiscale simulations in the field of controlled release from lipid mesophases, thus paving the way for the development of an invaluable complement to experimental assays.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pharmaceutics15020573/s1. Raw data used to generate the plots.

Author Contributions

Conceptualization, writing—original draft preparation, supervision, S.A.; investigation, analysis, writing—review and editing, A.M.B. and S.A. All authors have read and agreed to the published version of the manuscript.

Funding

The project that gave rise to these results received the support of a fellowship from “la Caixa” Foundation (ID 100010434) and from the European Union’s Horizon research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 847648. The fellowship code is LCF/BQ/PI20/11760019.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The files containing the data used to generate the plots are available as Supplementary Materials.

Acknowledgments

We thank Thomas Bourke for insightful discussions inspiring this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Briggs, J.; Chung, H.; Caffrey, M. The temperature-composition phase diagram and mesophase structure characterization of the monoolein/water system. J. Phys. II 1996, 6, 723–751. [Google Scholar]
  2. Templer, R.; Seddon, J.; Warrender, N.; Syrykh, A.; Huang, Z.; Winter, R.; Erbes, J. Inverse bicontinuous cubic phases in 2: 1 fatty acid/phosphatidylcholine mixtures. The effects of chain length, hydration, and temperature. J. Phys. Chem. B 1998, 102, 7251–7261. [Google Scholar]
  3. Tyler, A.I.; Barriga, H.M.; Parsons, E.S.; McCarthy, N.L.; Ces, O.; Law, R.V.; Seddon, J.M.; Brooks, N.J. Electrostatic swelling of bicontinuous cubic lipid phases. Soft Matter 2015, 11, 3279–3286. [Google Scholar] [PubMed]
  4. Mezzenga, R.; Meyer, C.; Servais, C.; Romoscanu, A.I.; Sagalowicz, L.; Hayward, R.C. Shear rheology of lyotropic liquid crystals: A case study. Langmuir 2005, 21, 3322–3333. [Google Scholar]
  5. Negrini, R.; Mezzenga, R. Diffusion, molecular separation, and drug delivery from lipid mesophases with tunable water channels. Langmuir 2012, 28, 16455–16462. [Google Scholar] [PubMed]
  6. Barauskas, J.; Johnsson, M.; Tiberg, F. Self-assembled lipid superstructures: Beyond vesicles and liposomes. Nano Lett. 2005, 5, 1615–1619. [Google Scholar]
  7. Turner, D.R.; Wang, Z.G.; Gruner, S.; Mannock, D.; Mcelhaney, R. Structural study of the inverted cubic phases of di-dodecyl alkyl-β-D-glucopyranosyl-rac-glycerol. J. Phys. II 1992, 2, 2039–2063. [Google Scholar] [CrossRef]
  8. Salvati Manni, L.; Assenza, S.; Duss, M.; Vallooran, J.J.; Juranyi, F.; Jurt, S.; Zerbe, O.; Landau, E.M.; Mezzenga, R. Soft biomimetic nanoconfinement promotes amorphous water over ice. Nat. Nanotechnol. 2019, 14, 609–615. [Google Scholar]
  9. Landau, E.M.; Rosenbusch, J.P. Lipidic cubic phases: A novel concept for the crystallization of membrane proteins. Proc. Natl. Acad. Sci. USA 1996, 93, 14532–14535. [Google Scholar]
  10. Mulet, X.; Boyd, B.J.; Drummond, C.J. Advances in drug delivery and medical imaging using colloidal lyotropic liquid crystalline dispersions. J. Colloid Interface Sci. 2013, 393, 1–20. [Google Scholar]
  11. Mezzenga, R.; Schurtenberger, P.; Burbidge, A.; Michel, M. Understanding foods as soft materials. Nat. Mater. 2005, 4, 729–740. [Google Scholar] [CrossRef] [PubMed]
  12. Fong, C.; Le, T.; Drummond, C.J. Lyotropic liquid crystal engineering–ordered nanostructured small molecule amphiphile self-assembly materials by design. Chem. Soc. Rev. 2012, 41, 1297–1322. [Google Scholar] [PubMed]
  13. Speziale, C.; Salvati Manni, L.; Manatschal, C.; Landau, E.M.; Mezzenga, R. A macroscopic H+ and Cl− ions pump via reconstitution of EcClC membrane proteins in lipidic cubic mesophases. Proc. Natl. Acad. Sci. USA 2016, 113, 7491–7496. [Google Scholar] [PubMed]
  14. Assenza, S.; Mezzenga, R. Soft condensed matter physics of foods and macronutrients. Nat. Rev. Phys. 2019, 1, 551–566. [Google Scholar] [CrossRef]
  15. Yao, Y.; Zhou, T.; Färber, R.; Grossner, U.; Floudas, G.; Mezzenga, R. Designing cryo-enzymatic reactions in subzero liquid water by lipidic mesophase nanoconfinement. Nat. Nanotechnol. 2021, 16, 802–810. [Google Scholar]
  16. Boyd, B.J.; Clulow, A.J. The influence of lipid digestion on the fate of orally administered drug delivery vehicles. Biochem. Soc. Trans. 2021, 49, 1749–1761. [Google Scholar]
  17. Negrini, R.; Mezzenga, R. pH-responsive lyotropic liquid crystals for controlled drug delivery. Langmuir 2011, 27, 5296–5303. [Google Scholar] [CrossRef]
  18. Martiel, I.; Baumann, N.; Vallooran, J.J.; Bergfreund, J.; Sagalowicz, L.; Mezzenga, R. Oil and drug control the release rate from lyotropic liquid crystals. J. Control. Release 2015, 204, 78–84. [Google Scholar]
  19. Zabara, A.; Mezzenga, R. Controlling molecular transport and sustained drug release in lipid-based liquid crystalline mesophases. J. Control. Release 2014, 188, 31–43. [Google Scholar]
  20. Negrini, R.; Fong, W.K.; Boyd, B.J.; Mezzenga, R. pH-responsive lyotropic liquid crystals and their potential therapeutic role in cancer treatment. Chem. Commun. 2015, 51, 6671–6674. [Google Scholar] [CrossRef]
  21. Nazaruk, E.; Miszta, P.; Filipek, S.; Gorecka, E.; Landau, E.M.; Bilewicz, R. Lyotropic cubic phases for drug delivery: Diffusion and sustained release from the mesophase evaluated by electrochemical methods. Langmuir 2015, 31, 12753–12761. [Google Scholar] [CrossRef] [PubMed]
  22. Clogston, J.; Craciun, G.; Hart, D.; Caffrey, M. Controlling release from the lipidic cubic phase by selective alkylation. J. Control. Release 2005, 102, 441–461. [Google Scholar] [CrossRef] [PubMed]
  23. Lee, K.W.; Nguyen, T.H.; Hanley, T.; Boyd, B.J. Nanostructure of liquid crystalline matrix determines in vitro sustained release and in vivo oral absorption kinetics for hydrophilic model drugs. Int. J. Pharm. 2009, 365, 190–199. [Google Scholar] [PubMed]
  24. Phan, S.; Fong, W.K.; Kirby, N.; Hanley, T.; Boyd, B.J. Evaluating the link between self-assembled mesophase structure and drug release. Int. J. Pharm. 2011, 421, 176–182. [Google Scholar] [CrossRef] [PubMed]
  25. Zabara, A.; Negrini, R.; Onaca-Fischer, O.; Mezzenga, R. Perforated Bicontinuous Cubic Phases with pH-Responsive Topological Channel Interconnectivity. Small 2013, 9, 3602–3609. [Google Scholar] [CrossRef]
  26. Meikle, T.G.; Yao, S.; Zabara, A.; Conn, C.E.; Drummond, C.J.; Separovic, F. Predicting the release profile of small molecules from within the ordered nanostructured lipidic bicontinuous cubic phase using translational diffusion coefficients determined by PFG-NMR. Nanoscale 2017, 9, 2471–2478. [Google Scholar]
  27. Ghanbari, R.; Assenza, S.; Saha, A.; Mezzenga, R. Diffusion of polymers through periodic networks of lipid-based nanochannels. Langmuir 2017, 33, 3491–3498. [Google Scholar]
  28. Assenza, S.; Mezzenga, R. Curvature and bottlenecks control molecular transport in inverse bicontinuous cubic phases. J. Chem. Phys. 2018, 148, 054902. [Google Scholar] [CrossRef]
  29. Ghanbari, R.; Assenza, S.; Mezzenga, R. The interplay of channel geometry and molecular features determines diffusion in lipidic cubic phases. J. Chem. Phys. 2019, 150, 094901. [Google Scholar]
  30. Ghanbari, R.; Assenza, S.; Zueblin, P.; Mezzenga, R. Impact of Molecular Partitioning and Partial Equilibration on the Estimation of Diffusion Coefficients from Release Experiments. Langmuir 2019, 35, 5663–5671. [Google Scholar] [CrossRef]
  31. Salvati Manni, L.; Duss, M.; Assenza, S.; Boyd, B.J.; Landau, E.M.; Fong, W.K. Enzymatic hydrolysis of monoacylglycerols and their cyclopropanated derivatives: Molecular structure and nanostructure determine the rate of digestion. J. Colloid Interface Sci. 2021, 588, 767–775. [Google Scholar] [PubMed]
  32. Chang, C.; Meikle, T.G.; Drummond, C.J.; Yang, Y.; Conn, C.E. Comparison of cubosomes and liposomes for the encapsulation and delivery of curcumin. Soft Matter 2021, 17, 3306–3313. [Google Scholar] [CrossRef] [PubMed]
  33. Rajesh, S.; Zhai, J.; Drummond, C.J.; Tran, N. Novel pH-Responsive Cubosome and Hexosome Lipid Nanocarriers of SN-38 Are Prospective for Cancer Therapy. Pharmaceutics 2022, 14, 2175. [Google Scholar] [CrossRef] [PubMed]
  34. Dubbeldam, D.; Snurr, R. Recent developments in the molecular modeling of diffusion in nanoporous materials. Mol. Simul. 2007, 33, 305–325. [Google Scholar] [CrossRef]
  35. Krishna, R. Diffusion in porous crystalline materials. Chem. Soc. Rev. 2012, 41, 3099–3118. [Google Scholar] [PubMed]
  36. Bagchi, B. Water dynamics in the hydration layer around proteins and micelles. Chem. Rev. 2005, 105, 3197–3219. [Google Scholar] [PubMed]
  37. Laage, D.; Elsaesser, T.; Hynes, J.T. Water dynamics in the hydration shells of biomolecules. Chem. Rev. 2017, 117, 10694–10725. [Google Scholar]
  38. Bourg, I.C.; Steefel, C.I. Molecular dynamics simulations of water structure and diffusion in silica nanopores. J. Phys. Chem. C 2012, 116, 11556–11564. [Google Scholar]
  39. Hande, V.R.; Chakrabarty, S. How Far Is “Bulk Water” from Interfaces? Depends on the Nature of the Surface and What We Measure. J. Phys. Chem. B 2022, 126, 1125–1135. [Google Scholar]
  40. Vallooran, J.J.; Assenza, S.; Mezzenga, R. Spatiotemporal control of enzyme-induced crystallization under lyotropic liquid crystal nanoconfinement. Angew. Chem. 2019, 131, 7367–7371. [Google Scholar]
  41. Kasim, N.A.; Whitehouse, M.; Ramachandran, C.; Bermejo, M.; Lennernäs, H.; Hussain, A.S.; Junginger, H.E.; Stavchansky, S.A.; Midha, K.K.; Shah, V.P.; et al. Molecular properties of WHO essential drugs and provisional biopharmaceutical classification. Mol. Pharm. 2004, 1, 85–96. [Google Scholar] [PubMed]
  42. Takagi, T.; Ramachandran, C.; Bermejo, M.; Yamashita, S.; Yu, L.X.; Amidon, G.L. A provisional biopharmaceutical classification of the top 200 oral drug products in the United States, Great Britain, Spain, and Japan. Mol. Pharm. 2006, 3, 631–643. [Google Scholar] [CrossRef] [PubMed]
  43. Eason, T.; Ramirez, G.; Clulow, A.J.; Salim, M.; Boyd, B.J. Revisiting the Dissolution of Praziquantel in Biorelevant Media and the Impact of Digestion of Milk on Drug Dissolution. Pharmaceutics 2022, 14, 2228. [Google Scholar] [PubMed]
  44. Kim, J.; Lu, W.; Qiu, W.; Wang, L.; Caffrey, M.; Zhong, D. Ultrafast hydration dynamics in the lipidic cubic phase: Discrete water structures in nanochannels. J. Phys. Chem. B 2006, 110, 21994–22000. [Google Scholar] [CrossRef]
  45. Helfrich, W. Effect of thermal undulations on the rigidity of fluid membranes and interfaces. J. Phys. 1985, 46, 1263–1268. [Google Scholar] [CrossRef]
  46. Angelov, B.; Ollivon, M.; Angelova, A. X-ray diffraction study of the effect of the detergent octyl glucoside on the structure of lamellar and nonlamellar lipid/water phases of use for membrane protein reconstitution. Langmuir 1999, 15, 8225–8234. [Google Scholar]
  47. Nademi, Y.; Amjad Iranagh, S.; Yousefpour, A.; Mousavi, S.Z.; Modarress, H. Molecular dynamics simulations and free energy profile of Paracetamol in DPPC and DMPC lipid bilayers. J. Chem. Sci. 2014, 126, 637–647. [Google Scholar]
  48. Zhou, T.; Vallooran, J.J.; Assenza, S.; Szekrenyi, A.; Clapés, P.; Mezzenga, R. Efficient asymmetric synthesis of carbohydrates by aldolase nano-confined in lipidic cubic mesophases. ACS Catal. 2018, 8, 5810–5815. [Google Scholar]
  49. Ma, Y.; Zhu, C.; Ma, P.; Yu, K. Studies on the diffusion coefficients of amino acids in aqueous solutions. J. Chem. Eng. Data 2005, 50, 1192–1196. [Google Scholar]
  50. Longsworth, L. Diffusion measurements, at 25, of aqueous solutions of amino acids, peptides and sugars. J. Am. Chem. Soc. 1953, 75, 5705–5709. [Google Scholar]
  51. Ribeiro, A.C.; Barros, M.C.; Verissimo, L.M.; Lobo, V.M.; Valente, A.J. Binary diffusion coefficients for aqueous solutions of l-aspartic acid and its respective monosodium salt. J. Solut. Chem. 2014, 43, 83–92. [Google Scholar] [CrossRef]
  52. Ribeiro, A.C.; Rodrigo, M.; Barros, M.C.; Verissimo, L.M.; Romero, C.; Valente, A.J.; Esteso, M.A. Mutual diffusion coefficients of L-glutamic acid and monosodium L-glutamate in aqueous solutions at T = 298.15 K. J. Chem. Thermodyn. 2014, 74, 133–137. [Google Scholar]
  53. Umecky, T.; Ehara, K.; Omori, S.; Kuga, T.; Yui, K.; Funazukuri, T. Binary diffusion coefficients of aqueous phenylalanine, tyrosine isomers, and aminobutyric acids at infinitesimal concentration and temperatures from (293.2 to 333.2) K. J. Chem. Eng. Data 2013, 58, 1909–1917. [Google Scholar]
  54. Mendes, F.S.; Cruz, C.E.; Martins, R.N.; Ramalho, J.P.P.; Martins, L.F. On the diffusion of ketoprofen and ibuprofen in water: An experimental and theoretical approach. J. Chem. Thermodyn. 2022, 2022, 106955. [Google Scholar] [CrossRef]
  55. Edwards, L. The dissolution and diffusion of aspirin in aqueous media. Trans. Faraday Soc. 1951, 47, 1191–1210. [Google Scholar] [CrossRef]
  56. Ribeiro, A.C.; Barros, M.C.; Veríssimo, L.M.; Santos, C.I.; Cabral, A.M.; Gaspar, G.D.; Esteso, M.A. Diffusion coefficients of paracetamol in aqueous solutions. J. Chem. Thermodyn. 2012, 54, 97–99. [Google Scholar] [CrossRef]
  57. Antognini, L.M.; Assenza, S.; Speziale, C.; Mezzenga, R. Quantifying the transport properties of lipid mesophases by theoretical modelling of diffusion experiments. J. Chem. Phys. 2016, 145, 084903. [Google Scholar]
  58. Anderson, D.M.; Wennerstroem, H. Self-diffusion in bicontinuous cubic phases, L3 phases, and microemulsions. J. Phys. Chem. 1990, 94, 8683–8694. [Google Scholar] [CrossRef]
  59. Doi, M.; Edwards, S.F.; Edwards, S.F. The Theory of Polymer Dynamics; Oxford University Press: Oxford, UK, 1988; Volume 73. [Google Scholar]
  60. Van Kampen, N.G. Stochastic Processes in Physics and Chemistry; Elsevier: Amsterdam, The Netherlands, 1992; Volume 1. [Google Scholar]
  61. Marsaglia, G.; Tsang, W.W. The ziggurat method for generating random variables. J. Stat. Softw. 2000, 5, 1–7. [Google Scholar] [CrossRef]
  62. Martin, Y.C. Quantitative Drug Design: A Critical Introduction; CRC Press: Boca Raton, FL, USA, 2010. [Google Scholar]
  63. MacCallum, J.L.; Bennett, W.D.; Tieleman, D.P. Distribution of amino acids in a lipid bilayer from computer simulations. Biophys. J. 2008, 94, 3393–3404. [Google Scholar] [CrossRef]
  64. Yue, Z.; Li, C.; Voth, G.A.; Swanson, J.M. Dynamic protonation dramatically affects the membrane permeability of drug-like molecules. J. Am. Chem. Soc. 2019, 141, 13421–13433. [Google Scholar] [PubMed]
  65. Koirala, R.P.; Bhusal, H.P.; Khanal, S.P.; Adhikari, N.P. Effect of temperature on transport properties of cysteine in water. AIP Adv. 2020, 10, 025122. [Google Scholar]
  66. Rodrigo, M.M.; Valente, A.J.; Barros, M.C.; Verissimo, L.M.; Romero, C.; Esteso, M.A.; Ribeiro, A.C. Mutual diffusion coefficients of L-lysine in aqueous solutions. J. Chem. Thermodyn. 2014, 74, 227–230. [Google Scholar] [CrossRef]
  67. Heller, W.T. Small-Angle Neutron Scattering for Studying Lipid Bilayer Membranes. Biomolecules 2022, 12, 1591. [Google Scholar] [PubMed]
  68. Johansson, A.C.; Lindahl, E. Position-resolved free energy of solvation for amino acids in lipid membranes from molecular dynamics simulations. Proteins Struct. Funct. Bioinform. 2008, 70, 1332–1344. [Google Scholar] [CrossRef] [PubMed]
  69. Zhu, F.; Hummer, G. Theory and simulation of ion conduction in the pentameric GLIC channel. J. Chem. Theory Comput. 2012, 8, 3759–3768. [Google Scholar] [CrossRef] [PubMed]
  70. Vallooran, J.J.; Bolisetty, S.; Mezzenga, R. Macroscopic alignment of lyotropic liquid crystals using magnetic nanoparticles. Adv. Mater. 2011, 23, 3932–3937. [Google Scholar] [PubMed]
  71. Boggara, M.B.; Krishnamoorti, R. Partitioning of nonsteroidal antiinflammatory drugs in lipid membranes: A molecular dynamics simulation study. Biophys. J. 2010, 98, 586–595. [Google Scholar] [PubMed]
  72. Keyvanfard, M.; Hatami, M.; Gupta, V.K.; Agarwal, S.; Sadeghifar, H.; Khalilzadeh, M.A. Liquid phase analysis of methyldopa in the presence of tyrosine using electrocatalytic effect of a catechol derivative at a surface of NiO nanoparticle modified carbon paste electrode. J. Mol. Liq. 2017, 230, 290–294. [Google Scholar] [CrossRef]
  73. Fagerholm, U.; Lennernäs, H. Experimental estimation of the effective unstirred water layer thickness in the human jejunum, and its importance in oral drug absorption. Eur. J. Pharm. Sci. 1995, 3, 247–253. [Google Scholar] [CrossRef]
  74. Lu, Y.; Li, M. Simultaneous rapid determination of the solubility and diffusion coefficients of a poorly water-soluble drug based on a novel UV imaging system. J. Pharm. Sci. 2016, 105, 131–138. [Google Scholar] [PubMed]
  75. Vilt, M.E.; Ho, W.W. Supported liquid membranes with strip dispersion for the recovery of Cephalexin. J. Membr. Sci. 2009, 342, 80–87. [Google Scholar] [CrossRef]
  76. Tajik, S.; Beitollahi, H.; Aflatoonian, M.R.; Mohtat, B.; Aflatoonian, B.; Shoaie, I.S.; Khalilzadeh, M.A.; Ziasistani, M.; Zhang, K.; Jang, H.W.; et al. Fabrication of magnetic iron oxide-supported copper oxide nanoparticles (Fe3O4/CuO): Modified screen-printed electrode for electrochemical studies and detection of desipramine. RSC Adv. 2020, 10, 15171–15178. [Google Scholar] [CrossRef]
  77. Beitollahi, H.; Hamzavi, M.; Torkzadeh-Mahani, M. Electrochemical determination of hydrochlorothiazide and folic acid in real samples using a modified graphene oxide sheet paste electrode. Mater. Sci. Eng. C 2015, 52, 297–305. [Google Scholar] [CrossRef]
  78. Pyka, A. Lipophilicity investigations of ibuprofen. J. Liquid Chromatogr. Related Technol.® 2009, 32, 723–731. [Google Scholar] [CrossRef]
  79. Barros, M.C.; Ribeiro, A.C.; Esteso, M.A.; Lobo, V.M.; Leaist, D.G. Diffusion of levodopa in aqueous solutions of hydrochloric acid at 25 °C. J. Chem. Thermodyn. 2014, 72, 44–47. [Google Scholar] [CrossRef]
  80. Ye, F.; Yaghmur, A.; Jensen, H.; Larsen, S.W.; Larsen, C.; Østergaard, J. Real-time UV imaging of drug diffusion and release from Pluronic F127 hydrogels. Eur. J. Pharm. Sci. 2011, 43, 236–243. [Google Scholar]
  81. Radak, B.K.; Chipot, C.; Suh, D.; Jo, S.; Jiang, W.; Phillips, J.C.; Schulten, K.; Roux, B. Constant-pH molecular dynamics simulations for large biomolecular systems. J. Chem. Theory Comput. 2017, 13, 5933–5944. [Google Scholar] [CrossRef]
  82. Johansson, A.C.; Lindahl, E. Titratable amino acid solvation in lipid membranes as a function of protonation state. J. Phys. Chem. B 2009, 113, 245–253. [Google Scholar]
  83. Enkavi, G.; Javanainen, M.; Kulig, W.; Róg, T.; Vattulainen, I. Multiscale simulations of biological membranes: The challenge to understand biological phenomena in a living substance. Chem. Rev. 2019, 119, 5607–5774. [Google Scholar]
  84. Prajapati, J.D.; Kleinekathöfer, U.; Winterhalter, M. How to enter a bacterium: Bacterial porins and the permeation of antibiotics. Chem. Rev. 2021, 121, 5158–5192. [Google Scholar] [PubMed]
Figure 1. General features of the theoretical model. (a) sketch of a repeating unit in a lamellar mesophase and definition of the main parameters; molecules are not at scale; (b) representative periodic potential of mean force U ( z ) corresponding to the system in (a). In this case, U ( z ) corresponds to a hydrophilic molecule ( Δ U > 0 ) with low affinity for the lipid heads ( Δ U b > 0 ). The orange circles correspond to the free energy extracted from a control simulation to determine the optimal timestep, as described in the Methods; (c) representative periodic, position-dependent diffusion coefficient D ( z ) corresponding to the system in (a), smoothly changing from D l i p within the lipid phase to D w a t in the water phase far from the lipid/water interface.
Figure 1. General features of the theoretical model. (a) sketch of a repeating unit in a lamellar mesophase and definition of the main parameters; molecules are not at scale; (b) representative periodic potential of mean force U ( z ) corresponding to the system in (a). In this case, U ( z ) corresponds to a hydrophilic molecule ( Δ U > 0 ) with low affinity for the lipid heads ( Δ U b > 0 ). The orange circles correspond to the free energy extracted from a control simulation to determine the optimal timestep, as described in the Methods; (c) representative periodic, position-dependent diffusion coefficient D ( z ) corresponding to the system in (a), smoothly changing from D l i p within the lipid phase to D w a t in the water phase far from the lipid/water interface.
Pharmaceutics 15 00573 g001
Figure 2. Representative evolution of mean-square displacement as a function of time for Δ U b = 4 (a) and Δ U b = 8 (b), while Δ U = 0 and D l i p = D w a t . In both panels, filled and empty symbols correspond to M S D and M S D , respectively, as reported in the legends. The dashed lines correspond to the formula 4 D w a t t . In the insets, the potential of mean force for each case is reported.
Figure 2. Representative evolution of mean-square displacement as a function of time for Δ U b = 4 (a) and Δ U b = 8 (b), while Δ U = 0 and D l i p = D w a t . In both panels, filled and empty symbols correspond to M S D and M S D , respectively, as reported in the legends. The dashed lines correspond to the formula 4 D w a t t . In the insets, the potential of mean force for each case is reported.
Pharmaceutics 15 00573 g002
Figure 3. (a) Dependence of long-time diffusion coefficients on the size of the barriers Δ U b , for Δ U = 0 and D l i p = D w a t . The horizontal dashed line corresponds to D w a t = 0.7 nm2/ns. The red continuous line is obtained by fitting the values for Δ U b > 5 k B T via an exponential function. In the inset, we report a representative potential of mean force U ( z ) for this study, corresponding to Δ U b = 4 k B T ; (b) same as (a), but focusing on varying Δ U in systems with Δ U b = 0 and D l i p = D w a t . In the inset, we report U ( z ) for Δ U = 4 k B T .
Figure 3. (a) Dependence of long-time diffusion coefficients on the size of the barriers Δ U b , for Δ U = 0 and D l i p = D w a t . The horizontal dashed line corresponds to D w a t = 0.7 nm2/ns. The red continuous line is obtained by fitting the values for Δ U b > 5 k B T via an exponential function. In the inset, we report a representative potential of mean force U ( z ) for this study, corresponding to Δ U b = 4 k B T ; (b) same as (a), but focusing on varying Δ U in systems with Δ U b = 0 and D l i p = D w a t . In the inset, we report U ( z ) for Δ U = 4 k B T .
Pharmaceutics 15 00573 g003
Figure 4. Effective diffusion coefficients in a system with position-dependent D ( z ) and no external potential. The red continuous line is the prediction from Equation (12). The blue dashed line corresponds to a fit of the simulation data by a power law.
Figure 4. Effective diffusion coefficients in a system with position-dependent D ( z ) and no external potential. The red continuous line is the prediction from Equation (12). The blue dashed line corresponds to a fit of the simulation data by a power law.
Pharmaceutics 15 00573 g004
Figure 5. Dependence of the perpendicular long-time diffusion coefficient D on hydropathy ( Δ U ) and affinity for lipid heads ( Δ U b ), implemented according to the corresponding potential of mean force U ( z ) (Figure 1a); (a) is obtained by assuming D l i p = D w a t , while (b) considers a more realistic diffusion profile with D l i p = 0.09 D w a t (Figure 1c). The contours correspond to the indicated values of D / D w a t .
Figure 5. Dependence of the perpendicular long-time diffusion coefficient D on hydropathy ( Δ U ) and affinity for lipid heads ( Δ U b ), implemented according to the corresponding potential of mean force U ( z ) (Figure 1a); (a) is obtained by assuming D l i p = D w a t , while (b) considers a more realistic diffusion profile with D l i p = 0.09 D w a t (Figure 1c). The contours correspond to the indicated values of D / D w a t .
Pharmaceutics 15 00573 g005
Figure 6. (a) Dependence of the parallel long-time diffusion coefficient D on hydropathy ( Δ U ) and affinity for lipid heads ( Δ U b ), implemented according to the corresponding potential of mean force U ( z ) (Figure 1a). The diffusion profile corresponds to Figure 1c with D l i p = 0.09 D w a t ; (b) dependence of D / D w a t on log P for the simulation data (golden stars) and for an extended set of systems (purple triangles), for which the lateral diffusion coefficient was computed according to Equation (13). The extended systems with log P < 0 were obtained by considering Δ U = Δ U b = 20 k B T , 200 k B T , 2000 k B T ; the extended systems with log P > 0 were obtained by setting Δ U = Δ U b = 10 k B T , 9 k B T , 8 k B T , 7 k B T , 6 k B T . The red continuous line corresponds to Equation (14).
Figure 6. (a) Dependence of the parallel long-time diffusion coefficient D on hydropathy ( Δ U ) and affinity for lipid heads ( Δ U b ), implemented according to the corresponding potential of mean force U ( z ) (Figure 1a). The diffusion profile corresponds to Figure 1c with D l i p = 0.09 D w a t ; (b) dependence of D / D w a t on log P for the simulation data (golden stars) and for an extended set of systems (purple triangles), for which the lateral diffusion coefficient was computed according to Equation (13). The extended systems with log P < 0 were obtained by considering Δ U = Δ U b = 20 k B T , 200 k B T , 2000 k B T ; the extended systems with log P > 0 were obtained by setting Δ U = Δ U b = 10 k B T , 9 k B T , 8 k B T , 7 k B T , 6 k B T . The red continuous line corresponds to Equation (14).
Pharmaceutics 15 00573 g006
Figure 7. (a) Potential of mean force for 16 amino acids as computed in Ref. [63]. The plots were colour-coded according to the physico-chemical properties of the side chain: gold ↔ charged, red ↔ polar and green ↔ apolar; (b) perpendicular diffusion coefficient D for charged (filled gold stars), polar (filled red circles) and apolar (filled green squares) residues as a function of the difference U m a x U m i n between maximum and minimum height of the corresponding potential of mean force. In the inset, the full range of U m a x U m i n is considered to include the case of Arg, for which U m a x U m i n 32 k B T . The dashed grey line is a best fit via an exponential decay; (c) parallel diffusion coefficient D for charged (empty gold stars), polar (empty red circles) and apolar (empty green squares) residues as a function of the logarithm of the partition coefficient log P . The dashed grey line corresponds to the theoretical prediction according to Equation (14).
Figure 7. (a) Potential of mean force for 16 amino acids as computed in Ref. [63]. The plots were colour-coded according to the physico-chemical properties of the side chain: gold ↔ charged, red ↔ polar and green ↔ apolar; (b) perpendicular diffusion coefficient D for charged (filled gold stars), polar (filled red circles) and apolar (filled green squares) residues as a function of the difference U m a x U m i n between maximum and minimum height of the corresponding potential of mean force. In the inset, the full range of U m a x U m i n is considered to include the case of Arg, for which U m a x U m i n 32 k B T . The dashed grey line is a best fit via an exponential decay; (c) parallel diffusion coefficient D for charged (empty gold stars), polar (empty red circles) and apolar (empty green squares) residues as a function of the logarithm of the partition coefficient log P . The dashed grey line corresponds to the theoretical prediction according to Equation (14).
Pharmaceutics 15 00573 g007
Figure 8. Approximation error on D eff when employing Equation (15) for the toy model (a) and the amino-acids simulations (b). The error is computed as 100 · 1 D eff , pred / D eff , sim , where D eff , pred is the predicted value according to Equation (15), and D eff , sim is computed directly from the simulations.
Figure 8. Approximation error on D eff when employing Equation (15) for the toy model (a) and the amino-acids simulations (b). The error is computed as 100 · 1 D eff , pred / D eff , sim , where D eff , pred is the predicted value according to Equation (15), and D eff , sim is computed directly from the simulations.
Pharmaceutics 15 00573 g008
Table 1. Estimated values of effective diffusion coefficient for release of various drugs from a lamellar mesophase at 43 °C with geometric parameters chosen from the literature. The values of D w a t were obtained by renormalizing experimental values obtained at different temperatures via the Stokes–Einstein equation, as discussed in the Methods [57]. Experimental values of D w a t and log P were taken from Refs. [42,54,72,73,74,75,76,77,78,79,80].
Table 1. Estimated values of effective diffusion coefficient for release of various drugs from a lamellar mesophase at 43 °C with geometric parameters chosen from the literature. The values of D w a t were obtained by renormalizing experimental values obtained at different temperatures via the Stokes–Einstein equation, as discussed in the Methods [57]. Experimental values of D w a t and log P were taken from Refs. [42,54,72,73,74,75,76,77,78,79,80].
Name log P D wat (nm2/ns) D eff (nm2/ns)
Cephalexin−0.670.700.27
Hydrochlorothiazide−0.151.690.43
Levodopa0.000.950.21
Piroxicam0.290.850.14
Methyldopa0.391.140.18
Paracetamol0.461.060.15
Antipyrine1.011.040.11
Carbamazepine2.931.130.10
Ketoprofen3.310.670.06
Desipramine3.940.460.04
Ibuprofen3.990.770.07
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bosch, A.M.; Assenza, S. Interplay of Hydropathy and Heterogeneous Diffusion in the Molecular Transport within Lamellar Lipid Mesophases. Pharmaceutics 2023, 15, 573. https://doi.org/10.3390/pharmaceutics15020573

AMA Style

Bosch AM, Assenza S. Interplay of Hydropathy and Heterogeneous Diffusion in the Molecular Transport within Lamellar Lipid Mesophases. Pharmaceutics. 2023; 15(2):573. https://doi.org/10.3390/pharmaceutics15020573

Chicago/Turabian Style

Bosch, Antonio M., and Salvatore Assenza. 2023. "Interplay of Hydropathy and Heterogeneous Diffusion in the Molecular Transport within Lamellar Lipid Mesophases" Pharmaceutics 15, no. 2: 573. https://doi.org/10.3390/pharmaceutics15020573

APA Style

Bosch, A. M., & Assenza, S. (2023). Interplay of Hydropathy and Heterogeneous Diffusion in the Molecular Transport within Lamellar Lipid Mesophases. Pharmaceutics, 15(2), 573. https://doi.org/10.3390/pharmaceutics15020573

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop