High-resolution polypeptide structure in a lamellar phase lipid environment from solid state NMR derived orientational constraints

Antiviral activity of the host defense peptide piscidin 1: investigating a membrane-mediated mode of action

Outbreaks of viral diseases are on the rise, fueling the search for antiviral therapeutics that act on a broad range of viruses while remaining safe to human host cells. In this research, we leverage the finding that the plasma membranes of host cells and the lipid bilayers surrounding enveloped viruses differ in lipid composition. We feature Piscidin 1 (P1), a cationic host defense peptide (HDP) that has antimicrobial effects and membrane activity associated with its N-terminal region where a cluster of aromatic residues and copper-binding motif reside. While few HDPs have demonstrated antiviral activity, P1 acts in the micromolar range against several enveloped viruses that vary in envelope lipid composition. Notably, it inhibits HIV-1, a virus that has an envelope enriched in cholesterol, a lipid associated with higher membrane order and stability. Here, we first document through plaque assays that P1 boasts strong activity against SARS-CoV-2, which has an envelope low in cholesterol. Second, we extend previous studies done with homogeneous bilayers and devise cholesterol-containing zwitterionic membranes that contain the liquid disordered (Ld; low in cholesterol) and ordered (Lo, rich in cholesterol) phases. Using dye leakage assays and cryo-electron microscopy on vesicles, we show that P1 has dramatic permeabilizing capability on the Lo/Ld, an effect matched by a strong ability to aggregate, fuse, and thin the membranes. Differential scanning calorimetry and NMR experiments demonstrate that P1 mixes the lipid content of vesicles and alters the stability of the Lo. Structural studies by NMR indicate that P1 interacts with the Lo/Ld by folding into an α-helix that lies parallel to the membrane surface. Altogether, these results show that P1 is more disruptive to phase-separated than homogenous cholesterol-containing bilayers, suggesting an ability to target domain boundaries. Overall, this multi-faceted research highlights how a peptide that interacts strongly with membranes through an aromatic-rich N-terminal motif disrupt viral envelope mimics. This represents an important step towards the development of novel peptides with broad-spectrum antiviral activity.

17O Solid-State NMR Spectroscopy of Lipid Membranes

Despite the limitations posed by poor sensitivity, studies have reported the unique advantages of 17O based NMR spectroscopy to study systems existing in liquid, solid, or semisolid states. 17O NMR studies have exploited the remarkable sensitivity of quadrupole coupling and chemical shift anisotropy tensors to the local environment in the characterization of a variety of intra- and intermolecular interactions and motion. Recent studies have considerably expanded the use of 17O NMR to study dynamic intermolecular interactions associated with some of the challenging biological systems under magic angle spinning (MAS) and aligned conditions. The very fast relaxing nature of 17O has been well utilized in cellular and in vivo MRS (magnetic resonance spectroscopy) and MRI (magnetic resonance imaging) applications. The main focus of this Review is to highlight the new developments in the biological solids with a detailed discussion for a few selected examples including membrane proteins and nanodiscs. In addition to the unique benefits and limitations, the remaining challenges to overcome, and the impacts of higher magnetic fields and sensitivity enhancement techniques are discussed.

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.

Solid-state NMR spectroscopy for structural studies of polypeptides and lipids in extended physiological membranes

Solid-state NMR is a quickly developing technique that allows one to obtain structural information at atomic resolution in extended lipid bilayers in a rather unique manner. Two approaches have been developed for membrane proteins and peptides namely magic angle sample spinning and the use of uniaxially oriented membrane samples. The state-of-the-art of both approaches will be introduced and the perspectives of solid-state NMR spectroscopy in the context of other structural biology techniques, pressing biomedical questions and membrane biophysics will be discussed.

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.

Orexin A, an amphipathic α-helical neuropeptide involved in pleiotropic functions in the nervous and immune systems: Synthetic approach and biophysical studies of the membrane-bound state

This research reports on the membrane interactions of orexin A (OXA), an α-helical and amphipathic neuropeptide that contains 33 residues and two disulfide bonds in the N-terminal region. OXA, which activates the orexins 1 and 2 receptors in neural and immune cell membranes, has essential pleiotropic physiological effects, including at the levels of arousal, sleep/wakefulness, energy balance, neuroprotection, lipid signaling, the inflammatory response, and pain. As a result, the orexin system has become a prominent target to treat diseases such as sleep disorders, drug addiction, and inflammation. While the high-resolution structure of OXA has been investigated in water and bound to micelles, there is a lack of information about its conformation bound to phospholipid membranes and its receptors. NMR is a powerful method to investigate peptide structures in a membrane environment. To facilitate the NMR structural studies of OXA exposed to membranes, we present a novel synthetic scheme, leading to the production of isotopically-labeled material at high purity. A receptor activation assay shows that the ¹⁵N-labeled peptide is biologically active. Biophysical studies are performed using surface plasmon resonance, circular dichroism, and NMR to investigate the interactions of OXA with phospholipid bilayers. The results demonstrate a strong interaction between the peptide and phospholipids, an increase in α-helical content upon membrane binding, and an in-plane orientation of the C-terminal region critical to function. This new knowledge about structure-activity relationships in OXA could inspire the design of novel therapeutics that leverage the anti-inflammatory and neuro-protective functions of OXA, and therefore could help address neuroinflammation, a major issue associated with neurological disorders such as Alzheimer's disease.

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.

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.

A theoretical assessment of structure determination of multi-span membrane proteins by oriented sample solid-state NMR spectroscopy

Oriented sample solid state NMR (OS-ssNMR) spectroscopy allows direct determination of the structure and topology of membrane proteins reconstituted into aligned lipid bilayers. While OS-ssNMR theoretically has no upper size limit, its application to multi-span membrane proteins has not been established since most studies have been restricted to single or dual span proteins and peptides. Here, we present a critical assessment of the application of this method to multi-span membrane proteins. We used molecular dynamics simulations to back-calculate [¹⁵N-¹H] separated local field (SLF) spectra from a G protein-coupled receptor (GPCR) and show that fully resolved spectra can be obtained theoretically for a multi-span membrane protein with currently achievable resonance linewidths.

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.

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.

Bilayer thickness determines the alignment of model polyproline helices in lipid membranes

Our understanding of protein folds relies fundamentally on the set of secondary structures found in the proteomes. Yet, there also exist intriguing structures and motifs that are underrepresented in natural biopolymeric systems. One example is the polyproline II helix, which is usually considered to have a polar character and therefore does not form membrane spanning sections of membrane proteins. In our work, we have introduced specially designed polyproline II helices into the hydrophobic membrane milieu and used 19F NMR to monitor the helix alignment in oriented lipid bilayers. Our results show that these artificial hydrophobic peptides can adopt several different alignment states. If the helix is shorter than the thickness of the hydrophobic core of the membrane, it is submerged into the bilayer with its long axis parallel to the membrane plane. The polyproline helix adopts a transmembrane alignment when its length exceeds the bilayer thickness. If the peptide length roughly matches the lipid thickness, a coexistence of both states is observed. We thus show that the lipid thickness plays a determining role in the occurrence of a transmembrane polyproline II helix. We also found that the adaptation of polyproline II helices to hydrophobic mismatch is in some notable aspects different from α-helices. Finally, our results prove that the polyproline II helix is a competent structure for the construction of transmembrane peptide segments, despite the fact that no such motif has ever been reported in natural systems.

Mechanistic Landscape of Membrane-Permeabilizing Peptides

Membrane permeabilizing peptides (MPPs) are as ubiquitous as the lipid bilayer membranes they act upon. Produced by all forms of life, most membrane permeabilizing peptides are used offensively or defensively against the membranes of other organisms. Just as nature has found many uses for them, translational scientists have worked for decades to design or optimize membrane permeabilizing peptides for applications in the laboratory and in the clinic ranging from antibacterial and antiviral therapy and prophylaxis to anticancer therapeutics and drug delivery. Here, we review the field of membrane permeabilizing peptides. We discuss the diversity of their sources and structures, the systems and methods used to measure their activities, and the behaviors that are observed. We discuss the fact that "mechanism" is not a discrete or a static entity for an MPP but rather the result of a heterogeneous and dynamic ensemble of structural states that vary in response to many different experimental conditions. This has led to an almost complete lack of discrete three-dimensional active structures among the thousands of known MPPs and a lack of useful or predictive sequence-structure-function relationship rules. Ultimately, we discuss how it may be more useful to think of membrane permeabilizing peptides mechanisms as broad regions of a mechanistic landscape rather than discrete molecular processes.

YidC: Evaluating the Importance of the Native Environment

In this issue of Structure, Baker et al. (2018) take advantage of recent technological breakthroughs in solid-state NMR spectroscopy and electron cryogenic tomography to characterize structural and functional differences between reconstituted YidC in E. coli lipids and YidC as overexpressed in the E. coli inner cellular membrane.

Structure determination of supra-molecular assemblies by solid-state NMR: Practical considerations

In the cellular environment, biomolecules assemble in large complexes which can act as molecular machines. Determining the structure of intact assemblies can reveal conformations and inter-molecular interactions that are only present in the context of the full assembly. Solid-state NMR (ssNMR) spectroscopy is a technique suitable for the study of samples with high molecular weight that allows the atomic structure determination of such large protein assemblies under nearly physiological conditions. This review provides a practical guide for the first steps of studying biological supra-molecular assemblies using ssNMR. The production of isotope-labeled samples is achievable via several means, which include recombinant expression, cell-free protein synthesis, extraction of assemblies directly from cells, or even the study of assemblies in whole cells in situ. Specialized isotope labeling schemes greatly facilitate the assignment of chemical shifts and the collection of structural data. Advanced strategies such as mixed, diluted, or segmental subunit labeling offer the possibility to study inter-molecular interfaces. Detailed and practical considerations are presented with respect to first setting up magic-angle spinning (MAS) ssNMR experiments, including the selection of the ssNMR rotor, different methods to best transfer the sample and prepare the rotor, as well as common and robust procedures for the calibration of the instrument. Diagnostic spectra to evaluate the resolution and sensitivity of the sample are presented. Possible improvements that can reduce sample heterogeneity and improve the quality of ssNMR spectra are reviewed.

Gramicidin Lateral Distribution in Phospholipid Membranes: Fluorescence Phasor Plots and Statistical Mechanical Model

When using small mole fraction increments to study gramicidins in phospholipid membranes, we found that the phasor dots of intrinsic fluorescence of gramicidin D and gramicidin A in dimyristoyl-sn-glycero-3-phosphocholine (DMPC) unilamellar and multilamellar vesicles exhibit a biphasic change with peptide content at 0.143 gramicidin mole fraction. To understand this phenomenon, we developed a statistical mechanical model of gramicidin/DMPC mixtures. Our model assumes a sludge-like mixture of fluid phase and aggregates of rigid clusters. In the fluid phase, gramicidin monomers are randomly distributed. A rigid cluster is formed by a gramicidin dimer and DMPC molecules that are condensed to the dimer, following particular stoichiometries (critical gramicidin mole fractions, Xcr including 0.143). Rigid clusters form aggregates in which gramicidin dimers are regularly distributed, in some cases, even to superlattices. At Xcr, the size of cluster aggregates and regular distributions reach a local maximum. Before a similar model was developed for cholesterol/DMPC mixtures (Sugar and Chong (2012) J. Am. Chem. Soc. 134, 1164⁻1171) and here the similarities and differences are discussed between these two models.

Photonic band-gap resonators for high-field/high-frequency EPR of microliter-volume liquid aqueous samples

High-field EPR provides significant advantages for studying structure and dynamics of molecular systems possessing an unpaired electronic spin. However, routine use of high-field EPR in biophysical research, especially for aqueous biological samples, is still facing substantial technical difficulties stemming from high dielectric millimeter wave (mmW) losses associated with non-resonant absorption by water and other polar molecules. The strong absorbance of mmW's by water also limits the penetration depth to just fractions of mm or even less, thus making fabrication of suitable sample containers rather challenging. Here we describe a radically new line of high Q-factor mmW resonators that are based on forming lattice defects in one-dimensional photonic band-gap (PBG) structures composed of low-loss ceramic discs of λ/4 in thickness and having alternating dielectric constants. A sample (either liquid or solid) is placed within the E = 0 node of the standing mm wave confined within the defect. A resonator prototype has been built and tested at 94.3 GHz. The resonator performance is enhanced by employing ceramic nanoporous membranes as flat sample holders of controllable thickness and tunable effective dielectric constant. The experimental Q-factor of an empty resonator was  ≈ 420. The Q-factor decreased slightly to  ≈ 370 when loaded with a water-containing nanoporous disc of 50 μm in thickness. The resonator has been tested with a number of liquid biological samples and demonstrated about tenfold gain in concentration sensitivity vs. a high-Q cylindrical TE₀₁₂-type cavity. Detailed HFSS Ansys simulations have shown that the resonator structure could be further optimized by properly choosing the thickness of the aqueous sample and employing metallized surfaces. The PBG resonator design is readily scalable to higher mmW frequencies and is capable of accommodating significantly larger sample volumes than previously achieved with either Fabry-Perot or cylindrical resonators.

Advances in instrumentation and methodology for solid-state NMR of biological assemblies

Many advances in instrumentation and methodology have furthered the use of solid-state NMR as a technique for determining the structures and studying the dynamics of molecules involved in complex biological assemblies. Solid-state NMR does not require large crystals, has no inherent size limit, and with appropriate isotopic labeling schemes, supports solving one component of a complex assembly at a time. It is complementary to cryo-EM, in that it provides local, atomic-level detail that can be modeled into larger-scale structures. This review focuses on the development of high-field MAS instrumentation and methodology; including probe design, benchmarking strategies, labeling schemes, and experiments that enable the use of quadrupolar nuclei in biomolecular NMR. Current challenges facing solid-state NMR of biological assemblies and new directions in this dynamic research area are also discussed.

1H magic-angle spinning NMR evolves as a powerful new tool for membrane proteins

Building on a decade of continuous advances of the community, the recent development of very fast (60 kHz and above) magic-angle spinning (MAS) probes has revolutionised the field of solid-state NMR. This new spinning regime reduces the ¹H-¹H dipolar couplings, so that direct detection of the larger magnetic moment available from ¹H is now possible at high resolution, not only in deuterated molecules but also in fully-protonated substrates. Such capabilities allow rapid "fingerprinting" of samples with a ten-fold reduction of the required sample amounts with respect to conventional approaches, and permit extensive, robust and expeditious assignment of small-to-medium sized proteins (up to ca. 300 residues), and the determination of inter-nuclear proximities, relative orientations of secondary structural elements, protein-cofactor interactions, local and global dynamics. Fast MAS and ¹H detection techniques have nowadays been shown to be applicable to membrane-bound systems. This paper reviews the strategies underlying this recent leap forward in sensitivity and resolution, describing its potential for the detailed characterization of membrane proteins.

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.

Study of (Cyclic Peptide)-Polymer Conjugate Assemblies by Small-Angle Neutron Scattering

We present a fundamental study into the self-assembly of (cyclic peptide)-polymer conjugates as a versatile supramolecular motif to engineer nanotubes with defined structure and dimensions, as characterised in solution using small-angle neutron scattering (SANS). This work demonstrates the ability of the grafted polymer to stabilise and/or promote the formation of unaggregated nanotubes by the direct comparison to the unconjugated cyclic peptide precursor. This ideal case permitted a further study into the growth mechanism of self-assembling cyclic peptides, allowing an estimation of the cooperativity. Furthermore, we show the dependency of the nanostructure on the polymer and peptide chemical functionality in solvent mixtures that vary in the ability to compete with the intermolecular associations between cyclic peptides and ability to solvate the polymer shell.

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.

Solid-state NMR methods for oriented membrane proteins

Oriented-sample solid-state NMR represents one of few experimental methods capable of characterising the membrane-bound conformation of proteins in the cell membrane. Since the technique was developed 25 years ago, the technique has been applied to study the structure of helix bundle membrane proteins and antimicrobial peptides, characterise protein-lipid interactions, and derive information on dynamics of the membrane anchoring of membrane proteins. We will review the major developments in various aspects of oriented-sample solid-state NMR, including sample-preparation methods, pulse sequences, theory required to interpret the experiments, perspectives for and guidelines to new experiments, and a number of representative applications.

Structural dynamics and conformational equilibria of SERCA regulatory proteins in membranes by solid-state NMR restrained simulations

Solid-state NMR spectroscopy is emerging as a powerful approach to determine structure, topology, and conformational dynamics of membrane proteins at the atomic level. Conformational dynamics are often inferred and quantified from the motional averaging of the NMR parameters. However, the nature of these motions is difficult to envision based only on spectroscopic data. Here, we utilized restrained molecular dynamics simulations to probe the structural dynamics, topology and conformational transitions of regulatory membrane proteins of the calcium ATPase SERCA, namely sarcolipin and phospholamban, in explicit lipid bilayers. Specifically, we employed oriented solid-state NMR data, such as dipolar couplings and chemical shift anisotropy measured in lipid bicelles, to refine the conformational ensemble of these proteins in lipid membranes. The samplings accurately reproduced the orientations of transmembrane helices and showed a significant degree of convergence with all of the NMR parameters. Unlike the unrestrained simulations, the resulting sarcolipin structures are in agreement with distances and angles for hydrogen bonds in ideal helices. In the case of phospholamban, the restrained ensemble sampled the conformational interconversion between T (helical) and R (unfolded) states for the cytoplasmic region that could not be observed using standard structural refinements with the same experimental data set. This study underscores the importance of implementing NMR data in molecular dynamics protocols to better describe the conformational landscapes of membrane proteins embedded in realistic lipid membranes.

Membrane protein structural validation by oriented sample solid-state NMR: diacylglycerol kinase

The validation of protein structures through functional assays has been the norm for many years. Functional assays perform this validation for water-soluble proteins very well, but they need to be performed in the same environment as that used for the structural analysis. This is difficult for membrane proteins that are often structurally characterized in detergent environments, although functional assays for these proteins are most frequently performed in lipid bilayers. Because the structure of membrane proteins is known to be sensitive to the membrane mimetic environment, such functional assays are appropriate for validating the protein construct, but not the membrane protein structure. Here, we compare oriented sample solid-state NMR spectral data of diacylglycerol kinase previously published with predictions of such data from recent structures of this protein. A solution NMR structure of diacylglycerol kinase has been obtained in detergent micelles and three crystal structures have been obtained in a monoolein cubic phase. All of the structures are trimeric with each monomer having three transmembrane and one amphipathic helices. However, the solution NMR structure shows typical perturbations induced by a micelle environment that is reflected in the predicted solid-state NMR resonances from the structural coordinates. The crystal structures show few such perturbations, especially for the wild-type structure and especially for the monomers that do not have significant crystal contacts. For these monomers the predicted and observed data are nearly identical. The thermostabilized constructs do show more perturbations, especially the A41C mutation that introduces a hydrophilic residue into what would be the middle of the lipid bilayer inducing additional hydrogen bonding between trimers. These results demonstrate a general technique for validating membrane protein structures with minimal data obtained from membrane proteins in liquid crystalline lipid bilayers by oriented sample solid-state NMR.

Idealizing ion channel recordings by a jump segmentation multiresolution filter

Based on a combination of jump segmentation and statistical multiresolution analysis for dependent data, a new approach called J-SMURF to idealize ion channel recordings has been developed. It is model-free in the sense that no a-priori assumptions about the channel’s characteristics have to be made; it thus complements existing methods which assume a model for the channel's dynamics, like hidden Markov models. The method accounts for the effect of an analog filter being applied before the data analysis, which results in colored noise, by adapting existing muliresolution statistics to this situation. J-SMURF’s ability to denoise the signal without missing events even when the signal-to-noise ratio is low is demonstrated on simulations as well as on ion current traces obtained from gramicidin A channels reconstituted into solvent-free planar membranes. When analyzing a newly synthesized acylated system of a fatty acid modified gramicidin channel, we are able to give statistical evidence for unknown gating characteristics such as subgating.

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.

Efficient and stable reconstitution of the ABC transporter BmrA for solid-state NMR studies

We present solid-state NMR sample preparation and first 2D spectra of the Bacillus subtilis ATP-binding cassette (ABC) transporter BmrA, a membrane protein involved in multidrug resistance. The homodimeric 130-kDa protein is a challenge for structural characterization due to its membrane-bound nature, size, inherent flexibility and insolubility. We show that reconstitution of this protein in lipids from Bacillus subtilis at a lipid-protein ratio of 0.5 w/w allows for optimal protein insertion in lipid membranes with respect to two central NMR requirements, high signal-to-noise in the spectra and sample stability over a time period of months. The obtained spectra point to a well-folded protein and a highly homogenous preparation, as witnessed by the narrow resonance lines and the signal dispersion typical for the expected secondary structure distribution of BmrA. This opens the way for studies of the different conformational states of the transporter in the export cycle, as well as on interactions with substrates, via chemical-shift fingerprints and sequential resonance assignments.

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.

Helical membrane protein conformations and their environment

Evidence that membrane proteins respond conformationally and functionally to their environment is growing. Structural models, by necessity, have been characterized in preparations where the protein has been removed from its native environment. Different structural methods have used various membrane mimetics that have recently included lipid bilayers as a more native-like environment. Structural tools applied to lipid bilayer-embedded integral proteins are informing us about important generic characteristics of how membrane proteins respond to the lipid environment as compared with their response to other nonlipid environments. Here, we review the current status of the field, with specific reference to observations of some well-studied α-helical membrane proteins, as a starting point to aid the development of possible generic principles for model refinement.

Solid state NMR strategy for characterizing native membrane protein structures

Unlike water soluble proteins, the structures of helical transmembrane proteins depend on a very complex environment. These proteins sit in the midst of dramatic electrical and chemical gradients and are often subject to variations in the lateral pressure profile, order parameters, dielectric constant, and other properties. Solid state NMR is a collection of tools that can characterize high resolution membrane protein structure in this environment. Indeed, prior work has shown that this complex environment significantly influences transmembrane protein structure. Therefore, it is important to characterize such structures under conditions that closely resemble its native environment. Researchers have used two approaches to gain protein structural restraints via solid state NMR spectroscopy. The more traditional approach uses magic angle sample spinning to generate isotropic chemical shifts, much like solution NMR. As with solution NMR, researchers can analyze the backbone chemical shifts to obtain torsional restraints. They can also examine nuclear spin interactions between nearby atoms to obtain distances between atomic sites. Unfortunately, for membrane proteins in lipid preparations, the spectral resolution is not adequate to obtain complete resonance assignments. Researchers have developed another approach for gaining structural restraints from membrane proteins: the use of uniformly oriented lipid bilayers, which provides a method for obtaining high resolution orientational restraints. When the bilayers are aligned with respect to the magnetic field of the NMR spectrometer, researchers can obtain orientational restraints in which atomic sites in the protein are restrained relative to the alignment axis. However, this approach does not allow researchers to determine the relative packing between helices. By combining the two approaches, we can take advantage of the information acquired from each technique to minimize the challenges and maximize the quality of the structural results. By combining the distance, torsional, and orientational restraints, we can characterize high resolution membrane protein structure in native-like lipid bilayer environments.

Solid state NMR and protein-protein interactions in membranes

Solid state NMR spectroscopy has evolved rapidly in recent years into an excellent tool for the characterization of membrane proteins and their complexes. In the past few years it has also become clear that the structure of membrane proteins, especially helical membrane proteins is determined, in part, by the membrane environment. Therefore, the modeling of this environment by a liquid crystalline lipid bilayer for solid state NMR has generated a unique tool for the characterization of native conformational states, local and global dynamics, and high-resolution structure for these proteins. Protein-protein interactions can also benefit from this solid state NMR capability to characterize membrane proteins in a native-like environment. These complexes take the form of oligomeric structures and hetero-protein interactions both with water-soluble proteins and other membrane proteins.

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.

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.

Glycines: role in α-helical membrane protein structures and a potential indicator of native conformation

Among the growing number of membrane protein structures in the Protein Data Bank, there are many transmembrane domains that appear to be native-like; at the same time, there are others that appear to have less than complete native-like character. Hence, there is an increasing need for validation tools that distinguish native-like from non-native-like structures. Membrane mimetics used in protein structural characterizations differ in numerous physicochemical properties from native membranes and provide many opportunities for introducing non-native-like features into membrane protein structures. One possible approach for validating membrane protein structures is based on the use of glycine residues in transmembrane domains. Here, we have reviewed the membrane protein structure database and identified a set of benchmark proteins that appear to be native-like. In these structures, conserved glycine residues rarely face the lipid interstices, and many of them participate in close helix-helix packing. Glycine-based validation allowed the identification of non-native-like features in several membrane proteins and also shows the potential for verifying the native-like character for numerous other membrane protein structures.

Ion fluxes through nanopores and transmembrane channels

We introduce an implicit solvent Molecular Dynamics approach for calculating ionic fluxes through narrow nanopores and transmembrane channels. The method relies on a dual-control-volume grand-canonical molecular dynamics (DCV-GCMD) simulation and the analytical solution for the electrostatic potential inside a cylindrical nanopore recently obtained by Levin [Europhys. Lett. 76, 163 (2006)]. The theory is used to calculate the ionic fluxes through an artificial transmembrane channel which mimics the antibacterial gramicidin A channel. Both current-voltage and current-concentration relations are calculated under various experimental conditions. We show that our results are comparable to the characteristics associated to the gramicidin A pore, especially the existence of two binding sites inside the pore and the observed saturation in the current-concentration profiles.

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. 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.