Molecular Dynamics Simulations of Membrane Permeability

Molecular Picture of the Transient Nature of Lipid Rafts

In biological membranes, lipid rafts are now thought to be transient and nanoscopic. However, the mechanism responsible for these nanoscopic assemblies remains poorly understood, even in the case of model membranes. As a result, it has proven extremely challenging to probe the physicochemical properties of lipid rafts at the molecular level. Here, we use all-atom molecular dynamics (MD) simulations and inelastic X-ray scattering (IXS), an intrinsically nanoscale technique, to directly probe the energy transfer and collective short-wavelength dynamics (phonons) of biologically relevant model membranes. We show that the nanoscale propagation of stress in lipid rafts takes place in the form of collective motions made up of longitudinal (compression waves) and transverse (shear waves) molecular vibrations. Importantly, we provide a molecular picture for the so-called van der Waals mediated "force from lipid" [Anishkin, A. et al. Proc. Natl. Acad. Sci. U.S.A. 2014, 111, 7898], a key parameter for the ionic channel mechano-transduction and the mechanism for the lipid transfer of molecular level stress [Aponte-Santamarı́a, C. et al. J. Am. Chem. Soc. 2017, 139, 13588]. Specifically, we describe how lipid rafts are formed and maintained through the propagation of molecular stress, lipid raft rattling dynamics, and a relaxation process. Eventually, the rafts dissipate through the self-diffusion of lipids making up the rafts. We also show that the molecular stress and viscoelastic properties of transient lipid rafts can be modulated through the use of hydrophobic biomolecules such as melatonin and tryptophan. Ultimately, the herein proposed mechanism describing the molecular interactions for the formation and dissolution of lipid rafts may offer insights as to how lipid rafts enable biological function.

Multi-phase Boltzmann weighting: accounting for local inhomogeneity in molecular simulations of water-octanol partition coefficients in the SAMPL6 challenge

Accurately computing partition coefficients is a pivotal part of drug discovery. Specifically, octanol-water partition coefficients can provide information into hydrophobicity of drug-like molecules, as well as a de facto representation of membrane permeability. However, one challenge facing the computation of partition coefficients is the need to encapsulate various microscopic environments. These include areas of largely bulk solvent (i.e., either water or octanol) or regions where octanol is saturated with water or areas of higher salt concentration. Also, tautomeric effects require consideration. Thus, we present a Boltzmann weighting approach that incorporates transfer free energies across varying microscopic media, as well as varying tautomeric state, to compute partition coefficients in the SAMPL6 challenge.

Design of Active Interfaces Using Responsive Molecular Components

Responsive interfaces are interfaces that show a defined and reversible change in physical properties in response to external stimuli. Typically, responsive interfaces result from the immobilization of responsive molecular components at the interface that translate a nanoscale signal into a macroscopic effect. Responsive interfaces can also be obtained if the topology of the interface can be reversibly changed using an external stimulus. As the surface of any material is its connection to the environment, responsive interfaces provide opportunities for interactive materials which are not only able to change properties upon demand, but also sense their environment and act autonomously. The application of responsive molecular components at interfaces, however, requires chemical and physical compatibility with the material surface of interest, posing a challenge not least in the retention of the responsive functionality. The state of the art in "active" interfaces which display responsive wettability, permeability, or adhesion is discussed, with a particular emphasis on microscale and nanoscale patterning since patterned interfaces can give rise to unique material properties. Finally, perspectives in the development of responsive interfaces, as well as promising approaches for bypassing the most prominent challenges are discussed.

Computational Approaches in Membrane Science and Engineering

Computational modelling and simulation form a consolidated branch in the multidisciplinary field of membrane science and technology [...]

Interaction of Camptothecin with Model Cellular Membranes

Accurate and efficient prediction of drug partitioning in model membranes is of significant interest to the pharmaceutical industry. Herein, we utilize advanced sampling methods, specifically, the adaptive biasing force methodology to calculate the potential of mean force for a model hydrophobic anticancer drug, camptothecin (CPT), across three model interfaces. We consider an octanol bilayer, a thick octanol/water interface, and a model 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/water interface. We characterize the enthalpic and entropic contributions of the drug to the potential of mean force. We show that the rotational entropy of the drug is inversely related to the probability of hydrogen bond formation of the drug with the POPC membrane. In addition, in long-time microsecond simulations of a high concentration of CPT above the POPC membrane, we show that strong drug-drug aromatic interactions shift the spatial orientation of the drug with the membrane. Stacks of hydrophobic drugs form, allowing penetration of the drug just under the POPC head groups. These results imply that inhomogeneous membrane models need to take into account the effect of drug aggregation on the membrane environment.

Understanding and modelling the interactions of peptides with membranes: from partitioning to self-assembly

Atomic detail simulations are starting to reveal how flexible polypeptides interact with fluid lipid bilayers. These insights are transforming our understanding of one of the fundamental processes in biology: membrane protein folding and assembly. Advanced molecular dynamics (MD) simulation techniques enable accurate prediction of protein structure, folding pathways and assembly in microsecond-timescales. Such simulations show how membrane-active peptides self-assemble in cell membranes, revealing their binding, folding, insertion, and aggregation, while at the same time providing atomic resolution details of peptide-lipid interactions. Essential to the impact of simulations are experimental approaches that enable calibration and validation of the computational models and techniques. In this review, we summarize the current development of applying unbiased atomic detail MD simulations and the relation to experimental techniques, to study peptide folding and provide our perspective of the field.

Lipidomic and biophysical homeostasis of mammalian membranes counteracts dietary lipid perturbations to maintain cellular fitness

Proper membrane physiology requires maintenance of biophysical properties, which must be buffered from external perturbations. While homeostatic adaptation of membrane fluidity to temperature variation is a ubiquitous feature of ectothermic organisms, such responsive membrane adaptation to external inputs has not been directly observed in mammals. Here, we report that challenging mammalian membranes by dietary lipids leads to robust lipidomic remodeling to preserve membrane physical properties. Specifically, exogenous polyunsaturated fatty acids are rapidly incorporated into membrane lipids, inducing a reduction in membrane packing. These effects are rapidly compensated both in culture and in vivo by lipidome-wide remodeling, most notably upregulation of saturated lipids and cholesterol, resulting in recovery of membrane packing and permeability. Abrogation of this response results in cytotoxicity when membrane homeostasis is challenged by dietary lipids. These results reveal an essential mammalian mechanism for membrane homeostasis wherein lipidome remodeling in response to dietary lipid inputs preserves functional membrane phenotypes.

Proper membrane physiology requires maintenance of a narrow range of physicochemical properties, which must be buffered from external perturbations. Here, authors report lipidomic remodeling to preserve membrane physical properties upon exogenous polyunsaturated fatty acids exposure.

Optimization of hydrophobic nanoparticles to better target lipid rafts with molecular dynamics simulations

Due to different interactions between lipids and proteins, a plasma membrane can segregate into different membrane domains. Among them, ordered functional membrane domains are defined as "lipid rafts", which play key roles in many biological processes (e.g., signal transduction, endocytosis, etc.) in the cell. Hence, it will be of much biological significance to monitor and even regulate the dynamics of lipid rafts. In this work, we designed a ligand-modified spherical nanoparticle with coarse-grained molecular dynamics simulations, which can be encapsulated into the hydrophobic region of the lipid membrane and specifically target either raft or non-raft membrane domains. The preferred localization of the nanoparticle can be tuned by adjusting ligand hydrophobicity, length and density. Generally, more hydrophobic nanoparticles tend to target the raft domain, while less hydrophobic nanoparticles prefer the non-raft domain. Besides, ligand length and density jointly determine the exposure of nanoparticle cores and thus affect the roles of ligands in nanoparticles' final localization. Our results may provide insights into the experimental design of functional nanoparticles, targeting the lipid raft and regulating its dynamics.

Prediction of Passive Membrane Permeability by Semi-Empirical Method Considering Viscous and Inertial Resistances and Different Rates of Conformational Change and Diffusion

Membrane permeability is an important property of drugs in adsorption. Many prediction methods work well for small molecules, but the prediction of middle-molecule permeability is still difficult. In the present study, we modified a classical permeability model based on Fick's law to study passive membrane permeability. The model consisted of the distribution of solute from water to membrane and the diffusion of solute in each solvent. The diffusion coefficient is the inverse of the resistance, and we examined the inertial resistance in addition to the viscous resistance, the latter of which has been widely used in permeability prediction. Also, we examined three models changing the balance between the diffusion of solute in membrane and the conformational change of solute. The inertial resistance improved the prediction results in addition to the viscous resistance. The models worked well not only for small molecules but also for middle molecules, whose structures have more conformational freedom.

Permeability of membranes in the liquid ordered and liquid disordered phases

The functional significance of ordered nanodomains (or rafts) in cholesterol rich eukaryotic cell membranes has only begun to be explored. This study exploits the correspondence of cellular rafts and liquid ordered (Lo) phases of three-component lipid bilayers to examine permeability. Molecular dynamics simulations of Lo phase dipalmitoylphosphatidylcholine (DPPC), dioleoylphosphatidylcholine (DOPC), and cholesterol show that oxygen and water transit a leaflet through the DOPC and cholesterol rich boundaries of hexagonally packed DPPC microdomains, freely diffuse along the bilayer midplane, and escape the membrane along the boundary regions. Electron paramagnetic resonance experiments provide critical validation: the measured ratio of oxygen concentrations near the midplanes of liquid disordered (Ld) and Lo bilayers of DPPC/DOPC/cholesterol is 1.75 ± 0.35, in very good agreement with 1.3 ± 0.3 obtained from simulation. The results show how cellular rafts can be structurally rigid signaling platforms while remaining nearly as permeable to small molecules as the Ld phase.

Accurate Biomolecular Simulations Account for Electronic Polarization

In this perspective, we discuss where and how accounting for electronic many-body polarization affects the accuracy of classical molecular dynamics simulations of biomolecules. While the effects of electronic polarization are highly pronounced for molecules with an opposite total charge, they are also non-negligible for interactions with overall neutral molecules. For instance, neglecting these effects in important biomolecules like amino acids and phospholipids affects the structure of proteins and membranes having a large impact on interpreting experimental data as well as building coarse grained models. With the combined advances in theory, algorithms and computational power it is currently realistic to perform simulations with explicit polarizable dipoles on systems with relevant sizes and complexity. Alternatively, the effects of electronic polarization can also be included at zero additional computational cost compared to standard fixed-charge force fields using the electronic continuum correction, as was recently demonstrated for several classes of biomolecules.

Quantitative Characterization of Protein-Lipid Interactions by Free Energy Simulation between Binary Bilayers

Using a recently developed binary bilayer system (BBS) consisting of two patches of laterally contacting bilayers, umbrella sampling molecular dynamics (MD) simulations were performed for quantitative characterization of protein-lipid interactions. The BBS is composed of 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) with an embedded model membrane protein, a gramicidin A (gA) channel. The calculated free energy difference for the transfer of a gA channel from DLPC (hydrophobic thickness ≈ 21.5 Å) to DMPC (hydrophobic thickness ≈ 25.5 Å) bilayers, ΔG(DLPC → DMPC), is -2.2 ± 0.7 kcal/mol. This value appears at odds with the traditional view that the hydrophobic length of the gA channel is ∼22 Å. To understand this discrepancy, we first note that recent MD simulations by different groups have shown that lipid bilayer thickness profiles in the vicinity of a gA channel differ qualitatively from the deformation profile predicted from continuum elastic bilayer models. Our MD simulations at low and high gA:lipid molar ratios and different membrane compositions indicate that the gA channel's effective hydrophobic length is ∼26 Å. Using this effective hydrophobic length, ΔG(DLPC → DMPC) determined here is in excellent agreement with predictions based on continuum elastic models (-3.0 to -2.2 kcal/mol) where the bilayer deformation energy is approximated as a harmonic function of the mismatch between the channel's effective hydrophobic length and the hydrophobic thickness of the bilayer. The free energy profile for gA in the BBS includes a barrier at the interface between the two bilayers which can be attributed to the line tension at the interface between two bilayers with different hydrophobic thicknesses. This observation implies that translation of a peptide between two different regions of a cell membrane (such as between the liquid ordered and disordered phases) may include effects of a barrier at the interface in addition to the relative free energies of the species far from the interface. The BBS allows for direct transfer free energy calculations between bilayers without a need of a reference medium, such as bulk water, and thus provides an efficient simulation protocol for the quantitative characterization of protein-lipid interactions at all-atom resolution.

Dynamic Protonation Dramatically Affects the Membrane Permeability of Drug-like Molecules

Permeability (P_m) across biological membranes is of fundamental importance and a key factor in drug absorption, distribution, and development. Although the majority of drugs will be charged at some point during oral delivery, our understanding of membrane permeation by charged species is limited. The canonical model assumes that only neutral molecules partition into and passively permeate across membranes, but there is mounting evidence that these processes are also facile for certain charged species. However, it is unknown whether such ionizable permeants dynamically neutralize at the membrane surface or permeate in their charged form. To probe protonation-coupled permeation in atomic detail, we herein apply continuous constant-pH molecular dynamics along with free energy sampling to study the permeation of a weak base propranolol (PPL), and evaluate the impact of including dynamic protonation on P_m. The simulations reveal that PPL dynamically neutralizes at the lipid-tail interface, which dramatically influences the permeation free energy landscape and explains why the conventional model overestimates the assigned intrinsic permeability. We demonstrate how fixed-charge-state simulations can account for this effect, and propose a revised model that better describes pH-coupled partitioning and permeation. Our results demonstrate how dynamic changes in protonation state may play a critical role in the permeation of ionizable molecules, including pharmaceuticals and drug-like molecules, thus requiring a revision of the standard picture.

Ionizable amino lipid interactions with POPC: implications for lipid nanoparticle function

Lipid nanoparticles (LNPs) composed of ionizable cationic lipids are currently the leading systems for siRNA delivery in liver disease, with the major limitation of low siRNA release efficacy into the cytoplasm. Ionizable cationic lipids are known to be of critical importance in LNP structure and stability, siRNA entrapment, and endosomal disruption. However, their distribution inside the LNPs and their exact role in cytoplasmic delivery remain unclear. A recent study [Kulkarni et al., On the formation and morphology of lipid nanoparticles containing ionizable cationic lipids and siRNA, ACS Nano, 2018, 12(5), 4787-4795] on LNP-siRNA systems containing the ionizable lipid DLin-KC2-DMA (also known as KC2 with an apparent pKa of ca. 6.7) suggested that neutral KC2 segregates from other components and forms an amorphous oil droplet in the core of LNPs. In this paper, we present evidence supporting the model proposed by Kulkarni et al. We studied KC2 segregation in the presence of POPC using molecular dynamics simulation, deuterium NMR, SAXS, and cryo-TEM experiments, and found that neutral KC2 has a high tendency to separate from POPC dispersions. KC2 confinement, upon raising the pH during the formulation process, could result in rearrangement of the internal structure of LNPs. As interactions between cationic KC2 and anionic endosomal lipids are thought to be a key factor in cargo release, KC2 confinement inside the LNP may be responsible for the observed low release efficacy.

Glycyrrhizin-Assisted Transport of Praziquantel Anthelmintic Drug through the Lipid Membrane: An Experiment and MD Simulation

Praziquantel (PZQ) is one of the most widespread anthelmintic drugs. However, the frequent insufficient application of PZQ after oral administration is associated with its low solubility, penetration rate, and bioavailability. In the present study, the permeation of PZQ through a 1,2-dioleoyl- sn-glycero-3-phosphocholine (DOPC) membrane was investigated to probe glycyrrhizin-assisted transport. Glycyrrhizin (or glycyrrhizic acid, GA), a natural saponin, shows the ability to enhance the therapeutic activity of various drugs when it is used as a drug delivery system. However, the molecular mechanism of this effect is still under debate. In the present study, the transport rate was measured experimentally by a parallel artificial membrane permeation assay (PAMPA) and molecular dynamics (MD) simulation with DOPC lipid bilayers. The formation of the noncovalent supramolecular complex of PZQ with disodium salt of GA (Na2GA) in an aqueous solution was proved by the NMR relaxation technique. PAMPA experiments show a strong increase in the amount of the penetrating praziquantel molecules in comparison with a saturated aqueous solution of pure drug used as a control. MD simulation of PZQ penetration through the bilayer demonstrates an increase in permeability into the membrane in the presence of a glycyrrhizin molecule. A decrease in the free energy barrier in the middle of the lipid bilayer was obtained, associated with the hydrogen bond between PZQ and GA. Also, GA reduces the local bilayer surface resistance to penetration of PZQ by rearranging the surface lipid headgroups. This study clarifies the mechanism of increasing the drug's bioavailability in the presence of glycyrrhizin.

Interactions of Water and Alkanes: Modifying Additive Force Fields to Account for Polarization Effects

Atomistic biomolecular simulations predominantly utilize additive force fields (FF), where the electrostatic potential is modeled by fixed point charges. Among other consequences, the lack of polarizability in these models undermines the balance of hydrophilic/hydrophobic nonbonded interactions. Simulations of water/alkane systems using the TIP3P water model and CHARMM36 parameters reveal a 1 kcal/mol overestimate of the experimental transfer free energy of water to hexadecane; more recent optimized water models (SPC/E, TIP4P/2005, TIP4P-Ew, TIP3P-FB, TIP4P-FB, OPC, TIP4P-D) overestimate this transfer free energy by approximately 2 kcal/mol. In contrast, the polarizable SWM4-NDP and SWM6 water models reproduce experimental values to within statistical error. As an alternative to explicitly modeling polarizability, this paper develops an efficient automated workflow to optimize pair-specific Lennard-Jones parameters within an additive FF. Water/hexadecane is used as a prototype and the free energy of water transfer to hexadecane as a target. The optimized model yields quantitative agreement with the experimental transfer free energy and improves the water/hexadecane interfacial tension by 6%. Simulations of five different lipid bilayers show a strong increase of water permeabilities compared to the unmodified CHARMM36 lipid FF which consistently improves match with experiment: the order-of-magnitude underestimate for monounsaturated bilayers is rectified and the factor of 2.8-4 underestimate for saturated bilayers is turned into a factor of 1.5-3 overestimate. While agreement with experiment is decreased for the diffusion constant of water in hexadecane, alkane transfer free energies, and the bilayers' area per lipid, the method provides a permeant-specific route to achieve a wide range of heterogeneous observables via rapidly optimized pairwise parameters.

Reactions at Biomembrane Interfaces

Membranes surrounding the biological cell and its internal compartments host proteins that catalyze chemical reactions essential for the functioning of the cell. Rather than being a passive structural matrix that holds membrane-embedded proteins in place, the membrane can largely shape the conformational energy landscape of membrane proteins and impact the energetics of their chemical reaction. Here, we highlight the challenges in understanding how lipids impact the conformational energy landscape of macromolecular membrane complexes whose functioning involves chemical reactions including proton transfer. We review here advances in our understanding of how chemical reactions occur at membrane interfaces gleaned with both theoretical and experimental advances using simple protein systems as guides. Our perspective is that of bridging experiments with theory to understand general physicochemical principles of membrane reactions, with a long term goal of furthering our understanding of the role of the lipids on the functioning of complex macromolecular assemblies at the membrane interface.