Interaction of Antimicrobial Lipopeptides with Bacterial Lipid Bilayers.

The resistance of pathogens to traditional antibiotics is currently a global issue of enormous concern. As the discovery and development of new antibiotics become increasingly challenging, synthetic antimicrobial lipopeptides (AMLPs) are now receiving renewed attention as a new class of antimicrobial agents. In contrast to traditional antibiotics, AMLPs act by physically disrupting the cell membrane (rather than targeting specific proteins), thus reducing the risk of inducing bacterial resistance. In this study, we use microsecond-timescale atomistic molecular dynamics simulations to quantify the interaction of a short AMLP (C16-KKK) with model bacterial lipid bilayers. In particular, we investigate how fundamental transmembrane properties change in relation to a range of lipopeptide concentrations. A number of structural, mechanical, and dynamical features are found to be significantly altered in a non-linear fashion. At 10 mol% concentration, lipopeptides have a condensing effect on bacterial bilayers, characterized by a decrease in the area per lipid and an increase in the bilayer order. Higher AMLP concentrations of 25 and 40 mol% destabilize the membrane by disrupting the bilayer core structure, inducing membrane thinning and water leakage. Important transmembrane properties such as the lateral pressure and dipole potential profiles are also affected. Potential implications on membrane function and associated proteins are discussed.

Physical properties of model biological lipid bilayers: insights from all-atom molecular dynamics simulations.

The physical properties of lipid bilayers are sensitive to the specific type and composition of the lipids that make up the many different types of cell membranes. Studying model bilayers of representative heterogeneous compositions can provide key insights into membrane functionality. In this work, we use atomistic molecular dynamics simulations to characterize key properties in a number of bilayer membranes of varying composition. We first examine basic properties, such as lipid area, volume, and bilayer thickness, of simple, homogeneous bilayers comprising several lipid types, which are prevalent in biological membranes. Such lipids are then used in simulations of heterogeneous systems representative of bacterial, mammalian, and cancer membranes. Our analysis is especially focused on depth-dependent, transmembrane profiles; in particular, we calculate lateral pressure and dipole potential profiles, two fundamental properties which play key roles in a large number of biological functions.

Effects of lipid composition on membrane permeation.

Passive permeation through lipid membranes is an essential process in biology. In vivo membranes typically consist of mixtures of lamellar and nonlamellar lipids. Lamellar lipids are characterized by their tendency to form lamellar sheet-like structures, which are predominant in nature. Nonlamellar lipids, when isolated, instead form more geometrically complex nonlamellar phases. While mixed lamellar/nonlamellar lipid membranes tend to adopt the ubiquitous lamellar bilayer structure, the presence of nonlamellar lipids is known to have profound effects on key membrane properties, such as internal distributions of stress and elastic properties, which in turn may alter related biological processes. This work focuses on one such process, i.e., permeation, by utilising atomistic molecular dynamics simulations in order to obtain transfer free energy profiles, diffusion profiles and permeation coefficients for a series of thirteen small molecules and drugs. Each permeant is tested on two bilayer membranes of different lipid composition, i.e., purely lamellar and mixed lamellar/nonlamellar. Our results indicate that the presence of nonlamellar lipids reduces permeation for smaller molecules (molecular weight < 100) but facilitates it for the largest ones (molecular weight > 100). This work represents an advancement towards the development of more realistic in silico permeability assays, which may have a substantial future impact in the area of rational drug design.

Effects of High Pressure on Phospholipid Bilayers.

The response of lipid membranes to changes in external pressure is important for many biological processes, and it can also be exploited for technological applications. In this work, we employ all-atom molecular dynamics simulations to characterize the changes in the physical properties of phospholipid bilayers brought about by high pressure (1000 bar). In particular, we study how the response differs, in relation to different chain unsaturation levels, by comparing monounsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and biunsaturated dioleoyl-phosphatidylcholine (DOPC) bilayers. Various structural, mechanical, and dynamical features are found to be altered by the pressure increase in both bilayers. Notably, for most properties, including bilayer area and thickness, lipid order parameters, lateral pressure profile, and curvature frustration energy, we observe significantly more pronounced effects for monounsaturated POPC than biunsaturated DOPC. Possible biological implications of the results obtained are discussed, especially in relation to how different lipids can control the structure and function of membrane proteins.

Effects of Lipid Composition on Bilayer Membranes Quantified by All-Atom Molecular Dynamics.

Biological bilayer membranes typically contain varying amounts of lamellar and nonlamellar lipids. Lamellar lipids, such as dioleoylphosphatidylcholine (DOPC), are defined by their tendency to form the lamellar phase, ubiquitous in biology. Nonlamellar lipids, such as dioleoylphosphatidylethanolamine (DOPE), prefer instead to form nonlamellar phases, which are mostly nonbiological. However, nonlamellar lipids mix with lamellar lipids in biomembrane structures that remain overall lamellar. Importantly, changes in the lamellar vs nonlamellar lipid composition are believed to affect membrane function and modulate membrane proteins. In this work, we employ atomistic molecular dynamics simulations to quantify how a range of bilayer properties are altered by variations in the lamellar vs nonlamellar lipid composition. Specifically, we simulate five DOPC/DOPE bilayers at mixing ratios of 1/0, 3/1, 1/1, 1/3, and 0/1. We examine properties including lipid area and bilayer thickness, as well as the transmembrane profiles of electron density, lateral pressure, electric field, and dipole potential. While the bilayer structure is only marginally altered by lipid composition changes, dramatic effects are observed for the lateral pressure, electric field, and dipole potential profiles. Possible implications for membrane function are discussed.

Modeling Short-Term Maximum Individual Exposure from Airborne Hazardous Releases in Urban Environments. Part I: Validation of a Deterministic Model with Field Experimental Data.

The release of airborne hazardous substances in the atmosphere has a direct effect on human health as, during the inhalation, an amount of concentration is inserted through the respiratory system into the human body, which can cause serious or even irreparable damage in health. One of the key problems in such cases is the prediction of the maximum individual exposure. Current state of the art methods, which are based on the concentration cumulative distribution function and require the knowledge of the concentration variance and the intermittency factor, have limitations. Recently, authors proposed a deterministic approach relating maximum individual exposure to parameters such as the fluctuation intensity and the concentration integral time scale. The purpose of the first part of this study is to validate the deterministic approach with the extensive dataset of the MUST (Mock Urban Setting Test) field experiment. This dataset includes 81 trials, which practically cover various atmospheric conditions and stability classes and contains in total 4004 non-zero concentration sensor data with time resolutions of 0.01-0.02 s. The results strengthen the usefulness of the deterministic model in predicting short-term maximum individual exposure. Another important output is the estimation of the methodology uncertainty involved.

Direct Mixing of Atomistic Solutes and Coarse-Grained Water.

We present a new dual-resolution approach for coupling atomistic and coarse-grained models in molecular dynamics simulations of hydrated systems. In particular, a coarse-grained point dipolar water model is used to solvate molecules represented with standard all-atom force fields. A unique characteristic of our methodology is that the mixing of resolutions is direct, meaning that no additional or ad hoc scaling factors, intermediate regions, or extra sites are required. To validate the methodology, we compute the hydration free energy of 14 atomistic small molecules (analogs of amino acid side chains) solvated by the coarse-grained water. Remarkably, our predictions reproduce the experimental data as accurately as the predictions from state-of-the-art fully atomistic simulations. We also show that the hydration free energy of the coarse-grained water itself is in comparable or better agreement with the experimental value than the predictions from all but one of the most common multisite atomistic models. The coarse-grained water is then applied to solvate a typical atomistic protein containing both α-helix and ß-strand elements. Moreover, parallel tempering simulations are performed to investigate the folding free energy landscape of a representative α helical and a ß hairpin structure. For the simulations considered in this work, our dual-resolution method is found to be 3 to 6 times more computationally efficient than corresponding fully atomistic approaches.