The reconciliation between the experimental and calculated octanol-water partition coefficient of 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine using atomistic molecular dynamics: an open question

The octanol-water partition coefficient of 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) was investigated using atomistic molecular dynamics simulations via thermodynamic integration and multistate Bennett acceptance ratio methods. The GAFF and CHARMM36 force fields were used with six water models widely used in molecular dynamics simulations. The OPC4 water model provided the best agreement with the experimental octanol-water partition coefficient of DPPC using the two force fields. However, there is still plenty of room for improvement in water models with correct estimation of surface tension that uses better and suitable non-bonded interaction parameters between water-water and water-DPPC. The Gibbs free energy of transferring DPPC from octanol to water phase was calculated to be 19.8 ± 0.3 and 20.2 ± 0.3 kcal mol-1, giving a partition coefficient of 14.5 ± 0.4 and 14.8 ± 0.3 for the GAFF and CHARMM36 force fields, respectively. This study reinforces the importance of developing new water models that reproduce experimental surface tensions to reconcile the water-water and water-DPPC non-bonded interactions and the existing discrepancy between experimental measurements of amphiphilic molecules that are important in many areas of scientific applications and industry such as biophysics, surfactant, colloids, membranes, medicine, nanotechnology, and food and pharmaceutical industries, and so on. It raises two important open questions: Is the experimental octanol-water partition coefficient of DPPC reliable? Or is its calculation accurate using the OPC4 water model? With respect to the experimental measurements, there may be non-treated aspects such as the formation of aggregates in aqueous phase and limit of detection of the applied method. And, in the calculation, some effects are not possible to be considered in a correct way or viable time such as calculating quantum effects, sampling all conformations, considering phase transitions, and correctly evaluating the intermolecular forces to estimate an accurate surface tension.Communicated by Ramaswamy H. Sarma.

Molecular Dynamics Simulations of Self-Assembling Colloids in Fed-State Human Intestinal Fluids and Their Solubilization of Lipophilic Drugs

Bioavailability of oral drugs often depends on how soluble the active pharmaceutical ingredient is in the fluid present in the small intestine. For efficient drug discovery and development, computational tools are needed for estimating this drug solubility. In this paper, we examined human intestinal fluids collected in the fed state, with coarse-grained molecular dynamics simulations. The experimentally obtained concentrations in aspirated duodenal fluids from five healthy individuals were used in three simulation sets to evaluate the importance of the initial distribution of molecules and the presence of glycerides in the simulation box when simulating the colloidal environment of the human intestinal fluid. We observed self-assembly of colloidal structures of different types: prolate, elongated, and oblate micelles, and vesicles. Glycerides were important for the formation of vesicles, and their absence was shown to induce elongated micelles. We then simulated the impact of digestion and absorption on the different colloidal types. Finally, we looked at the solubilization of three model compounds of increasing lipophilicity (prednisolone, fenofibrate, and probucol) by calculating contact ratios of drug-colloid to drug-water. Our simulation results of colloidal interactions with APIs were in line with experimental solubilization data but showed a dissimilarity to solubility values when comparing fasted-/fed-state ratios between two of the APIs. This work shows that coarse-grained molecular dynamics simulation is a promising tool for investigation of the intestinal fluids, in terms of colloidal attributes and drug solubility.

The concentration-dependent effect of hydrocortisone on the structure of model lung surfactant monolayer by using an in silico approach

Understanding the adsorption mechanism of corticosteroids in the lung surfactant requires the knowledge of corticosteroid molecular interactions with lung surfactant monolayer (LSM). We employed coarse-grained molecular dynamics simulation to explore the action of hydrocortisone on an LSM comprised of a phospholipid, cholesterol and surfactant protein. The structural and dynamical morphology of the lung surfactant monolayer at different surface tensions were investigated to assess the monolayer compressibility. The simulations were also conducted at the two extreme ends of breathing cycles: exhalation (0 mN m^-1 surface tension) and inhalation (20 mN m^-1 surface tension). The impact of surface tension and hydrocortisone concentration on the monolayer compressibility and stability are significant, resulting the monolayer expansion at higher surface tension. However, at low surface tension, the highly compressed monolayer induces monolayer instability in the presence of the drug due to the accumulation of surfactant protein and drug. The constant area per lipid simulation results demonstrate that the surface pressure-area isotherms show a decrease in area-per-lipid with increased drug concentration. The drug-induced expansion causes considerable instability in the monolayer after a specific drug concentration is attained at inhalation breathing condition, whereas, for exhalation breathing, the monolayer gets more compressed, causing the LSM to collapse. The monolayer collapse occurs for inhalation due to the higher drug concentration, whereas for exhalation due to the accumulation of surfactant proteins and drugs. The findings from this study will aid in enhancing the knowledge of molecular interactions of corticosteroid drugs with lung surfactants to treat respiratory diseases.

Lung surfactant negatively affects the photodynamic inactivation of bacteria-in vitro and molecular dynamic simulation analyses

In the context of the rapid increase of antibiotic-resistant infections, in particular of pneumonia, antimicrobial photodynamic therapy (aPDT), the microbiological application of photodynamic therapy (PDT), comes in as a promising treatment alternative since the induced damage and resultant death are not dependent on a specific biomolecule or cellular pathway. The applicability of aPDT using the photosensitizer indocyanine green with infrared light has been successfully demonstrated for different bacterial agents in vitro, and the combination of pulmonary delivery using nebulization and external light activation has been shown to be feasible. However, there has been little progress in obtaining sufficient in vivo efficacy results. This study reports the lung surfactant as a significant suppressor of aPDT in the lungs. In vitro, the clinical surfactant Survanta® reduced the aPDT effect of indocyanine green, Photodithazine®, bacteriochlorin-trizma, and protoporphyrin IX against Streptococcus pneumoniae. The absorbance and fluorescence spectra, as well as the photobleaching profile, suggested that the decrease in efficacy is not a result of singlet oxygen quenching, while a molecular dynamics simulation showed an affinity for the polar head groups of the surfactant phospholipids that likely impacts uptake of the photosensitizers by the bacteria. Methylene blue is the exception, likely because its high water solubility confers a higher mobility when interacting with the surfactant layer. We propose that the interaction between lung surfactant and photosensitizer must be taken into account when developing pulmonary aPDT protocols.

Concentration-Dependent Effect of the Steroid Drug Prednisolone on a Lung Surfactant Monolayer

The lung surfactant monolayer (LSM) is the main barrier for particles entering the lung, including steroid drugs used to treat lung diseases. The present study combines Langmuir experiments and coarse-grained (CG) molecular dynamics simulations to investigate the concentration-dependent effect of steroid drug prednisolone on the structure and morphology of a model LSM. The surface pressure-area isotherms for the Langmuir monolayers reveal a concentration-dependent decrease in area per lipid (APL). Results from simulations at a fixed surface tension, representing inhalation and exhalation conditions, suggest that at high drug concentrations, prednisolone induces a collapse of the LSM, which is likely caused by the inability of the drug to diffuse into the bilayer. Overall, the monolayer is most susceptible to drug-induced collapse at surface tensions representing exhalation conditions. The presence of cholesterol also exacerbates the instability. The findings of this investigation might be helpful for better understanding the interaction between steroid drug prednisolone and lung surfactants in relation to off-target effects.

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.

Molecular basis of transport of surface functionalised gold nanoparticles to pulmonary surfactant

Ligands like alkanethiol (e.g. dodecanethiol, hexadecanethiol, etc.) and polymers (e.g. poly(vinyl pyrrolidone), polyethylene glycol-thiol) capped to the gold nanoparticles (AuNPs) are widely used in biomedical field as drug carriers and as promising materials for probing and manipulating cellular processes. Ligand functionalised AuNPs are known to interact with the pulmonary surfactant (PS) monolayer once reaching the alveolar region. Therefore, it is crucial to understand the interaction between AuNPs and PS monolayers. Using coarse-grained molecular dynamics simulations, the effect of ligand density, and ligand length have been studied for two classes of ligands on a PS model monolayer consisting of DPPC, POPG, cholesterol and SP-B (mini-peptide). The ligands considered in this study are alkanethiol and polyethylene glycol (PEG) thiol as examples of hydrophobic and hydrophilic ligands, respectively. It was observed that the interaction between AuNPs and PS changes the biophysical properties of PS monolayer in compressed and expanded states. The AuNPs with hydrophilic ligand, can penetrate through the monolayer more easily, while the AuNPs with hydrophobic ligand are embedded in the monolayer and participated in deforming the monolayer structure particularly the monolayer in the compressed state. The bare AuNPs hinder to lower the monolayer surface tension value at the interface, however introducing ligand to the bare AuNPs or increasing the ligand length and density have an impact of lowering of monolayer surface tension to a minor extent. The simulation results guide the design of ligand protected NPs as drug carriers and can identify the nanoparticles' potential side effects on lung surfactant.

Molecular-level observations of the behavior of ligand functionalised gold nanoparticles with a lipid monolayers.

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.

Pulmonary surfactant: a unique biomaterial with life-saving therapeutic applications

Pulmonary surfactant is a complex lipoprotein mixture secreted into the alveolar lumen by type 2 pneumocytes, which is composed by tens of different lipids (approximately 90% of its entire mass) and surfactant proteins (approximately 10% of the mass). It is crucially involved in maintaining lung homeostasis by reducing the values of alveolar liquid surface tension close to zero at end-expiration, thereby avoiding the alveolar collapse, and assembling a chemical and physical barrier against inhaled pathogens. A deficient amount of surfactant or its functional inactivation is directly linked to a wide range of lung pathologies, including the neonatal respiratory distress syndrome. This paper reviews the main biophysical concepts of surfactant activity and its inactivation mechanisms, and describes the past, present and future roles of surfactant replacement therapy, focusing on the exogenous surfactant preparations marketed worldwide and new formulations under development. The closing section describes the pulmonary surfactant in the context of drug delivery. Thanks to its peculiar composition, biocompatibility, and alveolar spreading capability, the surfactant may work not only as a shuttle to the branched anatomy of the lung for other drugs but also as a modulator for their release, opening to innovative therapeutic avenues for the treatment of several respiratory diseases.

Computational Studies of Lipid-Wrapped Gold Nanoparticle Transport Through Model Lung Surfactant Monolayers

Colloidal nanoparticles, such as gold nanoparticles (AuNPs), are promising materials for the delivery of hydrophilic drugs via the pulmonary route. The inhaled nanoparticle drug carriers primarily deposit in lung alveoli and interact with the alveolar surface known as lung surfactants. Therefore, it is vital to understand the interactions of nanocarriers with the surfactant layer. To understand the interactions at the molecular level, here we simulated model lung surfactant monolayers with phospholipid (PL)-wrapped AuNPs at the vacuum-water interface using coarse-grained molecular dynamics simulations. The PL-wrapped AuNPs quickly adsorbed into the surfactant layer, altered the structural properties of the monolayer, and at high concentrations initiated the compressed monolayer to collapse/buckle. Among the surfactant monolayer lipid components, cholesterol adsorbed to the AuNPs preferentially over PL species. The position of the adsorbed PL-AuNPs within the monolayer, and subsequent monolayer perturbation, vary depending on the monolayer phase, monolayer composition, and species of PL used as a ligand. Information provided by these molecular dynamic simulations helps to rationalize why some colloidal nanoparticles work better as nanocarriers than others and aid the design of new ones, to avoid biological toxicity and improve efficacy for pulmonary drug delivery.

Recent Open Issues in Coarse Grained Force Fields

This viewpoint intends to show recent open issues of using coarse grained models in molecular dynamics simulation. It reviews the current knowledge of the comparison between experimental and simulation data of structural and physical chemical properties that depend on the hydrophilic and hydrophobic behavior of the molecule.

Molecular Dynamics Simulations on Interindividual Variability of Intestinal Fluids: Impact on Drug Solubilization

Efficient delivery of oral drugs is dependent on their solubility in human intestinal fluid, a complex and dynamic fluid that contains colloidal structures composed of small molecules. These structures solubilize poorly water-soluble compounds, increasing their apparent solubility, and possibly their bioavailability. In this study, we conducted coarse-grained molecular dynamics simulations with data from duodenal fluid samples previously acquired from five healthy volunteers. In these simulations, we observed the self-assembly of mixed micelles of bile salts, phospholipids, and free fatty acids. The micelles were ellipsoids with a size range of 4–7 nm. Next, we investigated micelle affinities of three model drugs. The affinities in our simulation showed the same trend as literature values for the solubility enhancement of drugs in human intestinal fluids. This type of simulations is useful for studies of events and interactions taking place in the small intestinal fluid.

Penetration of antimicrobial peptides in a lung surfactant model

Molecular dynamics simulations were successfully performed to understand the absorption mechanism of antimicrobial peptides LL-37, CATH-2, and SMAP-29 in a lung surfactant model. The antimicrobial peptides quickly penetrate in the lung surfactant model in dozens or hundreds nanoseconds, but they electrostatically interact with the lipid polar heads during the simulation time of 2 μs. This electrostatic interaction should be the explanation for the inactivation of the antimicrobial peptides when co-administrated with lung surfactant. As they strongly interact with the lipid polar heads of the lung surfactant, there is no positive charge available on the antimicrobial peptide to attack the negatively charged bacteria membrane. In order to avoid the interaction of peptides with the lipid polar heads, sodium cholate was used to form nanoparticles which act as an absorption enhancer of all antimicrobial peptides used in this investigation. The nanoparticles of 150 molecules of sodium cholate with one peptide were inserted on the top of the lung surfactant model. The nanoparticles penetrated into the lung surfactant model, spreading the sodium cholate molecules around the lipid polar heads. The sodium cholate molecules seem to protect the peptides from the interaction with the lipid polar heads, leaving them free to be delivered to the water phase. The penetration of peptides alone or even the peptide nanoparticles with sodium cholate do not collapse the lung surfactant model, indicating to be a promisor drug delivery system to the lung. The implications of this finding are that antimicrobial peptides may only be co-administered with an absorption enhancer such as sodium cholate into lung surfactant in order to avoid inactivation of their antimicrobial activity.

Interaction of testosterone-based compounds with dodecyl sulphate monolayers at the air–water interface

The encapsulation of testosterone enanthate into a sodium dodecyl sulphate monolayer.