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.     

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.     

The gramicidin A channel: The energy profile calculated for Na+ in the presence of water with inclusion of the flexibility of the ethanolamine tail

The effect of the conformational freedom of the ethanolamine tail of gramicidin A on the energy profile for the transfer of Na⁺, computed in the presence of water, shows an appreciable lowering of the minimum at 10.5 Å, and a splitting of the entrance barrier. The deep energy region at the channel mouth remains, however, the deepest one and contains a site of strong interaction of the ion with the Trp 11 carbonyl, at about 12 Å from the center.     

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.     

Modeling of alpha-MSH conformations with implicit solvent.

A conformational search for the most probable structures of the hormone alpha-MSH in aqueous solution was performed in order to help determine the structural features necessary for biological activity. The free-energy surface was modeled using methods from integral equation theory, and high-temperature molecular dynamics was used to enhance conformational sampling. Families of low free-energy structures have been found. The minimum energy structure shows a stable beta-turn conformation in the putative message region that is stabilized by a salt bridge between Glu5 and Lys11. The orientation of the side chains reflects the amphiphilic nature of the peptide, and a close interaction between the side chains of the His6, Phe7 and Trp9 was observed. Several structural features observed in the minimum energy structure agree well with experimental results. The conformational features led to a hypothesis of a receptor-hormone interaction model in which the hydrophobic side chains of Phe7 and Trp9 interact with the transmembrane portion of the human melanocortin (MC1) receptor. Also, the positively charged side chain of Arg8 and the imidazole side chain of His6 may interact with the negatively charged portions of the receptor which may even be on the receptor's extracellular loops.     

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).     

Energy profiles in the acetylcholine receptor (AChR) channel. The MII-helix model and the role of the remaining helices

It is demonstrated by theoretical computations that no favorable energy profile for cation transfer can be obtained in a model of the AChR channel constructed with the sole five MII helices of the inner wall. A favorable profile is obtained upon including the effect of the remaining helices of the five subunits. The decisive role, for the exit of the ion, of the charged residues situated at the N-terminal of the MII segments, established before, is underlined further. The role of the other elements of the channel wall (peptide carbonyl oxygens, hydrocarbon residues and polar side chains) is analyzed.     

Interaction of K+ ion with the solvated gramicidin A transmembrane channel.

Using Urry's gramicidin A (GA) atomic coordinates and ab into calculations, the interaction energies of a K+ ion with GA are examined. From these energies the values of the fitting parameters are obtained for 6-12-1 atom-atom pair potentials. The potential of the GA channel as experienced by the ion is analyzed in detail. An energy profile of the K+ ion in the GA channel is obtained by analyzing iso-energy maps. Using Monte Carlo simulations, the energy profiles of the K+ ion with the solvated GA channel are analyzed and the hydration structures in the presence of the K+ ion are studied.     

Energy Profile of Cs+ in Gramicidin A in the Presence of Water. Problem of the Ion Selectivity of the Channel

Abstract

The effect of water present at the mouth and inside the channel of Gramicidin A on the energy profile calculated for a caesium ion is determined. The total optimal interaction energy computed for the system GA-Cs⁺-(22 waters) leads to an energy profile characterized by a deep minimum at 11Å followed by an entrance energy barrier of 7 Kcal/mol expanding until 9 Å from the center. After this point, a second minimum less deep than the previous one is observed, itself followed by a central barrier. The shape of the profile at the entrance is governed by the balance between the progressive desolvation process of the ion and the increase of favorable hydrogen bond interactions implying both the water molecules and GA. The comparison of this energy profile with that obtained in vacuo shows that the presence of water molecules does not modify the pathway of the ion which, owing to its size, is constrained essentially to remain on the channel axis. The comparison Na⁺ versus Cs⁺ indicates that although the phenomena involved are globally the same, differences between the two profiles appear due firstly to the difference in the affinity of the two ions for water and secondly to their respective size. This last difference implies that the number of water molecules present in the interior of the channel during the cation progression is reduced roughly by one in the case of caesium.

The desolvation barrier computed for Cs⁺ is half the corresponding value for Na⁺, a result in agreement with the observed selectivity.     

Structure, dynamics, and selectivity of the sodium channel blocker mu-conotoxin SIIIA

mu-SIIIA, a novel mu-conotoxin from Conus striatus, appeared to be a selective blocker of tetrodotoxin-resistant sodium channels in frog preparations. It also exhibited potent analgesic activity in mice, although its selectivity profile against mammalian sodium channels remains unknown. We have determined the structure of mu-SIIIA in aqueous solution and characterized its backbone dynamics by NMR and its functional properties electrophysiologically. Consistent with the absence of hydroxyprolines, mu-SIIIA adopts a single conformation with all peptide bonds in the trans conformation. The C-terminal region contains a well-defined helix encompassing residues 11-16, while residues 3-5 in the N-terminal region form a helix-like turn resembling 3 10-helix. The Trp12 and His16 side chains are close together, as in the related conotoxin mu-SmIIIA, but Asn2 is more distant. Dynamics measurements show that the N-terminus and Ser9 have larger-magnitude motions on the subnanosecond time scale, while the C-terminus is more rigid. Cys4, Trp12, and Cys13 undergo significant conformational exchange on microsecond to millisecond time scales. mu-SIIIA is a potent, nearly irreversible blocker of Na V1.2 but also blocks Na V1.4 and Na V1.6 with submicromolar potency. The selectivity profile of mu-SIIIA, including poor activity against the cardiac sodium channel, Na V1.5, is similar to that of the closely related mu-KIIIA, suggesting that the C-terminal regions of both are critical for blocking neuronal Na V1.2. The structural and functional characterization described in this paper of an analgesic mu-conotoxin that targets neuronal subtypes of mammalian sodium channels provides a basis for the design of novel analogues with an improved selectivity profile.     

Energy profile of Cs+ in gramicidin A in the presence of water. Problem of the ion selectivity of the channel

The effect of water present at the mouth and inside the channel of Gramicidin A on the energy profile calculated for a caesium ion is determined. The total optimal interaction energy computed for the system GA-Cs+-(22 waters) leads to an energy profile characterized by a deep minimum at 11 A followed by an entrance energy barrier of 7 Kcal/mol expanding until 9 A from the center. After this point, a second minimum less deep than the previous one is observed, itself followed by a central barrier. The shape of the profile at the entrance is governed by the balance between the progressive desolvation process of the ion and the increase of favorable hydrogen bond interactions implying both the water molecules and GA. The comparison of this energy profile with that obtained in vacuo shows that the presence of water molecules does not modify the pathway of the ion which, owing to its size, is constrained essentially to remain on the channel axis. The comparison Na+ versus Cs+ indicates that although the phenomena involved are globally the same, differences between the two profiles appear due firstly to the difference in the affinity of the two ions for water and secondly to their respective size. This last difference implies that the number of water molecules present in the interior of the channel during the cation progression is reduced roughly by one in the case of caesium. The desolvation barrier computed for Cs+ is half the corresponding value for Na+, a result in agreement with the observed selectivity.     

Analysis of the ion transfer through the channel of 9,11,13,15-phenylalanylgramicidin A

The behavior of an analogue of gramicidin A in which all four tryptophanyl residues are substituted by phenylalanyl and which shows a strong voltage effect on the single channel conductance is analyzed on the basis of a 'three-barrier--two-site' model. It is shown that in the gramicidin family the side chains of some amino acids, in spite of their location, which point outside the channel can play a major role in the binding of ions in the channel and thus can significantly modify the energy profile of the channel.     

The structure,cation binding,transport, and conductance of Gly15-gramicidin A incorporated into SDS micelles and PC/PG vesicles

To further investigate the effect of single amino acid substitution on the structure and function of the gramicidin channel, an analogue of gramicidin A (GA) has been synthesized in which Trp(15) is replaced by Gly in the critical aqueous interface and cation binding region. The structure of Gly(15)-GA incorporated into SDS micelles has been determined using a combination of 2D-NMR spectroscopy and molecular modeling. Like the parent GA, Gly(15)-GA forms a dimeric channel composed of two single-stranded, right-handed beta(6.3)-helices joined by hydrogen bonds between their N-termini. The replacement of Trp(15) by Gly does not have a significant effect on backbone structure or side chain conformations with the exception of Trp(11) in which the indole ring is rotated away from the channel axis. Measurement of the equilibrium binding constants and Delta G for the binding of monovalent cations to GA and Gly(15)-GA channels incorporated into PC vesicles using (205)Tl NMR spectroscopy shows that monovalent cations bind much more weakly to the Gly(15)-GA channel entrance than to GA channels. Utilizing the magnetization inversion transfer NMR technique, the transport of Na(+) ions through GA and Gly(15)-GA channels incorporated into PC/PG vesicles has been investigated. The Gly(15) substitution produces an increase in the activation enthalpy of transport and thus a significant decrease in the transport rate of the Na(+) ion is observed. The single-channel appearances show that the conducting channels have a single, well-defined structure. Consistent with the NMR results, the single-channel conductances are reduced by 30% and the lifetimes by 70%. It is concluded that the decrease in cation binding, transport, and conductance in Gly(15)-GA results from the removal of the Trp(15) dipole and, to a lesser extent, the change in orientation of Trp(11).     

Molecular structure of the cell-attachment protein of reovirus: correlation of computer-processed electron micrographs with sequence-based predictions.

The receptor-recognition interaction that initiates reovirus infection is mediated by the sigma 1 protein, located at the vertices of the icosahedral virion. We have applied computer-based image-averaging techniques to electron micrographs of negatively stained preparations of sigma 1 purified from virions (serotype 2 Jones). Combining these results with inferences based on the amino acid sequence has led to a molecular model in which the overall folding of the chains is described; its conformation embodies motifs, coiled-coil alpha-helices and nodular multichain elements rich in beta-sheets, previously detected in the corresponding proteins of other viruses, but with some novel variations. Sigma 1 is a filamentous lollipop-shaped molecule with an overall length of approximately 48 nm; it has a flexible "tail," approximately 40 nm long by 4 to 6 nm wide, terminating at its distal end in a globular "head," approximately 9.5 nm in diameter. The purified protein is a tetramer (4 by 50 kilodaltons) consisting of two similarly oriented dimers bonded side by side and in register. For each chain, a cluster of hydrophobic residues at its amino terminus resides at the proximal end of the tail; next, an alpha-helical domain (residues 25 to 172) participates in a two-chained coiled coil, 22 nm long, with two such coiled coils pairing laterally to form the proximal half of the tail. The remainder of the tail (residues 173 to approximately 316) is less uniform in width and is expected to be rich in beta-sheet; the interdimer bonding is evidently sustained through this portion of the molecule. Finally, the globular head consists of the carboxy-terminal domains (which contain the receptor-binding sites) folded into compact globular conformations; in appropriate side views, the head is resolved into two subunits, presumably contributed by the respective dimers. This model for how the four sigma 1 polypeptide chains are threaded in parallel through the fiber is supported by the observed match between an empirical curvature profile, which identifies the locations of relatively flexible sites along the tail, and the flexibility profile predicted on the basis of the model. Appraisal of the interactions that stabilize the coiled coils suggests that (i) the alpha-helices are individually only marginally stable, a property that may be of significance with regard to the retracted conformation in which sigma 1 is accommodated in the intact virion, and (ii) the predominant interactions between the two coiled coils are likely to involve hydrogen bonding between patches of uncharged residues.     

Use of synthetic gramicidins in the determination of channel structure and mechanism

The syntheses of (1-13C) Trp9 gramicidin A (GA), (1-13C) Trp11 GA, (1-13C) Trp13 GA, (1-13C) Trp15 GA, and D . Leu2 GA are verified by means of high performance liquid chromatography, carbon-13 nuclear magnetic resonance, circular dichroism and characterization of transport properties. The use of these and other synthetic gramicidins is discussed in terms of determining ion binding sites within the channel, helix sense of the channel, the basis of monovalent vs divalent cation selectivity, and means of modulating channel conductance.     

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.     

The planar conformation of a strained proline ring: a QM/MM study

QM and QM/MM energy calculations have been carried out on an atomic resolution structure of liganded triosephosphate isomerase (TIM) that has an active site proline (Pro168) in a planar conformation. The origin of the planarity of this proline has been identified. Steric interactions between the atoms of the proline ring and a tyrosine ring (Tyr166) on one side of the proline prevent the ring from adopting the up pucker (chi1 is approximately -30 degrees), while the side chain of a nearby alanine (Ala171) forbids the down pucker (chi1 is approximately +30 degrees). To obtain a proline conformation that is in agreement with the experimentally observed planar state, a quantum system of sufficient size is required and should at least include the nearby side chains of Tyr166, Ala171, and Glu129 to provide enough stabilization. It is argued that the current force fields for structure optimization do not describe strained protein fragments correctly. The proline is part of a catalytic loop that closes upon ligand binding. Comparison of the proline conformation in different TIM X-ray structures, indicates that in the closed conformation of TIM the proline is planar or nearly planar, while in the open conformation it is down puckered. This suggests that the planarity possibly plays a role in the overall catalytic cycle of TIM, presumable acting as a reservoir of energy that becomes available upon loop opening.     

The binding site of sodium in the gramicidin A channel: comparison of molecular dynamics with solid-state NMR data.

The location of the main binding site for sodium in the gramicidin A (GA) channel was investigated with molecular dynamics simulations, using an atomic model of the channel embedded in a fully hydrated dimyristoyl phosphatidycholine (DMPC) bilayer. Twenty-four separate simulations in which a sodium was restrained at different locations along the channel axis were generated. The results are compared with carbonyl 13C chemical shift anisotropy solid-state NMR experimental data previously obtained with oriented GA:DMPC samples. Predictions are made for other solid-state NMR properties that could be observed experimentally. The combined information from experiment and simulation strongly suggests that the main binding sites for sodium are near the channel's mouth, approximately 9.2 A from the center of the dimer channel. The 13C chemical shift anisotropy of Leu10 is the most affected by the presence of a sodium ion in the binding site. In the binding site, the sodium ion is lying off-axis, making contact with two carbonyl oxygens and two single-file water molecules. The main channel ligand is provided by the carbonyl group of the Leu10-Trp11 peptide linkage, which exhibits the largest deviation from the ion-free channel structure. Transient contacts with the carbonyl group of Val8 and Trp15 are also present. The influence of the tryptophan side chains on the channel conductance is examined based on the current information about the binding site.     

Side-chain ionization states in a potassium channel

KcsA is a bacterial K+ channel that is gated by pH. Continuum dielectric calculations on the crystal structure of the channel protein embedded in a low dielectric slab suggest that side chains E71 and D80 of each subunit, which lie adjacent to the selectivity filter region of the channel, form a proton-sharing pair in which E71 is neutral (protonated) and D80 is negatively charged at pH 7. When K+ ions are introduced into the system at their crystallographic positions the pattern of proton sharing is altered. The largest perturbation is for a K+ ion at site S3, i.e., interacting with the carbonyls of T75 and V76. The presence of multiple K+ ions in the filter increases the probability of E71 being ionized and of D80 remaining neutral (i.e., protonated). The ionization states of the protein side chains influence the potential energy profile experienced by a K+ ion as it is translated along the pore axis. In particular, the ionization state of the E71-D80 proton-sharing pair modulates the shape of the potential profile in the vicinity of the selectivity filter. Such reciprocal effects of ion occupancy on side-chain ionization states, and of side-chain ionization states on ion potential energy profiles will complicate molecular dynamics simulations and related studies designed to calculate ion permeation energetics.     

Conduction properties of KcsA measured using brownian dynamics with flexible carbonyl groups in the selectivity filter

In the narrow segment of an ion conducting pathway, it is likely that a permeating ion influences the positions of the nearby atoms that carry partial or full electronic charges. Here we introduce a method of incorporating the motion of charged atoms lining the pore into Brownian dynamics simulations of ion conduction. The movements of the carbonyl groups in the selectivity filter of the KcsA channel are calculated explicitly, allowing their bond lengths, bond angles, and dihedral angels to change in response to the forces acting upon them. By systematically changing the coefficients of bond stretching and of angle bending, the carbon and oxygen atoms can be made to fluctuate from their fixed positions by varying mean distances. We show that incorporating carbonyl motion in this way does not alter the mechanism of ion conduction and only has a small influence on the computed current. The slope conductance of the channel increases by approximately 25% when the root mean-square fluctuations of the carbonyl groups are increased from 0.01 to 0.61 A. The energy profiles and the number of resident ions in the channel remain unchanged. The method we utilized here can be extended to allow the movement of glutamate or aspartate side chains lining the selectivity filters of other ionic channels.