Morten Ø. Jensen

Mechanism of Voltage Gating in Potassium Channels

Open and Shut Case Voltage-sensing domains (VSDs) control the activity of voltage-gated ion channels to regulate the ion flow that underlies nerve conduction. Structural and biophysical studies have provided insight into voltage gating; however, understanding has been hindered by the lack of a crystal structure of a fully closed state. Starting from a structure of an open conducting state, a voltage-gated K+ channel, Jensen et al. (p. 229) used all-atom molecular dynamics simulations to show the conformational changes involved in switching to the closed, nonconducting state. Additional simulations revealed the major steps of channel activation. The computational determination of a closed state may guide development of drugs to treat channelopathies associated with this resting state. Molecular dynamics simulations show how a voltage-gated channel closes. The mechanism of ion channel voltage gating—how channels open and close in response to voltage changes—has been debated since Hodgkin and Huxley’s seminal discovery that the crux of nerve conduction is ion flow across cellular membranes. Using all-atom molecular dynamics simulations, we show how a voltage-gated potassium channel (KV) switches between activated and deactivated states. On deactivation, pore hydrophobic collapse rapidly halts ion flow. Subsequent voltage-sensing domain (VSD) relaxation, including inward, 15-angstrom S4-helix motion, completes the transition. On activation, outward S4 motion tightens the VSD–pore linker, perturbing linker–S6-helix packing. Fluctuations allow water, then potassium ions, to reenter the pore; linker-S6 repacking stabilizes the open pore. We propose a mechanistic model for the sodium/potassium/calcium voltage-gated ion channel superfamily that reconciles apparently conflicting experimental data.

Energetics of glycerol conduction through aquaglyceroporin GlpF

Aquaglyceroporin GlpF selectively conducts water and linear polyalcohols, such as glycerol, across the inner membrane of Escherichia coli. We report steered molecular dynamics simulations of glycerol conduction through GlpF, in which external forces accelerate the transchannel conduction in a manner that preserves the intrinsic conduction mechanism. The simulations reveal channel-glycerol hydrogen bonding interactions and the stereoselectivity of the channel. Employing Jarzynskis identity between free energy and irreversible work, we reconstruct the potential of mean force along the conduction pathway through a time series analysis of molecular dynamics trajectories. This potential locates binding sites and barriers inside the channel; it also reveals a low energy periplasmic vestibule suited for efficient uptake of glycerol from the environment.

Identification of two distinct inactive conformations of the β2-adrenergic receptor reconciles structural and biochemical observations

Fully understanding the mechanisms of signaling proteins such as G protein-coupled receptors (GPCRs) will require the characterization of their conformational states and the pathways connecting those states. The recent crystal structures of the β2- and β1-adrenergic receptors in a nominally inactive state constituted a major advance toward this goal, but also raised new questions. Although earlier biochemical observations had suggested that these receptors possessed a set of contacts between helices 3 and 6, known as the ionic lock, which was believed to form a molecular switch for receptor activation, the crystal structures lacked these contacts. The unexpectedly broken ionic lock has raised questions about the true conformation(s) of the inactive state and the role of the ionic lock in receptor activation and signaling. To address these questions, we performed microsecond-timescale molecular dynamics simulations of the β2-adrenergic receptor (β2AR) in multiple wild-type and mutant forms. In wild-type simulations, the ionic lock formed reproducibly, bringing the intracellular ends of helices 3 and 6 together to adopt a conformation similar to that found in inactive rhodopsin. Our results suggest that inactive β2AR exists in equilibrium between conformations with the lock formed and the lock broken, whether or not the cocrystallized ligand is present. These findings, along with the formation of several secondary structural elements in the β2AR loops during our simulations, may provide a more comprehensive picture of the inactive state of the β-adrenergic receptors, reconciling the crystal structures with biochemical studies.

Principles of conduction and hydrophobic gating in K+ channels

We present the first atomic-resolution observations of permeation and gating in a K+ channel, based on molecular dynamics simulations of the Kv1.2 pore domain. Analysis of hundreds of simulated permeation events revealed a detailed conduction mechanism, resembling the Hodgkin–Keynes “knock-on” model, in which translocation of two selectivity filter–bound ions is driven by a third ion; formation of this knock-on intermediate is rate determining. In addition, at reverse or zero voltages, we observed pore closure by a novel “hydrophobic gating” mechanism: A dewetting transition of the hydrophobic pore cavity—fastest when K+ was not bound in selectivity filter sites nearest the cavity—caused the open, conducting pore to collapse into a closed, nonconducting conformation. Such pore closure corroborates the idea that voltage sensors can act to prevent pore collapse into the intrinsically more stable, closed conformation, and it further suggests that molecular-scale dewetting facilitates a specific biological function: K+ channel gating. Existing experimental data support our hypothesis that hydrophobic gating may be a fundamental principle underlying the gating of voltage-sensitive K+ channels. We suggest that hydrophobic gating explains, in part, why diverse ion channels conserve hydrophobic pore cavities, and we speculate that modulation of cavity hydration could enable structural determination of both open and closed channels.

A conserved protonation-dependent switch controls drug binding in the Abl kinase

In many protein kinases, a characteristic conformational change (the “DFG flip”) connects catalytically active and inactive conformations. Many kinase inhibitors—including the cancer drug imatinib—selectively target a specific DFG conformation, but the function and mechanism of the flip remain unclear. Using long molecular dynamics simulations of the Abl kinase, we visualized the DFG flip in atomic-level detail and formulated an energetic model predicting that protonation of the DFG aspartate controls the flip. Consistent with our models predictions, we demonstrated experimentally that the kinetics of imatinib binding to Abl kinase have a pH dependence that disappears when the DFG aspartate is mutated. Our model suggests a possible explanation for the high degree of conservation of the DFG motif: that the flip, modulated by electrostatic changes inherent to the catalytic cycle, allows the kinase to access flexible conformations facilitating nucleotide binding and release.

Microsecond Molecular Dynamics Simulation Shows Effect of Slow Loop Dynamics on Backbone Amide Order Parameters of Proteins

A molecular-level understanding of the function of a protein requires knowledge of both its structural and dynamic properties. NMR spectroscopy allows the measurement of generalized order parameters that provide an atomistic description of picosecond and nanosecond fluctuations in protein structure. Molecular dynamics (MD) simulation provides a complementary approach to the study of protein dynamics on similar time scales. Comparisons between NMR spectroscopy and MD simulations can be used to interpret experimental results and to improve the quality of simulation-related force fields and integration methods. However, apparent systematic discrepancies between order parameters extracted from simulations and experiments are common, particularly for elements of noncanonical secondary structure. In this paper, results from a 1.2 micros explicit solvent MD simulation of the protein ubiquitin are compared with previously determined backbone order parameters derived from NMR relaxation experiments [Tjandra, N.; Feller, S. E.; Pastor, R. W.; Bax, A. J. Am. Chem. Soc. 1995, 117, 12562-12566]. The simulation reveals fluctuations in three loop regions that occur on time scales comparable to or longer than that of the overall rotational diffusion of ubiquitin and whose effects would not be apparent in experimentally derived order parameters. A coupled analysis of internal and overall motion yields simulated order parameters substantially closer to the experimentally determined values than is the case for a conventional analysis of internal motion alone. Improved agreement between simulation and experiment also is encouraging from the viewpoint of assessing the accuracy of long MD simulations.

The Mechanism of Glycerol Conduction in Aquaglyceroporins

BACKGROUND The E. coli glycerol facilitator, GlpF, selectively conducts glycerol and water, excluding ions and charged solutes. The detailed mechanism of the glycerol conduction and its relationship to the characteristic secondary structure of aquaporins and to the NPA motifs in the center of the channel are unknown. RESULTS Molecular dynamics simulations of GlpF reveal spontaneous glycerol and water conduction driven, on a nanosecond timescale, by thermal fluctuations. The bidirectional conduction, guided and facilitated by the secondary structure, is characterized by breakage and formation of hydrogen bonds for which water and glycerol compete. The conduction involves only very minor changes in the protein structure, and cooperativity between the GlpF monomers is not evident. The two conserved NPA motifs are strictly linked together by several stable hydrogen bonds and their asparagine side chains form hydrogen bonds with the substrates passing the channel in single file. CONCLUSIONS A complete conduction of glycerol through the GlpF was deduced from molecular dynamics simulations, and key residues facilitating the conduction were identified. The nonhelical parts of the two half-membrane-spanning segments expose carbonyl groups towards the channel interior, establishing a curve-linear pathway. The conformational stability of the NPA motifs is important in the conduction and critical for selectivity. Water and glycerol compete in a random manner for hydrogen bonding sites in the protein, and their translocations in single file are correlated. The suggested conduction mechanism should apply to the whole family.

Exploring atomic resolution physiology on a femtosecond to millisecond timescale using molecular dynamics simulations.

Discovering the functional mechanisms of biological systems frequently requires information that challenges the spatial and temporal resolution limits of current experimental techniques. Recent dramatic methodological advances have made all-atom molecular dynamics (MD) simulations an ever more

Dynamic control of slow water transport by aquaporin 0: Implications for hydration and junction stability in the eye lens

Aquaporin 0 (AQP0), the most abundant membrane protein in mammalian lens fiber cells, not only serves as the primary water channel in this tissue but also appears to mediate the formation of thin junctions between fiber cells. AQP0 is remarkably less water permeable than other aquaporins, but the structural basis and biological significance of this low permeability remain uncertain, as does the permeability of the protein in a reported junctional form. To address these issues, we performed molecular dynamics (MD) simulations of water transport through membrane-embedded AQP0 in both its (octameric) junctional and (tetrameric) nonjunctional forms. From our simulations, we measured an osmotic permeability for the nonjunctional form that agrees with experiment and found that the distinct dynamics of the conserved, lumen-protruding side chains of Tyr-23 and Tyr-149 modulate water passage, accounting for the slow permeation. The junctional and nonjunctional forms conducted water equivalently, in contrast to a previous suggestion based on static crystal structures that water conduction is lost on junction formation. Our analysis suggests that the low water permeability of AQP0 may help maintain the mechanical stability of the junction. We hypothesize that the structural features leading to low permeability may have evolved in part to allow AQP0 to form junctions that both conduct water and contribute to the organizational structure of the fiber cell tissue and microcirculation within it, as required to maintain transparency of the lens.

Atomic-level simulation of current–voltage relationships in single-file ion channels

The difficulty in characterizing ion conduction through membrane channels at the level of individual permeation events has made it challenging to elucidate the mechanistic principles underpinning this fundamental physiological process. Using long, all-atom simulations enabled by special-purpose hardware, we studied K+ permeation across the KV1.2/2.1 voltage-gated potassium channel. At experimentally accessible voltages, which include the physiological range, the simulated permeation rate was substantially lower than the experimentally observed rate. The current–voltage relationship was also nonlinear but became linear at much higher voltages. We observed permeation consistent with a “knock-on” mechanism at all voltages. At high voltages, the permeation rate was in accordance with our previously reported KV1.2 pore-only simulations, after the simulated voltages from the previous study were recalculated using the correct method, new insight into which is provided here. Including the voltage-sensing domains in the simulated channel brought the linear current–voltage regime closer to the experimentally accessible voltages. The simulated permeation rate, however, still underestimated the experimental rate, because formation of the knock-on intermediate occurred too infrequently. Reducing the interaction strength between the ion and the selectivity filter did not increase conductance. In complementary simulations of gramicidin A, similar changes in interaction strength did increase the observed permeation rate. Permeation nevertheless remained substantially below the experimental value, largely because of infrequent ion recruitment into the pore lumen. Despite the need to apply large voltages to simulate the permeation process, the apparent voltage insensitivity of the permeation mechanism suggests that the direct simulation of permeation at the single-ion level can provide fundamental physiological insight into ion channel function. Notably, our simulations suggest that the knock-on permeation mechanisms in KV1.2 and KcsA may be different.