Structural basis of dimerization and nucleic acid binding of human DBHS proteins NONO and PSPC1.

The Drosophila behaviour/human splicing (DBHS) proteins are a family of RNA/DNA binding cofactors liable for a range of cellular processes. DBHS proteins include the non-POU domain-containing octamer-binding protein (NONO) and paraspeckle protein component 1 (PSPC1), proteins capable of forming combinatorial dimers. Here, we describe the crystal structures of the human NONO and PSPC1 homodimers, representing uncharacterized DBHS dimerization states. The structures reveal a set of conserved contacts and structural plasticity within the dimerization interface that provide a rationale for dimer selectivity between DBHS paralogues. In addition, solution X-ray scattering and accompanying biochemical experiments describe a mechanism of cooperative RNA recognition by the NONO homodimer. Nucleic acid binding is reliant on RRM1, and appears to be affected by the orientation of RRM1, influenced by a newly identified 'ß-clasp' structure. Our structures shed light on the molecular determinants for DBHS homo- and heterodimerization and provide a basis for understanding how DBHS proteins cooperatively recognize a broad spectrum of RNA targets.

Influence of force field choice on the conformational landscape of rat and human islet amyloid polypeptide.

Human islet amyloid polypeptide (hIAPP) is a naturally occurring, intrinsically disordered protein (IDP) whose abnormal aggregation into toxic soluble oligomers and insoluble amyloid fibrils is a pathological feature in type-2 diabetes. Rat IAPP (rIAPP) differs from hIAPP by only six amino acids yet has a reduced tendency to aggregate or form fibrils. The structures of the monomeric forms of IAPP are difficult to characterize due to their intrinsically disordered nature. Molecular dynamics simulations can provide a detailed characterization of the monomeric forms of rIAPP and hIAPP in near-physiological conditions. In this work, the conformational landscapes of rIAPP and hIAPP as a function of secondary structure content were predicted using well-tempered bias exchange metadynamics simulations. Several combinations of commonly used biomolecular force fields and water models were tested. The predicted conformational preferences of both rIAPP and hIAPP are typical of IDPs, exhibiting dominant random coil structures but showing a low propensity for transient α-helical conformations. Predicted nuclear magnetic resonance Cα chemical shifts reveal different preferences with each force field towards certain conformations, with AMBERff99SBnmr2/TIP4Pd showing the best agreement with the experiment. Comparisons of secondary structure content demonstrate residue-specific differences between hIAPP and rIAPP that may reflect their different aggregation propensities.

Biochemical interactions between the Atm1-like transporter from Novosphingobium aromaticivorans and heavy metals.

Novosphingobium aromaticivorans has the ability to survive in harsh environments by virtue of its suite of iron-containing oxygenases that biodegrade an astonishing array of aromatic compounds. It is also resistant to heavy metals through Atm1, an ATP-binding cassette protein that mediates active efflux of heavy metals conjugated to glutathione. However, Atm1 orthologues in higher organisms have been implicated in the intracellular transport of organic iron complexes. Our hypothesis suggests that the ability of Atm1 to remove heavy metals is related to the need for regulated iron handling in N. aromaticivorans to support high oxygenase activity. Here we provide the first data demonstrating a direct interaction between an iron-porphyrin compound (hemin) and NaAtm1. Hemin displayed considerably higher binding affinity and lower EC50 to stimulate ATP hydrolysis by Atm1 than Ag-GSH, GSSG or GSH, established substrates of the transporter. Co-incubation of NaAtm1 and hemin with Ag-GSH in ATPase assays revealed a non-competitive interaction, indicating distinct binding sites on NaAtm1 and this property was reinforced using molecular docking analysis. Our data suggests that NaAtm1 has considerable versatility in transporting organic conjugates of metals and that this versatility enables it to play roles in detoxification processes for toxic metals and in homeostasis of iron. The ability to play these distinct roles is enabled by the plasticity of the substrate binding site within the central cavity of NaAtm1.

The role of ion-lipid interactions and lipid packing in transient defects caused by phenolic compounds.

The transient disruption of membranes for the passive permeation of ions or small molecules is a complex process relevant to understanding physiological processes and biotechnology applications. Phenolic compounds are widely studied for their antioxidant and antimicrobial properties, and some of these activities are based on the interactions of the phenolic compound with membranes. Ions are ubiquitous in cells and are known to alter the structure of phospholipid bilayers. Yet, ion-lipid interactions are usually ignored when studying the membrane-altering properties of phenolic compounds. This study aims to assess the role of Ca2+ ions on the membrane-disrupting activity of two phenolic acids and to highlight the role of local changes in lipid packing in forming transient defects or pores. Results from tethered bilayer lipid membrane electrical impedance spectroscopy experiments showed that Ca2+ significantly reduces membrane disruption by caffeic acid methyl ester and caffeic acid. As phenolic acids are known metal chelators, we used UV-vis and fluorescence spectroscopy to exclude the possibility that Ca2+ interferes with membrane disruption by binding to the phenolic compound and subsequently preventing membrane binding. Molecular dynamics simulations showed that Ca2+ but not caffeic acid methyl ester or caffeic acid increases lipid packing in POPC bilayers. The combined data confirm that Ca2+ reduces the membrane-disrupting activity of the phenolic compounds, and that Ca2+-induced changes to lipid packing govern this effect. We discuss our data in the context of ion-induced pores and transient defects and how lipid packing affects membrane disruption by small molecules.

Concentration-Dependent Effect of the Steroid Drug Prednisolone on a Lung Surfactant Monolayer.

The lung surfactant monolayer (LSM) is the main barrier for particles entering the lung, including steroid drugs used to treat lung diseases. The present study combines Langmuir experiments and coarse-grained (CG) molecular dynamics simulations to investigate the concentration-dependent effect of steroid drug prednisolone on the structure and morphology of a model LSM. The surface pressure-area isotherms for the Langmuir monolayers reveal a concentration-dependent decrease in area per lipid (APL). Results from simulations at a fixed surface tension, representing inhalation and exhalation conditions, suggest that at high drug concentrations, prednisolone induces a collapse of the LSM, which is likely caused by the inability of the drug to diffuse into the bilayer. Overall, the monolayer is most susceptible to drug-induced collapse at surface tensions representing exhalation conditions. The presence of cholesterol also exacerbates the instability. The findings of this investigation might be helpful for better understanding the interaction between steroid drug prednisolone and lung surfactants in relation to off-target effects.

Calcium Ion Binding at the Lipid-Water Interface Alters the Ion Permeability of Phospholipid Bilayers.

Calcium ions (Ca2+) play a fundamental role in membrane-associated physiological processes. Ca2+ can also significantly modulate the physicochemical properties of phospholipid bilayers, but whether this occurs at physiologically relevant concentrations is difficult to determine because of the uncertainty in the reported affinity of Ca2+ for phospholipid bilayers. In this article, we determine the apparent affinity of Ca2+ for zwitterionic phospholipid bilayers using tethered bilayer lipid membranes (tBLMs) used in conjunction with swept-frequency electrical impedance spectroscopy (EIS). We report that Ca2+ binds to phospholipid bilayers at physiologically relevant concentrations and modulates membrane permeability. We present direct experimental evidence that this effect is governed by specific interactions with select lipid headgroup moieties, which is supported by data from molecular dynamics (MD) simulations. This is the first reported use of tBLM/EIS to estimate cation-membrane affinity. Combined with MD simulations, this technique provides a novel methodology to elucidate the molecular details of cation-membrane interactions at the water-phospholipid interface.

Phenolic compounds alter the ion permeability of phospholipid bilayers via specific lipid interactions.

This study aims to understand the role of specific phenolic-lipid interactions in the membrane-altering properties of phenolic compounds. We combine tethered lipid bilayer (tBLM) electrical impedance spectroscopy (EIS) with all-atom molecular dynamics (MD) simulations to study the membrane interactions of six phenolic compounds: caffeic acid methyl ester, caffeic acid, 3,4 dihydroxybenzoic acid, chlorogenic acid, syringic acid and p-coumaric acid. tBLM/EIS experiments showed that caffeic acid methyl ester, caffeic acid and 3,4 dihydroxybenzoic acid significantly increase the permeability of phospholipid bilayers to Na+ ions. In contrast, chlorogenic acid, syringic acid and p-coumaric acid showed no effect. Experiments with lipids lacking the phosphate group show a significant decrease in the membrane-altering effects indicating that specific phenolic-lipid interactions are critical in altering ion permeability. MD simulations confirm that compounds that alter ion permeability form stable interactions with the phosphate oxygen. In contrast, inactive phenolic compounds are superficially bound to the membrane surface and primarily interact with interfacial water. Our combined results show that compounds with similar structures can have very different effects on ion permeability in membranes. These effects are governed by specific interactions at the water-lipid interface and show no correlation with lipophilicity. Furthermore, none of the compounds alter the overall structure of the phospholipid bilayer as determined by area per lipid and order parameters. Based on data from this study and previous findings, we propose that phenolic compounds can alter membrane ion permeability by causing local changes in lipid packing that subsequently reduce the energy barrier for ion-induced pores.

The unusual conformation of cross-strand disulfide bonds is critical to the stability of ß-hairpin peptides.

The cross-strand disulfides (CSDs) found in ß-hairpin antimicrobial peptides (ß-AMPs) show a unique disulfide geometry that is characterized by unusual torsion angles and a short Cα-Cα distance. While the sequence and disulfide bond connectivity of disulfide-rich peptides is well studied, much less is known about the disulfide geometry found in CSDs and their role in the stability of ß-AMPs. To address this, we solved the nuclear magnetic resonance (NMR) structure of the ß-AMP gomesin (Gm) at 278, 298, and 310 K, examined the disulfide bond geometry of over 800 disulfide-rich peptides, and carried out extensive molecular dynamics (MD) simulation of the peptides Gm and protegrin. The NMR data suggests Cα-Cα distances characteristic for CSDs are independent of temperature. Analysis of disulfide-rich peptides from the Protein Data Bank revealed that right-handed and left-handed rotamers are equally likely in CSDs. The previously reported preference for right-handed rotamers was likely biased by restricting the analysis to peptides and proteins solved using X-ray crystallography. Furthermore, data from MD simulations showed that the short Cα-Cα distance is critical for the stability of these peptides. The unique disulfide geometry of CSDs poses a challenge to biomolecular force fields and to retain the stability of ß-hairpin fold over long simulation times, restraints on the torsion angles might be required.

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.

Revisiting the Interaction of Melittin with Phospholipid Bilayers: The Effects of Concentration and Ionic Strength.

Melittin is an anti-microbial peptide (AMP) and one of the most studied membrane-disrupting peptides. There is, however, a lack of accurate measurements of the concentration-dependent kinetics and affinity of binding of melittin to phospholipid membranes. In this study, we used surface plasmon resonance spectroscopy to determine the concentration-dependent effect on the binding of melittin to 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) bilayers in vesicles. Three concentration ranges were considered, and when combined, covered two orders of magnitudes (0.04 µM to 8 µM), corresponding to concentrations relevant to the membrane-disrupting and anti-microbial activities of melittin. Binding kinetics data were analysed using a 1:1 Langmuir-binding model and a two-state reaction model. Using in-depth quantitative analysis, we characterised the effect of peptide concentration, the addition of NaCl at physiological ionic strength and the choice of kinetic binding model on the reliability of the calculated kinetics and affinity of binding parameters. The apparent binding affinity of melittin for POPC bilayers was observed to decrease with increasing peptide/lipid (P/L) ratio, primarily due to the marked decrease in the association rate. At all concentration ranges, the two-state reaction model provided a better fit to the data and, thus, a more reliable estimate of binding affinity. Addition of NaCl significantly reduced the signal response during the association phase; however, no substantial effect on the binding affinity of melittin to the POPC bilayers was observed. These findings based on POPC bilayers could have important implications for our understanding of the mechanism of action of melittin on more complex model cell membranes of higher physiological relevance.

The Role of Structure and Biophysical Properties in the Pleiotropic Effects of Statins.

Statins are a class of drugs used to lower low-density lipoprotein cholesterol and are amongst the most prescribed medications worldwide. Most statins work as a competitive inhibitor of 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase (HMGR), but statin intolerance from pleiotropic effects have been proposed to arise from non-specific binding due to poor enzyme-ligand sensitivity. Yet, research into the physicochemical properties of statins, and their interactions with off-target sites, has not progressed much over the past few decades. Here, we present a concise perspective on the role of statins in lowering serum cholesterol levels, and how their reported interactions with phospholipid membranes offer a crucial insight into the mechanism of some of the more commonly observed pleiotropic effects of statin administration. Lipophilicity, which governs hepatoselectivity, is directly related to the molecular structure of statins, which dictates interaction with and transport through membranes. The structure of statins is therefore a clinically important consideration in the treatment of hypercholesterolaemia. This review integrates the recent biophysical studies of statins with the literature on the physiological effects and provides new insights into the mechanistic cause of statin pleiotropy, and prospective means of understanding the cholesterol-independent effects of statins.

Design, Synthesis and Biological Evaluation of Biphenylglyoxamide-Based Small Molecular Antimicrobial Peptide Mimics as Antibacterial Agents.

There has been an increasing interest in the development of antimicrobial peptides (AMPs) and their synthetic mimics as a novel class of antibiotics to overcome the rapid emergence of antibiotic resistance. Recently, phenylglyoxamide-based small molecular AMP mimics have been identified as potential leads to treat bacterial infections. In this study, a new series of biphenylglyoxamide-based small molecular AMP mimics were synthesised from the ring-opening reaction of N-sulfonylisatin bearing a biphenyl backbone with a diamine, followed by the conversion into tertiary ammonium chloride, quaternary ammonium iodide and guanidinium hydrochloride salts. Structure-activity relationship studies of the analogues identified the octanesulfonyl group as being essential for both Gram-positive and Gram-negative antibacterial activity, while the biphenyl backbone was important for Gram-negative antibacterial activity. The most potent analogue was identified to be chloro-substituted quaternary ammonium iodide salt 15c, which possesses antibacterial activity against both Gram-positive (MIC against Staphylococcus aureus = 8 µM) and Gram-negative bacteria (MIC against Escherichia coli = 16 µM, Pseudomonas aeruginosa = 63 µM) and disrupted 35% of pre-established S. aureus biofilms at 32 µM. Cytoplasmic membrane permeability and tethered bilayer lipid membranes (tBLMs) studies suggested that 15c acts as a bacterial membrane disruptor. In addition, in vitro toxicity studies showed that the potent compounds are non-toxic against human cells at therapeutic dosages.

A potential new, stable state of the E-cadherin strand-swapped dimer in solution.

E-cadherin is a transmembrane glycoprotein that facilitates inter-cellular adhesion in the epithelium. The ectodomain of the native structure is comprised of five repeated immunoglobulin-like domains. All E-cadherin crystal structures show the protein in one of three alternative conformations: a monomer, a strand-swapped trans homodimer and the so-called X-dimer, which is proposed to be a kinetic intermediate to forming the strand-swapped trans homodimer. However, previous studies have indicated that even once the trans strand-swapped dimer is formed, the complex is highly dynamic and the E-cadherin monomers may reorient relative to each other. Here, molecular dynamics simulations have been used to investigate the stability and conformational flexibility of the human E-cadherin trans strand-swapped dimer. In four independent, 100 ns simulations, the dimer moved away from the starting structure and converged to a previously unreported structure, which we call the Y-dimer. The Y-dimer was present for over 90% of the combined simulation time, suggesting that it represents a stable conformation of the E-cadherin dimer in solution. The Y-dimer conformation is stabilised by interactions present in both the trans strand-swapped dimer and X-dimer crystal structures, as well as additional interactions not found in any E-cadherin dimer crystal structures. The Y-dimer represents a previously unreported, stable conformation of the human E-cadherin trans strand-swapped dimer and suggests that the available crystal structures do not fully capture the conformations that the human E-cadherin trans homodimer adopts in solution.

Widespread convergence in toxin resistance by predictable molecular evolution.

The question about whether evolution is unpredictable and stochastic or intermittently constrained along predictable pathways is the subject of a fundamental debate in biology, in which understanding convergent evolution plays a central role. At the molecular level, documented examples of convergence are rare and limited to occurring within specific taxonomic groups. Here we provide evidence of constrained convergent molecular evolution across the metazoan tree of life. We show that resistance to toxic cardiac glycosides produced by plants and bufonid toads is mediated by similar molecular changes to the sodium-potassium-pump (Na(+)/K(+)-ATPase) in insects, amphibians, reptiles, and mammals. In toad-feeding reptiles, resistance is conferred by two point mutations that have evolved convergently on four occasions, whereas evidence of a molecular reversal back to the susceptible state in varanid lizards migrating to toad-free areas suggests that toxin resistance is maladaptive in the absence of selection. Importantly, resistance in all taxa is mediated by replacements of 2 of the 12 amino acids comprising the Na(+)/K(+)-ATPase H1-H2 extracellular domain that constitutes a core part of the cardiac glycoside binding site. We provide mechanistic insight into the basis of resistance by showing that these alterations perturb the interaction between the cardiac glycoside bufalin and the Na(+)/K(+)-ATPase. Thus, similar selection pressures have resulted in convergent evolution of the same molecular solution across the breadth of the animal kingdom, demonstrating how a scarcity of possible solutions to a selective challenge can lead to highly predictable evolutionary responses.

The Biological and Biophysical Properties of the Spider Peptide Gomesin.

This review summarises the current knowledge of Gomesin (Gm), an 18-residue long, cationic anti-microbial peptide originally isolated from the haemocytes of the Brazilian tarantula Acanthoscurria gomesiana. The peptide shows potent cytotoxic activity against clinically relevant microbes including Gram-positive and Gram-negative bacteria, fungi, and parasites. In addition, Gm shows in-vitro and in-vivo anti-cancer activities against several human and murine cancers. The peptide exerts its cytotoxic activity by permeabilising cell membranes, but the underlying molecular mechanism of action is still unclear. Due to its potential as a therapeutic agent, the structure and membrane-binding properties, as well as the leakage and cytotoxic activities of Gm have been studied using a range of techniques. This review provides a summary of these studies, with a particular focus on biophysical characterisation studies of peptide variants that have attempted to establish a structure-activity relationship. Future studies are still needed to rationalise the binding affinity and cell-type-specific selectivity of Gm and its variants, while more pre-clinical studies are required to develop Gm into a therapeutically useful peptide.

Characterisation of the Structure and Oligomerisation of Islet Amyloid Polypeptides (IAPP): A Review of Molecular Dynamics Simulation Studies.

Human islet amyloid polypeptide (hIAPP) is a naturally occurring, intrinsically disordered protein whose abnormal aggregation into amyloid fibrils is a pathological feature in type 2 diabetes, and its cross-aggregation with amyloid beta has been linked to an increased risk of Alzheimer's disease. The soluble, oligomeric forms of hIAPP are the most toxic to ß-cells in the pancreas. However, the structure of these oligomeric forms is difficult to characterise because of their intrinsic disorder and their tendency to rapidly aggregate into insoluble fibrils. Experimental studies of hIAPP have generally used non-physiological conditions to prevent aggregation, and they have been unable to describe its soluble monomeric and oligomeric structure at physiological conditions. Molecular dynamics (MD) simulations offer an alternative for the detailed characterisation of the monomeric structure of hIAPP and its aggregation in aqueous solution. This paper reviews the knowledge that has been gained by the use of MD simulations, and its relationship to experimental data for both hIAPP and rat IAPP. In particular, the influence of the choice of force field and water models, the choice of initial structure, and the configurational sampling method used, are discussed in detail. Characterisation of the solution structure of hIAPP and its mechanism of oligomerisation is important to understanding its cellular toxicity and its role in disease states, and may ultimately offer new opportunities for therapeutic interventions.

Interaction of Tarantula Venom Peptide ProTx-II with Lipid Membranes Is a Prerequisite for Its Inhibition of Human Voltage-gated Sodium Channel NaV1.7.

ProTx-II is a disulfide-rich peptide toxin from tarantula venom able to inhibit the human voltage-gated sodium channel 1.7 (hNaV1.7), a channel reported to be involved in nociception, and thus it might have potential as a pain therapeutic. ProTx-II acts by binding to the membrane-embedded voltage sensor domain of hNaV1.7, but the precise peptide channel-binding site and the importance of membrane binding on the inhibitory activity of ProTx-II remain unknown. In this study, we examined the structure and membrane-binding properties of ProTx-II and several analogues using NMR spectroscopy, surface plasmon resonance, fluorescence spectroscopy, and molecular dynamics simulations. Our results show a direct correlation between ProTx-II membrane binding affinity and its potency as an hNaV1.7 channel inhibitor. The data support a model whereby a hydrophobic patch on the ProTx-II surface anchors the molecule at the cell surface in a position that optimizes interaction of the peptide with the binding site on the voltage sensor domain. This is the first study to demonstrate that binding of ProTx-II to the lipid membrane is directly linked to its potency as an hNaV1.7 channel inhibitor.

Membrane-binding properties of gating modifier and pore-blocking toxins: Membrane interaction is not a prerequisite for modification of channel gating.

Many venom peptides are potent and selective inhibitors of voltage-gated ion channels, including channels that are validated therapeutic targets for treatment of a wide range of human diseases. However, the development of novel venom-peptide-based therapeutics requires an understanding of their mechanism of action. In the case of voltage-gated ion channels, venom peptides act either as pore blockers that bind to the extracellular side of the channel pore or gating modifiers that bind to one or more of the membrane-embedded voltage sensor domains. In the case of gating modifiers, it has been debated whether the peptide must partition into the membrane to reach its binding site. In this study, we used surface plasmon resonance, fluorescence spectroscopy and molecular dynamics to directly compare the lipid-binding properties of two gating modifiers (µ-TRTX-Hd1a and ProTx-I) and two pore blockers (ShK and KIIIA). Only ProTx-I was found to bind to model membranes. Our results provide further evidence that the ability to insert into the lipid bilayer is not a requirement to be a gating modifier. In addition, we characterised the surface of ProTx-I that mediates its interaction with neutral and anionic phospholipid membranes and show that it preferentially interacts with anionic lipids.

Spider peptide toxin HwTx-IV engineered to bind to lipid membranes has an increased inhibitory potency at human voltage-gated sodium channel hNaV1.7.

The human voltage-gated sodium channel sub-type 1.7 (hNaV1.7) is emerging as an attractive target for the development of potent and sub-type selective novel analgesics with increased potency and fewer side effects than existing therapeutics. HwTx-IV, a spider derived peptide toxin, inhibits hNaV1.7 with high potency and is therefore of great interest as an analgesic lead. In the current study we examined whether engineering a HwTx-IV analogue with increased ability to bind to lipid membranes would improve its inhibitory potency at hNaV1.7. This hypothesis was explored by comparing HwTx-IV and two analogues [E1PyrE]HwTx-IV (mHwTx-IV) and [E1G,E4G,F6W,Y30W]HwTx-IV (gHwTx-IV) on their membrane-binding affinity and hNaV1.7 inhibitory potency using a range of biophysical techniques including computational analysis, NMR spectroscopy, surface plasmon resonance, and fluorescence spectroscopy. HwTx-IV and mHwTx-IV exhibited weak affinity for lipid membranes, whereas gHwTx-IV showed improved affinity for the model membranes studied. In addition, activity assays using SH-SY5Y neuroblastoma cells expressing hNaV1.7 showed that gHwTx-IV has increased activity at hNaV1.7 compared to HwTx-IV. Based on these results we hypothesize that an increase in the affinity of HwTx-IV for lipid membranes is accompanied by improved inhibitory potency at hNaV1.7 and that increasing the affinity of gating modifier toxins to lipid bilayers is a strategy that may be useful for improving their potency at hNaV1.7.

The effect of hydronium ions on the structure of phospholipid membranes.

This work seeks to identify the mechanisms by which hydronium ions (H3O+) modulate the structure of phospholipid bilayers by studying the interactions of H3O+ with phospholipids at the molecular level. For this, we carried out multiple microsecond-long unrestrained molecular dynamics (MD) simulations of a POPC bilayer at different H3O+ concentrations. The results show that H3O+ accumulates at the membrane surface where it displaces water and forms strong and long-lived hydrogen bonds with the phosphate and carbonyl oxygens in phospholipids. This results in a concentration-dependent reduction of the area per lipid and an increase in bilayer thickness. This study provides an important molecular-level insight into the mechanism of how H3O+ modulates the structure of biological membranes and is a critical step towards a better understanding of the effect of low pH on mammalian and bacterial membranes.