Simulations of lipid monolayers

The accumulation of methylated porphyrins in dormant cells of Mycolicibacterium smegmatis is accompanied by a decrease in membrane fluidity and an impede of the functioning of the respiratory chain

Transition of Mycolicibacterium smegmatis (Msm) and Mycobacterium tuberculosis to dormancy in vitro is accompanied by an accumulation of free methylated forms of porphyrins (tetramethyl coproporphyrin - TMC) localized in the cell wall of dormant bacteria. A study of the fluorescence anisotropy of BODIPY based fluorescent probes on individual cell level using confocal microscope revealed significant changes in this parameter for BODIPY FL C16 from 0.05 to 0.22 for vegetative and dormant Msm cells correspondingly. Similarly, the increase of TMC concentration in vegetative Msm cells grown in the presence of 5-aminolevulinic acid (a known inducer of porphyrin synthesis) resulted in an increase of BODIPY FL C16 anisotropy. These changes in TMC concentration and membrane fluidity were accompanied by an inhibition of the activity of the respiratory chain measured by oxygen consumption and a reduction of the DCPIP redox acceptor. During the first 8 h of the reactivation of the dormant Msm cells, the porphyrin content and probe fluorescent anisotropy returned to the level for vegetative bacteria. We suggested that upon transition to dormancy, an accumulation of TMC in membranes leads to a decrease in membrane fluidity, resulting in an inhibition of the respiratory chain activity. However, direct interactions of TMC with membrane bound enzymes cannot also be excluded. This, in turn, may result in the down regulation of many metabolic energy-dependent reactions as a part of mechanisms accompanying the transition to a hypometabolic state of mycobacteria.

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.

The lung surfactant activity probed with molecular dynamics simulations

The surface of pulmonary alveolar subphase is covered with a mixture of lipids and proteins. This lung surfactant plays a crucial role in lung functioning. It shows a complex phase behavior which can be altered by the interaction with third molecules such as drugs or pollutants. For studying multicomponent biological systems, it is of interest to couple experimental approach with computational modelling yielding atomic-scale information. Simple two, three, or four-component model systems showed to be useful for getting more insight in the interaction between lipids, lipids and proteins or lipids and proteins with drugs and impurities. These systems were studied theoretically using molecular dynamic simulations and experimentally by means of the Langmuir technique. A better understanding of the structure and behavior of lung surfactants obtained from this research is relevant for developing new synthetic surfactants for efficient therapies, and may contribute to public health protection.

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.

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.

Effect of alcohol on the phase separation in model membranes

The discovery of coexisting liquid-ordered and liquid-disordered phases in multicomponent lipid bilayers has received widespread attention due to its potential relevance for biological systems. One of the many open questions is how the presence of additional components affects the nature of the coexisting phases. Of particular interest is the addition of alcohols because their anesthetic properties may arise from modulating bilayer behavior. We use coarse-grained Molecular Dynamics simulations to gain insight into the partitioning preferences of linear n-alcohols into ordered and disordered bilayers alongside their effects on local membrane structure. We find that alcohols cause only small changes to membrane composition alongside a lack of significant effects on membrane thickness and lipid tail order. Cholesterol and n-alcohol trans-bilayer motion is measured and found to be near or within the range of previous atomistic results. The cholesterol flip-flop rates increase with both n-alcohol length and concentration for octanol, dodecanol, and hexadecanol, indicating a decrease in lipid order. Umbrella sampling simulations of removing cholesterol from tertiary membranes find no significant difference with or without n-alcohols at various concentrations. Simulations of a phase-separated bilayer show that octanol preferentially partitions into the liquid-disordered phase in a ratio of approximately 3:1 over the liquid-ordered phase. Furthermore, partition coefficients of alcohol in single-phase membranes show a preference for longer alcohols (dodecanol and hexadecanol) to partition preferentially into the liquid-ordered phase, while decreasing the length of the alcohol reverses this trend. Our work tests experimental results while also investigating the ability for coarse-grained MARTINI simulations to capture minute differences in model membrane spatial arrangements on the nanoscale level.

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.

The selective interactions of cationic tetra-p-guanidinoethylcalix[4]arene with lipid membranes: theoretical and experimental model studies

Behavior of cationic tetra-p-guanidinoethylcalix[4]arene (CX1) and its building block, p-guanidinoethylphenol (mCX1) in model monolayer lipid membranes was investigated using all atom molecular dynamics simulations and surface pressure measurements. Members of two classes of lipids were taken into account: zwitterionic 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and anionic 1,2-dimyristoyl-sn-glycero-3-phospho-l-serine sodium salt (DMPS) as models of eukaryotic and bacterial cell membranes, respectively. It was demonstrated that CX1 and mCX1 accumulate near the negatively charged DMPS monolayers. The adsorption to neutral monolayers was negligible. In contrast to mCX1, CX1 penetrated into the hydrophobic part of the monolayer. The latter effect, which is possible due to a flip-flop inversion of the CX1 orientation in the lipid layer compared to the aqueous phase, may be responsible for its antibacterial activity.