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

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.

A lattice relaxation algorithm for three-dimensional Poisson-Nernst-Planck theory with application to ion transport through the gramicidin A channel

A lattice relaxation algorithm is developed to solve the Poisson-Nernst-Planck (PNP) equations for ion transport through arbitrary three-dimensional volumes. Calculations of systems characterized by simple parallel plate and cylindrical pore geometries are presented in order to calibrate the accuracy of the method. A study of ion transport through gramicidin A dimer is carried out within this PNP framework. Good agreement with experimental measurements is obtained. Strengths and weaknesses of the PNP approach are discussed.

Noncontact dipole effects on channel permeation. I. Experiments with (5F-indole)Trp13 gramicidin A channels

Gramicidin A (gA), with four Trp residues per monomer, has an increased conductance compared to its Phe replacement analogs. When the dipole moment of the Trp13 side chain is increased by fluorination at indole position 5 (FgA), the conductance is expected to increase further. gA and FgA conductances to Na+, K+, and H+ were measured in planar diphytanoylphosphatidylcholine (DPhPC) or glycerylmonoolein (GMO) bilayers. In DPhPC bilayers, Na+ and K+ conductances increased upon fluorination, whereas in GMO they decreased. The low ratio in the monoglyceride bilayer was not reversed in GMO-ether bilayers, solvent-inflated or -deflated bilayers, or variable fatty acid chain monoglyceride bilayers. In both GMO and DPhPC bilayers, fluorination decreased conductance to H+ but increased conductance in the mixed solution, 1 M KCl at pH 2.0, where K+ dominates conduction. Eadie-Hofstee plot slopes suggest similar destabilization of K+ binding in both lipids. Channel lifetimes were not affected by fluorination in either lipid. These observations indicate that fluorination does not change the rotameric conformation of the side chain. The expected difference in the rate-limiting step for transport through channels in the two bilayers qualitatively explains all of the above trends.

Heterodimer formation and crystal nucleation of gramicidin D

The linear pentadecapeptide antibiotic gramicidin D is a heterogeneous mixture of six components. Precise refinements of three-dimensional structures of naturally occurring gramicidin D in crystals obtained from methanol, ethanol, and n-propanol demonstrate the unexpected presence of stable left-handed antiparallel double-helical heterodimers that vary with the crystallization solvent. The side chains of Trp residues in the three structures exhibit sequence-specific patterns of conformational preference. Tyr substitution for Trp at position 11 appears to favor beta ribbon formation and stabilization of the antiparallel double helix that acts as a template for gramicidin folding and nucleation of different crystal forms. The fact that a minor component in a heterogeneous mixture influences aggregation and crystal nucleation has potential applications to other systems in which anomalous behavior is exhibited by aggregation of apparently homogeneous materials, such as the enigmatic behavior of prion proteins.

The conducting form of gramicidin A is a right-handed double-stranded double helix

The linear pentadecapeptide antibiotic, gramicidin D, is a naturally occurring product of Bacillus brevis known to form ion channels in synthetic and natural membranes. The x-ray crystal structures of the right-handed double-stranded double-helical dimers (DSDH) reported here agree with 15N-NMR and CD data on the functional gramicidin D channel in lipid bilayers. These structures demonstrate single-file ion transfer through the channels. The results also indicate that previous crystal structure reports of a left-handed double-stranded double-helical dimer in complex with Cs+ and K+ salts may be in error and that our evidence points to the DSDH as the major conformer responsible for ion transport in membranes.

Free energy profiles for H+ conduction along hydrogen-bonded chains of water molecules

The molecular mechanism for proton conduction along hydrogen-bonded chains, or "proton wires," is studied with free energy simulations. The complete transport of a charge along a proton wire requires two complementary processes: 1) translocation of an excess proton (propagation of an ionic defect), and 2) reorientation of the hydrogen-bonded chain (propagation of a bonding defect). The potential of mean force profile for these two steps is computed in model systems comprising a single-file chain of nine dissociable and polarizable water molecules represented by the PM6 model of Stillinger and co-workers. Results of molecular dynamics simulations with umbrella sampling indicate that the unprotonated chain is preferably polarized, and that the inversion of its total dipole moment involves an activation free energy of 8 kcal/mol. In contrast, the rapid translocation of an excess H+ across a chain extending between two spherical solvent droplets is an activationless process. These results suggest that the propagation of a bonding defect constitutes a limiting step for the passage of several protons along single-file chains of water molecules, whereas the ionic translocation may be fast enough to occur within the lifetime of transient hydrogen-bonded water chains in biological membranes.

Recent Advances in the High Resolution Structures of Bacterial Channels: Gramicidin A

Gramicidin is a polypeptide antibiotic which forms dimeric channels specific for the transport of monovalent cations across membranes. It adopts several different conformations, most notably double helical (pore) and helical dimer (channels) forms, which have very different structural and functional characteristics. This review focuses on recent high resolution structure determinations of both the pore and channel forms of the molecule by X-ray crystallographic and/or NMR spectroscopic techniques. It discusses the structural consequences of binding ions and the location of ion binding sites and how the structures are related to the conductance properties of the molecule. This relatively simple molecule is probably the best characterized ion channel (both structurally and functionally) and has, to date, been the principal proving-ground for many of our ideas about the molecular nature of ion conduction in membranes. Copyright 1998 Academic Press.

Models and simulations of ion channels and related membrane proteins

The past year has seen major advances in our understanding of ion channels, resulting from molecular dynamics simulations and modelling studies. Simulations of gramicidin have revealed that proton conduction along a water wire is limited by the dynamics of water reorientation. Plausible models are now available for a number of other channels, including alamethicin, the influenza A virus M2 protein, and the pore domains of the nicotinic acetylcholine receptor and Kv channels. Molecular dynamics simulations and continuum calculations have revealed some of the subtleties of the interactions between transmembrane helices and their lipid bilayer environment.

Using theory and simulation to understand permeation and selectivity in ion channels

It is clear that the function of ion channels must flow from their structure. With recent advances in computational power and methodology, it appears feasible to correlate structure to ion channel permeation at an atomistically detailed level of description. The overall strategy is to structure the calculations in a hierarchy, ranging from coarse-grained thermodynamic and kinetic descriptions to fine-grained molecular dynamics descriptions with atomic detail. Each level of description is connected to the others by appropriate statistical mechanical theory. The coarse-grained descriptions can be correlated directly with electrophysiological experiment. The fine-grained descriptions are used to parameterize the coarse-grained descriptions and to describe the permeation process at the most detailed level. This strategy has so far had varying degrees of success. It has successfully described water permeation through lipid bilayers and gramicidin channels. It has revealed the essential events of ion permeation through gramicidin channels at an atomistically detailed level. The role of channel protein motions in permeation has been elucidated. However, it appears that force fields used to describe molecular dynamics must be refined further to achieve completely accurate predictions of the permeation of such small ions as sodium. Channels with more complex structure and more multiion occupancy than gramicidin pose major computational challenges with respect to sampling protein conformations and ion distributions involved in the permeation process. Possible approaches to meeting these challenges are discussed.

High-resolution polypeptide structure in a lamellar phase lipid environment from solid state NMR derived orientational constraints

BACKGROUND

Solid-state nuclear magnetic resonance (NMR) spectroscopy provides novel structural constraints from uniformly aligned samples. These orientational constraints orient specific atomic sites with respect to the magnetic field direction and the unique molecular axis of alignment. Solid-state NMR is uniquely and ideally suited for providing such structural constraints on polypeptides and proteins in a lamellar phase lipid environment. Membrane protein structure represents a great challenge for structural biologists; a new approach for characterizing high resolution three-dimensional structure in such an environment is needed.

RESULTS

The optimal use of orientational constraints for defining three-dimensional structures is demonstrated with the elucidation of the gramicidin A channel structure at high resolution. Initial structures are refined against both the experimental constraints and the CHARMM energy using a novel simulated-annealing protocol to define torsion angle solutions with an error bar of approximately +/- 5 degrees.

CONCLUSIONS

This analysis results in the determination of a high-resolution, time averaged structure of gramicidin A obtained in a lipid bilayer environment above the gel-to-liquid crystalline phase transition temperature. It is demonstrated that solid-state NMR can be used to establish polypeptide, and potentially protein, structures in such an environment. Furthermore, this high-resolution structure is demonstrated to provide new insights into polypeptide function. For the gramicidin A channel the roles of the indole groups that facilitate ion transport and details of the cation solvation environment provided by the amide oxygens are characterized.

Influence of the membrane potential on the free energy of an intrinsic protein

A modified Poisson-Boltzmann equation is developed from statistical mechanical considerations to describe the influence of the transmembrane potential on macromolecular systems. Using a Green's function formalism, the electrostatic free energy of a protein associated with the membrane is expressed as the sum of three terms: a contribution from the energy required to charge the system's capacitance, a contribution corresponding to the interaction of the protein charges with the membrane potential, and a contribution corresponding to a voltage-independent reaction field free energy. The membrane potential, which is due to the polarization interface, is calculated in the absence of the protein charges, whereas the reaction field is calculated in the absence of transmembrane potential. Variations in the capacitive energy associated with typical molecular processes are negligible under physiological conditions. The formulation of the theory is closely related to standard algorithms used to solve the Poisson-Boltzmann equation and only small modifications to current source codes are required for its implementation. The theory is illustrated by examining the voltage-dependent membrane insertion of a simple polyalanine alpha-helix and by computing the electrostatic potential across a 60-A-diameter sphere meant to represent a large intrinsic protein.