Amino-acid solvation structure in transmembrane helices from molecular dynamics simulations

Characterizations of a novel peptide encoded by a circular RNA using in-silico analyses

CircRNAs have gained importance in recent times due to their involvement in gene regulation and also in the prognosis of cancer. Generally, the circRNA directly interact with miRNA or RNA binding proteins to exert their action, but some of them can be translated. These translated peptides often participate in the regulation of cellular processes. The circPPP1R12A translated peptide has been shown to influence the functioning of the Mst pathway. The Mst signaling is noteworthy for its role in the process of development, but it also has a function as a regulator of apoptosis, which is significant for regulation of cancer. Overexpression of this novel peptide deactivates the Mst signaling to induce the expression of the proliferative oncogene, Yap. Its molecular interaction with the molecules in the Mst pathway is hitherto unknown. In this short report we present our findings from in-silico studies the plausible structure of the peptide through bioinformatics and dynamics simulation studies. This is the first such report on the structure of the novel peptide encoded by circPPP1R12A, which could be important to predict in future its molecular interactions to understand its functionality.

Mapping temperature-dependent conformational change in the voltage-sensing domain of an engineered heat-activated K+ channel

Temperature-dependent regulation of ion channel activity is critical for a variety of physiological processes ranging from immune response to perception of noxious stimuli. Our understanding of the structural mechanisms that underlie temperature sensing remains limited, in part due to the difficulty of combining high-resolution structural analysis with temperature stimulus. Here, we use NMR to compare the temperature-dependent behavior of Shaker potassium channel voltage sensor domain (WT-VSD) to its engineered temperature sensitive (TS-VSD) variant. Further insight into the molecular basis for temperature-dependent behavior is obtained by analyzing the experimental results together with molecular dynamics simulations. Our studies reveal that the overall secondary structure of the engineered TS-VSD is identical to the wild-type channels except for local changes in backbone torsion angles near the site of substitution (V369S and F370S). Remarkably however, these structural differences result in increased hydration of the voltage-sensing arginines and the S4-S5 linker helix in the TS-VSD at higher temperatures, in contrast to the WT-VSD. These findings highlight how subtle differences in the primary structure can result in large-scale changes in solvation and thereby confer increased temperature-dependent activity beyond that predicted by linear summation of solvation energies of individual substituents.

Aromatic interactions with membrane modulate human BK channel activation

Large-conductance potassium (BK) channels are transmembrane (TM) proteins that can be synergistically and independently activated by membrane voltage and intracellular Ca²⁺. The only covalent connection between the cytosolic Ca²⁺ sensing domain and the TM pore and voltage sensing domains is a 15-residue ‘C-linker’. To determine the linker’s role in human BK activation, we designed a series of linker sequence scrambling mutants to suppress potential complex interplay of specific interactions with the rest of the protein. The results revealed a surprising sensitivity of BK activation to the linker sequence. Combining atomistic simulations and further mutagenesis experiments, we demonstrated that nonspecific interactions of the linker with membrane alone could directly modulate BK activation. The C-linker thus plays more direct roles in mediating allosteric coupling between BK domains than previously assumed. Our results suggest that covalent linkers could directly modulate TM protein function and should be considered an integral component of the sensing apparatus.

Dissecting the thermodynamic contributions of the charged residues in the membrane anchoring of Bcl-xl C-terminal domain

The C-terminal helix of the Bcl-xl is known to initiate the membrane insertion of the protein by anchoring into the mitochondrial outer membrane. The C-terminal charged residues of that helix, R232 and K233, are reported to have an important structural role in the process of that insertion. The present work provides a quantitative understanding of the thermodynamic contribution of these residues on the membrane insertion energy-profile, calculated from the Adaptive Biasing Force based MD simulations of 2.67 μs altogether. Interestingly, the effect of the single neutralizing mutations at the C-terminus, i.e. K233A or R232A, is easily tolerated by the peptide without impacting the nature of insertion energy-profile, indicating the efficiency of one positively charged residue to drive the insertion. Whereas a double mutant, i.e. R232A and K233A, makes a significant impact on the energy-profile by destabilizing the membrane-associated states, as well as the membrane-embedded states. The finding provides molecular-level mechanistic insight. The water-mediated interaction formed by the peptide polar side chains within the bilayer core is found to modulate the membrane response during peptide insertion and that subsequently regulates the insertion mechanism. Mutation of the C-terminal residues eventually alters such a cascade of interactions that results in an insertion through energetically more expensive pathway. Since any one of the positively charged residues at the terminal is critical to ensure the membrane insertion, it appears that the natural selection of 'two' instead of 'one' charged residue is redundant in the context of membrane anchoring but may be important for other biochemical events.

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.

Direct Measurement of the Effect of Cholesterol and 6-Ketocholestanol on the Membrane Dipole Electric Field Using Vibrational Stark Effect Spectroscopy Coupled with Molecular Dynamics Simulations

Biological membranes are heterogeneous structures with complex electrostatic profiles arising from lipids, sterols, membrane proteins, and water molecules. We investigated the effect of cholesterol and its derivative 6-ketocholestanol (6-kc) on membrane electrostatics by directly measuring the dipole electric field (F⃗d) within lipid bilayers containing cholesterol or 6-kc at concentrations of 0-40 mol% through the vibrational Stark effect (VSE). We found that adding low concentrations of cholesterol, up to ∼10 mol %, increases F⃗d, while adding more cholesterol up to 40 mol% lowers F⃗d. In contrast, we measured a monotonic increase in F⃗d as 6-kc concentration increased. We propose that this membrane electric field is affected by multiple factors: the polarity of the sterol molecules, the reorientation of the phospholipid dipole due to sterol, and the impact of the sterol on hydrogen bonding with surface water. We used molecular dynamics simulations to examine the distribution of phospholipids, sterol, and helix in bilayers containing these sterols. At low concentrations, we observed clustering of sterols near the vibrational probe whereas at high concentrations, we observed spatial correlation between the positions of the sterol molecules. This work demonstrates how a one-atom difference in a sterol changes the physicochemical and electric field properties of the bilayer.

Influence of glutamic acid residues and pH on the properties of transmembrane helices

Negatively charged side chains are important for the function of particular ion channels and certain other membrane proteins. To investigate the influence of single glutamic acid side chains on helices that span lipid-bilayer membranes, we have employed GWALP23 (acetyl-GGALW⁵LALALALALALALW¹⁹LAGA-amide) as a favorable host peptide framework. We substituted individual Leu residues with Glu residues (L12E or L14E or L16E) and incorporated specific ²H-labeled alanine residues within the core helical region or near the ends of the sequence. Solid-state ²H NMR spectra reveal little change for the core labels in GWALP23-E12, -E14 and -E16 over a pH range of 4 to 12.5, with the spectra being broader for samples in DOPC compared to DLPC bilayers. The spectra for samples with deuterium labels near the helix ends on alanines 3 and 21 show modest pH-dependent changes in the extent of unwinding of the helix terminals in DLPC and DOPC bilayers. The combined results indicate minor overall responses of these transmembrane helices to changes in pH, with the most buried residue E12 showing no pH dependence. While the Glu residues E14 and E16 may have high pK_a values in the lipid bilayer environment, it is also possible that a paucity of helix response is masking the pK_a values. Interestingly, when E16 is present, spectral changes at high pH report significant local unwinding of the core helix. Our results are consistent with the expectation that buried carboxyl groups aggressively hold their protons and/or waters of hydration.

The Single Transmembrane Segment of Minimal Sensor DesK Senses Temperature via a Membrane-Thickness Caliper

Thermosensors detect temperature changes and trigger cellular responses crucial for survival at different temperatures. The thermosensor DesK is a transmembrane (TM) histidine kinase which detects a decrease in temperature through its TM segments (TMS). Here, we address a key issue: how a physical stimulus such as temperature can be converted into a cellular response. We show that the thickness of Bacillus lipid membranes varies with temperature and that such variations can be detected by DesK with great precision. On the basis of genetic studies and measurements of in vitro activity of a DesK construct with a single TMS (minimal sensor DesK [MS-DesK]), reconstituted in liposomes, we propose an interplay mechanism directed by a conserved dyad, phenylalanine 8-lysine 10. This dyad is critical to anchor the only transmembrane segment of the MS-DesK construct to the extracellular water-lipid interphase and is required for the transmembrane segment of MS-DesK to function as a caliper for precise measurement of membrane thickness. The data suggest that positively charged lysine 10, which is located in the hydrophobic core of the membrane but is close to the water-lipid interface, pulls the transmembrane region toward the water phase to localize its charge at the interface. Nevertheless, the hydrophobic residue phenylalanine 8, located at the N-terminal extreme of the TMS, has a strong tendency to remain in the lipid phase, impairing access of lysine 10 to the water phase. The outcome of this interplay is a fine-tuned sensitivity to membrane thickness that elicits conformational changes that favor different signaling states of the protein.

IMPORTANCE

The ability to sense and respond to extracellular signals is essential for cell survival. One example is the cellular response to temperature variation. How do cells "sense" temperature changes? It has been proposed that the bacterial thermosensor DesK acts as a molecular caliper measuring membrane thickness variations that would occur as a consequence of temperature changes and activates a pathway to restore membrane fluidity at low temperature. Here, we demonstrated that membrane thickness variations do occur at physiological temperatures by directly measuring Bacillus lipid membrane thickness. We also dissected the N-terminal sensing motif of MS-DesK at the molecular-biophysical level and found that the dyad phenylalanine-lysine at the water-lipid phase is critical for achievement of a fine-tuned sensitivity to temperature.

Naratriptan aggregation in lipid bilayers: perspectives from molecular dynamics simulations

In order to understand the interaction between naratriptan and a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidyl-choline (POPC), we carried out molecular dynamics simulations. The simulations were performed considering neutral and protonated ionization states, starting from different initial conditions. At physiological pH, the protonated state of naratriptan is predominant. It is expected that neutral compounds could have larger membrane partition than charged compounds. However, for the specific case of triptans, it is difficult to study neutral species in membranes experimentally, making computer simulations an interesting tool. When the naratriptan molecules were originally placed in water, they partitioned between the bilayer/water interface and water phase, as has been described for similar compounds. From this condition, the drugs displayed low access to the hydrophobic environment, with no significant effects on bilayer organization. The molecules anchored in the interface, due mainly to the barrier function of the polar and oriented lipid heads. On the other hand, when placed inside the bilayer, both neutral and protonated naratriptan showed self-aggregation in the lipid tail environment. In particular, the protonated species exhibited a pore-like structure, dragging water through this environment. Graphical Abstract Different behaviour of Naratriptan and Sumatriptan, when the drugs were originally placed in the lipid core.

Novel pH-Sensitive Cyclic Peptides

A series of cyclic peptides containing a number of tryptophan (W) and glutamic acid (E) residues were synthesized and evaluated as pH-sensitive agents for targeting of acidic tissue and pH-dependent cytoplasmic delivery of molecules. Biophysical studies revealed the molecular mechanism of peptides action and localization within the lipid bilayer of the membrane at high and low pHs. The symmetric, c[(WE)₄WC], and asymmetric, c[E₄W₅C], cyclic peptides translocated amanitin, a polar cargo molecule of similar size, across the lipid bilayer and induced cell death in a pH- and concentration-dependent manner. Fluorescently-labelled peptides were evaluated for targeting of acidic 4T1 mammary tumors in mice. The highest tumor to muscle ratio (5.6) was established for asymmetric cyclic peptide, c[E₄W₅C], at 24 hours after intravenous administration. pH-insensitive cyclic peptide c[R₄W₅C], where glutamic acid residues (E) were replaced by positively charged arginine residues (R), did not exhibit tumor targeting. We have introduced a novel class of cyclic peptides, which can be utilized as a new pH-sensitive tool in investigation or targeting of acidic tissue.

C-terminal tail insertion of Bcl-xL in membrane occurs via partial unfolding and refolding cycle associating microsolvation

Insertion of the Bcl-x_L C-terminal into the mitochondrial outer membrane is found to be guided by the partial unfolding–refolding cycle, assisted by micro-solvation.

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.

Marginally hydrophobic transmembrane α-helices shaping membrane protein folding

Cells have developed an incredible machinery to facilitate the insertion of membrane proteins into the membrane. While we have a fairly good understanding of the mechanism and determinants of membrane integration, more data is needed to understand the insertion of membrane proteins with more complex insertion and folding pathways. This review will focus on marginally hydrophobic transmembrane helices and their influence on membrane protein folding. These weakly hydrophobic transmembrane segments are by themselves not recognized by the translocon and therefore rely on local sequence context for membrane integration. How can such segments reside within the membrane? We will discuss this in the light of features found in the protein itself as well as the environment it resides in. Several characteristics in proteins have been described to influence the insertion of marginally hydrophobic helices. Additionally, the influence of biological membranes is significant. To begin with, the actual cost for having polar groups within the membrane may not be as high as expected; the presence of proteins in the membrane as well as characteristics of some amino acids may enable a transmembrane helix to harbor a charged residue. The lipid environment has also been shown to directly influence the topology as well as membrane boundaries of transmembrane helices-implying a dynamic relationship between membrane proteins and their environment.

Measurement of the membrane dipole electric field in DMPC vesicles using vibrational shifts of p-cyanophenylalanine and molecular dynamics simulations

The magnitude of the membrane dipole field was measured using vibrational Stark effect (VSE) shifts of nitrile oscillators placed on the unnatural amino acid p-cyanophenylalanine (p-CN-Phe) added to a peptide sequence at four unique positions. These peptides, which were based on a repeating alanine-leucine motif, intercalated into small unilamellar DMPC vesicles which formed an α-helix as confirmed by circular dichroic (CD) spectroscopy. Molecular dynamics simulations of the membrane-intercalated helix containing two of the nitrile probes, one near the headgroup region of the lipid (αLAX(25)) and one buried in the interior of the bilayer (αLAX(16)), were used to examine the structure of the nitrile with respect to the membrane normal, the assumed direction of the dipole field, by quantifying both a small tilt of the helix in the bilayer and conformational rotation of the p-CN-Phe side chain at steady state. Vibrational absorption energies of the nitrile oscillator at each position showed a systematic blue shift as the nitrile was stepped toward the membrane interior; for several different concentrations of peptide, the absorption energy of the nitrile located in the middle of the bilayer was ∼3 cm(-1) greater than that of the nitrile closest to the surface of the membrane. Taken together, the measured VSE shifts and nitrile orientations within the membrane resulted in an absolute magnitude of 8-11 MV/cm for the dipole field, at the high end of the range of possible values that have been accumulated from a variety of indirect measurements. Implications for this are discussed.

Interaction of glycine, lysine, proline and histidine with dipalmitoylphosphatidylcholine lipid bilayers: a theoretical and experimental study

The interaction of unblocked glycine, lysine, proline and histidine with a DPPC bilayer was assessed using extensive atomistic molecular dynamics simulations.

Modeling the structure of SARS 3a transmembrane protein using a minimum unfavorable contact approach

3a is an accessory protein from SARS coronavirus that is known to play a significant role in the proliferation of the virus by forming tetrameric ion channels. Although the monomeric units are known to consist of three transmembrane (TM) domains, there are no solved structures available for the complete monomer. The present study proposes a structural model for the transmembrane region of the monomer by employing our previously tested approach, which predicts potential orientations of TM α-helices by minimizing the unfavorable contact surfaces between the different TM domains. The best model structure comprising all three α-helices has been subjected to MD simulations to examine its quality. The TM bundle was found to form a compact and stable structure with significant intermolecular interactions. The structural features of the proposed model of 3a account for observations from previous experimental investigations on the activity of the protein. Further analysis indicates that residues from the TM2 and TM3 domains are likely to line the pore of the ion channel, which is in good agreement with a recent experimental study. In the absence of an experimental structure for the protein, the proposed structure can serve as a useful model for inferring structure-function relationships about the protein.

Graphical Abstract

Side-chain to main-chain hydrogen bonding controls the intrinsic backbone dynamics of the amyloid precursor protein transmembrane helix

Many transmembrane helices contain serine and/or threonine residues whose side chains form intrahelical H-bonds with upstream carbonyl oxygens. Here, we investigated the impact of threonine side-chain/main-chain backbonding on the backbone dynamics of the amyloid precursor protein transmembrane helix. This helix consists of a N-terminal dimerization region and a C-terminal cleavage region, which is processed by γ-secretase to a series of products. Threonine mutations within this transmembrane helix are known to alter the cleavage pattern, which can lead to early-onset Alzheimer's disease. Circular dichroism spectroscopy and amide exchange experiments of synthetic transmembrane domain peptides reveal that mutating threonine enhances the flexibility of this helix. Molecular dynamics simulations show that the mutations reduce intrahelical amide H-bonding and H-bond lifetimes. In addition, the removal of side-chain/main-chain backbonding distorts the helix, which alters bending and rotation at a diglycine hinge connecting the dimerization and cleavage regions. We propose that the backbone dynamics of the substrate profoundly affects the way by which the substrate is presented to the catalytic site within the enzyme. Changing this conformational flexibility may thus change the pattern of proteolytic processing.

Model for coupled insertion and folding of membrane-spanning proteins

Current understanding of the forces directing the folding of integral membrane proteins is very limited compared to the detailed picture available for water-soluble proteins. While mechanistic studies of the folding process in vitro have been conducted for only a small number of membrane proteins, the available evidence indicates that their folding process is thermodynamically driven like that of soluble proteins. In vivo, however, the majority of integral membrane proteins are installed in membranes by dedicated machinery, suggesting that the cellular systems may act to facilitate and regulate the spontaneous physical process of folding. Both the in vitro folding process and the in vivo pathway must navigate an energy landscape dominated by the energetically favorable burial of hydrophobic segments in the membrane interior and the opposition to folding due to the need for passage of polar segments across the membrane. This manuscript describes a simple, exactly solvable model which incorporates these essential features of membrane protein folding. The model is used to compare the folding time under conditions which depict both the in vitro and in vivo pathways. It is proposed that the cellular complexes responsible for insertion of membrane proteins act by lowering the energy barrier for passage of polar regions through the membrane, thereby allowing the chain to more rapidly achieve the folded state.

Implicit membrane treatment of buried charged groups: application to peptide translocation across lipid bilayers

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

Partitioning of amino acids into a model membrane: capturing the interface

Energetics of protein side chain partitioning between aqueous solution and cellular membranes is of fundamental importance for correctly capturing the membrane binding and specific protein–lipid interactions of peripheral membrane proteins. We recently reported a highly mobile membrane mimetic (HMMM) model that accelerates lipid dynamics by modeling the membrane interior partially as a fluid organic solvent while retaining a literal description of the lipid head groups and the beginning of the tails. While the HMMM has been successfully applied to study spontaneous insertion of a number of peripheral proteins into membranes, a quantitative characterization of the energetics of membrane–protein interactions in HMMM membranes has not been performed. We report here the free energy profiles for partitioning of 10 protein side chain analogues into a HMMM membrane. In the interfacial and headgroup regions of the membrane, the side chain free energy profiles show excellent agreement with profiles previously reported for conventional membranes with full-tail lipids. In regions where the organic solvent is prevalent, the increased dipole and fluidity of the solvent generally result in a less accurate description, most notably overstabilization of aromatic and polar amino acids. As an additional measure of the ability of the HMMM model to describe membrane–protein interactions, the water-to-membrane interface transfer energies were analyzed and found to be in agreement with the previously reported experimental and computational hydrophobicity scales. We discuss strengths and weaknesses of HMMM in describing protein–membrane interactions as well as further development of model membranes.

Characterization of the water defect at the HIV-1 gp41 membrane spanning domain in bilayers with and without cholesterol using molecular simulations

The membrane spanning domain (MSD) of human immunodeficiency virus 1 (HIV-1) envelope glycoprotein gp41 is important for fusion and infection. We used molecular dynamics (MD) simulations (3.4 μs total) to relate membrane and peptide properties that lead to water solvation of the α-helical gp41 MSD's midspan arginine in pure dipalmitoylphosphatidylcholine (DPPC) and in 50/50 DPPC/cholesterol membranes. We find that the midspan arginine is solvated by water that penetrates the inner leaflet, leading to a so-called water defect. The water defect is surprisingly robust across initial conditions and membrane compositions, but the presence of cholesterol modulates its behavior in several key ways. In the cholesterol-containing membranes, fluctuations in membrane thickness and water penetration depth are localized near the midspan arginine, and the MSD helices display a tightly regulated tilt angle. In the cholesterol-free membranes, thickness fluctuations are not as strongly correlated to the peptide position and tilt angles vary significantly depending on protein position relative to boundaries between domains of differing thickness. Cholesterol in an HIV-1 viral membrane is required for infection. Therefore, this work suggests that the colocalized water defect and membrane thickness fluctuations in cholesterol-containing viral membranes play an important role in fusion by bringing the membrane closer to a stability limit that must be crossed for fusion to occur.

K(+) channels: function-structural overview

Potassium channels are particularly important in determining the shape and duration of the action potential, controlling the membrane potential, modulating hormone secretion, epithelial function and, in the case of those K(+) channels activated by Ca(2+), damping excitatory signals. The multiplicity of roles played by K(+) channels is only possible to their mammoth diversity that includes at present 70 K(+) channels encoding genes in mammals. Today, thanks to the use of cloning, mutagenesis, and the more recent structural studies using x-ray crystallography, we are in a unique position to understand the origins of the enormous diversity of this superfamily of ion channels, the roles they play in different cell types, and the relations that exist between structure and function. With the exception of two-pore K(+) channels that are dimers, voltage-dependent K(+) channels are tetrameric assemblies and share an extremely well conserved pore region, in which the ion-selectivity filter resides. In the present overview, we discuss in the function, localization, and the relations between function and structure of the five different subfamilies of K(+) channels: (a) inward rectifiers, Kir; (b) four transmembrane segments-2 pores, K2P; (c) voltage-gated, Kv; (d) the Slo family; and (e) Ca(2+)-activated SK family, SKCa.

Effects of Phospholipid Composition on the Transfer of a Small Cationic Peptide Across a Model Biological Membrane

The transfer of a lysine amino acid analogue across phospholipid membrane models was investigated using molecular-dynamics simulations. The evolution of the protonation state of this small peptide as a function of its position inside the membrane was studied by determining the local pKa by means of free-energy calculations. Permeability and mean-first-passage time were evaluated and showed that the transfer occurs on the submillisecond time scale. Comparative studies were conducted to evaluate changes in the pKa arising from differences in the phospholipid chemical structure. We compared, hence, the effect of an ether vs an ester linkage of the lipid headgroup as well as linear vs branched lipid tails. The study reveals that protonated lysine residues can be buried further inside an ether lipid membrane than an ester lipid membrane, while branched lipids are found to stabilize less the charged form compared to their unbranched lipid chain counterparts.

Ser/Thr motifs in transmembrane proteins: conservation patterns and effects on local protein structure and dynamics

We combined systematic bioinformatics analyses and molecular dynamics simulations to assess the conservation patterns of Ser and Thr motifs in membrane proteins, and the effect of such motifs on the structure and dynamics of α-helical transmembrane (TM) segments. We find that Ser/Thr motifs are often present in β-barrel TM proteins. At least one Ser/Thr motif is present in almost half of the sequences of α-helical proteins analyzed here. The extensive bioinformatics analyses and inspection of protein structures led to the identification of molecular transporters with noticeable numbers of Ser/Thr motifs within the TM region. Given the energetic penalty for burying multiple Ser/Thr groups in the membrane hydrophobic core, the observation of transporters with multiple membrane-embedded Ser/Thr is intriguing and raises the question of how the presence of multiple Ser/Thr affects protein local structure and dynamics. Molecular dynamics simulations of four different Ser-containing model TM peptides indicate that backbone hydrogen bonding of membrane-buried Ser/Thr hydroxyl groups can significantly change the local structure and dynamics of the helix. Ser groups located close to the membrane interface can hydrogen bond to solvent water instead of protein backbone, leading to an enhanced local solvation of the peptide.

Cellular uptake of nanoparticles by membrane penetration: a study combining confocal microscopy with FTIR spectroelectrochemistry

It is well-known that nanomaterials are capable of entering living cells, often by utilizing the cells' endocytic mechanisms. Passive penetration of the lipid bilayer may, however, occur as an alternative process. Here we have focused on the passive transport of small nanoparticles across the plasma membranes of red blood cells, which are incapable of endocytosis. By using fluorescence microscopy, we have observed that zwitterionic quantum dots penetrate through the cell membranes so that they can be found inside the cells. The penetration-induced structural changes of the lipid bilayer were explored by surface-enhanced infrared absorption spectroscopy and electrochemistry studies of model membranes prepared on solid supports with lipid compositions identical to those of red blood cell membranes. A detailed analysis of the infrared spectra revealed a markedly enhanced flexibility of the lipid bilayers in the presence of nanoparticles. The electrochemistry data showed that the overall membrane structure remained intact; however, no persistent holes were formed in the bilayers.

Membrane protein misassembly in disease

Helix-helix interactions play a central role in the folding and assembly of integral α-helical membrane proteins and are fundamentally dictated by the amino acid sequence of the TM domain. It is not surprising then that missense mutations that target these residues are often linked to disease. In this review, we focus on the molecular mechanisms through which missense mutations lead to aberrant folding and/or assembly of these proteins, and then discuss pharmacological approaches that may potentially mitigate or reverse the negative effects of these mutations. Improving our understanding of how missense mutations affect the interactions between TM α-helices will increase our capability to develop effective therapeutic approaches to counter the misassembly of these proteins and, ultimately, disease. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Outer membrane phospholipase A in phospholipid bilayers: a model system for concerted computational and experimental investigations of amino acid side chain partitioning into lipid bilayers

Understanding the forces that stabilize membrane proteins in their native states is one of the contemporary challenges of biophysics. To date, estimates of side chain partitioning free energies from water to the lipid environment show disparate values between experimental and computational measures. Resolving the disparities is particularly important for understanding the energetic contributions of polar and charged side chains to membrane protein function because of the roles these residue types play in many cellular functions. In general, computational free energy estimates of charged side chain partitioning into bilayers are much larger than experimental measurements. However, the lack of a protein-based experimental system that uses bilayers against which to vet these computational predictions has traditionally been a significant drawback. Moon & Fleming recently published a novel hydrophobicity scale that was derived experimentally by using a host-guest strategy to measure the side chain energetic perturbation due to mutation in the context of a native membrane protein inserted into a phospholipid bilayer. These values are still approximately an order of magnitude smaller than computational estimates derived from molecular dynamics calculations from several independent groups. Here we address this discrepancy by showing that the free energy differences between experiment and computation become much smaller if the appropriate comparisons are drawn, which suggests that the two fields may in fact be converging. In addition, we present an initial computational characterization of the Moon & Fleming experimental system used for the hydrophobicity scale: OmpLA in DLPC bilayers. The hydrophobicity scale used OmpLA position 210 as the guest site, and our preliminary results demonstrate that this position is buried in the center of the DLPC membrane, validating its usage in the experimental studies. We further showed that the introduction of charged Arg at position 210 is well tolerated in OmpLA and that the DLPC bilayers accommodate this perturbation by creating a water dimple that allows the Arg side chain to remain hydrated. Lipid head groups visit the dimple and can hydrogen bond with Arg, but these interactions are transient. Overall, our study demonstrates the unique advantages of this molecular system because it can be interrogated by both computational and experimental practitioners, and it sets the stage for free energy calculations in a system for which there is unambiguous experimental data. This article is part of a Special Issue entitled: Membrane protein structure and function.

Positions of polar amino acids alter interactions between transmembrane segments and detergents

α-Helical transmembrane (TM) segments in membrane proteins are comprised primarily of hydrophobic amino acids that accommodate insertion from water into the nonpolar membrane bilayer. In many such segments, however, polar residues are also present for structural or functional reasons. These latter residues impair the local favorable acyl interactions required for solvation by hydrophobic media such as phospholipids in native bilayers or detergents used for in vitro characterization. Using a series of Lys-tagged designed TM-like peptides (typified by KK-YAAAIAAIAWAIAAIAAAIAA-KKK) in which single-Asn residue substitutions (from Ile or Ala) were made successively from the center of the hydrophobic region toward the C-terminus, we demonstrate that polar residues strongly alter the nature of the interaction between TM segments and the solvating detergent. Through the application of sodium dodecyl sulfate-polyacrylamide gel electrophoresis, circular dichroism spectroscopy, and tryptophan fluorescence, we observed drastic differences in the structures of the detergent-peptide complexes that contain relatively minor sequence differences. For example, the blue shift of the Trp fluorescence (indicating local detergent solvation at this location) differs by as much as ~10 nm depending upon the position of a single Asn substitution in an otherwise identical segment. The overall results suggest that polar point mutations occurring in a biological membrane will elicit comparable effects, placing a significant refolding burden on the local protein structure and potentially leading to disease states through altered protein--lipid interactions in membrane proteins.

Residue-specific side-chain packing determines the backbone dynamics of transmembrane model helices

The transmembrane domains (TMDs) of membrane-fusogenic proteins contain an overabundance of β-branched residues. In a previous effort to systematically study the relation among valine content, fusogenicity, and helix dynamics, we developed model TMDs that we termed LV-peptides. The content and position of valine in LV-peptides determine their fusogenicity and backbone dynamics, as shown experimentally. Here, we analyze their conformational dynamics and the underlying molecular forces using molecular-dynamics simulations. Our study reveals that backbone dynamics is correlated with the efficiency of side-chain to side-chain van der Waals packing between consecutive turns of the helix. Leu side chains rapidly interconvert between two rotameric states, thus favoring contacts to its i±3 and i±4 neighbors. Stereochemical restraints acting on valine side chains in the α-helix force both β-substituents into an orientation where i,i±3 interactions are less favorable than i,i±4 interactions, thus inducing a local packing deficiency at VV3 motifs. We provide a quantitative molecular model to explain the relationship among chain connectivity, side-chain mobility, and backbone flexibility. We expect that this mechanism also defines the backbone flexibility of natural TMDs.

Helical integrity and microsolvation of transmembrane domains from Flaviviridae envelope glycoproteins

Charged and polar amino acids in the transmembrane domains of integral membrane proteins can be crucial for protein function and also promote helix-helix association or protein oligomerization. Yet, our current understanding is still limited on how these hydrophilic amino acids are efficiently translocated from the Sec61/SecY translocon into the cell membrane during the biogenesis of membrane proteins. In hepatitis C virus, the putative transmembrane segments of envelope glycoproteins E1 and E2 were suggested to heterodimerize via a Lys-Asp ion-pair in the host endoplasmic reticulum. Therefore in this work, we carried out molecular dynamic simulations in explicit lipid bilayer and solvent environment to explore the stability of all possible bridging ion-pairs using the model of H-segment helix dimers. We observed that, frequently, several water molecules penetrated from the interface into the membrane core to stabilize the charged and polar pairs. The hydration time and amount of water molecules in the membrane core depended on the position of the charged residues as well as on the type of ion-pairs. Similar microsolvation events were observed in simulations of the putative E1-E2 transmembrane helix dimers. Simulations of helix monomers from other members of the Flaviviridae family suggest that these systems show similar behaviors. Thus this study illustrates the important contribution of water microsolvation to overcome the unfavorable energetic cost of burying charged and polar amino acids in membrane lipid bilayers. Also, it emphasizes the novel role of bridging charged or polar interactions stabilized by water molecules in the hydrophobic lipid bilayer core that has an important biological function for helix dimerization in several envelope glycoproteins from the family of Flaviviridae viruses.

Determining peptide partitioning properties via computer simulation

The transfer of polypeptide segments into lipid bilayers to form transmembrane helices represents the crucial first step in cellular membrane protein folding and assembly. This process is driven by complex and poorly understood atomic interactions of peptides with the lipid bilayer environment. The lack of suitable experimental techniques that can resolve these processes both at atomic resolution and nanosecond timescales has spurred the development of computational techniques. In this review, we summarize the significant progress achieved in the last few years in elucidating the partitioning of peptides into lipid bilayer membranes using atomic detail molecular dynamics simulations. Indeed, partitioning simulations can now provide a wealth of structural and dynamic information. Furthermore, we show that peptide-induced bilayer distortions, insertion pathways, transfer free energies, and kinetic insertion barriers are now accurate enough to complement experiments. Further advances in simulation methods and force field parameter accuracy promise to turn molecular dynamics simulations into a powerful tool for investigating a wide range of membrane active peptide phenomena.

Why are polar residues within the membrane core evolutionary conserved?

Here, we present a study of polar residues within the membrane core of alpha-helical membrane proteins. As expected, polar residues are less frequent in the membrane than expected. Further, most of these residues are buried within the interior of the protein and are only rarely exposed to lipids. However, the polar groups often border internal water filled cavities, even if the rest of the sidechain is buried. A survey of their functional roles in known structures showed that the polar residues are often directly involved in binding of small compounds, especially in channels and transporters, but other functions including proton transfer, catalysis, and selectivity have also been attributed to these proteins. Among the polar residues histidines often interact with prosthetic groups in photosynthetic- and oxidoreductase-related proteins, whereas prolines often are required for conformational changes of the proteins. Indeed, the polar residues in the membrane core are more conserved than other residues in the core, as well as more conserved than polar residues outside the membrane. The reason is twofold; they are often (i) buried in the interior of the protein and (ii) directly involved in the function of the proteins. Finally, a method to identify which polar residues are present within the membrane core directly from protein sequences was developed. Applying the method to the set of all human membrane proteins the prediction indicates that polar residues were most frequent among active transporter proteins and GPCRs, whereas infrequent in families with few transmembrane regions, such as non-GPCR receptors.

Electrostatics of deformable lipid membranes

It was recently demonstrated that significant local deformations of biological membranes take place due to the fields of charged peptides and ions, challenging the standard model of membrane electrostatics. The ability of ions to retain their immediate hydration environment, combined with the lack of sensitivity of permeability to ion type or even ion pairs, led us to question the extent to which hydration energetics and electrostatics control membrane ion permeation. Using the arginine analog methyl-guanidinium as a test case, we find that although hydrocarbon electronic polarizability causes dramatic changes in ion solvation free energy, as well as a significant change (approximately 0.4 V) in the membrane dipole potential, little change in membrane permeation energetics occurs. We attribute this to compensation of solvation terms from polar and polarizable nonpolar components within the membrane, and explain why the dipole potential is not fully sensed in terms of the locally deformed bilayer interface. Our descriptions provide a deeper understanding of the translocation process and allow predictions for poly-ions, ion pairs, charged lipids, and lipid flip-flop. We also report simulations of large hydrophobic-ion-like membrane defects and the ionophore valinomycin, which exhibit little membrane deformation, as well as hydrophilic defects and the ion channel gramicidin A, to provide parallels to membranes deformed by unassisted ion permeation.

Contribution of charged and polar residues for the formation of the E1-E2 heterodimer from Hepatitis C Virus

The transmembrane domains of the envelope glycoprotein E1 and E2 have crucial multifunctional roles in the biogenesis of hepatitis C virus. We have performed molecular dynamics simulations to investigate a structural model of the transmembrane segments of the E1-E2 heterodimer. The simulations support the key role of the Lys370-Asp728 ion pair for mediating the E1-E2 heterodimerization. In comparison to these two residues, the simulation results also reveal the differential effect of the conserved Arg730 residue that has been observed in experimental studies. Furthermore, we discovered the formation of inter-helical hydrogen bonds via Asn367 that stabilize dimer formation. Simulations of single and double mutants further demonstrate the importance of the ion-pair and polar interactions between the interacting helix monomers. The conformation of the E1 fragment in the simulation of the E1-E2 heterodimer is in close agreement with an NMR structure of the E1 transmembrane segment. The proposed model of the E1-E2 heterodimer supports the postulated cooperative insertion of both helices by the translocon complex into the bilayer.

Reduced voltage sensitivity in a K+-channel voltage sensor by electric field remodeling

Propagation of the nerve impulse relies on the extreme voltage sensitivity of Na(+) and K(+) channels. The transmembrane movement of four arginine residues, located at the fourth transmembrane segment (S4), in each of their four voltage-sensing domains is mostly responsible for the translocation of 12 to 13 e(o) across the transmembrane electric field. Inserting additional positively charged residues between the voltage-sensing arginines in S4 would, in principle, increase voltage sensitivity. Here we show that either positively or negatively charged residues added between the two most external sensing arginines of S4 decreased voltage sensitivity of a Shaker voltage-gated K(+)-channel by up to approximately 50%. The replacement of Val363 with a charged residue displaced inwardly the external boundaries of the electric field by at least 6 A, leaving the most external arginine of S4 constitutively exposed to the extracellular space and permanently excluded from the electric field. Both the physical trajectory of S4 and its electromechanical coupling to open the pore gate seemed unchanged. We propose that the separation between the first two sensing charges at resting is comparable to the thickness of the low dielectric transmembrane barrier they must cross. Thus, at most a single sensing arginine side chain could be found within the field. The conserved hydrophobic nature of the residues located between the voltage-sensing arginines in S4 may shape the electric field geometry for optimal voltage sensitivity in voltage-gated ion channels.

Characterization of membrane-protein interactions for the leucine transporter from Aquifex aeolicus by molecular dynamics calculations

Multinanosecond molecular dynamics (MD) simulations have been employed to characterize the interaction of an integral membrane protein (IMP), the leucine transmitter from Aquifex aeolicus (Yamashita et al., Nature 2005, 437, 215-223), with hydrated lipid bilayer membranes in their physiologically relevant liquid crystalline phases. Analysis of the MD trajectories for dimyristoyl phosphatidylcholine (DMPC), 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC), and 1-palmitoyl-2-oleoyl phosphatidylethanolamine (POPE) focused on the contacts between aromatic and basic side chains of the IMP with the lipid head groups and water. Structural fluctuations of the IMP were investigated as well as the contact dynamics of neighboring lipids. In characterizing the IMP-membrane systems, the behaviors of the protein's cytoplasmic and periplasmic parts are considered separately. All three lipid membranes show a rather similar overall level of association with the IMP. However, for DMPC there is a better matching of the membrane core to the hydrophobic transmembrane portion of the IMP. The closed cytoplasmic end of the IMP exhibits a higher degree of association with lipids than the more open periplasmic end, an observation which correlates with the more compact structure and a slower dynamics of surrounding lipids.

Membrane insertion of marginally hydrophobic transmembrane helices depends on sequence context

In mammalian cells, most integral membrane proteins are initially inserted into the endoplasmic reticulum membrane by the so-called Sec61 translocon. However, recent predictions suggest that many transmembrane helices (TMHs) in multispanning membrane proteins are not sufficiently hydrophobic to be recognized as such by the translocon. In this study, we have screened 16 marginally hydrophobic TMHs from membrane proteins of known three-dimensional structure. Indeed, most of these TMHs do not insert efficiently into the endoplasmic reticulum membrane by themselves. To test if loops or TMHs immediately upstream or downstream of a marginally hydrophobic helix might influence the insertion efficiency, insertion of marginally hydrophobic helices was also studied in the presence of their neighboring loops and helices. The results show that flanking loops and nearest-neighbor TMHs are sufficient to ensure the insertion of many marginally hydrophobic helices. However, for at least two of the marginally hydrophobic helices, the local interactions are not enough, indicating that post-insertional rearrangements are involved in the folding of these proteins.

Protein contents in biological membranes can explain abnormal solvation of charged and polar residues

Transmembrane helices are generally believed to insert into membranes based on their hydrophobicity. Nevertheless, there are important exceptions where polar residues have great functional importance, for instance the S4 helix of voltage-gated ion channels. It has been shown experimentally that insertion can be accomplished by hydrophobic counterbalance, predicting an arginine insertion cost of only 2.5 kcal/mol, compared with 14.9 kcal/mol in cyclohexane. Previous simulations of pure bilayers have produced values close to the pure hydrocarbon, which has lead to spirited discussion about the experimental conditions. Here, we have performed computer simulations of models better mimicking biological membranes by explicitly including protein helices at mass fractions from 15% to 55%, as well as an actual translocon. This has a striking effect on the solvation free energy of arginine. With some polar residues present, the solvation cost comes close to experimental observation at approximately 30% mass fraction, and negligible at 40%. In the presence of a translocon in the membrane, the cost of inserting arginine next to the lateral gate can be as low as 3-5 kcal/mol. The effect is mainly due to the extra helices making it easier to retain hydration water. These results offer a possible explanation for the discrepancy between the in vivo hydrophobicity scale and computer simulations and highlight the importance of the high protein contents in membranes. Although many membrane proteins are stable in pure bilayers, such simplified models might not be sufficiently accurate for insertion of polar or charged residues in biological membranes.

Tilt and azimuthal angles of a transmembrane peptide: a comparison between molecular dynamics calculations and solid-state NMR data of sarcolipin in lipid membranes

We report molecular dynamics simulations in the explicit membrane environment of a small membrane-embedded protein, sarcolipin, which regulates the sarcoplasmic reticulum Ca-ATPase activity in both cardiac and skeletal muscle. In its monomeric form, we found that sarcolipin adopts a helical conformation, with a computed average tilt angle of 28 +/- 6 degrees and azymuthal angles of 66 +/- 22 degrees, in reasonable accord with angles determined experimentally (23 +/- 2 degrees and 50 +/- 4 degrees, respectively) using solid-state NMR with separated-local-field experiments. The effects of time and spatial averaging on both (15)N chemical shift anisotropy and (1)H/(15)N dipolar couplings have been analyzed using short-time averages of fast amide out-of-plane motions and following principal component dynamic trajectories. We found that it is possible to reproduce the regular oscillatory patterns observed for the anisotropic NMR parameters (i.e., PISA wheels) employing average amide vectors. This work highlights the role of molecular dynamics simulations as a tool for the analysis and interpretation of solid-state NMR data.