Dataset Information

Surface viscosities of lipid bilayers determined from equilibrium molecular dynamics simulations.

ABSTRACT: Lipid membrane viscosity is critical to biological function. Bacterial cells grown in different environments alter their lipid composition in order to maintain a specific viscosity, and membrane viscosity has been linked to the rate of cellular respiration. To understand the factors that determine the viscosity of a membrane, we ran equilibrium all-atom simulations of single component lipid bilayers and calculated their viscosities. The viscosity was calculated via a Green-Kubo relation, with the stress-tensor autocorrelation function modeled by a stretched exponential function. By simulating a series of lipids at different temperatures, we establish the dependence of viscosity on several aspects of lipid chemistry, including hydrocarbon chain length, unsaturation, and backbone structure. Sphingomyelin is found to have a remarkably high viscosity, roughly 20 times that of DPPC. Furthermore, we find that inclusion of the entire range of the dispersion interaction increases viscosity by up to 140%. The simulated viscosities are similar to experimental values obtained from the rotational dynamics of small chromophores and from the diffusion of integral membrane proteins but significantly lower than recent measurements based on the deformation of giant vesicles.

SUBMITTER: Fitzgerald JE

PROVIDER: S-EPMC10111272 | biostudies-literature | 2023 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Surface viscosities of lipid bilayers determined from equilibrium molecular dynamics simulations.

Fitzgerald James E JE Venable Richard M RM Pastor Richard W RW Lyman Edward R ER

Biophysical journal 20230204 6

Lipid membrane viscosity is critical to biological function. Bacterial cells grown in different environments alter their lipid composition in order to maintain a specific viscosity, and membrane viscosity has been linked to the rate of cellular respiration. To understand the factors that determine the viscosity of a membrane, we ran equilibrium all-atom simulations of single component lipid bilayers and calculated their viscosities. The viscosity was calculated via a Green-Kubo relation, with th ...[more]

PMID: 36739477

Similar Datasets

Project description:Human β defensin type 3 (hBD-3) is a small cationic cysteine-rich peptide. It has a broad spectrum of antimicrobial activities. However, at high concentrations, it also shows hemolytic activity by interrupting red blood cells. To understand the selectivity of hBD-3 disrupting cell membranes, investigating the capability of hBD-3 translocating through different membranes is important. Since hBD-3 in the analogue form in which all three pairs of disulfide bonds are broken has similar antibacterial activities to the wild-type, this project investigates the structure and dynamics of an hBD-3 analogue in monomer, dimer, and tetramer forms through both zwitterionic and negatively charged lipid bilayers using molecular dynamics (MD) simulations. One tetramer structure of hBD-3 was predicted by running all-atom MD simulations on hBD-3 in water at a high concentration, which was found to be stable in water during 400 ns all-atom simulations based on root-mean-squared deviation, root-mean-squared fluctuation, buried surface area, and binding interaction energy calculations. After that, hBD-3 in different forms was placed inside different membranes, and then steered MD simulation was conducted to pull the hBD-3 out of the membrane along the z-direction to generate different configurational windows to set up umbrella-sampling (US) simulations. Because extensive sampling is important to obtain accurate free energy barriers, coarse-grained US MD simulations were performed in each window. Based on the long-term simulation result, membrane thinning was found near hBD-3 in different lipid bilayers and in different hBD-3 oligomer systems. By calculating the root-mean-squared deviation of the z-coordinate of hBD-3 molecules, rotation of the oligomer inside the bilayer and stretching of the oligomer structure along the z-direction were observed. Although reorientation of lipid heads toward the hBD-3 tetramer was observed based on the density profile calculation, the order parameter calculation shows that hBD-3 disrupts 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS) lipids more significantly and makes it less ordered than on 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipids. Calculating the free energy of hBD-3 through different lipid bilayers, it was found that generally hBD-3 encounters a lower energy barrier through negatively charged lipid membranes than the zwitterionic membrane. hBD-3 in different forms needs to overcome a lower energy barrier crossing the combined POPC+POPS bilayer through the POPS leaflet than through the POPC leaflet. Besides that, the potential of mean force result suggests that hBD-3 forms an oligomer translocating negatively charged lipid membranes at a low concentration. This study supplied new insight into the antibacterial mechanism of hBD-3 through different membranes.

Dataset Information

Surface viscosities of lipid bilayers determined from equilibrium molecular dynamics simulations.

Publications

Surface viscosities of lipid bilayers determined from equilibrium molecular dynamics simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets