Search

This paper reports on a simulation of a gramicidin channel inserted into a fluid phase DMPC bilayer with 100 lipid molecules. Two lipid molecules per leaflet were removed to insert the gramicidin, so the resulting preparation had 96 lipid molecules and 3209 water molecules. Constant surface tension boundary conditions were employed. Like previous simulations with a lower lipid/gramicidin ratio (Woolf, T. B., and B. Roux. 1996. Proteins: Struct., Funct., Genet. 24:92-114), it is found that tryptophan-water hydrogen bonds are more common than tryptophan-phospholipid hydrogen bonds. However, one of the tryptophan NH groups entered into an unusually long-lived hydrogen bonding pattern with two glycerol oxygens of one of the phospholipid molecules. Comparisons were made between the behavior of the lipids adjacent to the channel with those farther away. It was found that hydrocarbon chains of lipids adjacent to the channel had higher-order parameters than those farther away. The thickness of the lipid bilayer immediately adjacent to the channel was greater than it was farther away. In general, the lipids adjacent to the membrane had similar orientations to those seen by Woolf and Roux, while those farther away had similar orientations to those pertaining before the insertion of the gramicidin. A corollary to this observation is that the thickness of the hydrocarbon region adjacent to the gramicidin was much thicker than what other studies have identified as the "hydrophobic length" of the gramicidin channel.

A gramicidin channel in a fluid phase DMPC bilayer with excess lipid and water has been simulated. By use of the formal correspondence between diffusion and random walk, a permeability for water through the channel was calculated, and was found to agree closely with the experimental results of Rosenberg and Finkelstein (Rosenberg, P.A., and A. Finkelstein. 1978. J. Gen. Physiol. 72:327-340; 341-350) for permeation of water through gramicidin in a phospholipid membrane. By using fluctuation analysis, components of resistance to permeation were computed for movement through the channel interior, for the transition step at the channel mouth where the water molecule solvation environment changes, and for the process of diffusion up to the channel mouth. The majority of the resistance to permeation appears to occur in the transition step at the channel mouth. A significant amount is also due to structurally based free energy barriers within the channel. Only small amounts are due to local friction within the channel or to diffusive resistance for approaching the channel mouth.

This paper reports on a coupled approach to determining the structure of the gramicidin A ion channel, utilizing solid state nuclear magnetic resonance (NMR) of isotopically labeled gramicidin channels aligned parallel to the magnetic field direction, and molecular dynamics (MD). MD computations using an idealized right-handed beta-helix as a starting point produce a refined molecular structure that is in excellent agreement with atomic resolution solid state NMR data. The data provided by NMR and MD are complementary to each other. When applied in a coordinated manner they provide a powerful approach to structure determination in molecular systems not readily amenable to x-ray diffraction.

Molecular dynamics simulations have been done on a system consisting of the polypeptide membrane channel former gramicidin, plus water molecules in the channel and caps of waters at the two ends of the channel. In the absence of explicit simulation of the surrounding membrane, the helical form of the channel was maintained by artificial restraints on the peptide motion. The characteristic time constant of the artificial restraint was varied to assess the effect of the restraints on the channel structure and water motions. Time-correlation analysis was done on the motions of individual channel waters and on the motions of the center of mass of the channel waters. It is found that individual water molecules confined in the channel execute higher frequency motions than bulk water, for all degrees of channel peptide restraint. The center-of-mass motion of the chain of channel waters (which is the motion that is critical for transmembrane transport, due to the mandatory single filing of water in the channel) does not exhibit these higher frequency motions. The mobility of the water chain is dramatically reduced by holding the channel rigid. Thus permeation through the channel is not like flow through a rigid pipe; rather permeation is facilitated by peptide motion. For the looser restraints we used, the mobility of the water chain was not very much affected by the degree of restraint. Depending on which set of experiments is considered, the computed mobility of our water chain in the flexible channel is four to twenty times too high to account for the experimentally measured resistance of the gramicidin channel to water flow. From this result it appears likely that the peptide motions of an actual gramicidin channel embedded in a lipid membrane may be more restrained than in our flexible channel model, and that these restraints may be a significant modulator of channel permeability. For the completely rigid channel model the "trapping" of the water molecules in preferred positions throughout the molecular dynamics run precludes a reasonable assessment of mobility, but it seems to be quite low.

This paper shows how Brownian motion theory can be used to analyze features of individual ion trajectories in channels as calculated by molecular dynamics, and that its use permits more precise determinations of diffusion coefficients than would otherwise be possible. We also show how a consideration of trajectories of single particles can distinguish between effects due to the magnitude of the diffusion coefficient and effects due to barriers and wells in the potential profile, effects which can not be distinguished by consideration of average fluxes.

The electrodiffusion equations were solved for the one-ion channel both by the analytical method due to Levitt and also by Brownian dynamic simulations. For both types of calculations equilibration of ion distribution between the bath and the ends of the channel was assumed. Potential profiles were found that give good fits to published data on Na+ permeation of gramicidin channels. The data were best fit by profiles that have no relative energy maximum at the mouth of the channel. This finding suggests that alignment of waters or channel charged groups inside the channel in response to an ion's approach may provide an energetically favorable situation for entry sufficient to overcome the energy required for removing bulk waters of hydration. An alternative possibility is that the barrier to ion entry is situated outside the region restricted to single-ion occupancy. Replacement of valine with more polar amino acids at the No. 1 location was found to correspond to a deepening of the potential minima near the channel mouths, an increase in height of the central barrier to ion translocation across the channel, and possibly a reduction in the mobility of the ion-water complex in the channel. The Levitt theory was extended to calculate passage times for ions to cross the channel and the blocking effects of ions that entered the channel but didn't cross. These quantities were also calculated by the Brownian dynamics method.

Theoretical studies of ion channels address several important questions. The mechanism of ion transport, the role of water structure, the fluctuations of the protein channel itself, and the influence of structural changes are accessible from these studies. In this paper, we have carried out a 70-ps molecular dynamics simulation on a model structure of gramicidin A with channel waters. The backbone of the protein has been analyzed with respect to the orientation of the carbonyl and the amide groups. The results are in conformity with the experimental NMR data. The structure of water and the hydrogen bonding network are also investigated. It is found that the water molecules inside the channel act as a collective chain; whereas the conformation in which all the waters are oriented with the dipoles pointing along the axis of the channel is a preferred one, others are also accessed during the dynamics simulation. A collective coordinate involving the channel waters and some of the hydrogen bonding peptide partners is required to describe the transition of waters from one configuration to the other.