Molecular dynamics simulations of model trans-membrane peptides in lipid bilayers: a systematic investigation of hydrophobic mismatch

Exploration of Lipid Bilayer Mechanical Properties Using Molecular Dynamics simulation

Important biological structures known for their exceptional mechanical qualities, lipid bilayers are essential to many cellular functions. Fluidity, elasticity, permeability, stiffness, tensile strength, compressibility, shear viscosity, line tension, and curvature elasticity are some of the fundamental characteristics affecting their behavior. The purpose of this review is to examine these characteristics in more detail by molecular dynamics simulation, elucidating their importance and the elements that lead to their appearance in lipid bilayers. Comprehending the mechanical characteristics of lipid bilayers is critical for creating medications, drug delivery systems, and biomaterials that interact with biological membranes because it allows one to understand how these materials respond to different stresses and deformations. The influence of mechanical characteristics on important lipid bilayer properties is examined in this review. The mechanical properties of lipid bilayers were clarified through the use of molecular dynamics simulation analysis techniques, including bilayer thickness, stress-strain analysis, lipid bilayer area compressibility, membrane bending rigidity, and time- or ensemble-averaged the area per lipid evaluation. We explain the significance of molecular dynamics simulation analysis methods, providing important new information about the stability and dynamic behavior of the bilayer. In the end, we hope to use molecular dynamics simulation to provide a comprehensive understanding of the mechanical properties and behavior of lipid bilayers, laying the groundwork for further studies and applications. Taken together, careful investigation of these mechanical aspects deepens our understanding of the adaptive capacities and functional roles of lipid bilayers in biological environments.

Lipoteichoic Acids Are Essential for Pneumococcal Colonization and Membrane Integrity

INTRODUCTION

The hydrophilic, polymeric chain of the lipoteichoic acid (LTA) of the Gram-positive pathobiont Streptococcus pneumoniae is covalently linked to the glycosylglycerolipid α-d-glucopyranosyl-(1,3)-diacylglycerol by the LTA ligase TacL, leading to its fixation in the cytoplasmic membrane. Pneumococcal LTA, sharing identical repeating units with the wall teichoic acids (WTA), is dispensable for normal growth but required for full virulence in invasive infections.

METHODS

Mutants deficient in TacL and complemented strains constructed were tested for their growth, resistance against oxidative stress, and susceptibility against antimicrobial peptides. Further, the membrane fluidity of pneumococci, their capability to adhere to lung epithelial cells, and virulence in a Galleria mellonella as well as intranasal mouse infection model were assessed.

RESULTS

In the present study, we indicate that LTA is already indispensable for pneumococcal adherence to human nasopharyngeal cells and colonization in an intranasal mouse infection model. Mutants deficient for TacL did not show morphological defects. However, our analysis of pneumococcal membranes in different serotypes showed an altered membrane fluidity and surface protein abundance of lipoproteins in mutants deficient for LTA but not WTA. These mutants had a decreased membrane fluidity, exhibited higher amounts of lipoproteins, and showed an increased susceptibility to antimicrobial peptides. In complemented mutant strains, this defect was fully restored.

CONCLUSION

Taken together, LTA is crucial for colonization and required to effectively protect pneumococci from innate immune defence mechanisms by maintaining the membrane integrity.

Energetics of the Transmembrane Peptide Sorting by Hydrophobic Mismatch

Hydrophobic mismatch between a lipid membrane and embedded transmembrane peptides or proteins plays a role in their lateral localization and function. Earlier studies have resolved numerous mechanisms through which the peptides and membrane proteins adapt to mismatch, yet the energetics of lateral sorting due to hydrophobic mismatch have remained elusive due to the lack of suitable computational or experimental protocols. Here, we pioneer a molecular dynamics simulation approach to study the sorting of peptides along a membrane thickness gradient. Peptides of different lengths tilt and diffuse along the membrane to eliminate mismatch with a rate directly proportional to the magnitude of mismatch. We extract the 2-dimensional free energy profiles as a function of local thickness and peptide orientation, revealing the relative contributions of sorting and tilting, and suggesting their thermally accessible regimes. Our approach can readily be applied to study other membrane systems of biological interest where hydrophobic mismatch, or membrane thickness in general, plays a role.

Sensing membrane voltage by reorientation of dipolar transmembrane peptides

Membrane voltage plays a vital role in the behavior and functions of the lipid bilayer membrane. For instance, it regulates the exchange of molecules across the membrane through transmembrane proteins such as ion channels. In this paper, we study the membrane voltage-sensing mechanism, which entails the reorientation of α-helices with a change in the membrane voltage. We consider a helix having a large electrical macrodipole embedded in a lipid bilayer as a model system. We performed extensive molecular dynamics simulations to study the effect of variation of membrane voltage on the tilt angle of peptides and ascertain the optimal parameters for designing such a voltage-sensing peptide. A theoretical model for the system is also developed to investigate the interplay of competing effects of hydrophobic mismatch and dipole-electric field coupling on the tilt of the peptide and further explore the parameter space. This work opens the possibility for the design and fabrication of artificial dipolar membrane voltage-sensing elements for biomedical applications.

Functional regulation of aquaporin dynamics by lipid bilayer composition

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics are examined. We demonstrate that SoPIP2;1's structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

Studying Conformational Properties of Transmembrane Domain of KCNE3 in a Lipid Bilayer Membrane Using Molecular Dynamics Simulations

KCNE3 is a single-pass integral membrane protein that regulates numerous voltage-gated potassium channel functions such as KCNQ1. Previous solution NMR studies suggested a moderate degree of curved α-helical structure in the transmembrane domain (TMD) of KCNE3 in lyso-myristoylphosphatidylcholine (LMPC) micelles and isotropic bicelles with the residues T71, S74 and G78 situated along the concave face of the curved helix. During the interaction of KCNE3 and KCNQ1, KCNE3 pushes its transmembrane domain against KCNQ1 to lock the voltage sensor in its depolarized conformation. A cryo-EM study of KCNE3 complexed with KCNQ1 in nanodiscs suggested a deviation of the KCNE3 structure from its independent structure in isotropic bicelles. Despite the biological significance of KCNE3 TMD, the conformational properties of KCNE3 are poorly understood. Here, all atom molecular dynamics (MD) simulations were utilized to investigate the conformational dynamics of the transmembrane domain of KCNE3 in a lipid bilayer containing a mixture of POPC and POPG lipids (3:1). Further, the effect of the interaction impairing mutations (V72A, I76A and F68A) on the conformational properties of the KCNE3 TMD in lipid bilayers was investigated. Our MD simulation results suggest that the KCNE3 TMD adopts a nearly linear α helical structural conformation in POPC-POPG lipid bilayers. Additionally, the results showed no significant change in the nearly linear α-helical conformation of KCNE3 TMD in the presence of interaction impairing mutations within the sampled time frame. The KCNE3 TMD is more stable with lower flexibility in comparison to the N-terminal and C-terminal of KCNE3 in lipid bilayers. The overall conformational flexibility of KCNE3 also varies in the presence of the interaction-impairing mutations. The MD simulation data further suggest that the membrane bilayer width is similar for wild-type KCNE3 and KCNE3 containing mutations. The Z-distance measurement data revealed that the TMD residue site A69 is close to the lipid bilayer center, and residue sites S57 and S82 are close to the surfaces of the lipid bilayer membrane for wild-type KCNE3 and KCNE3 containing interaction-impairing mutations. These results agree with earlier KCNE3 biophysical studies. The results of these MD simulations will provide complementary data to the experimental outcomes of KCNE3 to help understand its conformational dynamic properties in a more native lipid bilayer environment.

Electron paramagnetic resonance spectroscopic characterization of the human KCNE3 protein in lipodisq nanoparticles for structural dynamics of membrane proteins

One of the major challenges in solubilization of membrane proteins is to find the optimal physiological environment for their biophysical studies. EPR spectroscopy is a powerful biophysical technique for studying the structural and dynamic properties of macromolecules. However, the challenges in the membrane protein sample preparation and flexible motion of the spin label limit the utilization of EPR spectroscopy to a majority of membrane protein systems in a physiological membrane-bound state. Recently, lipodisq nanoparticles or styrene-maleic acid copolymer-lipid nanoparticles (SMALPs) have emerged as a membrane mimetic system for investigating the structural studies of membrane proteins. However, its detail characterization for membrane protein studies is still poorly understood. Recently, we characterized the potassium channel membrane protein KCNQ1 voltage sensing domain (KCNQ1-VSD) and KCNE1 reconstituted into lipodisq nanoparticles using EPR spectroscopy. In this study, the potassium channel accessory protein KCNE3 containing flexible N- and C-termini was encapsulated into proteoliposomes and lipodisq nanoparticles and characterized for studying its structural and dynamic properties using nitroxide based site-directed spin labeling EPR spectroscopy. CW-EPR lineshape analysis data indicated an increase in spectral line broadenings with the addition of the styrene-maleic acid (SMA) polymer which approaches close to the rigid limit providing a homogeneous stabilization of the protein-lipid complex. Similarly, EPR DEER measurements indicated an enhanced quality of distance measurements with an increase in the phase memory time (Tm) values upon incorporation of the sample into lipodisq nanoparticles, when compared to proteoliposomes. These results agree with the solution NMR structural structure of the KCNE3 and EPR studies of other membrane proteins in lipodisq nanoparticles. This study along with our earlier studies will provide the reference characterization data that will provide benefit to the membrane protein researchers for studying structural dynamics of challenging membrane proteins.

Functional Regulation of Aquaporin Dynamics by Lipid Bilayer Composition

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics were examined. We demonstrate that SoPIP2;1s structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

PEO-b-PBD Diblock Copolymers Induce Packing Defects in Lipid/Hybrid Membranes and Improve Insertion Rates of Natively Folded Peptides

Hybrid membranes assembled from biological lipids and synthetic polymers are a promising scaffold for the reconstitution and utilization of membrane proteins. Recent observations indicate that inclusion of small fractions of polymer in lipid membranes can improve protein folding and function, but the exact structural and physical changes a given polymer sequence imparts on a membrane often remain unclear. Here, we use all-atom molecular dynamics simulations to study the structure of hybrid membranes assembled from DOPC phospholipids and PEO-b-PBD diblock copolymers. We verified our computational model using new and existing experimental data and obtained a detailed picture of the polymer conformations in the lipid membrane that we can relate to changes in membrane elastic properties. We find that inclusion of low polymer fractions induces transient packing defects into the membrane. These packing defects act as insertion sites for two model peptides, and in this way, small amounts of polymer content in lipid membranes can lead to large increases in peptide insertion rates. Additionally, we report the peptide conformational space in both pure lipid and hybrid membranes. Both membranes support similar alpha helical peptide structures, exemplifying the biocompatibility of hybrid membranes.

Structure and dynamics of the SARS-CoV-2 envelope protein monomer

Coronaviruses, especially severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), present an ongoing threat to human wellbeing. Consequently, elucidation of molecular determinants of their function and interaction with the host is an important task. Whereas some of the coronaviral proteins are extensively characterized, others remain understudied. Here, we use molecular dynamics simulations to analyze the structure and dynamics of the SARS-CoV-2 envelope (E) protein (a viroporin) in the monomeric form. The protein consists of the hydrophobic α-helical transmembrane domain (TMD) and amphiphilic α-helices H2 and H3, connected by flexible linkers. We show that TMD has a preferable orientation in the membrane, while H2 and H3 reside at the membrane surface. Orientation of H2 is strongly influenced by palmitoylation of cysteines Cys40, Cys43, and Cys44. Glycosylation of Asn66 affects the orientation of H3. We also observe that the monomeric E protein both generates and senses the membrane curvature, preferably localizing with the C-terminus at the convex regions of the membrane; the protein in the pentameric form displays these properties as well. Localization to curved regions may be favorable for assembly of the E protein oligomers, whereas induction of curvature may facilitate the budding of the viral particles. The presented results may be helpful for a better understanding of the function of the coronaviral E protein and viroporins in general, and for overcoming the ongoing SARS-CoV-2 pandemic.

Spiers Memorial Lecture: Analysis and de novo design of membrane-interactive peptides

Membrane-peptide interactions play critical roles in many cellular and organismic functions, including protection from infection, remodeling of membranes, signaling, and ion transport. Peptides interact with membranes in a variety of ways: some associate with membrane surfaces in either intrinsically disordered conformations or well-defined secondary structures. Peptides with sufficient hydrophobicity can also insert vertically as transmembrane monomers, and many associate further into membrane-spanning helical bundles. Indeed, some peptides progress through each of these stages in the process of forming oligomeric bundles. In each case, the structure of the peptide and the membrane represent a delicate balance between peptide-membrane and peptide-peptide interactions. We will review this literature from the perspective of several biologically important systems, including antimicrobial peptides and their mimics, α-synuclein, receptor tyrosine kinases, and ion channels. We also discuss the use of de novo design to construct models to test our understanding of the underlying principles and to provide useful leads for pharmaceutical intervention of diseases.

Nanostructural Characterization of Cardiolipin-Containing Tethered Lipid Bilayers Adsorbed on Gold and Silicon Substrates for Protein Incorporation

A key to the development of lipid membrane-based devices is a fundamental understanding of how the molecular structure of the lipid bilayer membrane is influenced by the type of lipids used to build the membrane. This is particularly important when membrane proteins are included in these devices since the precise lipid environment affects the ability to incorporate membrane proteins and their functionality. Here, we used neutron reflectometry to investigate the structure of tethered bilayer lipid membranes and to characterize the incorporation of the NhaA sodium proton exchanger in the bilayer. The lipid membranes were composed of two lipids, dioleoyl phosphatidylcholine and cardiolipin, and were adsorbed on gold and silicon substrates using two different tethering architectures based on functionalized oligoethylene glycol molecules of different lengths. In all of the investigated samples, the addition of cardiolipin caused distinct structural rearrangement including crowding of ethylene glycol groups of the tethering molecules in the inner head region and a thinning of the lipid tail region. The incorporation of NhaA in the tethered bilayers following two different protocols is quantified, and the way protein incorporation modulates the structural properties of these membranes is detailed.

Extension of the CAVS model to the simulation of helical peptides in a membrane environment

Considering the effect of peptide insertion on the dipole potential of the lipid membrane, we extend the CAVS coarse-grained (CG) model to the simulation of helical peptides in a membrane environment. In this approach, the CG scheme for a peptide backbone is similar to the treatment in the united-atom model, while we treated the side chain of an amino acid by grouping 1-3 heavy atoms into a CG unit. The CAVS CG force field for peptides is optimized by reproducing the experimental results for the backbone (φ, ψ) distribution and predicting the PMF profiles of transferring organic molecules in a lipid bilayer membrane obtained from all-atom simulations. The CAVS simulation of a helical peptide in a phosphatidylcholine (PC) lipid bilayer revealed that the insertion of a peptide increases the dipole potential of the PC lipid bilayer, in which the peptide and its neutralized ions make a significant contribution. Finally, we carried out the CAVS simulation for five different helical peptides in the PC lipid bilayer to explore the behavior of peptide tilt, showing excellent agreement with the all-atom simulations. Our work suggests that the peptide tilt should relieve the deformation stress from the lipid bilayer, and the peptide aggregation could reduce the peptide tilt by resisting the deformation stress from the surrounding lipids.

The role of the envelope protein in the stability of a coronavirus model membrane against an ethanolic disinfectant

Ethanol is highly effective against various enveloped viruses and can disable the virus by disintegrating the protective envelope surrounding it. The interactions between the coronavirus envelope (E) protein and its membrane environment play key roles in the stability and function of the viral envelope. By using molecular dynamics simulation, we explore the underlying mechanism of ethanol-induced disruption of a model coronavirus membrane and, in detail, interactions of the E-protein and lipids. We model the membrane bilayer as N-palmitoyl-sphingomyelin and 1-palmitoyl-2-oleoylphosphatidylcholine lipids and the coronavirus E-protein. The study reveals that ethanol causes an increase in the lateral area of the bilayer along with thinning of the bilayer membrane and orientational disordering of lipid tails. Ethanol resides at the head-tail region of the membrane and enhances bilayer permeability. We found an envelope-protein-mediated increase in the ordering of lipid tails. Our simulations also provide important insights into the orientation of the envelope protein in a model membrane environment. At ∼25 mol. % of ethanol in the surrounding ethanol-water phase, we observe disintegration of the lipid bilayer and dislocation of the E-protein from the membrane environment.

Improved Protocol to Tackle the pH Effects on Membrane-Inserting Peptides

Many important biological pathways rely on membrane-interacting peptides or proteins, which can alter the biophysical properties of the cell membrane by simply adsorbing to its surface to undergo a full insertion process. To study these phenomena with atomistic detail, model peptides have been used to refine the current computational methodologies. Improvements have been made with force-field parameters, enhanced sampling techniques to obtain faster sampling, and the addition of chemical-physical properties, such as pH, whose influence dramatically increases at the water/membrane interface. The pH (low) insertion peptide (pHLIP) is a peptide that inserts across a membrane bilayer depending on the pH due to the presence of a key residue (Asp14) whose acidity-induced protonation triggers the whole process. The complex nature of these peptide/membrane interactions resulted in sampling limitations of the protonation and configurational space albeit using state-of-the-art methods such as the constant-pH molecular dynamics. To address this issue and circumvent those limitations, new simulations were performed with our newly developed pH-replica exchange method using wild-type (wt)-pHLIP in different 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine membrane sizes. This technique provided enhanced sampling and allowed for the calculation of more complete Asp14 pK_a profiles. The conformational heterogeneity derived from strong electrostatic interactions between Asp14 and the lipid phosphate groups was identified as the source of most pK_a variability. In spite of these persistent and harder-to-equilibrate phosphate interactions, the pK_a values at deeper regions (6.0-6.2) still predicted the experimental pK of insertion (6.0) since the electrostatic perturbation decays as the residue inserts further into the membrane. We also observed that reducing the system size leads to membrane deformations where it increasingly loses the ability to accommodate the pHLIP-induced perturbations. This indicates that large membrane patches, such as 256 or even 352 lipids, are needed to obtain stable and more realistic pHLIP/membrane systems. These results strengthen our method pK_a predictive and analytical capabilities to study the intricate play of electrostatic effects of the peptide/membrane interface, granting confidence for future applications in similar systems.

Effect of pH and temperature on the aggregation behaviour of dirhamnolipid biosurfactant. An experimental and molecular dynamics study

HYPOTHESIS

Pseudomonas aeruginosa dirhamnolipid (diRL) has been shown to form aggregates of different size and structure, under various conditions. Due to the presence of a carboxyl group in the molecule, it is expected that pH would strongly affect this aggregation behaviour. In addition, preliminary observations of temperature-induced changes in the states of aggregation of diRL supported the need of further investigation.

EXPERIMENTS

A systematic experimental study, using differential scanning calorimetry (DSC), small-angle Xray diffraction (SAXD), and Fourier-transform infrared spectroscopy (FTIR), has been carried out to characterize pH and temperature driven changes in the aggregation behavior of diRL biosurfactant. Molecular dynamics (MD) simulations, supported by the experimental results, allowed depicting molecular details on formation of diRL membranes and other aggregated structures under various physicochemical conditions.

FINDINGS

DiRL could adopt fairly organized multilayered structures (membranes) at low pH and temperature, which became highly disordered upon increasing either of these parameters. The effect of pH on the gauche/all-trans conformer ratio of the diRL acyl chains was not of significance, whereas temperature-induced effects were observed. For the first time it is described that diRL underwent an endothermic thermotropic transition with T_c = 34 °C as observed by DSC, at pH 4.5 (protonated diRL), but not at pH 7.4 (unprotonated diRL). FTIR confirmed these findings, showing a significant additional disordering of the all-trans acyl chains upon increasing temperature around that same value in the protonated form, an effect not observed for the dissociated form of the biosurfactant. In addition, at pH 7.4, changing temperature did not modify the hydration state of the polar moiety of diRL, whereas at pH 4.5 a significant decrease in the hydration state around 34 °C took place. SAXD data showed that protonated diRL formed multilayered structures at 20 °C, which converted into poorly correlated layers at 50 °C. MD simulations supported these findings, showing that the membrane-like structures formed by protonated diRL at 20 °C became unstable at higher temperatures, tending to form other structures, which could be micelles or other type of layered structures, whereas the negatively charged form of diRL organized in micelle-type aggregates in the whole range of temperature under study.

Pressure control using stochastic cell rescaling

Molecular dynamics simulations require barostats to be performed at a constant pressure. The usual recipe is to employ the Berendsen barostat first, which displays a first-order volume relaxation efficient in equilibration but results in incorrect volume fluctuations, followed by a second-order or a Monte Carlo barostat for production runs. In this paper, we introduce stochastic cell rescaling, a first-order barostat that samples the correct volume fluctuations by including a suitable noise term. The algorithm is shown to report volume fluctuations compatible with the isobaric ensemble and its anisotropic variant is tested on a membrane simulation. Stochastic cell rescaling can be straightforwardly implemented in the existing codes and can be used effectively in both equilibration and production phases.

Characterization of the Human KCNQ1 Voltage Sensing Domain (VSD) in Lipodisq Nanoparticles for Electron Paramagnetic Resonance (EPR) Spectroscopic Studies of Membrane Proteins

Membrane proteins are responsible for conducting essential biological functions that are necessary for the survival of living organisms. In spite of their physiological importance, limited structural information is currently available as a result of challenges in applying biophysical techniques for studying these protein systems. Electron paramagnetic resonance (EPR) spectroscopy is a very powerful technique to study the structural and dynamic properties of membrane proteins. However, the application of EPR spectroscopy to membrane proteins in a native membrane-bound state is extremely challenging due to the complexity observed in inhomogeneity sample preparation and the dynamic motion of the spin label. Detergent micelles are very popular membrane mimetics for membrane proteins due to their smaller size and homogeneity, providing high-resolution structure analysis by solution NMR spectroscopy. However, it is important to test whether the protein structure in a micelle environment is the same as that of its membrane-bound state. Lipodisq nanoparticles or styrene-maleic acid copolymer-lipid nanoparticles (SMALPs) have been introduced as a potentially good membrane-mimetic system for structural studies of membrane proteins. Recently, we reported on the EPR characterization of the KCNE1 membrane protein having a single transmembrane incorporated into lipodisq nanoparticles. In this work, lipodisq nanoparticles were used as a membrane mimic system for probing the structural and dynamic properties of the more complicated membrane protein system human KCNQ1 voltage sensing domain (Q1-VSD) having four transmembrane helices using site-directed spin-labeling EPR spectroscopy. Characterization of spin-labeled Q1-VSD incorporated into lipodisq nanoparticles was carried out using CW-EPR spectral line shape analysis and pulsed EPR double-electron electron resonance (DEER) measurements. The CW-EPR spectra indicate an increase in spectral line broadening with the addition of the styrene-maleic acid (SMA) polymer which approaches close to the rigid limit providing a homogeneous stabilization of the protein-lipid complex. Similarly, EPR DEER measurements indicated a superior quality of distance measurement with an increase in the phase memory time (T_m) values upon incorporation of the sample into lipodisq nanoparticles when compared to proteoliposomes. These results are consistent with the solution NMR structural studies on the Q1-VSD. This study will be beneficial for researchers working on investigating the structural and dynamic properties of more complicated membrane protein systems using lipodisq nanoparticles.

Active S2168 and inactive S21IRS pinholin interact differently with the lipid bilayer: A 31P and 2H solid state NMR study

Pinholins are a family of lytic membrane proteins responsible for the lysis of the cytosolic membrane in host cells of double stranded DNA bacteriophages. Protein-lipid interactions have been shown to influence membrane protein topology as well as its function. This work investigated the interactions of pinholin with the phospholipid bilayer while in active and inactive confirmations to elucidate the different interactions the two forms have with the bilayer. Pinholin incorporated into deuterated DMPC-d₅₄ lipid bilayers, along with ³¹P and ²H solid state NMR (SS-NMR) spectroscopy were used to probe the protein-lipid interactions with the phosphorus head group at the surface of the bilayer while interactions with the ²H nuclei were used to study the hydrophobic core. A comparison of the ³¹P chemical shift anisotropy (CSA) values of the active S²¹68 pinholin and inactive S²¹IRS pinholin indicated stronger head group interactions for the pinholin in its active form when compared to that of the inactive form supporting the model of a partially externalized peripheral transmembrane domain (TMD) of the active S²¹68 instead of complete externalized TMD1 as suggested by Ahammad et al. JPC B 2019. The ²H quadrupolar splitting analysis showed a decrease in spectral width for both forms of the pinholin when compared to the empty bilayers at all temperatures. In this case the decrease in the spectral width of the inactive S²¹IRS form of the pinholin showed stronger interactions with the acyl chains of the bilayer. The presence of the inactive form's additional TMD within the membrane was supported by the loss of peak resolution observed in the ²H NMR spectra.

Insights Into Protein S-Palmitoylation in Synaptic Plasticity and Neurological Disorders: Potential and Limitations of Methods for Detection and Analysis

S-palmitoylation (S-PALM) is a lipid modification that involves the linkage of a fatty acid chain to cysteine residues of the substrate protein. This common posttranslational modification (PTM) is unique among other lipid modifications because of its reversibility. Hence, like phosphorylation or ubiquitination, it can act as a switch that modulates various important physiological pathways within the cell. Numerous studies revealed that S-PALM plays a crucial role in protein trafficking and function throughout the nervous system. Notably, the dynamic turnover of palmitate on proteins at the synapse may provide a key mechanism for rapidly changing synaptic strength. Indeed, palmitate cycling on postsynaptic density-95 (PSD-95), the major postsynaptic density protein at excitatory synapses, regulates the number of synaptic α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) and thus affects synaptic transmission. Accumulating evidence suggests a relationship between impairments in S-PALM and severe neurological disorders. Therefore, determining the precise levels of S-PALM may be essential for understanding the ways in which this PTM is regulated in the brain and controls synaptic dynamics. Protein S-PALM can be characterized using metabolic labeling methods and biochemical tools. Both approaches are discussed herein in the context of specific methods and their advantages and disadvantages. This review clearly shows progress in the field, which has led to the development of new, more sensitive techniques that enable the detection of palmitoylated proteins and allow predictions of potential palmitate binding sites. Unfortunately, one significant limitation of these approaches continues to be the inability to use them in living cells.

Effects of Hydrophilic Residues and Hydrophobic Length on Flip-Flop Promotion by Transmembrane Peptides

Peptide-induced phospholipid flip-flop (scrambling) was evaluated using transmembrane model peptides in which the central residue was substituted with various amino acid residues (sequence: Ac-GKK(LA) _nXW(LA) _nLKKA-CONH₂). Peptides with a strongly hydrophilic residue (X = Q, N, or H) had higher scramblase activity than that of other peptides, and the activity was also dependent on the length of the peptides. Peptides with a hydrophobic stretch of 17 residues showed high flip-promotion propensity, whereas those of 21 and 25 residues did not, suggesting that membrane thinning under negative mismatch conditions promotes the flipping. Interestingly, a hydrophobic stretch of 19 residues intensively promoted phospholipid scrambling and membrane leakage. The distinctive characteristics of the peptide were ascribed by long-term molecular dynamics simulation to the arrangement of central glutamine and terminal four lysine residues on the same side of the helix. The combination of simulated and experimental data enables understanding of the mechanisms by which transmembrane helices, and ultimately unidentified scramblases in biomembranes, cause lipid scrambling.

Cell penetrating peptide modulation of membrane biomechanics by Molecular dynamics

The efficacy of a pharmaceutical treatment is often countered by the inadequate membrane permeability, that prevents drugs from reaching their specific intracellular targets. Cell penetrating peptides (CPPs) are able to route across cells' membrane various types of cargo, including drugs and nanoparticles. However, CPPs internalization mechanisms are not yet fully understood and depend on a wide variety of aspects. In this contest, the entry of a CPP into the lipid bilayer might induce molecular conformational changes, including marked variations on membrane's mechanical properties. Understanding how the CPP does influence the mechanical properties of cells membrane is crucial to design, engineer and improve new and existing penetrating peptides. Here, all atom Molecular Dynamics (MD) simulations were used to investigate the interaction between different types of CPPs embedded in a lipid bilayer of dioleoyl phosphatidylcholine (DOPC). In a greater detail, we systematically highlighted how CPP properties are responsible for modulating the membrane bending modulus. Our findings highlighted the CPP hydropathy strongly correlated with penetration of water molecules in the lipid bilayer, thus supporting the hypothesis that the amount of water each CPP can route inside the membrane is modulated by the hydrophobic and hydrophilic character of the peptide. Water penetration promoted by CPPs leads to a local decrease of the lipid order, which emerges macroscopically as a reduction of the membrane bending modulus.

Molecular Modeling of Multidrug Properties of Resistance Nodulation Division (RND) Transporters

Efflux pumps of the resistance nodulation division (RND) superfamily are among the major contributors to intrinsic and acquired multidrug resistance in Gram-negative bacteria. Structural information on AcrAB-TolC and MexAB-OprM, major efflux pumps of Escherichia coli and Pseudomonas aeruginosa respectively, boosted intensive research aimed at understanding the molecular mechanisms ruling the active extrusion processes. In particular, several studies were devoted to the understanding of the determinants behind the extraordinary broad specificity of the RND transporters AcrB and MexB. In this chapter, we discuss the ever-growing role computational methods have been playing in deciphering key structural and dynamical features of these transporters and of their interaction with substrates and inhibitors. We further discuss and illustrate examples from our lab of how molecular docking, homology modeling, all-atom molecular dynamics simulations and in silico free energy estimations can all together give precious insights into the processes of recognition and extrusion of substrates, as well as on the possible inhibition strategies.

The role of hydrophobic matching on transmembrane helix packing in cells

Folding and packing of membrane proteins are highly influenced by the lipidic component of the membrane. Here, we explore how the hydrophobic mismatch (the difference between the hydrophobic span of a transmembrane protein region and the hydrophobic thickness of the lipid membrane around the protein) influences transmembrane helix packing in a cellular environment. Using a ToxRED assay in Escherichia coli and a Bimolecular Fluorescent Complementation approach in human-derived cells complemented by atomistic molecular dynamics simulations we analyzed the dimerization of Glycophorin A derived transmembrane segments. We concluded that, biological membranes can accommodate transmembrane homo-dimers with a wide range of hydrophobic lengths. Hydrophobic mismatch and its effects on dimerization are found to be considerably weaker than those previously observed in model membranes, or under in vitro conditions, indicating that biological membranes (particularly eukaryotic membranes) can adapt to structural deformations through compensatory mechanisms that emerge from their complex structure and composition to alleviate membrane stress. Results based on atomistic simulations support this view, as they revealed that Glycophorin A dimers remain stable, despite of poor hydrophobic match, using mechanisms based on dimer tilting or local membrane thickness perturbations. Furthermore, hetero-dimers with large length disparity between their monomers are also tolerated in cells, and the conclusions that one can draw are essentially similar to those found with homo-dimers. However, large differences between transmembrane helices length hinder the monomer/dimer equilibrium, confirming that, the hydrophobic mismatch has, nonetheless, biologically relevant effects on helix packing in vivo.

How Does the P7C3-Series of Neuroprotective Small Molecules Prevent Membrane Disruption?

Molecular dynamics (MD) simulations are conducted to suggest a mechanism of action for the aminopropyl dibromocarbazole derivative (P7C3) small molecule, which protects neurons from apoptotic cell death. At first, the influence of embedded Aβ₄₂ stacks on the structure of membrane is studied. Then, the effect of P7C3 molecules on the Aβ₄₂ fibril enriched membrane and Aβ₄₂ fibril depleted membrane (when Aβ₄₂ fibrils are originally dissolved in the aqueous phase) are evaluated. Also, the formation of an amyloid ion channel in the Aβ₄₂ enriched membrane is examined by calculating deuterium order parameter, density profile, and surface thickness. For Aβ₄₂ in the fully inserted state, ion channel-like structures are formed. The presence of P7C3 molecules in this case just postpones membrane destruction but could not prevent pore formation. In contrast, when both Aβ₄₂ and P7C3 molecules are embedded in the aqueous solution, the P7C3 molecules are self-assembled at membrane/ionic aqueous solution interface and prevent the precipitation and deposition of Aβ₄₂ fibrils into the membrane.

The molecular behavior of a single β-amyloid inside a dipalmitoylphosphatidylcholine bilayer at three different temperatures: An atomistic simulation study: Aβ interaction with DPPC: Atomistic simulation

The behavior of a single Aβ40 molecule within a dipalmitoylphosphatidylcholine (DPPC) bilayer was studied by all-atom molecular dynamics simulations. The effect of membrane structure was investigated on Aβ40 behavior, secondary structure, and insertion depth. Simulations were performed at three temperatures (323, 310, and 300 K) to probe three different bilayer fluidities. Results show that at all above temperatures, the peptide contains two short helices, coil, bend, and turn structures. At 300 K, the peptide contains a region with β structure in C-terminal region. Our results also show that Aβ decreases the bilayer thickness and the order of lipids in its vicinity which leads to water insertion into the bilayer and concomitant increase in the local fluidity. The peptide remains embedded in the bilayer at all temperatures, and become inserted into the bilayer up to several residues at 323 and 310 K. At 310 and 300 K, the dominant interaction energy between Aβ and bilayer changes from electrostatic to van der Waals. It can be proposed that at higher temperatures (e.g., 323 K), Lys28 and the C-terminal region of the peptide play the role of two anchors that keep Aβ inside the top leaflet. This study demonstrates that Aβ molecule can perturb the integrity of cellular membranes. Proteins 2017; 85:1298-1310. © 2017 Wiley Periodicals, Inc.

Localization and Ordering of Lipids Around Aquaporin-0: Protein and Lipid Mobility Effects

Hydrophobic matching, lipid sorting, and protein oligomerization are key principles by which lipids and proteins organize in biological membranes. The Aquaporin-0 channel (AQP0), solved by electron crystallography (EC) at cryogenic temperatures, is one of the few protein-lipid complexes of which the structure is available in atomic detail. EC and room-temperature molecular dynamics (MD) of dimyristoylglycerophosphocholine (DMPC) annular lipids around AQP0 show similarities, however, crystal-packing and temperature might affect the protein surface or the lipids distribution. To understand the role of temperature, lipid phase, and protein mobility in the localization and ordering of AQP0-lipids, we used MD simulations of an AQP0-DMPC bilayer system. Simulations were performed at physiological and at DMPC gel-phase temperatures. To decouple the protein and lipid mobility effects, we induced gel-phase in the lipids or restrained the protein. We monitored the lipid ordering effects around the protein. Reducing the system temperature or inducing lipid gel-phase had a marginal effect on the annular lipid localization. However, restraining the protein mobility increased the annular lipid localization around the whole AQP0 surface, resembling EC. The distribution of the inter-phosphate and hydrophobic thicknesses showed that stretching of the DMPC annular layer around AQP0 surface is the mechanism that compensates the hydrophobic mismatch in this system. The distribution of the local area-per-lipid and the acyl-chain order parameters showed particular fluid- and gel-like areas that involved several lipid layers. These areas were in contact with the surfaces of higher and lower protein mobility, respectively. We conclude that the AQP0 surfaces induce specific fluid- and gel-phase prone areas. The presence of these areas might guide the AQP0 lipid sorting interactions with other membrane components, and is compatible with the squared array oligomerization of AQP0 tetramers separated by a layer of annular lipids.

Heterogeneous dielectric generalized Born model with a van der Waals term provides improved association energetics of membrane-embedded transmembrane helices

The heterogeneous dielectric generalized Born (HDGB) implicit membrane formalism is extended by the addition of a van der Waals dispersion term to better describe the nonpolar components of the free energy of solvation. The new model, termed HDGBvdW, improves the energy estimates in the hydrophobic interior of the membrane, where polar and charged species are rarely found and nonpolar interactions become significant. The implicit van der Waals term for the membrane environment extends the model from Gallicchio et al. (J. Comput. Chem. 2004, 25, 479) by combining separate contributions from each of the membrane components. The HDGBvdW model is validated with a series of test cases ranging from membrane insertion and pair association profiles of amino acid side chain analogs and transmembrane helices. Overall, the HDGBvdW model leads to increased agreement with explicit membrane simulation results and experimental data. © 2016 Wiley Periodicals, Inc.

Crystal Structure of the Human Cannabinoid Receptor CB1

Cannabinoid receptor 1 (CB₁) is the principal target of Δ⁹-tetrahydrocannabinol (THC), a psychoactive chemical from Cannabis sativa with a wide range of therapeutic applications and a long history of recreational use. CB₁ is activated by endocannabinoids and is a promising therapeutic target for pain management, inflammation, obesity, and substance abuse disorders. Here, we present the 2.8 Å crystal structure of human CB₁ in complex with AM6538, a stabilizing antagonist, synthesized and characterized for this structural study. The structure of the CB₁-AM6538 complex reveals key features of the receptor and critical interactions for antagonist binding. In combination with functional studies and molecular modeling, the structure provides insight into the binding mode of naturally occurring CB₁ ligands, such as THC, and synthetic cannabinoids. This enhances our understanding of the molecular basis for the physiological functions of CB₁ and provides new opportunities for the design of next-generation CB₁-targeting pharmaceuticals.

Folding and insertion thermodynamics of the transmembrane WALP peptide

The anchor of most integral membrane proteins consists of one or several helices spanning the lipid bilayer. The WALP peptide, GWW(LA)n (L)WWA, is a common model helix to study the fundamentals of protein insertion and folding, as well as helix-helix association in the membrane. Its structural properties have been illuminated in a large number of experimental and simulation studies. In this combined coarse-grained and atomistic simulation study, we probe the thermodynamics of a single WALP peptide, focusing on both the insertion across the water-membrane interface, as well as folding in both water and a membrane. The potential of mean force characterizing the peptide's insertion into the membrane shows qualitatively similar behavior across peptides and three force fields. However, the Martini force field exhibits a pronounced secondary minimum for an adsorbed interfacial state, which may even become the global minimum-in contrast to both atomistic simulations and the alternative PLUM force field. Even though the two coarse-grained models reproduce the free energy of insertion of individual amino acids side chains, they both underestimate its corresponding value for the full peptide (as compared with atomistic simulations), hinting at cooperative physics beyond the residue level. Folding of WALP in the two environments indicates the helix as the most stable structure, though with different relative stabilities and chain-length dependence.

Identification of dual ligands targeting angiotensin II type 1 receptor and peroxisome proliferator-activated receptor-γ by core hopping of telmisartan

It has been reported previously that some angiotensin II receptor blockers not only antagonize angiotensin II type 1 receptor (AT₁R), but also exert stimulation in peroxisome proliferator-activated receptor γ (PPARγ) partial activation, among which telmisartan displays the best. Telmisartan has been tested as a bifunctional ligand with antihypertensive and hypoglycemic activity. Aiming at more potent leads with selective AT₁R antagonism and PPARγ partial agonism, the three parts of telmisartan including the distal benzimidazole ring, the biphenyl moiety, and the carboxylic acid group experienced modification by core hopping method in our study. The central benzimidazole ring, however, remained intact considering its great affinity toward AT₁R and PPARγ. We utilized computational techniques for the sake of details on the binding interactions and conformational stability. Standard precision docking analysis and absorption, distribution, metabolism, excretion, and toxicity prediction received 10 molecules with higher Glide scores, similar interactions, and improved pharmacokinetic profiles compared to telmisartan. Comp#91 with highest scores for AT₁R (-11.92 kcal/mol) and PPARγ (-13.88 kcal/mol) exhibited excellent binding modes and pharmacokinetic parameters. Molecular dynamics trajectories on best docking pose of comp#91 confirmed the docking results and verified the conformational stability with both receptors throughout the course of 20-ns simulations. Thus, comp#91 could be identified as a promising lead in the development of dual AT₁R antagonist and PPARγ partial agonist against hypertension and type 2 diabetes.

Modeling Interactions between Multicomponent Vesicles and Antimicrobial Peptide-Inspired Nanoparticles

We examine the interaction between peptide-inspired nanoparticles, or nanopins, and multicomponent vesicles using the dissipative particle dynamics simulation technique. We study the role of nanopin architecture and cholesterol concentration on the binding of the nanopins to the lipid bilayer, their insertion, and postembedding self-organization. We find the insertion to be triggered by enthalpically unfavorable interactions between the hydrophilic solvent and the lipophilic components of the nanopins. The nanopins are observed to form aggregates in solution, insert into the bilayer, and disassemble into the individual nanopins following the insertion process. We examine factors that influence the orientation of the nanopins in the host vesicle. We report the length of the hydrophilic segment of the nanopins to regulate their orientation within the clusters before the embedding process and in the bilayer, after the postinsertion disassembly of the aggregates. The orientation angle distribution for a given nanopin architecture is found to be driven by energy minimization. In addition, higher concentration of cholesterol is observed to constrain the orientation of the nanopins. We also report thermal fluctuations to induce transverse diffusion of nanopins with specific architectures. The incidence of transverse diffusion is observed to decrease with the concentration of cholesterol. Our results can provide guidelines for designing peptide-inspired nanoparticles or macromolecules that can interface with living cells to serve as sensors for applications in medicine, sustainability, and energy.

NMR Structural Studies of Antimicrobial Peptides: LPcin Analogs

Lactophoricin (LPcin), a component of proteose peptone (113-135) isolated from bovine milk, is a cationic amphipathic antimicrobial peptide consisting of 23 amino acids. We designed a series of N- or C-terminal truncated variants, mutated analogs, and truncated mutated analogs using peptide-engineering techniques. Then, we selected three LPcin analogs of LPcin-C8 (LPcin-YK1), LPcin-T2WT6W (LPcin-YK2), and LPcin-T2WT6W-C8 (LPcin-YK3), which may have better antimicrobial activities than LPcin, and successfully expressed them in E. coli with high yield. We elucidated the 3D structures and topologies of the three LPcin analogs in membrane environments by conducting NMR structural studies. We investigated the purity of the LPcin analogs and the α-helical secondary structures by performing (1)H-(15)N 2D HSQC and HMQC-NOESY liquid-state NMR spectroscopy using protein-containing micelle samples. We measured the 3D structures and tilt angles in membranes by conducting (15)N 1D and 2D (1)H-(15)N SAMMY type solid-state NMR spectroscopy with an 800 MHz in-house-built (1)H-(15)N double-resonance solid-state NMR probe with a strip-shield coil, using protein-containing large bicelle samples aligned and confirmed by molecular-dynamics simulations. The three LPcin analogs were found to be curved α-helical structures, with tilt angles of 55-75° for normal membrane bilayers, and their enhanced activities may be correlated with these topologies.

The drug binding sites and transport mechanism of the RND pumps from Mycobacterium tuberculosis: Insights from molecular dynamics simulations

RND permease superfamily drug efflux pumps are involved in multidrug transport and are attractive to study them for therapeutic purpose. In previous work we have classified 14 members of MmpL proteins belong to RND superfamily from Mycobacterium tuberculosis (Mtb) within its families [Sandhu P. and Akhter Y., 2015. Int. J. Med. Microbiol., 305:413-423]. In this study, structures of these proteins are homology modelled. The drug binding sites and channels are identified using local micro-stereochemistry and charge densities. Potential transport mechanism based on differential structural behaviour in the absence and on the binding of drug molecules is explained using the molecular dynamics simulation results. Our studies show two potential drug binding sites positioned at opposite ends of the transport tunnel leading from cytoplasmic to the periplasmic space across MmpL5 trimer. The drug binding have effects on the structural conformation of the protein leading to molecular-scale peristaltic movements. The free binding energy calculations reveal that the subsequent binding events are interdependent and may have implications on transport mechanism. Two drug binding sites and a continuous channel in the RND pump have been reported. The proposed ligand binding mechanism shows peristaltic movements in the channel leading to the drug efflux. This study would be helpful in understanding the molecular basis of drugs resistance in Mtb.

Protein-fluctuation-induced water-pore formation in ion channel voltage-sensor translocation across a lipid bilayer membrane

We have applied a combined fluorescence microscopy and single-ion-channel electric current recording approach, correlating with molecular dynamics (MD) simulations, to study the mechanism of voltage-sensor domain translocation across a lipid bilayer. We use the colicin Ia ion channel as a model system, and our experimental and simulation results show the following: (1) The open-close activity of an activated colicin Ia is not necessarily sensitive to the amplitude of the applied cross-membrane voltage when the cross-membrane voltage is around the resting potential of excitable membranes; and (2) there is a significant probability that the activation of colicin Ia occurs by forming a transient and fluctuating water pore of ∼15 Å diameter in the lipid bilayer membrane. The location of the water-pore formation is nonrandom and highly specific, right at the insertion site of colicin Ia charged residues in the lipid bilayer membrane, and the formation is intrinsically associated with the polypeptide conformational fluctuations and solvation dynamics. Our results suggest an interesting mechanistic pathway for voltage-sensitive ion channel activation, and specifically for translocation of charged polypeptide chains across the lipid membrane under a transmembrane electric field: the charged polypeptide domain facilitates the formation of hydrophilic water pore in the membrane and diffuses through the hydrophilic pathway across the membrane; i.e., the charged polypeptide chain can cross a lipid membrane without entering into the hydrophobic core of the lipid membrane but entirely through the aqueous and hydrophilic environment to achieve a cross-membrane translocation. This mechanism sheds light on the intensive and fundamental debate on how a hydrophilic and charged peptide domain diffuses across the biologically inaccessible high-energy barrier of the hydrophobic core of a lipid bilayer: The peptide domain does not need to cross the hydrophobic core to move across a lipid bilayer.

Potential of mean force analysis of the self-association of leucine-rich transmembrane α-helices: difference between atomistic and coarse-grained simulations

Interaction of transmembrane (TM) proteins is important in many biological processes. Large-scale computational studies using coarse-grained (CG) simulations are becoming popular. However, most CG model parameters have not fully been calibrated with respect to lateral interactions of TM peptide segments. Here, we compare the potential of mean forces (PMFs) of dimerization of TM helices obtained using a MARTINI CG model and an atomistic (AT) Berger lipids-OPLS/AA model (AT(OPLS)). For helical, tryptophan-flanked, leucine-rich peptides (WL15 and WALP15) embedded in a parallel configuration in an octane slab, the AT(OPLS) PMF profiles showed a shallow minimum (with a depth of approximately 3 kJ/mol; i.e., a weak tendency to dimerize). A similar analysis using the CHARMM36 all-atom model (AT(CHARMM)) showed comparable results. In contrast, the CG analysis generally showed steep PMF curves with depths of approximately 16-22 kJ/mol, suggesting a stronger tendency to dimerize compared to the AT model. This CG > AT discrepancy in the propensity for dimerization was also seen for dilauroylphosphatidylcholine (DLPC)-embedded peptides. For a WL15 (and WALP15)/DLPC bilayer system, AT(OPLS) PMF showed a repulsive mean force for a wide range of interhelical distances, in contrast to the attractive forces observed in the octane system. The change from the octane slab to the DLPC bilayer also mitigated the dimerization propensity in the CG system. The dimerization energies of CG (AALALAA)3 peptides in DLPC and dioleoylphosphatidylcholine bilayers were in good agreement with previous experimental data. The lipid headgroup, but not the length of the lipid tails, was a key causative factor contributing to the differences between octane and DLPC. Furthermore, the CG model, but not the AT model, showed high sensitivity to changes in amino acid residues located near the lipid-water interface and hydrophobic mismatch between the peptides and membrane. These findings may help interpret CG and AT simulation results on membrane proteins.

Probing Structural Dynamics and Topology of the KCNE1 Membrane Protein in Lipid Bilayers via Site-Directed Spin Labeling and Electron Paramagnetic Resonance Spectroscopy

KCNE1 is a single transmembrane protein that modulates the function of voltage-gated potassium channels, including KCNQ1. Hereditary mutations in the genes encoding either protein can result in diseases such as congenital deafness, long QT syndrome, ventricular tachyarrhythmia, syncope, and sudden cardiac death. Despite the biological significance of KCNE1, the structure and dynamic properties of its physiologically relevant native membrane-bound state are not fully understood. In this study, the structural dynamics and topology of KCNE1 in bilayered lipid vesicles was investigated using site-directed spin labeling (SDSL) and electron paramagnetic resonance (EPR) spectroscopy. A 53-residue nitroxide EPR scan of the KCNE1 protein sequence including all 27 residues of the transmembrane domain (45-71) and 26 residues of the N- and C-termini of KCNE1 in lipid bilayered vesicles was analyzed in terms of nitroxide side-chain motion. Continuous wave-EPR spectral line shape analysis indicated the nitroxide spin label side-chains located in the KCNE1 TMD are less mobile when compared to the extracellular region of KCNE1. The EPR data also revealed that the C-terminus of KCNE1 is more mobile when compared to the N-terminus. EPR power saturation experiments were performed on 41 sites including 18 residues previously proposed to reside in the transmembrane domain (TMD) and 23 residues of the N- and C-termini to determine the topology of KCNE1 with respect to the 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (POPG) lipid bilayers. The results indicated that the transmembrane domain is indeed buried within the membrane, spanning the width of the lipid bilayer. Power saturation data also revealed that the extracellular region of KCNE1 is solvent-exposed with some of the portions partially or weakly interacting with the membrane surface. These results are consistent with the previously published solution NMR structure of KCNE1 in micelles.