Convergence of Free Energy Profile of Coumarin in Lipid Bilayer

Protonation State of a Bioactive Compound Regulates Its Release from Lamellar Gel-Phase Bilayers

Lamellar gel networks (LGNs) in personal care or pharmaceutical lotions and creams provide an opaque cream appearance and a creamy texture to these products. Within the LGNs, the lamellar gel (Lβ) phase composed of regularly spaced bilayers of surfactants and long-chain fatty alcohols is predominately responsible for the unique rheological properties of the LGNs. To extend the shelf life of LGN-containing products, bioactive compounds with antimicrobial properties are often incorporated into the formulation. However, how the protonation state of the bioactive compounds regulates their release from the Lβ-phase bilayers is currently unknown. Using molecular dynamics simulations, we found that the protonated (neutral) form of cinnamic acid, a common antimicrobial food additive, has a retention ratio higher than that of its deprotonated (charged) counterpart in the Lβ-phase bilayer. From free energy calculations, we determined that not only is the protonated molecule more stable in the hydrophobic interior of the bilayer but also the formation of hydrogen-bonded dimers significantly enhances its stability within the bilayer. Thus, the protonation state has a profound impact on bioavailability of the compounds. Our results also highlight the importance of considering possible oligomeric states of molecules when performing calculations to estimate the permeability of molecules within various bilayers.

Permeation of a Homologous Series of NBD-Labeled Fatty Amines through Lipid Bilayers: A Molecular Dynamics Study

Permeation through biomembranes is ubiquitous for drugs to reach their active sites. Asymmetry of the cell plasma membrane (PM) has been described as having an important role in this process. Here we describe the interaction of a homologous series of 7-nitrobenz-2-oxa-1,3-diazol-4-yl (NBD)-labeled amphiphiles (NBD-Cn, n = 4 to 16) with lipid bilayers of different compositions (1-palmitoyl, 2-oleoyl-sn-glycero-3-phosphocholine (POPC):cholesterol (1:1) and palmitoylated sphingomyelin (SpM):cholesterol (6:4)), including an asymmetric bilayer. Both unrestrained and umbrella sampling (US) simulations (at varying distances to the bilayer center) were carried out. The free energy profile of NBD-Cn at different depths in the membrane was obtained from the US simulations. The behavior of the amphiphiles during the permeation process was described regarding their orientation, chain elongation, and H-bonding to lipid and water molecules. Permeability coefficients were also calculated for the different amphiphiles of the series, using the inhomogeneous solubility-diffusion model (ISDM). Quantitative agreement with values obtained from kinetic modeling of the permeation process could not be obtained. However, for the longer, and more hydrophobic amphiphiles, the variation trend along the homologous series was qualitatively better matched by the ISDM when the equilibrium location of each amphiphile was taken as reference (ΔG = 0), compared to the usual choice of bulk water.

Multiscale Modeling Unravels the Influence of Biomembranes on the Photochemical Properties of Embedded Anti-Oxidative Polyphenolic and Phenanthroline Chelating Dyes

The embedding of caffeate methyl ester, the flavonoids luteolin and quercetin, and the o-phenanthroline and neocuproine in a liquid disordered lipid bilayer has been studied through extensive atomistic calculations. The location and the orientation of these bio-active antioxidants are explained and analyzed. While the two phenanthrolines strongly associate with the lipid tail region, the other three compounds are rather found among the head groups. The simulations showcase conformational changes of the flavonoids. Through the use of a hybrid quantum mechanics-molecular mechanics scheme and supported by a profound benchmarking of the electronic excited-state method for these compounds, the influence of the anisotropic environment on the compounds' optical properties is analyzed. Influences of surrounding water molecules and of the polar parts of the lipids on the transition dipole moments and excited-state dipole moments are weighted with respect to a change in conformation. The current study highlights the importance of the mapping of molecular interactions in model membranes and pinpoints properties, which can be biomedically used to discriminate and detect different lipid environments.

Interactions between Rhodamine Dyes and Model Membrane Systems—Insights from Molecular Dynamics Simulations

Background: rhodamines are dyes widely used as fluorescent tags in cell imaging, probing of mitochondrial membrane potential, and as P-glycoprotein model substrates. In all these applications, detailed understanding of the interaction between rhodamines and biomembranes is fundamental. Methods: we combined atomistic molecular dynamics (MD) simulations and fluorescence spectroscopy to characterize the interaction between rhodamines 123 and B (Rh123 and RhB, respectively) and POPC bilayers. Results: while the xanthene moiety orients roughly parallel to the membrane plane in unrestrained MD simulations, variations on the relative position of the benzoic ring (below the xanthene for Rh123, above it for RhB) were observed, and related to the structure of the two dyes and their interactions with water and lipids. Subtle distinctions were found among different ionization forms of the probes. Experimentally, RhB displayed a lipid/water partition coefficient more than two orders of magnitude higher than Rh123, in agreement with free energy profiles obtained from umbrella sampling MD. Conclusions: this work provided detailed insights on the similarities and differences in the behavior of bilayer-inserted Rh123 and RhB, related to the structure of the probes. The much higher affinity of RhB for the membranes increases the local concentration and explains its higher apparent affinity for P-glycoprotein reconstituted in model membranes.

Substrate Access Mechanism in a Novel Membrane-Bound Phospholipase A of Pseudomonas aeruginosa Concordant with Specificity and Regioselectivity

PlaF is a cytoplasmic membrane-bound phospholipase A₁ from Pseudomonas aeruginosa that alters the membrane glycerophospholipid (GPL) composition and fosters the virulence of this human pathogen. PlaF activity is regulated by a dimer-to-monomer transition followed by tilting of the monomer in the membrane. However, how substrates reach the active site and how the characteristics of the active site tunnels determine the activity, specificity, and regioselectivity of PlaF for natural GPL substrates have remained elusive. Here, we combined unbiased and biased all-atom molecular dynamics (MD) simulations and configurational free-energy computations to identify access pathways of GPL substrates to the catalytic center of PlaF. Our results map out a distinct tunnel through which substrates access the catalytic center. PlaF variants with bulky tryptophan residues in this tunnel revealed decreased catalysis rates due to tunnel blockage. The MD simulations suggest that GPLs preferably enter the active site with the sn-1 acyl chain first, which agrees with the experimentally demonstrated PLA₁ activity of PlaF. We propose that the acyl chain-length specificity of PlaF is determined by the structural features of the access tunnel, which results in favorable free energy of binding of medium-chain GPLs. The suggested egress route conveys fatty acid (FA) products to the dimerization interface and, thus, contributes to understanding the product feedback regulation of PlaF by FA-triggered dimerization. These findings open up opportunities for developing potential PlaF inhibitors, which may act as antibiotics against P. aeruginosa.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Transmembrane penetration mechanism of cyclic pollutants inspected by molecular dynamics and metadynamics: the case of morpholine, phenol, 1,4-dioxane and oxane

The presence of industrially produced chemicals in water is often not monitored, while their passive transport and accumulation can cause serious damage in living cells. Molecular dynamics simulations are an effective way to understand the mechanism of the action of these pollutants. In this paper, the passive membrane transport of 1,4-dioxane, phenol, oxane and morpholine was investigated and analyzed thoroughly from structural and energetic points of view. Free energy profiles for pollutant and water penetration into the bilayer were obtained from well-tempered metadynamics (WT-MD) simulations and a mass density-based approach. It was found that all four investigated compounds can penetrate biological membranes and affect the free energy profile of water penetration. Out of the investigated species, oxane has the thermodynamically most preferred position in the bilayer center, leading to a lower free energy barrier of water molecules by 3 kJ mol-1, resulting in 5 times more water molecules in the bilayer center. The concentration dependence of free energy was tested at two different phenol concentrations using WT-MD, and it was found that the higher phenol concentration lowers the main barrier by 3 kJ mol-1. Density-based free energy calculations were found to reproduce the results of WT-MD within the limits of chemical accuracy.

Qualitative similarities and differences in visual object representations between brains and deep networks

Deep neural networks have revolutionized computer vision, and their object representations across layers match coarsely with visual cortical areas in the brain. However, whether these representations exhibit qualitative patterns seen in human perception or brain representations remains unresolved. Here, we recast well-known perceptual and neural phenomena in terms of distance comparisons, and ask whether they are present in feedforward deep neural networks trained for object recognition. Some phenomena were present in randomly initialized networks, such as the global advantage effect, sparseness, and relative size. Many others were present after object recognition training, such as the Thatcher effect, mirror confusion, Weber’s law, relative size, multiple object normalization and correlated sparseness. Yet other phenomena were absent in trained networks, such as 3D shape processing, surface invariance, occlusion, natural parts and the global advantage. These findings indicate sufficient conditions for the emergence of these phenomena in brains and deep networks, and offer clues to the properties that could be incorporated to improve deep networks.

Deep neural networks are widely considered as good models for biological vision. Here, we describe several qualitative similarities and differences in object representations between brains and deep networks that elucidate when deep networks can be considered good models for biological vision and how they can be improved.

Selecting Collective Variables and Free-Energy Methods for Peptide Translocation across Membranes

The selective permeability of cellular membranes is a crucial property for controlled transport into and out of cells. Molecules that can bypass the cellular machinery and spontaneously translocate across membranes could be used as therapeutics or drug carriers. Peptides are a prominent class of such molecules, which include natural and man-developed antimicrobial and cell-penetrating peptides. However, the necessary peptide properties for translocation remain elusive. Computer simulations could uncover these properties once we have a good collective variable (CV) that accurately describes the translocation process. Here, we developed a new CV, which includes a description of peptide insertion, local membrane deformation, and peptide internal degrees of freedom related to its charged groups. By comparison of CVs, we demonstrated that all these components are necessary for an accurate description of peptide translocation. Moreover, the advantages and disadvantages of three common methods for free-energy calculations with our CV were evaluated using the MARTINI coarse-grained model: umbrella sampling, umbrella sampling with replica exchange, and metadynamics. The developed CV leads to the reliable and effective calculation of the free energy of peptide translocation, and thus, it could be useful in the design of spontaneously translocating peptides.

Binding of warfarin differently affects the thermal behavior and chain packing of anionic, zwitterionic and cationic lipid membranes

Warfarin is a coumarin derivative drug widely used for its anticoagulant properties. The interaction of warfarin with fully hydrated lipid bilayers has been studied by combining differential scanning calorimetry, spectrophotometry, electron spin resonance of chain-labelled lipids and molecular docking. Bilayers formed by lipids with different chemico-physical properties were considered, namely dimyristoyl-phosphatidylcholine (DMPC), dimyristoyl-phosphatidylglycerol (DMPG), and dioleoyltrimethyl-ammoniumpropane (DOTAP). We observed in all cases the binding of warfarin in proximity of the surface of the bilayers, leading to a variety of distinct effects on key molecular properties of the membranes. The drug associates with the lipid bilayers in the deprotonated open chain form, with an association constant similar for DMPC and DMPG (1.27·10⁴ and 2.82·10⁴ M^-1, respectively) and lower for DOTAP (0.46·10⁴ M^-1). In DMPC bilayers, which are zwitterionic and with saturated symmetrical chains, warfarin at 10 mol% suppresses the pre-transition, slightly stabilizes the fluid state and reduces the cooperativity of the main transition. Moreover, it alters the lateral packing density of the chain segments close to the polar/apolar interface at any temperature through the gel phase. In anionic DMPG bilayers, the drug slightly perturbs the thermotropic phase behavior, and at 10 mol% markedly loosens the compact gel phase packing of the first chain segments. In cationic DOTAP bilayers, possessing unsaturated acyl chains, the drug induces a slightly higher degree of order and motional restriction in the outer hydrocarbon region in the frozen state. In all cases, as a surface adsorbed molecule, warfarin does not affect the segmental chain order and dynamics for temperatures in the fluid phase. The overall results provide an outline of the action of warfarin on membranes formed by lipids of different types.

The influence of lipid membranes on fluorescent probes' optical properties

BACKGROUND

Organic fluorophores embedded in lipid bilayers can nowadays be described by a multiscale computational approach. Combining different length and time scales, a full characterization of the probe localization and optical properties led to novel insight into the effect of the environments.

SCOPE OF REVIEW

Following an introduction on computational advancements, three relevant probes are reviewed that delineate how a multiscale approach can lead to novel insight into the probes' (non) linear optical properties. Attention is paid to the quality of the theoretical description of the optical techniques.

MAJOR CONCLUSIONS

Computation can assess a priori novel probes' optical properties and guide the analysis and interpretation of experimental data in novel studies. The properties can be used to gain information on the phase and condition of the surrounding biological environment.

GENERAL SIGNIFICANCE

Computation showed that a canonical view on some of the probes should be revisited and adapted.

The cut-off effect of n-alcohols in lipid rafts: A lipid-dependent phenomenon☆

n-Aliphatic alcohols act as anesthetics only up to a certain chain length, beyond which its biological activity disappears. This is known as the 'cut-off' phenomenon. Although the most accepted explanation is based on action sites in membrane proteins, it is not well understood why alcohols alter their functions. The structural dependence of these protein receptors to lipid domains known as 'lipid rafts', suggests a new approach to tackle the puzzling phenomenon. In this work, by performing molecular dynamic simulations (MDS) to explore the lipid role, we provide relevant molecular details about the membrane-alcohol interaction at the cut-off point regime. Since the high variability of the cut-off points found on protein receptors in neurons may be a consequence of differences in the lipid composition surrounding such proteins, our results could have a clear-cut importance.

The Secret Lives of Fluorescent Membrane Probes as Revealed by Molecular Dynamics Simulations

Fluorescent probes have been employed for more than half a century to study the structure and dynamics of model and biological membranes, using spectroscopic and/or microscopic experimental approaches. While their utilization has led to tremendous progress in our knowledge of membrane biophysics and physiology, in some respects the behavior of bilayer-inserted membrane probes has long remained inscrutable. The location, orientation and interaction of fluorophores with lipid and/or water molecules are often not well known, and they are crucial for understanding what the probe is actually reporting. Moreover, because the probe is an extraneous inclusion, it may perturb the properties of the host membrane system, altering the very properties it is supposed to measure. For these reasons, the need for independent methodologies to assess the behavior of bilayer-inserted fluorescence probes has been recognized for a long time. Because of recent improvements in computational tools, molecular dynamics (MD) simulations have become a popular means of obtaining this important information. The present review addresses MD studies of all major classes of fluorescent membrane probes, focusing in the period between 2011 and 2020, during which such work has undergone a dramatic surge in both the number of studies and the variety of probes and properties accessed.

Effect of Protein Flexibility from Coarse-Grained Elastic Network Parameterizations on the Calculation of Free Energy Profiles of Ligand Binding

The characterization of the affinity and binding mechanism of specific molecules to a protein active site is scientifically and industrially relevant for many applications. In principle, this information can be obtained using molecular dynamics (MD) simulations by calculating the free energy profile of the process. However, this is a computationally demanding calculation. Currently, coarse-grained (CG) force fields are very well implemented for MD simulations of biomolecular systems. These computationally efficient force fields are a major advantage to the study of large model systems and/or those requiring long simulation times. The Martini model is currently one of the most popular CG force fields for these systems. For the specific case of protein simulations, to correctly maintain the macromolecular three-dimensional structure, the Martini model needs to include an elastic network (EN). In this work, the effect of protein flexibility, as induced by three EN models compatible with the Martini force field, was tested on the calculation of free energy profiles for protein-ligand binding. The EN models used were ElNeDyn, GoMartini, and GEN. The binding of triolein (TOG) and triacetin (TAG) to a lipase protein (thermomyces lanuginosa lipase-TLL) was used as a case study. The results show that inclusion of greater flexibility in the CG parameterization of proteins is of high importance in the calculation of the free energy profiles of protein-ligand systems. However, care must be taken in order to avoid unjustified large protein deformations. In addition, due to molecular flexibility there may be no absolute need for the center of the ligand to reach the center of the protein-binding site. The calculation of the energy profile to a distance of about 0.5 nm from the active site center can be sufficient to differentiate the affinity of different ligands to a protein.

Atomistic Molecular Dynamics Simulations of Propofol and Fentanyl in Phosphatidylcholine Lipid Bilayers

Atomistic molecular dynamics (MD) and steered MD simulations in combination with umbrella sampling methodology were utilized to study the general anesthetic propofol and the opioid analgesic fentanyl and their interaction with lipid bilayers, which is not yet fully understood. These molecules were inserted into two different fully hydrated phospholipid bilayers, namely, dioleoylphosphatidylcholine (DOPC) and dipalmitoylphosphatidylcholine (DPPC), to investigate the effects that these drugs have on the bilayer. We determined the role of the lipid chain length and saturation on the behavior of the two drugs. Pure, fully hydrated DOPC and DPPC bilayers were also simulated, and the results were in excellent agreement with the experimental values. Various structural and mechanical properties of each system, such as the area per lipid, area compressibility modulus, order parameter, lateral lipid diffusion, hydrogen bonds, and radial distribution functions, have been calculated to assess how the drug molecules affect the different bilayers. From the calculated results, we show that fentanyl and propofol generally follow similar trends in each bilayer but adopt different favorable positions close to the headgroup/chain interface at the carbonyl groups. Propofol was shown to selectively form hydrogen bonds at the carbonyl carbon in each bilayer, whereas fentanyl interacts with water molecules at the headgroup interface. From the calculated free-energy profiles, we determined that both molecules show a preference for the low-density, low-order acyl chain region of the bilayers and both significantly preferred the DOPC bilayer with propofol and fentanyl having energy minima at -6.66 and -43.07 kcal mol^-1, respectively. This study suggests that different chain lengths and levels of saturation directly affect the properties of these two important molecules, which are seen to work together to control anesthesia in surgical applications.

Molecular simulation as a computational pharmaceutics tool to predict drug solubility, solubilization processes and partitioning

In this review we will discuss how computational methods, and in particular classical molecular dynamics simulations, can be used to calculate solubility of pharmaceutically relevant molecules and systems. To the extent possible, we focus on the non-technical details of these calculations, and try to show also the added value of a more thorough and detailed understanding of the solubilization process obtained by using computational simulations. Although the main focus is on classical molecular dynamics simulations, we also provide the reader with some insights into other computational techniques, such as the COSMO-method, and also discuss Flory-Huggins theory and solubility parameters. We hope that this review will serve as a valuable starting point for any pharmaceutical researcher, who has not yet fully explored the possibilities offered by computational approaches to solubility calculations.

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.

Atomistic Picture of Fluorescent Probes with Hydrocarbon Tails in Lipid Bilayer Membranes: An Investigation of Selective Affinities and Fluorescent Anisotropies in Different Environmental Phases

By reverting to spectroscopy, changes in the biological environment of a fluorescent probe can be monitored and the presence of various phases of the surrounding lipid bilayer membranes can be detected. However, it is currently not always clear in which phase the probe resides. The well-known orange 1,1'-dioctadecyl-3,3,3',3'-tetramethylindodicarbo-cyanine perchlorate (DiI-C18(5)) fluorophore, for instance, and the new, blue BODIPY (4,4-difluoro-4-bora-3 a,4 a-diaza- s-indacene) derivative were experimentally seen to target and highlight identical parts of giant unilamellar vesicles of various compositions, comprising mixtures of dipalmitoylphosphatidylcholine (DPPC), dioleoylphosphatidylcholine (DOPC), sphingomyelin (SM), and cholesterol (Chol). However, it was not clear which of the coexisting membrane phases were visualized (Bacalum et al., Langmuir. 2016, 32, 3495). The present study addresses this issue by utilizing large-scale molecular dynamics simulations and the z-constraint method, which allows evaluating Gibbs free-energy profiles. The current calculations give an indication why, at room temperature, both BODIPY and DiI-C18(5) probes prefer the gel (S_o) phase in DOPC/DPPC (2:3 molar ratio) and the liquid-ordered (L_o) phase in DOPC/SM/Chol (1:2:1 molar ratio) mixtures. This study highlights the important differences in orientation and location and therefore in efficiency between the probes when they are used in fluorescence microscopy to screen various lipid bilayer membrane phases. Dependent on the lipid composition, the angle between the transition-state dipole moments of both probes and the normal to the membrane is found to deviate clearly from 90°. It is seen that the DiI-C18(5) probe is located in the headgroup region of the SM/Chol mixture, in close contact with water molecules. A fluorescence anisotropy study also indicates that DiI-C18(5) gives rise to a distinctive behavior in the SM/Chol membrane compared to the other considered membranes. The latter behavior has not been seen for the studied BODIPY probe, which is located deeper in the membrane.

Metronidazole within phosphatidylcholine lipid membranes: New insights to improve the design of imidazole derivatives

Metronidazole is a imidazole derivative with antibacterial and antiprotozoal activity. Despite its therapeutic efficacy, several studies have been developing new imidazole derivatives with lower toxicity. Considering that drug-membrane interactions are key factors for drugs pharmacokinetic and pharmacodynamic properties, the aim of this work is to provide new insights into the structure-toxicity relationship of metronidazole within phosphatidylcholine membranes. For that purpose, lipid membrane models (liposomes and monolayers) composed of dipalmitoylphosphatidylcholine were used. Experimental techniques (determination of partition coefficients and Langmuir isotherm measurements) were combined with molecular dynamics simulations. Different pHs and lipid phases were evaluated to enable a better extrapolation for in vivo conditions. The partition of metronidazole depends on the pH and on the biphasic system (octanol/water or DPPC/water system). At pH 1.2, metronidazole is hydrophilic. At pH 7.4, metronidazole disturbs the order and the packing of phospholipids. For this toxic effect, the hydroxyl group of the side chain of metronidazole is crucial by interacting with the water embedded in the membrane and with the phosphate group and the apolar chains of phospholipids.

Diffusional interaction behavior of NSAIDs in lipid bilayer membrane using molecular dynamics (MD) simulation: Aspirin and Ibuprofen

In this research, for the first time, molecular dynamics (MD) method was used to simulate aspirin and ibuprofen at various concentrations and in neutral and charged states. Effects of the concentration (dosage), charge state, and existence of an integral protein in the membrane on the diffusion rate of drug molecules into lipid bilayer membrane were investigated on 11 systems, for which the parameters indicating diffusion rate and those affecting the rate were evaluated. Considering the diffusion rate, a suitable score was assigned to each system, based on which, analysis of variance (ANOVA) was performed. By calculating the effect size of the indicative parameters and total scores, an optimum system with the highest diffusion rate was determined. Consequently, diffusion rate controlling parameters were obtained: the drug-water hydrogen bond in protein-free systems and protein-drug hydrogen bond in the systems containing protein.

Entropy and Polarity Control the Partition and Transportation of Drug-like Molecules in Biological Membrane

Partition and transportation of drug in the plasma membrane of a mammalian cell are the prerequisite for its function on target protein. Therefore, comprehensive understanding of the physicochemical properties and mechanism behind these complex phenomena is crucial in pharmaceutical research. By using the state-of-art molecular simulations with polarization effect implicitly or explicitly included, we studied the permeation behavior of 2-aminoethoxydiphenyl borate (2-APB), a broad-spectrum modulator for a number of membrane proteins. We showed that the protonation state and therefore the polarity of the drug is critical for its partition, and that the drug is likely to switch between different protonation states along its permeation pathway. By changing the degrees of freedom, protonation further affects the thermodynamic of the permeation pathway of 2-APB, leading to different entropic contributions. A survey on 54 analog structures with similar backbone to 2-APB showed that delicate balance between entropy and polarity plays an important role in drugs’ potency.

Dynamics and energetics of the mammalian phosphatidylinositol transfer protein phospholipid exchange cycle

Phosphatidylinositol-transfer proteins (PITPs) regulate phosphoinositide signaling in eukaryotic cells. The defining feature of PITPs is their ability to exchange phosphatidylinositol (PtdIns) molecules between membranes, and this property is central to PITP-mediated regulation of lipid signaling. However, the details of the PITP-mediated lipid exchange cycle remain entirely obscure. Here, all-atom molecular dynamics simulations of the mammalian StART-like PtdIns/phosphatidylcholine (PtdCho) transfer protein PITPα, both on membrane bilayers and in solvated systems, informed downstream biochemical analyses that tested key aspects of the hypotheses generated by the molecular dynamics simulations. These studies provided five key insights into the PITPα lipid exchange cycle: (i) interaction of PITPα with the membrane is spontaneous and mediated by four specific protein substructures; (ii) the ability of PITPα to initiate closure around the PtdCho ligand is accompanied by loss of flexibility of two helix/loop regions, as well as of the C-terminal helix; (iii) the energy barrier of phospholipid extraction from the membrane is lowered by a network of hydrogen bonds between the lipid molecule and PITPα; (iv) the trajectory of PtdIns or PtdCho into and through the lipid-binding pocket is chaperoned by sets of PITPα residues conserved throughout the StART-like PITP family; and (v) conformational transitions in the C-terminal helix have specific functional involvements in PtdIns transfer activity. Taken together, these findings provide the first mechanistic description of key aspects of the PITPα PtdIns/PtdCho exchange cycle and offer a rationale for the high conservation of particular sets of residues across evolutionarily distant members of the metazoan StART-like PITP family.

Revisiting Partition in Hydrated Bilayer Systems

We discuss potential-of-mean force convergence issues involving the methods for generating starting structures of umbrella sampling simulations for the study of the water-bilayer permeation of drug-like compounds. We provide a proof of concept of the benefits of potential-energy modification methods, in particular the flooding technique, to enhance the sampling and improve the convergence of the potential-of-mean force. In addition, we also discuss how these potential modifiers could introduce the necessary bias to induce the correct description of other relevant internal degrees of freedom for the compounds, thus providing a more accurate description of the permeation process. Effective methodological approximations are taken into consideration, more specifically the inclusion of atomic polarization, and experimental data are used in the validation of computed values whenever possible. The parametrization of the atomic point charges is always a relevant issue in biomolecular simulations. We discuss the influence of a molecule's conformation-dependent atomic point charges and their impact on the predicted potential-of-mean force.

Partitioning into Colloidal Structures of Fasted State Intestinal Fluid Studied by Molecular Dynamics Simulations

We performed molecular dynamics (MD) simulations to obtain insights into the structure and molecular interactions of colloidal structures present in fasted state intestinal fluid. Drug partitioning and interaction were studied with a mixed system of the bile salt taurocholate (TCH) and 1,2-dilinoleoyl-sn-glycero-3-phosphocholine (DLiPC). Spontaneous aggregation of TCH and DLiPC from unconstrained MD simulations at the united-atom level using the Berger/Gromos54A7 force fields demonstrated that intermolecular hydrogen bonding between TCH molecules was an important factor in determining the overall TCH and DLiPC configuration. In bilayered systems, these intermolecular hydrogen bonds resulted in embedded transmembrane TCH clusters. Free energy simulations using the umbrella sampling technique revealed that the stability of these transmembrane TCH clusters was superior when they consisted of 3 or 4 TCH per bilayer leaflet. All-atom simulations using the Slipids/GAFF force fields showed that the TCH embedded in the bilayer decreased the energy barrier to penetrate the bilayer (ΔG_pen) for water, ethanol, and carbamazepine, but not for the more lipophilic felodipine and danazol. This suggests that diffusion of hydrophilic to moderately lipophilic molecules through the bilayer is facilitated by the embedded TCH molecules. However, the effect of embedded TCH on the overall lipid/water partitioning was significant for danazol, indicating that the incorporation of TCH plays a crucial role for the partitioning of lipophilic solutes into e.g. lipidic vesicles existing in fasted state intestinal fluids. To conclude, the MD simulations revealed important intermolecular interactions in lipidic bilayers, both between the bile components themselves and with the drug molecules.

Two Relations to Estimate Membrane Permeability Using Milestoning

Prediction of passive permeation rates of solutes across lipid bilayers is important to drug design, toxicology, and other biological processes such as signaling. The inhomogeneous solubility-diffusion (ISD) equation is traditionally used to relate the position-dependent potential of mean force and diffusivity to the permeability coefficient. The ISD equation is derived via the Smoluchowski equation and assumes overdamped system dynamics. It has been suggested that the complex membrane environment may exhibit more complicated damping conditions. Here we derive a variant of the inhomogeneous solubility diffusion equation as a function of the mean first passage time (MFPT) and show how milestoning, a method that can estimate kinetic quantities of interest, can be used to estimate the MFPT of membrane crossing and, by extension, the permeability coefficient. We further describe a second scheme, agnostic to the damping condition, to estimate the permeability coefficient from milestoning results or other methods that compute a probability of membrane crossing. The derived relationships are tested using a one-dimensional Langevin dynamics toy system confirming that the presented theoretical methods can be used to estimate permeabilities given simulation and milestoning results.

Characterization of molecular association of poly(2-oxazoline)s-based micelles with various epoxides and diols via the Flory-Huggins theory: a molecular dynamics simulation approach

The hydrolytic kinetic resolution (HKR) of epoxides has been performed in a shell-crosslinked micellar (SCM) nanoreactor consisting of amphiphilic triblock copolymers based on poly(2-oxazline)s polymer derivatives with attached Co(iii)-salens to the micelle core. To investigate the effect of the molecular interaction of reactant/product molecules with the SCM nanoreactor on the rate of HKR, we calculated the Flory-Huggins interaction parameters (χ) using the molecular dynamics simulation method. For this, the blend systems were constructed with various compositions such as 15, 45, and 70 wt% of the reactant/product molecules with respect to the polymers such as poly(2-methyl-2-oxazoline) (PMOX), poly(2-(3-butinyl)2-oxazoline) (PBOX), and poly(methyl-3-oxazol-2-yl)pentanoate with Co(iii)-salen (PSCoX). From the χ parameters, we demonstrate that the miscibility of reactants/products with polymers has a strong correlation with the experimental reaction rate of the HKR: phenyl glycidyl ether (Reac-OPh) > epoxyhexane (Reac-C4) > styrene oxide (Reac-Ph) > epichlorohydrin (Reac-Cl). To validate this finding, we also conducted the potential of mean force analysis using steered molecular dynamics simulation for the molecular displacement of Reac-Cl and Reac-OPh through PMOX and PSCoX, revealing that the free energy reduction was greater when Reac-OPh molecule enters the polymer phase compared to Reac-Cl, which agrees with the findings from the χ parameters calculations.

Structural perturbation of a dipalmitoylphosphatidylcholine (DPPC) bilayer by warfarin and its bolaamphiphilic analogue: A molecular dynamics study

Compounds with nominally similar biological activity may exhibit differential toxicity due to differences in their interactions with cell membranes. Many pharmaceutical compounds are amphiphilic and can be taken up by phospholipid bilayers, interacting strongly with the lipid-aqueous interface whether or not subsequent permeation through the bilayer is possible. Bolaamphiphilic compounds, which possess two hydrophilic ends and a hydrophobic linker, can likewise undergo spontaneous uptake by bilayers. While membrane-spanning bolaamphiphiles can stabilize membranes, small molecules with this characteristic have the potential to create membrane defects via disruption of bilayer structure and dynamics. When compared to the amphiphilic therapeutic anticoagulant, warfarin, the bolaamphiphilic analogue, brodifacoum, exhibits heightened toxicity that goes beyond superior inhibition of the pharmacological target enzyme. We explore, herein, the consequences of anticoagulant accumulation in a dipalmitoylphosphatidylcholine (DPPC) bilayer. Coarse-grained molecular dynamics simulations reveal that permeation of phospholipid bilayers by brodifacoum causes a disruption of membrane barrier function that is driven by the bolaamphiphilic nature and size of this molecule. We find that brodifacoum partitioning into bilayers causes membrane thinning and permeabilization and promotes lipid flip-flop - phenomena that are suspected to play a role in triggering cell death. These phenomena are either absent or less pronounced in the case of the less toxic, amphiphilic compound, warfarin.

Simulation-Based Approaches for Determining Membrane Permeability of Small Compounds

Predicting the rate of nonfacilitated permeation of solutes across lipid bilayers is important to drug design, toxicology, and signaling. These rates can be estimated using molecular dynamics simulations combined with the inhomogeneous solubility-diffusion model, which requires calculation of the potential of mean force and position-dependent diffusivity of the solute along the transmembrane axis. In this paper, we assess the efficiency and accuracy of several methods for the calculation of the permeability of a model DMPC bilayer to urea, benzoic acid, and codeine. We compare umbrella sampling, replica exchange umbrella sampling, adaptive biasing force, and multiple-walker adaptive biasing force for the calculation of the transmembrane PMF. No definitive advantage for any of these methods in their ability to predict the membrane permeability coefficient Pm was found, provided that a sufficiently long equilibration is performed. For diffusivities, a Bayesian inference method was compared to a generalized Langevin method, both being sensitive to chosen parameters and the slow relaxation of membrane defects. Agreement within 1.5 log units of the computed Pm with experiment is found for all permeants and methods. Remaining discrepancies can likely be attributed to limitations of the force field as well as slowly relaxing collective movements within the lipid environment. Numerical calculations based on model profiles show that Pm can be reliably estimated from only a few data points, leading to recommendations for calculating Pm from simulations.

Permeability across lipid membranes

Molecular permeation through lipid membranes is a fundamental biological process that is important for small neutral molecules and drug molecules. Precise characterization of free energy surface and diffusion coefficients along the permeation pathway is required in order to predict molecular permeability and elucidate the molecular mechanisms of permeation. Several recent technical developments, including improved molecular models and efficient sampling schemes, are illustrated in this review. For larger penetrants, explicit consideration of multiple collective variables, including orientational, conformational degrees of freedom, are required to be considered in addition to the distance from the membrane center along the membrane normal. Although computationally demanding, this method can provide significant insights into the molecular mechanisms of permeation for molecules of medical and pharmaceutical importance. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

A Blue-Light-Emitting BODIPY Probe for Lipid Membranes

Here we describe a new BODIPY-based membrane probe (1) that provides an alternative to dialkylcarbocyanine dyes, such as DiI-C18, that can be excited in the blue spectral region. Compound 1 has unbranched octadecyl chains at the 3,5-positions and a meso-amino function. In organic solvents, the absorption and emission maxima of 1 are determined mainly by solvent acidity and dipolarity. The fluorescence quantum yield is high and reaches 0.93 in 2-propanol. The fluorescence decays are well fitted with a single-exponential in pure solvents and in small and giant unilamellar vesicles (GUV) with a lifetime of ca. 4 ns. Probe 1 partitions in the same lipid phase as DiI-C18(5) for lipid mixtures containing sphingomyelin and for binary mixtures of dipalmitoylphosphatidylcholine (DPPC) and dioleoylphosphatidylcholine (DOPC). The lipid phase has no effect on the fluorescence lifetime but influences the fluorescence anisotropy. The translational diffusion coefficients of 1 in GUVs and OLN-93 cells are of the same order as those reported for DiI-C18. The directions of the absorption and emission transition dipole moments of 1 are calculated to be parallel. This is reflected in the high steady-state fluorescence anisotropy of 1 in high ordered lipid phases. Molecular dynamic simulations of 1 in a model of the DOPC bilayer indicate that the average angle of the transition moments with respect to membrane normal is ca. 70°, which is comparable with the value reported for DiI-C18.

Sampling errors in free energy simulations of small molecules in lipid bilayers

Free energy simulations are a powerful tool for evaluating the interactions of molecular solutes with lipid bilayers as mimetics of cellular membranes. However, these simulations are frequently hindered by systematic sampling errors. This review highlights recent progress in computing free energy profiles for inserting molecular solutes into lipid bilayers. Particular emphasis is placed on a systematic analysis of the free energy profiles, identifying the sources of sampling errors that reduce computational efficiency, and highlighting methodological advances that may alleviate sampling deficiencies. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Synergism of antioxidant action of vitamins E, C and quercetin is related to formation of molecular associations in biomembranes

Vitamins E, C and polyphenols (flavonoids and non-flavonoids) are major natural antioxidants capable of preventing damage generated by oxidative stress. Here we show the capacity of these antioxidants to form non-covalent association within lipid bilayers close to the membrane/cytosol interface. Antioxidant regeneration is significantly enhanced in these complexes.

Liquid–liquid extraction of alkali cations by 18-crown-6: complexation and interface crossing studied by MD and PMF simulations

The 18C6/M⁺Pic⁻complexes form and adsorb “right at the nano-interface” where 18C6 prefers the K⁺guest.

From dioxin to dioxin congeners: understanding the differences in hydrophobic aggregation in water and absorption into lipid membranes by means of atomistic simulations

Molecular dynamics simulations indicate that dioxin likely accumulates in cell membranes more than its congeners (anthracene, tetrahydrodioxin, 3,3′,5,5′-tetrachlorobiphenyl, and 1,2-dihydroxytetrahydrodibenzo-p-dioxin).

Molecular simulation of nonfacilitated membrane permeation

This is a review. Non-electrolytic compounds typically cross cell membranes by passive diffusion. The rate of permeation is dependent on the chemical properties of the solute and the composition of the lipid bilayer membrane. Predicting the permeability coefficient of a solute is important in pharmaceutical chemistry and toxicology. Molecular simulation has proven to be a valuable tool for modeling permeation of solutes through a lipid bilayer. In particular, the solubility-diffusion model has allowed for the quantitative calculation of permeability coefficients. The underlying theory and computational methods used to calculate membrane permeability are reviewed. We also discuss applications of these methods to examine the permeability of solutes and the effect of membrane composition on permeability. The application of coarse grain and polarizable models is discussed. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Light-Driven Ca(2+) Ion Pump: How Does It Work?

Work done by Bennett et al. [ Nature 2002 , 420 , 398 - 401 ] demonstrated that Ca(2+) ions can be actively transported through a lipid bilayer membrane by an artificial photosynthetic machine. However, details of the pump process, such as the oxidation state of the shuttle molecule and stoichiometry of the shuttle-ion complex, are not fully understood, which hinders the development of ion pumps of this type with higher efficiency. In this study, we combine all atom molecular dynamics simulations and quantum mechanics calculations to estimate the time scale of the shuttle-ion complex diffusion process and charge transfer step. We find that the process of shuttle-ion complex diffusion across the lipid bilayer membrane is the rate-limiting step, with a time scale of seconds to minutes. Other processes such as charge transfer between the redox reaction center and the shuttle molecule have picoseconds time scales. We also show that a shuttle-ion complex with 2:1 stoichiometry ratio has a lower energy barrier across the lipid membrane than other choices of complexes. The calculations show that the Ca(2+) ion is likely to be shuttled by a semiquinone type of shuttle molecule as this has the lowest free energy barrier across the lipid bilayer membrane, the fewest electrons transferred in the redox cycle, and it does not generate (or require) proton flow. Estimates of ion flow rates are consistent with measured values.