Structure of gramicidin a in a lipid bilayer environment determined using molecular dynamics simulations and solid-state NMR data

Advancements in Engineering Planar Model Cell Membranes: Current Techniques, Applications, and Future Perspectives

Cell membranes are crucial elements in living organisms, serving as protective barriers and providing structural support for cells. They regulate numerous exchange and communication processes between cells and their environment, including interactions with other cells, tissues, ions, xenobiotics, and drugs. However, the complexity and heterogeneity of cell membranes-comprising two asymmetric layers with varying compositions across different cell types and states (e.g., healthy vs. diseased)-along with the challenges of manipulating real cell membranes represent significant obstacles for in vivo studies. To address these challenges, researchers have developed various methodologies to create model cell membranes or membrane fragments, including mono- or bilayers organized in planar systems. These models facilitate fundamental studies on membrane component interactions as well as the interactions of membrane components with external agents, such as drugs, nanoparticles (NPs), or biomarkers. The applications of model cell membranes have extended beyond basic research, encompassing areas such as biosensing and nanoparticle camouflage to evade immune detection. In this review, we highlight advancements in the engineering of planar model cell membranes, focusing on the nanoarchitectonic tools used for their fabrication. We also discuss approaches for incorporating challenging materials, such as proteins and enzymes, into these models. Finally, we present our view on future perspectives in the field of planar model cell membranes.

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.

Dimeric Tubulin Modifies Mechanical Properties of Lipid Bilayer, as Probed Using Gramicidin A Channel

Using the gramicidin A channel as a molecular probe, we show that tubulin binding to planar lipid membranes changes the channel kinetics-seen as an increase in the lifetime of the channel dimer-and thus points towards modification of the membrane's mechanical properties. The effect is more pronounced in the presence of non-lamellar lipids in the lipid mixture used for membrane formation. To interpret these findings, we propose that tubulin binding redistributes the lateral pressure of lipid packing along the membrane depth, making it closer to the profile expected for lamellar lipids. This redistribution happens because tubulin perturbs the lipid headgroup spacing to reach the membrane's hydrophobic core via its amphiphilic α-helical domain. Specifically, it increases the forces of repulsion between the lipid headgroups and reduces such forces in the hydrophobic region. We suggest that the effect is reciprocal, meaning that alterations in lipid bilayer mechanics caused by membrane remodeling during cell proliferation in disease and development may also modulate tubulin membrane binding, thus exerting regulatory functions. One of those functions includes the regulation of protein-protein interactions at the membrane surface, as exemplified by VDAC complexation with tubulin.

Calcium-gated potassium channel blockade via membrane-facing fenestrations

Quaternary ammonium blockers were previously shown to bind in the pore to block both open and closed conformations of large-conductance calcium-activated potassium (BK and MthK) channels. Because blocker entry was assumed through the intracellular entryway (bundle crossing), closed-pore access suggested that the gate was not at the bundle crossing. Structures of closed MthK, a Methanobacterium thermoautotrophicum homolog of BK channels, revealed a tightly constricted intracellular gate, leading us to investigate the membrane-facing fenestrations as alternative pathways for blocker access directly from the membrane. Atomistic free energy simulations showed that intracellular blockers indeed access the pore through the fenestrations, and a mutant channel with narrower fenestrations displayed no closed-state TPeA block at concentrations that blocked the wild-type channel. Apo BK channels display similar fenestrations, suggesting that blockers may use them as access paths into closed channels. Thus, membrane fenestrations represent a non-canonical pathway for selective targeting of specific channel conformations, opening novel ways to selectively drug BK channels.

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.

Peptide-induced membrane elastic deformations decelerate gramicidin dimer-monomer equilibration

Gramicidin A (gA) is a hydrophobic pentadecapeptide readily incorporating into a planar bilayer lipid membrane (BLM), thereby inducing a large macroscopic current across the BLM. This current results from ion-channel formation due to head-to-head transbilayer dimerization of gA monomers with rapidly established monomer-dimer equilibrium. Any disturbance of the equilibrium, e.g., by sensitized photoinactivation of a portion of gA monomers, causes relaxation toward a new equilibrium state. According to previous studies, the characteristic relaxation time of the gA-mediated electric current decreases as the current increases upon elevating the gA concentration in the membrane. Here, we report data on the current relaxation kinetics for gA analogs with N-terminal valine replaced by glycine or tyrosine. Surprisingly, the relaxation time increased rather than decreased upon elevation of the total membrane conductance induced by these gA analogs, thus contradicting the classical kinetic scheme. We developed a general theoretical model that accounts for lateral interaction of monomers and dimers mediated by membrane elastic deformations. The modified kinetic scheme of the gramicidin dimerization predicts the reverse dependence of the relaxation time on membrane conductance for gA analogs, with a decreased dimerization constant that is in a good agreement with our experimental data. The equilibration process may be also modulated by incorporation of other peptides ("impurities") into the lipid membrane.

Collective dynamics in lipid membranes containing transmembrane peptides

Biological membranes are composed of complex mixtures of lipids and proteins that influence each other's structure and function. The biological activities of many channel-forming peptides and proteins are known to depend on the material properties of the surrounding lipid bilayer. However, less is known about how membrane-spanning channels affect the lipid bilayer properties, and in particular, their collective fluctuation dynamics. Here we use neutron spin echo spectroscopy (NSE) to measure the collective bending and thickness fluctuation dynamics in dimyristoylphosphatidylcholine (di 14 : 0 PC, DMPC) lipid membranes containing two different antimicrobial peptides, alamethicin (Ala) and gramicidin (gD). Ala and gD are both well-studied antimicrobial peptides that form oligomeric membrane-spanning channels with different structures. At low concentrations, the peptides did not have a measurable effect on the average bilayer structure, yet significantly changed the collective membrane dynamics. Despite both peptides forming transmembrane channels, they had opposite effects on the relaxation time of the collective bending fluctuations and associated effective bending modulus, where gD addition stiffened the membrane while Ala addition softened the membrane. Meanwhile, the lowest gD concentrations enhanced the collective thickness fluctuation dynamics, while the higher gD concentrations and all studied Ala concentrations dampened these dynamics. The results highlight the synergy between lipids and proteins in determining the collective membrane dynamics and that not all peptides can be universally treated as rigid bodies when considering their effects on the lipid bilayer fluctuations.

Characterizing the Water Wire in the Gramicidin Channel Found by Monte Carlo Sampling Using Continuum Electrostatics and in Molecular Dynamics Trajectories with Conventional or Polarizable Force Fields

Water molecules play a key role in all biochemical processes. They help define the shape of proteins, and they are reactant or product in many reactions and are released as ligands are bound. They facilitate the transfer of protons through transmembrane proton channel, pump and transporter proteins. Continuum electrostatics (CE) force fields used by program Multiconformation CE (MCCE) capture electrostatic interactions in biomolecules with an implicit solvent, which captures the averaged solvent water equilibrium properties. Hybrid CE methods can use explicit water molecules within the protein surrounded by implicit solvent. These hybrid methods permit the study of explicit hydrogen bond networks within the protein and allow analysis of processes such as proton transfer reactions. Yet hybrid CE methods have not been rigorously tested. Here, we present an explicit treatment of water molecules in the Gramicidin A (gA) channel using MCCE and compare the resulting distributions of water molecules and key hydration features against those obtained with explicit solvent Molecular Dynamics (MD) simulations with the nonpolarizable CHARMM36 and polarizable Drude force fields. CHARMM36 leads to an aligned water wire in the channel characterized by a large absolute net water dipole moment; the MCCE and Drude analysis lead to a small net dipole moment as the water molecules change orientation within the channel. The correct orientation is not as yet known, so these calculations identify an open question.

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.

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.

Membrane interactions in drug delivery: Model cell membranes and orthogonal techniques

To generate the desired effect in the human body, the active pharmaceutical ingredient usually needs to interact with a receptor located on the cell membrane or inside the cell. Thus, understanding membrane interactions is of great importance when it comes to the development and testing of new drug molecules or new drug delivery systems. Nowadays, there is a tremendous selection of both model cell membranes and of techniques that can be used to characterize interactions between selected model cell membranes and a drug molecule, an excipient, or a drug delivery system. Having such a wide selection of model cell membranes and techniques available makes it sometimes challenging to select the optimal combination for a specific study. Furthermore, it is difficult to compare results obtained using different model cell membranes and techniques, and not all in vitro studies translate as well to an estimation of the in vivo biological activity or understanding of mode of action. This review provides an overview of the available lipid bilayer-based model cell membranes and of the most widely employed techniques for studying membrane interactions. Finally, the need for employing complimentary characterization techniques in order to acquire more reliable and in-depth information is highlighted.

A Computational Pipeline to Predict Cardiotoxicity: From the Atom to the Rhythm

RATIONALE

Drug-induced proarrhythmia is so tightly associated with prolongation of the QT interval that QT prolongation is an accepted surrogate marker for arrhythmia. But QT interval is too sensitive a marker and not selective, resulting in many useful drugs eliminated in drug discovery.

OBJECTIVE

To predict the impact of a drug from the drug chemistry on the cardiac rhythm.

METHODS AND RESULTS

In a new linkage, we connected atomistic scale information to protein, cell, and tissue scales by predicting drug-binding affinities and rates from simulation of ion channel and drug structure interactions and then used these values to model drug effects on the hERG channel. Model components were integrated into predictive models at the cell and tissue scales to expose fundamental arrhythmia vulnerability mechanisms and complex interactions underlying emergent behaviors. Human clinical data were used for model framework validation and showed excellent agreement, demonstrating feasibility of a new approach for cardiotoxicity prediction.

CONCLUSIONS

We present a multiscale model framework to predict electrotoxicity in the heart from the atom to the rhythm. Novel mechanistic insights emerged at all scales of the system, from the specific nature of proarrhythmic drug interaction with the hERG channel, to the fundamental cellular and tissue-level arrhythmia mechanisms. Applications of machine learning indicate necessary and sufficient parameters that predict arrhythmia vulnerability. We expect that the model framework may be expanded to make an impact in drug discovery, drug safety screening for a variety of compounds and targets, and in a variety of regulatory processes.

Simulating ion channel activation mechanisms using swarms of trajectories

Atomic-level studies of protein activity represent a significant challenge as a result of the complexity of conformational changes occurring on wide-ranging timescales, often greatly exceeding that of even the longest simulations. A prime example is the elucidation of protein allosteric mechanisms, where localized perturbations transmit throughout a large macromolecule to generate a response signal. For example, the conversion of chemical to electrical signals during synaptic neurotransmission in the brain is achieved by specialized membrane proteins called pentameric ligand-gated ion channels. Here, the binding of a neurotransmitter results in a global conformational change to open an ion-conducting pore across the nerve cell membrane. X-ray crystallography has produced static structures of the open and closed states of the proton-gated GLIC pentameric ligand-gated ion channel protein, allowing for atomistic simulations that can uncover changes related to activation. We discuss a range of enhanced sampling approaches that could be used to explore activation mechanisms. In particular, we describe recent application of an atomistic string method, based on Roux's "swarms of trajectories" approach, to elucidate the sequence and interdependence of conformational changes during activation. We illustrate how this can be combined with transition analysis and Brownian dynamics to extract thermodynamic and kinetic information, leading to understanding of what controls ion channel function. © 2019 Wiley Periodicals, Inc.

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.

Helix formation and stability in membranes

In this article we review current understanding of basic principles for the folding of membrane proteins, focusing on the more abundant alpha-helical class. Membrane proteins, vital to many biological functions and implicated in numerous diseases, fold into their active conformations in the complex environment of the cell bilayer membrane. While many membrane proteins rely on the translocon and chaperone proteins to fold correctly, others can achieve their functional form in the absence of any translation apparatus or other aides. Nevertheless, the spontaneous folding process is not well understood at the molecular level. Recent findings suggest that helix fraying and loop formation may be important for overall structure, dynamics and regulation of function. Several types of membrane helices with ionizable amino acids change their topology with pH. Additionally we note that some peptides, including many that are rich in arginine, and a particular analogue of gramicidin, are able passively to translocate across cell membranes. The findings indicate that a final protein structure in a lipid-bilayer membrane is sequence-based, with lipids contributing to stability and regulation. While much progress has been made toward understanding the folding process for alpha-helical membrane proteins, it remains a work in progress. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Permeating disciplines: Overcoming barriers between molecular simulations and classical structure-function approaches in biological ion transport

Ion translocation across biological barriers is a fundamental requirement for life. In many cases, controlling this process-for example with neuroactive drugs-demands an understanding of rapid and reversible structural changes in membrane-embedded proteins, including ion channels and transporters. Classical approaches to electrophysiology and structural biology have provided valuable insights into several such proteins over macroscopic, often discontinuous scales of space and time. Integrating these observations into meaningful mechanistic models now relies increasingly on computational methods, particularly molecular dynamics simulations, while surfacing important challenges in data management and conceptual alignment. Here, we seek to provide contemporary context, concrete examples, and a look to the future for bridging disciplinary gaps in biological ion transport. This article is part of a Special Issue entitled: Beyond the Structure-Function Horizon of Membrane Proteins edited by Ute Hellmich, Rupak Doshi and Benjamin McIlwain.

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.

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.

Determinants of Curvature-Sensing Behavior for MARCKS-Fragment Peptides

It is increasingly recognized that membrane curvature plays an important role in various cellular activities such as signaling and trafficking, as well as key issues involving health and disease development. Thus, curvature-sensing peptides are essential to the study and detection of highly curved bilayer structures. The effector domain of myristoylated alanine-rich C-kinase substrate (MARCKS-ED) has been demonstrated to have curvature-sensing ability. Research of the MARCKS-ED has further revealed that its Lys and Phe residues play an essential role in how MARCKS-ED detects and binds to curved bilayers. MARCKS-ED has the added property of being a lower-molecular-weight curvature sensor, which offers advantages in production. With that in mind, this work investigates peptide-sequence-related factors that influence curvature sensing and explores whether peptide fragments of even shorter length can function as curvature sensors. Using both experimental and computational methods, we studied the curvature-sensing capabilities of seven fragments of MARCKS-ED. Two of the longer fragments were designed from approximately the two halves of the full-length peptide whereas the five shorter fragments were taken from the central stretch of MARCKS-ED. Fully atomistic molecular dynamics simulations show that the fragments that remain bound to the bilayer exhibit interactions with the bilayer similar to that of the full-length MARCKS-ED peptide. Fluorescence enhancement and anisotropy assays, meanwhile, reveal that five of the MARCKS fragments possess the ability to sense membrane curvature. Based on the sequences of the curvature-sensing fragments, it appears that the ability to sense curvature involves a balance between the numbers of positively charged residues and hydrophobic anchoring residues. Together, these findings help crystallize our understanding of the molecular mechanisms underpinning the curvature-sensing behaviors of peptides, which will prove useful in the design of future curvature sensors.

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.

Estimation of Potentials of Mean Force from Nonequilibrium Pulling Simulations Using Both Minh-Adib Estimator and Weighted Histogram Analysis Method

Nonequilibrium pulling simulations have been a useful approach for investigating a variety of physical and biological problems. The major target in the simulations is to reconstruct reliable potentials of mean force (PMFs) or unperturbed free-energy profiles for quantitatively addressing both equilibrium mechanistic properties and contributions from nonequilibrium processes. While several current nonequilibrium methods were shown to be accurate in computing free-energy profiles in systems with relatively simple dynamics, they have proved to be unsuitable in more complicated systems. To extend the applicability of nonequilibrium sampling, we demonstrate a novel method that combines Minh-Adib's bidirectional estimator with nonlinear WHAM equations to reconstruct and assess PMFs from relatively fast pulling trajectories. We test the method in a one-dimensional model system and in a system of an antibiotic gramicidin-A (gA) channel, which is considered a significant challenge for nonequilibrium sampling. We identify key parameters for efficiently performing pulling simulations to improve and ensure the convergence and accuracy of estimated PMFs. We show that a few pulling trajectories of a relatively fast pulling speed v = 10 Å/ns can return a fair estimate of the PMF of a single potassium ion in gA.

How transmembrane peptides insert and orientate in biomembranes: a combined experimental and simulation study

After the synthesis of transmembrane peptides/proteins (TMPs), their insertion into a lipid bilayer is a fundamental biophysical process.

Curvature sensing MARCKS-ED peptides bind to membranes in a stereo-independent manner

Membrane curvature and lipid composition plays a critical role in interchanging of matter and energy in cells. Peptide curvature sensors are known to activate signaling pathways and promote molecular transport across cell membranes. Recently, the 25-mer MARCKS-ED peptide, which is derived from the effector domain of the myristoylated alanine-rich C kinase substrate protein, has been reported to selectively recognize highly curved membrane surfaces. Our previous studies indicated that the naturally occurring L-MARCKS-ED peptide could simultaneously detect both phosphatidylserine and curvature. Here, we demonstrate that D-MARCKS-ED, composed by unnatural D-amino acids, has the same activities as its enantiomer, L-MARCKS-ED, as a curvature and lipid sensor. An atomistic molecular dynamics simulation suggests that D-MARCKS-ED may change from linear to a boat conformation upon binding to the membrane. Comparable enhancement of fluorescence intensity was observed between D- and L-MARCKS-ED peptides, indicating similar binding affinities. Meanwhile, circular dichroism spectra of D- and L-MARCKS-ED are almost symmetrical both in the presence and absence of liposomes. These results suggest similar behavior of artificial D- and natural L-MARCKS-ED peptides when binding to curved membranes. Our studies may contribute to further understanding of how MARCKS-ED senses membrane curvature as well as provide a new direction to develop novel membrane curvature probes.

Biophysical investigations with MARCKS-ED: dissecting the molecular mechanism of its curvature sensing behaviors

Curved membranes are a common and important attribute in cells. Protein and peptide curvature sensors are known to activate signaling pathways, initiate vesicle budding, trigger membrane fusion, and facilitate molecular transport across cell membranes. Nonetheless, there is little understanding how these proteins and peptides achieve preferential binding of different membrane curvatures. The current study is to elucidate specific factors required for curvature sensing. As a model system, we employed a recently identified peptide curvature sensor, MARCKS-ED, derived from the effector domain of the myristoylated alanine-rich C kinase substrate protein, for these biophysical investigations. An atomistic molecular dynamics (MD) simulation suggested an important role played by the insertion of the Phe residues within MARCKS-ED. To test these observations from our computational simulations, we performed electron paramagnetic resonance (EPR) studies to determine the insertion depth of MARCKS-ED into differently curved membrane bilayers. Next, studies with varied lipid compositions revealed their influence on curvature sensing by MARCKS-ED, suggesting contributions from membrane fluidity, rigidity, as well as various lipid structures. Finally, we demonstrated that the curvature sensing by MARCKS-ED is configuration independent. In summary, our studies have shed further light to the understanding of how MARCKS-ED differentiates between membrane curvatures, which may be generally applicable to protein curvature sensing behavior.

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.

The influenza fusion peptide promotes lipid polar head intrusion through hydrogen bonding with phosphates and N-terminal membrane insertion depth

Influenza infection requires fusion between the virus envelope and a host cell endosomal membrane. The influenza hemagglutinin fusion peptide (FP) is essential to viral membrane fusion. It was recently proposed that FPs would fuse membranes by increasing lipid tail protrusion, a membrane fusion transition state. The details of how FPs induce lipid tail protrusion, however, remain to be elucidated. To decipher the molecular mechanism by which FPs promote lipid tail protrusion, we performed molecular dynamics simulations of the wild-type (WT) FP, fusogenic mutant F9A, and nonfusogenic mutant W14A in model bilayers. This article presents the peptide-lipid interaction responsible for lipid tail protrusion and a related lipid perturbation, polar head intrusion, where polar heads are sunk under the membrane surface. The backbone amides from the four N-terminal peptide residues, deeply inserted in the membrane, promoted both perturbations through H bonding with lipid phosphates. Polar head intrusion correlated with peptides N-terminal insertion depth and activity: the N-termini of WT and F9A were inserted deeper into the membrane than nonfusogenic W14A. Based on these results, we propose that FP-induced polar head intrusion would complement lipid tail protrusion in catalyzing membrane fusion by reducing repulsions between juxtaposed membranes headgroups. The presented model provides a framework for further research on membrane fusion and influenza antivirals.

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.

An experimental study of the solvent-dependent self-assembly/disassembly and conformer preferences of gramicidin A

The solvent dependence of self-assembly/disassembly kinetics and conformer preferences of the gramicidin A (GA) dimer is investigated using a combination of techniques, viz., electrospray ionization-ion mobility-mass spectrometry (IM-MS), collision-induced dissociation (CID), and hydrogen/deuterium exchange (HDX)-MS. IM-MS measurements reveal that there are possibly three distinct GA dimeric species, detected as sodium ion adduct ions [2GA + 2Na](2+), and these are assigned as the parallel β-helix, antiparallel β-helix, and head-to-head dimer. The monomerization kinetics and equilibrium abundances of the dimer ions depend upon solvent polarity. The antiparallel β-helix was the thermodynamically preferred species in less polar solvents. HDX measurements and collision-induced dissociation (CID) of the intermediate complex confirm the well-protected dimer geometry with strong intermolecular hydrogen bonds. This combined IM-HDX-CID methodology provides a comprehensive view of GA self-assembly/disassembly in low dielectric solutions, showing its potential utility in solving solution-phase protein self-assembly/disassembly kinetics and providing structural information of the multimers at the same time.

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.

Supramolecular structure of membrane-associated polypeptides by combining solid-state NMR and molecular dynamics simulations

Elemental biological functions such as molecular signal transduction are determined by the dynamic interplay between polypeptides and the membrane environment. Determining such supramolecular arrangements poses a significant challenge for classical structural biology methods. We introduce an iterative approach that combines magic-angle spinning solid-state NMR spectroscopy and atomistic molecular dynamics simulations for the determination of the structure and topology of membrane-bound systems with a resolution and level of accuracy difficult to obtain by either method alone. Our study focuses on the Shaker B ball peptide that is representative for rapid N-type inactivating domains of voltage-gated K(+) channels, associated with negatively charged lipid bilayers.

The determinants of hydrophobic mismatch response for transmembrane helices

Hydrophobic mismatch arises from a difference in the hydrophobic thickness of a lipid membrane and a transmembrane protein segment, and is thought to play an important role in the folding, stability and function of membrane proteins. We have investigated the possible adaptations that lipid bilayers and transmembrane α-helices undergo in response to mismatch, using fully-atomistic molecular dynamics simulations totaling 1.4 μs. We have created 25 different tryptophan-alanine-leucine transmembrane α-helical peptide systems, each composed of a hydrophobic alanine-leucine stretch, flanked by 1-4 tryptophan side chains, as well as the β-helical peptide dimer, gramicidin A. Membrane responses to mismatch include changes in local bilayer thickness and lipid order, varying systematically with peptide length. Adding more flanking tryptophan side chains led to an increase in bilayer thinning for negatively mismatched peptides, though it was also associated with a spreading of the bilayer interface. Peptide tilting, bending and stretching were systematic, with tilting dominating the responses, with values of up to ~45° for the most positively mismatched peptides. Peptide responses were modulated by the number of tryptophan side chains due to their anchoring roles and distributions around the helices. Potential of mean force calculations for local membrane thickness changes, helix tilting, bending and stretching revealed that membrane deformation is the least energetically costly of all mismatch responses, except for positively mismatched peptides where helix tilting also contributes substantially. This comparison of energetic driving forces of mismatch responses allows for deeper insight into protein stability and conformational changes in lipid membranes.

The influenza fusion peptide adopts a flexible flat V conformation in membranes

Knowledge about the influenza fusion peptide (FP) membrane insertion mode is crucial for understanding its fusogenic mechanism. NMR and electron paramagnetic resonance experiments showed that in micelles, the FP inserted as a fixed-angle inverted V. In membranes, however, it was shown to insert as a straight α-helix (by molecular-dynamics simulations) and to adopt multiple kinked conformations (by solid-state NMR). In this work we performed explicit-solvent molecular-dynamics simulations of the influenza FP, and its F9A and W14A mutants, in POPC membranes. The Hα1 chemical shifts predicted from the molecular-dynamics structures are in excellent agreement with the experimental values obtained for the three peptides. The peptide orientation and conformations observed from the simulations lead to a flexible flat-V model in which the peptide lies almost flat on the membrane surface and alternates between kinked and straight-helix conformations.

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers

Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4-5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations-except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.

Single-file water in nanopores

Water molecules confined to pores with sub-nanometre diameters form single-file hydrogen-bonded chains. In such nanoscale confinement, water has unusual physical properties that are exploited in biology and hold promise for a wide range of biomimetic and nanotechnological applications. The latter can be realized by carbon and boron nitride nanotubes which confine water in a relatively non-specific way and lend themselves to the study of intrinsic properties of single-file water. As a consequence of strong water-water hydrogen bonds, many characteristics of single-file water are conserved in biological and synthetic pores despite differences in their atomistic structures. Charge transport and orientational order in water chains depend sensitively on and are mainly determined by electrostatic effects. Thus, mimicking functions of biological pores with apolar pores and corresponding external fields gives insight into the structure-function relation of biological pores and allows the development of technical applications beyond the molecular devices found in living systems. In this Perspective, we revisit results for single-file water in apolar pores, and examine the similarities and the differences between these simple systems and water in more complex pores.

Linear rate-equilibrium relations arising from ion channel-bilayer energetic coupling

Linear rate-equilibrium (RE) relations, also known as linear free energy relations, are widely observed in chemical reactions, including protein folding, enzymatic catalysis, and channel gating. Despite the widespread occurrence of linear RE relations, the principles underlying the linear relation between changes in activation and equilibrium energy in macromolecular reactions remain enigmatic. When examining amphiphile regulation of gramicidin channel gating in lipid bilayers, we noted that the gating process could be described by a linear RE relation with a simple geometric interpretation. This description is possible because the gating process provides a well-understood reaction, in which structural changes in a bilayer-embedded model protein can be studied at the single-molecule level. It is thus possible to obtain quantitative information about the energetics of the reaction transition state and its position on a spatial coordinate. It turns out that the linear RE relation for the gramicidin monomer-dimer reaction can be understood, and the quantitative relation between changes in activation energy and equilibrium energy can be interpreted, by considering the effects of amphiphiles on the changes in bilayer elastic energy associated with channel gating. We are not aware that a similar simple mechanistic explanation of a linear RE relation has been provided for a chemical reaction in a macromolecule. RE relations generally should be useful for examining how amphiphile-induced changes in bilayer properties modulate membrane protein folding and function, and for distinguishing between direct (e.g., due to binding) and indirect (bilayer-mediated) effects.

Gramicidin A backbone and side chain dynamics evaluated by molecular dynamics simulations and nuclear magnetic resonance experiments. II: nuclear magnetic resonance experiments

Motional properties are important for understanding protein function and are accessible to NMR relaxation measurements. The goal of this study is to investigate the internal dynamics occurring in gramicidin A (gA) channels in order to provide benchmark experimental data for comparison with the results of molecular dynamics simulations. We therefore synthesized several (15)N isotope-enriched gA samples, covering all backbone residues as well as the Trp indole side chains for NMR relaxation experiments. On the basis of the (15)N NMR spectra for labeled gA samples incorporated in sodium dodecylsulfate (SDS) micelles, we determined T(1), T(2), and heteronuclear NOE values for backbone and indole (15)NH groups. The results indicate that the SDS-incorporated gA channel is a constrained structure without an especially "floppy" region. The NMR observables, particularly those for backbone groups, are predicted well by the molecular dynamics simulations in the accompanying article (DOI 10.1021/jp200904d ).

Gramicidin A backbone and side chain dynamics evaluated by molecular dynamics simulations and nuclear magnetic resonance experiments. I: molecular dynamics simulations

Gramicidin A (gA) channels provide an ideal system to test molecular dynamics (MD) simulations of membrane proteins. The peptide backbone lines a cation-selective pore, and due to the small channel size, the average structure and extent of fluctuations of all atoms in the peptide will influence ion permeation. This raises the question of how well molecular mechanical force fields used in MD simulations and potential of mean force (PMF) calculations can predict structure and dynamics as well as ion permeation. To address this question, we undertook a comparative study of nuclear magnetic resonance (NMR) observables predicted by fully atomistic MD simulations on a gA dimer embedded in a sodium dodecyl sulfate (SDS) micelle with measurements of the gA dimer backbone and tryptophan side chain dynamics using solution-state (15)N NMR on gA dimers in SDS micelles (Vostrikov, V. V.; Gu, H.; Ingólfsson, H. I.; Hinton, J. F.; Andersen, O. S.; Roux, B.; Koeppe, R. E., II. J. Phys. Chem. B2011, DOI 10.1021/jp200906y , accompanying article). This comparison enables us to examine the robustness of the MD simulations done using different force fields as well as their ability to predict important features of the gA channel. We find that MD is able to predict NMR observables, including the generalized order parameters (S(2)), the (15)N spin-lattice (T(1)) and spin-spin (T(2)) relaxation times, and the (1)H-(15)N nuclear Overhauser effect (NOE), with remarkable accuracy. To examine further how differences in the force fields can affect the channel conductance, we calculated the PMF for K(+) and Na(+) permeation through a gA channel in a dimyristoylphosphatidylcholine (DMPC) bilayer. In this case, we find that MD is less successful in quantitatively predicting the single-channel conductance.

Energetics of double-ion occupancy in the gramicidin A channel

To understand the energetics of double-ion occupancy in gramicidin A (gA) channels, the 2D potential of mean force (PMF) is calculated for two ions at different positions along the channel axis. The cross sections of this 2D PMF are compared with available one-ion PMFs to highlight the effect of one ion on the permeation dynamics of the other. It is found that, if the first ion stays on one side in the channel, the second ion has to pass over an additional barrier to move into the outer binding site. At the same time, both outer and inner binding sites for the second ion become shallower than those in the one-ion PMF. The calculated ion-ion repulsion for a doubly occupied channel is about 2 kcal/mol, in good agreement with previous experimental estimates. The number of water molecules inside the channel and their dipole moment are calculated to interpret the energetics of double-ion occupancy. As the first ion moves into the outer binding site and then further into the channel, the oxygen atoms of the single-file water column in the channel are oriented to point toward the ion. The observed dipole moment distribution of a singly occupied channel has only one sharp peak, and the water alignment is essentially perfect once the ion is in the inner binding site. For this reason, there is an energy penalty to accommodate a second ion at the opposite end of the channel.

The structural and topological analysis of membrane-associated polypeptides by oriented solid-state NMR spectroscopy: established concepts and novel developments

Solid-state NMR spectroscopy is a powerful technique for the investigation of membrane-associated peptides and proteins as well as their interactions with lipids, and a variety of conceptually different approaches have been developed for their study. The technique is unique in allowing for the high-resolution investigation of liquid disordered lipid bilayers representing well the characteristics of natural membranes. Whereas magic angle solid-state NMR spectroscopy follows approaches that are related to those developed for solution NMR spectroscopy the use of static uniaxially oriented samples results in angular constraints which also provide information for the detailed analysis of polypeptide structures. This review introduces this latter concept theoretically and provides a number of examples. Furthermore, ongoing developments combining solid-state NMR spectroscopy with information from solution NMR spectroscopy and molecular modelling as well as exploratory studies using dynamic nuclear polarization solid-state NMR will be presented.

Characterization of horizontal lipid bilayers as a model system to study lipid phase separation

Artificial lipid membranes are widely used as a model system to study single ion channel activity using electrophysiological techniques. In this study, we characterize the properties of the artificial bilayer system with respect to its dynamics of lipid phase separation using single-molecule fluorescence fluctuation and electrophysiological techniques. We determined the rotational motions of fluorescently labeled lipids on the nanosecond timescale using confocal time-resolved anisotropy to probe the microscopic viscosity of the membrane. Simultaneously, long-range mobility was investigated by the lateral diffusion of the lipids using fluorescence correlation spectroscopy. Depending on the solvent used for membrane preparation, lateral diffusion coefficients in the range D(lat) = 10-25 mum(2)/s and rotational diffusion coefficients ranging from D(rot) = 2.8 - 1.4 x 10(7) s(-1) were measured in pure liquid-disordered (L(d)) membranes. In ternary mixtures containing saturated and unsaturated phospholipids and cholesterol, liquid-ordered (L(o)) domains segregated from the L(d) phase at 23 degrees C. The lateral mobility of lipids in L(o) domains was around eightfold lower compared to those in the L(d) phase, whereas the rotational mobility decreased by a factor of 1.5. Burst-integrated steady-state anisotropy histograms, as well as anisotropy imaging, were used to visualize the rotational mobility of lipid probes in phase-separated bilayers. These experiments and fluorescence correlation spectroscopy measurements at different focal diameters indicated a heterogeneous microenvironment in the L(o) phase. Finally, we demonstrate the potential of the optoelectro setup to study the influence of lipid domains on the electrophysiological properties of ion channels. We found that the electrophysiological activity of gramicidin A (gA), a well-characterized ion-channel-forming peptide, was related to lipid-domain partitioning. During liquid-liquid phase separation, gA was largely excluded from L(o) domains. Simultaneously, the number of electrically active gA dimers increased due to the increased surface density of gA in the L(d) phase.