Nourddine Hadrioui,
Mohammed Lemaalem*,
Abdelali Derouiche and
Hamid Ridouane
Laboratoire de Physique des Polymères et Phénomènes Critiques Sciences Faculty Ben M’Sik, Hassan II University, P.O. Box 7955, Casablanca, Morocco. E-mail: mohammedlemaalem@gmail.com
First published on 27th February 2020
This work deals with a quantitative investigation of the physical properties of pulmonary surfactant near melting temperature. To this end, we make use of molecular dynamics simulations, using the MARTINI coarse-grained model, for determining the physical properties of the system, such as the potential energy, the specific heat, the microstructure, the diffusion laws, and the elastic properties of the surfactant. The microstructure is studied by computation of the radial-distribution-function upon varying the distance between constituents (lipid molecules or proteins). The diffusion phenomenon is investigated by determination of the mean-squared-displacement and the time dependent velocity-autocorrelation-function for various values of temperature. We show that the dynamics of lipids and proteins exhibit a subdiffusion regime (slow movement) due to the cage effect within pulmonary surfactant. From the obtained mean-squared-displacement, we get the values of the self-diffusion-coefficients and the anomalous exponents at different temperatures close to the melting temperature. For the mathematical description of the cage effect, we make use of the scale relations in terms of the waiting time probability distribution. The last study is concerned with determination of the dependence of the lateral stress upon the strain of pulmonary surfactant, which is found to be linear, and from which we deduce the lateral-elastic-modulus.
Pulmonary surfactant is a complex structure, composed of 90% lipids and 10% proteins, and is strikingly comparable across various species including humans.6 The majority of the surfactant lipids are phospholipids (80–90%), of which phosphatidylcholine (PC) is quantitatively the most important, accounting for 70–80% of the total lipids.7 Four surfactant proteins have so far been identified in pulmonary surfactant,8 namely, SP-A, SP-B, SP-C and SP-D, which have different functions. They are either exclusively lung associated or predominately found in the lung.8–10 The main function of the surfactant is to reduce the surface tension at the air–liquid interface in the alveoli, in a manner that depends on the alveolar surface area. The reduction of surface tension at the very extensive alveolar surface promotes lung expansion on inspiration and prevents lung collapse on expiration. Thus, the surfactant plays a major role in the pressure–volume characteristics of the lung and in gas exchange. However, the presence of the proteins, particularly SP-B and SP-C, is absolutely required for interfacial adsorption, film stability, and respreading capacities.11,12 SP-B and SP-C are small hydrophobic polypeptides that play crucial roles in the formation and stabilization of pulmonary surfactant films.13 In contrast to SP-B, the surfactant protein SP-C is not absolutely required for survival. However, targeted deletion of the SP-C gene in mice leads to the development of severe progressive pneumonitis depending on the genetic background.14,15 Experimental data suggest that SP-C would be required to maintain the association of the surfactant complexes (reservoirs) with the interfacial film in the most compressed states, i.e., those reached at the end of expiration.9 The SP-C-promoted attachment would then facilitate the reinsertion of the surface-active molecules from the reservoirs, with the probable critical participation of SP-B during the re-expansion.11,16–19 Some authors have suggested18,20 the existence of SP-C/cholesterol complexes within the bilayer membranes; however, both direct evidence and a model explaining the role of SP-C/cholesterol interactions in surfactant function are lacking. It has been proposed that cholesterol could contribute in preventing SP-C from aggregating, leading to better miscibility of the protein into the monolayers.20 On the other hand, a previous study suggested that concomitant changes in the concentrations of cholesterol and SP-C may occur in the surfactant of heterothermic animals, and therefore, the roles of these two molecules in the surfactant may be connected.21
The present study is motivated by the interest in pulmonary surfactant in many disciplines such as biological research and as an exogenous surfactant treatment to neonates in medicine. The aim is to rely on the physical properties of the pulmonary surfactant components to understand its physiological roles.
In recent years, molecular simulations have become an excellent tool for studying the physicochemical behavior of biological membranes, especially with the advancement in computational tools.22–24 Numerous studies in relation to pulmonary surfactant have been accomplished using simulation methods.
In this context, Wang and Frenkel have proposed25 a simple model of the amphiphilic bilayer, in order to reproduce the elastic properties, using Monte Carlo simulations. They have shown that such a model gives a good description of the elastic properties and agrees with experiments. Also, they have observed that the solid phase exists only at low temperature, whereas at ambient and high temperatures, the bilayer membrane is in the liquid state. Also, Freites et al. have performed26 an all-atom MD simulation on a model system, composed of a monolayer of palmitic acid and a surfactant protein B peptide SP-B. They have shown that the palmitic acid/SP-B monolayer is stable in a wide range of lipid specific areas in regions of the monolayer containing tilted and untilted condensed phases. S. Baoukina and D. P. Tieleman27 presented a review of recent simulation studies on the properties of lipid–protein interactions, and the effects of nanoparticles on lung surfactant. In another work, Q. Hu et al.28 investigated the interaction between nanoparticles and the pulmonary surfactant system in alveolar fluid, using coarse-grain molecular dynamics simulations. Their simulation results revealed that, when brought into contact with pulmonary surfactant, silver and polystyrene nanoparticles are immediately covered with a biomolecular corona composed of both lipids and proteins. In another recently published paper, K. Dutta et al.29 studied the paradoxical bactericidal effects of hydrophobic lung surfactant proteins and their peptide mimics using a molecular Trojan horse liposome. This study has alarming implications in using high dosage exogenous surfactant for treatment of respiratory distress syndrome, genetic knock-out abnormalities associated with these proteins, and the novel roles played by SP-B/C as bactericidal agents.
In the present work, we are interested in the study of several physical aspects of a pulmonary surfactant system near melting temperature. For the sake of simplification, we proposed a system composed of SP-B and SP-C integral proteins, and a phospholipid bilayer, with water as an explicit solvent at the interfaces. To determine the physical properties of the surfactant, we used MD simulations using the MARTINI coarse-grained force-field. These simulations allow us to determine the melting temperature, structure (through the radial-distribution-function), dynamics (subdiffusion laws via the mean-squared-displacement and velocity-autocorrelation-function), and mechanical properties of pulmonary surfactant by the lateral-elastic-modulus.
Our results show the importance of the phospholipids and integral proteins in the physiological role of pulmonary surfactant. The analysis of the simulation results suggests that, at human body temperature, the surfactant is in the gel phase, due to the presence of congested proteins and locally condensed lipids. These constituents are also characterized by slow diffusion, called subdiffusion, and significant stiffness against stretching and compression during the breathing cycle.
The proteins which exist in the structure of pulmonary surfactant are SP-A, SP-B, SP-C, and SP-D. These proteins can be divided into two groups: small hydrophobic proteins (SP-B and SP-C) and larger hydrophilic proteins (SP-A and SP-D).
Many experimental studies have been conducted to study the three-dimensional structure of SP-B as an atomic representation,39,40 giving rise to an overview of the SP-B conformation and disposition in pulmonary surfactant. According to these studies, SP-B is a hydrophobic protein consisting of 79 amino acids, with molecular weight of about 9 kDa. It is in the form of dimers and trimers (respectively two and three repetitive units). Because of its hydrophobicity, it will preferentially interact with the lipid-tails, as has been shown using the electron paramagnetic resonance technique.11 One of the most important roles of SP-B is the stabilization of lipid structures at the water–liquid interface. Thus, the presence of SP-B consolidates the lipid structures, and at the same time, improves their mechanical properties during the cyclic compression-expansion movements of the respiratory system. SP-C consists of 35 amino acids, and has a molecular weight between 3 and 5 kDa. Its three dimensional structure was determined by Johansson et al.41 using the nuclear magnetic resonance technique.
The structure showed a valyl-rich α-helix formed by amino acids. It was also well matched with the analysis of the secondary structure by circular dichroism and Fourier transform infrared spectroscopy. During the respiratory cycle, SP-C is required to maintain the combination of the pulmonary surfactant with the interfacial film at the end of exhalation. Then, SP-C facilitates the re-insertion of the surfactant molecules into the reservoir during re-expansion.16–19 The integral proteins SP-B and SP-C account for approximately 2% of the total surfactant weight.
SP-A and SP-D, which constitute approximately 6% of the total surfactant-weight,13 are surfactant-associated water-soluble proteins. These proteins play an important role in host defense but have little effect on the biophysical properties of the surfactant.42 We note that pulmonary surfactant can contain other proteins accounting for approximately 2% of the total surfactant weight.
In this study, we are interested in the biophysical properties of pulmonary surfactant. Thus, the proposed model for pulmonary surfactant used in this study consists only of phospholipids and integral proteins (SP-B and SP-C).
The two hydrophobic surfactant proteins, SP-B and SP-C, were constructed using their all-atom models obtained from the Protein Data Bank. The all-atom models were then transferred to the corresponding coarse-grained models with the martinize.py script provided by the MARTINI force field. This transformation gives a coarse-grained representation similar to the α-helix model. This model is adopted in several recent reports.13,27–29 In the MARTINI model, each amino acid is represented by a single bead that can be polar (sL), non-polar (sH) or neutral (tN). The sL and sH beads are mutually repulsive and differ only in the parameterization of their respective dihedral potentials. The coarse-grained protein scheme has also been well discussed in the literature by Monticelli et al.44 and in other work by Bellesia et al.45
(1) |
Fig. 1(a) illustrates pulmonary surfactant in the gel phase, for T = 293 K, which is below the melting temperature (about 320 K). But, above this temperature, the bilayer membrane begins to fluctuate. Beyond the melting temperature, as shown in Fig. 1(b), at T = 333 K, the bilayer membrane is in the liquid phase. Therefore, the transition from the gel phase to the liquid phase occurs at the melting temperature, TM.
Fig. 1 The simulated pulmonary surfactant model in an orthogonal simulation box of lengths lx = ly = 425 Å and lz = 100 Å, at two different temperatures: (a) T = 293 K and (b) T = 333 K. The structure was generated using the Moltemplate package (https://moltemplate.org/),24 and visualized using the Ovito software (https://ovito.org).23 The lipid tails C1 are represented in white, the different heads of each lipid are presented in yellow (Na), purple (Qa) and blue (Q0, P4), the water beads are presented in turquoise, and the protein beads are presented in black (sH), cyan (sL) and green (sN). More detailed information about the MARTINI coarse-grained representation and parametrization of the interaction potentials of lipids, proteins (amino acids) and water can be found on the MARTINI home page (http://cgmartini.nl/). |
In Fig. 2, we report the variation of the specific heat (solid line + triangles in blue) and potential energy (solid line + spheres in black) with temperature, calculated using a heating rate of 1011 K s−1. At fixed surfactant area, we find that the potential energy increases with increasing temperature, and its variation becomes steeper beyond the melting temperature. In particular, the energy exhibits a jump at this temperature, which is a signature of a first-order phase transition. The specific heat exhibits a peak at the melting temperature, which equals 320 K, marking the possibility of the appearance of a solid phase at physiological temperatures. Several experimental studies have been carried out on a pure DPPC bilayer and revealed that TM = 314 K.51,52 We note that the insertion of the proteins leads to an increase of the melting temperature, TM, since their presence has a tendency to rigidify the bilayer membrane. We expect that TM increases with the protein fraction, Φ, as: TM = T0M + AΦx, with exponent x > 1. Here, T0M represents the melting temperature relative to surfactant free from proteins.
Fig. 2 Potential energy (Ep) and specific heat (CV) of pulmonary surfactant versus temperature, for a heating rate of 1011 K s−1. |
From the literature, Yan et al. conducted an experimental study to determine the melting temperature of pulmonary surfactant.53 Thus, the authors found that TM is significantly higher in the lungs than in the lipid bilayer. Their results indicate that the change in pulmonary mechanics should occur in the range of 321–328 K.53 Our results are in good agreement with the experimental results discussed above. In addition, the present study provides a clear overview of other physical properties at temperatures near the melting point, including physiological temperatures. Also, it should be noted that the physical properties of pulmonary surfactant are still not completely understood and are different from those of simple lipid bilayers.
(2) |
Fig. 4 (a) RDF of the group of proteins, and (b) RDF of the group of pulmonary surfactant constituents (phospholipid molecules and proteins). |
Fig. 3(a) and (b) indicate that the phospholipid molecules exist in a condensed liquid phase, which is characterized by the presence of short-range order (for the closest neighbors). This liquid behavior has been shown experimentally54 for long-range lipids (heads and tails). Fig. 4(a) suggests that the proteins exhibit a congested structure, close to a solid, which manifests itself in the form of intense and narrow peaks (for the first and second closest neighbors) in the RDF. The position of the first peak is the average distance between the lipid head groups (about 5 Å), which are inside pores formed by the integral proteins, while the second peak accounts for the pore radius (about 9.5 Å), in agreement with the literature.55,56
Fig. 4(b) shows the RDF of the surfactant molecule group (phospholipid molecules and proteins). This figure reveals that the structure of pulmonary surfactant is almost identical to the lipid structure. This means that the inclusion of the proteins, with a small percentage of (10%), has a minor influence on the static structure of the surfactant. The computation of RDFs for temperatures close to the melting point suggests that there is no remarkable influence on the surfactant structure. These findings may be useful for theorists who are interested in the determination of the effective interactions between the proteins which are included in the bilayer membrane, where the latter may be considered a two dimensional continuum fluid. They can, therefore, consider the proteins as nanoparticles in a two dimensional solvent, but this ignores the fact that the proteins have slow dynamics in comparison with lipids, and a solid-like structure.
(3) |
In the expression above, rn(t) represents the two dimensional time position of a constituent n, which is a protein or a phospholipid molecule. Here, t denotes time and t0 is the initial time at which the tracer (random walker or target) starts moving. In our case, the tracer is a phospholipid molecule or a protein.
Due to thermal agitation, the proteins and phospholipids exhibit lateral diffusion. Their motion is hampered by the presence of other neighboring entities. As a consequence, a given tracer (a protein or a phospholipid molecule) exhibits an anomalous diffusion that we are interested in. In comparison to normal or Brownian diffusion, the anomalous diffusion is slower, due to the extreme difficulty with which a particle moves in a complex structure. This kind of diffusion is usually referred to as subdiffusion, and it is characterized by an MSD which behaves as
W(t) = 2dDαtα, 0 < α < 1. | (4) |
logW(t) = log(2dDα) + αlogt. | (5) |
This equality means that the anomalous exponent α is simply the slope of the logarithm of MSD, logW, as a function of logt, and the generalized-diffusion-coefficient, Dα, can be determined by extrapolation to the vertical axis.
We emphasize that subdiffusion is a feature of crowded systems, where the trajectories of the mobile constituents are strongly correlated. Notice that the above scaling relation is valid for long times, which are beyond some characteristic times depending on the specific details of the diffusion process and the structure of the host medium. In general, a particle is said to be subdiffusive if the condition W(t)/t → 0, for t → +∞, is satisfied (very slow diffusion). This explains why the anomalous exponent α must be in the interval 0 < α < 1.
Another interesting dynamic quantity to consider is the velocity-autocorrelation-function (VACF), denoted as cvv(t). The latter is simply the time-average: cvv(t) = 〈v(t) × v(0)〉, where v(t) is the velocity of the tracer at time t. We recall that the MSD and VACF are connected by67
(6) |
According to the asymptotic time behavior in eqn (4), the VACF then scales as
cvv(t) = 2Dαα(α − 1)tα−2, 0 < α < 1, | (7) |
Quantitatively, the cage effect can be studied through the waiting time probability,69 denoted as P(τ). For the reasoning, we can assimilate the membrane surface to a regular two dimensional lattice on which the tracer moves at random and may diffuse from time to time, with Nc mobile traps forming a cage (a surface on the bilayer membrane) that confines this tracer. The latter has to wait for a time τ on each site before performing the next jump. In fact, this waiting time is a random variable independently chosen for each new jump according to a distribution, P(τ). The problem dealt with here is called a “continuous time random walk” in the literature.70
After a characteristic long time, tc, the tracer leaves its cage. Physically, this time is simply the stay duration in the cage formed by the neighbors. It is related to the mean waiting time, 〈τ〉, calculated with the probability, P(τ), by the following simple formula70
tc = Nc〈τ〉. | (8) |
When the distribution, P(τ), is board, it behaves as70
P(τ) ∼ τ−(1+α), τ < tc. | (9) |
At long times, that is for τ > tc, this distribution falls exponentially. Here, 0 < α < 1 is the anomalous exponent, characterizing the long-time behavior of the MSD, according to relation (4). Using relation (9), we find that the mean waiting time is given by
(9a) |
tc ∼ τ0N1/αc. | (10) |
It is noted that both the tracer and the mobile traps diffuse normally among themselves and anomalously with the molecules of the host liquid (the bilayer membrane).
We emphasize that the above scaling-relation can be rapidly recovered by equating the linear-size of the confinement space (cage), aN1/2c, to the square-root of the MSD, taken at t = tc. Using eqn (4) for the MSD yields formula (10).
Now, if the mobile traps form an aggregate of fractal dimension, dF, we find that, in this case, the stay duration scales with Nc as
tc ∼ τ0N2/αdcF. | (10a) |
Notice that, in order to get this scaling-relation, we have equated the linear size of the confinement space (aggregate), aN1/dcF, to the square-root of the MSD, at t = tc. The above behavior may be applied to proteins with high density, where these are aggregated due to the attractive van der Waals forces between them. Finally, notice that, if the tracer diffuses with the surrounding traps in Euclidean space, we would have dF = 2, and then, the scaling-relation (10) is recovered.
The next step consists of determining the dynamic characteristics of the constituents (lipid molecules or proteins) within the bilayer membranes, using MD simulations.
(11) |
Fig. 5 Lateral mean-squared-displacements of the phospholipid molecules versus time: (a) MSDs and (b) log–log plots of MSDs. |
Fig. 6 Lateral mean-squared-displacements of the proteins with time: (a) MSDs and (b) log–log plots of MSDs. |
Of course, this characteristic time is temperature dependent. Its temperature dependence will be discussed below. Tables 1 and 2 summarize the computed anomalous exponents and the associated generalized-diffusion-coefficients for the lipids and proteins for various values of temperature (near the melting point) and time-regimes R1 and R2. These tables clearly show that the anomalous exponent does not depend on temperature (around the melting point), while the generalized-diffusion-coefficient does. But, the values of these two quantities are naturally different for the time-regimes R1 and R2. We notice that the generalized-diffusion-coefficient increases with increasing temperature. Since this coefficient is essentially the inverse of the medium viscosity, we expect that this coefficient follows the following Arrhenius-like law:
Dα(T) = Aαe−Bα/T. | (12) |
log(Dα(T)) = log(Aα) − Bα/T, | (12a) |
Aα1 ≃ 8.21 × 10−7, Bα1 ≃ 289.94 K, (R1), | (12b) |
Aα2 ≃ 36.5 × 10−7, Bα2 ≃ 631.43 K, (R2) | (12c) |
Aα1 ≃ 10.00 × 10−8, Bα1 ≃ 256.88 K, (R1), | (12c) |
Aα2 ≃ 130.40 × 10−8, Bα2 ≃ 759.12 K, (R2) | (12d) |
Regime | T (K) | α | Dα (cm2/tα) |
---|---|---|---|
R1 | 293 | 0.37 | 3.07 × 10−7 |
303 | 0.37 | 3.14 × 10−7 | |
313 | 0.37 | 3.23 × 10−7 | |
323 | 0.37 | 3.35 × 10−7 | |
333 | 0.37 | 3.45 × 10−7 | |
R2 | 293 | 0.53 | 4.24 × 10−7 |
303 | 0.53 | 4.5 × 10−7 | |
313 | 0.53 | 4.85 × 10−7 | |
323 | 0.53 | 5.2 × 10−7 | |
333 | 0.53 | 5.45 × 10−7 |
Regime | T (K) | α | Dα (cm2/tα) |
---|---|---|---|
R1 | 293 | 0.19 | 4.18 × 10−8 |
303 | 0.19 | 4.31 × 10−8 | |
313 | 0.19 | 4.44 × 10−8 | |
323 | 0.19 | 4.53 × 10−8 | |
333 | 0.19 | 4.65 × 10−8 | |
R2 | 293 | 0.65 | 9.50 × 10−8 |
303 | 0.65 | 1.07 × 10−7 | |
313 | 0.65 | 1.16 × 10−7 | |
323 | 0.65 | 1.22 × 10−7 | |
333 | 0.65 | 1.31 × 10−7 |
In a recent series of reports, Matti Javanainen, Ralf Metzler and their coworkers71–77 studied the anomalous dynamics and subdiffusion in protein–lipid mixture bilayers, using all-atom and coarse-grained MD simulations and single-particle tracking experiments. For a pure lipid bilayer with different molar fractions of incorporated proteins, the authors observed anomalous diffusion of lipids with a scaling exponent α ≃ 0.8 below a crossover time of τc ≃ 10 ns,71,72 which is a long subdiffusion time. Also, they found that the calculated diffusion coefficient of lipids is of the order of 10−7 cm2 s−1 and decreases proportionally to the protein molar fraction. In this work, we found that the scaling exponent for the phospholipid mixture in pulmonary surfactant is α ≃ 0.53. This value is slightly lower than the value found in the work of R. Metzler et al. (α ≃ 0.8), and this can be explained by the fact that the variety of lipids with different tail lengths and the proteins included in the surfactant may form cage-like environments.78 Thus, this heterogeneity causes slower and more pronounced anomalous diffusion of lipids in pulmonary surfactant. The calculated diffusion coefficient of lipids is also found to be in good agreement with previous MD simulations and experimental work.71,72,79
According to Matti Javanainen et al.,71,72 the protein–lipid mixture in the crowded case corresponding to a molar fraction of 1:50, which is close to the pulmonary surfactant protein–lipid molar fraction (10:100), is characterized by slow diffusion with a scaling exponent α ≃ 0.65 and a diffusion coefficient of the order of 10−8 cm2 s−1. From the present study, the subdiffusion exponent of the proteins in pulmonary surfactant was found to be α = 0.65. These results are in the order of those found by Matti Javanainen et al.71,72 Also, the obtained diffusion coefficient of proteins, which is in the order of 10−8 cm2 s−1, agrees with the experimental findings.80–82 Thus, anomalous diffusion is an inherent property of both lipids and proteins in pulmonary surfactant. The time scales of the anomalous diffusion coefficients depend on temperature, following an Arrhenius law, whereas the anomalous diffusion exponents do not. It should be noted that the generalized-diffusion-coefficients, Dα, of proteins are less than those of phospholipid molecules. Thus, the proteins diffuse slowly compared to the lipid molecules in pulmonary surfactant.
According to Tables 1 and 2, invoked above, the crossover-time, τ*, decreases on increasing the temperature toward the melting point. Also, it follows the Arrhenius law:
(13) |
Since α2 > α1 and Bα2 > Bα1, τ*(T) is a decreasing function of temperature. Explicitly, we have the closer forms
τ*(T) = 1.91 × 10−4 × e1924.20/T, phospholipid molecules, | (13a) |
τ*(T) = 7.70 × 10−3 × e943.44/T, proteins. | (13b) |
Table 3 summarizes the crossover-times, τ*, for the lipid molecules and proteins, for three values of temperature close to the melting point. This table indicates that the second time-regime, R2, is reached later by the proteins, in comparison with the lipids.
Molecules | T (K) | τ* (ns) |
---|---|---|
Phospholipid molecules | 293 | 0.1358 |
313 | 0.0893 | |
333 | 0.0617 | |
Proteins | 293 | 0.1927 |
313 | 0.1568 | |
333 | 0.1308 |
From an initial configuration balanced in the NPT-ensemble, with P = 1 bar and T = 300 K, we perform a uniaxial tensile strain MD simulation. We find that the stress increases linearly with the strain (Fig. 9), and a linear fit to the data gives the elastic-constant Ka = 10 N m−2. This computed value agrees with previous experimental work on bilayer membranes.87,88
In Fig. 9, we show the stress–strain plot of pulmonary surfactant, with a strain rate of 10−5 fs−1. During respiratory cycles, the elasticity-modulus describes qualitatively the effectiveness of pulmonary surfactant, that is to say, its ability to oppose the decrease of the volume of the pulmonary cavity. In fact, this decreases the work of the respiratory muscles and then prevents alveolar collapse. The presence of the lipids in the gel phase induces a phase separation, which is translated by the presence of the condensed domains. Pulmonary surfactant contains two phases (condensed and expanded) and has an elasticity-modulus between those of the condensed and expanded phases. Indeed, the presence of the lipids, which fluidify and thus reduce the elasticity-modulus of the pulmonary surfactant in the pulmonary alveoli compared to the DPPC monolayer, is a compromise between the rigidity necessary for good efficiency and the fluidity necessary for good adsorption at the air–liquid interface. These two conditions are a prerequisite for the proper functioning of pulmonary surfactant under certain dynamic conditions (respiratory cycles).
This journal is © The Royal Society of Chemistry 2020 |