Interactions of the M2delta segment of the acetylcholine receptor with lipid bilayers: a continuum-solvent model study

Investigating the Secondary Structure of Membrane Peptides Utilizing Multiple 2H-Labeled Hydrophobic Amino Acids via Electron Spin Echo Envelope Modulation (ESEEM) Spectroscopy

An electron spin echo envelope modulation (ESEEM) approach was used to probe local secondary structures of membrane proteins and peptides. This ESEEM method detects dipolar couplings between ²H-labeled nuclei on the side chains of an amino acid (Leu or Val) and a strategically placed nitroxide spin-label in the proximity up to 8 Å. ESEEM spectra patterns for different samples correlate directly to the periodic structural feature of different secondary structures. Since this pattern can be affected by the side chain length and flexibility of the ²H-labeled amino acid used in the experiment, it is important to examine several different hydrophobic amino acids (d₃ Ala, d₈ Val, d₈ Phe) utilizing this ESEEM approach. In this work, a series of ESEEM data were collected on the AChR M2δ membrane peptide to build a reference for the future application of this approach for various biological systems. The results indicate that, despite the relative intensity and signal-to-noise level, all amino acids share a similar ESEEM modulation pattern for α-helical structures. Thus, all commercially available ²H-labeled hydrophobic amino acids can be utilized as probes for the further application of this ESEEM approach. Also, the ESEEM signal intensities increase as the side chain length gets longer or less rigid. In addition, longer side chain amino acids had a larger ²H ESEEM FT peak centered at the ²H Larmor frequency for the i ± 4 sample when compared to the corresponding i ± 3 sample. For shorter side chain amino acids, the ²H ESEEM FT peak intensity ratio between i ± 4 and i ± 3 was not well-defined.

Probing the Secondary Structure of Membrane Peptides Using (2)H-Labeled d(10)-Leucine via Site-Directed Spin-Labeling and Electron Spin Echo Envelope Modulation Spectroscopy

Previously, we reported an electron spin echo envelope modulation (ESEEM) spectroscopic approach for probing the local secondary structure of membrane proteins and peptides utilizing (2)H isotopic labeling and site-directed spin-labeling (SDSL). In order to probe the secondary structure of a peptide sequence, an amino acid residue (i) side chain was (2)H-labeled, such as (2)H-labeled d10-Leucine, and a cysteine residue was strategically placed at a subsequent nearby position (denoted as i + 1 to i + 4) to which a nitroxide spin label was attached. In order to fully access and demonstrate the feasibility of this new ESEEM approach with (2)H-labeled d10-Leu, four Leu residues within the AChR M2δ peptide were fully mapped out using this ESEEM method. Unique (2)H-ESEEM patterns were observed with the (2)H-labeled d10-Leu for the AChR M2δ α-helical model peptide. For proteins and peptides with an α-helical secondary structure, deuterium modulation can be clearly observed for i ± 3 and i ± 4 samples, but not for i ± 2 samples. Also, a deuterium peak centered at the (2)H Larmor frequency of each i ± 4 sample always had a significantly higher intensity than the corresponding i + 3 sample. This unique feature can be potentially used to distinguish an α-helix from a π-helix or 310-helix. Moreover, (2)H modulation depth for ESEEM samples on Leu10 were significantly enhanced which was consistent with a kinked or curved structural model of the AChR M2δ peptide as suggested by previous MD simulations and NMR experiments.

Refining the treatment of membrane proteins by coarse-grained models

Obtaining a quantitative description of the membrane proteins stability is crucial for understanding many biological processes. However the advance in this direction has remained a major challenge for both experimental studies and molecular modeling. One of the possible directions is the use of coarse-grained models but such models must be carefully calibrated and validated. Here we use a recent progress in benchmark studies on the energetics of amino acid residue and peptide membrane insertion and membrane protein stability in refining our previously developed coarse-grained model (Vicatos et al., Proteins 2014;82:1168). Our refined model parameters were fitted and/or tested to reproduce water/membrane partitioning energetics of amino acid side chains and a couple of model peptides. This new model provides a reasonable agreement with experiment for absolute folding free energies of several β-barrel membrane proteins as well as effects of point mutations on a relative stability for one of those proteins, OmpLA. The consideration and ranking of different rotameric states for a mutated residue was found to be essential to achieve satisfactory agreement with the reference data.

Solid-state NMR (31)P paramagnetic relaxation enhancement membrane protein immersion depth measurements

Paramagnetic relaxation enhancement (PRE) is a widely used approach for measuring long-range distance constraints in biomolecular solution NMR spectroscopy. In this paper, we show that (31)P PRE solid-state NMR spectroscopy can be utilized to determine the immersion depth of spin-labeled membrane peptides and proteins. Changes in the (31)P NMR PRE times coupled with modeling studies can be used to describe the spin-label position/amino acid within the lipid bilayer and the corresponding helical tilt. This method provides valuable insight on protein-lipid interactions and membrane protein structural topology. Solid-state (31)P NMR data on the 23 amino acid α-helical nicotinic acetylcholine receptor nAChR M2δ transmembrane domain model peptide followed predicted behavior of (31)P PRE rates of the phospholipid headgroup as the spin-label moves from the membrane surface toward the center of the membrane. Residue 11 showed the smallest changes in (31)P PRE (center of the membrane), while residue 22 shows the largest (31)P PRE change (near the membrane surface), when compared to the diamagnetic control M2δ sample. This PRE SS-NMR technique can be used as a molecular ruler to measure membrane immersion depth.

Time-resolved EPR immersion depth studies of a transmembrane peptide incorporated into bicelles

The reduction in EPR signal intensity of nitroxide spin-labels by ascorbic acid has been measured as a function of time to investigate the immersion depth of the spin-labeled M2δ AChR peptide incorporated into a bicelle system utilizing EPR spectroscopy. The corresponding decay curves of n-DSA (n=5, 7, 12, and 16) EPR signals have been used to (1) calibrate the depth of the bicelle membrane and (2) establish a calibration curve for measuring the depth of spin-labeled transmembrane peptides. The kinetic EPR data of CLS, n-DSA (n=5, 7, 12, and 16), and M2δ AChR peptide spin-labeled at Glu-1 and Ala-12 revealed excellent exponential and linear fits. For a model M2δ AChR peptide, the depth of immersion was calculated to be 5.8Å and 3Å for Glu-1, and 21.7Å and 19Å for Ala-12 in the gel-phase (298K) and L(α)-phases (318K), respectively. The immersion depth values are consistent with the pitch of an α-helix and the structural model of M2δ AChR incorporated into the bicelle system is in a good agreement with previous studies. Therefore, this EPR time-resolved kinetic technique provides a new reliable method to determine the immersion depth of membrane-bound peptides, as well as, explore the structural characteristics of the M2δ AChR peptide.

Capturing spontaneous partitioning of peripheral proteins using a biphasic membrane-mimetic model

Membrane binding of peripheral proteins, mediated by specialized anchoring domains, is a crucial step for their biological function. Computational studies of membrane insertion, however, have proven challenging and largely inaccessible, due to the time scales required for the complete description of the process, mainly caused by the slow diffusion of the lipid molecules composing the membrane. Furthermore, in many cases, the nature of the membrane "anchor", i.e., the part of the protein that inserts into the membrane, is also unknown. Here, we address some of these issues by developing and employing a simplified representation of the membrane by a biphasic solvent model which we demonstrate can be used efficiently to capture and describe the process of hydrophobic insertion of membrane anchoring domains in all-atom molecular dynamics simulations. Applying the model, we have studied the insertion of the anchoring domain of a coagulation protein (the GLA domain of human protein C), starting from multiple initial configurations varying with regard to the initial orientation and height of the protein with respect to the membrane. In addition to efficiently and consistently identifying the "keel" region as the hydrophobic membrane anchor, within a few nanoseconds each configuration simulated showed a convergent height (2.20 ± 1.04 Å) and angle with respect to the interface normal (23.37 ± 12.48°). We demonstrate that the model can produce the same results as those obtained from a full representation of a membrane, in terms of both the depth of penetration and the orientation of the protein in the final membrane-bound form with an order of magnitude decrease in the required computational time compared to previous models, allowing for a more exhaustive search for the correct membrane-bound configuration.

Interaction of an antimicrobial peptide with membranes: experiments and simulations with NKCS

We used Monte Carlo simulations and biophysical measurements to study the interaction of NKCS, a derivative of the antimicrobial peptide NK-2, with a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) membrane. The simulations showed that NKCS adsorbed on the membrane surface and the dominant conformation featured two amphipathic helices connected by a hinge region. We designed two mutants in the hinge to investigate the interplay between helicity and membrane affinity. Simulations with a Leu-to-Pro substitution showed that the helicity and membrane affinity of the mutant (NKCS-[LP]) decreased. Two Ala residues were added to NKCS to produce a sequence that is compatible with a continuous amphipathic helix structure (NKCS-[AA]), and the simulations showed that the mutant adsorbed on the membrane surface with a particularly high affinity. The circular dichroism spectra of the three peptides also showed that NKCS-[LP] is the least helical and NKCS-[AA] is the most. However, the activity of the peptides, determined in terms of their antimicrobial potency and influence on the temperature of the transition of the lipid to hexagonal phase, displayed a complex behavior: NKCS-[LP] was the least potent and had the smallest influence on the transition temperature, and NKCS was the most potent and had the largest effect on the temperature.

A generalized born implicit-membrane representation compared to experimental insertion free energies

An implicit-membrane representation based on the generalized Born theory of solvation has been developed. The method was parameterized against the water-to-cyclohexane insertion free energies of hydrophobic side-chain analogs. Subsequently, the membrane was compared with experimental data from translocon inserted polypeptides and validated by comparison with an independent dataset of six membrane-associated peptides and eight integral membrane proteins of known structure and orientation. Comparison of the insertion energy of alpha-helical model peptides with the experimental values from the biological hydrophobicity scale of Hessa et al. gave a correlation of 93% with a mean unsigned error of 0.64 kcal/mol, when charged residues were ignored. The membrane insertion energy was found to be dependent on residue position. This effect is particularly pronounced for charged and polar residues, which strongly prefer interfacial locations. All integral membrane proteins investigated orient and insert correctly into the implicit-membrane model. Remarkably, the membrane model correctly predicts a partially inserted configuration for the monotopic membrane protein cyclooxygenase, matching experimental and theoretical predictions. To test the applicability and usefulness of the implicit-membrane method, molecular simulations of influenza A M2 as well as the glycophorin A dimer were performed. Both systems remain structurally stable and integrated into the membrane.

Properties of integral membrane protein structures: derivation of an implicit membrane potential

Distributions of each amino acid in the trans-membrane domain were calculated as a function of the membrane normal using all currently available alpha-helical membrane protein structures with resolutions better than 4 A. The results were compared with previous sequence- and structure-based analyses. Calculation of the average hydrophobicity along the membrane normal demonstrated that the protein surface in the membrane domain is in fact much more hydrophobic than the protein core. While hydrophobic residues dominate the membrane domain, the interfacial regions of membrane proteins were found to be abundant in the small residues glycine, alanine, and serine, consistent with previous studies on membrane protein packing. Charged residues displayed nonsymmetric distributions with a preference for the intracellular interface. This effect was more prominent for Arg and Lys resulting in a direct confirmation of the positive inside rule. Potentials of mean force along the membrane normal were derived for each amino acid by fitting Gaussian functions to the residue distributions. The individual potentials agree well with experimental and theoretical considerations. The resulting implicit membrane potential was tested on various membrane proteins as well as single trans-membrane alpha-helices. All membrane proteins were found to be at an energy minimum when correctly inserted into the membrane. For alpha-helices both interfacial (i.e. surface bound) and inserted configurations were found to correspond to energy minima. The results demonstrate that the use of trans-membrane amino acid distributions to derive an implicit membrane representation yields meaningful residue potentials.

Evaluating tilt angles of membrane-associated helices: comparison of computational and NMR techniques

A computational method to calculate the orientation of membrane-associated alpha-helices with respect to a lipid bilayer has been developed. It is based on a previously derived implicit membrane representation, which was parameterized using the structures of 46 alpha-helical membrane proteins. The method is validated by comparison with an independent data set of six transmembrane and nine antimicrobial peptides of known structure and orientation. The minimum energy orientations of the transmembrane helices were found to be in good agreement with tilt and rotation angles known from solid-state NMR experiments. Analysis of the free-energy landscape found two types of minima for transmembrane peptides: i), Surface-bound configurations with the helix long axis parallel to the membrane, and ii), inserted configurations with the helix spanning the membrane in a perpendicular orientation. In all cases the inserted configuration also contained the global energy minimum. Repeating the calculations with a set of solution NMR structures showed that the membrane model correctly distinguishes native transmembrane from nonnative conformers. All antimicrobial peptides investigated were found to orient parallel and bind to the membrane surface, in agreement with experimental data. In all cases insertion into the membrane entailed a significant free-energy penalty. An analysis of the contributions of the individual residue types confirmed that hydrophobic residues are the main driving force behind membrane protein insertion, whereas polar, charged, and aromatic residues were found to be important for the correct orientation of the helix inside the membrane.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.