Noncontact dipole effects on channel permeation. V. Computed potentials for fluorinated gramicidin

Experimental and theoretical calculations indicate that the dipole moment of the four Trp side chains in gramicidin A (gA) channels modify channel conductance through long-range electrostatic interactions. Electrostatic ion/side-chain interaction energies along the channel were computed with CHARMM using ab initio atom charges for native and 4-, 5-, or 6-fluorinated Trp side chains. The bulk water reaction to the polar side chains was included using the method of images as implemented by, and channel waters in idealized structures were included. Ion/Trp interaction energies were approximately -0.6 kcal/mol throughout the channel for all four of the native Trp pairs. Channel waters produced a modest reduction in the magnitude of interactions, essentially offsetting images representing the bulk water outside the channel. The effects of side-chain fluorination depended on ring position and, to a lesser extent, residue number. Compared with native Trp, 5-fluorination reduces the translocation barrier with minor effects on the exit barrier. In contrast, 6-fluorination primarily reduces exit barrier. 4-Fluorination produces a more complex double-well energy profile. Effects of measured side-chain movements resulting from fluorination or change in lipid bilayer were negligible whereas thermal side chain librations cause large effects, especially in the region of the ion-binding sites.     

DFT study of alpha- and beta-D-mannopyranose at the B3LYP/6-311++G** level

Thirty-five conformations of alpha- and beta-d-mannopyranose, the C-2 substituted epimer of glucopyranose, were geometry optimized using the density functional (B3LYP), and basis set (6-311++G**). Full geometry optimization was performed on the hydroxymethyl rotamers (gg/gt/tg) and an analytical hessian program was used to calculate the harmonic vibrational frequencies, zero point energy, enthalpy, and entropy. The lowest energy conformation investigated is the beta-tg in the (4)C(1) chair conformation. The in vacuo calculations showed little energetic preference for either the alpha or beta anomer for mannopyranose in the (4)C(1) chair conformation. Results are compared to similar glucopyranose calculations in vacuo where the alpha anomer is approximately 1kcal/mol lower in electronic energy than the beta anomer. In the case of the generally higher energy (1)C(4) chair conformations, one low-energy, low-entropy beta-gg-(1)C(4) chair conformation was identified that is within approximately 1.4kcal/mol of the lowest energy (4)C(1) conformation of mannopyranose. Other (1)C(4) chair conformations in our investigation are approximately 2.9-7.9kcal/mol higher in overall energy. Many of the (3,O)B, B(3,O), (1,4)B, and B(1,4) boat forms passed through transitions without barriers to (1)S(3), (5)S(1), (1)S(5) skew forms with energies between approximately 3.6 and 8.9kcal/mol higher in energy than the lowest energy conformation of mannopyranose. Boat forms were found that remained stable upon gradient optimization. As with glucopyranose, the orientation and interaction of the hydroxy groups make a significant contribution to the conformation/energy relationship in vacuo.     

Assignment of side-chain conformation using adiabatic energy mapping, free energy perturbation, and molecular dynamic simulations

NMR spectroscopic analysis of the C-terminal Kunitz domain fragment (alpha3(VI)) from the human alpha3-chain of type VI collagen has revealed that the side chain of Trp21 exists in two unequally populated conformations. The major conformation (M) is identical to the conformation observed in the X-ray crystallographic structure, while the minor conformation (m) cannot structurally be resolved in detail by NMR due to insufficient NOE data. In the present study, we have applied: (1) rigid and adiabatic mapping, (2) free energy simulations, and (3) molecular dynamic simulations to elucidate the structure of the m conformer and to provide a possible pathway of the Trp21 side chain between the two conformers. Adiabatic energy mapping of conformations of the Trp21 side chain obtained by energy minimization identified two energy minima: One corresponding to the conformation of Trp21 observed in the X-ray crystallographic structure and solution structure of alpha3(VI) (the M conformation) and the second corresponding to the m conformation predicted by NMR spectroscopy. A transition pathway between the M and m conformation is suggested. The free-energy difference between the two conformers obtained by the thermodynamic integration method is calculated to 1.77+/-0.7 kcal/mol in favor of the M form, which is in good agreement with NMR results. Structural and dynamic properties of the major and minor conformers of the alpha3(VI) molecule were investigated by molecular dynamic. Essential dynamics analysis of the two resulting 800 ps trajectories reveals that when going from the M to the m conformation only small, localized changes in the protein structure are induced. However, notable differences are observed in the mobility of the binding loop (residues Thr13-Ile18), which is more flexible in the m conformation than in the M conformation. This suggests that the reorientation of Trp2 might influence the inhibitory activity against trypsin, despite the relative large distance between the binding loop and Trp21.     

DFT study of alpha- and beta-D-allopyranose at the B3LYP/6-311++G ** level of theory

One hundred and two conformations of alpha- and beta-D-allopyranose, the C-3 substituted epimer of glucopyranose, were geometry optimized using the density functional, B3LYP, and the basis set, 6-311++G **. Full geometry optimization was performed on different ring geometries and on the hydroxymethyl rotamers (gg/gt/tg). Analytically derived Hessians were used to calculate zero point energy, enthalpy, and entropy. The lowest energy and free energy conformation found is the alpha-tg(g-)-4C1-c conformation, which is only slightly higher in electronic (approximately 0.2 kcal/mol) and free energy than the lowest energy alpha-D-glucopyranose. The in vacuo calculations showed a small (approximately 0.3 kcal/mol) energetic preference for the alpha- over the beta-anomer for allopyranose in the 4C1 conformation, whereas in the 1C4 conformation a considerable (approximately 1.6 kcal/mol) energetic preference for the beta- over the alpha-anomer for allopyranose was encountered. The results are compared to previous aldohexose calculations in vacuo. Boat and skew forms were found that remained stable upon gradient optimization although many starting boat conformations moved to other skew forms upon optimization. As found for glucose, mannose, and galactose the orientation and interaction of the hydroxyl groups make the most significant contributions to the conformation/energy relationship in vacuo. A comparison of different basis sets and density functionals is made in the Discussion section, confirming the appropriateness of the level of theory used here.     

New developments of the electrostatically driven monte carlo method: Test on the membrane-bound portion of melittin

The electrostatically driven Monte Carlo (EDMC) method has been greatly improved by adding a series of new features, including a procedure for cluster analysis of the accepted conformations. This information is used to guide the search for the global energy minimum. Alternative procedures for generating perturbed conformations to sample the conformational space were also included. These procedures enhance the efficiency of the method by generating a larger number of low-energy conformations. The improved EDMC method has been used to explore the conformational space of a 20-residue polypeptide chain whose sequence corresponds to the membrane-bound portion of melittin. The ECEPP/3 (Empirical Conformational Energy Program for Peptides) algorithm was used to describe the conformational energy of the chain. After an exhaustive search involving 14 independent runs, the lowest energy conformation (LEC) (−91.0 kcal/mol) of the entire study was encountered in four of the runs, while conformations higher in energy by no more than 1.8 kcal/mol were found in the remaining runs with the exception of one of them (run 8). The LEC is identical to the conformation found recently by J. Lee, H.A. Scheraga, and S. Rackovsky [(1998) “Conformational Analysis of the 20-Residue Membrane-Bound Portion of Melittin by Conformational Space Annealing,” Biopolymers, Vol. 46, pp. 103–115] as the lowest energy conformation obtained in their study using the conformational space annealing method. These results suggest that this conformation corresponds to the global energy minimum of the ECEPP/3 potential function for this specific sequence; it also appears to be the conformation of lowest free energy. © 1998 John Wiley & Sons, Inc. Biopoly 46: 117–126, 1998     

B3LYP/6-311++G** study of alpha- and beta-D-glucopyranose and 1,5-anhydro-D-glucitol: 4C1 and 1C4 chairs, (3,O)B and B(3,O) boats, and skew-boat conformations

Geometry optimization, at the B3LYP/6-311++G** level of theory, was carried out on 4C1 and 1C4 chairs, (3,O)B and B(3,O) boats, and skew-boat conformations of alpha- and beta-D-glucopyranose. Similar calculations on 1,5-anhydro-D-glucitol allowed examination of the effect of removal of the 1-hydroxy group on the energy preference of the hydroxymethyl rotamers. Stable minimum energy boat conformers of glucose were found, as were stable skew boats, all having energies ranging from approximately 4-15 kcal/mol above the global energy 4C1 chair conformation. The 1C4 chair electronic energies were approximately 5-10 kcal/mol higher than the 4C1 chair, with the 1C4 alpha-anomers being lower in energy than the beta-anomers. Zero-point energy, enthalpy, entropy, and relative Gibbs free energies are reported at the harmonic level of theory. The alpha-anomer 4C1 chair conformations were found to be approximately 1 kcal/mol lower in electronic energy than the beta-anomers. The hydroxymethyl gt conformation was of lowest electronic energy for both the alpha- and beta-anomers. The glucose alpha/beta anomer ratio calculated from the relative free energies is 63/37%. From a numerical Hessian calculation, the tg conformations were found to be approximately 0.4-0.7 kcal/mol higher in relative free energy than the gg or gt conformers. Transition-state barriers to rotation about the C-5-C-6 bond were calculated for each glucose anomer with resulting barriers to rotation of approximately 3.7-5.8 kcal/mol. No energy barrier was found for the path between the alpha-gt and alpha-gg B(3,O) boat forms and the equivalent 4C1 chair conformations. The alpha-tg conformation has an energy minimum in the 1S3 twist form. Other boat and skew-boat forms are described. The beta-anomer boats retained their starting conformations, with the exception of the beta-tg-(3,O)B boat that moved to a skew form upon optimization.     

Tryptophan-47 rotational isomerization in variant-3 scorpion neurotoxin. A combination thermodynamic perturbation and umbrella sampling study.

A combination thermodynamic perturbation and umbrella sampling study predicts two free energy wells for the rotational isomerization of the variant-3 scorpion neurotoxin tryptophan-47 indole side chain. One well has the indole side chain in the crystallographic orientation; the other has the indole rotated approximately 220 degrees to form a new conformation with a relative free energy of 3 +/- 2 kcal/mol. The activation barrier is 8.5 kcal/mol from the crystallographic well, from which transition state theory predicts a rate of escape of 2 x 10(5) s-1. Correlations in the displacements of side chains neighboring tryptophan-47 and the isomerization reaction coordinate last up to 20 ps. Favorable conditions of experimental verification are discussed.     

The effect of the amino-acid side chains on the energy profiles for ion transport in the gramicidin A channel

Computations on the energy profiles for Na+ in the gramicidin A (GA) channel have been extended by introducing the effect, previously neglected, of the amino acid side chains of GA, fixed in their most stable conformations. The calculations have been performed in two approximations: 1) with the ethanolamine tail fixed in its most stable conformation, 2) with the tail allowed to optimize its conformation upon the progression of the ion. In both approximations the overall shape of the energy profile is very similar to that obtained in the absence of the side chains. One observes, however, a general lowering of the profile upon the adjunction of the side chains. The analysis of the factors responsible for this energy lowering indicates that it is due essentially to the electrostatic and polarisation components of the interaction which interplay differently, however, in the different parts of the channel. A particular role is attributed in this respect to the tryptophan residues of GA. The role of the 4 tryptophans present, Trp 15, 13, 11 and 9, is individualized by stripping of one of them at a time. The strongest effect on the energy deepening is due to Trp 13 and is particularly prominent in the entrance zone at 14.5A from the center of the channel. The result indicates the possibility of investigating theoretically the effect on the energy profiles of the substitution of the "natural" side chain by others.     

Conformational analysis of the 20-residue membrane-bound portion of melittin by conformational space annealing

The conformational space of the 20-residue membrane-bound portion of melittin has been investigated extensively with the conformational space annealing (CSA) method and the ECEPP/3 (Empirical Conformational Energy Program for Peptides) algorithm. Starting from random conformations, the CSA method finds that there are at least five different classes of conformations, within 4 kcal/mol, which have distinct backbone structures. We find that the lowest energy conformation of this peptide from previous investigations is not the global minimum-energy conformation (GMEC); but it belongs to the second lowest energy class of the five classes found here. In four independent runs, one conformation is found repeatedly as the lowest energy conformation of the peptide (two of the four lowest energy conformations are identical; the other two have essentially identical backbone conformations but slightly different side-chain conformations). We propose this conformation, whose energy is lower than that found previously by 1.9 kcal/mol, as the GMEC of the ECEPP/3 force field. The structure of the proposed GMEC is less helical and more compact than the previous one. It appears that the CSA method can find several classes of conformations of a 20-residue peptide starting from random conformations utilizing only its amino acid sequence information. The proposed GMEC has also been found with a modified electrostatically driven Monte Carlo method [D. R. Ripoll, A. Liwo, and H.A. Scheraga (1998) “New Developments of the Electrostatically Driven Monte Carlo Method: Test on the Membrane-Bound Portion of Melittin,” Biopolymers, Vol. 46, pp. 117–126]. © 1998 John Wiley & Sons, Inc. Biopoly 46: 103–115, 1998     

Molecular dynamics simulations of Trp side-chain conformational flexibility in the gramicidin A channel

Gramicidin A (gA) is prototypical peptide antibiotic and a model ion channel former. Configured in the solid-state NMR beta(6.5)-helix channel conformation, gA was subjected to 1-ns molecular dynamics (MD) gas phase simulations using the all-atom charmm22 force field to ascertain the conformational stability of the Trp side chains as governed by backbone and neighboring side-chain contacts. Three microcanonical trajectories were computed using different initial atomic velocities for each of twenty different initial structures. For each set, one of the four Trp side chains in each monomer was initially positioned in one of the five non-native conformations (A. E. Dorigo et al., Biophysical Journal, 1999, Vol. 76, 1897-1908), the other Trps being positioned in the native state, o1. In three additional control simulations, all Trps were initiated in the native conformation. After equilibration, constraints were removed and subsequent conformational changes of the initially constrained Trp were measured. The chi(1) was more flexible than chi(2.1). The energetically optimal orientation, o1 (Dorigo et al., 1999), was the most stable in all four Trp positions (9, 11, 13, 15) and remained unchanged for the entire 1 ns simulation in 19 of 24 trials. Changes in chi(1) from each of the 5 suboptimal states occur readily. Two of the non-native conformations reverted readily to o1, whereas the other three converted to an intermediate state, i2. There were frequent interconversions between i2 and o1. We speculate that experimentally observed Trp stability is caused by interactions with the lipid-water interface, and that stabilization of one of the suboptimal conformations in gA, such as i2, by lipid headgroups could produce a secondary, metastable conformational state. This could explain recent experimental studies of differences in the channel conductance dispersity between gA and a Trp-to-Phe gA analog, gramicidin M (gM, J. C. Markham et al., Biochimica et Biophysica Acta, 2001, Vol. 1513, 185-192).     

Conformational analysis of d-tubocurarine: Implications for minimal dimensions of its binding site within ion channels

All the minimum-energy conformations of d-tubocurarine were calculated by the method of molecular mechanics. The energy was minimized from 413 closed forms of the 18-member ring. The set of minimum-energy conformations includes 10 forms with energies less than 6 kcal/mol from the most stable one. Among the four lowest minimum-energy conformations, two forms correspond to those known from X-ray studies, whereas two conformations were not detected experimentally earlier. The flexibility of d-tubocurarine was estimated by calculating six paths of interconversion between the four lowest minimum-energy conformations. Using a molecular graphics technique, it was found that the most extended minimum-energy conformation of d-tubocurarine may fit in an ion channel of a rectangular profile of 8.7 × 11.2 Å, while one tetrahydroisoquinoline head may fit a profile as small as 6.9 × 11.0 Å. A possible model of d-tubocurarine location within the ion channel of the neuronal nicotinic acetylcholine receptor is suggested.     

The Effect of the Amino-Acid Side Chains on the Energy Profiles for Ion Transport in the Gramicidin A Channel

Abstract

Computations on the energy profiles for Na⁺ in the gramicidin A (GA) channel have been extended by introducing the effect, previously neglected, of the amino acid side chains of GA, fixed in their most stable conformations. The calculations have been performed in two approximations: 1) with the ethanolamine tail fixed in its most stable conformation, 2) with the tail allowed to optimize its conformation upon the progression of the ion. In both approximations the overall shape of the energy profile is very similar to that obtained in the absence of the side chains. One observes, however, a general lowering of the profile upon the adjunction of the side chains. The analysis of the factors responsible for this energy lowering indicates that it is due essentially to the electrostatic and polarisation components of the interaction which interplay differently, however, in the different parts of the channel. A particular role is attributed in this respect to the tryptophan residues of GA. The role of the 4 tryptophans present, Trp 15, 13, 11 and 9, is individualized by stripping of one of them at a time. The strongest effect on the energy deepening is due to Trp 13 and is particularly prominent in the entrance zone at 14.5Å from the center of the channel. The result indicates the possibility of investigating theoretically the effect on the energy profiles of the substitution of the “natural” side chain by others.     

Conformational energy analysis of melanostatin

Conformational energy calculations were carried out on the hypothalamic hormone melanostatin, a tripeptide with the primary structure H-L-Pro-L-Leu-Gly-NH2. The calculated lowest energy conformation was a type II beta bend, very similar to that reported in an X-ray crystal study. This conformation, however, was only one of 109 low-energy structures (less than or equal to 3 kcal/mol above the global minimum), indicating that the molecule in solution exists as an ensemble of conformations and is very flexible, in agreement with relaxation data from n.m.r. measurements. A statistical analysis yielded an average end-to-end distance of 6.8 A and a bend probability of 0.62, suggesting that, in nonpolar solvents, bend structures predominate within the statistical ensemble. The statistical analysis, however, also yielded a probability of only 0.11 for the occurrence of a 4 leads to 1 hydrogen bond. Hence, the calculations show that, although bend conformations predominate, bends would be difficult to observe in solution if the experiments were designed only to detect 4 leads to 1 hydrogen bonds.     

Studying excited states of proteins by NMR spectroscopy.

Protein structure is inherently dynamic, with function often predicated on excursions from low to higher energy conformations. For example, X-ray studies of a cavity mutant of T4 lysozyme, L99A, show that the cavity is sterically inaccessible to ligand, yet the protein is able to bind substituted benzenes rapidly. We have used novel relaxation dispersion NMR techniques to kinetically and thermodynamically characterize a transition between a highly populated (97%, 25 degrees C) ground state conformation and an excited state that is 2.0 kcal mol(-1) higher in free energy. A temperature-dependent study of the rates of interconversion between ground and excited states allows the separation of the free energy change into enthalpic (Delta H = 7.1 kcal mol(-1)) and entropic (T Delta S = 5.1 kcal mol(-1), 25 degrees C) components. The residues involved cluster about the cavity, providing evidence that the excited state facilitates ligand entry.     

Survey of conformational role of ester bonds in a cyclic depsipeptide. Study on cyclo(-L-Ala-L-Hmb-)2 by energy calculation and NMR spectroscopy.

The effect of ester bond on the conformation of peptide molecule was studied by designing and synthesizing a model tetradepsipeptide cyclo(-L-Ala-L-Hmb-)2 and by analyzing the conformation both theoretically and experimentally. Theoretical analysis showed that both ester and peptide bonds in the calculated low-energy conformations within 3 kcal/mol of the global minimum take a trans but distorted configuration. The distortion is larger in ester bonds than in peptide bonds. Further, the four carbonyls project from one side of the plane of the cyclic backbone, whereas the side chains project from the other side. These results are consistent with the experimental results obtained by NMR measurement; first, the coupling constant deduced from 1H-NMR species in DMSO-d6 is consistent with the dihedral angles of the calculated low-energy conformations; second, results of NOE measurement can reproduce the calculated configuration of the carbonyls and side chains. From the consistency between theoretical and experimental results, it is concluded that this model tetradepsipeptide takes an all-trans backbone conformation in solution and this backbone conformation is stabilized by large distortion in the ester bond, which compensates the strain resulted from the 12-membered cyclic backbone structure consisting only of L-residues.     

Flexibility of bovine pancreatic trypsin inhibitor

The native conformation of a protein may be expressed in terms of the dihedral angles, phi's and psi's for the backbone, and kappa's for the side chains, for a given geometry (bond lengths and bond angles). We have developed a method to obtain the dihedral angles for a low-energy structure of a protein, starting with the X-ray structure; it is applied here to examine the degree of flexibility of bovine pancreatic trypsin inhibitor. Minimization of the total energy of the inhibitor (including nonbonded, electrostatic, torsional, hydrogen bonding, and disulfide loop energies) yields a conformation having a total energy of -221 kcal/mol and a root mean square deviation between all atoms of the computed and experimental structures of 0.63 A. The optimal conformation is not unique, however, there being at least two other conformations of low-energy (-222 and -220 kcal/mol), which resemble the experimental one (root mean square deviations of 0.66 and 0.64 A, respectively). These three conformations are located in different positions in phi, psi space, i.e., with a total deviation of 81 degrees, 100 degrees and 55 degrees from each other (with a root mean square deviation of several degrees per dihedral angle from each other). The nonbonded energies of the backbones, calculated along lines in phi, psi space connecting these three conformations, are all negative, without any intervening energy barriers (on an energy contour map in the phi, psi plane). Side chains were attached at several representative positions in this plane, and the total energy was minimized by varying the kappa's. The energies were of approximately the same magnitude as the previous ones, indicating that the conformation of low energy is flexible to some extent in a restricted region of phi, psi space. Interestingly, the difference delta phi i+1 in phi i+1 for the (i + 1)th residue from one conformation to another is approximately the same as -delta psi i for the ith residue; i.e., the plane of the peptide group between the ith and (i + 1)th residues re-orient without significant changes in the positions of the other atoms. The flexibility of the orientations of the planes of the peptide groups is probably coupled in a cooperative manner to the flexibility of the positions of the backbone and side-chain atoms.     

Monte Carlo loop refinement and virtual screening of the thyroid-stimulating hormone receptor transmembrane domain

Metropolis Monte Carlo (MMC) loop refinement has been performed on the three extracellular loops (ECLs) of rhodopsin and opsin-based homology models of the thyroid-stimulating hormone receptor transmembrane domain, a class A type G protein-coupled receptor. The Monte Carlo sampling technique, employing torsion angles of amino acid side chains and local moves for the six consecutive backbone torsion angles, has previously reproduced the conformation of several loops with known crystal structures with accuracy consistently less than 2?Å. A grid-based potential map, which includes van der Waals, electrostatics, hydrophobic as well as hydrogen-bond potentials for bulk protein environment and the solvation effect, has been used to significantly reduce the computational cost of energy evaluation. A modified sigmoidal distance-dependent dielectric function has been implemented in conjunction with the desolvation and hydrogen-bonding terms. A long high-temperature simulation with 2?kcal/mol repulsion potential resulted in extensive sampling of the conformational space. The slow annealing leading to the low-energy structures predicted secondary structure by the MMC technique. Molecular docking with the reported agonist reproduced the binding site within 1.5?Å. Virtual screening performed on the three lowest structures showed that the ligand-binding mode in the inter-helical region is dependent on the ECL conformations.     

DFT study of alpha- and beta-D-galactopyranose at the B3LYP/6-311++G** level of theory

Forty-one conformations of alpha- and beta-d-galactopyranose were geometry optimized using the B3LYP density functional and 6-311++G** basis set. Full geometry optimization was performed on different ring geometries and different hydroxymethyl rotamers (gg/gt/tg). Analytically derived Hessians were used to calculate zero point energy, enthalpy, and entropy. The lowest energy and free-energy conformation found is the alpha-gg-(4)C(1)-c chair conformation, which is of lower electronic and free energy than the lowest energy alpha-d-glucopyranose conformer because of favorable hydrogen-bonding interactions. The in vacuo calculations showed considerable ( approximately 2.2kcal/mol) energetic preference for the alpha over the beta anomer for galactopyranose in both the (4)C(1) and (1)C(4) chair conformations. Results are compared to glucopyranose and mannopyranose calculations in vacuo. Boat and skew-boat forms were found that remained stable upon gradient optimization, although many starting conformations moved to other boat forms upon optimization. As with glucopyranose and mannopyranose, the orientation and interaction of the hydroxyl groups make the most significant contributions to the conformation-energy relationship in vacuo.     

Crystal structure and ligand-binding studies of a screened peptide complexed with streptavidin.

The thermodynamic binding parameters and crystal structure for streptavidin-peptide complexes where the peptide sequences were obtained by random screening methods are reported. The affinities between streptavidin and two heptapeptides were determined by titrating calorimetric methods [Phe-Ser-His-Pro-Gln-Asn-Thr, Ka = 7944 (+/- 224) M-1, delta G degrees = -5.32 (+/- 0.01) kcal/mol, and delta H degrees = -19.34 (+/- 0.48) kcal/mol; His-Asp-His-Pro-Gln-Asn-Leu, Ka = 3542 (+/- 146) M-1, delta G degrees = -4.84 (+/- 0.03) kcal/mol, and delta H degrees = -19.00 (+/- 0.64) kcal/mol]. The crystal structure of streptavidin complexed with one of these peptides has been determined at 2.0-A resolution. The peptide (Phe-Ser-His-Pro-Gln-Asn-Thr) binds in a turn conformation with the histidine, proline, and glutamine side chains oriented inward at the biotin-binding site. A water molecule is immobilized between the histidine and glutamine side chains of the peptide and an aspartic acid side chain of the protein. Although some of the residues that participate in binding biotin also interact with the screened peptide, the peptide adopts an alternate method of utilizing binding determinants in the biotin-binding site of streptavidin.     

Noncontact dipole effects on channel permeation. III. Anomalous proton conductance effects in gramicidin

Proton transport on water wires, of interest for many problems in membrane biology, is analyzed in side-chain analogs of gramicidin A channels. In symmetrical 0.1 N HCl solutions, fluorination of channel Trp(11), Trp-(13), or Trp(15) side chains is found to inhibit proton transport, and replacement of one or more Trps with Phe enhances proton transport, the opposite of the effects on K(+) transport in lecithin bilayers. The current-voltage relations are superlinear, indicating that some membrane field-dependent process is rate limiting. The interfacial dipole effects are usually assumed to affect the rate of cation translocation across the channel. For proton conductance, however, water reorientation after proton translocation is anticipated to be rate limiting. We propose that the findings reported here are most readily interpreted as the result of dipole-dipole interactions between channel waters and polar side chains or lipid headgroups. In particular, if reorientation of the water column begins with the water nearest the channel exit, this hypothesis explains the negative impact of fluorination and the positive impact of headgroup dipole on proton conductance.