1,721,084 research outputs found
Atypical mechanism of conduction in potassium channels.
Potassium channels can conduct passively K+ ions with rates of up to approximately 10(8) ions per second at physiological conditions, and they are selective to these species by a factor of 10(4) over Na+ ions. Ion conduction has been proposed to involve transitions between 2 main states, with 2 or 3 K+ ions occupying the selectivity filter separated by an intervening water molecule. The largest free energy barrier of such a process was reported to be of the order of 2-3 kcal mol(-1). Here, we present an alternative mechanism for conduction of K+ in potassium channels where site vacancies are involved, and we propose that coexistence of several ion permeation mechanisms is energetically possible. Conduction can be described as a more anarchic phenomenon than previously characterized by the concerted translocations of K+-water-K+
Dynamics, energetics, and selectivity of the low-K+ KcsA channel structure.
Potassium channels are a diverse family of integral membrane proteins through which K(+) can pass selectively. There is ongoing debate about the nature of conformational changes associated with the opening/closing and conductive/nonconductive states of potassium channels. The channels partly exert their function by varying their conductance through a mechanism known as C-type inactivation. Shortly after the activation of K(+) channels, their selectivity filter stops conducting ions at a rate that depends on various stimuli. The molecular mechanism of C-type inactivation has not been fully understood yet. However, the X-ray structure of the KcsA channel obtained in the presence of low K(+) concentration is thought to be representative of a K(+) channel in the C-type inactivated state. Here, extensive, fully atomistic molecular dynamics and free-energy simulations of the low-K(+) KcsA structure in an explicit lipid bilayer are performed to evaluate the stability of this structure and the selectivity of its binding sites. We find that the low-K(+) KcsA structure is stable on the timescale of the molecular dynamics simulations performed, and that ions preferably remain in S1 and S4. In the absence of ions, the selectivity filter evolves toward an asymmetric architecture, as already observed in other computations of the high-K(+) structure of KcsA and KirBac. The low-K(+) KcsA structure is not permeable by Na(+), K(+), or Rb(+), and the selectivity of its binding sites is different from that of the high-K(+) structure
Gating at the selectivity filter of ion channels that conduct Na+ and K+ ions.
The NaK channel is a cation selective channel with similar permeability for K(+) and Na(+). The available crystallographic structure of wild-type (WT) NaK is usually associated with a conductive state of the channel. Here, potential of mean force for complete conduction events of Na(+) and K(+) ions through NaK show that: i), large energy barriers prevent the passage of ions through the WT NaK structure, ii), the barriers are correlated to the presence of a hydrogen bond between Asp-66 and Asn-68, and iii), the structure of NaK mutated to mimic cyclic nucleotide-gated channels conducts Na(+) and K(+). These results support the hypothesis that the filter of cation selective channels can adopt at least two different structures: a conductive one, represented by the x-ray structures of the NaK-CNG chimeras, and a closed one, represented by the x-ray structures of the WT NaK
On ionic conduction in potassium channels.
In ref. 1, we presented an alternative mechanism for conduc- tion of K+ in K+ channels where site vacancies are involved, and we proposed that coexistence of several ion permeation mechanisms is energetically possible. Specifically, we found that conduction can be described as a more anarchic phe- nomenon than previously habitually characterized by the con- certed translocations of K+-water-K+. This alternative pathway entails the possible presence of vacancies, with neither K+ nor water molecules in certain sites; sometimes, ions can even be found at adjacent binding sites. Moreover, we sug- gested that this mechanism is likely to be just one example among a plethora of alternative configurations and conduction pathways that ions and water may adopt during permeation, and that it can be viewed as a perturbation to the long- standing accepted mechanism involving neat, organized ion–water fluxe
Critical Assessment of Common Force Fields for Molecular Dynamics Simulations of Potassium Channels
For the last two decades, the KcsA K+ channel has served as a case study to understand how potassium channels operate at the atomic scale, and molecular dynamics simulations have contributed significantly to the current knowledge of the atomic mechanisms of conduction, inactivation, and gating in this family of membrane proteins. Currently, microsecond trajectories are becoming the new standard in the field, and consequently, it is critical to assess and compare the performance of the classical force fields ordinarily used in simulations of biological systems as well as to quantitatively assess the agreement with experimental data for trajectories of this order of magnitude. To that extent, we performed classical molecular dynamics simulations with CHARMM36 and AMBER-ff14sb force fields using atomic models based on the experimental structure of the KcsA channel in the open/conductive state, at conditions of ionic concentrations and membrane potentials resembling the ones adopted in experiments. In simulations using the CHARMM force field, the experimental open/conductive structure of the channel exhibited high conformational plasticity and fast collapse toward an occluded state. In contrast, in an identical set of simulations using the AMBER force field, no major deviations from the experimental structure were recorded. Force field development is a complex process in which many approximations are typically required and adopted. The results presented here provide additional motivation to further improve the existing models and, crucially, alert practitioners about limitations
Selectivity and permeation of alkali metal ions in K +-channels
Ion conduction in K +-channels is usually described in terms of concerted movements of K + progressing in a single file through a narrow pore. Permeation is driven by an incoming ion knocking on those ions already inside the protein. A fine-tuned balance between high-affinity binding and electrostatic repulsive forces between permeant ions is needed to achieve efficient conduction. While K +-channels are known to be highly selective for K + over Na +, some K + channels conduct Na + in the absence of K +. Other ions are known to permeate K +-channels with a more moderate preference and unusual conduction features. We describe an extensive computational study on ion conduction in K +-channels rendering free energy profiles for the translocation of three different alkali ions and some of their mixtures. The free energy maps for Rb + translocation show at atomic level why experimental Rb + conductance is slightly lower than that of K +. In contrast to K + or Rb +, external Na + block K + currents, and the sites where Na + transport is hindered are characterized. Translocation of K +/Na + mixtures is energetically unfavorable owing to the absence of equally spaced ion-binding sites for Na +, excluding Na + from a channel already loaded with K +. © 2011 Elsevier Ltd. All rights reserved
Insights into the sliding movement of the lac repressor nonspecifically bound to DNA.
The Lac repressor finds its DNA binding sequences with an association rate 2 orders of magnitude higher than what is expected for a random diffusive process. This experimental data stimulated numerous theoretical and experimental studies, which led to the facilitated diffusion model. In facilitated diffusion, the Lac repressor binds nonspecifically to DNA. This nonspecific binding is followed by an exploration of the DNA molecule in a reduced space. Single-molecule imaging confirmed that the Lac repressor may move along the DNA molecule; however, it is still under debate whether the LacI movement proceeds through sliding, with a continuous close contact between the protein and DNA, or through hopping between adjacent binding sites. We have investigated the one-dimensional sliding movement of the Lac repressor along nonspecific DNA by full-atomistic molecular dynamics simulations and free-energy calculations based on the umbrella sampling technique. The computed free-energy profile along a helical trajectory was periodic, with periodicity equal to the distance between successive nucleotides and an energy barrier between successive minima of 8.7 +/- 0.7 kcal/mol. The results from the molecular simulations were subsequently used in a Langevin dynamics framework to estimate the diffusion coefficient of the Lac repressor sliding along nonspecific DNA. The computed diffusion coefficient is close to the lower limit of the experimental range
Effects of the Protonation States of the EEEE Motif of a Bacterial Na+-Channel on Conduction and Pore Structure
Permeation of water through the KcsA K+ channel
Previous studies have reported that the KcsA potassium channel has an osmotic permeability coefficient of 4.8 x 10(-12) cm3/s, giving it a significantly higher osmotic permeability coefficient than that of some membrane channels specialized in water transport. This high osmotic permeability is proposed to occur when the channel is depleted of potassium ions, the presence of which slow down the water permeation process. The atomic structure of the potassium-depleted KcsA channel and the mechanisms of water permeation have not been well characterized so far. Here, all-atom molecular dynamics simulations, in conjunction with an umbrella sampling strategy and a nonequilibrium approach to simulate pressure gradients are employed to illustrate the permeation of water in the absence of ions through the KcsA K+ channel. Equilibrium molecular dynamics simulations (95 ns combined total length) identified a possible structure of the potassium-depleted KcsA channel, and umbrella sampling calculations (160 ns combined total length) revealed that this structure is not permeable by water molecules moving along the channel axis. The simulation of a pressure gradient across the channel (30 ns combined total length) identified an alternative permeation pathway with a computed osmotic permeability of approximately (2.7 +/- 0.9) x 10(-13) cm3/s. Water fluxes along this pathway did not proceed through collective water motions or transitions to vapor state. All of the major results of this study were robust against variations in a wide set of simulation parameters (force field, water model, membrane model, and channel conformation)
DNA-recognition process described by MD simulations of the lactose repressor protein on a specific and a non-specific DNA sequence.
The lactose repressor protein may bind DNA in two possible configurations: a specific one, if the DNA sequence corresponds to a binding site, and a non-specific one otherwise. To find its target sequences, the lactose repressor first binds non-specifically to DNA, and subsequently, it rapidly searches for a binding site. Atomic structures of non-specific and specific complexes are available from crystallographic and nuclear magnetic resonance experiments. However, what remains unknown is a detailed description of the steps that transform the non-specific complex into the specific one. Here, how the protein first recognizes its binding site has been studied using molecular dynamics simulations. The picture that emerges is that of a protein that is as mobile when interacting with non-specific DNA sequences as when free in solution. This high degree of mobility allows the protein to rapidly sample different DNA sequences. In contrast, when the protein encounters a binding site, the configuration ensemble collapses, and the protein sliding movements along the DNA sequence become scarce. The binding energies in the specific and non-specific complexes were analysed using the Molecular Mechanics Poisson Boltzmann Surface Area approach. These results represent a first step towards a throughout characterization of the DNA-recognition process
- …
