Dataset Information

Efficient Simulations of Solvent Asymmetry Across Lipid Membranes Using Flat-Bottom Restraints.

ABSTRACT: The routinely employed periodic boundary conditions complicate molecular simulations of physiologically relevant asymmetric lipid membranes together with their distinct solvent environments. Therefore, separating the extracellular fluid from its cytosolic counterpart has often been performed using a costly double-bilayer setup. Here, we demonstrate that the lipid membrane and solvent asymmetry can be efficiently modeled with a single lipid bilayer by applying an inverted flat-bottom potential to ions and other solute molecules, thereby restraining them to only interact with the relevant leaflet. We carefully optimized the parameters of the suggested method so that the results obtained using the flat-bottom and double-bilayer approaches become mutually indistinguishable. Then, we apply the flat-bottom approach to lipid bilayers with various compositions and solvent environments, covering ions and cationic peptides to validate the approach in a realistic use case. We also discuss the possible limitations of the method as well as its computational efficiency and provide a step-by-step guide on how to set up such simulations in a straightforward manner.

SUBMITTER: Biriukov D

PROVIDER: S-EPMC10537000 | biostudies-literature | 2023 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Efficient Simulations of Solvent Asymmetry Across Lipid Membranes Using Flat-Bottom Restraints.

Biriukov Denys D Javanainen Matti M

Journal of chemical theory and computation 20230831 18

The routinely employed periodic boundary conditions complicate molecular simulations of physiologically relevant asymmetric lipid membranes together with their distinct solvent environments. Therefore, separating the extracellular fluid from its cytosolic counterpart has often been performed using a costly double-bilayer setup. Here, we demonstrate that the lipid membrane and solvent asymmetry can be efficiently modeled with a single lipid bilayer by applying an inverted flat-bottom potential to ...[more]

PMID: 37651714

Similar Datasets

Project description:Despite considerable advances in computing power, atomistic simulations under nonperiodic boundary conditions, with Coulombic electrostatic interactions and in systems large enough to reduce finite-size associated errors in thermodynamic quantities to within the thermal energy, are still not affordable. As a result, periodic boundary conditions, systems of microscopic size and effective electrostatic interaction functions are frequently resorted to. Ensuing artifacts in thermodynamic quantities are nowadays routinely corrected a posteriori, but the underlying configurational sampling still descends from spurious forces. The present study addresses this problem through the introduction of on-the-fly corrections to the physical forces during an atomistic molecular dynamics simulation. Two different approaches are suggested, where the force corrections are derived from special potential energy terms. In the first approach, the solvent-generated electrostatic potential sampled at a given atom site is restrained to a target value involving corrections for electrostatic artifacts. In the second approach, the long-range regime of the solvent polarization around a given atom site is restrained to the Born polarization, i.e., the solvent polarization corresponding to the ideal situation of a macroscopic system under nonperiodic boundary conditions and governed by Coulombic electrostatic interactions. The restraints are applied to the explicit-water simulation of a hydrated sodium ion, and the effect of the restraints on the structural and energetic properties of the solvent is illustrated. Furthermore, by means of the calculation of the charging free energy of a hydrated sodium ion, it is shown how the electrostatic potential restraint translates into the on-the-fly consideration of the corresponding free-energy correction terms. It is discussed how the restraints can be generalized to situations involving several solute particles. Although the present study considers a very simple system only, it is an important step toward the on-the-fly elimination of finite-size and approximate-electrostatic artifacts during atomistic molecular dynamics simulations.

Project description:Interfacing nanoelectronic devices with cell membranes can enable multiplexed detection of fundamental biological processes (such as signal transduction, electrophysiology, and import/export control) even down to the single ion channel level, which can lead to a variety of applications in pharmacology and clinical diagnosis. Therefore, it is necessary to understand and control the chemical and electrical interface between the device and the lipid bilayer membrane. Here, we develop a simple bottom-up approach to assemble tethered bilayer lipid membranes (tBLMs) on silicon wafers and glass slides, using a covalent tether attachment chemistry based on silane functionalization, followed by step-by-step stacking of two other functional molecular building blocks (oligo-poly(ethylene glycol) (PEG) and lipid). A standard vesicle fusion process was used to complete the bilayer formation. The monolayer synthetic scheme includes three well-established chemical reactions: self-assembly, epoxy-amine reaction, and EDC/NHS cross-linking reaction. All three reactions are facile and simple and can be easily implemented in many research labs, on the basis of common, commercially available precursors using mild reaction conditions. The oligo-PEG acts as the hydrophilic spacer, a key role in the formation of a homogeneous bilayer membrane. To explore the broad applicability of this approach, we have further demonstrated the formation of tBLMs on three common classes of (nano)electronic biosensor devices: indium-tin oxide-coated glass, silicon nanoribbon devices, and high-density single-walled carbon nanotubes (SWNT) networks on glass. More importantly, we incorporated alemethicin into tBLMs and realized the real-time recording of single ion channel activity with high sensitivity and high temporal resolution using the tBLMs/SWNT network transistor hybrid platform. This approach can provide a covalently bonded lipid coating on the oxide layer of nanoelectronic devices, which will enable a variety of applications in the emerging field of nanoelectronic interfaces to electrophysiology.

Dataset Information

Efficient Simulations of Solvent Asymmetry Across Lipid Membranes Using Flat-Bottom Restraints.

Publications

Efficient Simulations of Solvent Asymmetry Across Lipid Membranes Using Flat-Bottom Restraints.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets