Benzohydroxamate and nitrobenzohydroxamate affect membrane order: Correlations between spectroscopic and molecular dynamics to approach tuberculosis.

This work correlates the effects of benzohydroxamate (BH) and nitrobenzohydroxamate (NBH) anions in two membrane models which may be used for anti-tuberculosis (anti-TB) spectroscopic studies and/or computational studies. Firstly, the BH and NBH influence in the physico-chemical properties of soy asolectin (ASO)-based large multilamellar vesicles (MLVs) were evaluated by spectroscopic and calorimetric studies. In parallel, the BH and NBH interaction with a Mycobacterium tuberculosis (Mtb) inner membrane model, composed of phosphatidyl-myo-inositol-dimannoside (PIM2), was investigated by molecular dynamics (MD) simulations. Spectroscopic data showed a localization of BH close to the lipid phosphate group, while NBH was found close to the choline region. The BH ordered the ASO choline, phosphate and carbonyl regions and disrupted the acyl methylenes, reducing the membrane packing of the lipid hydrophobic region. On the other hand, NBH showed an ordering effect in all the lipid groups (polar, interface and hydrophobic ones). By MD studies, it was found that NBH enhanced the stability of the PIM2 membrane more than BH, while also being positioned closer to its mannosyl oxygens. As in ASO MLVs, BH was localized close to the PIM2 phosphate group and disrupted its acyl chains. However, higher values of lateral diffusion were observed for NBH than BH. Despite this, BH and NBH increased the membrane thickness by 35 %, which suggests a global ordering effect of both drugs. Findings of this work reinforce the accordance and complementarity between MLVs based on ASO and the PIM2 MD model results to study the drug effects in Mtb membrane properties.

The disordering effect of SARMs on a biomembrane model.

From medicine to sport, selective androgen receptor modulators (SARMs) have represented promising applications. The ability of SARMs to selectively interact with the androgen receptor (AR) indicates that this kind of molecule can interfere with numerous physiological and pathological processes controlled by the AR regulatory mechanism. However, critical concerns in relation to safety and potential side effects of SARMs remain under discussion and investigation. SARMs, being hydrophobic/organic compounds, can be subjected to hydrophobic interactions. In this perspective, we hypothesize that SARMs interact with lipid membranes, producing significant physical and chemical changes that could be associated with several effects that SARMs represent in biological systems. In this context, the effect of SARMs on lipid membranes mediated by non-specific interactions is little explored. Here, we report significant information related to the changes that ostarine, ligandrol, andarine, and cardarine produce in the thermodynamic properties of a lipid biomembrane model. Physical changes and chemical interactions of the systems were investigated by differential scanning calorimetry (DSC), dynamic light scattering (DLS), attenuated total reflectance-Fourier transform infrared spectroscopy (ATR-FTIR), and theoretical calculations implementing density functional theory (DFT). We demonstrate that ostarine, ligandrol, andarine, and cardarine can strongly interact with a lipid biomembrane model composed of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), and accordingly, these molecules can be incorporated into the polar/hydrophobic regions of the lipid bilayer. By employing theoretical calculations, we gained insights into the possible electrostatic interactions between SARMs and phospholipid molecules, enhancing our understanding of the driving forces behind the interactions of SARMs with lipid membranes. Overall, this investigation provides relevant knowledge related to the biophysical-chemical effects that SARMs produce in biomembrane models and could be of practical reference for promising applications of SARMs in medicine and sport.

The Secret Ballet Inside Multivesicular Bodies.

Lipid bilayers possess the capacity for self-assembly due to the amphipathic nature of lipid molecules, which have both hydrophobic and hydrophilic regions. When confined, lipid bilayers exhibit astonishing versatility in their forms, adopting diverse shapes that are challenging to observe through experimental means. Exploiting this adaptability, lipid structures motivate the development of bio-inspired mechanomaterials and integrated nanobio-interfaces that could seamlessly merge with biological entities, ultimately bridging the gap between synthetic and biological systems. In this work, we demonstrate how, in numerical simulations of multivesicular bodies, a fascinating evolution unfolds from an initial semblance of order toward states of higher entropy over time. We observe dynamic rearrangements in confined vesicles that reveal unexpected limit shapes of distinct geometric patterns. We identify five structures as the basic building blocks that systematically repeat under various conditions of size and composition. Moreover, we observe more complex and less frequent shapes that emerge in confined spaces. Our results provide insights into the dynamics of multivesicular systems, offering a richer understanding of how confined lipid bodies spontaneously self-organize.

Vesicle protrusion induced by antimicrobial peptides suggests common carpet mechanism for short antimicrobial peptides.

Short-cationic alpha-helical antimicrobial peptides (SCHAMPs) are promising candidates to combat the growing global threat of antimicrobial resistance. They are short-sequenced, selective against bacteria, and have rapid action by destroying membranes. A full understanding of their mechanism of action will provide key information to design more potent and selective SCHAMPs. Molecular Dynamics (MD) simulations are invaluable tools that provide detailed insights into the peptide-membrane interaction at the atomic- and meso-scale level. We use atomistic and coarse-grained MD to look into the exact steps that four promising SCHAMPs-BP100, Decoralin, Neurokinin-1, and Temporin L-take when they interact with membranes. Following experimental set-ups, we explored the effects of SCHAMPs on anionic membranes and vesicles at multiple peptide concentrations. Our results showed all four peptides shared similar binding steps, initially binding to the membrane through electrostatic interactions and then flipping on their axes, dehydrating, and inserting their hydrophobic moieties into the membrane core. At higher concentrations, fully alpha-helical peptides induced membrane budding and protrusions. Our results suggest the carpet mode of action is fit for the description of SCHAMPs lysis activity and discuss the importance of large hydrophobic residues in SCHAMPs design and activity.

Interaction of L-phenylalanine with carbonyl groups in mixed lipid membranes.

The interaction of L-Phe with the membrane components, i.e., lipids and proteins, has been discussed in the current literature due to the interest to understand the effect of single amino acids in relation to the formation of amyloid aggregates. In the present work, it is shown that L-Phe interacts with 9:1 DMPC (1,2-dimyristoyl-sn-glycero-3 phosphocholine)/DPPC (1,2-dipalmitoyl-sn-glycero-3 phosphocholine) mixtures but not in the 1:9 one. An important observation is that the interaction disappears when DPPC is replaced by diether PC (2-di-O-hexadecyl-sn-glycero-3-phosphocholine) a lipid lacking carbonyl groups (CO). This denotes that CO groups may interact specifically with L-Phe in accordance with the appearance of a new peak observed by Infrared spectroscopy (FTIR-ATR). The interaction of L-Phe affects the compressibility pattern of the 9:1 DMPC/DPPC mixture which is congruent with the changes observed by Raman spectra. The specific interaction of L-Phe with CO, propagates to phosphate and choline groups in this particular mixture as analyzed by FTIR-ATR spectroscopy and is absent when DMPC is dopped with diether PC.

Myomerger Induces Membrane Hemifusion and Regulates Fusion Pore Expansion.

Membrane fusion is a crucial mechanism in a wide variety of important events in cell biology from viral infection to exocytosis. However, despite many efforts and much progress, cell-cell fusion has remained elusive to our understanding. Along the life of the fusion pore, large conformational changes take place from the initial lipid bilayer bending, passing through the hemifusion intermediates, and ending with the formation of the first nascent fusion pore. In this sense, computer simulations are an ideal technique for describing such complex lipid remodeling at the molecular level. In this work, we studied the role played by the muscle-specific membrane protein Myomerger during the formation of the fusion pore. We have conducted µs length atomistic and coarse-grained molecular dynamics, together with free-energy calculations using ad hoc collective variables. Our results show that Myomerger favors the hemifusion diaphragm-stalk transition, reduces the nucleation-expansion energy difference, and promotes the formation of nonenlarging fusion pores.

Aß Amyloid Fibers Drastically Alter the Topography and Mechanical Properties of Lipid Membranes.

Alzheimer's disease (AD) is related to the fibrillation of the Aß peptides at neuronal membranes, a process that depends on the lipid composition and may impart different physical states to the membrane. In the present work, we study the properties of the Aß peptide when mixed with a zwitterionic lipid (DMPC), using the Langmuir monolayer technique as an approach to control membrane physical conditions. First, we build on previous characterizations of pure Aß monolayers and observe that, in addition to high shear, these films present a pronounced compressional hysteresis. When Aß is assembled with DMPC in a binary film, the resulting membranes become heterogeneous, with a peptide-enriched phase distributed in a network-like pattern, and they exhibit a lateral transition that depends on the Aß content. At lower peptide proportions, the films segregate into two well-defined phases: one consisting of lipids and another enriched with peptides. The reflectivity of these phases differs from that obtained for pure Aß films. Thus, the formed fibers effectively cover most of the interface area and remain stable at higher pressures (from 20 to 30 mN m-1 depending on Aß content) compared to pure peptide films (17 mN m-1). Furthermore, such structures induce a compressional hysteresis in the film, similar to that of pure peptide films (which is nonexistent in the pure lipid monolayer), even at low peptide proportions. We claim that the mechanical properties at the interface are governed by the size of the fibril-like structures. Based on the low molar fractions and surface packing at which these phenomena were observed, we postulate that as a consequence of peptide intermolecular interactions, Aß may have drastic effects on the molecular arrangement and mechanical properties of a lipid membrane.

Human HSP70-escort protein 1 (hHep1) interacts with negatively charged lipid bilayers and cell membranes.

Human Hsp70-escort protein 1 (hHep1) is a cochaperone that assists in the function and stability of mitochondrial HSPA9. Similar to HSPA9, hHep1 is located outside the mitochondria and can interact with liposomes. In this study, we further investigated the structural and thermodynamic behavior of interactions between hHep1 and negatively charged liposomes, as well as interactions with cellular membranes. Our results showed that hHep1 interacts peripherally with liposomes formed by phosphatidylserine and cardiolipin and remains partially structured, exhibiting similar affinities for both. In addition, after being added to the cell membrane, recombinant hHep1 was incorporated by cells in a dose-dependent manner. Interestingly, the association of HSPA9 with hHep1 improved the incorporation of these proteins into the lipid bilayer. These results demonstrated that hHep1 can interact with lipids also present in the plasma membrane, indicating roles for this cochaperone outside of mitochondria.

Out-of-plane deformability and its coupling with electrostatics in biomembranes.

Cell membranes are quasi-bidimensional soft systems formed by multipoles in an ordered array that can be polarized in an electric field. Consequently, electrostatic potentials emerge inside membranes, and membranes respond to external electric fields. From a mechanical perspective, membranes can be easily compressed-expanded, laterally deformed, and curved. Bending is particularly easy, and this kind of deformation translates to changes in the relative positions of the negative and positive charges, leading to strain gradient-induced polarization. Conversely, an external electric field gradient will exert a bending stress that translates to mechanical membrane deformation. These phenomena are described through membrane flexoelectricity. Here, we describe this property in lipid bilayers and cell membranes and summarize the studies in the field with emphasis on the effects promoted by membrane asymmetry.

Molecular dynamics insights on temperature and pressure effects on electroporation.

Electroporation is a cell-level phenomenon caused by an ionic imbalance in the membrane, being of great relevance in various fields of knowledge. A dependence of the pore formation kinetics on the environmental conditions (temperature and pressure) of the cell membrane has already been reported, but further clarification regarding how these variables affect the pore formation/resealing dynamics and the transport of molecules through the membrane is still lacking. The objective of the present study was to investigate the temperature (288-348 K) and pressure (1-5000 atm) effects on the electroporation kinetics using coarse-grained molecular dynamics simulations. Results shown that the time for pore formation and resealing increased with pressure and decreased with temperature, whereas the maximum pore radius increased with temperature and decreased with pressure. This behavior influenced the ion migration through the bilayer, and the higher ionic mobility was obtained in the 288 K/1000 atm simulations, i.e., a combination of low temperature and (not excessively) high pressure. These results were used to discuss some experimental observations regarding the extraction of intracellular compounds applying this technique. This study contributes to a better understanding of electroporation under different thermodynamic conditions and to an optimal selection of processing parameters in practical applications which exploit this phenomenon.

Interaction of Spike protein and lipid membrane of SARS-CoV-2 with Ursodeoxycholic acid, an in-silico analysis.

Numerous repositioned drugs have been sought to decrease the severity of SARS-CoV-2 infection. It is known that among its physicochemical properties, Ursodeoxycholic Acid (UDCA) has a reduction in surface tension and cholesterol solubilization, it has also been used to treat cholesterol gallstones and viral hepatitis. In this study, molecular docking was performed with the SARS-CoV-2 Spike protein and UDCA. In order to confirm this interaction, we used Molecular Dynamics (MD) in "SARS-CoV-2 Spike protein-UDCA". Using another system, we also simulated MD with six UDCA residues around the Spike protein at random, naming this "SARS-CoV-2 Spike protein-6UDCA". Finally, we evaluated the possible interaction between UDCA and different types of membranes, considering the possible membrane conformation of SARS-CoV-2, this was named "SARS-CoV-2 membrane-UDCA". In the "SARS-CoV-2 Spike protein-UDCA", we found that UDCA exhibits affinity towards the central region of the Spike protein structure of - 386.35 kcal/mol, in a region with 3 alpha helices, which comprises residues from K986 to C1032 of each monomer. MD confirmed that UDCA remains attached and occasionally forms hydrogen bonds with residues R995 and T998. In the presence of UDCA, we observed that the distances between residues atoms OG1 and CG2 of T998 in the monomers A, B, and C in the prefusion state do not change and remain at 5.93 ± 0.62 and 7.78 ± 0.51 Å, respectively, compared to the post-fusion state. Next, in "SARS-CoV-2 Spike protein-6UDCA", the three UDCA showed affinity towards different regions of the Spike protein, but only one of them remained bound to the region between the region's heptad repeat 1 and heptad repeat 2 (HR1 and HR2) for 375 ps of the trajectory. The RMSD of monomer C was the smallest of the three monomers with a value of 2.89 ± 0.32, likewise, the smallest RMSF was also of the monomer C (2.25 ± 056). In addition, in the simulation of "SARS-CoV-2 membrane-UDCA", UDCA had a higher affinity toward the virion-like membrane; where three of the four residues remained attached once they were close (5 Å, to the centre of mass) to the membrane by 30 ns. However, only one of them remained attached to the plasma-like membrane and this was in a cluster of cholesterol molecules. We have shown that UDCA interacts in two distinct regions of Spike protein sequences. In addition, UDCA tends to stay bound to the membrane, which could potentially reduce the internalization of SARS-CoV-2 in the host cell.

Myristoylation and its effects on the human Golgi Reassembly and Stacking Protein 55.

GRASP55 is a myristoylated protein localized in the medial/trans-Golgi faces and involved in the Golgi structure maintenance and the regulation of unconventional secretion pathways. It is believed that GRASP55 achieves its main functionalities in the Golgi organization by acting as a tethering factor. When bound to the lipid bilayer, its orientation relative to the membrane surface is restricted to determine its proper trans-oligomerization. Despite the paramount role of myristoylation in GRASP function, the impact of such protein modification on the membrane-anchoring properties and the structural organization of GRASP remains elusive. Here, an optimized protocol for the myristoylation in E. coli of the membrane-anchoring domain of GRASP55 is presented. The biophysical properties of the myristoylated/non-myristoylated GRASP55 GRASP domain were characterized in a membrane-mimicking micellar environment. Although myristoylation did not cause any impact on the protein's secondary structure, according to our circular dichroism data, it had a significant impact on the protein's thermal stability and solubility. Electrophoresis of negatively charged liposomes incubated with the two GRASP55 constructions showed different electrophoretic mobility for the myristoylated anchored protein only, thus demonstrating that myristoylation is essential for the biological membrane anchoring. Molecular dynamics simulations were used to further explore the anchoring process in determining the restricted orientation of GRASPs in the membrane.

Cholesterol plays a decisive role in tetraspanin assemblies during bilayer deformations.

The tetraspanin family plays key roles in many physiological processes, such as, tumour invasion, cell motility, virus infection, cell attachment and entry. Tetraspanins function as molecular scaffolds organized in microdomains with interesting downstream cellular consequences. However, despite their relevance in human physiology, the precise mechanisms of their various functions remain elusive. In particular, the full-length CD81 tetraspanin has interesting cholesterol-related properties that modulate its activity in cells. In this work, we study the opening transition of CD81 under different conditions. We propose that such conformational change is a collaborative process enhanced by simultaneous interactions between multiple identical CD81 tetraspanins. With molecular dynamics simulations we describe the crucial role of a ternary lipid bilayer with cholesterol in CD81 conformational dynamics, observing two emergent properties: first, clusters of CD81 collectively segregate one tetraspanin while favouring one opening transition, second, cumulative cholesterol sequestering by CD81 tetraspanins inhibits large membrane deformations due to local density variations.

Peptide VSAK maintains tissue glucose uptake and attenuates pro-inflammatory responses caused by LPS in an experimental model of the systemic inflammatory response syndrome: a PET study.

The present investigation using Positron Emission Tomography shows how peptide VSAK can reduce the detrimental effects produced by lipopolysaccharides in Dutch dwarf rabbits, used to develop the Systemic Inflammatory Response Syndrome (SIRS). Animals concomitantly treated with lipopolysaccharides (LPS) and peptide VSAK show important protection in the loss of radiolabeled-glucose uptake observed in diverse organs when animals are exclusively treated with LPS. Treatment with peptide VSAK prevented the onset of changes in serum levels of glucose and insulin associated with the establishment of SIRS and the insulin resistance-like syndrome. Treatment with peptide VSAK also allowed an important attenuation in the circulating levels of pro-inflammatory molecules in LPS-treated animals. As a whole, our data suggest that peptide VSAK might be considered as a candidate in the development of new therapeutic possibilities focused on mitigating the harmful effects produced by lipopolysaccharides during the course of SIRS.

A certain proportion of docosahexaenoic acid tends to revert structural and dynamical effects of cholesterol on lipid membranes.

This work investigates how docosahexaenoic acid (DHA) modifies the effect of Cholesterol (Chol) on the structural and dynamical properties of dipalmitoylphosphatidylcholine (DPPC) membrane. We employ low-cost and non-invasive methods: zeta potential (ZP), conductivity, density, and ultrasound velocity, complemented by molecular dynamics simulations. Our studies reveal that 30% of DHA added to the DPPC-Chol system tends to revert Chol action on a model lipid bilayer. Results obtained in this work shed light on the effect of polyunsaturated fatty acids - particularly DHA - on lipid membranes, with potential preventive applications in many diseases, e.g. neuronal as, Alzheimer's disease, and viral, as Covid-19.

The Cytotoxic Effect of α-Synuclein Aggregates.

Parkinson's disease is a neurodegenerative disorder involving a functional protein, α-synuclein, whose primary function is related to vesicle trafficking. However, α-synuclein is prone to form aggregates, and these inclusions, known as Lewy bodies, are the hallmark of Parkinson's disease. α-synuclein can alter its conformation and acquire aggregating capacity, forming aggregates containing ß-sheets. This protein's pathogenic importance is based on its ability to form oligomers that impair synaptic transmission and neuronal function by increasing membrane permeability and altering homeostasis, generating a deleterious effect over cells. First, we establish that oligomers interfere with the mechanical properties of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) membrane, as demonstrated by nanoindentation curves. In contrast, nanoindentation revealed that the α-synuclein monomer's presence leads to a much more resistant lipid bilayer. Moreover, the oligomers' interaction with cell membranes can promote lactate dehydrogenase (LDH) release, suggesting the activation of cytotoxic events.

Physicochemical-guided design of cathelicidin-derived peptides generates membrane active variants with therapeutic potential.

The spread of multi-drug resistance and the slow pace at which antibiotics come onto the market are undermining our ability to treat human infections, leading to high mortality rates. Aiming to overcome this global crisis, antimicrobial peptides are considered promising alternatives to counter bacterial infections with multi-drug resistant bacteria. The cathelicidins comprise a well-studied class of AMPs whose members have been used as model molecules for sequence modifications, aiming at enhanced biological activities and stability, along with reduced toxic effects on mammalian cells. Here, we describe the antimicrobial activities, modes of action and structural characterization of two novel cathelicidin-like peptides, named BotrAMP14 and CrotAMP14, which were re-designed from snake batroxicidin and crotalicidin, respectively. BotrAMP14 and CrotAMP14 showed broad-spectrum antibacterial activity against susceptible microorganisms and clinical isolates with minimal inhibitory concentrations ranging from 2-35.1 µM. Moreover, both peptides had low cytotoxicity against Caco-2 cells in vitro. In addition, in vivo toxicity against Galleria mellonella moth larvae revealed that both peptides led to>76% larval survival after 144 h. Microscopy studies suggest that BotrAMP14 and CrotAMP14 destabilize E. coli membranes. Furthermore, circular dichroism and molecular dynamics simulations indicate that, in a membrane-like environment, both peptides adopt α-helical structures that interact with bilayer phospholipids through hydrogen bonds and electrostatic interaction. Thus, we concluded that BotrAMP14 and CrotAMP14 are helical membrane active peptides, with similar antibacterial properties but lower cytotoxicity than the larger parent peptides batroxicidin and crotalicidin, having advantages for drug development strategies.

Order-disorder skewness in alpha-synuclein: a key mechanism to recognize membrane curvature.

Currently, membrane curvature is understood as an active mechanism to control cells spatial organization and activity. Protein processes involved in sensing and generating curvature are therefore of major interest. In this work, we have studied α-synuclein interactions with a model lipid bilayer, inducing curvature in a controlled manner and describing protein responses at molecular level. We show that the intrinsically disordered region of α-synuclein binds to the bilayer as an acknowledgment to the induced curvature, a mechanism used by the interacting protein-membrane assembly to relieve free energy. We have calculated free energies for bending the bilayer with α-synuclein adsorbed on the surface and we have established the crucial role of the intrinsically disordered region, suggesting that a dynamic order/disorder interplay takes place as the bilayer reorganizes to bend.

Localization of Annexin A6 in Matrix Vesicles During Physiological Mineralization.

Annexin A6 (AnxA6) is the largest member of the annexin family of proteins present in matrix vesicles (MVs). MVs are a special class of extracellular vesicles that serve as a nucleation site during cartilage, bone, and mantle dentin mineralization. In this study, we assessed the localization of AnxA6 in the MV membrane bilayer using native MVs and MV biomimetics. Biochemical analyses revealed that AnxA6 in MVs can be divided into three distinct groups. The first group corresponds to Ca2+-bound AnxA6 interacting with the inner leaflet of the MV membrane. The second group corresponds to AnxA6 localized on the surface of the outer leaflet. The third group corresponds to AnxA6 inserted in the membrane's hydrophobic bilayer and co-localized with cholesterol (Chol). Using monolayers and proteoliposomes composed of either dipalmitoylphosphatidylcholine (DPPC) to mimic the outer leaflet of the MV membrane bilayer or a 9:1 DPPC:dipalmitoylphosphatidylserine (DPPS) mixture to mimic the inner leaflet, with and without Ca2+, we confirmed that, in agreement with the biochemical data, AnxA6 interacted differently with the MV membrane. Thermodynamic analyses based on the measurement of surface pressure exclusion (πexc), enthalpy (ΔH), and phase transition cooperativity (Δt1/2) showed that AnxA6 interacted with DPPC and 9:1 DPPC:DPPS systems and that this interaction increased in the presence of Chol. The selective recruitment of AnxA6 by Chol was observed in MVs as probed by the addition of methyl-ß-cyclodextrin (MßCD). AnxA6-lipid interaction was also Ca2+-dependent, as evidenced by the increase in πexc in negatively charged 9:1 DPPC:DPPS monolayers and the decrease in ΔH in 9:1 DPPC:DPPS proteoliposomes caused by the addition of AnxA6 in the presence of Ca2+ compared to DPPC zwitterionic bilayers. The interaction of AnxA6 with DPPC and 9:1 DPPC:DPPS systems was distinct even in the absence of Ca2+ as observed by the larger change in Δt1/2 in 9:1 DPPC:DPPS vesicles as compared to DPPC vesicles. Protrusions on the surface of DPPC proteoliposomes observed by atomic force microscopy suggested that oligomeric AnxA6 interacted with the vesicle membrane. Further work is needed to delineate possible functions of AnxA6 at its different localizations and ways of interaction with lipids.

Evolution of Bonding during the Insertion of Anionic Ibuprofen into Model Cell Membranes.

Descriptors of chemical bonding derived from five different analysis tools based on quantum mechanics (natural charges, electron density differences, atoms in molecules (AIM), natural bond orbitals (NBO), and non-covalent interactions (NCI) index) consistently afford a picture of a wall of weak, non-covalent intermolecular interactions separating anionic Ibuprofen from the environment. This wall, arising from the cumulative effect of a multitude of individual weak charge transfer interactions to the interstitial region between fragments, stabilizes the drug at all equilibrium positions in the free energy profile for its insertion into model cell membranes. The formal charge in anionic Ibuprofen strengthens all intermolecular interactions, having a particularly strong effect in the network of water to water hydrogen bonds in the solvent. Electron redistribution during the insertion process leads to a sensible reduction of electron delocalization in both the -CO2- group and the aromatic ring of Ibuprofen. Here, we conclusively show that, despite their purely classical origin, randomly chosen configurations from molecular dynamics simulations provide deep insight into the purely quantum nature of bonding interactions.