Calculation of membrane bending rigidity using field-theoretic umbrella sampling

Determining the Bending Rigidity of Free-Standing Planar Phospholipid Bilayers

We describe a method to determine membrane bending rigidity from capacitance measurements on large area, free-standing, planar, biomembranes. The bending rigidity of lipid membranes is an important biological mechanical property that is commonly optically measured in vesicles, but difficult to quantify in a planar, unsupported system. To accomplish this, we simultaneously image and apply an electric potential to free-standing, millimeter area, planar lipid bilayers composed of DOPC and DOPG phospholipids to measure the membrane Young's (elasticity) modulus. The bilayer is then modeled as two adjacent thin elastic films to calculate bending rigidity from the electromechanical response of the membrane to the applied field. Using DOPC, we show that bending rigidities determined by this approach are in good agreement with the existing work using neutron spin echo on vesicles, atomic force spectroscopy on supported lipid bilayers, and micropipette aspiration of giant unilamellar vesicles. We study the effect of asymmetric calcium concentration on symmetric DOPC and DOPG membranes and quantify the resulting changes in bending rigidity. This platform offers the ability to create planar bilayers of controlled lipid composition and aqueous ionic environment, with the ability to asymmetrically alter both. We aim to leverage this high degree of compositional and environmental control, along with the capacity to measure physical properties, in the study of various biological processes in the future.

Determination of Elastic Parameters of Lipid Membranes with Molecular Dynamics: A Review of Approaches and Theoretical Aspects

Lipid membranes are abundant in living organisms, where they constitute a surrounding shell for cells and their organelles. There are many circumstances in which the deformations of lipid membranes are involved in living cells: fusion and fission, membrane-mediated interaction between membrane inclusions, lipid-protein interaction, formation of pores, etc. In all of these cases, elastic parameters of lipid membranes are important for the description of membrane deformations, as these parameters determine energy barriers and characteristic times of membrane-involved phenomena. Since the development of molecular dynamics (MD), a variety of in silico methods have been proposed for the determination of elastic parameters of simulated lipid membranes. These MD methods allow for the consideration of details unattainable in experimental techniques and represent a distinct scientific field, which is rapidly developing. This work provides a review of these MD approaches with a focus on theoretical aspects. Two main challenges are identified: (i) the ambiguity in the transition from the continuum description of elastic theories to the discrete representation of MD simulations, and (ii) the determination of intrinsic elastic parameters of lipid mixtures, which is complicated due to the composition-curvature coupling effect.

Human Learning for Molecular Simulations: The Collective Variables Dashboard in VMD

The Collective Variables Dashboard is a software tool for real-time, seamless exploration of molecular structures and trajectories in a customizable space of collective variables. The Dashboard arises from the integration of the Collective Variables Module (also known as Colvars) with the visualization software VMD, augmented with a fully discoverable graphical interface offering interactive workflows for the design and analysis of collective variables. Typical use cases include a priori design of collective variables for enhanced sampling and free energy simulations as well as analysis of any type of simulation or collection of structures in a collective variable space. A combination of those cases commonly occurs when preliminary simulations, biased or unbiased, reveal that an optimized set of collective variables is necessary to improve sampling in further simulations. Then the Dashboard provides an efficient way to intuitively explore the space of likely collective variables, validate them on existing data, and use the resulting collective variable definitions directly in further biased simulations using the Collective Variables Module. Visualization of biasing energies and forces is proposed to help analyze or plan biased simulations. We illustrate the use of the Dashboard on two applications: discovering coordinates to describe ligand unbinding from a protein binding site and designing volume-based variables to bias the hydration of a transmembrane pore.

How does curvature affect the free-energy barrier of stalk formation? Small vesicles vs apposing, planar membranes

Using molecular simulations of POPC lipids in conjunction with the calculation of the Minimum Free-Energy Path (MFEP), we study the effect of strong membrane curvature on the formation of the first fusion intermediate—the stalk between a vesicle and its periodic image. We find that the thermodynamic stability of this hourglass-shaped, hydrophobic connection between two vesicles is largely increased by the strong curvature of small vesicles, whereas the intrinsic barrier to form a stalk, i.e., associated with dimple formation and lipid tails protrusions, is similar to the case of two, apposing, planar membranes. A significant reduction of the barrier of stalk formation, however, stems from the lower dehydration free energy that is required to bring highly curved vesicle into a distance, at which stalk formation may occur, compared to the case of apposing, planar membranes.

How proteins open fusion pores: insights from molecular simulations

Fusion proteins can play a versatile and involved role during all stages of the fusion reaction. Their roles go far beyond forcing the opposing membranes into close proximity to drive stalk formation and fusion. Molecular simulations have played a central role in providing a molecular understanding of how fusion proteins actively overcome the free energy barriers of the fusion reaction up to the expansion of the fusion pore. Unexpectedly, molecular simulations have revealed a preference of the biological fusion reaction to proceed through asymmetric pathways resulting in the formation of, e.g., a stalk-hole complex, rim-pore, or vertex pore. Force-field based molecular simulations are now able to directly resolve the minimum free-energy path in protein-mediated fusion as well as quantifying the free energies of formed reaction intermediates. Ongoing developments in Graphics Processing Units (GPUs), free energy calculations, and coarse-grained force-fields will soon gain additional insights into the diverse roles of fusion proteins.

Density Field Thermodynamic Integration (DFTI): A "Soft" Approach to Calculate the Free Energy of Surfactant Self-Assemblies

Thermodynamic integration is one of the most established methods to quantify excess free energies between different metastable states. Excess intermolecular interactions in surfactant assemblies are on the scale of the energy of thermal fluctuations. Therefore, these materials can be deformed and topologically altered via relatively small mechanical stresses. It is thus intuitive to design reaction paths and associated order parameters that exploit the "soft" nature of these materials to mechanically rather than alchemically morph surfactant assemblies from state to state. Here, we propose a novel method coined "density field thermodynamic integration" (DFTI) that adopts the universality and transferability of alchemical methods while simultaneously exploiting the soft excess interactions between surfactant molecules. DFTI was designed for a rapid quantification of the free energy differences between different metastable structures in soft fluid materials. The DFTI method uses an external field coupled to the local density to mechanically morph the system between metastable states of interest. Here, we explored the capability of the DFTI method to swiftly and accurately calculate free energy differences between states. To this aim, we studied two different coarse-grained lipidic surfactant systems: (i) a fusion stalk and (ii) a worm-like micelle. Our results illustrate that DFTI can provide an efficient, versatile, and rather reliable method to calculate the free energy differences between surfactant assemblies.

Direct Derivation of Free Energies of Membrane Deformation and Other Solvent Density Variations From Enhanced Sampling Molecular Dynamics

We report a methodology to calculate the free energy of a shape transformation in a lipid membrane directly from a molecular dynamics simulation. The bilayer need not be homogeneous or symmetric and can be atomically detailed or coarse grained. The method is based on a collective variable that quantifies the similarity between the membrane and a set of predefined density distributions. Enhanced sampling of this "Multi-Map" variable re-shapes the bilayer and permits the derivation of the corresponding potential of mean force. Calculated energies thus reflect the dynamic interplay of atoms and molecules, rather than postulated effects. Evaluation of deformations of different shape, amplitude, and range demonstrates that the macroscopic bending modulus assumed by the Helfrich-Canham model is increasingly unsuitable below the 100-Å scale. In this range of major biological significance, direct free-energy calculations reveal a much greater plasticity. We also quantify the stiffening effect of cholesterol on bilayers of different composition and compare with experiments. Lastly, we illustrate how this approach facilitates analysis of other solvent reorganization processes, such as hydrophobic hydration. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

Thermodynamically reversible paths of the first fusion intermediate reveal an important role for membrane anchors of fusion proteins

Biological membrane fusion proceeds via an essential topological transition of the two membranes involved. Known players such as certain lipid species and fusion proteins are generally believed to alter the free energy and thus the rate of the fusion reaction. Quantifying these effects by theory poses a major challenge since the essential reaction intermediates are collective, diffusive and of a molecular length scale. We conducted molecular dynamics simulations in conjunction with a state-of-the-art string method to resolve the minimum free-energy path of the first fusion intermediate state, the so-called stalk. We demonstrate that the isolated transmembrane domains (TMDs) of fusion proteins such as SNARE molecules drastically lower the free energy of both the stalk barrier and metastable stalk, which is not trivially explained by molecular shape arguments. We relate this effect to the local thinning of the membrane (negative hydrophobic mismatch) imposed by the TMDs which favors the nearby presence of the highly bent stalk structure or prestalk dimple. The distance between the membranes is the most crucial determinant of the free energy of the stalk, whereas the free-energy barrier changes only slightly. Surprisingly, fusion enhancing lipids, i.e., lipids with a negative spontaneous curvature, such as PE lipids have little effect on the free energy of the stalk barrier, likely because of its single molecular nature. In contrast, the lipid shape plays a crucial role in overcoming the hydration repulsion between two membranes and thus rather lowers the total work required to form a stalk.

Numerical algorithms for solving self-consistent field theory reversely for block copolymer systems

Besides dictating the equilibrium phase diagram, the rugged free-energy landscape of AB block copolymers gives rise to a multitude of non-equilibrium phenomena. Self-consistent field theory (SCFT) can be employed to calculate the mean-field free energy, F [ ϕ A t a r g e t ] , of a non-equilibrium unstable state that is characterized by a given spatial density distribution, ϕ A t a r g e t , in the incompressible system. Such a free-energy functional is the basis of describing the structure formation by dynamic SCFT techniques or the identification of minimum free-energy paths via the string method. The crucial step consists in computing the external potential fields that generate the given density distribution in the corresponding system of non-interacting copolymers, i.e., the potential-to-density relation employed in equilibrium SCFT calculations has to be inverted (reverse SCFT calculation). We describe, generalize, and evaluate the computational efficiency of two different numerical algorithms for this reverse SCFT calculation-the Debye-function algorithm based on the structure factor and the field-theoretic umbrella-potential (FUP) algorithm. In contrast to the Debye-function algorithm, the FUP algorithm only yields the exact mean-field values of the given target densities in the limit of a strong umbrella potential, and we devise a two-step variant of the FUP algorithm that significantly mitigates this issue. For Gaussian copolymers, the Debye-function algorithm is more efficient for highly unstable states that are far away from the equilibrium, whereas the improved FUP algorithm outperforms the Debye-function algorithm closer to metastable states and is easily transferred to more complex molecular architectures.

Bridging molecular simulation models and elastic theories for amphiphilic membranes

The Single Chain Mean Field theory is used to link coarse-grained models of amphiphilic molecules with analytical models for membrane elasticity, where phenomenological parameters are deduced from explicit molecular models and force fields. We estimate the elastic constants based on the free energy of the amphiphilic bilayer in planar and cylindrical geometries on the example of four amphiphilic molecules that differ in length and stiffness. We study how these variations affect the equilibrium bilayer structure, the equilibrium free energy, and the elastic constants. Bending rigidities are obtained within the typical range of experimental values for phospholipid membranes in a liquid state.

Bending Lipid Bilayers: A Closed-Form Collective Variable for Effective Free-Energy Landscapes in Quantitative Biology

Curvature-related processes are of major importance during protein-membrane interactions. The illusive simplicity of membrane reshaping masks a complex molecular process crucial for a wide range of biological functions like fusion, endo- and exocytosis, cell division, cytokinesis, and autophagy. To date, no functional expression of a reaction coordinate capable of biasing molecular dynamics simulations to produce membrane curvature has been reported. This represents a major drawback given that the adequate identification of proper collective variables to enhance sampling is fundamental for restrained dynamics techniques. In this work, we present a closed-form equation of a collective variable that induces bending in lipid bilayers in a controlled manner, allowing for straightforward calculation of free energy landscapes of important curvature-related events, using standard methods such as umbrella sampling and metadynamics. As a direct application of the collective variable, we calculate the bending free energies of a ternary lipid bilayer in the presence and the absence of a Bin/Amphiphysin/Rvs domain with an N-terminal amphipathic helix (N-BAR), a well-known peripheral membrane protein that induces curvature.

Exploiting Lipid Permutation Symmetry to Compute Membrane Remodeling Free Energies

A complete physical description of membrane remodeling processes, such as fusion or fission, requires knowledge of the underlying free energy landscapes, particularly in barrier regions involving collective shape changes, topological transitions, and high curvature, where Canham-Helfrich (CH) continuum descriptions may fail. To calculate these free energies using atomistic simulations, one must address not only the sampling problem due to high free energy barriers, but also an orthogonal sampling problem of combinatorial complexity stemming from the permutation symmetry of identical lipids. Here, we solve the combinatorial problem with a permutation reduction scheme to map a structural ensemble into a compact, nondegenerate subregion of configuration space, thereby permitting straightforward free energy calculations via umbrella sampling. We applied this approach, using a coarse-grained lipid model, to test the CH description of bending and found sharp increases in the bending modulus for curvature radii below 10 nm. These deviations suggest that an anharmonic bending term may be required for CH models to give quantitative energetics of highly curved states.