Conformational changes due to vicinal glycosylation: The branched α-l-Rhap(1–2)[β-d-Galp(1–3)]-β-d-Glc1-OMe trisaccharide compared with its parent disaccharides*

Conformations of the α-l -Rhap(1-2)-β-d -Glc1-OMe and β-d -Galp(1-3)-β-d -Glc1-OMe disaccharides and the branched title trisaccharide were examined in DMSO-d₆ solution by ¹H-nmr. The distance mapping procedure was based on rotating frame nuclear Overhauser effect (NOE) constraints involving C- and O-linked protons, and hydrogen-bond constraints manifested by the splitting of the OH nmr signals for partially deuteriated samples. An “isotopomer-selected NOE” method for the unequivocal identification of mutually hydrogen-bonded hydroxyl groups was suggested. The length of hydrogen bonds thus detected is considered the only one motionally nonaveraged nmr-derived constraint. Molecular mechanics and molecular dynamics methods were used to model the conformational properties of the studied oligosaccharides. Complex conformational search, relying on a regular Φ,Ψ-grid based scanning of the conformational space of the selected glycosidic linkage, combined with simultaneous modeling of different allowed orientations of the pendant groups and the third, neighboring sugar residue, has been carried out. Energy minimizations were performed for each member of the Φ,Ψ grid generated set of conformations. Conformational clustering has been done to group the minimized conformations into families with similar values of glycosidic torsion angles. Several stable syn and anti conformations were found for the 1→2 and 1→3 bonds in the studied disaccharides. Vicinal glycosylation affected strongly the occupancy of conformational states in both branches of the title trisaccharide. The preferred conformational family of the trisaccharide (with average Φ,Ψ values of 38°, 17° for the 1→2 and 48°, 1° for the 1→3 bond, respectively) was shown by nmr to be stabilized by intramolecular hydrogen bonding between the nonbonded Rha and Gal residues. © 1998 John Wiley & Sons, Inc. Biopoly 46: 417–432, 1998     

A Monte Carlo method for generating structures of short single-stranded DNA sequences

A Monte Carlo method has been developed for generating the conformations of short single-stranded DNAs from arbitrary starting states. The chain conformers are constructed from energetically favorable arrangements of the constituent mononucleotides. Minimum energy states of individual dinucleotide monophosphate molecules are identified using a torsion angle minimizer. The glycosyl and acyclic backbone torsions of the dimers are allowed to vary, while the sugar rings are held fixed in one of the two preferred puckered forms. A total of 108 conformationally distinct states per dimer are considered in this first stage of minimization. The torsion angles within 5 kcal/mole of the global minimum in the resulting optimized states are then allowed to vary by ±10° in an effort to estimate the breadth of the different local minima. The energies of a total of 2187 (3⁷) angle combinations are examined per local conformational minimum. Finally, the energies of all dinucleotide conformers are scaled so that the populations of differently puckered sugar rings in the theoretical sample match those found in nmr solution studies. This last step is necessitated by limitations in the theoretical methods to predict DNA sugar puckering accurately. The conformer populations of the individual acyclic torsion angles in the composite dimer ensembles are found to be in good agreement with the distributions of backbone conformations deduced from nmr coupling constants and the frequencies of glycosyl conformations in x-ray crystal structures, suggesting that the low energy states are reasonable. The low energy dimer forms (consisting of 150–325 conformational states per dimer step) are next used as variables in a Monte Carlo algorithm, which generates the conformations of single-stranded d(CX_nG) chains, where X = A, T and n = 3, 4, 5. The oligonucleotides are built sequentially from the 5′ end of the chain using random numbers to select the conformations of overlapping dimer units. The simulations are very fast, involving a total of 10⁶ conformations per chain sequence. The potential errors in the buildup procedure are minimized by taking advantage of known rotational interdependences in the sugar–phosphate backbone. The distributions of oligonucleotide conformations are examined in terms of the magnitudes, positions, and orientations of the end-to-end vectors of the chains. The differences in overall flexibility and extension of the oligomers are discussed in terms of the conformations of the constituent dinucleotide steps, while the general methodology is discussed and compared with other nucleic acid model building techniques. © 1993 John Wiley & Sons, Inc.     

Theoretical studies on the conformation of β-D-N-acetyl neuraminic acid (sialic acid)

The possible conformations of sialic acid were analysed using semi-empirical potential functions. The solid state conformation has approx. 0.2 kcal/mol higher energy than the minimum energy conformation. These studies suggest that in solution sialic acid may exist preponderantly in two different conformations which differ in the orientation of the terminal hydroxymethyl group of glycerol side-chain. The present model is consistent with ¹H- and ¹³C-NMR data, but differs from the earlier models.     

DFTr studies of five- and six-residue cyclic-β(1→4) cellulosic molecules

Density functional (DFT) conformational in vacuo studies of cellobiose have shown that ?_H‐anti conformations are low in energy relative to the syn forms, while the ψ_H‐anti forms are higher in energy. Further, as the cellulosic fragments became larger than a disaccharide and new hydrogen bonding interactions between multiple residues become available, stable low energy ?_H‐anti, and ψ_H‐anti cellulosic structures became possible. To test the stability of cyclic anti‐conformations, a number of β‐linked five‐ and six‐residue molecules were created and then energy optimized in solvent (water, n‐heptane) using the implicit solvation method COSMO at the B3LYP level of theory. The created symmetric cyclic structures were without distortion. Upon optimization some cyclic conformations were found to be of low energy when compared with linear five‐ and six‐residue chains, after correcting the energy for the exclusion of a water molecule upon cyclization. It was also obvious from the hydrogen bonding network formed above and below the plane of the cyclic structure that these structures could exhibit strong synergistic tendencies. The conformational energy preferences for clockwise “c” and counter‐clockwise “r” hydroxyl groups and preference for the hydroxymethyl rotamers is described. Because these structures contain energetically unfavorable flipped conformations in water, that is, dihedral angles of ～180°/0° or ～0°/180° in ?_H/ψ_H, it is clear that the synthesis of these compounds will be challenging. © 2012 Wiley Periodicals, Inc. Biopolymers 97:568–576, 2012.     

The influence of nearest-neighboring amino acid residues on aspects of secondary structure of proteins. Attempts to locate -helices and -sheets

The influence of nearest-neighbor pairs of amino acids (n ? 1) and (n + 1) on the conformation of amino acid (n) in proteins has been studied. From experimental data on eleven proteins of known three dimensional structures, our definition of an α-helical domain in the Φ,Ψ plot has been reexamined and found to be satisfactory. On the same principle, a regular β-sheet domain has been delineated. We then revised our 20 × 20 table of frequencies of occurrences of various conformations tabulating three values: α-helical, β-sheet, and neither. These frequencies were then used to locate the helixbreaking positions in cytochrome b₅, papain, thermolysin, and calcium-binding protein. In conjuction with the helical wheel method, they were useful for predicting the locations of most α-helical segments. Similarly the β-sheet breaking positions in papain were located and most of the β-sheets found by X-ray diffraction were close to or between these positions. Data on β-sheets are extremely sparse so that extensive tests were not possible. The application of this method to abnormal hemoglobins suggested possible distortions of helices and in several instances correlated with abnormal properties of the hemoglobins and association with disease. The variable region of human immunoglobin heavy chains was found to have a very low α-helical content though β-sheet structures might exist.     

Conformational energy calculations of enzyme-substrate complexes of lysozyme. I. Energy minimization of monosaccharide and oligosaccharide inhibitors and substrates of lysozyme

Conformational energy calculations were performed on monosaccharide and oligosaccharide inhibitors and substrates of lysozyme to examine the preferred conformations of these molecules. A grid-search method was used to locate all of the low-energy conformational regions for N-acetyl-β-D -glycosamine (NAG), and energy minimization was then carried out in each of these regions. Three stable positions for the N-acetyl group have ben located, in two of which the plane of the amide unit is normal to the mean plane of the pyranosyl ring. Nine local energy minima were located for the —CH₂OH group. The positions of the two vicinal cis —OH groups are determined predominantly by interactions with either the —CH₂OH or the N-acetyl group. The most stable conformations of β-N-acetylmuramic acid (NAM) were determined from the study of the low-energy conformations of NAG. In the two stable orientations for the D -lactic acid side chain, the O—C—C′ plane (C′ being the carbon atom of the terminal carboxyl group) was found to be normal to the mean plane of the pyranosyl ring. The low-energy positions for the COOH group of NAM are determined mainly by interactions with neighboring groups. The conformational preferences of the α-anomers of NAG and NAM were also explored. The calculated conformation of the N-acetyl group for α-NAG was quite close to that determined by X-ray analysis. Two of the three lowest energy conformations of α-NAM are similar to the corresponding conformations of the β-anomer. A third low-energy structure, which has a hydrogen bond from the NH of the N-acetyl group to the C?O of the lactic acid group, corresponds very closely to the X-ray structure of this molecule. The preferred conformations of the disaccharides NAG–NAG, NAM–NAG and NAG–NAM were also investigated. Two preferred orientations of the reducing pyranosyl ring relative to the nonreducing ring were found for all of these disaccharides, both of which are close to the extended conformation. In one of these conformations, a hydrogen bond can form between the OH group attached to C₃ of the reducing sugar and the ring oxygen of the preceding residue. Each conformation can be stabilized further by a hydrogen bond between the CH₂OH (donor) of residue i + 1 and the C?O of residue i (acceptor). The interactions that determine conformations for all oligosaccharides containing both NAG and NAM are shown to be exclusively intraresidue and nearest neighbor interactions, so that it is possible to predict all stable conformations of oligosaccharides containing NAG and NAM in any sequence.     

DFTMD studies of β-cellobiose: conformational preference using implicit solvent

Previous DFT in vacuo studies on the conformational preferences for cellobiose showed that upon optimization the φ_H-anti conformations were of lower energy than the syn forms. Upon optimization using an implicit solvation method, COSMO, the syn or observed form was still not predicted to be of lower energy than the φ_H-anti form, even though optimization after addition of several explicit water molecules did show a relative energy difference favoring the syn form. In order to examine the predictive ability of COSMO on this carbohydrate, constant energy dynamics, DFTMD, simulations were carried out on low energy syn and φ_H-anti conformations with and without COSMO included during the dynamics. The resulting analysis confirmed that when COSMO is included in the dynamics, the syn conformations become energetically favored over the φ_H-anti forms suggesting that both solvent and entropy play roles in dictating the solution conformation of cellobiose. Analysis of the dynamic runs includes distributions of selected dihedral angles versus time, conformational transitions, and populations of some quasi-planar, boat, skew forms during the simulations.     

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.     

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.     

Structural basis of hierarchical multiple substates of a protein. I: Introduction

A computer experiment of protein dynamics is carried out, which consists of two steps: (1) A Monte Carlo simulation of thermal fluctuations in the native state of a globular protein, bovine pancreatic trypsin inhibitor; and (2) a simulation of the quick freezing of fluctuating conformations into energy minima by minimization of the energy of a number of conformations sampled in the Monte Carlo simulation. From the analysis of results of the computer experiment is obtained the following picture of protein dynamics: multiple energy minima exist in the native state, and they are distributed in clusters in the conformational space. The dynamics has a hierarchical structure which has at least two levels. In the first level, dynamics is restricted within one of the clusters of minima. In the second, transitions occur among the clusters. Local parts of a protein molecule, side chains and local main chain segments, can take multiple locally stable conformations in the native state. Many minima result from combinations of these multiple local conformations. The hierarchical structure in the dynamics comes from interactions among the local parts. Protein molecules have two types of flexibility, each associated with elastic and plastic deformations, respectively.     

The chbG Gene of the Chitobiose (chb) Operon of Escherichia coli Encodes a Chitooligosaccharide Deacetylase

The chb operon of Escherichia coli is involved in the utilization of the β-glucosides chitobiose and cellobiose. The function of chbG (ydjC), the sixth open reading frame of the operon that codes for an evolutionarily conserved protein is unknown. We show that chbG encodes a monodeacetylase that is essential for growth on the acetylated chitooligosaccharides chitobiose and chitotriose but is dispensable for growth on cellobiose and chitosan dimer, the deacetylated form of chitobiose. The predicted active site of the enzyme was validated by demonstrating loss of function upon substitution of its putative metal-binding residues that are conserved across the YdjC family of proteins. We show that activation of the chb promoter by the regulatory protein ChbR is dependent on ChbG, suggesting that deacetylation of chitobiose-6-P and chitotriose-6-P is necessary for their recognition by ChbR as inducers. Strains carrying mutations in chbR conferring the ability to grow on both cellobiose and chitobiose are independent of chbG function for induction, suggesting that gain of function mutations in ChbR allow it to recognize the acetylated form of the oligosaccharides. ChbR-independent expression of the permease and phospho-β-glucosidase from a heterologous promoter did not support growth on both chitobiose and chitotriose in the absence of chbG, suggesting an additional role of chbG in the hydrolysis of chitooligosaccharides. The homologs of chbG in metazoans have been implicated in development and inflammatory diseases of the intestine, indicating that understanding the function of E. coli chbG has a broader significance.     

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.     

Chitosan Supplementation and Fecal Fat Excretion in Men

Objective : Few weight loss supplements are clinically tested for efficacy, yet their proliferation continues. Chitosan‐based supplements are sold as fat trappers and fat magnets. They purportedly block fat absorption and cause weight loss without food restriction. We quantified the in vivo effect of a chitosan product on fat absorption. Research Methods and Procedures : Participants (n = 15) consumed five meals per day for 12 days. Energy intake was not restricted. Participants consumed no supplements during a 4‐day control period and two capsules five times per day (4.5 g chitosan/d), 30 minutes before each meal, during a 4‐day supplement period. All feces were collected from days 2 to 12. Oral charcoal markers permitted division of the feces into two periods. The two fecal pools were analyzed for fat content. Results : Participants were male, 26.3 ± 5.9 years old, BMI of 25.6 ± 2.3 kg/m². Subjects consumed 133 ± 23 g of fat/d and 12.91 ± 1.79 MJ/d (3084 ± 427 kcal/d). Individual meals averaged 26.3 ± 9.3 g of fat. With chitosan supplementation at 10 capsules/day, fecal fat excretion increased by 1.1 ± 1.8 g/d (p = 0.02), from 6.1 ± 1.2 to 7.2 ± 1.8 g/d. Discussion : The effect of chitosan on fat absorption is clinically negligible. Far from being a fat trapper, at 0.11 ± 0.18 g of fat trapped per 0.45‐g capsule or 1.1 g (9.9 kcal) fat trapped per day, this product would have no significant effect on energy balance. The fat trapping claims associated with chitosan are unsubstantiated.     

The flexibility of the LeaLex Tumor Associated Antigen central fragment studied by systematic and stochastic searches as well as dynamic simulations

The solution conformational behavior of the Tumor-Associated Carbohydrate Antigen Le^aLe^x central fragment: methyl α-l-fucopyranosyl-(1→4)-2-acetamido-2-deoxy-β-d-glucopyranosyl-(1→3)-β-d-galactopyranoside was studied using three computational methods: a rigid systematic search as implemented in Sybyl, a stochastic search as implemented in MOE2004, and dynamics simulations using the SANDER module of AMBER9. Our results illustrate the complementarity of these methods to identify energetically relevant conformations and flexible linkages. In particular, the β-GlcNAc-(1→3)-Gal linkage was shown to be extremely flexible adopting a wide range of orientations around two energy minima. The modeling results were validated by comparison of theoretical distances, derived from the simulations, with experimental measurements obtained from 1D selective ROESY buildup curves on the synthetic fragment.     

Properties of Chitosanase from Bacillus cereus S1

Chitosanase from Bacillus cereus S1 was purified, and the enzymatic properties were investigated. The molecular weight was estimated to 45,000 on SDS-PAGE. Optimum pH was about 6, and stable pH in the incubation at 40°C for 60 min was 6–11. This chitosanase was stable in alkaline side. Optimum temperature was around 60°C, and enzyme activity was relatively stable below 60°C. The degradations of colloidal chitosan and carboxymethyl cellulose (CMC) were about 30 and 20% relative to the value of soluble chitosan, respectively, but colloidal chitin and crystalline cellulose were not almost hydrolyzed. On the other hand, S1 chitosanase adsorbed on colloidal chitin completely and by about 50% also on crystalline cellulose, in contrast to colloidal chitosan, which it did not adsorb. S1 chitosanase finally hydrolyzed 100% N-deacetylated chitosan (soluble state) to chitobiose (27.2%), chitotriose (40.6%), and chitotetraose (32.2%). In the hydrolysis of various chitooligosaccharides, chitobiose and chitotriose were not hydrolyzed, and chitotetraose was hydrolyzed to chitobiose. Chitobiose and chitotriose were released from chitopentaose and chitohexaose. From this specificity, it was hypothesized that the active site of S1 chitosanase recognized more than two glucosamine residues posited in both sides against splitting point for glucosamine polymer. Received: 8 June 1999 / Accepted: 20 July 1999     

Enzymatic hydrolysis of cellulose to glucose lung immobilized β-glucosidase

The production of sugars by enzymatic hydrolysis of cellulose is a multistep process which includes conversion of the intermediate cellobiose to glucose by β-glucosidase. Aside from its role as an intermediate, cellobiose inhibits the endoglucanase components of typical cellulase enzyme systems. Because these enzyme systems often contain insufficient concentrations of β-glucosidase to prevent accumulation of inhibitory cellobiose, this research investigated the use of supplemental immobilized β-glucosidase to increase yield of glucose. Immobilized β-glucosidase from Aspergillus phoenicis was produced by sorption at controlled-pore alumina with about 90% activity retention. The product lost only about 10% of the original activity during an on-stream reaction period of 500 hr with cellobiose as substrate; maximum activity occurred near pH 3.5 and the apparent activation energy was about 11 kcal/mol. The immobilized β-glucosidase was used together with Trichoderma reesei cellulase to hydrolyze cellulosic materials, such as Solka Floc, corn stove and exploded wood. Increased yields of glucose and greater conversions of cellobiose of glucose were observed when the reaction systems contained supplemental immobilized β-glucosidase.     

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.     

Noncontact dipole effects on channel permeation. II. Trp conformations and dipole potentials in gramicidin A

The four Trp dipoles in the gramicidin A (gA) channel modulate channel conductance, and their side chain conformations should therefore be important, but the energies of different conformations are unknown. A conformational search for the right-handed helix based on molecular mechanics in vacuo yielded 46 conformations within 20 kcal/mol of the lowest energy conformation. The two lowest energy conformations correspond to the solid-state and solution-state NMR conformations, suggesting that interactions within the peptide determine the conformation. For representative conformations, the electrostatic potential of the Trp side chains on the channel axis was computed. A novel application of the image-series method of. Biophys. J. 9:1160-1170) was introduced to simulate the polarization of bulk water by the Trp side chains. For the experimentally observed structures, the CHARm toph19 potential energy (PE) of a cation in the channel center is -1.65 kcal/mol without images. With images, the PE is -1.9 kcal/mol, demonstrating that the images further enhance the direct dipole effect. Nonstandard conformations yielded less favorable PEs by 0.4-1.1 kcal/mol.