Hydration of biologically relevant tetramethylammonium cation by neutron scattering and molecular dynamics.

Neutron scattering and molecular dynamics studies were performed on a concentrated aqueous tetramethylammonium (TMA) chloride solution to gain insight into the hydration shell structure of TMA, which is relevant for understanding its behavior in biological contexts of, e.g., properties of phospholipid membrane headgroups or interactions between DNA and histones. Specifically, neutron diffraction with isotopic substitution experiments were performed on TMA and water hydrogens to extract the specific correlation between hydrogens in TMA (HTMA) and hydrogens in water (HW). Classical molecular dynamics simulations were performed to help interpret the experimental neutron scattering data. Comparison of the hydration structure and simulated neutron signals obtained with various force field flavors (e.g. overall charge, charge distribution, polarity of the CH bonds and geometry) allowed us to gain insight into how sensitive the TMA hydration structure is to such changes and how much the neutron signal can capture them. We show that certain aspects of the hydration, such as the correlation of the hydrogen on TMA to hydrogen on water, showed little dependence on the force field. In contrast, other correlations, such as the ion-ion interactions, showed more marked changes. Strikingly, the neutron scattering signal cannot discriminate between different hydration patterns. Finally, ab initio molecular dynamics was used to examine the three-dimensional hydration structure and thus to benchmark force field simulations. Overall, while neutron scattering has been previously successfully used to improve force fields, in the particular case of TMA we show that it has only limited value to fully determine the hydration structure, with other techniques such as ab initio MD being of a significant help.

Why do polyarginines adsorb at neutral phospholipid bilayers and polylysines do not? An insight from density functional theory calculations and molecular dynamics simulations.

Adsorption of cell-penetrating peptides (CPPs) at cellular membranes is the first and necessary step for their subsequent translocation across cellular membranes into the cytosol. It has been experimentally shown that CPPs rich in arginine (Arg) amino acid penetrate across phospholipid bilayers more effectively than their lysine (Lys) rich counterparts. In this work, we aim to understand the differences in the first translocation step, adsorption of Arg9 and Lys9 peptides at fully hydrated neutral phosphatidylcholine (PC) and phosphatidylethanolamine (PE) lipid bilayers and evaluate in detail the energetics of the process using molecular dynamics (MD) simulations and free energy calculations of adsorption of the single peptide. We show that the adsorption of Arg9 is energetically feasible, with the free energy of adsorption being ∼-5.0 kcal mol-1 at PC and ∼-5.5 kcal mol-1 at PE bilayers. In contrast, adsorption of Lys9 is not observed at PC bilayers, and their adsorption at PE bilayers is very weak, being ∼-0.5 kcal mol-1. We show by energy decomposition and analysis of peptide hydration along the membrane that significantly stronger electrostatic interactions of Arg9 with lipid phosphate groups, together with the greater loss of peptide hydration (and in turn stronger hydrophobic interactions) along the membrane translocation path, are the main driving factors governing the adsorption of Arg-rich peptides at neutral lipid bilayers in contrast to Lys-rich peptides. Finally, we also compare the energetics in lipid/bilayer systems with the density functional theory (DFT) calculations of the corresponding model systems in the continuum water model and reveal the energetic differences in different environments.

FA Sliding as the Mechanism for the ANT1-Mediated Fatty Acid Anion Transport in Lipid Bilayers.

Mitochondrial adenine nucleotide translocase (ANT) exchanges ADP for ATP to maintain energy production in the cell. Its protonophoric function in the presence of long-chain fatty acids (FA) is also recognized. Our previous results imply that proton/FA transport can be best described with the FA cycling model, in which protonated FA transports the proton to the mitochondrial matrix. The mechanism by which ANT1 transports FA anions back to the intermembrane space remains unclear. Using a combined approach involving measurements of the current through the planar lipid bilayers reconstituted with ANT1, site-directed mutagenesis and molecular dynamics simulations, we show that the FA anion is first attracted by positively charged arginines or lysines on the matrix side of ANT1 before moving along the positively charged protein-lipid interface and binding to R79, where it is protonated. We show that R79 is also critical for the competitive binding of ANT1 substrates (ADP and ATP) and inhibitors (carboxyatractyloside and bongkrekic acid). The binding sites are well conserved in mitochondrial SLC25 members, suggesting a general mechanism for transporting FA anions across the inner mitochondrial membrane.

Ionic Strength and Solution Composition Dictate the Adsorption of Cell-Penetrating Peptides onto Phosphatidylcholine Membranes.

Adsorption of arginine-rich positively charged peptides onto neutral zwitterionic phosphocholine (PC) bilayers is a key step in the translocation of those potent cell-penetrating peptides into the cell interior. In the past, we have shown both theoretically and experimentally that polyarginines adsorb to the neutral PC-supported lipid bilayers in contrast to polylysines. However, comparing our results with previous studies showed that the results often do not match even at the qualitative level. The adsorption of arginine-rich peptides onto 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) may qualitatively depend on the actual experimental conditions where binding experiments have been performed. In this work, we systematically studied the adsorption of R9 and K9 peptides onto the POPC bilayer, aided by molecular dynamics (MD) simulations and fluorescence cross-correlation spectroscopy (FCCS) experiments. Using MD simulations, we tested a series of increasing peptide concentrations, in parallel with increasing Na+ and Ca2+ salt concentrations, showing that the apparent strength of adsorption of R9 decreases upon the increase of peptide or salt concentration in the system. The key result from the simulations is that the salt concentrations used experimentally can alter the picture of peptide adsorption qualitatively. Using FCCS experiments with fluorescently labeled R9 and K9, we first demonstrated that the binding of R9 to POPC is tighter by almost 2 orders of magnitude compared to that of K9. Finally, upon the addition of an excess of either Na+ or Ca2+ ions with R9, the total fluorescence correlation signal is lost, which implies the unbinding of R9 from the PC bilayer, in agreement with our predictions from MD simulations.

Molecular Dynamics Simulations of Mitochondrial Uncoupling Protein 2.

Molecular dynamics (MD) simulations of uncoupling proteins (UCP), a class of transmembrane proteins relevant for proton transport across inner mitochondrial membranes, represent a complicated task due to the lack of available structural data. In this work, we use a combination of homology modelling and subsequent microsecond molecular dynamics simulations of UCP2 in the DOPC phospholipid bilayer, starting from the structure of the mitochondrial ATP/ADP carrier (ANT) as a template. We show that this protocol leads to a structure that is impermeable to water, in contrast to MD simulations of UCP2 structures based on the experimental NMR structure. We also show that ATP binding in the UCP2 cavity is tight in the homology modelled structure of UCP2 in agreement with experimental observations. Finally, we corroborate our results with conductance measurements in model membranes, which further suggest that the UCP2 structure modeled from ANT protein possesses additional key functional elements, such as a fatty acid-binding site at the R60 region of the protein, directly related to the proton transport mechanism across inner mitochondrial membranes.

ANT1 Activation and Inhibition Patterns Support the Fatty Acid Cycling Mechanism for Proton Transport.

Adenine nucleotide translocase (ANT) is a well-known mitochondrial exchanger of ATP against ADP. In contrast, few studies have shown that ANT also mediates proton transport across the inner mitochondrial membrane. The results of these studies are controversial and lead to different hypotheses about molecular transport mechanisms. We hypothesized that the H+-transport mediated by ANT and uncoupling proteins (UCP) has a similar regulation pattern and can be explained by the fatty acid cycling concept. The reconstitution of purified recombinant ANT1 in the planar lipid bilayers allowed us to measure the membrane current after the direct application of transmembrane potential ΔΨ, which would correspond to the mitochondrial states III and IV. Experimental results reveal that ANT1 does not contribute to a basal proton leak. Instead, it mediates H+ transport only in the presence of long-chain fatty acids (FA), as already known for UCPs. It depends on FA chain length and saturation, implying that FA's transport is confined to the lipid-protein interface. Purine nucleotides with the preference for ATP and ADP inhibited H+ transport. Specific inhibitors of ATP/ADP transport, carboxyatractyloside or bongkrekic acid, also decreased proton transport. The H+ turnover number was calculated based on ANT1 concentration determined by fluorescence correlation spectroscopy and is equal to 14.6 ± 2.5 s-1. Molecular dynamic simulations revealed a large positively charged area at the protein/lipid interface that might facilitate FA anion's transport across the membrane. ANT's dual function-ADP/ATP and H+ transport in the presence of FA-may be important for the regulation of mitochondrial membrane potential and thus for potential-dependent processes in mitochondria. Moreover, the expansion of proton-transport modulating drug targets to ANT1 may improve the therapy of obesity, cancer, steatosis, cardiovascular and neurodegenerative diseases.

Self-association of a highly charged arginine-rich cell-penetrating peptide.

Small-angle X-ray scattering (SAXS) measurements reveal a striking difference in intermolecular interactions between two short highly charged peptides-deca-arginine (R10) and deca-lysine (K10). Comparison of SAXS curves at high and low salt concentration shows that R10 self-associates, while interactions between K10 chains are purely repulsive. The self-association of R10 is stronger at lower ionic strengths, indicating that the attraction between R10 molecules has an important electrostatic component. SAXS data are complemented by NMR measurements and potentials of mean force between the peptides, calculated by means of umbrella-sampling molecular dynamics (MD) simulations. All-atom MD simulations elucidate the origin of the R10-R10 attraction by providing structural information on the dimeric state. The last two C-terminal residues of R10 constitute an adhesive patch formed by stacking of the side chains of two arginine residues and by salt bridges formed between the like-charge ion pair and the C-terminal carboxyl groups. A statistical analysis of the Protein Data Bank reveals that this mode of interaction is a common feature in proteins.

Arginine "Magic": Guanidinium Like-Charge Ion Pairing from Aqueous Salts to Cell Penetrating Peptides.

It is a textbook knowledge that charges of the same polarity repel each other. For two monovalent ions in the gas phase at a close contact this repulsive interaction amounts to hundreds of kilojoules per mole. In aqueous solutions, however, this Coulomb repulsion is strongly attenuated by a factor equal to the dielectric constant of the medium. The residual repulsion, which now amounts only to units of kilojoules per mole, may be in principle offset by attractive interactions. Probably the smallest cationic pair, where a combination of dispersion and cavitation forces overwhelms the Coulomb repulsion, consists of two guanidinium ions in water. Indeed, by a combination of molecular dynamics with electronic structure calculations and electrophoretic, as well as spectroscopic, experiments, we have demonstrated that aqueous guanidinium cations form (weakly) thermodynamically stable like-charge ion pairs. The importance of pairing of guanidinium cations in aqueous solutions goes beyond a mere physical curiosity, since it has significant biochemical implications. Guanidinium chloride is known to be an efficient and flexible protein denaturant. This is due to the ability of the orientationally amphiphilic guanidinium cations to disrupt various secondary structural motifs of proteins by pairing promiscuously with both hydrophobic and hydrophilic groups, including guanidinium-containing side chains of arginines. The fact that the cationic guanidinium moiety forms the dominant part of the arginine side chain implies that the like-charge ion pairing may also play a role for interactions between peptides and proteins. Indeed, arginine-arginine pairing has been frequently found in structural protein databases. In particular, when strengthened by a presence of negatively charged glutamate, aspartate, or C-terminal carboxylic groups, this binding motif helps to stabilize peptide or protein dimers and is also found in or near active sites of several enzymes. The like-charge pairing of the guanidinium side-chain groups may also hold the key to the understanding of the arginine "magic", that is, the extraordinary ability of arginine-rich polypeptides to passively penetrate across cellular membranes. Unlike polylysines, which are also highly cationic but lack the ease in crossing membranes, polyarginines do not exhibit mutual repulsion. Instead, they accumulate at the membrane, weaken it, and might eventually cross in a concerted, "train-like" manner. This behavior of arginine-rich cell penetrating peptides can be exploited when devising smart strategies how to deliver in a targeted way molecular cargos into the cell.

Revisited Mechanism of Reaction between a Model Lysine Amino Acid Side Chain and 4-Hydroxynonenal in Different Solvent Environments.

We revisit the mechanism of reaction between a model lysine side chain and reactive aldehyde 4-hydroxynonenal in different solvents with an increasing water content. We show by model organic reactions and qualitative spectrometric analysis that a nonpolar pyrrole adduct is dominantly formed in non-aqueous solvents dichloromethane and acetonitrile. On the other hand, in aqueous acetonitrile and neat water, other polar products are also isolated, including Michael adducts, hemiacetal adducts, and pyridinium salt adducts, at the same time as the ratio of nonpolar products to polar products is decreasing. The experiments are supported by detailed quantum chemical calculations of the reaction mechanism with different computational setups showing that the pyrrole adduct is the most thermodynamically stable product compared to Michael adducts and hemiacetal adducts and also indicating that water molecules released along the reaction pathway are catalyzing reaction steps involving proton transfer. Finally, we also identify the mechanism of the pyridinium salt adduct that is formed only in aqueous solutions.

Calculation of apparent pKa values of saturated fatty acids with different lengths in DOPC phospholipid bilayers.

We performed all-atom molecular dynamics simulations and calculated free energy profiles and apparent pKa values for neutral and anionic forms of single myristic (C14:0), palmitic (C16:0) and stearic (C18:0) fatty acid embedded in a DOPC bilayer and explicit water solvent. We showed that the neutral forms of the fatty acids are stabilized inside the bilayer by hydrogen bonding of a fatty acid carboxylic group with DOPC phosphate and carbonyl groups. In contrast to the neutral form, the anionic forms of the fatty acids are shifted towards the water-membrane interface and are instead stabilized by hydrogen bonding to interfacial water. By using umbrella sampling simulations, we calculated free energies of stabilization and revealed that the free energy of stabilization inside the bilayer increases with the chain length for both the neutral and deprotonated forms. On the other hand, the free energies of flip-flop of both the neutral and anionic forms are constant upon the prolongation of the fatty acid. Based on the free energy curves, we also calculated apparent fatty acid pKa,app values in the bilayer, which are 7.0, 7.2 and 6.3 for myristic, palmitic and stearic acid and are increased by several pKa units compared to the corresponding pKa values in water. By further analysis of the calculated curves we found that spontaneous protonation of fatty acid anions takes place in the bilayer interior at ca. 1.4 nm from the bilayer center for all studied fatty acids.

Mechanism of Long-Chain Free Fatty Acid Protonation at the Membrane-Water Interface.

Long-chain free fatty acids (FFAs) play an important role in several physiological and pathological processes such as lipid fusion, adjustments of membrane permeability and fluidity, and the regulation of enzyme and protein activities. FFA-facilitated membrane proton transport (flip-flop) and FFA-dependent proton transport by membrane proteins (e.g., mitochondrial uncoupling proteins) are governed by the difference between FFA's intrinsic pKa value and the pH in the immediate membrane vicinity. Thus far, a quantitative understanding of the process has been hampered, because the pKa value shifts upon moving the FFA from the aqueous solution into the membrane. For the same FFA, pKa values between 5 and 10.5 were reported. Here, we systematically evaluated the dependence of pKa values on chain length and number of double bonds by measuring the ζ-potential of liposomes reconstituted with FFA at different pH values. The experimentally obtained intrinsic pKa values (6.25, 6.93, and 7.28 for DOPC membranes) increased with FFA chain length (C16, C18, and C20), indicating that the hydrophobic energy of transfer into the bilayer is an important pKa determinant. The observed pKa decrease in DOPC with increasing number of FFA double bonds (7.28, 6.49, 6.16, and 6.13 for C20:0, C20:1, C20:2, and C20:4, respectively) is in line with a decrease in transfer energy. Molecular dynamic simulations revealed that the ionized carboxylic group of the FFAs occupied a fixed position in the bilayer independent of chain length, underlining the importance of Born energy. We conclude that pKa is determined by the interplay between the energetic costs for 1) burying the charged moiety into the lipid bilayer and 2) transferring the hydrophobic protonated FFA into the bilayer.

Photocyclization of Tetra- and Pentapeptides Containing Adamantylphthalimide and Phenylalanines: Reaction Efficiency and Diastereoselectivity.

A series of tetrapeptides and pentapeptides was synthesized bearing a phthalimide chromophore at the N-terminus. The C-terminus of the peptides was strategically substituted with an amino acid, Phe, Phe(OMe), or Phe(OMe)2 characterized by different oxidation potentials. The photochemical reactivity of the peptides was investigated by preparative irradiation and isolation of photoproducts, as well as with laser flash photolysis. Upon photoexcitation, the peptides undergo photoinduced electron transfer (PET) and decarboxylation, followed by diastereoselective cyclization with the retention of configuration for tetrapeptides or inversion of configuration for pentapeptides. Molecular dynamics (MD) simulations and NOE experiments enabled assignment of the stereochemistry of the cyclic peptides. MD simulations of the linear peptides disclosed conformational reasons for the observed diastereoselectivity, being due to the peptide backbone spatial orientation imposed by the Phe amino acids. The photochemical efficiency for the decarboxylation and cyclization is not dependent on the peptide length, but it depends on the oxidation potential of the amino acid at the C-terminus. The results described herein are particularly important for the rational design of efficient photochemical reactions for the preparation of cyclic peptides with the desired selectivity.

Synthesis of Glycomimetics by Diastereoselective Passerini Reaction.

We describe the utilization of bis-isopropylidene-protected d-fructose-derived aldehyde in the Passerini reaction with various acids and isocyanides. A library of densely functionalized glycomimetics bearing up to 3 carbohydrate units was obtained in high yields and diastereoselectivities. The configuration of the newly formed stereocenter was determined and the diastereoselectivity was rationalized by DFT calculations.

Coarse-grained model of titrating peptides interacting with lipid bilayers.

Molecular-level computer simulations of peptide aggregation, translocation, and protonation at and in biomembranes are impeded by the large time and length scales involved. We present a computationally efficient, coarse-grained, and solvent-free model for the interaction between lipid bilayers and peptides. The model combines an accurate description of mechanical membrane properties with a new granular representation of the dielectric mismatch between lipids and the aqueous phase. All-atom force fields can be easily mapped onto the coarse-grained model, and parameters for coarse-grained monopeptides accurately extrapolate to membrane permeation free energies for the corresponding dipeptides and tripeptides. Acid-base equilibria of titratable amino acid residues are further studied using a constant-pH ensemble, capturing protonation state changes upon membrane translocation. Important differences between histidine, lysine, and arginine are observed, which are in good agreement with experimental observations.

Anomalous surface behavior of hydrated guanidinium ions due to ion pairing.

Surface affinity of aqueous guanidinium chloride (GdmCl) is compared to that of aqueous tetrapropylammonium chloride (TPACl) upon addition of sodium chloride (NaCl) or disodium sulfate (Na2SO4). The experimental results have been acquired using the surface sensitive technique X-ray photoelectron spectroscopy on a liquid jet. Molecular dynamics simulations have been used to produce radial distribution functions and surface density plots. The surface affinities of both TPA+ and Gdm+ increase upon adding NaCl to the solution. With the addition of Na2SO4, the surface affinity of TPA+ increases, while that of Gdm+ decreases. From the results of MD simulations it is seen that Gdm+ and SO42- ions form pairs. This finding can be used to explain the decreased surface affinity of Gdm+ when co-dissolved with SO42- ions. Since SO42- ions avoid the surface due to the double charge and strong water interaction, the Gdm+-SO42- ion pair resides deeper in the solutions' bulk than the Gdm+ ions. Since TPA+ does not form ion pairs with SO42-, the TPA+ ions are instead enriched at the surface.

Reaction Mechanism of Covalent Modification of Phosphatidylethanolamine Lipids by Reactive Aldehydes 4-Hydroxy-2-nonenal and 4-Oxo-2-nonenal.

4-Hydroxy-2-nonenal (HNE) and 4-oxo-2-nonenal (ONE) are biologically important reactive aldehydes formed during oxidative stress in phospholipid bilayers. They are highly reactive species due to presence of several reaction centers and can react with amino acids in peptides and proteins, as well as phosphoethanolamine (PE) lipids, thus modifying their biological activity. The aim of this work is to study in a molecular detail the reactivity of HNE and ONE toward PE lipids in a simplified system containing only lipids and reactive aldehydes in dichloromethane as an inert solvent. We use a combination of quantum chemical calculations, 1H NMR measurements, FT-IR spectroscopy, and mass spectrometry experiments and show that for both reactive aldehydes two types of chemical reactions are possible: formation of Michael adducts and Schiff bases. In the case of HNE, an initially formed Michael adduct can also undergo an additional cyclization step to a hemiacetal derivative, whereas no cyclization occurs in the case of ONE and a Michael adduct is identified. A Schiff base product initially formed when HNE is added to PE lipid can also further cyclize to a pyrrole derivative in contrast to ONE, where only a Schiff base product is isolated. The suggested reaction mechanism by quantum-chemical calculations is in a qualitative agreement with experimental yields of isolated products and is also additionally investigated by 1H NMR measurements, FT-IR spectroscopy, and mass spectrometry experiments.

The Antiferromagnetic Spin Coupling in Non-Kekulé Acenes-Impressive Polyradical Character Revealed by High-Level Multireference Methods.

Complete active space (CASSCF) and multireference (MR-CISD(Q) and MR-AQCC) calculations were performed for non-Kekulé analogues of acenes, dimethylenepolycyclobutadienes, with lengths of up to eight cyclobutadiene (CBD) units. Multireference calculations predict that the most stable energy state of the system is either triplet (if there is an odd number of CBD units) or singlet (if there is an even number of CBD units) due to antiferromagnetic spin coupling, which thus violates Hund's rule in larger molecules. We also show an impressive polyradical character in the system that increases with the size of the molecule, as witnessed by more than eleven unpaired electrons in the singlet state of the molecule with eight CBD units. Together with the small energy gap between singlet and higher multiplicity energy states even above the triplet state, this demonstrates the exceptional polyradical properties of these π-conjugated oligomeric chains.

Ab initio study of the reaction of ozone with bromide ion.

Surface level ozone destruction in polar environments may be initiated by oxidation of bromide ions by ozone, ultimately leading to Br2 production. Ab initio calculations are used to support the development of atmospheric chemistry models, but errors can occur in study of the bromide-ozone reaction due to inappropriate treatment of the many-electron species and the ionic nature of the reaction. In this work, a high level ab initio study is used to take into account the electronic correlation and the polarization effects. Our results show three possible pathways for the reaction. In particular, we find that this process, though endothermic on the singlet spin state surface, can be energetically feasible on the triplet surface. The triplet surface can be reached through photoexcitation of ozone or by the spin crossing of the potential energy surface. Because this process is known to occur in the dark, it may be that it occurs after intersystem crossing to a triplet surface. This paper also provides a starting point calibration for any future ab initio calculation studies of the bromide-ozone reaction, from the gas to the condensed phase.

Does fluoride disrupt hydrogen bond network in cationic lipid bilayer? Time-dependent fluorescence shift of Laurdan and molecular dynamics simulations.

Time-dependent fluorescence shift (TDFS) of Laurdan embedded in phospholipid bilayers reports on hydration and mobility of the phospholipid acylgroups. Exchange of H2O with D2O prolongs the lifetime of lipid-water and lipid-water-lipid interactions, which is reflected in a significantly slower TDFS kinetics. Combining TDFS measurements in H2O and D2O hydrated bilayers with atomistic molecular dynamics (MD) simulations provides a unique tool for characterization of the hydrogen bonding at the acylgroup level of lipid bilayers. In this work, we use this approach to study the influence of fluoride anions on the properties of cationic bilayers composed of trimethylammonium-propane (DOTAP). The results obtained for DOTAP are confronted with those for neutral phosphatidylcholine (DOPC) bilayers. Both in DOTAP and DOPC H2O/D2O exchange prolongs hydrogen-bonding lifetime and does not disturb bilayer structure. These results are confirmed by MD simulations. TDFS experiments show, however, that for DOTAP this effect is cancelled in the presence of fluoride ions. We interpret these results as evidence that strongly hydrated fluoride is able to steal water molecules that bridge lipid carbonyls. Consequently, when attracted to DOTAP bilayer, fluoride disrupts the local hydrogen-bonding network, and the differences in TDFS kinetics between H2O and D2O hydrated bilayers are no longer observed. A distinct behavior of fluoride is also evidenced by MD simulations, which show different lipid-ion binding for Cl(-) and F(-).

Protonation of palmitic acid embedded in DPPC lipid bilayers obscures detection of ripple phase by FTIR spectroscopy.

The transformation of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid bilayers from the gel (Lß') to the fluid (Lα) phase involves an intermediate ripple (Pß') phase forming a few degrees below the main transition temperature (Tm). While the exact cause of bilayer rippling is still debated, the presence of amphiphilic molecules, pH, and lipid bilayer architecture are all known to influence (pre)transition behavior. In particular, fatty acid chains interact with hydrophobic lipid tails, while the carboxylic groups simultaneously participate in proton transfer with interfacial water in the polar lipid region which is controlled by the pH of the surrounding aqueous medium. The molecular-level variations in the DPPC ripple phase in the presence of 2% palmitic acid (PA) were studied at pH levels 4.0, 7.3, and 9.1, where PA is fully protonated, partially protonated, or fully deprotonated. Bilayer thermotropic behavior was investigated by differential scanning calorimetry (DSC) and Fourier-transform infrared (FTIR) spectroscopy which agreed in their characterization of (pre)transition at pH of 9.1, but not at pH 4.0 and especially not at 7.3. Owing to the different insertion depths of protonated and deprotonated PA, along with the ability of protonated PA to undergo flip-flop in the bilayer, these two forms of PA show a different hydration pattern in the interfacial water layer. Finally, these results demonstrated the hitherto undiscovered potential of FTIR spectroscopy in the detection of the events occurring at the surface of lipid bilayers that obscure the low-cooperativity phase transition explored in this work.