Structures of gramicidins A, B, and C incorporated into sodium dodecyl sulfate micelles

Intrinsic Lipid Curvature and Bilayer Elasticity as Regulators of Channel Function: A Comparative Single-Molecule Study

Perturbations in bilayer material properties (thickness, lipid intrinsic curvature and elastic moduli) modulate the free energy difference between different membrane protein conformations, thereby leading to changes in the conformational preferences of bilayer-spanning proteins. To further explore the relative importance of curvature and elasticity in determining the changes in bilayer properties that underlie the modulation of channel function, we investigated how the micelle-forming amphiphiles Triton X-100, reduced Triton X-100 and the HII lipid phase promoter capsaicin modulate the function of alamethicin and gramicidin channels. Whether the amphiphile-induced changes in intrinsic curvature were negative or positive, amphiphile addition increased gramicidin channel appearance rates and lifetimes and stabilized the higher conductance states in alamethicin channels. When the intrinsic curvature was modulated by altering phospholipid head group interactions, however, maneuvers that promote a negative-going curvature stabilized the higher conductance states in alamethicin channels but destabilized gramicidin channels. Using gramicidin channels of different lengths to probe for changes in bilayer elasticity, we found that amphiphile adsorption increases bilayer elasticity, whereas altering head group interactions does not. We draw the following conclusions: first, confirming previous studies, both alamethicin and gramicidin channels are modulated by changes in lipid bilayer material properties, the changes occurring in parallel yet differing dependent on the property that is being changed; second, isolated, negative-going changes in curvature stabilize the higher current levels in alamethicin channels and destabilize gramicidin channels; third, increases in bilayer elasticity stabilize the higher current levels in alamethicin channels and stabilize gramicidin channels; and fourth, the energetic consequences of changes in elasticity tend to dominate over changes in curvature.

Screening for bilayer-active and likely cytotoxic molecules reveals bilayer-mediated regulation of cell function

A perennial problem encountered when using small molecules (drugs) to manipulate cell or protein function is to assess whether observed changes in function result from specific interactions with a desired target or from less specific off-target mechanisms. This is important in laboratory research as well as in drug development, where the goal is to identify molecules that are unlikely to be successful therapeutics early in the process, thereby avoiding costly mistakes. We pursued this challenge from the perspective that many bioactive molecules (drugs) are amphiphiles that alter lipid bilayer elastic properties, which may cause indiscriminate changes in membrane protein (and cell) function and, in turn, cytotoxicity. Such drug-induced changes in bilayer properties can be quantified as changes in the monomer↔dimer equilibrium for bilayer-spanning gramicidin channels. Using this approach, we tested whether molecules in the Pathogen Box (a library of 400 drugs and drug-like molecules with confirmed activity against tropical diseases released by Medicines for Malaria Venture to encourage the development of therapies for neglected tropical diseases) are bilayer modifiers. 32% of the molecules in the Pathogen Box were bilayer modifiers, defined as molecules that at 10 µM shifted the monomer↔dimer equilibrium toward the conducting dimers by at least 50%. Correlation analysis of the molecules' reported HepG2 cell cytotoxicity to bilayer-modifying potency, quantified as the shift in the gramicidin monomer↔dimer equilibrium, revealed that molecules producing <25% change in the equilibrium had significantly lower probability of being cytotoxic than molecules producing >50% change. Neither cytotoxicity nor bilayer-modifying potency (quantified as the shift in the gramicidin monomer↔dimer equilibrium) was well predicted by conventional physico-chemical descriptors (hydrophobicity, polar surface area, etc.). We conclude that drug-induced changes in lipid bilayer properties are robust predictors of the likelihood of membrane-mediated off-target effects, including cytotoxicity.

Delivery of antimicrobial peptides to model membranes by cubosome nanocarriers

Antimicrobial peptides (AMPs), which typically disrupt the bacterial wall prompting leakage or lysis of the cell, form a growing contingent in the arsenal against antibiotic resistant bacteria. The effectiveness of AMPs is, however, hampered by their low solubility, general chemical and physical instability, and short half-life in vivo. Lipid nanocarriers such as cubosomes are effective at encapsulating and protecting proteins while simultaneously showing promise in delivery applications. Here, the efficacy of cubosome mediated delivery of AMPs is evaluated by the in-situ surface characterization of model membranes with varying composition. The cubosomes were observed to initially fuse with the membranes, with subsequent membrane disruption observed after approximately 20 - 60 min. The time for the disruption was sensitive to the charge of the cubosome as well as the composition of the bilayer. More physiologically relevant bilayers including lipids with phospho-(1'-rac-glycerol) (PG) or phosphoethanolamine (PE) headgroups were more vulnerable than those of neat phosphocholine (PC). Notably, disruption to the bilayer occurred an order of magnitude faster for encapsulated AMP compared to free AMP.

Polarization Effects in Water-Mediated Selective Cation Transport across a Narrow Transmembrane Channel

Despite the progress in modeling complex molecular systems of ever-increasing complexity, a quantitatively accurate computational treatment of ion permeation through narrow membrane channels remains challenging. An important factor to reach this goal is induced electronic polarization, which is likely to impact the permeation rate of small ions through narrow molecular pores. In this work, we extended the recently developed polarizable force field based on the classical Drude oscillators to assess the role of induced polarization effects on the energetics of sodium and potassium ion transport across the gramicidin A (gA) ion channel. The inclusion of induced polarization lowers barriers present in 1D potential of mean force (PMF) for cation permeation by ∼50% compared to those obtained with the additive force field. Conductance properties calculated with 1D PMFs from Drude simulations are in better agreement with experimental results. Polarization of single-file water molecules and protein atoms forming the narrow pore has a direct impact on the free-energy barriers and cation-specific solid-state NMR chemical shifts. Sensitivity analysis indicates that small changes to water-channel interactions can alter the free energy barrier for ion permeation. These results, illustrating polarization effects present in the complex electrostatic environment of the gA channel, have broad implications for revising proposed mechanisms of ion permeation and selectivity in a variety of ion channels.

Atomistic Characterization of Gramicidin Channel Formation

We investigated gramicidin A (gA) subunit dimerization in lipid bilayers using microsecond-long replica-exchange umbrella sampling simulations, millisecond-long unbiased molecular dynamics simulations, and machine learning. Our simulations led to a dimer structure that is indistinguishable from the experimentally determined gA channel structures, with the two gA subunits joined by six hydrogen bonds (6HB). The simulations also uncovered two additional dimer structures, with different gA-gA stacking orientations that were stabilized by four or two hydrogen bonds (4HB or 2HB). When examining the temporal evolution of the dimerization, we found that two bilayer-inserted gA subunits can form the 6HB dimer directly, with no discernible intermediate states, as well as through paths that involve the 2HB and 4HB dimers.

Deciphering the Structure of the Gramicidin A Channel in the Presence of AOT Reverse Micelles in Pentane Using Molecular Dynamics Simulations

Structural studies of proteins and, in particular, integral membrane proteins (IMPs) using solution NMR spectroscopy approaches are challenging due to not only their inherent structural complexities but also the fact that they need to be solubilized in biomimetic environments (such as micelles), which enhances the slow molecular reorientation. To deal with these difficulties and increase the effective rate of molecular reorientation, the encapsulation of IMPs in the aqueous core of the reverse micelle (RM) dissolved in a low-viscosity solvent has been proven to be a viable approach. However, the effect of the reverse micelle (RM) environment on the IMP structure and function is little known. To gain insight into these aspects, this article presents a series of atomistic unconstrained molecular dynamics (MD) of a model ion channel (gramicidin A, gA) with RMs formed with anionic surfactant diacyl chain bis(2-ethylhexyl) sodium succinate (AOT) in pentane at a water-to-surfactant molar ratio (W₀) of 6. The simulations were carried out with different protocols and starting conditions for a total of 2.4 μs and were compared with other MDs used with the gA channel inserted in models of the SDS micelle or the DMPC membrane. We show here that in the presence of AOT RMs the gA dimer did not look like the "dumbbell-like" model anticipated by experiments, where the C-terminal parts of the gA are capped with two RMs and the rest of the dimer is protected from the oil solvent by the AOT acyl chains. In contrast, the MD simulations reveal that the AOT, Na⁺, and water formed two well-defined and elongated RMs attached to the C-terminal ends of the gA dimer, while the rest is in direct contact with the pentane. The initial β^6.3 secondary structure of the gA is well conserved and filled with 6-9 waters, as in SDS micelles or the DMPC membrane. Finally, the water movement inside the gA is strongly affected by the presence of RMs at each extremity, and no passage of water molecules through the gA channel is observed even after a long simulation period, whereas the opposite was found for gA in SDS and DMPC.

Anionic nanoparticle-induced perturbation to phospholipid membranes affects ion channel function

Understanding the mechanisms of nanoparticle interaction with cell membranes is essential for designing materials for applications such as bioimaging and drug delivery, as well as for assessing engineered nanomaterial safety. Much attention has focused on nanoparticles that bind strongly to biological membranes or induce membrane damage, leading to adverse impacts on cells. More subtle effects on membrane function mediated via changes in biophysical properties of the phospholipid bilayer have received little study. Here, we combine electrophysiology measurements, infrared spectroscopy, and molecular dynamics simulations to obtain insight into a mode of nanoparticle-mediated modulation of membrane protein function that was previously only hinted at in prior work. Electrophysiology measurements on gramicidin A (gA) ion channels embedded in planar suspended lipid bilayers demonstrate that anionic gold nanoparticles (AuNPs) reduce channel activity and extend channel lifetimes without disrupting membrane integrity, in a manner consistent with changes in membrane mechanical properties. Vibrational spectroscopy indicates that AuNP interaction with the bilayer does not perturb the conformation of membrane-embedded gA. Molecular dynamics simulations reinforce the experimental findings, showing that anionic AuNPs do not directly interact with embedded gA channels but perturb the local properties of lipid bilayers. Our results are most consistent with a mechanism in which anionic AuNPs disrupt ion channel function in an indirect manner by altering the mechanical properties of the surrounding bilayer. Alteration of membrane mechanical properties represents a potentially important mechanism by which nanoparticles induce biological effects, as the function of many embedded membrane proteins depends on phospholipid bilayer biophysical properties.

Assessing the Perturbing Effects of Drugs on Lipid Bilayers Using Gramicidin Channel-Based In Silico and In Vitro Assays

Partitioning of bioactive molecules, including drugs, into cell membranes may produce indiscriminate changes in membrane protein function. As a guide to safe drug development, it therefore becomes important to be able to predict the bilayer-perturbing potency of hydrophobic/amphiphilic drugs candidates. Toward this end, we exploited gramicidin channels as molecular force probes and developed in silico and in vitro assays to measure drugs' bilayer-modifying potency. We examined eight drug-like molecules that were found to enhance or suppress gramicidin channel function in a thick 1,2-dierucoyl-sn-glycero-3-phosphocholine (DC22:1PC) but not in thin 1,2-dioleoyl-sn-glycero-3-phosphocholine (DC18:1PC) lipid bilayer. The mechanism underlying this difference was attributable to the changes in gramicidin dimerization free energy by drug-induced perturbations of lipid bilayer physical properties and bilayer-gramicidin interactions. The combined in silico and in vitro approaches, which allow for predicting the perturbing effects of drug candidates on membrane protein function, have implications for preclinical drug safety assessment.

Functional stability of water wire-carbonyl interactions in an ion channel

Despite the well-characterized structural symmetry of the dimeric transmembrane antibiotic gramicidin A, we show that the symmetry is broken by selective hydrogen bonding between eight waters comprising a transmembrane water wire and a specific subset of the 26 pore-lining carbonyl oxygens of the gramicidin A channel. The ¹⁷O NMR spectroscopic resolution of the carbonyl resonances from the two subunits required the use of a world record high field magnet (35.2 T; 1,500 MHz for ¹H). Uniquely, this result documented the millisecond timescale stability of the water wire orientation within the gramicidin A pore that had been reported to have only subnanosecond stability. These ¹⁷O spectroscopic results portend wide applications in molecular biophysics and beyond.

Water wires are critical for the functioning of many membrane proteins, as in channels that conduct water, protons, and other ions. Here, in liquid crystalline lipid bilayers under symmetric environmental conditions, the selective hydrogen bonding interactions between eight waters comprising a water wire and a subset of 26 carbonyl oxygens lining the antiparallel dimeric gramicidin A channel are characterized by ¹⁷O NMR spectroscopy at 35.2 T (or 1,500 MHz for ¹H) and computational studies. While backbone ¹⁵N spectra clearly indicate structural symmetry between the two subunits, single site ¹⁷O labels of the pore-lining carbonyls report two resonances, implying a break in dimer symmetry caused by the selective interactions with the water wire. The ¹⁷O shifts document selective water hydrogen bonding with carbonyl oxygens that are stable on the millisecond timescale. Such interactions are supported by density functional theory calculations on snapshots taken from molecular dynamics simulations. Water hydrogen bonding in the pore is restricted to just three simultaneous interactions, unlike bulk water environs. The stability of the water wire orientation and its electric dipole leads to opposite charge-dipole interactions for K⁺ ions bound at the two ends of the pore, thereby providing a simple explanation for an ∼20-fold difference in K⁺ affinity between two binding sites that are ∼24 Å apart. The ¹⁷O NMR spectroscopy reported here represents a breakthrough in high field NMR technology that will have applications throughout molecular biophysics, because of the acute sensitivity of the ¹⁷O nucleus to its chemical environment.

Further Optimization and Validation of the Classical Drude Polarizable Protein Force Field

The CHARMM Drude-2013 polarizable force field (FF) was developed to include the explicit treatment of induced electronic polarizability, resulting in a more accurate description of the electrostatic interactions in molecular dynamics (MD) simulations. While the Drude-2013 protein FF has shown success in improving the folding properties of α-helical peptides and to reproduce experimental observables in simulations up to 1 μs, some limitations were noted regarding the stability of β-sheet structures in simulations longer than 100 ns as well as larger deviations from crystal structures in simulations of a number of proteins compared to the additive CHARMM36 protein FF. The origin of the instability has been identified and appears to be primarily due to overestimated atomic polarizabilities and induced dipole-dipole interactions on the Cβ, Cγ, and Cδ side chain atoms. To resolve this and other issues, a number of aspects of the model were revisited, resulting in Drude-2019 protein FF. Backbone parameters were optimized targeting the conformational properties of the (Ala)5 peptide in solution along with gas phase properties of the alanine dipeptide. Dipeptides that contain N-acetylated and N'-methylamidated termini, excluding Gly, Pro, and Ala, were used as models to optimize the atomic polarizabilities and Thole screening factors on selected Cβ, Cγ, and Cδ carbons by targeting quantum mechanical (QM) dipole moments and molecular polarizabilities. In addition, to obtain better conformational properties, side chain χ1 and χ2 dihedral parameters were optimized targeting QM data for the respective side chain dipeptide conformations as well as Protein Data Bank survey data based on the χ1, χ2 sampling from Hamiltonian replica-exchange MD simulations of (Ala)4-X-(Ala)4 in solution, where X is the amino acid of interest. Further improvements include optimizing nonbonded interactions between charged residues to reproduce QM interaction energies of the charged-protein model compounds and experimental osmotic pressures. Validation of the optimized Drude protein FF includes MD simulations of a collection of peptides and proteins including β-sheet structures, as well as transmembrane ion channels. Results showed that the updated Drude-2019 protein FF yields smaller overall root-mean-square differences of proteins as compared to the additive CHARMM36m and Drude-2013 FFs as well as similar or improved agreement with experimental NMR properties, allowing for long time scale simulation studies of proteins and more complex biomolecular systems in conjunction with the remainder of the Drude polarizable FF.

Targeting the Multiple Physiologic Roles of VDAC With Steroids and Hydrophobic Drugs

There is accumulating evidence that endogenous steroids and non-polar drugs are involved in the regulation of mitochondrial physiology. Many of these hydrophobic compounds interact with the Voltage Dependent Anion Channel (VDAC). This major metabolite channel in the mitochondrial outer membrane (MOM) regulates the exchange of ions and water-soluble metabolites, such as ATP and ADP, across the MOM, thus governing mitochondrial respiration. Proteomics and biochemical approaches together with molecular dynamics simulations have identified an impressively large number of non-polar compounds, including endogenous, able to bind to VDAC. These findings have sparked speculation that both natural steroids and synthetic hydrophobic drugs regulate mitochondrial physiology by directly affecting VDAC ion channel properties and modulating its metabolite permeability. Here we evaluate recent studies investigating the effect of identified VDAC-binding natural steroids and non-polar drugs on VDAC channel functioning. We argue that while many compounds are found to bind to the VDAC protein, they do not necessarily affect its channel functions in vitro. However, they may modify other aspects of VDAC physiology such as interaction with its cytosolic partner proteins or complex formation with other mitochondrial membrane proteins, thus altering mitochondrial function.

Molecular Mechanism for Gramicidin Dimerization and Dissociation in Bilayers of Different Thickness

Membrane protein functions can be altered by subtle changes in the host lipid bilayer physical properties. Gramicidin channels have emerged as a powerful system for elucidating the underlying mechanisms of membrane protein function regulation through changes in bilayer properties, which are reflected in the thermodynamic equilibrium distribution between nonconducting gramicidin monomers and conducting bilayer-spanning dimers. To improve our understanding of how subtle changes in bilayer thickness alter the gramicidin monomer and dimer distributions, we performed extensive atomistic molecular dynamics simulations and fluorescence-quenching experiments on gramicidin A (gA). The free-energy calculations predicted a nonlinear coupling between the bilayer thickness and channel formation. The energetic barrier inhibiting gA channel formation was sharply increased in the thickest bilayer (1,2-dierucoyl-sn-glycero-3-phosphocholine). This prediction was corroborated by experimental results on gramicidin channel activity in bilayers of different thickness. To further explore the mechanism of channel formation, we performed extensive unbiased molecular dynamics simulations, which allowed us to observe spontaneous gA dimer formation in lipid bilayers. The simulations revealed structural rearrangements in the gA subunits and changes in lipid packing, as well as water reorganization, that occur during the dimerization process. Together, the simulations and experiments provide new, to our knowledge, insights into the process and mechanism of gramicidin channel formation, as a prototypical example of the bilayer regulation of membrane protein function.

Quantitative Characterization of Protein-Lipid Interactions by Free Energy Simulation between Binary Bilayers

Using a recently developed binary bilayer system (BBS) consisting of two patches of laterally contacting bilayers, umbrella sampling molecular dynamics (MD) simulations were performed for quantitative characterization of protein-lipid interactions. The BBS is composed of 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) with an embedded model membrane protein, a gramicidin A (gA) channel. The calculated free energy difference for the transfer of a gA channel from DLPC (hydrophobic thickness ≈ 21.5 Å) to DMPC (hydrophobic thickness ≈ 25.5 Å) bilayers, ΔG(DLPC → DMPC), is -2.2 ± 0.7 kcal/mol. This value appears at odds with the traditional view that the hydrophobic length of the gA channel is ∼22 Å. To understand this discrepancy, we first note that recent MD simulations by different groups have shown that lipid bilayer thickness profiles in the vicinity of a gA channel differ qualitatively from the deformation profile predicted from continuum elastic bilayer models. Our MD simulations at low and high gA:lipid molar ratios and different membrane compositions indicate that the gA channel's effective hydrophobic length is ∼26 Å. Using this effective hydrophobic length, ΔG(DLPC → DMPC) determined here is in excellent agreement with predictions based on continuum elastic models (-3.0 to -2.2 kcal/mol) where the bilayer deformation energy is approximated as a harmonic function of the mismatch between the channel's effective hydrophobic length and the hydrophobic thickness of the bilayer. The free energy profile for gA in the BBS includes a barrier at the interface between the two bilayers which can be attributed to the line tension at the interface between two bilayers with different hydrophobic thicknesses. This observation implies that translation of a peptide between two different regions of a cell membrane (such as between the liquid ordered and disordered phases) may include effects of a barrier at the interface in addition to the relative free energies of the species far from the interface. The BBS allows for direct transfer free energy calculations between bilayers without a need of a reference medium, such as bulk water, and thus provides an efficient simulation protocol for the quantitative characterization of protein-lipid interactions at all-atom resolution.

Structural Insight into a Novel Formyltransferase and Evolution to a Nonribosomal Peptide Synthetase Tailoring Domain

Nonribosomal peptide synthetases (NRPSs) increase the chemical diversity of their products by acquiring tailoring domains. Linear gramicidin synthetase starts with a tailoring formylation (F) domain, which likely originated from a sugar formyltransferase (FT) gene. Here, we present studies on an Anoxybacillus kamchatkensis sugar FT representative of the prehorizontal gene transfer FT. Gene cluster analysis reveals that this FT acts on a UDP-sugar in a novel pathway for synthesis of a 7-formamido derivative of CMP-pseudaminic acid. We recapitulate the pathway up to and including the formylation step in vitro, experimentally demonstrating the role of the FT. We also present X-ray crystal structures of the FT alone and with ligands, which unveil contrasts with other structurally characterized sugar FTs and show close structural similarity with the F domain. The structures reveal insights into the adaptations that were needed to co-opt and evolve a sugar FT into a functional and useful NRPS domain.

Identification of Effective Dimeric Gramicidin-D Peptide as Antimicrobial Therapeutics over Drug Resistance: In-Silico Approach

Discovering and developing the antimicrobial peptides are recently focused on pharmaceutical firm, since they serve as complementary to antibiotics in prevailing over drug resistance by eliciting the disruption of microbial membrane. Still, there are lots of challenges to bring up the structurally stable and functionally efficient antimicrobial peptides. It is well known that gramicidin D is the prominent antimicrobial peptide that exists as g-AB, g-BC, and g-AC. This study analyzes the structural stability and the functional activity of hetero-dimeric double-stranded gramicidin-D peptides, thereby demonstrating its potent antimicrobial activity against antibiotic-resistant micro-organisms. To investigate the structural stability and functionality of gramicidin D, we performed static and dynamic analysis. Initially, we observed a maximum number of intermolecular interactions and membrane penetration in g-AB as compared to g-BC and g-AC. To substantiate further, the geometrical and thermodynamic parameters revealed the retention of maximum stability in g-AB than g-AC and g-BC. Thus, the conformational free energy and the binding free energy showed the variation among gramicidin-D peptides for the prediction of increased stability and functionality. In conclusion, g-AB peptide has definitely demonstrated adequate structural stability and functionality and this work will need to be considered in peptide-based drug discovery.

The effect of gramicidin inclusions on the local order of membrane components

We study the local effect of the antimicrobial peptide Gramicidin A on bilayers composed of lipids or surfactants using nuclear magnetic resonance spectroscopy and wide-angle X-ray scattering, techniques that probe the orientational and positional order of the alkyl chains, respectively. The two types of order vary with temperature and peptide concentration in complex ways which depend on the membrane composition, highlighting the subtlety of the interaction between inclusions and the host bilayer. The amplitude of the variation is relatively low, indicating that the macroscopic constants used to describe the elasticity of the bilayer are unlikely to change with the addition of peptide.

Exchange of Gramicidin between Lipid Bilayers: Implications for the Mechanism of Channel Formation

The canonical mechanism of gramicidin (gA) channel formation is transmembrane dimerization of nonconducting subunits that reside in opposite bilayer leaflets. The channels do not open and close; they appear and disappear due to subunit association and dissociation. Many different types of experiments support this monomer ↔ dimer mechanism. Recently, however, this mechanism was challenged, based on experiments with lipid vesicle-incorporated gA under conditions where vesicle fusion could be controlled. In these experiments, sustained channel activity was observed long after fusion had been terminated, which led to the proposal that gA single-channel current transitions result from closed-open transitions in long-lived bilayer-spanning dimers. This proposal is at odds with 40 years of experiments, but involves the key assumption that gA monomers do not exchange between bilayers. We tested the possibility of peptide exchange between bilayers using three different types of experiments. First, we demonstrated the exchange of gA between 1,2-dierucoyl-sn-glycero-3-phosphocholine (DC_22:1PC) or 1,2-dioleoyl-sn-glycero-3-phosphocholine (DC_18:1PC) lipid vesicles using a fluorescence assay for gA channel activity. Second, we added gA-free DC_22:1PC vesicles to both sides of planar DC_18:1PC bilayers preincubated with gA, which reduced channel activity up to 10-fold. Third, we added gA-containing DC_22:1PC vesicles to one or both sides of DC_18:1PC planar bilayers, which produced much higher channel activity when the gA-containing vesicles were added to both sides of the bilayer, as compared to one side only. All three types of experiments show that gA subunits can exchange between lipid bilayers. The exchange of subunits between bilayers thus is firmly established, which becomes a crucial consideration with respect to the mechanism of channel formation.

Classical Molecular Dynamics with Mobile Protons

An important limitation of standard classical molecular dynamics simulations is the inability to make or break chemical bonds. This restricts severely our ability to study processes that involve even the simplest of chemical reactions, the transfer of a proton. Existing approaches for allowing proton transfer in the context of classical mechanics are rather cumbersome and have not achieved widespread use and routine status. Here we reconsider the combination of molecular dynamics with periodic stochastic proton hops. To ensure computational efficiency, we propose a non-Boltzmann acceptance criterion that is heuristically adjusted to maintain the correct or desirable thermodynamic equilibria between different protonation states and proton transfer rates. Parameters are proposed for hydronium, Asp, Glu, and His. The algorithm is implemented in the program CHARMM and tested on proton diffusion in bulk water and carbon nanotubes and on proton conductance in the gramicidin A channel. Using hopping parameters determined from proton diffusion in bulk water, the model reproduces the enhanced proton diffusivity in carbon nanotubes and gives a reasonable estimate of the proton conductance in gramicidin A.

Delivering Transmembrane Peptide Complexes to the Gas Phase Using Nanodiscs and Electrospray Ionization

The gas-phase conformations of dimers of the channel-forming membrane peptide gramicidin A (GA), produced from isobutanol or aqueous solutions of GA-containing nanodiscs (NDs), are investigated using electrospray ionization-ion mobility separation-mass spectrometry (ESI-IMS-MS) and molecular dynamics (MD) simulations. The IMS arrival times measured for (2GA + 2Na)²⁺ ions from isobutanol reveal three different conformations, with collision cross-sections (Ω) of 683 Å² (conformation 1, C1), 708 Å² (C2), and 737 Å² (C3). The addition of NH₄CH₃CO₂ produced (2GA + 2Na)²⁺ and (2GA + H + Na)²⁺ ions, with Ω similar to those of C1, C2, and C3, as well as (2GA + 2H)²⁺, (2GA + 2NH₄)²⁺, and (2GA + H + NH₄)²⁺ ions, which adopt a single conformation with a Ω similar to that of C2. These results suggest that the nature of the charging agents, imparted by the ESI process, can influence dimer conformation in the gas phase. Notably, the POPC NDs produced exclusively (2GA + 2NH₄)²⁺ dimer ions; the DMPC NDs produced both (2GA + 2H)²⁺ and (2GA + 2NH₄)²⁺ dimer ions. While the Ω of (2GA + 2H)²⁺ is similar to that of C2, the (2GA + 2NH₄)²⁺ ions from NDs adopt a more compact structure, with a Ω of 656 Å². It is proposed that this compact structure corresponds to the ion conducting single stranded head-to-head helical GA dimer. These findings highlight the potential of NDs, combined with ESI, for transferring transmembrane peptide complexes directly from lipid bilayers to the gas phase. Graphical Abstract ᅟ.

Time-resolved measurements of an ion channel conformational change driven by a membrane phase transition

Using temperature-jump infrared spectroscopy, we are able to trigger a gel-to-fluid phase transition in lipid vesicles and monitor in real time how a membrane protein responds to structural changes in the membrane. The melting of lipid domains in 1,2-dimyristoyl-sn-glycero-3-phosphocholine vesicles is observed to occur in as fast as 50 ns, with a temperature dependence characteristic of critical slowing. Gramicidin D (gD) added to the membrane responds primarily to the change in thickness of the membrane on a timescale coincident with the membrane melting. Using structure-based spectral modeling, we assign the conformational changes to compression and rotation of a partially dissociated gD dimer. Free energy calculations indicate that the high rate is a result of near-barrierless diffusion on a protein energy landscape that is radically reshaped by membrane thinning. The structural changes associated with the phase transition are similar to the fluctuation modes of fluid phase membranes, highlighting the importance of understanding the dynamic nature of the membrane environment around proteins.

Characterizing Residue-Bilayer Interactions Using Gramicidin A as a Scaffold and Tryptophan Substitutions as Probes

Previous experiments have shown that the lifetime of a gramicidin A dimer channel (which forms from two nonconducting monomers) in a lipid bilayer is modulated by mutations of the tryptophan (Trp) residues at the bilayer-water interface. We explore this further using extensive molecular dynamics simulations of various gA dimer and monomer mutants at the Trp positions in phosphatidylcholine bilayers with different tail lengths. gA interactions with the surrounding bilayer are strongly modulated by mutating these Trp residues. There are three principal effects: eliminating residue hydrogen bonding ability (i.e., reducing the channel-monolayer coupling strength) reduces the extent of the bilayer deformation caused by the assembled dimeric channel; a residue's size and geometry affects its orientation, leading to different hydrogen bonding partners; and increasing a residue's hydrophobicity increases the depth of gA monomer insertion relative to the bilayer center, thereby increasing the lipid bending frustration.

New Continuum Approaches for Determining Protein-Induced Membrane Deformations

The influence of the membrane on transmembrane proteins is central to a number of biological phenomena, notably the gating of stretch activated ion channels. Conversely, membrane proteins can influence the bilayer, leading to the stabilization of particular membrane shapes, topological changes that occur during vesicle fission and fusion, and shape-dependent protein aggregation. Continuum elastic models of the membrane have been widely used to study protein-membrane interactions. These mathematical approaches produce physically interpretable membrane shapes, energy estimates for the cost of deformation, and a snapshot of the equilibrium configuration. Moreover, elastic models are much less computationally demanding than fully atomistic and coarse-grained simulation methodologies; however, it has been argued that continuum models cannot reproduce the distortions observed in fully atomistic molecular dynamics simulations. We suggest that this failure can be overcome by using chemically and geometrically accurate representations of the protein. Here, we present a fast and reliable hybrid continuum-atomistic model that couples the protein to the membrane. We show that the model is in excellent agreement with fully atomistic simulations of the ion channel gramicidin embedded in a POPC membrane. Our continuum calculations not only reproduce the membrane distortions produced by the channel but also accurately determine the channel's orientation. Finally, we use our method to investigate the role of membrane bending around the charged voltage sensors of the transient receptor potential cation channel TRPV1. We find that membrane deformation significantly stabilizes the energy of insertion of TRPV1 by exposing charged residues on the S4 segment to solution.

Gramicidin A Channel Formation Induces Local Lipid Redistribution II: A 3D Continuum Elastic Model

To change conformation, a protein must deform the surrounding bilayer. In this work, a three-dimensional continuum elastic model for gramicidin A in a lipid bilayer is shown to describe the sensitivity to thickness, curvature stress, and the mechanical properties of the lipid bilayer. A method is demonstrated to extract the gramicidin-lipid boundary condition from all-atom simulations that can be used in the three-dimensional continuum model. The boundary condition affects the deformation dramatically, potentially much more than typical variations in the material stiffness do as lipid composition is changed. Moreover, it directly controls the sensitivity to curvature stress. The curvature stress and hydrophobic surfaces of the all-atom and continuum models are found to be in excellent agreement. The continuum model is applied to estimate the enrichment of hydrophobically matched lipids near the channel in a mixture, and the results agree with single-channel experiments and extended molecular dynamics simulations from the companion article by Beaven et al. in this issue of Biophysical Journal.

Gramicidin A Channel Formation Induces Local Lipid Redistribution I: Experiment and Simulation

Integral membrane protein function can be modulated by the host bilayer. Because biological membranes are diverse and nonuniform, we explore the consequences of lipid diversity using gramicidin A channels embedded in phosphatidylcholine (PC) bilayers composed of equimolar mixtures of di-oleoyl-PC and di-erucoyl-PC (dC_18:1+dC_22:1, respectively), di-palmitoleoyl-PC and di-nervonoyl-PC (dC_16:1+dC_24:1, respectively), and di-eicosenoyl-PC (pure dC_20:1), all of which have the same average bilayer chain length. Single-channel lifetime experiments, molecular dynamics simulations, and a simple lipid compression model are used in tandem to gain insight into lipid redistribution around the channel, which partially alleviates the bilayer deformation energy associated with channel formation. The average single-channel lifetimes in the two-component bilayers (95 ± 10 ms for dC_18:1+dC_22:1 and 195 ± 20 ms for dC_16:1+dC_24:1) were increased relative to the single-component dC_20:1 control bilayer (65 ± 10 ms), implying lipid redistribution. Using a theoretical treatment of thickness-dependent changes in channel lifetimes, the effective local enrichment of lipids around the channel was estimated to be 58 ± 4% dC_18:1 and 66 ± 2% dC_16:1 in the dC_18:1+dC_22:1 and dC_16:1+dC_24:1 bilayers, respectively. 3.5-μs molecular dynamics simulations show 66 ± 2% dC_16:1 in the first lipid shell around the channel in the dC_16:1+dC_24:1 bilayer, but no significant redistribution (50 ± 4% dC_18:1) in the dC_18:1+dC_22:1 bilayer; these simulated values are within the 95% confidence intervals of the experimental averages. The strong preference for the better matching lipid (dC_16:1) near the channel in the dC_16:1+dC_24:1 mixture and lesser redistribution in the dC_18:1+dC_22:1 mixture can be explained by the energetic cost associated with compressing the lipids to match the channel's hydrophobic length.

Lipid and Peptide Diffusion in Bilayers: The Saffman-Delbrück Model and Periodic Boundary Conditions

The periodic Saffman-Delbrück (PSD) model, an extension of the Saffman-Delbrück model developed to describe the effects of periodic boundary conditions on the diffusion constants of lipids and proteins obtained from simulation, is tested using the coarse-grained Martini and all-atom CHARMM36 (C36) force fields. Simulations of pure Martini dipalmitoylphosphatidylcholine (DPPC) bilayers and those with one embedded gramicidin A (gA) dimer or one gA monomer with sizes ranging from 512 to 2048 lipids support the PSD model. Underestimates of D^∞ (the value of the diffusion constant for an infinite system) from the 512-lipid system are 35% for DPPC, 45% for the gA monomer, and 70% for the gA dimer. Simulations of all-atom DPPC and dioleoylphosphatidylcholine (DOPC) bilayers yield diffusion constants not far from experiment. However, the PSD model predicts that diffusion constants at the sizes of the simulation should underestimate experiment by approximately a factor of 3 for DPPC and 2 for DOPC. This likely implies a deficiency in the C36 force field. A Bayesian method for extrapolating diffusion constants of lipids and proteins in membranes obtained from simulation to infinite system size is provided.

LiCl solvation in N-methyl-acetamide (NMA) as a model for understanding Li(+) binding to an amide plane

The thermodynamics of ion solvation in non-aqueous solvents remains of great significance for understanding cellular transport and ion homeostasis for the design of novel ion-selective materials and applications in molecular pharmacology. Molecular simulations play pivotal roles in connecting experimental measurements to the microscopic structures of liquids. One of the most useful and versatile mimetic systems for understanding biological ion transport is N-methyl-acetamide (NMA). A plethora of theoretical studies for ion solvation in NMA have appeared recently, but further progress is limited by two factors. One is an apparent lack of experimental data on solubility and thermodynamics of solvation for a broad panel of 1 : 1 salts over an appropriate temperature and concentration range. The second concern is more substantial and has to do with the limitations hardwired in the additive (fixed charge) approximations used for most of the existing force-fields. In this submission, we report on the experimental evaluation of LiCl solvation in NMA over a broad range of concentrations and temperatures and compare the results with those of MD simulations with several additive and one polarizable force-field (Drude). By comparing our simulations and experimental results to density functional theory computations, we discuss the limiting factors in existing potential functions. To evaluate the possible implications of explicit and implicit polarizability treatments on ion permeation across biological channels, we performed potential of mean force (PMF) computations for Li(+) transport through a model narrow ion channel with additive and polarizable force-fields.

Computational Investigation of the Effect of Lipid Membranes on Ion Permeation in Gramicidin A

Membrane proteins are embedded in a lipid bilayer and interact with the lipid molecules in subtle ways. This can be studied experimentally by examining the effect of different lipid bilayers on the function of membrane proteins. Understanding the causes of the functional effects of lipids is difficult to dissect experimentally but more amenable to a computational approach. Here we perform molecular dynamics simulations and free energy calculations to study the effect of two lipid types (POPC and NODS) on the conductance of the gramicidin A (gA) channel. A larger energy barrier is found for the K⁺ potential of mean force in gA embedded in POPC compared to that in NODS, which is consistent with the enhanced experimental conductance of cations in gA embedded in NODS. Further analysis of the contributions to the potential energy of K⁺ reveals that gA and water molecules in gA make similar contributions in both bilayers but there are significant differences between the two bilayers when the lipid molecules and interfacial waters are considered. It is shown that the stronger dipole moments of the POPC head groups create a thicker layer of interfacial waters with better orientation, which ultimately is responsible for the larger energy barrier in the K⁺ PMF in POPC.

Gramicidin A Mutants with Antibiotic Activity against Both Gram-Positive and Gram-Negative Bacteria

Antimicrobial peptides (AMPs) have shown potential as alternatives to traditional antibiotics for fighting infections caused by antibiotic-resistant bacteria. One promising example of this is gramicidin A (gA). In its wild-type sequence, gA is active by permeating the plasma membrane of Gram-positive bacteria. However, gA is toxic to human red blood cells at similar concentrations to those required for it to exert its antimicrobial effects. Installing cationic side chains into gA has been shown to lower its hemolytic activity while maintaining the antimicrobial potency. In this study, we present the synthesis and the antibiotic activity of a new series of gA mutants that display cationic side chains. Specifically, by synthesizing alkylated lysine derivatives through reductive amination, we were able to create a broad selection of structures with varied activities towards Staphylococcus aureus and methicillin-resistant S. aureus (MRSA). Importantly, some of the new mutants were observed to have an unprecedented activity towards important Gram-negative pathogens, including Escherichia coli, Klebsiella pneumoniae and Psuedomonas aeruginosa.

Distinguishing gramicidin D conformers through two-dimensional infrared spectroscopy of vibrational excitons

Gramicidin D is a short peptide which dimerizes to form helical pores, adopting one of two conformations in the process. These conformations differ primarily in number of residues per turn and the hydrogen-bond registry between rungs of the helix. Using amide I 2D infrared (IR) and FTIR, we have demonstrated that it is possible to distinguish between the different conformers of gramicidin D in solution. We show that the spectra observed for this helical peptide bear no resemblance to the spectra of α- or 310-helices and that while the FTIR spectra appear similar to spectra of β-sheets, 2D IR reveals that the observed resonances arise from vibrational modes unlike those observed in β-sheets. We also present an idealized model which reproduces the experimental data with high fidelity. This model is able to explain the polarization-dependence of the experimental 2D IR data. Using this model, we show the coupling between the rungs of the helix dominates the spectra, and as a consequence of this, the number of residues per turn can greatly influence the amide I spectra of gramicidin D.

Gramicidin D enhances the antibacterial activity of fluoride

Fluoride is a toxic anion found in many natural environments. One of the major bacterial defenses against fluoride is the cell envelope, which limits passage of the membrane-impermeant fluoride anion. Accordingly, compounds that enhance the permeability of bacterial membranes to fluoride should also enhance fluoride toxicity. In this study, we demonstrate that the pore-forming antibiotic gramicidin D increases fluoride uptake in Bacillus subtilis and that the antibacterial activity of this compound is potentiated by fluoride. Polymyxin B, another membrane-targeting antibiotic with a different mechanism of action, shows no such improvement. These results, along with previous findings, indicate that certain compounds that destabilize bacterial cell envelopes can enhance the toxicity of fluoride.

Solid state NMR: The essential technology for helical membrane protein structural characterization

NMR spectroscopy of helical membrane proteins has been very challenging on multiple fronts. The expression and purification of these proteins while maintaining functionality has consumed countless graduate student hours. Sample preparations have depended on whether solution or solid-state NMR spectroscopy was to be performed - neither have been easy. In recent years it has become increasingly apparent that membrane mimic environments influence the structural result. Indeed, in these recent years we have rediscovered that Nobel laureate, Christian Anfinsen, did not say that protein structure was exclusively dictated by the amino acid sequence, but rather by the sequence in a given environment (Anfinsen, 1973) [106]. The environment matters, molecular interactions with the membrane environment are significant and many examples of distorted, non-native membrane protein structures have recently been documented in the literature. However, solid-state NMR structures of helical membrane proteins in proteoliposomes and bilayers are proving to be native structures that permit a high resolution characterization of their functional states. Indeed, solid-state NMR is uniquely able to characterize helical membrane protein structures in lipid environments without detergents. Recent progress in expression, purification, reconstitution, sample preparation and in the solid-state NMR spectroscopy of both oriented samples and magic angle spinning samples has demonstrated that helical membrane protein structures can be achieved in a timely fashion. Indeed, this is a spectacular opportunity for the NMR community to have a major impact on biomedical research through the solid-state NMR spectroscopy of these proteins.

Small-molecule photostabilizing agents are modifiers of lipid bilayer properties

Small-molecule photostabilizing or protective agents (PAs) provide essential support for the stability demands on fluorescent dyes in single-molecule spectroscopy and fluorescence microscopy. These agents are employed also in studies of cell membranes and model systems mimicking lipid bilayer environments, but there is little information about their possible effects on membrane structure and physical properties. Given the impact of amphipathic small molecules on bilayer properties such as elasticity and intrinsic curvature, we investigated the effects of six commonly used PAs--cyclooctatetraene (COT), para-nitrobenzyl alcohol (NBA), Trolox (TX), 1,4-diazabicyclo[2.2.2]octane (DABCO), para-nitrobenzoic acid (pNBA), and n-propyl gallate (nPG)--on bilayer properties using a gramicidin A (gA)-based fluorescence quench assay to probe for PA-induced changes in the gramicidin monomer↔dimer equilibrium. The experiments were done using fluorophore-loaded large unilamellar vesicles that had been doped with gA, and changes in the gA monomer↔dimer equilibrium were assayed using a gA channel-permeable fluorescence quencher (Tl⁺). Changes in bilayer properties caused by, e.g., PA adsorption at the bilayer/solution interface that alter the equilibrium constant for gA channel formation, and thus the number of conducting gA channels in the large unilamellar vesicle membrane, will be detectable as changes in the rate of Tl⁺ influx-the fluorescence quench rate. Over the experimentally relevant millimolar concentration range, TX, NBA, and pNBA, caused comparable increases in gA channel activity. COT, also in the millimolar range, caused a slight decrease in gA channel activity. nPG increased channel activity at submillimolar concentrations. DABCO did not alter gA activity. Five of the six tested PAs thus alter lipid bilayer properties at experimentally relevant concentrations, which becomes important for the design and analysis of fluorescence studies in cells and model membrane systems. We therefore tested combinations of COT, NBA, and TX; the combinations altered the fluorescence quench rate less than would be predicted assuming their effects on bilayer properties were additive. The combination of equimolar concentrations of COT and NBA caused minimal changes in the fluorescence quench rate.

Density-functional theory study of gramicidin A ion channel geometry and electronic properties

Understanding the mechanisms underlying ion channel function from the atomic-scale requires accurate ab initio modelling as well as careful experiments. Here, we present a density functional theory (DFT) study of the ion channel gramicidin A (gA), whose inner pore conducts only monovalent cations and whose conductance has been shown to depend on the side chains of the amino acids in the channel. We investigate the ground state geometry and electronic properties of the channel in vacuum, focusing on their dependence on the side chains of the amino acids. We find that the side chains affect the ground state geometry, while the electrostatic potential of the pore is independent of the side chains. This study is also in preparation for a full, linear scaling DFT study of gA in a lipid bilayer with surrounding water. We demonstrate that linear scaling DFT methods can accurately model the system with reasonable computational cost. Linear scaling DFT allows ab initio calculations with 10,000-100,000 atoms and beyond, and will be an important new tool for biomolecular simulations.

Atomic-level simulation of current-voltage relationships in single-file ion channels

The difficulty in characterizing ion conduction through membrane channels at the level of individual permeation events has made it challenging to elucidate the mechanistic principles underpinning this fundamental physiological process. Using long, all-atom simulations enabled by special-purpose hardware, we studied K⁺ permeation across the K_V1.2/2.1 voltage-gated potassium channel. At experimentally accessible voltages, which include the physiological range, the simulated permeation rate was substantially lower than the experimentally observed rate. The current–voltage relationship was also nonlinear but became linear at much higher voltages. We observed permeation consistent with a “knock-on” mechanism at all voltages. At high voltages, the permeation rate was in accordance with our previously reported K_V1.2 pore-only simulations, after the simulated voltages from the previous study were recalculated using the correct method, new insight into which is provided here. Including the voltage-sensing domains in the simulated channel brought the linear current–voltage regime closer to the experimentally accessible voltages. The simulated permeation rate, however, still underestimated the experimental rate, because formation of the knock-on intermediate occurred too infrequently. Reducing the interaction strength between the ion and the selectivity filter did not increase conductance. In complementary simulations of gramicidin A, similar changes in interaction strength did increase the observed permeation rate. Permeation nevertheless remained substantially below the experimental value, largely because of infrequent ion recruitment into the pore lumen. Despite the need to apply large voltages to simulate the permeation process, the apparent voltage insensitivity of the permeation mechanism suggests that the direct simulation of permeation at the single-ion level can provide fundamental physiological insight into ion channel function. Notably, our simulations suggest that the knock-on permeation mechanisms in K_V1.2 and KcsA may be different.

Dynamic Heterogeneous Dielectric Generalized Born (DHDGB): An implicit membrane model with a dynamically varying bilayer thickness

An extension to the heterogeneous dielectric generalized Born (HDGB) implicit membrane formalism is presented to allow dynamic membrane deformations in response to membrane-inserted biomolecules during molecular dynamic simulations. The flexible membrane is implemented through additional degrees of freedom that represent the membrane deformation at the contact points of a membrane-inserted solute with the membrane. The extra degrees of freedom determine the dielectric and non-polar solvation free energy profiles that are used to obtain the solvation free energy in the presence of the membrane and are used to calculate membrane deformation free energies according to an elastic membrane model. With the dynamic HDGB (DHDGB) model the membrane is able to deform in response to the insertion of charged molecules thereby avoiding the overestimation of insertion free energies with static membrane models. The DHDGB model also allows the membrane to respond to the insertion of membrane-spanning solutes with hydrophobic mismatch. The model is tested with the membrane insertion of amino acid side chain analogs, arginine-containing helices, the WALP23 peptide, and the gramicidin A channel.

Assessing smectic liquid-crystal continuum models for elastic bilayer deformations

For four decades, since W. Helfrich's pioneering study of smectic A liquid crystals in 1973, continuum elastic models (CEMs) have been employed as tools to understand the energetics of protein-induced lipid bilayer deformations. Among the assumptions underlying this use is that all relevant protein-lipid interactions can be included in the continuum representation of the protein-bilayer interactions through the physical parameters determined for protein-free bilayers and the choice of boundary conditions at the protein/bilayer interface. To better understand this assumption, we review the general structure of CEMs, examine how different choices of boundary conditions and physical moduli profiles alter the predicted bilayer thickness profiles around gramicidin A (gA) and mitochondrial voltage-dependent anion channels (VDAC), respectively, and compare these profiles with those obtained from all-atom molecular dynamics simulations. We find that the profiles differ qualitatively in the first lipid shell around the channels, indicating that the CEMs do not capture accurately the consequences of the protein-induced local changes in lipid bilayer dynamics. Therefore, one needs to be careful when interpreting the results of CEM-based analyses of lipid bilayer-membrane protein interactions.

Three-dimensional stress field around a membrane protein: atomistic and coarse-grained simulation analysis of gramicidin A

Using both atomistic and coarse-grained (CG) models, we compute the three-dimensional stress field around a gramicidin A (gA) dimer in lipid bilayers that feature different degrees of negative hydrophobic mismatch. The general trends in the computed stress field are similar at the atomistic and CG levels, supporting the use of the CG model for analyzing the mechanical features of protein/lipid/water interfaces. The calculations reveal that the stress field near the protein-lipid interface exhibits a layered structure with both significant repulsive and attractive regions, with the magnitude of the stress reaching 1000 bar in certain regions. Analysis of density profiles and stress field distributions helps highlight the Trp residues at the protein/membrane/water interface as mechanical anchors, suggesting that similar analysis is useful for identifying tension sensors in other membrane proteins, especially membrane proteins involved in mechanosensation. This work fosters a connection between microscopic and continuum mechanics models for proteins in complex environments and makes it possible to test the validity of assumptions commonly made in continuum mechanics models for membrane mediated processes. For example, using the calculated stress field, we estimate the free energy of membrane deformation induced by the hydrophobic mismatch, and the results for regions beyond the annular lipids are in general consistent with relevant experimental data and previous theoretical estimates using elasticity theory. On the other hand, the assumptions of homogeneous material properties for the membrane and a bilayer thickness at the protein/lipid interface being independent of lipid type (e.g., tail length) appear to be oversimplified, highlighting the importance of annular lipids of membrane proteins. Finally, the stress field analysis makes it clear that the effect of even rather severe hydrophobic mismatch propagates to only about two to three lipid layers, thus putting a limit on the range of cooperativity between membrane proteins in crowded cellular membranes.

Interactions of drugs and amphiphiles with membranes: modulation of lipid bilayer elastic properties by changes in acyl chain unsaturation and protonation

Poly-unsaturated fatty acids (PUFAs) alter the function of many membrane proteins, whereas monounsatured fatty acids generally are inert. We previously showed that docosahexaenoic acid (DHA) at pH 7 decreases the bilayer stiffness, consistent with an amphiphile-induced increase in elasticity, but not with a negative change in curvature; oleic acid (OA) was inert (Bruno, Koeppe and Andersen, Proc. Natl. Acad. Sci., 2007, 104, 9638-9643). To further explore how PUFAs and other amphiphiles may alter lipid bilayer properties, and thus membrane protein function, we examined how changes in acyl chain unsaturation and head group charge and size alter bilayer properties, as sensed by bilayer-spanning gramicidin A (gA) channels of different lengths. Compared to DHA, the neutral DHA-methyl ester has reduced effects on bilayer properties and 1-palmitoyl-2-docosahexaenoyl-phosphatidylcholine (PDPC) forms bilayers that are softer than dioleoylphosphatidylcholine (DOPC). The changes in channel function are larger for the short gA channels, indicating that changes in elasticity dominate over changes in curvature. We altered the fatty acid protonation by titration: docosahexaenoic acid (DHA) is more potent at pH 9 (relative to pH 7) and is inert at pH 4; OA, which was inert at pH 7, becomes a potent modifier of bilayer properties at pH 9. At both pH 7 and 9, DHA and OA produced larger changes in the lifetimes of the short gA channels, demonstrating that they increase lipid bilayer elasticity when deprotonated--though OA promotes the formation of inverted hexagonal phases at pH 7. The positively charged oleylamine (OAm), which has a small head-group and therefore should be a negative curvature promoter, inhibited gA channel function with similar reductions in the lifetimes of the short and long gA channels, indicating a curvature-dominated effect. Monitoring the single-channel conductance, we find that the negatively charged fatty acids increase the conductance by increasing the local negative charge around the channel, whereas the positively charged OAm has no effect. These results suggest that deprotonated fatty acids increase bilayer elasticity by reversibly adsorbing at the bilayer/solution interface.