Systematic evaluation of lecithin:cholesterol acyltransferase binding sites in apolipoproteins via peptide based nanodiscs: regulatory role of charged residues at positions 4 and 7.

Lecithin:cholesterol acyltransferase (LCAT) exhibits α-activity on high-density and ß-activity on low-density lipoproteins. However, the molecular determinants governing LCAT activation by different apolipoproteins remain elusive. Uncovering these determinants would offer the opportunity to design and explore advanced therapies against dyslipidemias. Here, we have conducted coarse-grained and all-atom molecular dynamics simulations of LCAT with nanodiscs made with α-helical amphiphilic peptides either derived from apolipoproteins A1 and E (apoA1 and apoE) or apoA1 mimetic peptide 22A that was optimized to activate LCAT. This study aims to explore what drives the binding of peptides to our previously identified interaction site in LCAT. We hypothesized that this approach could be used to screen for binding sites of LCAT in different apolipoproteins and would provide insights to differently localized LCAT activities. Our screening approach was able to discriminate apoA1 helixes 4, 6, and 7 as key contributors to the interaction with LCAT supporting the previous research data. The simulations provided detailed molecular determinants driving the interaction with LCAT: the formation of hydrogen bonds or salt bridges between peptides E4 or D4 and LCAT S236 or K238 residues. Additionally, salt bridging between R7 and D73 was observed, depending on the availability of R7. Expanding our investigation to diverse plasma proteins, we detected novel LCAT binding helixes in apoL1, apoB100, and serum amyloid A. Our findings suggest that the same binding determinants, involving E4 or D4 -S236 and R7-D73 interactions, influence LCAT ß-activity on low-density lipoproteins, where apoE and or apoB100 are hypothesized to interact with LCAT.

Development of Robust Cationic Light-Activated Thermosensitive Liposomes: Choosing the Right Lipids.

Extensive research has been conducted on cationic light-activated thermosensitive liposomes (CLTSLs) as a means for site-specific and controlled drug release; however, less attention has been given to the stability of these nanoparticles. Selecting the appropriate lipids is crucial for the development of a stable and responsive system. In this study, we investigated the impact of various lipids on the physical properties of cationic light-activated liposomes. Incorporating poly(ethylene glycol) PEG molecules resulted in uniform liposomes with low polydispersity index, while the addition of unsaturated lipid (DOTAP) resulted in extremely leaky liposomes, with almost 80% release in just 10 min of incubation at body temperature. Conversely, the inclusion of cholesterol in the formulation increased liposome stability too much and decreased their sensitivity to stimuli-responsive release, with only 14% release after 2 min of light exposure. To achieve stable and functional CLTSL, we substituted an equivalent amount of unsaturated lipid with a saturated lipid (DPTAP), resulting in stable liposomes at body temperature that were highly responsive to light, releasing 90% of their content in 10 s of light exposure. We also conducted two atomistic molecular dynamics simulations using lipid compositions with saturated and unsaturated lipids to investigate the effect of lipid composition on the dynamical properties of the liposomal lipid bilayer. Our findings suggest that the nature of lipids used to prepare liposomes significantly affects their properties, especially when the drug loading needs to be stable but triggered drug release properties are required at the same time. Selecting the appropriate lipids in the right amount is therefore essential for the preparation of liposomes with desirable properties.

Mechanistic Insights into the Activation of Lecithin-Cholesterol Acyltransferase in Therapeutic Nanodiscs Composed of Apolipoprotein A-I Mimetic Peptides and Phospholipids.

The mechanistic details behind the activation of lecithin-cholesterol acyltransferase (LCAT) by apolipoprotein A-I (apoA-I) and its mimetic peptides are still enigmatic. Resolving the fundamental principles behind LCAT activation will facilitate the design of advanced HDL-mimetic therapeutic nanodiscs for LCAT deficiencies and coronary heart disease and for several targeted drug delivery applications. Here, we have combined coarse-grained molecular dynamics simulations with complementary experiments to gain mechanistic insight into how apoA-Imimetic peptide 22A and its variants tune LCAT activity in peptide-lipid nanodiscs. Our results highlight that peptide 22A forms transient antiparallel dimers in the rim of nanodiscs. The dimerization tendency considerably decreases with the removal of C-terminal lysine K22, which has also been shown to reduce the cholesterol esterification activity of LCAT. In addition, our simulations revealed that LCAT prefers to localize to the rim of nanodiscs in a manner that shields the membrane-binding domain (MBD), αA-αA', and the lid amino acids from the water phase, following previous experimental evidence. Meanwhile, the location and conformation of LCAT in the rim of nanodiscs are spatially more restricted when the active site covering the lid of LCAT is in the open form. The average location and spatial dimensions of LCAT in its open form were highly compatible with the electron microscopy images. All peptide 22A variants studied here had a specific interaction site in the open LCAT structure flanked by the lid and MBD domain. The bound peptides showed different tendencies to form antiparallel dimers and, interestingly, the temporal binding site occupancies of the peptide variants affected their in vitro ability to promote LCAT-mediated cholesterol esterification.

Positive allosteric modulators of lecithin: Cholesterol acyltransferase adjust the orientation of the membrane-binding domain and alter its spatial free energy profile.

Lecithin:cholesterol acyltransferase protein (LCAT) promotes the esterification reaction between cholesterol and phospholipid-derived acyl chains. Positive allosteric modulators have been developed to treat LCAT deficiencies and, plausibly, also cardiovascular diseases in the future. The mechanism of action of these compounds is poorly understood. Here computational docking and atomistic molecular dynamics simulations were utilized to study the interactions between LCAT and the activating compounds. Results indicate that all drugs bind to the allosteric binding pocket in the membrane-binding domain in a similar fashion. The presence of the compounds in the allosteric site results in a distinct spatial orientation and sampling of the membrane-binding domain (MBD). The MBD's different spatial arrangement plausibly affects the lid's movement from closed to open state and vice versa, as suggested by steered molecular dynamics simulations.

Mechanistic Insight into How PEGylation Reduces the Efficacy of pH-Sensitive Liposomes from Molecular Dynamics Simulations.

Liposome-based drug delivery systems composed of DOPE stabilized with cholesteryl hemisuccinate (CHMS) have been proposed as a drug delivery mechanism with pH-triggered release as the anionic form (CHSa) is protonated (CHS) at reduced pH; PEGylation is known to decrease this pH sensitivity. In this manuscript, we set out to use molecular dynamics (MD) simulations with a model with all-atom resolution to provide insight into why incorporation of poly(ethyleneglycol) (PEG) into DOPE-CHMS liposomes reduces their pH sensitivity; we also address two additional questions: (1) How CHSa stabilizes DOPE bilayers into a lamellar conformation at a physiological pH of 7.4? and (2) how the change from CHSa to CHS at acidic pH triggers the destabilization of DOPE bilayers? We found that (A) CHSa stabilizes the DOPE lipid membrane by increasing the hydrophilicity of the bilayer surface, (B) when CHSa changes to CHS by pH reduction, DOPE bilayers are destabilized due to a reduction in bilayer hydrophilicity and a reduction in the area per lipid, and (C) PEG stabilizes DOPE bilayers into the lamellar phase, thus reducing the pH sensitivity of the liposomes by increasing the area per lipid through penetration into the bilayer, which is our main focus.

Rigorous Computational Study Reveals What Docking Overlooks: Double Trouble from Membrane Association in Protein Kinase C Modulators.

Increasing protein kinase C (PKC) activity is of potential therapeutic value. Its activation involves an interaction between the C1 domain and diacylglycerol (DAG) at intracellular membrane surfaces; DAG mimetics hold promise as new drugs. We previously developed the isophthalate derivative HMI-1a3, an effective but highly lipophilic (clogP = 6.46) DAG mimetic. Although a less lipophilic pyrimidine analog, PYR-1gP (clogP = 3.30), gave positive results in computational docking, it unexpectedly presented greatly diminished binding to PKC in vitro. Through more rigorous computational molecular modeling, we reveal that, unlike HMI-1a3, PYR-1gP forms an intramolecular hydrogen bond, which both obstructs binding and reorients PYR-1gP in the membrane in a fashion that prevents it from correctly accessing the PKC C1 domain. Our results highlight the great value of molecular dynamics simulations as a key component for the drug design process of ligands targeting weakly membrane-associated proteins, where simulation in the relevant membrane environment is crucial for obtaining biologically applicable results.

Interaction of lecithin:cholesterol acyltransferase with lipid surfaces and apolipoprotein A-I-derived peptides.

LCAT is an enzyme responsible for the formation of cholesteryl esters from unesterified cholesterol (UC) and phospholipid (PL) molecules in HDL particles. However, it is poorly understood how LCAT interacts with lipoproteins and how apoA-I activates it. Here we have studied the interactions between LCAT and lipids through molecular simulations. In addition, we studied the binding of LCAT to apoA-I-derived peptides, and their effect on LCAT lipid association-utilizing experiments. Results show that LCAT anchors itself to lipoprotein surfaces by utilizing nonpolar amino acids located in the membrane-binding domain and the active site tunnel opening. Meanwhile, the membrane-anchoring hydrophobic amino acids attract cholesterol molecules next to them. The results also highlight the role of the lid-loop in the lipid binding and conformation of LCAT with respect to the lipid surface. The apoA-I-derived peptides from the LCAT-activating region bind to LCAT and promote its lipid surface interactions, although some of these peptides do not bind lipids individually. The transfer free-energy of PL from the lipid bilayer into the active site is consistent with the activation energy of LCAT. Furthermore, the entry of UC molecules into the active site becomes highly favorable by the acylation of SER181.

The biophysical properties of ethanolamine plasmalogens revealed by atomistic molecular dynamics simulations.

Given the importance of plasmalogens in cellular membranes and neurodegenerative diseases, a better understanding of how plasmalogens affect the lipid membrane properties is needed. Here we carried out molecular dynamics simulations to study a lipid membrane comprised of ethanolamine plasmalogens (PE-plasmalogens). We compared the results to the PE-diacyl counterpart and palmitoyl-oleyl-phosphatidylcholine (POPC) bilayers. Results show that PE-plasmalogens form more compressed, thicker, and rigid lipid bilayers in comparison with the PE-diacyl and POPC membranes. The results also point out that the vinyl-ether linkage increases the ordering of sn-1 chain substantially and the ordering of the sn-2 chain to a minor extent. Further, the vinyl-ether linkage changes the orientation of the lipid head group, but it does not cause changes in the head group and glycerol backbone tilt angles with respect to the bilayer normal. The vinyl-ether linkage also packs the proximal regions of the sn-1 and sn-2 chains more closely together which also decreases the distance between the rest of the sn-1 and sn-2 chains.

A computational approach for studying antibody-antigen interactions without prior structural information: the anti-testosterone binding antibody as a case study.

Given the increasing exploitation of antibodies in different contexts such as molecular diagnostics and therapeutics, it would be beneficial to unravel the atomistic level properties of antibody-antigen complexes with the help of computational modeling. Thus, here we have studied the feasibility of computational tools to gather atomic scale information regarding the antibody-antigen complexes solely starting from an amino acid sequence. First, we constructed a homology model for the anti-testosterone binding antibody based on the knowledge based classification of complementary determining regions (CDRs) and implicit solvent molecular dynamics simulations. To further examine whether the generated homology model is suitable for studying antibody-antigen interactions, docking calculations were carried out followed by binding free-energy simulations. Our results indicate that with the antibody modeling approach presented here it is possible to construct accurate homology models for antibodies which correctly describes the antibody-antigen interactions, and produces absolute binding free-energies that are comparable with experimental values. In addition, our simulations suggest that the conformations of complementary determining regions (CDRs) may considerably change from the X-ray configuration upon solvation. In conclusion, here we have introduced an antibody modeling workflow that can be used in studying the interactions between antibody and antigen solely based on an amino acid sequence, which in turn provides novel opportunities to tune the properties of antibodies in different applications. Proteins 2017; 85:322-331. © 2016 Wiley Periodicals, Inc.

Interactions between Chloramphenicol, Carrier Polymers, and Bacteria-Implications for Designing Electrospun Drug Delivery Systems Countering Wound Infection.

Antibacterial drug-loaded electrospun nano- and microfibrous dressings are of major interest as novel topical drug delivery systems in wound care. In this study, chloramphenicol (CAM)-loaded polycaprolactone (PCL) and PCL/poly(ethylene oxide) (PEO) fiber mats were electrospun and characterized in terms of morphology, drug distribution, physicochemical properties, drug release, swelling, cytotoxicity, and antibacterial activity. Computational modeling together with physicochemical analysis helped to elucidate possible interactions between the drug and carrier polymers. Strong interactions between PCL and CAM together with hydrophobicity of the system resulted in much slower drug release compared to the hydrophilic ternary system of PCL/PEO/CAM. Cytotoxicity studies confirmed safety of the fiber mats to murine NIH 3T3 cells. Disc diffusion assay demonstrated that both fast and slow release fiber mats reached effective concentrations and had similar antibacterial activity. A biofilm formation assay revealed that both blank matrices are good substrates for the bacterial attachment and formation of biofilm. Importantly, prolonged release of CAM from drug-loaded fibers helps to avoid biofilm formation onto the dressing and hence avoids the treatment failure.

How anacetrapib inhibits the activity of the cholesteryl ester transfer protein? Perspective through atomistic simulations.

Cholesteryl ester transfer protein (CETP) mediates the reciprocal transfer of neutral lipids (cholesteryl esters, triglycerides) and phospholipids between different lipoprotein fractions in human blood plasma. A novel molecular agent known as anacetrapib has been shown to inhibit CETP activity and thereby raise high density lipoprotein (HDL)-cholesterol and decrease low density lipoprotein (LDL)-cholesterol, thus rendering CETP inhibition an attractive target to prevent and treat the development of various cardiovascular diseases. Our objective in this work is to use atomistic molecular dynamics simulations to shed light on the inhibitory mechanism of anacetrapib and unlock the interactions between the drug and CETP. The results show an evident affinity of anacetrapib towards the concave surface of CETP, and especially towards the region of the N-terminal tunnel opening. The primary binding site of anacetrapib turns out to reside in the tunnel inside CETP, near the residues surrounding the N-terminal opening. Free energy calculations show that when anacetrapib resides in this area, it hinders the ability of cholesteryl ester to diffuse out from CETP. The simulations further bring out the ability of anacetrapib to regulate the structure-function relationships of phospholipids and helix X, the latter representing the structural region of CETP important to the process of neutral lipid exchange with lipoproteins. Altogether, the simulations propose CETP inhibition to be realized when anacetrapib is transferred into the lipid binding pocket. The novel insight gained in this study has potential use in the development of new molecular agents capable of preventing the progression of cardiovascular diseases.

Far-red light-triggered cargo release from liposomes b ound to a photosensitizer-cellulose nanofiber hydrogel.

In our research we used the anionic nanofibrillar cellulose (ANFC) as a platform for far-red light-induced release of cargo from liposomes. In contrast to previous works, where photosensitizers are usually in the liposomal bilayers, we used a cellulose-binding dye. Our phthalocyanine derivative has been shown to bind very strongly to cellulose and cellulose nanofiber hydrogels, allowing us to place it outside of the liposomes. Both the sensitizer and cationic liposomes bind strongly to the ANFC after mixing, making the system easy to fabricate. Upon light activation, the photosensitizer generates reactive oxygen species (ROS) within the ANFC hydrogel, where the reactive oxygen species oxidize unsaturated lipids in the liposomal membrane, which makes the liposomes more permeable, resulting in on-demand cargo release. We were able to achieve ca. 70 % release of model hydrophilic cargo molecule calcein from the hydrogels with a relatively low dose of light (262 J/cm2) while employing the straightforward fabrication techniques. Our system was remarkably responsive to the far-red light (730 nm), enabling deep tissue penetration. Therefore, this very promising novel cellulose-immobilized photosensitizer liposomal platform could be used as a controlled drug delivery system, which can have applications in externally activated coatings or implants.

Interfacial properties of high-density lipoprotein-like lipid droplets with different lipid and apolipoprotein A-I compositions.

The surface properties of high-density lipoproteins (HDLs) are important because different enzymes bind and carry out their functions at the surface of HDL particles during metabolic processes. However, the surface properties of HDL and other lipoproteins are poorly known because they cannot be directly measured for nanoscale particles with contemporary experimental methods. In this work, we carried out coarse-grained molecular dynamics simulations to study the concentration of core lipids in the surface monolayer and the interfacial tension of droplets resembling HDL particles. We simulated lipid droplets composed of different amounts of phospholipids, cholesterol esters (CEs), triglycerides (TGs), and apolipoprotein A-Is. Our results reveal that the amount of TGs in the vicinity of water molecules in the phospholipid monolayer is 25-50% higher compared to the amount of CEs in a lipid droplet with a mixed core of an equal amount of TG and CE. In addition, the correlation time for the exchange of molecules between the core and the monolayer is significantly longer for TGs compared to CEs. This suggests that the chemical potential of TG is lower in the vicinity of aqueous phase but the free-energy barrier for the translocation between the monolayer and the core is higher compared to CEs. From the point of view of enzymatic modification, this indicates that TG molecules are more accessible from the aqueous phase. Further, our results point out that CE molecules decrease the interfacial tension of HDL-like lipid droplets whereas TG keeps it constant while the amount of phospholipids varies.

Lipid exchange mechanism of the cholesteryl ester transfer protein clarified by atomistic and coarse-grained simulations.

Cholesteryl ester transfer protein (CETP) transports cholesteryl esters, triglycerides, and phospholipids between different lipoprotein fractions in blood plasma. The inhibition of CETP has been shown to be a sound strategy to prevent and treat the development of coronary heart disease. We employed molecular dynamics simulations to unravel the mechanisms associated with the CETP-mediated lipid exchange. To this end we used both atomistic and coarse-grained models whose results were consistent with each other. We found CETP to bind to the surface of high density lipoprotein (HDL) -like lipid droplets through its charged and tryptophan residues. Upon binding, CETP rapidly (in about 10 ns) induced the formation of a small hydrophobic patch to the phospholipid surface of the droplet, opening a route from the core of the lipid droplet to the binding pocket of CETP. This was followed by a conformational change of helix X of CETP to an open state, in which we found the accessibility of cholesteryl esters to the C-terminal tunnel opening of CETP to increase. Furthermore, in the absence of helix X, cholesteryl esters rapidly diffused into CETP through the C-terminal opening. The results provide compelling evidence that helix X acts as a lid which conducts lipid exchange by alternating the open and closed states. The findings have potential for the design of novel molecular agents to inhibit the activity of CETP.

Rational design of a polysaccharide-based viral mimicry nanocomplex for potent gene silencing in inflammatory tissues.

Rational design and fabrication of small interfering RNA (siRNA) delivery system with simple production scheme, specific targeting capability, responsiveness to endogenous stimuli and potential multi-functionalities remains technically challenging. Herein, we screen and design a virus-mimicking polysaccharide nanocomplex that shows specific gene delivery capability in a selective subset of leukocytes. A virus-inspired poly (alkyl methacrylate-co-methacrylic acid) fragment was conjugated on barley ß-glucans (EEPG) to endow the nanocomplex with pH-dependent endosomal membrane destabilization capabilities, as confirmed both biologically and computationally. siRNA loaded EEPG nanocomplex is feasibly fabricated in a single-step manner, which exhibit efficient gene silencing efficacy towards Dectin-1+ monocytes/macrophages. The inherent targeting affinity and feasible gene silencing potency of EEPG nanocomplex are investigated in three independent murine inflammation models, including myocardial infarction, lung fibrosis and acute liver damage. Significant enhanced accumulation level of EEPG nanocomplex is observed in cardiac lesion site, indicating its exclusive targeting capability for ischemic heart diseases. As a proof of concept, siTGF-ß based gene therapy is confirmed in murine model with heart fibrosis. Overall, our findings suggest the designed EEPG nanocomplex is favorable for siRNA delivery, which might have translational potential as a versatile platform in inflammation-related diseases.

Expanding the Paradigm of Structure-Based Drug Design: Molecular Dynamics Simulations Support the Development of New Pyridine-Based Protein Kinase C-Targeted Agonists.

Protein kinase C (PKC) modulators hold therapeutic potential for various diseases, including cancer, heart failure, and Alzheimer's disease. Targeting the C1 domain of PKC represents a promising strategy; the available protein structures warrant the design of PKC-targeted ligands via a structure-based approach. However, the PKC C1 domain penetrates the lipid membrane during binding, complicating the design of drug candidates. The standard docking-scoring approach for PKC lacks information regarding the dynamics and the membrane environment. Molecular dynamics (MD) simulations with PKC, ligands, and membranes have been used to address these shortcomings. Previously, we observed that less computationally intensive simulations of just ligand-membrane interactions may help elucidate C1 domain-binding prospects. Here, we present the design, synthesis, and biological evaluation of new pyridine-based PKC agonists implementing an enhanced workflow with ligand-membrane MD simulations. This workflow holds promise to expand the approach in drug design for ligands targeted to weakly membrane-associated proteins.

Interfacial tension and surface pressure of high density lipoprotein, low density lipoprotein, and related lipid droplets.

Lipid droplets play a central role in energy storage and metabolism on a cellular scale. Their core is comprised of hydrophobic lipids covered by a surface region consisting of amphiphilic lipids and proteins. For example, high and low density lipoproteins (HDL and LDL, respectively) are essentially lipid droplets surrounded by specific proteins, their main function being to transport cholesterol. Interfacial tension and surface pressure of these particles are of great interest because they are related to the shape and the stability of the droplets and to protein adsorption at the interface. Here we use coarse-grained molecular-dynamics simulations to consider a number of related issues by calculating the interfacial tension in protein-free lipid droplets, and in HDL and LDL particles mimicking physiological conditions. First, our results suggest that the curvature dependence of interfacial tension becomes significant for particles with a radius of ∼5 nm, when the area per molecule in the surface region is <1.4 nm(2). Further, interfacial tensions in the used HDL and LDL models are essentially unaffected by single apo-proteins at the surface. Finally, interfacial tensions of lipoproteins are higher than in thermodynamically stable droplets, suggesting that HDL and LDL are kinetically trapped into a metastable state.

Role of neutral lipids in tear fluid lipid layer: coarse-grained simulation study.

Tear fluid lipid layer (TFLL) residing at the air-water interface of tears has been recognized to play an important role in the development of dry eye syndrome. Yet, the composition, structure, and mechanical properties of TFLL are only partly known. Here, we report results of coarse-grained simulations of a lipid layer comprising phospholipids, free fatty acids, cholesteryl esters, and triglycerides at the air-water interface to shed light on the properties of TFLL. We consider structural as well as dynamical properties of the lipid layer as a function of surface pressure. Simulations revealed that neutral lipids reside heterogeneously between phospholipids at relatively low pressures but form a separate hydrophobic phase with increasing surface pressure, transforming the initial lipid monolayer to a two-layered structure. When the model of TFLL was compared to a one-component phospholipid monolayer system, we found drastic differences in both structural and dynamical properties that explain the prominent role of neutral lipids as stabilizers of the TFLL. Based on our results, we suggest that neutral lipids are able to increase the stability of the TFLL by modulating its dynamical and structural behavior, which is important for the proper function of tear film.

Insights into the behavior of unsaturated diacylglycerols in mixed lipid bilayers in relation to protein kinase C activation-A molecular dynamics simulation study.

The lipid second messenger diacylglycerol (DAG) is known for its involvement in many types of cellular signaling, especially as an endogenous agonist for protein kinase C (PKC). Evidence has emerged that the degree of saturation of the DAG molecules can affect PKC activation. DAG molecules with different acyl chain saturation have not only been observed to induce varying extents of PKC activation, but also to express selectivity towards different PKC isozymes. Both qualities are important for precise therapeutic activation of PKC; understanding DAG behavior at the molecular level in different environments has much potential in the development of drugs to target PKC. We used molecular dynamics simulations to study the behavior of two different unsaturated DAG species in lipid environments with varying degrees of unsaturation. We focus on phosphatidylethanolamine (PE) instead of phosphatidylcholine (PC) to more accurately model the relevant biomembranes. The effect of cholesterol (CHOL) on these systems was also explored. We found that both high level of unsaturation in the acyl chains of the DAG species and presence of CHOL in the surrounding membrane increase DAG molecule availability at the lipid-water interface. This can partially explain the previously observed differences in PKC activation strength and specificity, the complete mechanism is, however, likely to be more complex. Our simulations coupled with the current understanding of lipids highlight the need for more simulations of biologically accurate lipid environments in order to determine the correct correlations between molecular mechanisms and biological behavior when studying PKC activation.

Membrane simulations mimicking acidic pH reveal increased thickness and negative curvature in a bilayer consisting of lysophosphatidylcholines and free fatty acids.

Phospholipids are key components of biological membranes and their lipolysis with phospholipase A(2) (PLA(2)) enzymes occurs in different cellular pH environments. Since no studies are available on the effect of pH on PLA(2)-modified phospholipid membranes, we performed 50-ns atomistic molecular dynamics simulations at three different pH conditions (pH 9.0, 7.5, and 5.5) using a fully PLA(2)-hydrolyzed phosphatidylcholine (PC) bilayer which consists solely of lysophosphatidylcholine and free fatty acid molecules. We found that a decrease in pH results in lateral squeezing of the membrane, i.e. in decreased surface area per headgroup. Thus, at the decreased pH, the lipid hydrocarbon chains had larger S(CD) order parameter values, and also enhanced membrane thickness, as seen in the electron density profiles across the membrane. From the lateral pressure profiles, we found that the values of spontaneous curvature of the two opposing monolayers became negative when the pH was decreased. At low pH, protonation of the free fatty acid headgroups reduces their mutual repulsion and accounts for the pH dependence of all the above-mentioned properties. The altered structural characteristics may significantly affect the overall surface properties of biomembranes in cellular vesicles, lipid droplets, and plasma lipoproteins, play an important role in membrane fission and fusion, and modify interactions between membrane lipids and the proteins embedded within them.