The protonation state of the Glu-71/Asp-80 residues in the KcsA potassium channel: a first-principles QM/MM molecular dynamics study

Fluoride Transport and Inhibition Across CLC Transporters

The Chloride Channel (CLC) family includes proton-coupled chloride and fluoride transporters. Despite their similar protein architecture, the former exchange two chloride ions for each proton and are inhibited by fluoride, whereas the latter efficiently transport one fluoride in exchange for one proton. The combination of structural, mutagenesis, and functional experiments with molecular simulations has pinpointed several amino acid changes in the permeation pathway that capitalize on the different chemical properties of chloride and fluoride to fine-tune protein function. Here we summarize recent findings on fluoride inhibition and transport in the two prototypical members of the CLC family, the chloride/proton transporter from Escherichia coli (CLC-ec1) and the fluoride/proton transporter from Enterococcus casseliflavus (CLCF-eca).

Unraveling of a Strongly Correlated Dynamical Network of Residues Controlling the Permeation of Potassium in KcsA Ion Channel

The complicated patterns of the single-channel currents in potassium ion channel KcsA are governed by the structural variability of the selectivity filter. A comparative analysis of the dynamics of the wild type KcsA channel and several of its mutants showing different conducting patterns was performed. A strongly correlated dynamical network of interacting residues is found to play a key role in regulating the state of the wild type channel. The network is centered on the aspartate D80 which plays the role of a hub by strong interacting via hydrogen bonds with residues E71, R64, R89, and W67. Residue D80 also affects the selectivity filter via its backbones. This network further compromises ions and water molecules located inside the channel that results in the mutual influence: the permeation depends on the configuration of residues in the network, and the dynamics of network's residues depends on locations of ions and water molecules inside the selectivity filter. Some features of the network provide a further understanding of experimental results describing the KcsA activity. In particular, the necessity of anionic lipids to be present for functioning the channel is explained by the interaction between the lipids and the arginine residues R64 and R89 that prevents destabilizing the structure of the selectivity filter.

The influence of membrane bilayer thickness on KcsA channel activity

Atomic resolution structures have provided significant insight into the gating and permeation mechanisms of various ion channels, including potassium channels. However, ion channels may also be regulated by numerous factors, including the physiochemical properties of the membrane in which they are embedded. For example, the matching of the bilayer's hydrophobic region to the hydrophobic external surface of the ion channel is thought to minimize the energetic penalty needed to solvate hydrophobic residues or exposed lipid tails. To understand the molecular basis of such regulation by hydrophobic matching requires examining channels in the presence of the lipid membrane. Here we examine the role of hydrophobic matching in regulating the activity of the model potassium channel, KcsA. ⁸⁶Rb⁺ influx assays and single-channel recordings indicate that the non-inactivating E71A KcsA channel is most active in thin bilayers (<diC18:1PC). Bilayer thickness affects the open probability of KcsA and not its unitary conductance. Molecular dynamics simulations indicate that the bilayer can sufficiently modify its dimensions to accommodate KcsA channels without major perturbations in the protein helical packing within the nanosecond timescale. Based on experimental results and MD simulations, we present a model in which bilayer thickness influences the stability of the open and closed conformations of the intracellular gate of KcsA, with minimal impact on the stability of the selectivity filter of the non-inactivating mutant, E71A.

CpHMD-Then-QM/MM Identification of the Amino Acids Responsible for the Anabaena Sensory Rhodopsin pH-Dependent Electronic Absorption Spectrum

Anabaena Sensory Rhodopsin (ASR), a microbial photoactive protein featuring the retinal chromophore in two different conformations, exhibits a pH-dependent electronic absorption spectrum. Using the recently developed CpHMD-then-QM/MM multiscale protocol applied to ASR embedded in a membrane model, the pH-induced changes in its maximum absorption wavelength have been reproduced and analyzed. While the acidic tiny red-shift is essentially correlated with the deprotonation of an aspartic acid located on the ASR extracellular side, the larger blue-shift experimentally reported at pH values larger than 5 involves a cluster of titrating residues sitting on the cytoplasmic side. The ASR pH-dependent spectrum is the consequence of the competitive stabilization of retinal ground and excited states by the protein electrostatic potential.

Chloride Ion Transport by the E. coli CLC Cl-/H+ Antiporter: A Combined Quantum-Mechanical and Molecular-Mechanical Study

We performed steered molecular dynamics (SMD) and umbrella sampling simulations of Cl^- ion migration through the transmembrane domain of a prototypical E. coli CLC Cl^-/H⁺ antiporter by employing combined quantum-mechanical (QM) and molecular-mechanical (MM) calculations. The SMD simulations revealed interesting conformational changes of the protein. While no large-amplitude motions of the protein were observed during pore opening, the side chain rotation of the protonated external gating residue Glu148 was found to be critical for full access of the channel entrance by Cl^-. Moving the anion into the external binding site (S_ext) induced small-amplitude shifting of the protein backbone at the N-terminal end of helix F. As Cl^- traveled through the pore, rigid-body swinging motions of helix R separated it from helix D. Helix R returned to its original position once Cl^- exited the channel. Population analysis based on polarized wavefunction from QM/MM calculations discovered significant (up to 20%) charge loss for Cl^- along the ion translocation pathway inside the pore. The delocalized charge was redistributed onto the pore residues, especially the functional groups containing π bonds (e.g., the Tyr445 side chain), while the charges of the H atoms coordinating Cl^- changed almost negligibly. Potentials of mean force computed from umbrella sampling at the QM/MM and MM levels both displayed barriers at the same locations near the pore entrance and exit. However, the QM/MM PMF showed higher barriers (~10 kcal/mol) than the MM PMF (~2 kcal/mol). Binding energy calculations indicated that the interactions between Cl^- and certain pore residues were overestimated by the semi-empirical PM3 Hamiltonian and underestimated by the CHARMM36 force fields, both of which were employed in the umbrella sampling simulations. In particular, CHARMM36 underestimated binding interactions for the functional groups containing π bonds, missing the stabilizations of the Cl^- ion due to electron delocalization. The results suggested that it is important to explore these quantum effects for accurate descriptions of the Cl^- transport.

H. R, Thamleena A. H, Joseph A, Ben A, V. U. K. Roles of different amino-acid residues towards binding and selective transport of K+ through KcsA K+-ion channel

Each amino acid in the selectivity filter plays a distinct role towards binding and transport of K⁺ ion through KcsA.

Protonation state of E71 in KcsA and its role for channel collapse and inactivation

The prototypical prokaryotic potassium channel KcsA alters its pore depending on the ambient potassium; at high potassium, it exists in a conductive form, and at low potassium, it collapses into a nonconductive structure with reduced ion occupancy. We present solid-state NMR studies of KcsA in which we test the hypothesis that an important channel-inactivation process, known as C-type inactivation, proceeds via a state similar to this collapsed state. We test this using an inactivation-resistant mutant E71A, and show that E71A is unable to collapse its pore at both low potassium and low pH, suggesting that the collapsed state is structurally similar to the inactivated state. We also show that E71A has a disordered selectivity filter. Using site-specific K(+) titrations, we detect a local change at E71 that is coupled to channel collapse at low K(+). To gain more insight into this change, we site specifically measure the chemical shift tensors of the side-chain carboxyls of E71 and its hydrogen bond partner D80, and use the tensors to assign protonation states to E71 and D80 at high K(+) and neutral pH. Our measurements show that E71 is protonated at pH 7.5 and must have an unusually perturbed pK(a) (> 7.5) suggesting that the change at E71 is a structural rearrangement rather than a protonation event. The results offer new mechanistic insights into why the widely used mutant KcsA-E71A does not inactivate and establish the ambient K(+) level as a means to populate the inactivated state of KcsA in a controlled way.

Transmembrane communication: general principles and lessons from the structure and function of the M2 proton channel, K⁺ channels, and integrin receptors

Signal transduction across biological membranes is central to life. This process generally happens through communication between different domains and hierarchical coupling of information. Here, we review structural and thermodynamic principles behind transmembrane (TM) signal transduction and discuss common themes. Communication between signaling domains can be understood in terms of thermodynamic and kinetic principles, and complex signaling patterns can arise from simple wiring of thermodynamically coupled domains. We relate this to functions of several signal transduction systems: the M2 proton channel from influenza A virus, potassium channels, integrin receptors, and bacterial kinases. We also discuss key features in the structural rearrangements responsible for signal transduction in these systems.

Membrane protein dynamics from femtoseconds to seconds

Membrane proteins play a key role in energy conversion, transport, signal recognition, transduction, and other fundamental biological processes. Despite considerable progress in experimental techniques, the determination of structure and dynamics of membrane proteins still represents a great challenge. Computer simulation methods are becoming an increasingly important tool not only in the interpretation of experiments but also in the prediction of membrane protein dynamics. In the present review, we give a brief introduction to molecular modeling techniques currently used to explore protein dynamics on time scales ranging from femtoseconds to microseconds. We then describe a few recent example applications of these techniques to membrane proteins. In conclusion, we also discuss some of the newest developments in simulation methodology that have the potential to further extend the time scale accessible to explore (membrane) protein dynamics.

Proton inhibition of unitary currents of vanilloid receptors

Protons, which are released during inflammation and injury, regulate many receptors and ion channels involved in pain transduction, including capsaicin channels (transient receptor potential vanilloid receptors 1). Whereas extracellular acidification both sensitizes and directly activates the channel, it also causes concomitant reduction of the unitary current amplitudes. Here, we investigate the mechanisms and molecular basis of this inhibitory effect of protons on channel conductance. Single-channel recordings showed that the unitary current amplitudes decreased with extracellular pH in a dose-dependent manner, consistent with a model in which protons bind to a site within the channel with an apparent pKa of ∼6. The inhibition was voltage dependent, ∼65% at −60 mV and 37% at +60 mV when pH was reduced from 7.4 to 5.5. The unitary current amplitudes reached saturation at [K⁺] ≥ 1 M, and notably the maximum amplitudes did not converge with different pHs, inconsistent with a blockade model based on surface charge screening or competitive inhibition of permeating ions. Mutagenesis experiments uncovered two acidic residues critical for proton inhibition, one located at the pore entrance and the other on the pore helix. Based on homology to the KcsA structure, the two acidic residues, along with another basic residue also on the pore helix, could form a triad interacting with each other through extensive hydrogen bonds and electrostatic contacts, suggesting that protons may mediate the interactions between the selectivity filter and pore helix, thereby altering the local structure in the filter region and consequently the conductance of the channel.

Computational analysis of membrane proteins: the largest class of drug targets

Given the key roles of integral membrane proteins as transporters and channels, it is necessary to understand their structures and, hence, mechanisms and regulation at the molecular level. Membrane proteins represent approximately 30% of all proteins of currently sequenced genomes. Paradoxically, however, only approximately 2% of crystal structures deposited in the protein data bank are of membrane proteins, and very few of these are at high resolution (better than 2A). The great disparity between our understanding of soluble proteins and our understanding of membrane proteins is because of the practical problems of working with membrane proteins - specifically, difficulties in expression, purification and crystallization. Thus, computational modeling has been utilized extensively to make crucial advances in understanding membrane protein structure and function.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Biomolecular simulation and modelling: status, progress and prospects

Molecular simulation is increasingly demonstrating its practical value in the investigation of biological systems. Computational modelling of biomolecular systems is an exciting and rapidly developing area, which is expanding significantly in scope. A range of simulation methods has been developed that can be applied to study a wide variety of problems in structural biology and at the interfaces between physics, chemistry and biology. Here, we give an overview of methods and some recent developments in atomistic biomolecular simulation. Some recent applications and theoretical developments are highlighted.

Non-Markovian noise mediated through anomalous diffusion within ion channels

It is evident from a wide range of experimental findings that ion channel gating is inherently stochastic. The issue of "memory effects" (diffusional retardation due to local changes in water viscosity) in ionic flow has been recently addressed using Brownian dynamics simulations. The results presented indicate such memory effects are negligible, unless the diffusional barrier is much higher than that of free solute. In this paper using differential stochastic methods we conclude that the Markovian property of exponential dwell times gives rise to a high barrier, resulting in diffusional memory effects that cannot be ignored in determining ionic flow through channels. We have addressed this question using a generalized Langevin equation that contains a combination of Markovian and non-Markovian processes with different time scales. This approach afforded the development of an algorithm that describes an oscillatory ionic diffusional sequence. The resulting oscillatory function behavior, with exponential decay, was obtained at the weak non-Markovian limit with two distinct time scales corresponding to the processes of ionic diffusion and drift. This will be analyzed further in future studies using molecular dynamics simulations. We propose that the rise of time scales and memory effects is related to differences of shear viscosity in the cytoplasm and extracellular matrix.

The molecular mechanism of toxin-induced conformational changes in a potassium channel: relation to C-type inactivation

Recently, a solid-state NMR study revealed that scorpion toxin binding leads to conformational changes in the selectivity filter of potassium channels. The exact nature of the conformational changes, however, remained elusive. We carried out all-atom molecular dynamics simulations that enabled us to cover the complete pathway of toxin approach and binding, and we validated our simulation results by using solid-state NMR data and electrophysiological measurements. Our structural model revealed a mechanism of cooperative toxin-induced conformational changes that accounts both for the signal changes observed in solid-state NMR and for the tight interaction between KcsA-Kv1.3 and Kaliotoxin. We show that this mechanism is structurally and functionally closely related to recovery from C-type inactivation. Furthermore, our simulations indicate heterogeneity in the binding modes of Kaliotoxin, which might serve to enhance its affinity for KcsA-Kv1.3 further by entropic stabilization.

Conformational changes in the selectivity filter of the open-state KcsA channel: an energy minimization study

Potassium channels switch between closed and open conformations and selectively conduct K(+) ions. There are at least two gates. The TM2 bundle at the intracellular site is the primary gate of KcsA, and rearrangements at the selectivity filter (SF) act as the second gate. The SF blocks ion flow via an inactivation process similar to C-type inactivation of voltage-gated K(+) channels. We recently generated the open-state conformation of the KcsA channel. We found no major, possibly inactivating, structural changes in the SF associated with this massive inner-pore rearrangement, which suggests that the gates might act independently. Here we energy-minimize the open state of wild-type and mutant KcsA, validating in silico structures of energy-minimized SFs by comparison with crystallographic structures, and use these data to gain insight into how mutation, ion depletion, and K(+) to Na(+) substitution influence SF conformation. Both E71 or D80 protonations/mutations and the presence/absence of protein-buried water molecule(s) modify the H-bonding network stabilizing the P-loops, spawning numerous SF conformations. We find that the inactivated state corresponds to conformations with a partially unoccupied or an entirely empty SF. These structures, involving modifications in all four P-loops, are stabilized by H-bonds between amide H and carbonyl O atoms from adjacent P-loops, which block ion passage. The inner portions of the P-loops are more rigid than the outer parts. Changes are localized to the outer binding sites, with innermost site S4 persisting in the inactivated state. Strong binding by Na(+) locally contracts the SF around Na(+), releasing ligands that do not participate in Na(+) coordination, and occluding the permeation pathway. K(+) selectivity primarily appears to arise from the inability of the SF to completely dehydrate Na(+) ions due to basic structural differences between liquid water and the "quasi-liquid" SF matrix.

Conformational dynamics of the KcsA potassium channel governs gating properties

K+ channels conduct and regulate K+ flux across the cell membrane. Several crystal structures and biophysical studies of tetrameric ion channels have revealed many of the structural details of ion selectivity and gating. A narrow pore lined with four arrays of carbonyl groups is responsible for ion selectivity, whereas a conformational change of the four inner transmembrane helices (TM2) is involved in gating. We used NMR to examine full-length KcsA, a prototypical K+ channel, in its open, closed and intermediate states. These studies reveal that at least two conformational states exist both in the selectivity filter and near the C-terminal ends of the TM2 helices. In the ion-conducting open state, we observed rapid structural exchange between two conformations of the filter, presumably of low and high K+ affinity, respectively. Such measurements of millisecond-timescale dynamics reveal the basis for simultaneous ion selection and gating.