Shedding light on the puzzle of drug-membrane interactions: Experimental techniques and molecular dynamics simulations

An emerging terpolymeric nanoparticle pore former as an internal recrystallization inhibitor of celecoxib in controlled release amorphous solid dispersion beads: Experimental studies and molecular dynamics analysis

Solid oral controlled release formulations feature numerous clinical advantages for drug candidates with adequate solubility and dissolution rate. However, most new chemical entities exhibit poor water solubility, and hence are exempt from such benefits. Although combining drug amorphization with controlled release formulation is promising to elevate drug solubility, like other supersaturating systems, the problem of drug recrystallization has yet to be resolved, particularly within the dosage form. Here, we explored the potential of an emerging, non-leachable terpolymer nanoparticle (TPN) pore former as an internal recrystallization inhibitor within controlled release amorphous solid dispersion (CRASD) beads comprising a poorly soluble drug (celecoxib) reservoir and insoluble polymer (ethylcellulose) membrane. Compared to conventional pore former, polyvinylpyrrolidone (PVP), TPN-containing membranes exhibited superior structural integrity, less crystal formation at the CRASD bead surface, and greater extent of celecoxib release. All-atom molecular dynamics analyses revealed that in the presence of TPN, intra-molecular bonding, crystal formation tendency, diffusion coefficient, and molecular flexibility of celecoxib were reduced, while intermolecular H-bonding was increased as compared to PVP. This work suggests that selection of a pore former that promotes prolonged molecular separation within a nanoporous controlled release membrane structure may serve as an effective strategy to enhance amorphicity preservation inside CRASD.

Aspirin Interacts with Cholesterol-Containing Membranes in a pH-Dependent Manner

Aspirin has been used for broad therapeutic treatment, including secondary prevention of cardiovascular disease associated with increased cholesterol levels. Aspirin and other nonsteroidal anti-inflammatory drugs have been shown to interact with lipid membranes and change their biophysical properties. In this study, mixed lipid model bilayers made from 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) or 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) comprising varying concentrations of cholesterol (10:1, 4:1, and 1:1 mole ratio of lipid:chol), prepared by the droplet interface bilayer method, were used to examine the effects of aspirin at various pH on transbilayer water permeability. The presence of aspirin increases the water permeability of POPC bilayers in a concentration-dependent manner, with a greater magnitude of increase at pH 3 compared to pH 7. In the presence of cholesterol, aspirin is similarly shown to increase water permeability; however, the extent of the increase depends on both the concentration of cholesterol and the pH, with the least pronounced enhancement in water permeability at high cholesterol levels at pH 7. A fusion of data from differential scanning calorimetry, confocal Raman microspectrophotometry, and interfacial tensiometric measurements demonstrates that aspirin can promote significant thermal, structural, and interfacial property perturbations in the mixed-lipid POPC or DOPC membranes containing cholesterol, indicating a disordering effect on the lipid membranes. Our findings suggest that aspirin fluidizes phosphocholine membranes in both cholesterol-free and cholesterol-enriched states and that the overall effect is greater when aspirin is in a neutral state. These results confer a deeper comprehension of the divergent effects of aspirin on biological membranes having heterogeneous compositions, under varying physiological pH and different cholesterol compositions, with implications for a better understanding of the gastrointestinal toxicity induced by the long term use of this important nonsteroidal anti-inflammatory molecule.

Concentration-dependent mechanism of the binding behavior of ibuprofen to the cell membrane: A molecular dynamic simulation study

Ibuprofen is a commonly used drug for treating headaches, pain, and fever. The lipid bilayer is the primary and most important interface for drugs to interact with biological systems. However, the molecular interactions between ibuprofen and the cell membrane are not well understood. Our findings suggest that the interactions between ibuprofen and the bilayer involve multiple steps and depend on the concentration of the drug. At low concentrations of ibuprofen, it can bind to the surface of the lipid bilayer. The electrostatic and vdW energies of IBU-lipid at 0 ns of the simulation were -22.5 ± 3.2 and -5.9 ± 1.2 kj.mol-1 Fig. 2. In the following, the vdW energy of the IBU-lipid was increased by around -134.6 ± 3.7 kj.mol-1 whereas the electrostatic energy of the IBU-lipid was significantly decreased. This binding is facilitated by electrostatic and vdW interactions between ibuprofen and the head group of lipids. In the second step, ibuprofen is inserted into the lipid bilayer and positioned at the interface between the bilayer and the aqueous phase. In high concentrations of ibuprofen, it moved to the central region of the lipid bilayer. At this concentration, the physical and structural properties of the cell membrane change significantly. Results from the radial distribution function analysis indicate that at low concentrations, ibuprofen molecules are situated close to the head groups of phosphate groups. However, at high concentrations of ibuprofen, these molecules move to the inner side of the lipid bilayer. In addition, our findings indicate that at low concentrations of ibuprofen, these molecules did not significantly alter the physical properties of the cell membrane. In contrast, at high concentrations of ibuprofen, the physical parameters of the hydrocarbon tails, such as thickness, fluidity, and order, changed dramatically. APL parameter for POPC membrane increased slightly to 0.60 and 0.63 nm2 in the presence of low and high concentrations of ibuprofen molecules. The three-step interaction between ibuprofen and the lipid bilayer involves several events, such as the movement of ibuprofen molecules towards the central region of the lipid bilayer and the deformation and alteration of the structural and stability properties of the cell membrane. These effects are observed only at high concentrations of ibuprofen. It appears that the side effects of ibuprofen overdose are related to changes in the properties of the cell membrane and, subsequently, the function of membrane-anchored target proteins.

Theoretical and Experimental Study of Molecular Interactions of Fluralaner with Lipid Membranes

Fluralaner is a relatively new insecticide belonging to the isoxazoline group, whose action mechanism involves the blocking of GABA_A-receptors in the insect nervous system. Because of its high hydrophobicity, fluralaner could bioaccumulate and reach toxic local concentrations. Since there are no data available about the penetration and persistence of isoxazolines in biological membranes, we intend to evaluate fluralaner permanence as a pollutant by using model membranes. We used experimental and in silico models to characterize the incorporation of fluralaner into the lipid phase at different packing states. We determined its impact in the membrane structure and organization. Our results confirm that fluralaner is capable of penetrating, holding, and accumulating in the lipid membrane and provide details on its precise location and orientation. These properties would allow fluralaner to reach high local concentrations in different membranes and organs, which could be dangerous for vertebrate organisms if its handling is not properly controlled.

Current Trends and Changes in Use of Membrane Molecular Dynamics Simulations within Academia and the Pharmaceutical Industry

There has been an almost exponential increase in the use of molecular dynamics simulations in basic research and industry over the last 5 years, with almost a doubling in the number of publications each year. Many of these are focused on neurological membranes, and biological membranes in general, applied to the medical industry. A smaller portion have utilized membrane simulations to answer more basic questions related to the function of specific proteins, chemicals or biological processes. This review covers some newer studies, alongside studies from the last two decades, to determine changes in the field. Some of these are basic, while others are more profound, such as multi-component embedded membrane machinery. It is clear that many facets of the discipline remain the same, while the focus on and uses of the technology are broadening in scope and utilization as a general research tool. Analysis of recent literature provides an overview of the current methodologies, covers some of the recent trends or advances and tries to make predictions of the overall path membrane molecular dynamics will follow in the coming years. In general, the overview presented is geared towards the general scientific community, who may wish to introduce the use of these methodologies in light of these changes, making molecular dynamic simulations more feasible for general scientific or medical research.

Characterisation of the Molecular Mechanism of Permeation of the Prodrug Me-5ALA across the Human Stratum Corneum Using Molecular Dynamics Simulations

The barrier imposed by the outer layer of the skin, the stratum corneum, creates an almost impermeable environment for exogenous substances. Few lipophilic drugs with low molecular mass can passively diffuse through this layer, highlighting the need to develop methods to enable the delivery of more drugs via the transdermal route. The prodrug approach involves modifying the structure of a drug molecule to enhance its permeability across the skin, but it is often difficult to predict how exactly changes in chemical structure affect permeation. This study uses molecular dynamics simulations to predict permeability values and adequately characterise the molecular mechanism of permeation of the prodrugs Me-5ALA and its parent compound 5ALA across a molecular model of the lipid bilayers of the human stratum corneum. The influence of increased hydrophobicity in Me-5ALA on its permeation revealed a reduction in hydrogen bonding capability that enables it to interact more favourably with the hydrophobic region of the bilayer and diffuse at a faster rate with less resistance, thus making it a better permeant compared to its more hydrophilic parent compound. This molecular simulation approach offers a promising route for the rational design of drug molecules that can permeate effectively across the stratum corneum.

Partition of the Reactive Species of the Suzuki-Miyaura Reaction between Aqueous and Micellar Environments

The Suzuki-Miyaura reaction between the aryl halide (1) and the phenyl boronic acid (2), in the presence of the palladium(0) complex (3) as catalyst, gives the cross-coupling product (4) in quantitative yield when performed in basic aqueous solution of the nonionic surfactant Kolliphor-EL (K-EL). The partition between the aqueous and micellar environments of the species of this reaction has been investigated by means of Molecular Dynamics (MD) simulations. Starting from the K-EL molecules dispersed in water, a micelle model has been generated by MD simulations, adopting the 2016H66 force field. Reagent and product species have been described with the same force field, once the reliability of this force field has been tested comparing the n-octanol/water partition free energies calculated from the MD and Free Energy Perturbation (FEP) method with those obtained from the quantum-mechanical SMD method. The potential of mean force for the transfer process between water and the micellar phase of the different species has been calculated by the MD simulations and the Umbrella Sampling (US) method. The overall picture that emerges from these results confirms that the molecular species involved in this reaction prefers the micellar environment and concentrates in different but close zones of the micelle. This supports the experimental evidence that the use of suitable surfactant agents promotes reactivity, allowing micelles to behave as nanoreactors in which reactive species are solubilized and enhance their local concentration.

Drug–Membrane Interaction as Revealed by Spectroscopic Methods: The Role of Drug Structure in the Example of Rifampicin, Levofloxacin and Rapamycin

We have investigated the nature of the interaction of small organic drug molecules with lipid membranes of various compositions. Using infrared spectroscopy and differential scanning calorimetry methods, we studied the role of the structure of the active molecule in interaction with the membrane using the example of dipalmitoylphosphatidylcholine (DPPC) and dipalmitoylphosphatidylcholine:cardiolipin (DPPC:CL) liposomes. We discovered the key role of the heterocycle in interaction with the polar part of the bilayer and the network of unsaturated bonds in interaction with the hydrophobic part. For rifampicin and levofloxacin, the main binding sites were phosphate and carbonyl groups of lipids, and in the case of anionic liposomes we found a slight penetration of rifampicin into the hydrophobic part of the bilayer. For rapamycin, experimental confirmation of the localization of the molecule in the region of fatty acid chains was obtained, and perturbation in the region of phosphate groups was demonstrated for the first time. The process of phase transition of liposomal forms of rifampicin and levofloxacin was studied. DPPC liposomes accelerate the phase transition when loaded with a drug. DPPC:CL liposomes are less susceptible to changes in the phase transition rate.

Effects of the RNA-Polymerase Inhibitors Remdesivir and Favipiravir on the Structure of Lipid Bilayers-An MD Study

The structure and dynamics of membranes are crucial to ensure the proper functioning of cells. There are some compounds used in therapeutics that show nonspecific interactions with membranes in addition to their specific molecular target. Among them, two compounds recently used in therapeutics against COVID-19, remdesivir and favipiravir, were subjected to molecular dynamics simulation assays. In these, we demonstrated that the compounds can spontaneously bind to model lipid membranes in the presence or absence of cholesterol. These findings correlate with the corresponding experimental results recently reported by our group. In conclusion, insertion of the compounds into the membrane is observed, with a mean position close to the phospholipid head groups.

Bioavailability by design — Vitamin D3 liposomal delivery vehicles

Vitamin D3 deficiency has serious health consequences, as demonstrated by its effect on severity and recovery after COVID-19 infection. Because of high hydrophobicity, its absorption and subsequent redistribution throughout the body are inherently dependent on the accompanying lipids and/or proteins. The effective oral vitamin D₃ formulation should ensure penetration of the mucus layer followed by internalization by competent cells. Isothermal titration calorimetry and computer simulations show that vitamin D₃ molecules cannot leave the hydrophobic environment, indicating that their absorption is predominantly driven by the digestion of the delivery vehicle. In the clinical experiment, liposomal vitamin D₃ was compared to the oily formulation. The results obtained show that liposomal vitamin D₃ causes a rapid increase in the plasma concentration of calcidiol. No such effect was observed when the oily formulation was used. The effect was especially pronounced for people with severe vitamin D₃ deficiency.

Review on the QM/MM Methodologies and Their Application to Metalloproteins

The multiscaling quantum mechanics/molecular mechanics (QM/MM) approach was introduced in 1976, while the extensive acceptance of this methodology started in the 1990s. The combination of QM/MM approach with molecular dynamics (MD) simulation, otherwise known as the QM/MM/MD approach, is a powerful and promising tool for the investigation of chemical reactions’ mechanism of complex molecular systems, drug delivery, properties of molecular devices, organic electronics, etc. In the present review, the main methodologies in the multiscaling approaches, i.e., density functional theory (DFT), semiempirical methodologies (SE), MD simulations, MM, and their new advances are discussed in short. Then, a review on calculations and reactions on metalloproteins is presented, where particular attention is given to nitrogenase that catalyzes the conversion of atmospheric nitrogen molecules N₂ into NH₃ through the process known as nitrogen fixation and the FeMo-cofactor.

Lipid-mediated antimicrobial resistance: a phantom menace or a new hope?

The proposition of a post-antimicrobial era is all the more realistic with the continued rise of antimicrobial resistance. The development of new antimicrobials is failing to counter the ever-increasing rates of bacterial antimicrobial resistance. This necessitates novel antimicrobials and drug targets. The bacterial cell membrane is an essential and highly conserved cellular component in bacteria and acts as the primary barrier for entry of antimicrobials into the cell. Although previously under-exploited as an antimicrobial target, the bacterial cell membrane is attractive for the development of novel antimicrobials due to its importance in pathogen viability. Bacterial cell membranes are diverse assemblies of macromolecules built around a central lipid bilayer core. This lipid bilayer governs the overall membrane biophysical properties and function of its membrane-embedded proteins. This mini-review will outline the mechanisms by which the bacterial membrane causes and controls resistance, with a focus on alterations in the membrane lipid composition, chemical modification of constituent lipids, and the efflux of antimicrobials by membrane-embedded efflux systems. Thorough insight into the interplay between membrane-active antimicrobials and lipid-mediated resistance is needed to enable the rational development of new antimicrobials. In particular, the union of computational approaches and experimental techniques for the development of innovative and efficacious membrane-active antimicrobials is explored.

Toxicity pathways of lipid metabolic disorders induced by typical replacement flame retardants via data-driven analysis, in silico and in vitro approaches

Endocrine-disrupting chemicals can interfere with hormone action via various pathways, thereby increasing the risk of adverse health outcomes. Organophosphorus ester (OPEs) retardants, a group of new emerging endocrine disruption chemicals, have been referred to as metabolism disruptors and reported to induce chronic health problems. However, the toxicity pathways were mainly focused on nuclear receptor signaling pathways. Significantly, the membrane receptor pathway (such as G protein-coupled estrogen receptor 1 (GPER) signaling pathway) had been gradually realized as the important role in respond more effective to lipid metabolism disorder than traditional nuclear receptors, whereas the detailed mechanism was unclear yet. Therefore, this study innovatively integrated the bibliometric analysis, in silico and in vitro approach to develop toxicity pathways for the mechanism interpretation. Bibliometric analysis found that the typical OPEs - triphenyl phosphate was a major concern of lipid metabolism abnormality. Results verified that TPP could damage the structures of cell membranes and exert an agonistic effect of GPER as the molecular initiating event. Then, the activated GPER could trigger the PI3K-Akt/NCOR1 and mTOR/S6K2/PPARα transduction pathways as key event 1 (KE1) and affect the process of lipid metabolism and synthesis (CPT1A, CPT2, SREBF2 and SCD) as KE2. As a result, these alterations led to lipid accumulation as adverse effect at cellular-levels. Furthermore, the potential outcomes (such as immunity damage, weight change and steatohepatitis) at high biological levels were expanded. These findings improved knowledge to deeply understand toxicity pathways of phosphorus flame retardants and then provided a theoretical basis for risk assessments.

Phosphatidylserine Exposed Lipid Bilayer Models for Understanding Cancer Cell Selectivity of Natural Compounds: A Molecular Dynamics Simulation Study

Development of drugs that are selectively toxic to cancer cells and safe to normal cells is crucial in cancer treatment. Evaluation of membrane permeability is a key metric for successful drug development. In this study, we have used in silico molecular models of lipid bilayers to explore the effect of phosphatidylserine (PS) exposure in cancer cells on membrane permeation of natural compounds Withaferin A (Wi-A), Withanone (Wi-N), Caffeic Acid Phenethyl Ester (CAPE) and Artepillin C (ARC). Molecular dynamics simulations were performed to compute permeability coefficients. The results indicated that the exposure of PS in cancer cell membranes facilitated the permeation of Wi-A, Wi-N and CAPE through a cancer cell membrane when compared to a normal cell membrane. In the case of ARC, PS exposure did not have a notable influence on its permeability coefficient. The presented data demonstrated the potential of PS exposure-based models for studying cancer cell selectivity of drugs.

Molecular modeling in cardiovascular pharmacology: Current state of the art and perspectives

Molecular modeling in pharmacology is a promising emerging tool for exploring drug interactions with cellular components. Recent advances in molecular simulations, big data analysis, and artificial intelligence (AI) have opened new opportunities for rationalizing drug interactions with their pharmacological targets. Despite the obvious utility and increasing impact of computational approaches, their development is not progressing at the same speed in different fields of pharmacology. Here, we review current in silico techniques used in cardiovascular diseases (CVDs), cardiological drug discovery, and assessment of cardiotoxicity. In silico techniques are paving the way to a new era in cardiovascular medicine, but their use somewhat lags behind that in other fields.

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.

The Plasma Membrane at the Cornerstone Between Flexibility and Adaptability: Implications for Saccharomyces cerevisiae as a Cell Factory

In the last decade, microbial-based biotechnological processes are paving the way toward sustainability as they implemented the use of renewable feedstocks. Nonetheless, the viability and competitiveness of these processes are often limited due to harsh conditions such as: the presence of feedstock-derived inhibitors including weak acids, non-uniform nature of the substrates, osmotic pressure, high temperature, extreme pH. These factors are detrimental for microbial cell factories as a whole, but more specifically the impact on the cell’s membrane is often overlooked. The plasma membrane is a complex system involved in major biological processes, including establishing and maintaining transmembrane gradients, controlling uptake and secretion, intercellular and intracellular communication, cell to cell recognition and cell’s physical protection. Therefore, when designing strategies for the development of versatile, robust and efficient cell factories ready to tackle the harshness of industrial processes while delivering high values of yield, titer and productivity, the plasma membrane has to be considered. Plasma membrane composition comprises diverse macromolecules and it is not constant, as cells adapt it according to the surrounding environment. Remarkably, membrane-specific traits are emerging properties of the system and therefore it is not trivial to predict which membrane composition is advantageous under certain conditions. This review includes an overview of membrane engineering strategies applied to Saccharomyces cerevisiae to enhance its fitness under industrially relevant conditions as well as strategies to increase microbial production of the metabolites of interest.

Improving the accuracy and convergence of drug permeation simulations via machine-learned collective variables

Understanding the permeation of biomolecules through cellular membranes is critical for many biotechnological applications, including targeted drug delivery, pathogen detection, and the development of new antibiotics. To this end, computer simulations are routinely used to probe the underlying mechanisms of membrane permeation. Despite great progress and continued development, permeation simulations of realistic systems (e.g., more complex drug molecules or biologics through heterogeneous membranes) remain extremely challenging if not intractable. In this work, we combine molecular dynamics simulations with transition-tempered metadynamics and techniques from the variational approach to conformational dynamics to study the permeation mechanism of a drug molecule, trimethoprim, through a multicomponent membrane. We show that collective variables (CVs) obtained from an unsupervised machine learning algorithm called time-structure based Independent Component Analysis (tICA) improve performance and substantially accelerate convergence of permeation potential of mean force (PMF) calculations. The addition of cholesterol to the lipid bilayer is shown to increase both the width and height of the free energy barrier due to a condensing effect (lower area per lipid) and increase bilayer thickness. Additionally, the tICA CVs reveal a subtle effect of cholesterol increasing the resistance to permeation in the lipid head group region, which is not observed when canonical CVs are used. We conclude that the use of tICA CVs can enable more efficient PMF calculations with additional insight into the permeation mechanism.

The Whole Is Bigger than the Sum of Its Parts: Drug Transport in the Context of Two Membranes with Active Efflux

Cell envelope plays a dual role in the life of bacteria by simultaneously protecting it from a hostile environment and facilitating access to beneficial molecules. At the heart of this ability lie the restrictive properties of the cellular membrane augmented by efflux transporters, which preclude intracellular penetration of most molecules except with the help of specialized uptake mediators. Recently, kinetic properties of the cell envelope came into focus driven on one hand by the urgent need in new antibiotics and, on the other hand, by experimental and theoretical advances in studies of transmembrane transport. A notable result from these studies is the development of a kinetic formalism that integrates the Michaelis-Menten behavior of individual transporters with transmembrane diffusion and offers a quantitative basis for the analysis of intracellular penetration of bioactive compounds. This review surveys key experimental and computational approaches to the investigation of transport by individual translocators and in whole cells, summarizes key findings from these studies and outlines implications for antibiotic discovery. Special emphasis is placed on Gram-negative bacteria, whose envelope contains two separate membranes. This feature sets these organisms apart from Gram-positive bacteria and eukaryotic cells by providing them with full benefits of the synergy between slow transmembrane diffusion and active efflux.

Recent Experimental Developments in Studying Passive Membrane Transport of Drug Molecules

The ability to measure the passive membrane permeation of drug-like molecules is of fundamental biological and pharmaceutical importance. Of significance, passive diffusion across the cellular membranes plays an effective role in the delivery of many pharmaceutical agents to intracellular targets. Hence, approaches for quantitative measurement of membrane permeability have been the topics of research for decades, resulting in sophisticated biomimetic systems coupled with advanced techniques. In this review, recent developments in experimental approaches along with theoretical models for quantitative and real-time analysis of membrane transport of drug-like molecules through mimetic and living cell membranes are discussed. The focus is on time-resolved fluorescence-based, surface plasmon resonance, and second-harmonic light scattering approaches. The current understanding of how properties of the membrane and permeant affect the permeation process is discussed.

Ibuprofen and the Phosphatidylcholine Bilayer: Membrane Water Permeability in the Presence and Absence of Cholesterol

The interactions between drugs and cell membranes can modulate the structural and physical properties of membranes. The resultant perturbations of the membrane integrity may affect the conformation of the proteins inserted within the membrane, disturbing the membrane-hosted biological functions. In this study, the droplet interface bilayer (DIB), a model cell membrane, is used to examine the effects of ibuprofen, a nonsteroidal anti-inflammatory drug (NSAID), on transbilayer water permeability, which is a fundamental membrane biophysical property. Our results indicate that the presence of neutral ibuprofen (pH 3) increases the water permeability of the lipid membranes composed of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). When cholesterol is present with the DOPC, however, the water permeability is not influenced by addition of ibuprofen, regardless of the cholesterol content in DOPC. Given the fact that cholesterol is generally considered to impact packing in the hydrocarbon chain regions, our findings suggest that a potential competition between opposing effects of ibuprofen molecules and cholesterol on the hydrocarbon core environment of the phospholipid assembly may influence the overall water transport phenomena. Results from confocal Raman microspectroscopy and interfacial tensiometry show that ibuprofen molecules induce substantial structural and dynamic changes in the DOPC lipid bilayer. These results, demonstrating that the presence of ibuprofen increases the water permeability of pure DOPC but not that of DOPC-cholesterol mixtures, provide insight into the differential effect of a representative NSAID on heterogeneous biological membranes, depending upon the local composition and structure, results which will signal increased understanding of the gastrointestinal damage and toxicity induced by these molecules.

Interface-Mediated Mechanism of Action-The Root of the Cytoprotective Effect of Immediate-Release Omeprazole

Omeprazole is usually administered under an enteric coating. However, there is a Food and Drug Administration-approved strategy that enables its release in the stomach. When locally absorbed, omeprazole shows a higher efficacy and a cytoprotective effect, whose mechanism was still unknown. Therefore, we aimed to assess the effect of the absorption route on the gastric mucosa. 2D and 3D models of dipalmitoylphosphatidylcholine (DPPC) at different pH values (5.0 and 7.4) were used to mimic different absorption conditions. Several experimental techniques, namely, fluorescence studies, X-ray scattering methodologies, and Langmuir monolayers coupled with microscopy, X-ray diffraction, and infrared spectroscopy techniques, were combined with molecular dynamics simulations. The results showed that electrostatic and hydrophobic interactions between omeprazole and DPPC rearranged the conformational state of DPPC. Omeprazole intercalates among DPPC molecules, promoting domain formation with untilted phospholipids. Hence, the local release of omeprazole enables its action as a phospholipid-like drug, which can reinforce and protect the gastric mucosa.

Molecular dynamics exploring of atmosphere components interacting with lung surfactant phospholipid bilayers

Sulfur dioxide (SO₂), nitrogen oxide (NO₂) and ozone (O₃) in the atmosphere are significantly correlated with various respiratory and cardiovascular diseases. High doses of each of these gases or a mixture can change the physical and chemical properties of the lung membrane, thus leading to an increased pulmonary vascular permeability and structural failure of the alveolar cell membrane. In the present study, detailed molecular dynamic (MD) modeling was applied to investigate the effects of SO₂, NO₂, O₃ and mixtures of these gases on the dipalmitoyl phosphatidylcholine (DPPC) phospholipid bilayer. The results showed that several key physical properties, including the mass density, lipid ordering parameter, lipid diffusion, and electrostatic potential of the cell membrane, have been changed by the binding of different compounds. This resulted in significant variations and more disorder in the DPPC bilayer. The multiple analyses of membrane properties proved the toxicity of NO₂, O₃, and SO₂ to the DPPC bilayer, providing a theoretical basis for the experimental phenomenon that SO₂, NO₂ and O₃ can cause lung cell apoptosis. For the single systems, the damage to DPPC bilayer caused by O₃ was more serious than NO₂ and SO₂. More importantly, the MD simulations using the mixtures of SO₂, NO₂, and O₃ showed a much greater decline of membrane fluidity and the aggravation of membrane damage than the single systems, indicating a synergistic effect when NO₂, SO₂, and O₃ coexisted in the atmosphere, which could lead to much more severe damage and greater toxicities to the lung.

Molecular Dynamics Simulations Reveal Membrane Interactions for Poorly Water-Soluble Drugs: Impact of Bile Solubilization and Drug Aggregation

Molecular transport mechanisms of poorly soluble hydrophobic drug compounds to lipid membranes were investigated using molecular dynamics (MD) simulations. The model compound danazol was used to investigate the mechanism(s) by which bile micelles delivered it to the membrane. The interactions between lipid membrane and pure drug aggregates-in the form of amorphous aggregates and nanocrystals-were also studied. Our simulations indicate that bile micelles formed in the intestinal fluid may facilitate danazol incorporation into cellular membranes through two different mechanisms. The micelle may be acting as: i) a shuttle that presents the danazol directly to the membrane or ii) an elevator that moves the solubilized danazol with it as the colloidal structure itself becomes incorporated and solubilized within the membrane. The elevator hypothesis was supported by complementary lipid monolayer adsorption experiments. In these experiments, colloidal structures formed with simulated intestinal fluid were observed to rapidly incorporate into the monolayer. Simulations of membrane interaction with drug aggregates showed that both the amorphous aggregates and crystalline nanostructures incorporated into the membrane. However, the amorphous aggregates solubilized more quickly than the nanocrystals into the membrane, thereby improving the danazol absorption.

Comparison of Different Chitosan Lipid Nanoparticles for Improved Ophthalmic Tetrandrine Delivery: Formulation, Characterization, Pharmacokinetic and Molecular Dynamics Simulation

In this study, three different chitosan, namely carboxymethyl chitosan (CMC), hydroxypropyl chitosan (HPC) and trimethyl chitosan (TMC) were used as cationic materials to prepare tetrandrine lipid nanoparticles (TET-LNPs) for the treatment of glaucoma. In vitro drug release and pre-corneal retention were used to select the optimal chitosan. In vitro drug release curves of three kinds of LNPs showed a sustained release and TMC-TET-LNPs were the slowest. Moreover, compared with CMC-TET-LNPs and HPC-TET-LNPs, TMC-TET-LNPs had longer corneal retention time. Afterwards, the characteristics of TMC-TET-LNPs were investigated. The ocular irritation study revealed no sign of irritation in rabbit eyes. The pharmacokinetic studies showed that the area under the curve of TMC-TET-LNPs was increased by 2.03 times than TET solution (p < 0.01). Furthermore, the drug biofilm interactions were evaluated by molecular dynamics (MD) simulation. In MD simulation, the strong hydrophobic group of TET interacted with the tail of POPC, making it hard to enter the hydrophobic region of the membrane, thereby restricting TET ocular bioavailability. The experiments and MD simulation indicated that TMC-TET-LNPs had great potential for ocular administration and MD simulation could predict transmembrane transport of drugs.

In Silico Study of Membrane Lipid Composition Regulating Conformation and Hydration of Influenza Virus B M2 Channel

The proton conduction of transmembrane influenza virus B M2 (BM2) proton channel is possibly mediated by the membrane environment, but the detailed molecular mechanism is challenging to determine. In this work, how membrane lipid composition regulates the conformation and hydration of BM2 channel is elucidated in silico. The appearance of several important hydrogen-bond networks has been discovered, as the addition of negatively charged lipid palmitoyloleoyl phosphatidylglycerol (POPG) and cholesterol reduces membrane fluidity and augments membrane rigidity. A more rigid membrane environment is beneficial to expand the channel, allow more water to enter the channel, promote channel hydration, and then even affect the proton conduction facilitated by the hydrated channel. Thus, membrane environment could be identified as an important influence factor of conformation and hydration of BM2. These findings can provide a unique perspective for understanding the mechanism of membrane lipid composition regulating conformation and hydration of BM2 and have important significance to the further study of anti-influenza virus B drugs.

Membrane-Dependent Binding and Entry Mechanism of Dopamine into Its Receptor

Synaptic neurotransmission has recently been proposed to function via either a membrane-independent or a membrane-dependent mechanism, depending on the neurotransmitter type. In the membrane-dependent mechanism, amphipathic neurotransmitters first partition to the lipid headgroup region and then diffuse along the membrane plane to their membrane-buried receptors. However, to date, this mechanism has not been demonstrated for any neurotransmitter-receptor complex. Here, we combined isothermal calorimetry measurements with a diverse set of molecular dynamics simulation methods to investigate the partitioning of an amphipathic neurotransmitter (dopamine) and the mechanism of its entry into the ligand-binding site. Our results show that the binding of dopamine to its receptor is consistent with the membrane-dependent binding and entry mechanism. Both experimental and simulation results showed that dopamine favors binding to lipid membranes especially in the headgroup region. Moreover, our simulations revealed a ligand-entry pathway from the membrane to the binding site. This pathway passes through a lateral gate between transmembrane alpha-helices 5 and 6 on the membrane-facing side of the protein. All in all, our results demonstrate that dopamine binds to its receptor by a membrane-dependent mechanism, and this is complemented by the more traditional binding mechanism directly through the aqueous phase. The results suggest that the membrane-dependent mechanism is common in other synaptic receptors, too.

Molecular Dynamics Simulation of Small Molecules Interacting with Biological Membranes

Cell membranes protect and compartmentalise cells and their organelles. The semi-permeable nature of these membranes controls the exchange of solutes across their structure. Characterising the interaction of small molecules with biological membranes is critical to understanding of physiological processes, drug action and permeation, and many biotechnological applications. This review provides an overview of how molecular simulations are used to study the interaction of small molecules with biological membranes, with a particular focus on the interactions of water, organic compounds, drugs and short peptides with models of plasma cell membrane and stratum corneum lipid bilayers. This review will not delve on other types of membranes which might have different composition and arrangement, such as thylakoid or mitochondrial membranes. The application of unbiased molecular dynamics simulations and enhanced sampling methods such as umbrella sampling, metadynamics and replica exchange are described using key examples. This review demonstrates how state-of-the-art molecular simulations have been used successfully to describe the mechanism of binding and permeation of small molecules with biological membranes, as well as associated changes to the structure and dynamics of these membranes. The review concludes with an outlook on future directions in this field.

A relationship between membrane permeation and partitioning of nitroaromatic explosives and their functional groups. A computational study

Nitroaromatic explosives, such as 2,4,6-trinitrotoluene, are representative aromatic compounds, which are generally highly toxic.

Lipid Head Group Parameterization for GROMOS 54A8: A Consistent Approach with Protein Force Field Description

Membranes are a crucial component of both bacterial and mammalian cells, being involved in signaling, transport, and compartmentalization. This versatility requires a variety of lipid species to tailor the membrane’s behavior as needed, increasing the complexity of the system. Molecular dynamics simulations have been successfully applied to study model membranes and their interactions with proteins, elucidating some crucial mechanisms at the atomistic detail and thus complementing experimental techniques. An accurate description of the functional interplay of the diverse membrane components crucially depends on the selected parameters that define the adopted force field. A coherent parameterization for lipids and proteins is therefore needed. In this work, we propose and validate new lipid head group parameters for the GROMOS 54A8 force field, making use of recently published parametrizations for key chemical moieties present in lipids. We make use additionally of a new canonical set of partial charges for lipids, chosen to be consistent with the parameterization of soluble molecules such as proteins. We test the derived parameters on five phosphocholine model bilayers, composed of lipid patches four times larger than the ones used in previous studies, and run 500 ns long simulations of each system. Reproduction of experimental data like area per lipid and deuterium order parameters is good and comparable with previous parameterizations, as well as the description of liquid crystal to gel-phase transition. On the other hand, the orientational behavior of the head groups is more realistic for this new parameter set, and this can be crucial in the description of interactions with other polar molecules. For that reason, we tested the interaction of the antimicrobial peptide lactoferricin with two model membranes showing that the new parameters lead to a weaker peptide–membrane binding and give a more realistic outcome in comparing binding to antimicrobial versus mammal membranes.

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.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Exploring the membrane toxicity of decabromodiphenyl ethane (DBDPE): Based on cell membranes and lipid membranes model

Decabromodiphenyl ethane (DBDPE) is widely used in industry as an alternative to the decabromodiphenyl ether (BDEs). The large-scale use of DBDPE could lead to rapid growth of the human accumulation level of DBDPE. However, the biophysics of accumulation of DBDPE in cell membranes, as one of determinants of DBDPE metabolism is not clear. In the present study, detailed observations of cell lactate dehydrogenase (LDH) and reactive oxygen species (ROS) levels measurements proved that the DBDPE exposure to cell could result in significant cell membrane damage by concentration-dependent manners. The fluorescence anisotropy analysis supported the evidence that high concentration DBDPE bound decreased membrane fluidity significantly. Besides it, a detailed molecular dynamic (MD) simulation was approached to investigate the effects of DBDPE on the DPPC (dipalmitoyl phosphatidylcholine) phospholipid bilayer, which was constructed as the model of cell membrane. The molecular dynamic simulation revealed that DBDPE molecules can easily enter the membrane from the aqueous phase. Under the concentration of a threshold, the DBDPE molecules tended to aggregate inside the DPPC bilayer and caused pore formation. The bound of high concentration of DBDPE could result in significant variations in DPPC bilayer with a less dense, more disorder and rougher layer. The knowledge about DBDPEs interactions with lipid membranes is fundamentally essential to understand the in vivo process of DBDPE and the physical basis for the toxicity of DBDPE in cell membranes.

Licofelone-DPPC Interactions: Putting Membrane Lipids on the Radar of Drug Development

(1) Background: Membrane lipids have been disregarded in drug development throughout the years. Recently, they gained attention in drug design as targets, but they are still disregarded in the latter stages. Thus, this study aims to highlight the relevance of considering membrane lipids in the preclinical phase of drug development. (2) Methods: The interactions of a drug candidate for clinical use (licofelone) with a membrane model system made of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) were evaluated by combining Langmuir isotherms, Brewster angle microscopy (BAM), polarization-modulation infrared reflection-absorption spectroscopy (PM-IRRAS), and grazing-incidence X-ray diffraction (GIXD) measurements. (3) Results: Licofelone caused the expansion of the DPPC isotherm without changing the lipid phase transition profile. Moreover, licofelone induced the reduction of DPPC packing density, while increasing the local order of the DPPC acyl chains. (4) Conclusions: The licofelone-induced alterations in the structural organization of phosphatidylcholine monolayers may be related to its pharmacological actions. Thus, the combination of studying drug-membrane interactions with the pharmacological characterization that occurs in the preclinical stage may gather additional information about the mechanisms of action and toxicity of drug candidates. Ultimately, the addition of this innovative step shall improve the success rate of drug development.

COSMOplex: self-consistent simulation of self-organizing inhomogeneous systems based on COSMO-RS

During the past 20 years, the efficient combination of quantum chemical calculations with statistical thermodynamics by the COSMO-RS method has become an important alternative to force-field based simulations for the accurate prediction of free energies of molecules in liquid systems. While it was originally restricted to homogeneous liquids, it later has been extended to the prediction of the free energy of molecules in inhomogeneous systems such as micelles, biomembranes, or liquid interfaces, but these calculations were based on external input about the structure of the inhomogeneous system. Here we report the rigorous extension of COSMO-RS to a self-consistent prediction of the structure and the free energies of molecules in self-organizing inhomogeneous systems. This extends the application range to many new areas, such as the prediction of micellar structures and critical micelle concentrations, finite loading effects in micelles and biomembranes, the free energies and structure of liquid interfaces, microemulsions, and many more related topics, which often are of great practical importance.

In silico screening of drug-membrane thermodynamics reveals linear relations between bulk partitioning and the potential of mean force

The partitioning of small molecules in cell membranes-a key parameter for pharmaceutical applications-typically relies on experimentally available bulk partitioning coefficients. Computer simulations provide a structural resolution of the insertion thermodynamics via the potential of mean force but require significant sampling at the atomistic level. Here, we introduce high-throughput coarse-grained molecular dynamics simulations to screen thermodynamic properties. This application of physics-based models in a large-scale study of small molecules establishes linear relationships between partitioning coefficients and key features of the potential of mean force. This allows us to predict the structure of the insertion from bulk experimental measurements for more than 400 000 compounds. The potential of mean force hereby becomes an easily accessible quantity-already recognized for its high predictability of certain properties, e.g., passive permeation. Further, we demonstrate how coarse graining helps reduce the size of chemical space, enabling a hierarchical approach to screening small molecules.

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.

Accessing lipophilicity of drugs with biomimetic models: A comparative study using liposomes and micelles

Lipophilicity is a physicochemical property of crucial importance in drug discovery and drug design. Biomimetic models, such as liposomes and micelles, constitute a valuable tool for the assessment of lipophilicity through the determination of partition coefficients (log K_p). However, the lack of standardization hampers the judgment about which model or method has the best and broadest passive drug permeation predictive capacity. This work provides a comparative analysis between the methodologies based on biomimetic models to determine the partition coefficient (log K_p). For that purpose, a set of reference substances preconized by the Organization for Economic Cooperation and Development (OECD) guidelines was used. The biomimetic models employed were liposomes and micelles composed by 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) and hexadecylphosphocholine (HePC), respectively. Both lipids were used as representative phospholipids of natural membranes. The partition coefficients between biomimetic models and aqueous phases were determined by derivative spectroscopy at physiological conditions (37 °C and pH 7.4). The partition coefficients obtained using biomimetic models are quite different and more reliable than the ones obtained using an octanol/water system. Comparing the performance of the two biomimetic models, micelles revealed to be suitable only for substances with high molar absorption coefficient and log K_p > 3, but in general liposomes are the best model for accessing lipophilicity of drugs. Furthermore, a comparison between experimental data and the partition coefficients determined by the computational method COSMOmic is also provided and discussed. As a final summarizing result, a decision tree is provided in order to guide the selection of a tool for assessing the lipophilicity of drugs.

Micellization and Partition Equilibria in Mixed Nonionic/Ionic Micellar Systems: Predictions with Molecular Models

In practical applications, surfactant solutions are mostly used in mixtures of nonionic and ionic surfactants because they have improved characteristics compared to those of single surfactant solutions. By adjusting the composition of the micelles and the pH value, the solubilization of solutes can be enhanced. Nevertheless, the partitioning of solutes between nonionic/ionic mixed micelles and the aqueous phase is studied to a much lesser extent than for single surfactant solutions. Theoretical methods to predict partition equilibria in mixed micelles are of interest for screening studies. For those, the composition of the mixed micelle has to be known. Here we investigate mixtures of TX-114 (Triton X-114), Brij35 (C12E23), SDS (sodium dodecyl sulfate), and CTAB (cetyltrimethylammonium bromide). First, to investigate the surfactant compositions in the micelles, molecular dynamics (MD) self-assembly simulations were applied. Thereafter, the predictive COSMO-RS model, which applies the pseudophase approach, and its extension to anisotropic systems termed COSMOmic were compared for the prediction of partition equilibria in mixed micelles, where various molar ratios of the surfactants were considered. It could be demonstrated that both methods are applicable and lead to reasonable predictions for neutral molecules. However, taking into account the three-dimensional structure of the micelle is beneficial because the calculations with COSMOmic are in better agreement with experimental results. Because the partitioning behavior of ionizable molecules in mixed micelles is of particular interest, the partitioning of ionized isovanillin in mixed Brij35/CTAB micelles at different micelle compositions was calculated with COSMOmic. Using a thermodynamic cycle, the position-dependent pK_a of isovanillin within the micelle is calculated on the basis of COSMOmic free energy profiles. As a result, the protolytic equilibrium of isovanillin within the micelles can be taken into account, which is crucial for the reliable prediction of partition coefficients.