Computer simulations of lung surfactant

A Review on Investigating the Interactions between Nanoparticles and the Pulmonary Surfactant Monolayer with Coarse-Grained Molecular Dynamics Method

Pulmonary drug delivery has garnered significant attention due to its targeted local lung action, minimal toxic side effects, and high drug utilization. However, the physicochemical properties of inhaled nanoparticles (NPs) used as drug carriers can influence their interactions with the pulmonary surfactant (PS) monolayer, potentially altering the fate of the NPs and impairing the biophysical function of the PS monolayer. Thus, the objective of this review is to summarize how the physicochemical properties of NPs affect their interactions with the PS monolayer. Initially, the definition and properties of NPs, as well as the composition and characteristics of the PS monolayer, are introduced. Subsequently, the coarse-grained molecular dynamics (CGMD) simulation method for studying the interactions between NPs and the PS monolayer is presented. Finally, the implications of the hydrophobicity, size, shape, surface charge, surface modification, and aggregation of NPs on their interactions with the PS monolayer and on the composition of biomolecular corona are discussed. In conclusion, gaining a deeper understanding of the effects of the physicochemical properties of NPs on their interactions with the PS monolayer will contribute to the development of safer and more effective nanomedicines for pulmonary drug delivery.

Analyzing Lipid Membrane Defects via a Coarse-Grained to Triangulated Surface Map: The Role of Lipid Order and Local Curvature in Molecular Binding

Lipid packing defects are known to serve as quantitative indicators for protein binding to lipid membranes. This paper presents a protocol for mapping molecular lipid detail onto a triangulated continuum leaflet representation. Besides establishing the desired forward counterpart to the existing inverse TS2CG map, this coarse-grained to triangulated surface (CG2TS) map enables straightforward extraction of the defect characteristics for any membrane geometry found in nature. We have applied our protocol to investigate the role of local curvature and varying lipid packing on the defect constant π. We find that the defect size is greatly influenced by both factors, arguing strongly against the usual assignment of a single defect constant in the case of more realistic membrane conditions. An important discovery is that lipids in the gel phase produce larger defects, or a higher π, in domains of high (local) curvature than the same lipid in a liquid phase of any curvature. This finding suggests that membranes featuring very ordered lipid packing can bind proteins via large defects in curved regions. Finally, we propose a route for estimating defect constants directly from the standard membrane properties. Identifying the precise role of composition, lipid (tail) order, and (local) curvature in defects for the irregular lipid structures that are (temporally) present in many biological processes is instrumental for obtaining fundamental insight as well as for a rational design of membrane binding targets.

Pulmonary Surfactant: A Mighty Thin Film

Pulmonary surfactant is a critical component of lung function in healthy individuals. It functions in part by lowering surface tension in the alveoli, thereby allowing for breathing with minimal effort. The prevailing thinking is that low surface tension is attained by a compression-driven squeeze-out of unsaturated phospholipids during exhalation, forming a film enriched in saturated phospholipids that achieves surface tensions close to zero. A thorough review of past and recent literature suggests that the compression-driven squeeze-out mechanism may be erroneous. Here, we posit that a surfactant film enriched in saturated lipids is formed shortly after birth by an adsorption-driven sorting process and that its composition does not change during normal breathing. We provide biophysical evidence for the rapid formation of an enriched film at high surfactant concentrations, facilitated by adsorption structures containing hydrophobic surfactant proteins. We examine biophysical evidence for and against the compression-driven squeeze-out mechanism and propose a new model for surfactant function. The proposed model is tested against existing physiological and pathophysiological evidence in neonatal and adult lungs, leading to ideas for biophysical research, that should be addressed to establish the physiological relevance of this new perspective on the function of the mighty thin film that surfactant provides.

Adsorption of pulmonary and exogeneous surfactants on SARS-CoV-2 spike protein

COVID-19 is transmitted by airborne particles containing virions of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Coronavirus virions represent nanoparticles enveloped by a lipid bilayer decorated by a "crown" of Spike protein protrusions. Virus transmission into the cells is induced by binding of Spike proteins with ACE2 receptors of alveolar epithelial cells. Active clinical search is ongoing for exogenous surfactants and biologically active chemicals capable of hindering virion-receptor binding. Here, we explore by using coarse-grained molecular dynamics simulations the physico-chemical mechanisms of adsorption of selected pulmonary surfactants, zwitterionic dipalmitoyl phosphatidyl choline and cholesterol, and exogeneous anionic surfactant, sodium dodecyl sulfate, on the S1-domain of the Spike protein. We show that surfactants form micellar aggregates that selectively adhere to the specific regions of the S1-domain that are responsible for binding with ACE2 receptors. We find distinctly higher cholesterol adsorption and stronger cholesterol-S1 interactions in comparison with other surfactants, that is consistent with the experimental observations of the effects of cholesterol on COVID-19 infection. Distribution of adsorbed surfactant along the protein residue chain is highly specific and inhomogeneous with preferential adsorption around specific amino acid sequences. We observe preferential adsorption of surfactants on cationic arginine and lysine residues in the receptor-binding domain (RBD) that play an important role in ACE2 binding and are present in higher amounts in Delta and Omicron variants, which may lead to blocking direct Spike-ACE2 interactions. Our findings of strong selective adhesion of surfactant aggregates to Spike proteins have important implications for informing clinical search for therapeutic surfactants for curing and preventing COVID-19 caused by SARS-CoV-2 and its variants.

Surfactant Proteins SP-B and SP-C in Pulmonary Surfactant Monolayers: Physical Properties Controlled by Specific Protein-Lipid Interactions

The lining of the alveoli is covered by pulmonary surfactant, a complex mixture of surface-active lipids and proteins that enables efficient gas exchange between inhaled air and the circulation. Despite decades of advancements in the study of the pulmonary surfactant, the molecular scale behavior of the surfactant and the inherent role of the number of different lipids and proteins in surfactant behavior are not fully understood. The most important proteins in this complex system are the surfactant proteins SP-B and SP-C. Given this, in this work we performed nonequilibrium all-atom molecular dynamics simulations to study the interplay of SP-B and SP-C with multicomponent lipid monolayers mimicking the pulmonary surfactant in composition. The simulations were complemented by z-scan fluorescence correlation spectroscopy and atomic force microscopy measurements. Our state-of-the-art simulation model reproduces experimental pressure-area isotherms and lateral diffusion coefficients. In agreement with previous research, the inclusion of either SP-B and SP-C increases surface pressure, and our simulations provide a molecular scale explanation for this effect: The proteins display preferential lipid interactions with phosphatidylglycerol, they reside predominantly in the lipid acyl chain region, and they partition into the liquid expanded phase or even induce it in an otherwise packed monolayer. The latter effect is also visible in our atomic force microscopy images. The research done contributes to a better understanding of the roles of specific lipids and proteins in surfactant function, thus helping to develop better synthetic products for surfactant replacement therapy used in the treatment of many fatal lung-related injuries and diseases.

Effect of Fluoxetine on the Surface Behavior of the Lipid Monolayers at Different Surface Pressures

Fluoxetine (FLX), used in the clinic to treat depression, is a well-known cationic amphiphilic antidepressant. However, there is a lack of research on the effect of FLX on the surface behavior of lipid monolayers under different surface pressures. In this study, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/CHOL (DPPC/POPC/CHOL) monolayers were prepared via the Langmuir method, and FLX was added to these monolayers under various surface pressures. The effect of FLX on the surface behavior of DPPC/POPC/CHOL monolayers under various surface pressures was studied using a combination of surface pressure-area isotherms, compressibility modulus-surface pressure curves, and atomic force microscope (AFM). The results showed that the effect of FLX on the lipid monolayers was different under different surface pressures. The interaction between FLX and lipid molecules was weak under low surface pressures, and FLX could easily intercalate between the lipid molecules to inhibit monolayer phase transition. The interaction between FLX and lipid molecules was enhanced and FLX tended to self-aggregate to reduce the monolayer stability when the surface pressure was high. This study lays the foundation for further studies on the interaction between FLX and lipid monolayers.

A Core-Shell Approach for Systematically Coarsening Nanoparticle-Membrane Interactions: Application to Silver Nanoparticles

The continuous release of engineered nanomaterial (ENM) into the environment may bring about health concerns following human exposure. One important source of ENMs are silver nanoparticles (NPs) that are extensively used as anti-bacterial additives. The introduction of ENMs into the human body can occur via ingestion, skin uptake or the respiratory system. Therefore, evaluating how NPs translocate over bio-membranes is essential in assessing their primary toxicity. Unfortunately, data regarding membrane-NP interaction is still scarce, as is theoretical and in silico insight into what governs adhesion and translocation for the most relevant NPs and membranes. Coarse-grained (CG) molecular descriptions have the potential to alleviate this situation, but are hampered by the absence of a direct link to NP materials and membrane adhesion mechanisms. Here, we interrogate the relationship between the most common NP representation at the CG level and the adhesion characteristics of a model lung membrane. We find that this representation for silver NPs is non-transferable, meaning that a proper CG representation for one size is not suited for other sizes. We also identify two basic types of primary adhesion-(partial) NPs wrapping by the membrane and NP insertion into the membrane-that closely relate to the overall NP hydrophobicity and significantly differ in terms of lipid coatings. The proven non-transferability of the standard CG representation with size forms an inspiration for introducing a core-shell model even for bare NPs that are uniform in composition. Using existing all-atom molecular dynamics (MD) data as a reference, we show that this extension does allow us to reproduce size-dependent NP adhesion properties and lipid responses to NP binding at the CG level. The subsequent CGMD evaluation for 10 nm Ag NPs provides new insight into membrane binding for relevant NP sizes and into the role of water in trapping NPs into defected mixed monolayer-bilayer states. This development will be instrumental for simulating NP-membrane adhesion towards more experimentally relevant length and time scales for particular NP materials.

Potential health risks of the interaction of microplastics and lung surfactant

Microplastics (MPs), as pollutants of environmental concern, are correlated with increased risk of various respiratory diseases. Nevertheless, whether or not MPs have adverse influences on the interfacial properties of lung surfactant (LS), and its effect on the generation of reactive oxygen species are poorly understood. In the present study, natural LS extracted from porcine lungs was used to investigate the interaction with polystyrene as a representative MPs. The results showed that the phase behavior, surface tension, and membrane structure of the LS were altered in the presence of polystyrene. Adsorption experiments demonstrated that in the mixed system of polystyrene and LS (the main active ingredients are phospholipids and proteins), adsorption of phospholipid components by polystyrene was notably higher than that of proteins. Moreover, polystyrene can accelerate the conversion between ascorbic acid and deoxyascorbic acid, thereby producing hydrogen peroxide (HOOH) in simulated lung fluid (containing LS) and further giving rise to an increase in the content of hydroxyl radicals (•OH). This work provides new insight into the potential hazard of MPs in human respiratory system, which is helpful for deeply understanding the unfavorable physicochemical effects of MPs exposure and the role of inhaled MPs on lung health.

Accurate Simulations of Lipid Monolayers Require a Water Model with Correct Surface Tension

Lipid monolayers provide our lungs and eyes their functionality and serve as proxy systems in biomembrane research. Therefore, lipid monolayers have been studied intensively including using molecular dynamics simulations, which are able to probe their lateral structure and interactions with, e.g., pharmaceuticals or nanoparticles. However, such simulations have struggled in describing the forces at the air-water interface. Particularly, the surface tension of water and long-range van der Waals interactions have been considered critical, but their importance in monolayer simulations has been evaluated only separately. Here, we combine the recent C36/LJ-PME lipid force field that includes long-range van der Waals forces with water models that reproduce experimental surface tensions to elucidate the importance of these contributions in monolayer simulations. Our results suggest that a water model with correct surface tension is necessary to reproduce experimental surface pressure-area isotherms and monolayer phase behavior. The latter includes the liquid expanded and liquid condensed phases, their coexistence, and the opening of pores at the correct area per lipid upon expansion. Despite these improvements of the C36/LJ-PME with certain water models, the standard cutoff-based CHARMM36 lipid model with the 4-point OPC water model still provides the best agreement with experiments. Our results emphasize the importance of using high-quality water models in applications and parameter development in molecular dynamics simulations of biomolecules.

Drug Meets Monolayer: Understanding the Interactions of Sterol Drugs with Models of the Lung Surfactant Monolayer Using Molecular Dynamics Simulations

The lung surfactant monolayer (LSM) is a thin layer of lipids and proteins that forms the air/water interface of the alveoli. The primary function of the LSM is to reduce the surface tension at the air/water interface during breathing. The LSM also forms the main biological barrier for any inhaled particles, including drugs, to treat lung diseases. Elucidating the mechanism by which these drugs bind to and absorb into the LSM requires a molecular-level understanding of any drug-induced changes to the morphology, structure, and phase changes of the LSM.Molecular dynamics simulations have been used extensively to study the structure and dynamics of the LSM. The monolayer is usually simulated in at least two states: the compressed state, mimicking exhalation, and the expanded state, mimicking inhalation. In this chapter, we provide detailed instructions on how to set up, run, and analyze coarse-grained MD simulations to study the concentration-dependent effect of a sterol drug on the LSM, both in the expanded and compressed state.

The role of size and nature in nanoparticle binding to a model lung membrane: an atomistic study

Understanding the uptake of nanoparticles (NPs) by different types of cellular membranes plays a pivotal role in the design of NPs for medical applications and in avoiding adverse effects that result in nanotoxicity. Yet, the role of key design parameters, such as the bare NP material, NP size and surface reactivity, and the nature of NP coatings, in membrane remodelling and uptake mechanisms is still very poorly understood, particularly towards the lower range of NP dimensions that are beyond the experimental imaging resolution. The same can be said about the role of a particular membrane composition. Here, we systematically employ biased and unbiased molecular dynamics simulations to calculate the binding energy for three bare materials (Ag/SiO₂/TiO₂) and three NP sizes (1/3/5 nm diameter) with a representative lung surfactant membrane, and to study their binding kinetics. The calculated binding energies show that irrespective of size, Ag nanoparticles bind very strongly to the bilayer, while the NPs made of SiO₂ or TiO₂ experience very low to no binding. The unbiased simulations provide insight into how the NPs and membrane affect each other in terms of the solvent-accessible surface area (SASA) of the NPs and the defect types and fluidity of the membrane. Using these systematic fine-grained results in coarsening procedures will pave the way for simulations considering NP sizes that are well beyond the membrane thickness, i.e. closer to experimental dimensions, for which different binding characteristics and more significant membrane remodelling are expected.

The steroid mometasone alters protein containing lung surfactant monolayers in a concentration-dependent manner

Mometasone is an investigational anti-inflammatory steroidal drug to treat inflammation via pulmonary administration. For steroid drugs to be effective they need to be adsorbed by lung surfactants, a thin monolayer at the air-water interface in alveoli that reduces surface tension. Information on the molecular-level interactions of the drug with lung surfactants is useful to understand the mechanism of adsorption. In this study, we use coarse-grained molecular dynamics simulation to understand the concentration-dependent effect of mometasone on a lung surfactant monolayer (LSM) composed of lipids and surfactant proteins, under two different breathing conditions (exhalation, at surface tension 0 mNm^-1 and inhalation, surface tension 20-25 mNm^-1). A series of fixed-APL and fixed-surface tension simulations were used to demonstrate that in the absence of drugs, the model LSM reproduces the surface tensions for the compressed and expanded states, as well as compressibility at different surface tensions. In-depth analysis of simulations of a LSM in the presence of five different drug concentrations shows that mometasone alters the structure and dynamics of the LSM in a concentration-dependent manner. Mometasone induces a collapse in the monolayer that is affected by the surfactant protein and surface tension. Overall, these findings suggest that the surfactant proteins, surface tension and drug concentration are all critical components affecting monolayer stability and drug adsorption. The outcomes of this study may be beneficial for a more in-depth understanding of how mometasone is adsorbed by lung surfactants.

Shape matters-the interaction of gold nanoparticles with model lung surfactant monolayers

The lung surfactant monolayer (LSM) forms the main biological barrier for any inhaled particles to enter our bloodstream, including gold nanoparticles (AuNPs) present as air pollutants and under investigation for use in biomedical applications. Understanding the interaction of AuNPs with lung surfactant can assist in understanding how AuNPs enter our lungs. In this study, we use coarse-grained molecular dynamics simulations to investigate the effect of four different shape D AuNPs (spherical, box, icosahedron and rod) on the structure and dynamics of a model LSM, with a particular focus on differences resulting from the shape of the AuNP. Monolayer-AuNP systems were simulated in two different states: the compressed state and the expanded state, representing inhalation and exhalation conditions, respectively. Our results indicate that the compressed state is more affected by the presence of the AuNPs than the expanded state. Our results show that in the compressed state, the AuNPs prevent the monolayer from reaching the close to zero surface tension required for normal exhalation. In the compressed state, all four nanoparticles (NPs) reduce the lipid order parameters and cause a thinning of the monolayer where the particles drag surfactant molecules into the water phase. Comparing the different properties shows no trend concerning which shape has the biggest effect on the monolayer, as shape-dependent effects vary among the different properties. Insights from this study might assist future work of how AuNP shapes affect the LSM during inhalation or exhalation conditions.

Simulated Breathing: Application of Molecular Dynamics Simulations to Pulmonary Lung Surfactant

In this review, we delve into the topic of the pulmonary surfactant (PS) system, which is present in the respiratory system. The total composition of the PS has been presented and explored, from the types of cells involved in its synthesis and secretion, down to the specific building blocks used, such as the various lipid and protein components. The lipid and protein composition varies across species and between individuals, but ultimately produces a PS monolayer with the same role. As such, the composition has been investigated for the ways in which it imposes function and confers peculiar biophysical characteristics to the system as a whole. Moreover, a couple of theories/models that are associated with the functions of PS have been addressed. Finally, molecular dynamic (MD) simulations of pulmonary surfactant have been emphasized to not only showcase various group’s findings, but also to demonstrate the validity and importance that MD simulations can have in future research exploring the PS monolayer system.

Modeling Gas-Liquid Interfaces by Dissipative Particle Dynamics: Adsorption and Surface Tension of Cetyl Trimethyl Ammonium Bromide at the Air-Water Interface

Adsorption of surfactants at gas-liquid interfaces that causes reduction in the surface tension is a classical problem in colloid and interface science with multiple practical applications in oil and gas recovery, separations, cosmetics, personal care, and biomedicine. Here, we develop an original coarse-grained model of the liquid-gas interface within the conventional dissipative particle dynamics (DPD) framework with the goal of quantitatively predicting the surface tension in the presence of surfactants. As a practical case-study example, we explore the adsorption of the cationic surfactant cetyl trimethyl ammonium bromide (CTAB) on the air-water interface. The gas phase is modeled as a DPD fluid composed of fictitious hard-core "gas" beads with exponentially decaying repulsive potentials to prevent penetration of the liquid phase components. A rigorous parametrization scheme is proposed based on matching the bulk and interfacial properties of water and octane taken as the reference compounds. Quantitative agreement between the simulated and experimental surface tension of CTAB solutions is found for a wide range of bulk surfactant concentrations (∼10^-3 to ∼1 mmol/L) with the reduction of the surface tension from ∼72 mN/m (pure water) to the limiting value of ∼37.5 mN/m at the critical micelle concentration. The gas phase DPD model with the proposed parametrization scheme can be extended and applied to modeling various gas-liquid interfaces with surfactant and lipid monolayers, such as bubble suspensions, foams, froths, etc.

A Mechanical Mechanism for Vitamin E Acetate in E-cigarette/Vaping-Associated Lung Injury

The outbreak of electronic-cigarette/vaping-associated lung injury (EVALI) has made thousands ill. This lung injury has been attributed to a physical interaction between toxicants from the vaping solution and the pulmonary surfactant. In particular, studies have implicated vitamin E acetate as a potential instigator of EVALI. Pulmonary surfactant is vital to proper respiration through the mechanical processes of adsorption and interface stability to achieve and maintain low surface tension at the air-liquid interface. Using neutron spin echo spectroscopy, we investigate the impact of vitamin E acetate on the mechanical properties of two lipid-only pulmonary surfactant mimics: pure 1,2-dipalmitoyl-sn-glycero-3-phosphocholine and a more comprehensive lipid mixture. It was found that increasing vitamin E acetate concentration nonlinearly increased membrane fluidity and area compressibility to a plateau. Softer membranes would promote adsorption to the air-liquid interface during inspiration as well as collapse from the interface during expiration. These findings indicate the potential for the failure of the pulmonary surfactant upon expiration, attributed to monolayer collapse. This collapse could contribute to the observed EVALI signs and symptoms, including shortness of breath and pneumonitis.

The role of disulfide bonds in a Solanum tuberosum saposin-like protein investigated using molecular dynamics

The Solanum tuberosum plant specific insert (StPSI) has a defensive role in potato plants, with the requirements of acidic pH and anionic lipids. The StPSI contains a set of three highly conserved disulfide bonds that bridge the protein’s helical domains. Removal of these bonds leads to enhanced membrane interactions. This work examined the effects of their sequential removal, both individually and in combination, using all-atom molecular dynamics to elucidate the role of disulfide linkages in maintaining overall protein tertiary structure. The tertiary structure was found to remain stable at both acidic (active) and neutral (inactive) pH despite the removal of disulfide linkages. The findings include how the dimer structure is stabilized and the impact on secondary structure on a residue-basis as a function of disulfide bond removal. The StPSI possesses an extensive network of inter-monomer hydrophobic interactions and intra-monomer hydrogen bonds, which is likely the key to the stability of the StPSI by stabilizing local secondary structure and the tertiary saposin-fold, leading to a robust association between monomers, regardless of the disulfide bond state. Removal of disulfide bonds did not significantly impact secondary structure, nor lead to quaternary structural changes. Instead, disulfide bond removal induces regions of amino acids with relatively higher or lower variation in secondary structure, relative to when all the disulfide bonds are intact. Although disulfide bonds are not required to preserve overall secondary structure, they may have an important role in maintaining a less plastic structure within plant cells in order to regulate membrane affinity or targeting.

Multi-walled carbon nanotubes (MWCNTs) transformed THP-1 macrophages into foam cells: Impact of pulmonary surfactant component dipalmitoylphosphatidylcholine

Pulmonary surfactant or its components can function as barriers toward nanomaterials (NMs) entering pulmonary systems. However, since pulmonary surfactant mainly consists of lipids, it may be necessary to investigate the effects of co-exposure to NMs and pulmonary surfactant or its components on lipid metabolism and related signaling pathways. Recently we found that multi-walled carbon nanotubes (MWCNTs) transformed THP-1 macrophages into lipid-laden foam cells via ER stress pathway. Here this study further investigated the impact of pulmonary surfactant component dipalmitoylphosphatidylcholine (DPPC) on this process. Up to 64 μg/mL hydroxylated or carboxylated MWCNTs induced lipid accumulation and IL-6 release in THP-1 macrophages, accompanying with increased oxidative stress and p-chop proteins (biomarker for ER stress). Incubation with 100 μg/mL DPPC led to MWCNT surface coating but did not significantly alter MWCNT internalization, lipid burden or IL-6 release. However, lipidomics indicated that DPPC altered lipid profliles in MWCNT-exposed cells. DPPC also led to a higher level of de novo lipogenesis regulator FASN in cells exposed to hydroxylated MWCNTs, as well as a higher level of p-chop and scavenger receptor MSR1 in cells exposed to carboxylated MWCNTs. Combined, DPPC did not significantly affect MWCNT-induced lipid accumulation but altered lipid components and ER stress in macrophages.

Physical properties of phospholipids and integral proteins and their biofunctional roles in pulmonary surfactant from molecular dynamics simulation

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.

This work deals with a quantitative investigation of the physical properties of pulmonary surfactant near melting temperature.

Nanoparticle translocation across the lung surfactant film regulated by grafting polymers

Nanoparticle-based pulmonary drug delivery has gained significant attention due to its ease of administration, increased bioavailability, and reduced side effects caused by a high systemic dosage. After being delivered into the deep lung, the inhaled nanoparticles first interact with the lung surfactant lining layer composed of phospholipids and surfactant proteins and then potentially cause the dysfunction of the lung surfactant. Conditioning the surface properties of nanoparticles with grafting polymers to avoid these side effects is of crucial importance to the efficiency and safety of pulmonary drug delivery. Herein, we perform coarse-grained molecular simulations to decipher the involved mechanism responsible for the translocation of the polymer-grafted Au nanoparticles across the lung surfactant film. The simulations illustrate that conditioning of the grafting polymers, including their length, terminal charge, and grafting density, can result in different translocation processes. Based on the energy analysis, we find that these discrepancies in translocation stem from the affinity of the nanoparticles with the lipid tails and heads and their contact with the proteins, which can be tuned by the surface polarity and surface charge of the nanoparticles. We further demonstrate that the interaction between the nanoparticles and the lung surfactant is related to the depletion of the lipids and proteins during translocation, which affects the surface tension of the surfactant film. The change in the surface tension in turn affects the nanoparticle translocation and the collapse of the surfactant film. These results can help understand the adverse effects of the nanoparticles on the lung surfactant film and provide guidance to the design of inhaled nanomedicines for improved permeability and targeting.

The role of SP-B1–25 peptides in lung surfactant monolayers exposed to gold nanoparticles

Lung surfactant monolayer’s (acts as the first line barrier for inhaled nanoparticles) components (lipids and peptides) rearrange themselves by the influence of exposed gold nanoparticles at various stages of the breathing cycle.

The Fats of Life: Using Computational Chemistry to Characterise the Eukaryotic Cell Membrane

Our current knowledge of the structural dynamics and complexity of lipid bilayers is still developing. Computational techniques, especially molecular dynamics simulations, have increased our understanding significantly as they allow us to model functions that cannot currently be experimentally resolved. Here we review available computational tools and techniques, the role of the major lipid species, insights gained into lipid bilayer structure and function from molecular dynamics simulations, and recent progress towards the computational modelling of the physiological complexity of eukaryotic lipid bilayers.

Molecular Ordering in Lipid Monolayers: An Atomistic Simulation

We report on atomistic simulations of DPPC lipid monolayers using the CHARMM36 lipid force field (and also the Slipid force field as a control case), combined with a four-point OPC water model. The entire two-phase region where domains of the "liquid-condensed" (LC) phase coexist with domains of the "liquid-expanded" (LE) phase has been explored. The simulations are long enough that the complete phase-transition stage, with two domains coexisting in the monolayer, is reached in all cases. Also, system sizes used are larger than those in previous works. As expected, domains of the minority phase are elongated, emphasizing the importance of anisotropic van der Waals and/or electrostatic dipolar interactions in the monolayer plane. The molecular structure is quantified in terms of distribution functions for the hydrocarbon chains and the PN dipoles. In contrast to previous work, where average distributions are calculated, distributions are here extracted for each of the coexisting phases by first identifying lipid molecules that belong to either LC or LE regions. In the case of the CHARMM36 force field, the three-dimensional distributions show that the average tilt angle of the chains with respect to the normal outward direction is (39.0 ± 0.1)° in the LC phase and (48.1 ± 0.5)° in the LC phase. In the case of the PN dipoles, the distributions indicate a tilt angle of (110.8 ± 0.5)° in the LC phase and (112.5 ± 0.5)° in the LE phase. These results are quantitatively different from those in previous works, which indicated a smaller normal component of the PN dipole. Also, the distributions of the monolayer-projected chains and PN dipoles have been calculated. Chain distributions peak along a particular direction in the LC domains, while they are uniform in the LE phase. Long-range ordering associated with the projected PN dipoles is absent in both phases. These results strongly suggest that LC domains do not exhibit dipolar ordering in the plane of the monolayer, the effect of these components being averaged out at short distances. Therefore, the only relevant component of the molecular dipoles, with regard to both intra- and long-range interdomain interactions, is normal to the monolayer. Also, the local orientation of chain projections is almost constant in LC domains and points in the direction along which domains are elongated, suggesting that the line tension driving the phase transition might be anisotropic with respect to the interfacial domain boundary.

All-Atom Molecular Dynamics Simulations of Dimeric Lung Surfactant Protein B in Lipid Multilayers

Although lung surfactant protein B (SP-B) is an essential protein that plays a crucial role in breathing, the details of its structure and mechanism are not well understood. SP-B forms covalent homodimers, and in this work we use all-atom molecular dynamics simulations to study dimeric SP-B’s structure and its behavior in promoting lipid structural transitions. Four initial system configurations were constructed based on current knowledge of SP-B’s structure and mechanism, and the protein maintained a helicity consistent with experiment in all systems. Several SP-B-induced lipid reorganization behaviors were observed, and regions of the protein particularly important for these activities included SP-B’s “central loop” and “hinge” regions. SP-B dimers with one subunit initially positioned in each of two adjacent bilayers appeared to promote close contact between two bilayers. When both subunits were initially positioned in the same bilayer, SP-B induced the formation of a defect in the bilayer, with water penetrating into the centre of the bilayer. Similarly, dimeric SP-B showed a propensity to interact with preformed interpores in the bilayer. SP-B dimers also promoted bilayer thinning and creasing. This work fleshes out the atomistic details of the dimeric SP-B structures and SP-B/lipid interactions that underlie SP-B’s essential functions.

Poly(amidoamine) Dendrimer as a Respiratory Nanocarrier: Insights from Experiments and Molecular Dynamics Simulations

Pulmonary drug delivery is superior to the systemic administration in treating lung diseases. An optimal respiratory nanocarrier should be able to efficiently and safely cross the pulmonary surfactant film, which serves as the first biological barrier for respiratory delivery and plays paramount roles in maintaining the proper mechanics of breathing. In this work, we focused on the interactions between poly(amidoamine) (PAMAM) dendrimers and a model pulmonary surfactant. With combined Langmuir monolayer experiments and coarse-grained molecular dynamics simulations, we studied the effect of environmental temperature, size, and surface property of PAMAM dendrimers (G3-OH, G3-NH₂, G5-OH, and G5-NH₂) on the dipalmitoylphosphatidylcholine (DPPC) monolayer. Our simulations indicated that the environmental temperature could significantly affect the influence of PAMAM dendrimers on the DPPC monolayer. Therefore, results obtained at room temperature cannot be directly applied to elucidate interactions at body temperature. Simulations at body temperature found that all tested PAMAM dendrimers can easily penetrate the lipid monolayer during the monolayer expansion process (mimicking "inhalation"), and the cationic PAMAM dendrimers (-NH₂) show promising penetration ability during the monolayer compression process (mimicking "expiration"). Larger PAMAM dendrimers (G5) adsorbed onto the lipid monolayer tend to induce structural collapse and inhibit normal phase transitions of the lipid monolayer. These adverse effects could be mitigated in the subsequent expansion-compression cycle. These findings suggest that the PAMAM dendrimer may be used as a potential respiratory drug nanocarrier.

Interaction of nano carbon particles and anthracene with pulmonary surfactant: The potential hazards of inhaled nanoparticles

Understanding the alteration of the air-liquid interfacial properties of pulmonary surfactant (PS) in the presence of nanoparticles (NPs) and polycyclic aromatic hydrocarbons (PAHs) is particularly important for pulmonary risk assessment. Here, we investigated the interaction of natural PS (extracted from pig's lungs) with nano carbon particles (NCPs) and anthracene as a representative PAH. Our results showed that PS exhibited a significant solubilization effect on anthracene. Solubilization experiment for the substructures of PS demonstrated that the mixed phospholipid components of PS played the primary role in the solubilization of PS for anthracene. Adsorption experiment indicated that in the mixed system of PS, NCPs, and anthracene, PS can inhibit the adsorption of anthracene on NCPs due to the solubilization, agglomeration, and competitive adsorption. In addition, the surface tension, phase behavior, and foaming ability of PS were obviously altered in the presence of NCPs. These findings indicate that the solubilization effect of PS on anthracene, the inhibitive effect of PS for the adsorption of anthracene on NCPs, and the alternation of air-liquid interfacial properties of PS containing NCPs may increase the pulmonary risk in the exposure of atmospheric environment containing both PAHs and NCPs.

Role of Lipid Coating in the Transport of Nanodroplets across the Pulmonary Surfactant Layer Revealed by Molecular Dynamics Simulations

Hydrophilic drugs can be delivered into lungs via nebulization for both local and systemic therapies. Once inhaled, ultrafine nanodroplets preferentially deposit in the alveolar region, where they first interact with the pulmonary surfactant (PS) layer, with nature of the interaction determining both efficiency of the pulmonary drug delivery and extent of the PS perturbation. Here, we demonstrate by molecular dynamics simulations the transport of nanodroplets across the PS layer being improved by lipid coating. In the absence of lipids, bare nanodroplets deposit at the PS layer to release drugs that can be directly translocated across the PS layer. The translocation is quicker under higher surface tensions but at the cost of opening pores that disrupt the ultrastructure of the PS layer. When the PS layer is compressed to lower surface tensions, the nanodroplet prompts collapse of the PS layer to induce severe PS perturbation. By coating the nanodroplet with lipids, the disturbance of the nanodroplet on the PS layer can be reduced. Moreover, the lipid-coated nanodroplet can be readily wrapped by the PS layer to form vesicular structures, which are expected to fuse with the cell membrane to release drugs into secondary organs. Properties of drug bioavailability, controlled drug release, and enzymatic tolerance in real systems could be improved by lipid coating on nanodroplets. Our results provide useful guidelines for the molecular design of nanodroplets as carriers for the pulmonary drug delivery.

Excessive Extracellular ATP Desensitizes P2Y2 and P2X4 ATP Receptors Provoking Surfactant Impairment Ending in Ventilation-Induced Lung Injury

Stretching the alveolar epithelial type I (AT I) cells controls the intercellular signaling for the exocytosis of surfactant by the AT II cells through the extracellular release of adenosine triphosphate (ATP) (purinergic signaling). Extracellular ATP is cleared by extracellular ATPases, maintaining its homeostasis and enabling the lung to adapt the exocytosis of surfactant to the demand. Vigorous deformation of the AT I cells by high mechanical power ventilation causes a massive release of extracellular ATP beyond the clearance capacity of the extracellular ATPases. When extracellular ATP reaches levels >100 μM, the ATP receptors of the AT II cells become desensitized and surfactant impairment is initiated. The resulting alteration in viscoelastic properties and in alveolar opening and collapse time-constants leads to alveolar collapse and the redistribution of inspired air from the alveoli to the alveolar ducts, which become pathologically dilated. The collapsed alveoli connected to these dilated alveolar ducts are subject to a massive strain, exacerbating the ATP release. After reaching concentrations >300 μM extracellular ATP acts as a danger-associated molecular pattern, causing capillary leakage, alveolar space edema, and further deactivation of surfactant by serum proteins. Decreasing the tidal volume to 6 mL/kg or less at this stage cannot prevent further lung injury.

Domain size polydispersity effects on the structural and dynamical properties in lipid monolayers with phase coexistence

In lipid monolayers with phase coexistence, domains of the liquid-condensed phase always present size polydispersity. However, very few theoretical works consider size distribution effects on the monolayer properties. Because of the difference in surface densities, domains have excess dipolar density with respect to the surrounding liquid expanded phase, originating a dipolar inter-domain interaction. This interaction depends on the domain area, and hence the presence of a domain size distribution is associated with interaction polydispersity. Inter-domain interactions are fundamental to understanding the structure and dynamics of the monolayer. For this reason, it is expected that polydispersity significantly alters monolayer properties. By means of Brownian dynamics simulations, we study the radial distribution function (RDF), the average mean square displacement and the average time-dependent self-diffusion coefficient, D(t), of lipid monolayers with normally distributed size domains. For this purpose, we vary the relevant system parameters, polydispersity and interaction strength, within a range of experimental interest. We also analyze the consequences of using a monodisperse model to determine the interaction strength from an experimental RDF. We find that polydispersity strongly affects the value of the interaction strength, which is greatly underestimated if polydispersity is not considered. However, within a certain range of parameters, the RDF obtained from a polydisperse model can be well approximated by that of a monodisperse model, by suitably fitting the interaction strength, even for 40% polydispersities. For small interaction strengths or small polydispersities, the polydisperse systems obtained from fitting the experimental RDF have an average mean square displacement and D(t) in good agreement with that of the monodisperse system.

Atomistic Model for Nearly Quantitative Simulations of Langmuir Monolayers

Lung surfactant and a tear film lipid layer are examples of biologically relevant macromolecular structures found at the air-water interface. Because of their complexity, they are often studied in terms of simplified lipid layers, the simplest example being a Langmuir monolayer. Given the profound biological significance of these lipid assemblies, there is a need to understand their structure and dynamics on the nanoscale, yet there are not many techniques able to provide this information. Atomistic molecular dynamics simulations would be a tool fit for this purpose; however, the simulation models suggested until now have been qualitative instead of quantitative. This limitation has mainly stemmed from the challenge to correctly describe the surface tension of water with simulation parameters compatible with other biomolecules. In this work, we show that this limitation can be overcome by using the recently introduced four-point OPC water model, whose surface tension for water is demonstrated to be quantitatively consistent with experimental data and which is also shown to be compatible with the commonly employed lipid models. We further establish that the approach of combining the OPC four-point water model with the CHARMM36 lipid force field provides nearly quantitative agreement with experiments for the surface pressure-area isotherm for POPC and DPPC monolayers, also including the experimentally observed phase coexistence in a DPPC monolayer. The simulation models reported in this work pave the way for nearly quantitative atomistic studies of lipid-rich biological structures at air-water interfaces.

Toxicity of ZnO nanoparticles (NPs) to A549 cells and A549 epithelium in vitro: Interactions with dipalmitoyl phosphatidylcholine (DPPC)

Once inhaled, nanoparticles (NPs) will first interact with lung surfactant system, which may influence the colloidal aspects of NPs and consequently the toxic potential of NPs to pulmonary cells. In this study, we investigated the effects of dipalmitoyl phosphatidylcholine (DPPC), the major component in lung surfactant, on stability and toxicity of ZnO NPs. The presence of DPPC increased the UV-vis spectra, hydrodynamic size, Zeta potential and dissolution rate of ZnO NPs, which indicates that DPPC might interact with NPs and affect the colloidal stability of NPs. Exposure to ZnO NPs induced cytotoxicity associated with increased intracellular Zn ions but not superoxide in A549 cells. In A549 epithelium model, exposure to ZnO NPs induced cytotoxicity and decreased the release of interleukin 6 (IL-6) without a significant effect on epithelial permeability rate. Co-exposure of A549 cells or A549 epithelium model to DPPC and ZnO NPs induced a higher release of lactate dehydrogenase (LDH) and interleukin-6 (IL-6) compared with the exposure of ZnO NPs alone. We concluded that the presence of DPPC could influence the colloidal stability of ZnO NPs and increase the damage of NPs to membrane probably due to the increased positive surface charge.

Unperturbed hydrocarbon chains and liquid phase bilayer lipid chains: a computer simulation study

In this work, the properties of saturated and unsaturated fatty acid acyl chains 16:0, 18:0, 18:1(n-9)cis, 18:2(n-6)cis, 18:3(n-3)cis, 18:4(n-3)cis, 18:5(n-3)cis, 20:4(n-6)cis, 20:5(n-3)cis and 22:6(n-3)cis in a bilayer liquid crystalline state and similar hydrocarbon chains (with CH[Formula: see text] terminal groups instead of C=O groups) in the unperturbed state characterised by a lack of long-range interaction were investigated. The unperturbed hydrocarbon chains were modelled by Monte Carlo simulations at temperature [Formula: see text] K; sixteen fully hydrated homogeneous liquid crystalline phosphatidylcholine bilayers containing these chains were studied by molecular dynamics simulations at the same temperature. To eliminate effects of the simulation parameters, the molecular dynamics and Monte Carlo simulations were carried out using the same structural data and force field coefficients. From these computer simulations, the average distances between terminal carbon atoms of the chains (end-to-end distances) were calculated and compared. The trends in the end-to-end distances obtained for the unperturbed chains were found to be qualitatively similar to those obtained for the same lipid chains in the bilayers. So, for understanding of a number of processes in biological membranes (e.g., changes in fatty acid composition caused by environmental changes such as temperature and pressure), it is possible to use, at least as a first approximation, the relationships between the structure and properties for unperturbed or isolated hydrocarbon chains.

Transport of nanoparticles across pulmonary surfactant monolayer: a molecular dynamics study

Three types of nanoparticles, including hydrophobic nanoparticles, hydrophilic nanoparticles, and hydrophilic nanoparticles coated with lipids, were found by our molecular dynamics simulations to be transported across the pulmonary surfactant monolayer, but via different pathways, which affect their subsequent interactions with target cell membranes.

Surface patterning of single-walled carbon nanotubes enhances their perturbation on a pulmonary surfactant monolayer: frustrated translocation and bilayer vesiculation

In contrast to SWCNTs with unique surface properties, the surface patterning of SWCNTs is found to enhance their perturbation on the pulmonary surfactantsviafrustrated translocation and destructive bilayer vesiculation.