Model for the conformational analysis of hydrated peptides. Effect of hydration on the conformational stability of the terminally blocked residues of the 20 naturally occurring amino acids

The effects of hydration are included in empirical conformational energy computations on oligopeptides by means of a modified hydration-shell model. Free energy terms are introduced to account for “specific hydration” due to water–solute hydrogen bonding and for “nonspecific hydration” describing the interaction of the solute with water molecules in a first-neighbor shell. The dielectric constant has been doubled (over the value used for calculations in the absence of water) to take into account the presence of solvent. Computations were carried out for the N-acetyl-N′-methylamides of the 20 naturally occurring amino acids. Conformational energy maps are compared with similar maps calculated in the absence of hydration. Minimum-energy conformations are located and compared with the corresponding minima for unhydrated peptides in terms of ordering with respect to potential energy, the dihedral angles at the minima, and the presence of intramolecular hydrogen bonds. The Boltzmann factors for various conformational regions are altered significantly on hydration in some cases. These changes can be explained in terms of differences in the hydration free energy terms for various conformations.     

Conformational study of angiotensin II

Conformational free energy calculations using an empirical potential ECEPP/3 (Empirical Conformational Energy Program for Peptides, Version 3) were carried out on angiotensin II (AII) of sequence Asp-Arg-Val-Tyr-Ile-His-Pro-Phe to find the stable conformations of the free state in the unhydrated and the hydrated states. A conformational analysis of the unhydrated state was carried out using the buildup procedure. The free energy calculation using the hydration shell model was also carried out to obtain the stable conformation of the hydrated state. The calculated stable conformations of AII in both states have a partially right-handed α-helical structure stabilized by short- and medium-range interactions. The similarity between the lowest free energy conformations of the unhydrated and hydrated states suggests that the hydration might not be important to stabilize the overall conformation of AII in a free state. The absence of any intramolecular interaction of the Tyr side chain suggests the possible interaction of this residue with the receptor. In this study, we found that the low free energy conformations contain both the parallel-plate and the perpendicular-plate geometries of the His and Phe rings, suggesting the coexistence of both conformations. © 1996 John Wiley & Sons, Inc.     

Influence of local interactions on protein structure. III. Conformational energy studies of N-acety-N′-methylamides of Gly-X and X-Gly dipeptides

Conformational energy calculations using an empirical conformational energy program for peptides (ECEPP) were carried out on 20 N-acetyl- N′-methylamides of Gly-X and X-Gly depeptides, where X = Ala, Asn, Asp, Gly, Phe, Ser, Thr, Tyr, Val, and Pro, and also of Leu-Gly. Each depeptde was found to have 25 or more low-energy minima, except Gly-Thr, which had only 11 low-energy minima because of the stable side chian-backbone hydrogen present in all low-energy conformation. As a group, the stble chain-backbone hydrogen bonds present in all low-energy conformations. As a group, the Gly-containing dipeptides were calculated in all low-energy prpensity for formation of bends than the Ala-containing depeptides. The X- Gly dipeptides were calculated to favor bends more than the Gly-X dipeptides, primarlly because of the high stability of the type II bend in X-Gly dipeptides. These results are in agreement with obseved occurrences of bends in the x-ray structures of globular proteins. The calculated conformation properties were found to be in good agreement with experimental results.     

Influence of local interactions on protein structure. I. Conformational energy studies of N-acetyl-N′-methylamides of pro-X and X-pro dipeptides

Conformational energy calculations using an Empirical Conformational Energy Program for Peptides (ECEPP) were carried out on the N-acetyl-N′-methylamides of Pro-X, where X = Ala, Asn, Asp, Gly, Leu, Phe, Ser, and Val, and of X-Pro, where X = Ala, Asn, Gly, and Pro. The conformational energy was minimized from starting conformations which included all combinations of low-energy single-residue minima and several standard bend structures. It was found that almost all resulting minima are combinations of low-energy single-residue minima, suggesting that intra residue interactions predominate in determining conformation. The calculations also indicate, however, that inter residue interactions can be important. In addition, librational entropy was found to influence the relative stabilities of some minima. Because of the existence of 10–100 low-energy minima for each dipeptide, the normalized statistical weight of an individual minimum rarely exceeds 0.3, suggesting that these dipeptides have considerable conformational flexibility and exist as statistical ensembles of low-energy structures. The propensity of each dipeptide to form bend conformations was calculated, and the results were compared with available experimental data. It was found that bends are favored in Pro-X dipeptides because ?_Pro is fixed by the pyrrolidine ring in a conformation which is frequently found in bends, but that bends are not favored in X-Pro dipeptides because interactions between the X residue and the pyrrolidine ring restrict the X residue to conformations which are not usually found in bends.     

Molecular dynamics simulations of helix folding: The effects of amino acid substitution

Molecular dynamics simulations were applied to helix folding of alanine-based synthetic peptides. A single alanine residue in the middle of the peptide was substituted with various nonpolar amino acids (leucine, isoleucine, valine, glycine or proline) to study the effect of the substitution. Unlike many other molecular dynamics simulations, nonhelical initial conformations were used in our simulations to study the folding process. An average solvent effect was included in the energy function to simplify the solvent calculation and to overcome the multiple minima problem. During the simulations, the peptides folded into helices in nanoseconds. Compact structures containing two helical segments were also observed. The calculated helical ratios of the peptides showed the same rank order as observed experimentally for the alanine-based peptides. Within a peptide, the helical ratio of each residue was calculated and a minimum was found near the center of the sequence for all peptides. The substitutions had different asymmetric effects on the helical ratios of the residues preceding and following the substitution site, indicating different helix capping preferences of the substituting amino acids. © 1997 John Wiley & Sons, Inc. Biopoly 42: 633–644, 1997     

Analysis of protein main-chain solvation as a function of secondary structure

We have analysed the hydration of main-chain carbonyl and amide groups in 24 high-resolution well-refined protein structures as a function of the secondary structure in which these polar groups occur. We find that main-chain atoms in beta-sheets are as hydrated as those in alpha-helices, with most interactions involving "free" amide and carbonyl groups that do not participate in secondary structure hydrogen bonds. The distributions of water molecules around these non-bonded carbonyl groups reflect specific steric interactions due to the local secondary structure. Approximately 20% and 4%, respectively of bonded carbonyl and amide groups interact with solvent. These include interactions with carbonyl groups on the exposed faces of alpha-helices that have been correlated previously with bending of the helix. Water molecules interacting with alpha-helices occur mainly at the amino and carbonyl termini of the helices, in which case the solvent sites maintain the hydrogen bonding by bridging between residues i and i-3 or i-4 at the amino terminus and between i and i+3 or i+4 at the carbonyl terminus. We also see a number of solvent-mediated Ncap and Ccap interactions. The water molecules interacting with beta-sheets occur mainly at the edges, in which case they extend the sheet structure, or at the ends of strands, in which case they extend the beta-ladder. In summary, the solvent networks appear to extend the hydrogen-bonding structure of the secondary structures. In beta-turns, which usually occur at the surface of a protein, exposed amide and carbonyl groups are often hydrated, especially close to glycine residues. Occasionally water molecules form a bridge between residues i and i+3 in the turn and this may provide extra stabilization.     

Design and synthesis of novel nonpolar host peptides for the determination of the 310- and α-helix compatibilities of α-amino acid buildig blocks: An assessment of α,α-disubstituted glycines

The present work describes three novel nonpolar host peptide sequences that provide a ready assessment of the 3₁₀- and α-helix compatibilities of natural and unnatural amino acids at different positions of small- to medium-size peptides. The unpolar peptides containing Ala, Aib, and a C-terminal p-iodoanilide group were designed in such a way that the peptides could be rapidly assembled in a modular fashion, were highly soluble in solvent mixtures of triflouroethanol and H₂O for CD- and two-dimensional (2D) nmr spectroscopic analyses, and showed excellent crystallinity suited for x-ray structure analysis. To validate our approach we synthesized 9-mer peptides 79a–96 (Table IV), 12-mer peptides 99–110c (Table V), and 10-mer peptides 120a–125d and 129–133 (Table VI and Scheme 8) incorporating a series of optically pure cyclic and open-chain (R)- and (S)-α,α-disubstituted glycines 1–10 (Figure 2). These amino acids are known to significantly modulate the conformations of small peptides. Based on x-ray structures of 9-mers 79a, 80, and 87 (Figures 4–7), 10-mers 124c, 131, and 132 (Figures 9–12), and 12-mer peptide 102b (Figure 13), CD spectra of all peptides recorded in acidic, neutral, and basic media and detailed 2D-nmr analyses of 9-mer peptide 86 and 12-mer 102b, several interesting conformational observations were made. Especially interesting results were obtained using the convex constraint CD analysis proposed by Fasman on 9-mer peptides 79a–d, 80, 81, 86, and 87, which allowed us to determine the relative content of 3₁₀- and α-helical conformations. These results were fully supported by the corresponding x-ray and 2D-nmr analyses. As a striking example we found that the (S)- and (R)-β-tetralin derived amino acids (R)- and (S)-1 show excellent α-helix stabilisation, more pronounced than Aib and Ala. These novel reference peptide sequences should help establish a scale for natural and unnatural amino acids concerning their intrinsic 3₁₀- and α-helix compatibilities at different positions of medium-sized peptides and thus improve our understanding in the folding processes of peptides. © 1997 John Wiley & Sons, Inc. Biopoly 42: 575–626, 1997     

An NMR study of conformations of substituted dipeptides in dodecylphosphocholine micelles: implications for drug transport

Efficient transport of intact drug (solute) across the intestinal epithelium is typically a requirement for good oral activity. In general, the membrane permeability of a solute is a complex function of its size, lipophilicity, hydrogen bond potential, charge, and conformation. In conjunction with theoretical/computational and in vitro drug transport studies, seven dipeptide (R(1)-D-Xaa-D-Phe-NHMe) homologues were each dissolved in a micellar d(38)-dodecylphosphocholine solvent system. In this homologous dipeptide series, factors such as size, lipophilicity, hydrogen-bond potential, and charge were either tightly controlled or well-characterized by other methods in order to investigate by nmr how conformational factors relate to transport. Nuclear Overhauser effect spectroscopy experiments and amide-NH-H(2)O chemical exchange rates showed that the five more lipophilic dipeptides were predominately associated with micelle, whereas the two less lipophilic analogues were not. Rotating frame nuclear Overhauser effect spectroscopy derived interproton distance restraints for each analogue, along with (3)J(HH)-derived dihedral restraints, were used in molecular dynamics/simulated annealing computations. Our results suggest that-other factors being equal-flexible dipeptides having a propensity to fold together nonpolar N- and C-terminal moieties allow greater segregation of polar and nonpolar domains and may possess enhanced transport characteristics. Dipeptides that were less flexible or that retained a less amphiphilic conformation did not have comparably enhanced transport characteristics. We suggest that these conformational/transport correlations may hold true for small, highly functionalized solutes (drugs) in general.     

Binary patterning of polar and nonpolar amino acids in the sequences and structures of native proteins

Protein sequences can be represented as binary patterns of polar (○) and nonpolar (?) amino acids. These binary sequence patterns are categorized into two classes: Class A patterns match the structural repeat of an idealized amphiphilic α-helix (3.6 residues per turn), and class B patterns match the structural repeat of an idealized amphiphilic β-strand (2 residues per turn). The difference between these two classes of sequence patterns has led to a strategy for de novo protein design based on binary patterning of polar and nonpolar amino acids. Here we ask whether similar binary patterning is incorporated in the sequences and structures of natural proteins. Analysis of the Protein Data Bank demonstrates the following. (1) Class A sequence patterns occur considerably more frequently in the sequences of natural proteins than would be expected at random, but class B patterns occur less often than expected. (2) Each pattern is found predominantly in the secondary structure expected from the binary strategy for protein design. Thus, class A patterns are found more frequently in α-helices than in β-strands, and class B patterns are found more frequently in β-strands than in α-helices. (3) Among the α-helices of natural proteins, the most commonly used binary patterns are indeed the class A patterns. (4) Among all β-strands in the database, the most commonly used binary patterns are not the expected class B patterns. (5) However, for solvent-exposed β-strands, the correlation is striking: All β-strands in the database that contain the class B patterns are exposed to solvent. (6) The bias of class A patterns for α-structure over β-structure and the bias of class B patterns for β-structure over α-structure are significant, not merely when compared to other binary patterns of polar (○) and nonpolar (?) amino acids, but also when compared to the full range of sequences in the database. The implications for the design of novel proteins are discussed.     

Criteria that discriminate between native proteins and incorrectly folded models

Various theoretical concepts, such as free energy potentials, electrostatic interaction potentials, atomic packing, solvent-exposed surface, and surface charge distribution, were tested for their ability to discriminate between native proteins and misfolded protein models. Misfolded models were constructed by introducing incorrect side chains onto polypeptide backbones: side chains of the alpha-helical hemerythrin were modeled on the beta-sheeted backbone of immunoglobulin VL domain, whereas those of the VL domain were similarly modeled on the hemerythrin backbone. CONGEN, a conformational space sampling program, was used to construct the side chains, in contrast to the previous work, where incorrect side chains were modeled in all trans conformations. Capability of the conformational search procedure to reproduce native conformations was gauged first by rebuilding (the correct) side chains in hemerythrin and the VL domain: constructs with r.m.s. differences from the x-ray side chains 2.2-2.4 A were produced, and many calculated conformations matched the native ones quite well. Incorrectly folded models were then constructed by the same conformational protocol applied to incorrect amino acid sequences. All CONGEN constructs, both correctly and incorrectly folded, were characterized by exceptionally small molecular surfaces and low potential energies. Surface charge density, atomic packing, and Coulomb formula-based electrostatic interactions of the misfolded structures and the correctly folded proteins were similar, and therefore of little interest for diagnosing incorrect folds. The following criteria clearly favored the native structures over the misfolded ones: 1) solvent-exposed side-chain nonpolar surface, 2) number of buried ionizable groups, and 3) empirical free energy functions that incorporate solvent effects.     

Helix propagation and N-cap propensities of the amino acids measured in alanine-based peptides in 40 volume percent trifluoroethanol.

The helix propagation and N-cap propensities of the amino acids have been measured in alanine-based peptides in 40 volume percent trifluoroethanol (40% TFE) to determine if this helix-stabilizing solvent uniformly affects all amino acids. The propensities in 40% TFE are compared with revised values of the helix parameters of alanine-based peptides in water. Revision of the propensities in water is the result of redefining the capping statistical weights and evaluating the helix nucleation constant with N-capping explicitly included in the helix-coil model. The propagation propensities of all amino acids increase in 40% TFE relative to water, but the increases are highly variable. In water, all beta-branched and beta-substituted amino acids are helix breakers. In 40% TFE, the propagation propensities of the nonpolar amino acids increase greatly, leaving charged and neutral polar, beta-substituted amino acids as helix breakers. Glycine and proline are strong helix breakers in both solvents. Free energy differences for helix propagation (delta delta G) between alanine and other nonpolar amino acids are twice as large in water as predicted from side-chain conformational entropies, but delta delta G values in 40% TFE are close to those predicted from side-chain entropies. This dependence of delta delta G on the solvent points to a specific role of water in determining the relative helix propensities of the nonpolar amino acids. The N-cap propensities converge toward a common value in 40% TFE, suggesting that differential solvation by water contributes to the diversity of N-cap values shown by the amino acids.     

Preferred conformations of RGDX tetrapeptides to inhibit the binding of fibrinogen to platelets

The conformational study on Arg-Gly-Asp (RGD)-containing tetrapeptides in the unhydrated and hydrated states has been carried out using the force field ECEPP/3 and the hydration shell model. The tetrapeptides studied here are H-RGDX-OH (X = Trp, Tyr, Phe, Leu, Val, Cys, Gln, and Ser), which show the inhibitory activity for binding of fibrinogen to platelets in the order of RGDW approximately equal to RGDY approximately equal to RGDF approximately equal to RGDL > RGDV > or = RGDC > or = RGDQ > or = RGDS. The backbone conformations with two C(7) backbone-to-backbone hydrogen bonds between Asp and Arg residues and between Xaa and Gly residues are in common most probable for the RGD sequence of RGDX tetrapeptides in the hydrated state. The dominant beta-turns for RGDX are found to be the types V' and IV at Gly-Asp and Asp-Xaa sequences, respectively, which are quite similar to the types II' and I (or II), respectively. However, it cannot be ruled out that the extended conformations are also remarkably feasible for RGDX tetrapeptides in water by peering the distributions of backbone conformations. These calculated results are consistent with the experimental results on RGD-containing proteins and conformationally constrained RGD-containing peptides. The reason why the RGDX becomes more potent as the side chain of the X residue is more hydrophobic may be ascribed to that the more hydrophobic is the residue X, the more populated are beta-turn structures for the Gly-Asp sequence. The hydrophobic side chain of X residue exposed to water is likely to interact with the hydrophobic region of receptor easily.     

Cations Form Sequence Selective Motifs within DNA Grooves via a Combination of Cation-Pi and Ion-Dipole/Hydrogen Bond Interactions

The fine conformational subtleties of DNA structure modulate many fundamental cellular processes including gene activation/repression, cellular division, and DNA repair. Most of these cellular processes rely on the conformational heterogeneity of specific DNA sequences. Factors including those structural characteristics inherent in the particular base sequence as well as those induced through interaction with solvent components combine to produce fine DNA structural variation including helical flexibility and conformation. Cation-pi interactions between solvent cations or their first hydration shell waters and the faces of DNA bases form sequence selectively and contribute to DNA structural heterogeneity. In this paper, we detect and characterize the binding patterns found in cation-pi interactions between solvent cations and DNA bases in a set of high resolution x-ray crystal structures. Specifically, we found that monovalent cations (Tl⁺) and the polarized first hydration shell waters of divalent cations (Mg²⁺, Ca²⁺) form cation-pi interactions with DNA bases stabilizing unstacked conformations. When these cation-pi interactions are combined with electrostatic interactions a pattern of specific binding motifs is formed within the grooves.     

Analysis of the stability of looped-out and stacked-in conformations of an adenine bulge in DNA using a continuum model for solvent and ions

A combination of conformational search, energy minimization, and energetic evaluation using a continuum solvent treatment has been employed to study the stability of various conformations of the DNA fragment d(CGCAGAA)/d(TTCGCG) containing a single adenine bulge. The extra-helical (looped-out) bulge conformation derived from a published x-ray structure and intra-helical (stacked bulge base) model structures partially based on nuclear magnetic resonance (NMR) data were used as start structures for the conformational search. Solvent-dependent contributions to the stability of the conformations were calculated from the solvent exposed molecular surface area and by using the finite difference Poisson-Boltzmann approach. Three classes (I-III) of bulge conformations with calculated low energies can be distinguished. The lowest-energy conformations were found in class I, corresponding to structures with the bulge base stacked between flanking helices, and class II, composed of structures forming a triplet of the bulge base and a flanking base pair. All extra-helical bulge structures, forming class III, were found to be less stable compared with the lowest energy structures of class I and II. The results are consistent with NMR data on an adenine bulge in the same sequence context indicating an intra-helical or triplet bulge conformation in solution. Although the total energies and total electrostatic energies of the low-energy conformations show only relatively modest variations, the energetic contributions to the stability were found to vary significantly among the classes of bulge structures. All intra-helical bulge structures are stabilized by a more favorable Coulomb charge-charge interaction but destabilized by a larger electrostatic reaction field contribution compared with all extra-helical and most triplet bulge structures. Van der Waals packing interactions and nonpolar surface-area-dependent contributions appear to favor triplet class II structures and to a lesser degree also the intra-helical stacked bulge conformations. The large conformational variation found for class III conformers might add a favorable entropic contribution to the stability of the extra-helical bulge form.     

Aib residues in peptaibiotics and synthetic sequences: analysis of nonhelical conformations

The alpha-aminoisobutyric (Aib) residue has generally been considered to be a strongly helicogenic residue as evidenced by its ability to promote helical folding in synthetic and natural sequences. Crystal structures of several peptide natural products, peptaibols, have revealed predominantly helical conformations, despite the presence of multiple helix-breaking Pro or Hyp residues. Survey of synthetic Aib-containing peptides shows a preponderance of 3(10)-, alpha-, and mixed 3(10)/alpha-helical structures. This review highlights the examples of Aib residues observed in nonhelical conformations, which fall 'primarily' into the polyproline II (P(II)) and fully extended regions of conformational space. The achiral Aib residue can adopt both left (alpha(L))- and right (alpha(R))-handed helical conformations. In sequences containing chiral amino acids, helix termination can occur by means of chiral reversal at an Aib residue, resulting in formation of a Schellman motif. Examples of Aib residues in unusual conformations are illustrated by surveying a database of Aib-containing crystal structures.     

Analysis on folding of misgurin using two-dimensional HP model

Misgurin is an antimicrobial peptide from the loach, while the hydrophobic-polar (HP) model is a way to study the folding conformations and native states in peptide and protein although several amino acids cannot be classified either hydrophobic or polar. Practically, the HP model requires extremely intensive computations, thus it has yet to be used widely. In this study, we use the two-dimensional HP model to analyze all possible folding conformations and native states of misgurin with conversion of natural amino acids according to the normalized amino acid hydrophobicity index as well as the shortest benchmark HP sequence. The results show that the conversion of misgurin into HP sequence with glycine as hydrophobic amino acid at pH 2 has 1212 folding conformations with the same native state of minimal energy -6; the conversion of glycine as polar amino acid at pH 2 has 13,386 folding conformations with three native states of minimal energy -5; the conversion of glycine as hydrophobic amino acid at pH 7 has 2538 folding conformations with three native states of minimal energy -5; and the conversion of glycine as polar amino acid at pH 7 has 12,852 folding conformations with three native states of minimal energy -4. Those native states can be ranked according to the normalized amino acid hydrophobicity index. The detailed discussions suggest two ways to modify misgurin.     

Utilization of amino acids and dipeptides by Lactobacillus plantarum from orange in nutritionally stressed conditions

Aims:  To investigate amino acid and dipeptide utilization by Lactobacillus plantarum N4 isolated from orange peel, in a nutritionally depleted medium based on MRS (Mann, Rogosa, Sharpe).
Methods and Results:  In MRS with 0·1 g l⁻¹ of meat extract and without peptone and yeast extract, growth increased when essential and stimulatory amino acids and nonessential amino acid were added to the medium. Replacement of the essential amino acid, leucine, and the nonessential amino acid, glycine, by leucyl-leucine (Leu-Leu) and/or glycyl-glycine (Gly-Gly) significantly enhanced growth. Essential amino acids were mainly consumed and the dipeptides were almost completely used at the end of growth. Leucine and glycine accumulated internally from the peptides were higher than from the free amino acids. Glucose utilization increased in the media containing dipeptides compared with the medium containing free amino acids.
Conclusions:  In a N-depleted medium, Leu-Leu and/or Gly-Gly were more effective than the respective amino acids in supporting growth of the micro-organism. The more efficient internal accumulation of glycine and especially leucine from dipeptides confirmed the ability of the strain to assimilate mainly complex nitrogen molecules rather than simple ones.
Significance and Impact of the Study:  The ability of Lact. plantarum N4 to efficiently use dipeptides could contribute to spoilage development in the natural medium of the organism, orange juice.     

Protein design simulations suggest that side-chain conformational entropy is not a strong determinant of amino acid environmental preferences

Loss of side-chain conformational entropy is an important force opposing protein folding and the relative preferences of the amino acids for being buried or solvent exposed may be partially determined by which amino acids lose more side-chain entropy when placed in the core of a protein. To investigate these preferences, we have incorporated explicit modeling of side-chain entropy into the protein design algorithm, RosettaDesign. In the standard version of the program, the energy of a particular sequence for a fixed backbone depends only on the lowest energy side-chain conformations that can be identified for that sequence. In the new model, the free energy of a single amino acid sequence is calculated by evaluating the average energy and entropy of an ensemble of structures generated by Monte Carlo sampling of amino acid side-chain conformations. To evaluate the impact of including explicit side-chain entropy, sequences were designed for 110 native protein backbones with and without the entropy model. In general, the differences between the two sets of sequences are modest, with the largest changes being observed for the longer amino acids: methionine and arginine. Overall, the identity between the designed sequences and the native sequences does not increase with the addition of entropy, unlike what is observed when other key terms are added to the model (hydrogen bonding, Lennard-Jones energies, and solvation energies). These results suggest that side-chain conformational entropy has a relatively small role in determining the preferred amino acid at each residue position in a protein.     

Hydrophobic tendencies of polar groups as a major force in molecular recognition

Proteins and nucleic acids are able to adopt their native conformation and perform their biological role only in the presence of water with which they actively interact in a mutually modifying way. Traditionally, hydrophobic effect has been considered to be the major factor stabilizing biopolymeric structures. However, solvent reorganization around polar groups is an event thermodynamically more unfavorable than solvent reorganization around nonpolar groups. Consequently, burial of polar groups with formation of complementary solute-solute hydrogen bonds out of contact with water is an energetically favorable process that also provides a major force driving macromolecular association and folding. In contrast to nonpolar groups, polar groups may form their complementary intra- or intersolute hydrogen bonds out of contact with water only provided that an appropriate solute structure has been formed with properly positioned hydrogen bond donors and acceptors. Formation of such structures is disfavored entropically and may not be possible due to steric reasons. However, the interior of a folded protein, alpha-helices and beta-sheets, double helical nucleic acid structures, and protein-ligand interfaces all provide rigid matrices where polar groups may form their complementary hydrogen bonds. For these structures, the inward drive of polar groups represents a considerable stabilizing factor.     

Conserved gating hinge in ligand- and voltage-dependent K+ channels

Ion channels open and close their pore in a process called gating. On the basis of crystal structures of two voltage-independent K(+) channels, KcsA and MthK, a conformational change for gating has been proposed whereby the inner helix bends at a glycine hinge point (gating hinge) to open the pore and straightens to close it. Here we ask if a similar gating hinge conformational change underlies the mechanics of pore opening of two eukaryotic voltage-dependent K(+) channels, Shaker and BK channels. In the Shaker channel, substitution of the gating hinge glycine with alanine and several other amino acids prevents pore opening, but the ability to open is recovered if a secondary glycine is introduced at an adjacent position. A proline at the gating hinge favors the open state of the Shaker channel as if by preventing inner helix straightening. In BK channels, which have two adjacent glycine residues, opening is significantly hindered in a graded manner with single and double mutations to alanine. These results suggest that K(+) channels, whether ligand- or voltage-dependent, open when the inner helix bends at a conserved glycine gating hinge.