Molecular Dynamics Simulations of a K Channel Model: Sensitivity to Changes in Ions, Waters, and Membrane Environment

Calcium interactions with Cx26 hemmichannel: Spatial association between MD simulations biding sites and variant pathogenicity

Connexinophaties are a collective of diseases related to connexin channels and hemichannels. In particular many Cx26 alterations are strongly associated to human deafness. Calcium plays an important role on this structures regulation. Here, using calcium as a probe, extensive atomistic Molecular Dynamics simulations were performed on the Cx26 hemichannel embedded in a lipid bilayer. Exploring different initial conditions and calcium concentration, simulation reached ∼4 μs. Several analysis were carried out in order to reveal the calcium distribution and localization, such as electron density profiles, density maps and distance time evolution, which is directly associated to the interaction energy. Specific amino acid interactions with calcium and their stability were capture within this context. Few of these sites such as, GLU42, GLU47, GLY45 and ASP50, were already suggested in the literature. Besides, we identified novel calcium biding sites: ASP2, ASP117, ASP159, GLU114, GLU119, GLU120 and VAL226. To the best of our knowledge, this is the first time that these sites are reported within this context. Furthermore, since various pathologies involving the Cx26 hemichannel are associated with pathogenic variants in the corresponding CJB2 gene, using ClinVar, we were able to spatially associate the 3D positions of the identified calcium binding sites within the framework of this work with reported pathogenic variants in the CJB2 gene. This study presents a first step on finding associations between molecular features and pathological variants of the Cx26 hemichannel.

Changing water affinity from hydrophobic to hydrophilic in hydrophobic channels

The behavior of water at hydrophobic interfaces can play a significant role in determining chemical reaction outcomes and physical properties. Carbon nanotubes and aluminophosphate materials have one-dimensional hydrophobic channels, which are entirely surrounded by hydrophobic interfaces. Unique water behavior was observed in such hydrophobic channels. In this article, changes in the water affinity in one-dimensional hydrophobic channels were assessed using water vapor adsorption isotherms at 303 K and grand canonical Monte Carlo simulations. Hydrophobic behavior of water adsorbed in channels wider than 3 nm was observed for both adsorption and desorption processes, owing to the hydrophobic environment. However, water showed hydrophilic properties in both adsorption and desorption processes in channels narrower than 1 nm. In intermediate-sized channels, the hydrophobic properties of water during the adsorption process were seen to transition to hydrophilic behavior during the desorption process. Hydrophilic properties in the narrow channels for both adsorption and desorption processes are a result of the relatively strong water-channel interactions (10-15 kJ mol(-1)). In the 2-3 nm channels, the water-channel interaction energy of 4-5 kJ mol(-1) was comparable to the thermal translational energy. The cohesive water interaction was approximately 35 kJ mol(-1), which was larger than the others. Thus, the water affinity change in the 2-3 nm channels for the adsorption and desorption processes was attributed to weak water-channel interactions and strong cohesive interactions. These results are inherently important to control the properties of water in hydrophobic environments.

Construction and validation of a homology model of the human voltage-gated proton channel hHV1

The topological similarity of voltage-gated proton channels (H_V1s) to the voltage-sensing domain (VSD) of other voltage-gated ion channels raises the central question of whether H_V1s have a similar structure. We present the construction and validation of a homology model of the human H_V1 (hH_V1). Multiple structural alignment was used to construct structural models of the open (proton-conducting) state of hH_V1 by exploiting the homology of hH_V1 with VSDs of K⁺ and Na⁺ channels of known three-dimensional structure. The comparative assessment of structural stability of the homology models and their VSD templates was performed using massively repeated molecular dynamics simulations in which the proteins were allowed to relax from their initial conformation in an explicit membrane mimetic. The analysis of structural deviations from the initial conformation based on up to 125 repeats of 100-ns simulations for each system reveals structural features consistently retained in the homology models and leads to a consensus structural model for hH_V1 in which well-defined external and internal salt-bridge networks stabilize the open state. The structural and electrostatic properties of this open-state model are compatible with proton translocation and offer an explanation for the reversal of charge selectivity in neutral mutants of Asp¹¹². Furthermore, these structural properties are consistent with experimental accessibility data, providing a valuable basis for further structural and functional studies of hH_V1. Each Arg residue in the S4 helix of hH_V1 was replaced by His to test accessibility using Zn²⁺ as a probe. The two outermost Arg residues in S4 were accessible to external solution, whereas the innermost one was accessible only to the internal solution. Both modeling and experimental data indicate that in the open state, Arg²¹¹, the third Arg residue in the S4 helix in hH_V1, remains accessible to the internal solution and is located near the charge transfer center, Phe¹⁵⁰.

Determination of the charge profile in the KcsA selectivity filter using ab initio calculations and molecular dynamics simulations

The charge profile of K+ and Na+ ions moving in a single file across the filter of the KcsA channel is determined using both molecular dynamics simulations and ab initio calculations. We show a strong correlation between the charge variation and the ion location resulting in a saw-tooth profile, which provides additional information on the influence of charge transfer on the permeation and selectivity of the channel.

Anomalous hydration shell order of Na+ and K+ inside carbon nanotubes

We performed molecular dynamics simulations of the hydration of Na+ and K+ in infinitely long single-walled armchair carbon nanotubes (CNTs) at 298 K. Simulation results indicate that the preferential orientation of water molecules in coordination shells of these two cations presents an anomalous change in the CNTs and causes a diameter-dependent variation for the interaction energy between the cation and water molecules in its coordination shell. In the five CNTs of this work, it is energetically favorable for confining a hydrated K+ inside the two narrow CNTs with diameters of 0.60 and 0.73 nm, whereas the situation is reverse inside the wide CNTs with diameters of 0.87, 1.0, and 1.28 nm. This finding is important for CNT applications in ionic systems that control the selectivity and the ionic flow.

Molecular dynamics simulations of Na+/Cl(-)-dependent neurotransmitter transporters in a membrane-aqueous system

We have performed molecular dynamics simulations of a homology model of the human serotonin transporter (hSERT) in a membrane environment and in complex with either the natural substrate 5-HT or the selective serotonin reuptake inhibitor escitalopram. We have also included a transporter homologue, the Aquifex aeolicus leucine transporter (LeuT), in our study to evaluate the applicability of a simple and computationally attractive membrane system. Fluctuations in LeuT extracted from simulations are in good agreement with crystallographic B factors. Furthermore, key interactions identified in the X-ray structure of LeuT are maintained throughout the simulations indicating that our simple membrane system is suitable for studying the transmembrane protein hSERT in complex with 5-HT or escitalopram. For these transporter complexes, only relatively small fluctuations are observed in the ligand-binding cleft. Specific interactions responsible for ligand recognition, are identified in the hSERT-5HT and hSERT-escitalopram complexes. Our findings are in good agreement with predictions from mutagenesis studies.

Is the mobility of the pore walls and water molecules in the selectivity filter of KcsA channel functionally important?

We performed in-depth analysis of the forces which act on the K(+) ions in the selectivity filter of the KcsA channel in order to estimate the relative importance of static and dynamic influence of the filter wall and water molecules on ion permeation and selectivity. The forces were computed using the trajectories of all-atom molecular dynamics simulations. It is shown that the dynamics of the selectivity filter contributes about 3% to the net force acting on the ions and can be neglected in the studies focused on the macroscopic properties of the channel, such as the current. Among the filter atoms, only the pore-forming carbonyl groups can be considered as dynamic in the studies of microscopic events of conduction, while the dynamic effects from all other atoms are negligible. We also show that the dynamics of the water molecules in the filter can not be neglected. The fluctuating forces from the water molecules can be as strong as net forces from the pore walls and can effectively drive the ions through the local energy barriers in the filter.

Molecular dynamics simulations of asymmetric NaCl and KCl solutions separated by phosphatidylcholine bilayers: potential drops and structural changes induced by strong Na+-lipid interactions and finite size effects

Differences of ionic concentrations across lipid bilayers are some of the primary energetic driving forces for cellular electrophysiology. While macroscopic models of asymmetric ionic solutions are well-developed, their connection to ion, water, and lipid interactions at the atomic scale are much more poorly understood. In this study, we used molecular dynamics to examine a system of two chambers of equal ionic strength, but differing amounts of NaCl and KCl, separated by a lipid bilayer. Our expectation was that the net electrostatic potential difference between the two chambers should be small or zero. Contrary to our expectation, a large potential difference (-70 mV) slowly evolved across the two water chambers over the course of our 172-ns simulation. This potential primarily originated from strong Na(+) binding to the carbonyls of the phosphatidylcholine lipids. This ion adsorption also led to significant structural and mechanical changes in the lipid bilayer. We discuss this surprising result in the context of indirect experimental evidence for Na(+) interaction with bilayers as well as potential caveats in current biomembrane simulation methodology, including force-field parameters and finite size effects.

Molecular dynamics simulations of potassium channels

Despite the complexity of ion-channels, MD simulations based on realistic all-atom models have become a powerful technique for providing accurate descriptions of the structure and dynamics of these systems, complementing and reinforcing experimental work. Successful multidisciplinary collaborations, progress in the experimental determination of three-dimensional structures of membrane proteins together with new algorithms for molecular simulations and the increasing speed and availability of supercomputers, have made possible a considerable progress in this area of biophysics. This review aims at highlighting some of the work in the area of potassium channels and molecular dynamics simulations where numerous fundamental questions about the structure, function, folding and dynamics of these systems remain as yet unresolved challenges.

Molecular dynamics simulations of inwardly rectifying (Kir) potassium channels: a comparative study

Inward rectifier potassium (Kir) channels regulate cell excitability and transport K+ ions across membranes. Homotetrameric models of three mammalian Kir channels (Kir1.1, Kir3.1, and Kir6.2) have been generated, using the KirBac3.1 transmembrane and rat Kir3.1 intracellular domain structures as templates. All three models have been explored by 10 ns molecular dynamics simulations in phospholipid bilayers. Analysis of the initial structures revealed conservation of potential lipid interaction residues (Trp/Tyr and Arg/Lys side chains near the lipid headgroup-water interfaces). Examination of the intracellular domains revealed key structural differences between Kir1.1 and Kir6.2 which may explain the difference in channel inhibition by ATP. The behavior of all three models in the MD simulations revealed that they have conformational stability similar to that seen for comparable simulations of, for example, structures derived from cryoelectron microscopy data. Local distortions of the selectivity filter were seen during the simulations, as observed in previous simulations of KirBac and in simulations and structures of KcsA. These may be related to filter gating of the channel. The intracellular hydrophobic gate does not undergo any substantial changes during the simulations and thus remains functionally closed. Analysis of lipid-protein interactions of the Kir models emphasizes the key role of the M0 (or "slide") helix which lies approximately parallel to the bilayer-water interface and forms a link between the transmembrane and intracellular domains of the channel.

Conformational sampling and dynamics of membrane proteins from 10-nanosecond computer simulations

In the current report, we provide a quantitative analysis of the convergence of the sampling of conformational space accomplished in molecular dynamics simulations of membrane proteins of duration in the order of 10 nanoseconds. A set of proteins of diverse size and topology is considered, ranging from helical pores such as gramicidin and small beta-barrels such as OmpT, to larger and more complex structures such as rhodopsin and FepA. Principal component analysis of the C(alpha)-atom trajectories was employed to assess the convergence of the conformational sampling in both the transmembrane domains and the whole proteins, while the time-dependence of the average structure was analyzed to obtain single-domain information. The membrane-embedded regions, particularly those of small or structurally simple proteins, were found to achieve reasonable convergence. By contrast, extra-membranous domains lacking secondary structure are often markedly under-sampled, exhibiting a continuous structural drift. This drift results in a significant imprecision in the calculated B-factors, which detracts from any quantitative comparison to experimental data. In view of such limitations, we suggest that similar analyses may be valuable in simulation studies of membrane protein dynamics, in order to attach a level of confidence to any biologically relevant observations.

Conformational dynamics of M2 helices in KirBac channels: helix flexibility in relation to gating via molecular dynamics simulations

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.

Membrane protein structure quality in molecular dynamics simulation

Our goal was to assess the relationship between membrane protein quality, output from protein quality checkers and output from molecular dynamics (MD) simulations. Membrane transport proteins are essential for a wide range of cellular processes. Structural features of integral membrane proteins are still under-explored due to experimental limitations in structure determination. Computational techniques can be used to exploit biochemical and medium resolution structural data, as well as sequence homology to known structures, and enable us to explore the structure-function relationships in several transmembrane proteins. The quality of the models produced is vitally important to obtain reliable predictions. An examination of the relationship between model stability in molecular dynamics (MD) simulations derived from RMSD (root mean squared deviation) and structure quality assessment from various protein quality checkers was undertaken. The results were compared to membrane protein structures, solved at various resolution, by either X-ray or electron diffraction techniques. The checking programs could predict the potential success of MD in making functional conclusions. MD stability was shown to be a good indicator for the quality of structures. The quality was also shown to be dependent on the resolution at which the structures were determined.

Computational analysis of mutations in the transmembrane region of Vpu from HIV-1

Vpu is an 81 amino acid integral membrane protein encoded by HIV-1. Its alpha-helical transmembrane (TM) domain (residues approximately 6-28) enhances virion release by oligomerizing into bundles and forming ion-conducting channels across the plasma membrane. Its cytoplasmic domain (residues approximately 29-81) is also alpha-helical and binds to the transmembrane protein CD4, inducing its degradation. Mutations within the TM domain have been found to abrogate enhanced particle release from the infected cell (Tiganos et al. Virology (1998) 251 96-107). A series of computational models of monomeric, pentameric and hexameric Vpu(1-31) mutants have been constructed, embedded in fully hydrated lipid bilayers and subjected to a 3 ns molecular dynamics (MD) simulation. None of the mutations has any destabilizing effect on the secondary and tertiary structure. One of the mutants, in which the position of a tryptophan residue within the TM domain is altered, is known not to induce CD4 degradation; an extended kinked model of this mutant has been generated (Vpu(1-52)IVW-k) and during subsequent MD simulations, the bend between the TM and a part of the cytoplasmic domain is found to unwind and a complex salt bridge involving Lys-37 is formed.

A strategy for integrative computational physiology

Organ function (the heart beat for example) can only be understood through knowledge of molecular and cellular processes within the constraints of structure-function relations at the tissue level. A quantitative modeling framework that can deal with these multiscale issues is described here under the banner of the International Union of Physiological Sciences Physiome Project.

Molecular dynamics simulation approaches to K channels: conformational flexibility and physiological function

Molecular modeling and simulations enable extrapolation for the structure of bacterial potassium channels to the function of their mammalian homologues. Molecular dynamics simulations have revealed the concerted single-file motion of potassium ions and water molecules through the selectivity filter of K channels and the role of filter flexibility in ion permeation and in "fast gating." Principal components analysis of extended K channel simulations suggests that hinge-bending of pore-lining M2 (or S6) helices plays a key role in K channel gating. Based on these and other simulations, a molecular model for gating of inward rectifier K channel gating is presented.

Targeted molecular dynamics of an open-state KcsA channel

Pore opening of KcsA channel is studied using targeted molecular dynamics simulations. Conformational changes of the protein are determined, starting from the crystallized refined 2.0 A structure (pdb 1K4C) determined in x-ray experiments and arriving to the open-state structure constructed on the basis of electron paramagnetic resonance spectroscopic data (pdb 1JQ1). Our results corroborate the essential role played by the terminal residues located on the transmembrane helices M2 which were not taken into account at that time. The aperture mechanism of the channel appears to be ziplike. A small constraint (approximately equal to 5 x 10(-2) kcal mol(-1) A(-2) per C(alpha)) applied to the terminal residues located on the intracellular side is sufficient to initialize the pore opening at the innermost part of the gate, but additional constraint must be applied to definitely complete the pore aperture. The open structure is proved to be a metastable state since releasing the constraint leads to another relaxed open conformation which seems to reach stability.

Filter flexibility and distortion in a bacterial inward rectifier K+ channel: simulation studies of KirBac1.1

The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (<0.1 nm) filter flexibility during the concerted motion of ions and water molecules within the filter, such local changes in conformation occurring on an approximately 1-ns timescale. In the absence of K(+) ions, the KirBac selectivity filter undergoes more substantial distortions. These resemble those seen in comparable simulations of other channels (e.g., KcsA and KcsA-based homology models) and are likely to lead to functional closure of the channel. This suggests filter distortions may provide a mechanism of K-channel gating in addition to changes in the hydrophobic gate formed at the intracellular crossing point of the M2 helices. The simulation data also provide evidence for interactions of the "slide" (pre-M1) helix of KirBac with phospholipid headgroups.

Role of water molecules in the KcsA protein channel by molecular dynamics calculations

Molecular dynamics simulations supported by electrostatic calculations have been conducted on the KcsA channel to determine the role of water molecules in the pore. Starting from the X-ray structure of the KcsA channel in its closed state at 2.0 angstroms resolution, the opening of the pore towards a conformation built on the basis of EPR results is studied. We show that water molecules act as a structural element for the K+ ions inside the filter and the hydrophobic cavity of the channel. In the filter, water tends to enhance the depth of the wells occupied by the K+ ions, while in the cavity there is a strong correlation between the water molecules and the cavity ion. As a consequence, the protein remains very stable in the presence of three K+ ions in the selectivity filter and one in the cavity. The analysis of the dynamics of water molecules in the cavity reveals preferred orientations of the dipoles along the pore axis, and a correlated behavior between this dipole orientation and the displacement of the K+ ion during the gating process.

Prediction of the mutation-induced change in thermodynamic stabilities of membrane proteins from free energy simulations

Comparative protein structure modeling and free energy perturbation simulation have been applied in a consecutive manner to investigate the mutation-induced stabilization of membrane proteins (MPs) in aqueous solution without knowledge of their three-dimensional structures. The calculated difference in protein solvation free energy between the wild type and a mutant compares well with their relative thermodynamic stabilities in solution. For monomeric MPs, a mutant reveals a higher stability than the wild type if the calculated solvation free energy indicates a favorable change. On the contrary, for oligomeric MPs the stability of a mutant increases as the solvation free energy of a mutated monomer becomes less favorable, indicating that the oligomeric MP mutant would be stabilized in solution due to the reduced desolvation cost for oligomerization. The present computational strategy is expected to find its way as a useful tool for assessing the relative stability of a mutant MP with respect to its wild type in solution.

Toward the active conformations of rhodopsin and the beta2-adrenergic receptor

Using sets of experimental distance restraints, which characterize active or inactive receptor conformations, and the X-ray crystal structure of the inactive form of bovine rhodopsin as a starting point, we have constructed models of both the active and inactive forms of rhodopsin and the beta2-adrenergic G-protein coupled receptors (GPCRs). The distance restraints were obtained from published data for site-directed crosslinking, engineered zinc binding, site-directed spin-labeling, IR spectroscopy, and cysteine accessibility studies conducted on class A GPCRs. Molecular dynamics simulations in the presence of either "active" or "inactive" restraints were used to generate two distinguishable receptor models. The process for generating the inactive and active models was validated by the hit rates, yields, and enrichment factors determined for the selection of antagonists in the inactive model and for the selection of agonists in the active model from a set of nonadrenergic GPCR drug-like ligands in a virtual screen using ligand docking software. The simulation results provide new insights into the relationships observed between selected biochemical data, the crystal structure of rhodopsin, and the structural rearrangements that occur during activation.

Coupled motions between pore and voltage-sensor domains: a model for Shaker B, a voltage-gated potassium channel

A high-resolution crystal structure of KvAP, an archeabacterial voltage-gated potassium (Kv) channel, complexed with a monoclonal Fab fragment has been recently determined. Based on this structure, a mechanism for the activation (opening) of Kv channels has been put forward. This mechanism has since been criticized, suggesting that the resolved structure is not representative of the family of voltage-gated potassium channels. Here, we propose a model of the transmembrane domain of Shaker B, a well-characterized Kv channel, built by homology modeling and docking calculations. In this model, the positively charged S4 helices are oriented perpendicular to the membrane and localized in the groove between segments S5 and S6 of adjacent subunits. The structure and the dynamics of the full atomistic model embedded in a hydrated lipid bilayer were investigated by means of two large-scale molecular dynamics simulations under transmembrane-voltage conditions known to induce, respectively, the resting state (closed) and the activation (opening) of voltage-gated channels. Upon activation, the model undergoes conformational changes that lead to an increase of the hydration of the charged S4 helices, correlated with an upward translation and a tilting of the latter, concurrently with movements of the S5 helices and the activation gate. Although small, these conformational changes ultimately result in an alteration of the ion-conduction pathway. Our findings support the transporter model devised by Bezanilla and collaborators, and further underline the crucial role played by internal hydration in the activation of the channel.

Molecular dynamics simulations of the lipid bilayer edge

Phospholipid bilayers have been intensively studied by molecular dynamics (MD) simulation in recent years. The properties of bilayer edges are important in determining the structure and stability of pores formed in vesicles and biomembranes. In this work, we use molecular dynamics simulation to investigate the structure, dynamics, and line tension of the edges of bilayer ribbons composed of pure dimyristoylphosphatidylcholine (DMPC) or palmitoyl-oleoylphosphatidylethanolamine (POPE). As expected, we observe a significant reorganization of lipids at and near the edges. The treatment of electrostatic effects is shown to have a qualitative impact on the structure and stability of the edge, and significant differences are observed in the dynamics and structure of edges formed by DMPC and palmitoyl-oleoylphosphatidylethanolamine. From the pressure anisotropy in the simulation box, we calculate a line tension of approximately 10-30 pN for the DMPC edge, in qualitative agreement with experimental estimates for similar lipids.

Homology modelling and molecular dynamics simulations: comparative studies of human aquaporin-1

The structures of the mammalian water transport protein Aqp1 and of its bacterial homologue GlpF enables us to test whether homology models can be used to explore relationships between structure, dynamics and function in mammalian transport proteins. Molecular dynamics simulations (totalling almost 40 ns) were performed starting from: the X-ray structure of Aqp1; a homology model of Aqp1 based on the GlpF structure; and intermediate resolution structures of Aqp1 derived from electron microscopy. Comparisons of protein RMSDs vs. time suggest that the homology models are of comparable conformational stability to the X-ray structure, whereas the intermediate resolution structures exhibit significant conformation drift. For simulations based on the X-ray structure and on homology models, the flexibility profile vs. residue number correlates well with the crystallographic B-values for each residue. In the simulations based on intermediate resolution structures, mobility of the highly conserved NPA loops is substantially higher than in the simulations based on the X-ray structure or the homology models. Pore radius profiles remained relatively constant in the X-ray and homology model simulations but showed substantial fluctuations (reflecting the higher NPA loop mobility) in the intermediate resolution simulations. The orientation of the dipoles of water molecules within the pore is of key importance in maintaining low proton permeability through Aqp1. This property seems to be quite robust to the starting model used in the simulation. These simulations suggest that homology models based on bacterial homologues may be used to derive functionally relevant information on the structural dynamics of mammalian transport proteins.

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.

Modeling the concentration-dependent permeation modes of the KcsA potassium ion channel

The potassium channel from Streptomyces lividans (KcsA) is an integral membrane protein with sequence similarity to all known potassium channels, particularly in the selectivity filter region. A recently proposed model for ion channels containing either n or (n-1) single-file ions in their selectivity filters [P. H. Nelson, J. Chem. Phys. 177, 11396 (2002)] is applied to published KcsA channel K+ permeation data that exhibit a high-affinity process at low concentrations and a low-affinity process at high concentrations [M. LeMasurier et al., J. Gen. Physiol. 118, 303 (2001)]. The kinetic model is shown to provide a reasonable first-order explanation for both the high- and low-concentration permeation modes observed experimentally. The low-concentration mode ([K+]<200 mM) has a 200-mV dissociation constant of 56 mM and a conductance of 88 pS. The high-concentration mode ([K+]>200 mM) has a 200-mV dissociation constant of 1100 mM and a conductance of 500 pS. Based on the permeation model, and x-ray analysis [J. H. Morais-Cabral et al., Nature (London) 414, 37 (2001)], it is suggested that the experimentally observed K+ permeation modes correspond to an n=3 mechanism at high concentrations and an n=2 mechanism at low concentrations. The ratio of the electrical dissociation distances for the high- and low-concentration modes is 3:2, also consistent with the proposed n=3 and n=2 modes. Model predictions for K+ channels that exhibit asymmetric current-voltage (I-V) curves are presented, and further validation of the kinetic model via molecular simulation and experiment is discussed. The qualitatively distinct I-V characteristics exhibited experimentally by Tl+, NH+4, and Rb+ ions at 100 mM concentration can also be explained using the model, but more extensive experimental tests are required for quantitative validation of the model predictions.

Potassium channel, ions, and water: simulation studies based on the high resolution X-ray structure of KcsA

Interactions of Na(+), K(+), Rb(+), and Cs(+) ions within the selectivity filter of a potassium channel have been investigated via multiple molecular dynamics simulations (total simulation time, 48 ns) based on the high resolution structure of KcsA, embedded in a phospholipid bilayer. As in simulations based on a lower resolution structure of KcsA, concerted motions of ions and water within the filter are seen. Despite the use of a higher resolution structure and the inclusion of four buried water molecules thought to stabilize the filter, this region exhibits a significant degree of flexibility. In particular, pronounced distortion of filter occurs if no ions are present within it. The two most readily permeant ions, K(+) and Rb(+), are similar in their interactions with the selectivity filter. In contrast, Na(+) ions tend to distort the filter by binding to a ring of four carbonyl oxygens. The larger Cs(+) ions result in a small degree of expansion of the filter relative to the x-ray structure. Cs(+) ions also appear to interact differently with the gate region of the channel, showing some tendency to bind within a predominantly hydrophobic pocket. The four water molecules buried between the back of the selectivity filter and the remainder of the protein show comparable mobility to the surrounding protein and do not exchange with water molecules within the filter or the central cavity. A preliminary comparison of the use of particle mesh Ewald versus cutoff protocols for the treatment of long-range electrostatics suggests some difference in the kinetics of ion translocation within the filter.

A prokaryotic glutamate receptor: homology modelling and molecular dynamics simulations of GluR0

GluR0 is a prokaryotic homologue of mammalian glutamate receptors that forms glutamate-activated, potassium-selective ion channels. The topology of its transmembrane (TM) domain is similar to that of simple potassium channels such as KcsA. Two plausible alignments of the sequence of the TM domain of GluR0 with KcsA are possible, differing in the region of the P helix. We have constructed homology models based on both alignments and evaluated them using 6 ns duration molecular dynamics simulations in a membrane-mimetic environment. One model, in which an insertion in GluR0 relative to KcsA is located in the loop between the M1 and P helices, is preferred on the basis of lower structural drift and maintenance of the P helix conformation during simulation. This model also exhibits inter-subunit salt bridges that help to stabilise the TM domain tetramer. During the simulation, concerted K(+) ion-water movement along the selectivity filter is observed, as is the case in simulations of KcsA. K(+) ion exit from the central cavity is associated with opening of the hydrophobic gate formed by the C-termini of the M2 helices. In the intact receptor the opening of this gate will be controlled by interactions with the extramembranous ligand-binding domains.

Effects of confinement on static and dynamical properties of water

Molecular-dynamics results on water confined in a silica pore are reviewed and discussed in connection with experiments performed on water in Vycor and with studies of water in contact with proteins. The properties of confined water are studied as a function of both temperature and hydration level. The interaction of water in the film close to the substrate with the silica atoms induces a strong distortion of the hydrogen bond network. At high hydration levels a double dynamical regime is observed. At low hydration an anomalous diffusion is found upon supercooling with a transition from a Brownian to a non-Brownian regime on approaching the substrate in agreement with results found in studies of water in contact with globular proteins.

Investigating lipid composition effects on the mechanosensitive channel of large conductance (MscL) using molecular dynamics simulations

Previous experimental work has shown that the functional properties of the mechanosensitive channel of large conductance (MscL) are affected by variations in lipid composition. Here, we utilize molecular dynamics simulations of Mycobacterium tuberculosis MscL to investigate such lipid composition effects on a molecular level. In particular, two sets of simulations were performed. In the first, trajectories using lipids with different headgroups (phosphatidylcholine and phosphatidylethanolamine) were compared. Protein-lipid interactions were clearly altered by the headgroup changes, leading to conformational differences in the C-terminal region of M. tuberculosis MscL. In the second set of simulations, lipid tails were gradually shortened, thinning the membrane over a molecular dynamics trajectory. These simulations showed evidence of hydrophobic matching between MscL and the lipid membrane, as previously proposed. For all simulations, protein-lipid interaction energies in the second transmembrane region were correlated to mutagenic data, emphasizing the importance of lipid interactions for proper MscL function.

Membrane protein dynamics versus environment: simulations of OmpA in a micelle and in a bilayer

The bacterial outer membrane protein OmpA is one of the few membrane proteins whose structure has been solved both by X-ray crystallography and by NMR. Crystals were obtained in the presence of detergent, and the NMR structure is of the protein in a detergent micelle. We have used 10 ns duration molecular dynamics simulations to compare the behaviour of OmpA in a detergent micelle and in a phospholipid bilayer. The dynamic fluctuations of the protein structure seem to be ca 1.5 times greater in the micelle environment than in the lipid bilayer. There are subtle differences between the nature of OmpA-detergent and OmpA-lipid interactions. As a consequence of the enhanced flexibility of the OmpA protein in the micellar environment, side-chain torsion angle changes are such as to lead to formation of a continuous pore through the centre of the OmpA molecule. This may explain the experimentally observed channel formation by OmpA.

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.

Filter flexibility in a mammalian K channel: models and simulations of Kir6.2 mutants

The single-channel conductance varies significantly between different members of the inward rectifier (Kir) family of potassium channels. Mutations at three sites in Kir6.2 have been shown to produce channels with reduced single-channel conductance, the largest reduction (to 40% of wild-type) being for V127T. We have used homology modeling (based on a KcsA template) combined with molecular dynamics simulations in a phosphatidycholine bilayer to explore whether changes in structural dynamics of the filter were induced by three such mutations: V127T, M137C, and G135F. Overall, 12 simulations of Kir6.2 models, corresponding to a total simulation time of 27 ns, have been performed. In these simulations we focused on distortions of the selectivity filter, and on the presence/absence of water molecules lying behind the filter, which form interactions with the filter and the remainder of the protein. Relative to the wild-type simulation, the V127T mutant showed significant distortion of the filter such that approximately 50% of the simulation time was spent in a closed conformation. While in this conformation, translocation of K(+) ions between sites S1 and S2 was blocked. The distorted filter conformation resembles that of the bacterial channel KcsA when crystallized in the presence of a low [K(+)]. This suggests filter distortion may be a possible general model for determining the conductance of K channels.

Structures of neat and hydrated 1-octanol from computer simulations

Molecular dynamics computer simulations of 1-octanol and its mixtures with water have been performed. The liquid is composed of regions enriched in either hydrocarbons or hydroxyl groups. In neat octanol, the hydroxyl groups form clusters of long, thin chains. Upon the addition of water, the clusters become longer and more spherical, forming a structure that can be described as consisting of "overlapping elongated inverse micelles". The structures of the mixtures obtained at different hydration levels are consistent with those of experimental diffraction studies of water/octanol mixtures and previous computer simulations of neat and water-saturated octanol. The saturation point of the model has been calculated using the cavity-bias particle insertion method. The solubility of water in octanol is slightly too low compared to experimental results, and suggestions for possible improvements to the force field are made.

Analysis and evaluation of channel models: simulations of alamethicin

Alamethicin is an antimicrobial peptide that forms stable channels with well-defined conductance levels. We have used extended molecular dynamics simulations of alamethicin bundles consisting of 4, 5, 6, 7, and 8 helices in a palmitoyl-oleolyl-phosphatidylcholine bilayer to evaluate and analyze channel models and to link the models to the experimentally measured conductance levels. Our results suggest that four helices do not form a stable water-filled channel and might not even form a stable intermediate. The lowest measurable conductance level is likely to correspond to the pentamer. At higher aggregation numbers the bundles become less symmetrical. Water properties inside the different-sized bundles are similar. The hexamer is the most stable model with a stability comparable with simulations based on crystal structures. The simulation was extended from 4 to 20 ns or several times the mean passage time of an ion. Essential dynamics analyses were used to test the hypothesis that correlated motions of the helical bundles account for high-frequency noise observed in open channel measurements. In a 20-ns simulation of a hexameric alamethicin bundle, the main motions are those of individual helices, not of the bundle as a whole. A detailed comparison of simulations using different methods to treat long-range electrostatic interactions (a twin range cutoff, Particle Mesh Ewald, and a twin range cutoff combined with a reaction field correction) shows that water orientation inside the alamethicin channels is sensitive to the algorithms used. In all cases, water ordering due to the protein structure is strong, although the exact profile changes somewhat. Adding an extra 4-nm layer of water only changes the water ordering slightly in the case of particle mesh Ewald, suggesting that periodicity artifacts for this system are not serious.

Potassium channels: structures, models, simulations

Potassium channels have been studied intensively in terms of the relationship between molecular structure and physiological function. They provide an opportunity to integrate structural and computational studies in order to arrive at an atomic resolution description of mechanism. We review recent progress in K channel structural studies, focussing on the bacterial channel KcsA. Structural studies can be extended via use of computational (i.e. molecular simulation) approaches in order to provide a perspective on aspects of channel function such as permeation, selectivity, block and gating. Results from molecular dynamics simulations are shown to be in good agreement with recent structural studies of KcsA in terms of the interactions of K(+) ions with binding sites within the selectivity filter of the channel, and in revealing the importance of filter flexibility in channel function. We discuss how the KcsA structure may be used as a template for developing structural models of other families of K channels. Progress in this area is explored via two examples: inward rectifier (Kir) and voltage-gated (Kv) potassium channels. A brief account of structural studies of ancillary domains and subunits of K channels is provided.

K(+) versus Na(+) ions in a K channel selectivity filter: a simulation study

Molecular dynamics simulations of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal significant differences in interactions of the selectivity filter with K(+) compared with Na(+) ions. K(+) ions and water molecules within the filter undergo concerted single-file motion in which they translocate between adjacent sites within the filter on a nanosecond timescale. In contrast, Na(+) ions remain bound to sites within the filter and do not exhibit translocation on a nanosecond timescale. Furthermore, entry of a K(+) ion into the filter from the extracellular mouth is observed, whereas this does not occur for a Na(+) ion. Whereas K(+) ions prefer to sit within a cage of eight oxygen atoms of the filter, Na(+) ions prefer to interact with a ring of four oxygen atoms plus two water molecules. These differences in interactions in the selectivity filter may contribute to the selectivity of KcsA for K(+) ions (in addition to the differences in dehydration energy between K(+) and Na(+)) and the block of KcsA by internal Na(+) ions. In our simulations the selectivity filter exhibits significant flexibility in response to changes in ion/protein interactions, with a somewhat greater distortion induced by Na(+) than by K(+) ions.