Anionic phospholipid interactions with the potassium channel KcsA: simulation studies

Molecular dynamics (MD) simulations have been used to unmask details of specific interactions of anionic phospholipids with intersubunit binding sites on the surface of the bacterial potassium channel KcsA. Crystallographic data on a diacyl glycerol fragment at this site were used to model phosphatidylethanolamine (PE), or phosphatidylglycerol (PG), or phosphatidic acid (PA) at the intersubunit binding sites. Each of these models of a KcsA-lipid complex was embedded in phosphatidyl choline bilayer and explored in a 20 ns MD simulation. H-bond analysis revealed that in terms of lipid-protein interactions PA > PG > PE and revealed how anionic lipids (PG and PA) bind to a site provided by two key arginine residues (R(64) and R(89)) at the interface between adjacent subunits. A 27 ns simulation was performed in which KcsA (without any lipids initially modeled at the R(64)/R(89) sites) was embedded in a PE/PG bilayer. There was a progressive specific increase over the course of the simulation in the number of H-bonds of PG with KcsA. Furthermore, two specific PG binding events at R(64)/R(89) sites were observed. The phosphate oxygen atoms of bound PG formed H-bonds to the guanidinium group of R(89), whereas the terminal glycerol H-bonded to R(64). Overall, this study suggests that simulations can help identify and characterize sites for specific lipid interactions on a membrane protein surface.     

Quantum calculations on water in the KcsA channel cavity with permeant and non-permeant ions

Different ions in the pore of the KcsA channel behave differently, and we relate this to their solvation. We show that the selectivity is dependent, in part, on the solvation in the cavity (sometimes referred to as the vestibule, it is the region containing water molecules between the intracellular gate and the selectivity filter at the extracellular end of the pore). We have shown earlier that potassium is more dependent at the upper end of the cavity region on solvation by the threonines there, while sodium ion has more water molecules as ligands. In addition, sodium ion is placed asymmetrically, while potassium is nearly exactly symmetric with respect to the four-fold symmetry of the channel. We have now extended these calculations to rubidium and cesium ions, and find that rubidium solvation resembles that of potassium (and both are permeant ions), while cesium resembles sodium (and both are non-permeant), in terms of the geometry of up to eight hydrating, and four non-hydrating, water molecules. In each case, a maximum of 12 water molecules are relevant to the calculation. The placement of the water molecules in the two cases is essentially the same as found from the electron density in the X-ray structure of Zhou and MacKinnon. For Na⁺ and K⁺, we show that energy decreases from bulk to the cavity to the lowest position in the selectivity filter (accurate energy could not be calculated for the heavier ions). A separate calculation shows that fixing the Na⁺ ion at the position of the K⁺ minimum, followed by re-optimization produced a significantly modified system, not something that could be produced by thermal fluctuations. Moving the K⁺ into the Na⁺ position in the upper cavity led to a small increase in energy, ≈ 3 k_BT, but was accompanied by large shifts in the positions of hydrating waters, which would create a major kinetic barrier. Therefore, thermal fluctuations could not invalidate the conclusions of the main calculations.     

Effects of conducting and blocking ions on the structure and stability of the potassium channel KcsA

This article reports on the interaction of conducting (K(+)) and blocking (Na(+)) monovalent metal ions with detergent-solubilized and lipid-reconstituted forms of the K(+) channel KcsA. Monitoring of the protein intrinsic fluorescence reveals that the two ions bind competitively to KcsA with distinct affinities (dissociation constants for the KcsA.K(+) and KcsA.Na(+) complexes of approximately 8 and 190 mm, respectively) and induce different conformations of the ion-bound protein. The differences in binding affinity as well as the higher K(+) concentration bathing the intracellular mouth of the channel, through which the cations gain access to the protein binding sites, should favor that only KcsA.K(+) complexes are formed under physiological-like conditions. Nevertheless, despite such prediction, it was also found that concentrations of Na(+) well below its dissociation constant and even in the presence of higher K(+) concentrations, cause a remarkable decrease in the protein thermal stability and facilitate thermal dissociation into subunits of the tetrameric KcsA, as concluded from the temperature dependence of the protein infrared spectra and from gel electrophoresis, respectively. These latter observations cannot be explained based on the occupancy of the binding sites from above and suggest that there must be additional ion binding sites, whose occupancy could not be detected by fluorescence and in which the affinity for Na(+) must be higher or at least similar to that of K(+). Moreover, cation binding as reported by means of fluorescence does not suffice to explain the large differences in free energy of stabilization involved in the formation of the KcsA.Na(+) and KcsA.K(+) complexes, which for the most part should arise from synergistic effects of the ion-mediated intersubunit interactions.     

Potassium permeation through the KcsA channel: a density functional study

We present a theoretical study on structural and electronic aspects of K+ permeation through the binding sites of the KcsA channel's selectivity filter. Density functional calculations are carried out on models taken from selected snapshots of a molecular dynamics simulation recently reported [FEBS Lett. 477 (2000) 37]. During the translocation process from one binding site to the other, the coordination number of the permeating K+ ion turns out to decrease and K+ ion polarizes significantly its ligands, backbone carbonyl groups and a water molecule. K+-induced polarization increases significantly at the transition state (TS) between the two binding sites. These findings suggest that polarization effects play a significant role in the microscopic mechanisms regulating potassium permeation.     

KcsA closed and open: modelling and simulation studies

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height ~5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of ~3.3 Å C RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid (<5 ns) exit of all three K⁺ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.     

Lipids in the structure,folding, and function of the KcsA K+ channel

Lipid molecules surround an ion channel in its native environment of cellular membranes. The importance of the lipid bilayer and the role of lipid protein interactions in ion channel structure and function are not well understood. Here we demonstrate that the bacterial potassium channel KcsA binds a negatively charged lipid molecule. We have defined the potential binding site of the lipid molecule on KcsA by X-ray crystallographic analysis of a complex of KcsA with a monoclonal antibody Fab fragment. We also demonstrate that lipids are required for the in vitro refolding of the KcsA tetramer from the unfolded monomeric state. The correct refolding of the KcsA tetramer requires lipids, but it is not dependent on negatively charged lipids as refolding takes place in the absence of such lipids. We confirm that the presence of negatively charged lipids is required for ion conduction through the KcsA potassium channel, suggesting that the lipid bound to KcsA is important for ion channel function.     

Structural models of the KtrB, TrkH, and Trk1,2 symporters based on the structure of the KcsA K(+) channel.

Three-dimensional computer modeling is used to further investigate the hypothesis forwarded in the accompanying paper of an evolutionary relationship between four related families of K(+) sympoter proteins and the superfamily of K(+) channel proteins. Atomic-scale models are developed for the transmembrane regions of one member from each of the three more distinct symporter families, i.e., a TrkH protein from Escherichia coli, a KtrB protein from Aquifex aeolicus, and a Trk1,2 protein from Schizosaccharomyces pombe. The portions of the four consecutive M1-P-M2 motifs in the symporters that can be aligned with K(+) channel sequences are modeled directly from the recently determined crystal structure of the KcsA K(+) channel from Streptomyces lividans. The remaining portions are developed using our previously accumulated theoretical modeling criteria and principles. Concurrently, the use of these criteria and principles is further supported by the now verified predictions of our previous K(+) channel modeling efforts and the degree to which they are satisfied by the known structure of the KcsA protein. Thus the observed ability of the portions of the symporter models derived from the KcsA crystal structure to also satisfy the theoretical modeling criteria provides additional support for an evolutionary link with K(+) channel proteins. Efforts to further satisfy the criteria and principles suggest that the symporter proteins from fungi and plants (i.e., Trk1,2 and HKT1) form dimeric and/or tetrameric complexes in the membrane. Furthermore, analysis of the atomic-scale models in relation to the sequence conservation within and between the protein families suggests structural details for previously proposed mechanisms for the linked symport of K(+) with Na(+) and H(+). Suggestions are also given for experiments to test these structures and hypotheses.     

Water and ions in a high resolution structure of B-DNA.

A detailed picture of hydration and counterion location in the B-DNA duplex d(GCGAATTCG) is presented. Detailed data have been obtained by single crystal x-ray diffraction at atomic resolution (0.89 A) in the presence of Mg(2+). The latter is the highest resolution ever obtained for a B-DNA oligonucleotide. Minor groove hydration is compared with that found in the Na(+) and Ca(2+) crystal forms of the related dodecamer d(CGCGAATTCGCG). High resolution data (1.45 A) of the Ca(2+) form obtained in our laboratory are used for that purpose. The central GAATTC has a very stable hydration spine identical in all cases, independent of duplex length and crystallization conditions (counterions, space group). However, the organization of the water molecules (tertiary and quaternary layers) associated with the central spine vary in each case.     

The influence of a membrane environment on the structure and stability of a prokaryotic potassium channel, KcsA

The lack of a membrane environment in membrane protein crystals is considered one of the major limiting factors to fully imply X-ray structural data to explain functional properties of ion channels [Gulbis, J.M. and Doyle, D. (2004) Curr. Opin. Struct. Biol. 14, 440-446]. Here, we provide infrared spectroscopic evidence that the structure and stability of the potassium channel KcsA and its chymotryptic derivative 1-125 KcsA reconstituted into native-like membranes differ from those exhibited by these proteins in detergent solution, the latter taken as an approximation of the mixed detergent-protein crystal conditions.     

Efficacy of external tetraethylammonium block of the KcsA potassium channel: Molecular and Brownian dynamics studies

Blockade of the KcsA potassium channel by externally applied tetraethylammonium is investigated using molecular dynamics calculations and Brownian dynamics simulations. In KcsA, the aromatic rings of four tyrosine residues located just external to the selectivity filter create an attractive energy well or a binding cage for a tetraethylammonium molecule. We first investigate the effects of re-orienting the four tyrosine residues such that the centers of the aromatic rings face the tetraethylammonium molecule directly. Then, we systematically move the residues inward in both orientations so that the radius of the binding cage formed by them becomes smaller. For each configuration, we construct a one-dimensional free energy profile by bringing in a tetraethylammonium molecule from the external reservoir toward the selectivity filter. The free energy profile is then converted to a one-dimensional potential energy profile, taking the available space between the tyrosine residues and the tetraethylammonium molecule into account. Incorporating this potential energy profile into the Brownian dynamics algorithm, we determine the conductance properties of the channel under various conditions, construct the current-tetraethylammonium-concentration curve and compare it with the experimentally determined inhibitory constant k_i for externally applied tetraethylammonium. We show that the experimentally determined binding affinity for externally applied tetraethylammonium can be replicated when each of the four tyrosine residues is moved inward by about 0.7 Å, irrespective of orientation of their aromatic rings.     

Molecular dynamics study of the KcsA channel at 2.0-A resolution: stability and concerted motions within the pore

The stability of the KcsA channel accommodating more than one ion in the pore has been studied with molecular dynamics. We have used the very last X-ray structure of the KcsA channel at 2.0-A resolution determined by Zhou et al. [Nature 414 (2001) 43]. In this channel, six of the seven experimentally evidenced sites have been considered. We show that the protein remains very stable in the presence of four K+ ions (three in the selectivity filter and one in the cavity). The locations and the respective distances of the different K+ ions and water molecules (W), calculated within our KWKWKK sequence, also fits well with the experimental observations. The analysis of the K+ ions and water molecules displacements shows concerted file motions on the simulated time scale (approximately 1 ns), which could act as precursor to the diffusion of K+ ions inside the channel. A simple one-dimensional dynamical model is used to interpret the concerted motions of the ions and water molecules in the pore leading ultimately to ion transfer.     

Molecular dynamics study of the KcsA channel at 2.0-Å resolution: stability and concerted motions within the pore

The stability of the KcsA channel accommodating more than one ion in the pore has been studied with molecular dynamics. We have used the very last X-ray structure of the KcsA channel at 2.0-Å resolution determined by Zhou et al. [Nature 414 (2001) 43]. In this channel, six of the seven experimentally evidenced sites have been considered. We show that the protein remains very stable in the presence of four K⁺ ions (three in the selectivity filter and one in the cavity). The locations and the respective distances of the different K⁺ ions and water molecules (W), calculated within our KWKWKK sequence, also fits well with the experimental observations. The analysis of the K⁺ ions and water molecules displacements shows concerted file motions on the simulated time scale (≈1 ns), which could act as precursor to the diffusion of K⁺ ions inside the channel. A simple one-dimensional dynamical model is used to interpret the concerted motions of the ions and water molecules in the pore leading ultimately to ion transfer.     

Extending the structure of an ABC transporter to atomic resolution: modeling and simulation studies of MsbA

Molecular modeling and simulation approaches have been use to generate a complete model of the prokaryotic ABC transporter MsbA from Escherichia coli, starting from the low-resolution structure-based Calpha trace (PDB code 1JSQ). MsbA is of some biomedical interest as it is homologous to mammalian transporters such as P-glycoprotein and TAP. The quality of the MsbA model is assessed using a combination of molecular dynamics simulations and static structural analysis. These results suggest that the approach adopted for MsbA may be of general utility for generating all atom models from low-resolution crystal structures of membrane proteins. Molecular dynamics simulations of the MsbA model inserted in a fully solvated octane slab (a membrane mimetic environment) reveal that while the monomer is relatively stable, the dimer is unstable and undergoes significant conformational drift on a nanosecond time scale. This suggests that the MsbA crystal dimer may not correspond to the MsbA dimer in vivo. An alternative model of the dimer is discussed in the context of available experimental data.     

The structure of the aquaporin-1 water channel: a comparison between cryo-electron microscopy and X-ray crystallography

Three different medium-resolution structures of the human water channel aquaporin-1 (AQP1) have been solved by cryo-electron microscopy (cryo-EM) during the last two years. Recently, the structure of the strongly related bovine AQP1 was solved by X-ray crystallography at higher resolution, allowing a validation of the original medium-resolution structures, and providing a good indication for the strengths and limitations of state of the art cryo-EM methods. We present a detailed comparison between the different models, which shows that overall, the structures are highly similar, deviating less than 2.5 A from each other in the helical backbone regions. The two original cryo-EM structures, however, also show a number of significant deviations from the X-ray structure, both in the backbone positions of the transmembrane helices and in the location of the amino acid side-chains facing the pore. In contrast, the third cryo-EM structure that included information from the X-ray structure of the homologous bacterial glycerol facilitator GlpF and that was subsequently refined against cryo-EM AQP1 data, shows a root mean square deviation of 0.9A from the X-ray structure in the helical backbone regions.     

Model studies of chromatin structure based on X-ray diffraction data.

Model calculations are presented in order to interpret the X-ray diffraction diagrams given by chromatin gels. It is shown that by taking into account the hydration of chromatin subunits, the problem of calculating the interference function in concentrated gels is greatly simplified. In this way it is spossible to fully interpret the influence of concentration on the position and intensity of the various rings present in the X-ray diffraction patterns. The possibilities and limitations of models based on spherical symmetry are also discussed. It is concluded that each chromatin subunit most likely contains three turns of DNA in each 200 base pairs segment surrounding a central protein core. With the method presented here it is possible to test if other models of chromatin based on different kinds of evidence are compatible with the X-ray diffraction data.     

X-ray studies of water structure in 2 Zn insulin crystal

The water structure of rhombohedral 2 Zn insulin crystal which contains about 280 water molecules and 0.55-0.60 mol citrate molecules per dimer has been studied by X-ray crystallographic refinement with 1.1 A resolution data. Atomic parameters of 83 fully occupied and 258 partially occupied water molecules and 0.3 mol of citrate were obtained. Full matrix least-squares method with isotropic temperature factor was used for the refinement of partially occupied water molecules. The water molecules in this crystal exist in one of the three states: fully occupied water, partially occupied water and water continuum, and a schematic model of water structure in protein crystal was proposed. The flexibility of water molecules is described.     

Modelling the effect of a GHz electric field on the dynamics of K+ ions in KcsA potassium channel

The dynamics of potassium ions in a KcsA channel, located within a stochastically fluctuating medium, is modelled via the application of the molecular dynamics simulation method. We investigate the effect of presence and absence of an applied electric field, either constant or periodic, on the dynamics of the channel. It is found that the ions undergo a hopping motion when the channel is exposed to a constant electric field of strength 0.03 V/nm. Furthermore, an alternating electric field in the GHz range, normally present in the daily environment and encountered by most biological systems, is applied to the channel, showing that in this frequency range, the rigidity of the atomic bonds of the filter is increased, which in turn disturbs the ionic passage rate through the filter. Consequently, in this frequency range, the application of electric fields may affect the function of such channels.     

Local and global structure of the monomeric subunit of the potassium channel KcsA probed by NMR

KcsA is a homotetrameric 68-kDa membrane-associated potassium channel which selectively gates the flux of potassium ions across the membrane. The channel is known to undergo a pH-dependent open-to-closed transition. Here we describe an NMR study of the monomeric subunit of the channel (KcsA^M), solubilized in SDS micelles. Chemical shift, solvent exchange, backbone ¹⁵N relaxation and residual dipolar coupling (RDC) data show the TM1 helix to remain intact, but the TM2 helix contains a distinct kink, which is subject to concentration-independent but pH-dependent conformational exchange on a microsecond time scale. The kink region, centered at G99, was previously implicated in the gating of the tetrameric KcsA channel. An RDC-based model of KcsA^M at acidic pH orients TM1 and the two helical segments of the kinked TM2 in a configuration reminiscent of the open conformation of the channel. Thus, the transition between states appears to be an inherent capability of the monomer, with the tetrameric assembly exerting a modulatory effect upon the transition which gives the channel its physiological gating profile.     

X-ray and neutron scattering studies of the hydration structure of alkali ions in concentrated aqueous solutions

The presence of ions in water provides a rich and varied environment in which many natural processes occur with important consequences in biology, geology and chemistry. This article will focus on the structural properties of ions in water and it will be shown how the 'difference' methods of neutron diffraction with isotopic substitution (NDIS) and anomalous X-ray diffraction (AXD) can be used to obtain direct information regarding the radial pair distribution functions of many cations and anions in solution. This information can subsequently be used to calculate coordination numbers and to determine ion-water conformation in great detail. As well as enabling comparisons to be made amongst ions in particular groups in the periodic table, such information can also be contrasted with results provided by molecular dynamics (MD) simulation techniques. To illustrate the power of these 'difference' methods, reference will be made to the alkali group of ions, all of which have been successfully investigated by the above methods, with the exception of the radioactive element francium. Additional comments will be made on how NDIS measurements are currently being combined with MD simulations to determine the structure around complex ions and molecules, many of which are common in biological systems.     

Surface structure and its dynamic rearrangements of the KcsA potassium channel upon gating and tetrabutylammonium blocking

KcsA is the first potassium channel for which the molecular structure was revealed. However, the high resolution structural information is limited to the transmembrane domain, and the dynamic picture of the full KcsA channel remains unsolved. We have developed a new approach to investigate the surface structure of proteins, and we applied this method to investigate the full length of the KcsA channel. Single-cysteine substitution was introduced into 25 sites, and specific reaction of these mutated channels to a bare surface of a flat gold plate was evaluated by surface plasmon resonance measurements. The surface plasmon resonance signals revealed the highest exposure for the mutant of the C-terminal end. When the gate of the KcsA channel is kept closed at pH 7.5, the extent of exposure showed periodic patterns for the consecutive sites located in the cytoplasmic (CP) and N-terminal domain. This suggests that these stretches take the alpha-helical structure. When the channel was actively gated at pH 4.0, many sites in the CP domain became exposed. Compared with the rigid structure in pH 7.5, these results indicate that the CP domain became loosely packed upon active gating. The C-terminal end of the M2 helix is a moving part of the gate, and it is exposed to the outer surface slightly at pH 4.0. By adding a channel blocker, tetrabutylammonium, the gate is further exposed. This suggests that in the active gating tetrabutylammonium keeps the gate open rather than being trapped in the central cavity.