Conformational dynamics of helix S6 from Shaker potassium channel: simulation studies

Prolines in transmembrane (TM) alpha-helices are believed to play an important structural and/or functional role in membrane proteins. At a structural level a proline residue distorts alpha-helical structure due to the loss of at least one stabilizing backbone hydrogen bond, and introduces flexibility in the helix that may result in substantial kink and swivel motions about the effective "hinge." At a functional level, for example in Kv channels, it is believed that proline-induced molecular hinges may have a direct role in gating, i.e., the conformational change linked to opening/closing the channel to movement of ions. In this article we study the conformational dynamics of the S6 TM helix from of the Kv channel Shaker, which possesses the motif PVP--a motif that is conserved in Kv channels. We perform multiple molecular dynamics simulations of single S6 helices in a membrane-mimetic environment in order to effectively map the kink-swivel conformational space of the protein, exploiting the ability of multiple simulations to achieve greater sampling. We show that the presence of proline locally perturbs the helix, disrupting local dihedral angles and producing local twist and unwinding in the region of the hinge--an effect that is relaxed with distance from the PVP motif. We furthermore show that motions about the hinge are highly anisotropic, reflecting a preferred region of kink-swivel conformation space that may have implications for the gating process.     

Visualization of conformational distribution of short to medium size segments in globular proteins and identification of local structural motifs

Analysis of the conformational distribution of polypeptide segments in a conformational space is the first step for understanding a principle of structural diversity of proteins. Here, we present a statistical analysis of protein local structures based on interatomic C(alpha) distances. Using principal component analysis (PCA) on the intrasegment C(alpha)-C(alpha) atomic distances, the conformational space of protein segments, which we call the protein segment universe, has been visualized, and three essential coordinate axes, suitable for describing the universe, have been identified. Three essential axes specified radius of gyration, structural symmetry, and separation of hairpin structures from other structures. Among the segments of arbitrary length, 6-22 residues long, the conservation of those axes was uncovered. Further application of PCA to the two largest clusters in the universe revealed local structural motifs. Although some of motifs have already been reported, we identified a possibly novel strand motif. We also showed that a capping box, which is one of the helix capping motifs, was separated into independent subclusters based on the C(alpha) geometry. Implications of the strand motif, which may play a role for protein-protein interaction, are discussed. The currently proposed method is useful for not only mapping the immense universe of protein structures but also identification of structural motifs.     

Sequence determinants of the capping box, a stabilizing motif at the N-termini of alpha-helices.

The capping box, a recurrent hydrogen bonded motif at the N-termini of alpha-helices, caps 2 of the initial 4 backbone amide hydrogen donors of the helix (Harper ET, Rose GD, 1993, Biochemistry 32:7605-7609). In detail, the side chain of the first helical residue forms a hydrogen bond with the backbone of the fourth helical residue and, reciprocally, the side chain of the fourth residue forms a hydrogen bond with the backbone of the first residue. We now enlarge the earlier definition of this motif to include an accompanying hydrophobic interaction between residues that bracket the capping box sequence on either side. The expanded box motif--in which 2 hydrogen bonds and a hydrophobic interaction are localized within 6 consecutive residues--resembles a glycine-based capping motif found at helix C-termini (Aurora R, Srinivasan R, Rose GD, 1994, Science 264:1126-1130).     

Comparison of the human and worm p53 structures suggests a way for enhancing stability

Maintaining the native conformation is essential for the proper function of tumor suppressor protein p53. However, p53 is a low-stability protein that can easily lose its function upon structural perturbations such as those resulting from missense mutations, leading to the development of cancer. Therefore, it is important to develop strategies to design stable p53 which still maintains its normal function. Here, we compare the stabilities of the human and worm p53 core domains using molecular dynamics simulations. We find that the worm p53 is significantly more stable than the human form. Detailed analysis of the structural fluctuations shows that the stability difference lies in the peripheral structural motifs that contrast in their structural features and flexibility. The most dramatic difference in stability originates from loop L1, from the turn between helix H1 and beta-strand S5, and from the turn that connects beta-strands S7 and S8. Structural analysis shows significant differences for these motifs between the two proteins. Loop L1 lacks secondary structure, and the turns between helix H1 and strand S5 and between strands S7 and S8 are much longer in the human form p53. On the basis of these differences, we designed a mutant by shortening the turn between strands S7 and S8 to enhance the stability. Surprisingly, this mutant was very stable when probed by molecular dynamics simulations. In addition, the stabilization was not localized in the turn region. Loop L1 was also significantly stabilized. Our results show that stabilizing peripheral structural motifs can greatly enhance the stability of the p53 core domain and therefore is likely to be a viable alternative in the development of stable p53. In addition, loop- or turn-related mutants with different stabilities may also be used to probe the relationship between function, a particular structural motif, and its flexibility.     

Molecular dynamics studies of a DNA-binding protein: 1. A comparison of the trp repressor and trp aporepressor aqueous simulations.

The results of two 30-ps molecular dynamics simulations of the trp repressor and trp aporepressor proteins are presented in this paper. The simulations were obtained using the AMBER molecular mechanical force field and in both simulations a 6-A shell of TIP3P waters surrounded the proteins. The trp repressor protein is a DNA-binding regulatory protein and it utilizes a helix-turn-helix (D helix-turn-E helix) motif to interact with DNA. The trp aporepressor, lacking two molecules of the L-tryptophan corepressor, cannot bind specifically to DNA. Our simulations show that the N- and C-termini and the residues in and near the helix-turn-helix motifs are the most mobile regions of the proteins, in agreement with the X-ray crystallographic studies. Our simulations also find increased mobility of the residues in the turn-D helix-turn regions of the proteins. We find the average distance separating the DNA-binding motifs to be larger in the repressor as compared to the aporepressor. In addition to examining the protein residue fluctuations and deviations with respect to X-ray structures, we have also focused on backbone dihedral angles and corepressor hydrogen-bonding patterns in this paper.     

Sequence and conformational preferences at termini of α‐helices in membrane proteins: Role of the helix environment

α‐helices are amongst the most common secondary structural elements seen in membrane proteins and are packed in the form of helix bundles. These α‐helices encounter varying external environments (hydrophobic, hydrophilic) that may influence the sequence preferences at their N and C‐termini. The role of the external environment in stabilization of the helix termini in membrane proteins is still unknown. Here we analyze α‐helices in a high‐resolution dataset of integral α‐helical membrane proteins and establish that their sequence and conformational preferences differ from those in globular proteins. We specifically examine these preferences at the N and C‐termini in helices initiating/terminating inside the membrane core as well as in linkers connecting these transmembrane helices. We find that the sequence preferences and structural motifs at capping (Ncap and Ccap) and near‐helical (N' and C') positions are influenced by a combination of features including the membrane environment and the innate helix initiation and termination property of residues forming structural motifs. We also find that a large number of helix termini which do not form any particular capping motif are stabilized by formation of hydrogen bonds and hydrophobic interactions contributed from the neighboring helices in the membrane protein. We further validate the sequence preferences obtained from our analysis with data from an ultradeep sequencing study that identifies evolutionarily conserved amino acids in the rat neurotensin receptor. The results from our analysis provide insights for the secondary structure prediction, modeling and design of membrane proteins. Proteins 2014; 82:3420–3436. © 2014 Wiley Periodicals, Inc.     

A conformational equilibrium in a protein fragment caused by two consecutive capping boxes: 1H-, 13C-NMR, and mutational analysis.

The conformational properties of an 18 residues peptide spanning the entire sequence, L1KTPA5QFDAD10ELRAA15MKG, of the first helix (A-helix) of domain 2 of annexin I, were thoroughly investigated. This fragment exhibits several singular features, and in particular, two successive potential capping boxes, T3xxQ6 and D8xxE11. The former corresponds to the native hydrogen bond network stabilizing the alpha helix N-terminus in the protein; the latter is a non-native capping box able to break the helix at residue D8, and is observed in the domain 2 partially folded state. Using 2D-NMR techniques, we showed that two main populations of conformers coexist in aqueous solution. The first corresponds to a single helix extending from T3 to K17. The second corresponds to a broken helix at residue Ds. Four mutants, T3A, F7A, D8A, and E11A, were designed to further analyze the role of key amino acids in the equilibrium between the two ensembles of conformers. The sensitivity of NMR parameters to account for the variations in the populations of conformers was evaluated for each peptide. Our data show the delta13Calpha chemical shift to be the most relevant parameter. We used it to estimate the population ratio in the various peptides between the two main ensembles of conformers, the full helix and the broken helix. For the WT, E11A, and F7A peptides, these ratios are respectively 35/65, 60/40, 60/40. Our results were compared to the data obtained from helix/coil transition algorithms.     

Pathways and intermediates in forced unfolding of spectrin repeats

Spectrin repeats are triple-helical coiled-coil domains found in many proteins that are regularly subjected to mechanical stress. We used atomic force microscopy technique and steered molecular dynamics simulations to study the behavior of a wild-type spectrin repeat and two mutants. The experiments indicate that spectrin repeats can form stable unfolding intermediates when subjected to external forces. In the simulations the unfolding proceeded via a variety of pathways. Stable intermediates were associated to kinking of the central helix close to a proline residue. A mutant stabilizing the central helix showed no intermediates in experiments, in agreement with simulation. Spectrin repeats may thus function as elastic elements, extendable to intermediate states at various lengths.     

Motifs of two small residues can assist but are not sufficient to mediate transmembrane helix interactions

Both experimental and statistical searches for specific motifs that mediate transmembrane helix-helix interactions showed that two glycine residues separated by three intervening residues (GxxxG) provide a framework for specific interactions. Further work suggested that other motifs of small residues can mediate the interaction of transmembrane domains, so that the AxxxA-motif could also drive strong interactions of alpha-helices in soluble proteins. Thus, all these data indicate that a motif of two small residues in a distance of four might be enough to provide a framework for transmembrane helix-helix interaction. To test whether GxxxG is equivalent to (small)xxx(small), we investigated the effect of a substitution of either of the two Gly residues in the glycophorin A GxxxG-motif by Ala or Ser using the recently developed GALLEX system. The results of this mutational study demonstrate that, while a replacement of either of the two Gly by Ala strongly disrupts GpA homo-dimerization, the mutation to Ser partly stabilizes a dimeric structure. We suggest that the Ser residue can form a hydrogen bond with a backbone carbonyl group of the adjacent helix stabilizing a preformed homo-dimer. While (small)xxx(small) serves as a useful clue, the context of adjacent side-chains is essential for stable helix interaction, so each case must be tested.     

Molecular dynamics simulations of isolated transmembrane helices of potassium channels

In the middle of the S6 helix in voltage-gated potassium channels there is a highly conserved Pro-Val-Pro motif, while the equivalent M2 helix of inward rectifier potassium channels contains a conserved glycine residue in a comparable position. The structural implications of these conserved motifs are of interest given the evidence that S6 and M2 are components of the lining of their respective pores. Multiple sequence alignment and TM helix prediction methods were used to define consensus regions for S6 and M2. Ensembles of 50 structures for each helix were generated by simulated annealing and restrained molecular dynamics. Time-dependent fluctuations of S6 and M2 were investigated by long time scale molecular dynamics simulations on representative members of each ensemble carried out in vacuo in the presence and absence of a hydrophobic potential that mimics a lipid bilayer. The results are discussed in terms of the structural basis of the kink in S6 and M2 and of a putative functional role for flexible helices as “molecular swivels.” © 1996 John Wiley & Sons, Inc.     

Probing the structural hierarchy and energy landscape of an RNA T-loop hairpin

The T-loop motif is an important recurrent RNA structural building block consisting of a U-turn sub-motif and a UA trans Watson–Crick/Hoogsteen base pair. In the presence of a hairpin stem, the UA non-canonical base pair becomes part of the UA-handle motif. To probe the hierarchical organization and energy landscape of the T-loop, we performed replica exchange molecular dynamics (REMD) simulations of the T-loop in isolation and as part of a hairpin. Our simulations reveal that the isolated T-loop adopts coil conformers stabilized by base stacking. The T-loop hairpin shows a highly rugged energy landscape featuring multiple local minima with a transition state for folding consisting of partially zipped states. The U-turn displays a high conformational flexibility both when the T-loop is in isolation and as part of a hairpin. On the other hand, the stability of the UA non-canonical base pair is enhanced in the presence of the UA-handle. This motif is apparently a key component for stabilizing the T-loop, while the U-turn is mostly involved in long-range interaction. Our results suggest that the stability and folding of small RNA motifs are highly dependent on local context.     

Role of the helix capping in the stability of the mouse prion (180-213) segment: investigation through molecular dynamics simulations.

Molecular dynamics simulation of the 180-213 segment, forming the B and C helices in the mouse prion protein, and of three mutants, where the capping box residues or the hydrophobic staple motif residues were selectively mutated, have been carried out. The results indicate that the wild type segment is stable over all the trajectory, whilst the mutants display different degrees of destabilization. In detail mutation of Asp202 brings to a rapid unfolding of helix C likely because of the concomitant loss of a hydrogen bond and of a negative charge able to stabilize the dipole in the first turn of the helix. A lower destabilizing effect is observed upon mutation Thr199. On the other hand mutation of Phe198 and Val203, the hydrophobic staple residues, brings to an incorrect orientation of the first helix relative to the second one due to a weakening of the hydrophobic interaction. The results confirm the importance of the presence of both motifs for the structural integrity of the isolated fragment and suggest that these residues may have a main role in the structural transition observed in the inherited human prion diseases.     

Ribosomal RNA kink-turn motif--a flexible molecular hinge

Ribosomal RNA K-turn motifs are asymmetric internal loops characterized by a sharp bend in the phosphodiester backbone resulting in "V" shaped structures, recurrently observed in ribosomes and showing a high degree of sequence conservation. We have carried out extended explicit solvent molecular dynamics simulations of selected K-turns, in order to investigate their intrinsic structural and dynamical properties. The simulations reveal an unprecedented dynamical flexibility of the K-turns around their X-ray geometries. The K-turns sample, on the nanosecond timescale, different conformational substates. The overall behavior of the simulations suggests that the sampled geometries are essentially isoenergetic and separated by minimal energy barriers. The nanosecond dynamics of isolated K-turns can be qualitatively considered as motion of two rigid helix stems controlled by a very flexible internal loop which then leads to substantial hinge-like motions between the two stems. This internal dynamics of K-turns is strikingly different for example from the bacterial 5S rRNA Loop E motif or BWYV frameshifting pseudoknot which appear to be rigid in the same type of simulations. Bistability and flexibility of K-turns was also suggested by several recent biochemical studies. Although the results of MD simulations should be considered as a qualitative picture of the K-turn dynamics due to force field and sampling limitations, the main advantage of the MD technique is its ability to investigate the region close to K-turn ribosomal-like geometries. This part of the conformational space is not well characterized by the solution experiments due to large-scale conformational changes seen in the experiments. We suggest that K-turns are well suited to act as flexible structural elements of ribosomal RNA. They can for example be involved in mediation of large-scale motions or they can allow a smooth assembling of the other parts of the ribosome.     

Phosphorylation-induced structural changes in the amyloid precursor protein cytoplasmic tail detected by NMR

The cytoplasmic tail of the amyloid precursor protein (APPc) interacts with several cellular factors implicated in intracellular signaling or proteolytic production of amyloid beta peptide found in senile plaques of Alzheimer's disease patients. APPc contains two threonine residues (654 and 668 relative to APP695, or 6 and 20 relative to APPc) and a serine residue (655 or 7, respectively) that are known to be phosphorylated in vivo and may play regulatory roles in these events. We show by solution NMR spectroscopy of a 49 residue cytoplasmic tail peptide (APP-C) that in all three cases, phosphorylation induces changes in backbone dihedral angles that can be attributed to formation of local hydrogen bonds between the phosphate group and nearby amide protons. Phosphorylation of S7 also induces chemical shift changes in the hydrophobic cluster (residues I8-V13), indicating additional medium-range effects. The most pronounced changes occur upon phosphorylation of T20, a neuron-specific phosphorylation site, where the N-terminal helix capping box previously characterized for this region is altered. Characterization of torsion angles and transient hydrogen bonds indicates that prolyl isomerization of the pThr-Pro peptide bond results from both destabilization of the N-terminal helix capping box and stabilization of the cis isomer by transient hydrogen bonds. The significant population of the cis isomer (9 %) present after phosphorylation of T20 suggests a potential role of selective recognition of cis versus trans isomers in response to phosphorylation of APP. Together, these structural changes indicate that phosphorylation may act as a conformational switch in the cytoplasmic tail of APP to alter specificity and affinity of binding to cytosolic partners, particularly in response to the abnormal phosphorylation events associated with Alzheimer's disease.     

Alpha helix capping in synthetic model peptides by reciprocal side chain–main chain interactions: Evidence for an N terminal “capping box”

A significant fraction of the amino acids in proteins are alpha helical in conformation. Alpha helices in globular proteins are short, with an average length of about twelve residues, so that residues at the ends of helices make up an important fraction of all helical residues. In the middle of a helix, H-bonds connect the NH and CO groups of each residue to partners four residues along the chain. At the ends of a helix, the H-bond potential of the main chain remains unfulfilled, and helix capping interactions involving bonds from polar side chains to the NH or CO of the backbone have been proposed and detected. In a study of synthetic helical peptides, we have found that the sequence Ser-Glu-Asp-Glu stabilizes the alpha helix in a series of helical peptides with consensus sequences. Following the report by Harper and Rose, which identifies SerXaaXaaGlu as a member of a class of common motifs at the N termini of alpha helices in proteins that they refer to as “capping boxes,” we have reexamined the side chain–main chain interactions in a varient sequence using ¹H NMR, and find that the postulated reciprocal side chain-backbone bonding between the first Ser and last Glu side chains and their peptide NH partners can be resolved: Deletion of two residues N terminal to the Ser-Glu-Asp-Glu sequence in these peptides has no effect on the initiation of helical structure, as defined by two-dimensional (2D) NMR experiments on this variant. Thus the capping box sequence Ser-Glu-Asp-Glu inhibits N terminal fraying of the N terminus of alpha helix in these peptides, and shows the side chain–main chain interactions proposed by Harper and Rose. It thus acts as a helix initiating signal. Since normal a helix cannot propagate beyond the N terminus of this structure, the box acts as a termination signal in this direction as well. © 1994 John Wiley & Sons, Inc.     

Structures of thermolabile mutants of human glutathione transferase P1-1

An N-capping box motif (Ser/Thr-Xaa-Xaa-Asp) is strictly conserved at the beginning of helix alpha6 in the core of virtually all glutathione transferases (GST) and GST-related proteins. It has been demonstrated that this local motif is important in determining the alpha-helical propensity of the isolated alpha6-peptide and plays a crucial role in the folding and stability of GSTs. Its removal by site-directed mutagenesis generated temperature-sensitive folding mutants unable to refold at physiological temperature (37 degrees C). In the present work, variants of human GSTP1-1 (S150A and D153A), in which the capping residues have been substituted by alanine, have been generated and purified for structural analysis. Thus, for the first time, temperature-sensitive folding mutants of an enzyme, expressed at a permissive temperature, have been crystallized and their three-dimensional structures determined by X-ray crystallography. The crystal structures of human pi class GST temperature-sensitive mutants provide a basis for understanding the structural origin of the dramatic effects observed on the overall stability of the enzyme at higher temperatures upon single substitution of a capping residue.     

An unexpected extended conformation for the third TPR motif of the peroxin PEX5 from Trypanosoma brucei

A number of helix-rich protein motifs are involved in a variety of critical protein-protein interactions in living cells. One of these is the tetratrico peptide repeat (TPR) motif that is involved, amongst others, in cell cycle regulation, chaperone function and post-translation modifications. So far, these helix-rich TPR motifs have always been observed to be a compact unit of two helices interacting with each other in antiparallel fashion. Here, we describe the structure of the first three TPR-motifs of the peroxin PEX5 from Trypanosoma brucei, the causative agent of sleeping sickness. Peroxins are proteins involved in peroxisome, glycosome and glyoxysome biogenesis. PEX5 is the receptor of the proteins targeted to these organelles by the "peroxisomal targeting signal-1", a C-terminal tripeptide called PTS-1. The first two of the three TPR-motifs of T. brucei PEX5 appear to adopt the canonical antiparallel helix hairpin structure. In contrast, the third TPR motif of PEX5 has a dramatically different conformation in our crystals: the two helices that were supposed to form a hairpin are folded into one single 44 A long continuous helix. Such a conformation has never been observed before for a TPR motif. This raises interesting questions including the potential functional importance of a "jack-knife" conformational change in TPR motifs.     

The intrinsic flexibility of the Kv voltage sensor and its implications for channel gating

Analysis of the crystal structures of the intact voltage-sensitive potassium channel KvAP (from Aeropyrum pernix) and Kv1.2 (from rat brain), along with the isolated voltage sensor (VS) domain from KvAP, raises the question of the exact nature of the voltage-sensing conformational change that triggers activation of Kv and related voltage-gated channels. Molecular dynamics simulations of the isolated VS of KvAP in a detergent micelle environment at two different temperatures (300 K and 368 K) have been used to probe the intrinsic flexibility of this domain on a tens-of-nanoseconds timescale. The VS contains a positively charged (S4) helix which is packed against a more hydrophobic S3 helix. The simulations at elevated temperature reveal an intrinsic flexibility/conformational instability of the S3a region (i.e., the C-terminus of the S3 helix). It is also evident that the S4 helix undergoes hinge bending and swiveling about its central I130 residue. The conformational instability of the S3a region facilitates the motion of the N-terminal segment of S4 (i.e., S4a). These simulations thus support a gating model in which, in response to depolarization, an S3b-S4a "paddle" may move relative to the rest of the VS domain. The flexible S3a region may in turn act to help restore the paddle to its initial conformation upon repolarization.     

The Molecular Mechanism of Ion-Dependent Gating in Secondary Transporters

LeuT-like fold Na-dependent secondary active transporters form a large family of integral membrane proteins that transport various substrates against their concentration gradient across lipid membranes, using the free energy stored in the downhill concentration gradient of sodium ions. These transporters play an active role in synaptic transmission, the delivery of key nutrients, and the maintenance of osmotic pressure inside the cell. It is generally believed that binding of an ion and/or a substrate drives the conformational dynamics of the transporter. However, the exact mechanism for converting ion binding into useful work has yet to be established. Using a multi-dimensional path sampling (string-method) followed by all-atom free energy simulations, we established the principal thermodynamic and kinetic components governing the ion-dependent conformational dynamics of a LeuT-like fold transporter, the sodium/benzyl-hydantoin symporter Mhp1, for an entire conformational cycle. We found that inward-facing and outward-facing states of Mhp1 display nearly the same free energies with an ion absent from the Na2 site conserved across the LeuT-like fold transporters. The barrier separating an apo-state from inward-facing or outward-facing states of the transporter is very low, suggesting stochastic gating in the absence of ion/substrate bound. In contrast, the binding of a Na2 ion shifts the free energy stabilizing the outward-facing state and promoting substrate binding. Our results indicate that ion binding to the Na2 site may also play a key role in the intracellular thin gate dynamics modulation by altering its interactions with the transmembrane helix 5 (TM5). The Potential of Mean Force (PMF) computations for a substrate entrance displays two energy minima that correspond to the locations of the main binding site S1 and proposed allosteric S2 binding site. However, it was found that substrate''s binds to the site S1 ∼5 kcal/mol more favorable than that to the site S2 for all studied bound combinations of ions and a substrate.     

Proline-induced hinges in transmembrane helices: possible roles in ion channel gating

A number of ion channels contain transmembrane (TM) alpha-helices that contain proline-induced molecular hinges. These TM helices include the channel-forming peptide alamethicin (Alm), the S6 helix from voltage-gated potassium (Kv) channels, and the D5 helix from voltage-gated chloride (CLC) channels. For both Alm and KvS6, experimental data implicate hinge-bending motions of the helix in an aspect of channel gating. We have compared the hinge-bending motions of these TM helices in bilayer-like environments by multi-nanosecond MD simulations in an attempt to describe motions of these helices that may underlie possible modes of channel gating. Alm is an alpha-helical channel-forming peptide, which contains a central kink associated with a Gly-x-x-Pro motif in its sequence. Simulations of Alm in a TM orientation for 10 ns in an octane slab indicate that the Gly-x-x-Pro motif acts as a molecular hinge. The S6 helix from Shaker Kv channels contains a Pro-Val-Pro motif. Modeling studies and recent experimental data suggest that the KvS6 helix may be kinked in the vicinity of this motif. Simulations (10 ns) of an isolated KvS6 helix in an octane slab and in a POPC bilayer reveal hinge-bending motions. A pattern-matching approach was used to search for possible hinge-bending motifs in the TM helices of other ion channel proteins. This uncovered a conserved Gly-x-Pro motif in TM helix D5 of CLC channels. MD simulations of a model of hCLC1-D5 spanning an octane slab suggest that this channel also contains a TM helix that undergoes hinge-bending motion. In conclusion, our simulations suggest a model in which hinge-bending motions of TM helices may play a functional role in the gating mechanisms of several different families of ion channels.