Symmetry-restrained molecular dynamics simulations improve homology models of potassium channels

Symmetric Molecular Dynamics

We derive a formulation of molecular dynamics that generates only symmetric configurations. We implement it for all 2D planar and 3D space groups. An atlas of 2D Lennard-Jones crystals under all planar groups is created with symmetric molecular dynamics.

Discrimination of Native-like States of Membrane Proteins with Implicit Membrane-based Scoring Functions

A scoring protocol based on implicit membrane-based scoring functions and a new protocol for optimizing the positioning of proteins inside the membrane was evaluated for its capacity to discriminate native-like states from misfolded decoys. A decoy set previously established by the Baker lab (Proteins: Struct., Funct., Genet. 2006, 62, 1010-1025) was used along with a second set that was generated to cover higher resolution models. The Implicit Membrane Model 1 (IMM1), IMM1 model with CHARMM 36 parameters (IMM1-p36), generalized Born with simple switching (GBSW), and heterogeneous dielectric generalized Born versions 2 (HDGBv2) and 3 (HDGBv3) were tested along with the new HDGB van der Waals (HDGBvdW) model that adds implicit van der Waals contributions to the solvation free energy. For comparison, scores were also calculated with the distance-scaled finite ideal-gas reference (DFIRE) scoring function. Z-scores for native state discrimination, energy vs root-mean-square deviation (RMSD) correlations, and the ability to select the most native-like structures as top-scoring decoys were evaluated to assess the performance of the scoring functions. Ranking of the decoys in the Baker set that were relatively far from the native state was challenging and dominated largely by packing interactions that were captured best by DFIRE with less benefit of the implicit membrane-based models. Accounting for the membrane environment was much more important in the second decoy set where especially the HDGB-based scoring functions performed very well in ranking decoys and providing significant correlations between scores and RMSD, which shows promise for improving membrane protein structure prediction and refinement applications. The new membrane structure scoring protocol was implemented in the MEMScore web server ( http://feiglab.org/memscore ).

Gating energetics of a voltage-dependent K+ channel pore domain

We used targeted molecular dynamics, informed by experimentally determined inter-atomic distances defining the pore region of open and closed states of the KvAP voltage-gated potassium channel, to generate a gating pathway of the pore domain in the absence of the voltage-sensing domains. We then performed umbrella sampling simulations along this pathway to calculate a potential of mean force that describes the free energy landscape connecting the closed and open conformations of the pore domain. The resulting energetic landscape displays three minima, corresponding to stable open, closed, and intermediate conformations with roughly similar stabilities. We found that the extent of hydration of the interior of the pore domain could influence the free energy landscape for pore opening/closing. © 2017 Wiley Periodicals, Inc.

A competing hydrophobic tug on L596 to the membrane core unlatches S4-S5 linker elbow from TRP helix and allows TRPV4 channel to open

We have some generalized physical understanding of how ion channels interact with surrounding lipids but few detailed descriptions on how interactions of particular amino acids with contacting lipids may regulate gating. Here we discovered a structure-specific interaction between an amino acid and inner-leaflet lipid that governs the gating transformations of TRPV4 (transient receptor potential vanilloid type 4). Many cation channels use a S4-S5 linker to transmit stimuli to the gate. At the start of TRPV4's linker helix is leucine 596. A hydrogen bond between the indole of W733 of the TRP helix and the backbone oxygen of L596 secures the helix/linker contact, which acts as a latch maintaining channel closure. The modeled side chain of L596 interacts with the inner lipid leaflet near the polar-nonpolar interface in our model-an interaction that we explored by mutagenesis. We examined the outward currents of TRPV4-expressing Xenopus oocyte upon depolarizations as well as phenotypes of expressing yeast cells. Making this residue less hydrophobic (L596A/G/W/Q/K) reduces open probability [Po; loss-of-function (LOF)], likely due to altered interactions at the polar-nonpolar interface. L596I raises Po [gain-of-function (GOF)], apparently by placing its methyl group further inward and receiving stronger water repulsion. Molecular dynamics simulations showed that the distance between the levels of α-carbons of H-bonded residues L596 and W733 is shortened in the LOFs and lengthened in the GOFs, strengthening or weakening the linker/TRP helix latch, respectively. These results highlight that L596 lipid attraction counteracts the latch bond in a tug-of-war to tune the Po of TRPV4.

The cytoplasmic cage domain of the mechanosensitive channel MscS is a sensor of macromolecular crowding

The cytoplasmic “cage” domain of the bacterial MscS channel senses macromolecular crowding to promote channel inactivation and prevent excessive loss of small osmolytes.

Cells actively regulate the macromolecular excluded volume of the cytoplasm to maintain the reciprocal fraction of free aqueous solution that is optimal for intracellular processes. However, the mechanisms whereby cells sense this critical parameter remain unclear. The mechanosensitive channel of small conductance (MscS channel), which is the major regulator of turgor in bacteria, mediates efflux of small osmolytes in response to increased membrane tension. At moderate sustained tensions produced by a decrease in external osmolarity, MscS undergoes slow adaptive inactivation; however, it inactivates abruptly in the presence of cytoplasmic crowding agents. To understand the mechanism underlying this rapid inactivation, we combined extrapolated and equilibrium molecular dynamics simulations with electrophysiological analyses of MscS mutants to explore possible transitions of MscS and generated models of the resting and inactivated states. Our models suggest that the coupling of the gate formed by TM3 helices to the peripheral TM1–TM2 pairs depends on the axial position of the core TM3 barrel relative to the TM1–TM2 shaft and the state of the associated hollow cytoplasmic domain (“cage”). They also indicate that the tension-driven inactivation transition separates the gate from the peripheral helices and promotes kinks in TM3s at G113 and that this conformation is stabilized by association of the TM3b segment with the β domain of the cage. We found that mutations destabilizing the TM3b–β interactions preclude inactivation and make the channel insensitive to crowding agents and voltage; mutations that strengthen this association result in a stable closed state and silent inactivation. Steered simulations showed that pressure exerted on the cage bottom in the inactivated state reduces the volume of the cage in the cytoplasm and at the same time increases the footprint of the transmembrane domain in the membrane, implying coupled sensitivity to both membrane tension and crowding pressure. The cage, therefore, provides feedback on the increasing crowding that disengages the gate and prevents excessive draining and condensation of the cytoplasm. We discuss the structural mechanics of cells surrounded by an elastic cell wall where this MscS-specific feedback mechanism may be necessary.

TRPV channel-mediated calcium transients in nociceptor neurons are dispensable for avoidance behaviour

Animals need to sense and react to potentially dangerous environments. TRP ion channels participate in nociception, presumably via Ca²⁺ influx, in most animal species. However, the relationship between ion permeation and animals’ nocifensive behaviour is unknown. Here we use an invertebrate animal model with relevance for mammalian pain. We analyse the putative selectivity filter of OSM-9, a TRPV channel, in osmotic avoidance behaviour of Caenorhabditis elegans. Using mutagenized OSM-9 expressed in the head nociceptor neuron, ASH, we study nocifensive behaviour and Ca²⁺ influx. Within the selectivity filter, M⁶⁰¹-F⁶⁰⁹, Y604G strongly reduces avoidance behaviour and eliminates Ca²⁺ transients. Y604F also abolishes Ca²⁺ transients in ASH, while sustaining avoidance behaviour, yet it disrupts behavioral plasticity. Homology modelling of the OSM-9 pore suggests that Y⁶⁰⁴ may assume a scaffolding role. Thus, aromatic residues in the OSM-9 selectivity filter are critical for pain behaviour and ion permeation. These findings have relevance for understanding evolutionary roots of mammalian nociception.

TRPs are calcium-permeable channels involved in the sensing of damaging stimuli but the relationship between calcium influx and pain behaviour has been elusive. Here the authors find that the TRP channel OSM-9 functions as an ion channel in vivo in C. elegans, and establish residues that are critical for worm pain-like behaviour.

Probing the energy landscape of activation gating of the bacterial potassium channel KcsA

The bacterial potassium channel KcsA, which has been crystallized in several conformations, offers an ideal model to investigate activation gating of ion channels. In this study, essential dynamics simulations are applied to obtain insights into the transition pathways and the energy profile of KcsA pore gating. In agreement with previous hypotheses, our simulations reveal a two phasic activation gating process. In the first phase, local structural rearrangements in TM2 are observed leading to an intermediate channel conformation, followed by large structural rearrangements leading to full opening of KcsA. Conformational changes of a highly conserved phenylalanine, F114, at the bundle crossing region are crucial for the transition from a closed to an intermediate state. 3.9 µs umbrella sampling calculations reveal that there are two well-defined energy barriers dividing closed, intermediate, and open channel states. In agreement with mutational studies, the closed state was found to be energetically more favorable compared to the open state. Further, the simulations provide new insights into the dynamical coupling effects of F103 between the activation gate and the selectivity filter. Investigations on individual subunits support cooperativity of subunits during activation gating.

Asymmetric ligand binding facilitates conformational transitions in pentameric ligand-gated ion channels

The anesthetic propofol inhibits the currents of the homopentameric ligand-gated ion channel GLIC, yet the crystal structure of GLIC with five propofol molecules bound symmetrically shows an open-channel conformation. To address this dilemma and determine if the symmetry of propofol binding sites affects the channel conformational transition, we performed a total of 1.5 μs of molecular dynamics simulations for different GLIC systems with propofol occupancies of 0, 1, 2, 3, and 5. GLIC without propofol binding or with five propofol molecules bound symmetrically, showed similar channel conformation and hydration status over multiple replicates of 100-ns simulations. In contrast, asymmetric binding to one, two or three equivalent sites in different subunits accelerated the channel dehydration, increased the conformational heterogeneity of the pore-lining TM2 helices, and shifted the lateral and radial tilting angles of TM2 toward a closed-channel conformation. The results differentiate two groups of systems based on the propofol binding symmetry. The difference between symmetric and asymmetric groups is correlated with the variance in the propofol-binding cavity adjacent to the hydrophobic gate and the force imposed by the bound propofol. Asymmetrically bound propofol produced greater variance in the cavity size that could further elevate the conformation heterogeneity. The force trajectory generated by propofol in each subunit over the course of a simulation exhibits an ellipsoidal shape, which has the larger component tangential to the pore. Asymmetric propofol binding creates an unbalanced force that expedites the channel conformation transitions. The findings from this study not only suggest that asymmetric binding underlies the propofol functional inhibition of GLIC, but also advocate for the role of symmetry breaking in facilitating channel conformational transitions.

Symmetry-restrained flexible fitting for symmetric EM maps

Many large biological macromolecules have inherent structural symmetry, being composed of a few distinct subunits, repeated in a symmetric array. These complexes are often not amenable to traditional high-resolution structural determination methods, but can be imaged in functionally relevant states using cryo-electron microscopy (cryo-EM). A number of methods for fitting atomic-scale structures into cryo-EM maps have been developed, including the molecular dynamics flexible fitting (MDFF) method. However, quality and resolution of the cryo-EM map are the major determinants of a method's success. In order to incorporate knowledge of structural symmetry into the fitting procedure, we developed the symmetry-restrained MDFF method. The new method adds to the cryo-EM map-derived potential further restraints on the allowed conformations of a complex during fitting, thereby improving the quality of the resultant structure. The benefit of using symmetry-based restraints during fitting, particularly for medium to low-resolution data, is demonstrated for three different systems.

Analyses of gating thermodynamics and effects of deletions in the mechanosensitive channel TREK-1: comparisons with structural models

TREK-1, a mechanosensitive K channel from the two-pore family (K(2)P), is involved in protective regulation of the resting potential in CNS neurons and other tissues. The structure of TREK-1 and the basis of its sensitivity to stretch and variety of lipid-soluble factors remain unknown. Using existing K channel structures as modeling templates, TREK-1 was envisioned as a two-fold symmetrical complex with the gate formed primarily by the centrally positioned TM2b helices of the second homologous repeat. Opening was modeled as a conical expansion of the barrel separating TM2b's accompanied by extension of TM2a helices with the cytoplasmic TM2a-TM1b connector. Seeking first experimental support to the models we have accomplished thermodynamic analysis of mouse TREK-1 gating and functional testing of several deletion mutants. The predicted increase of the channel in-plane area (ΔA) of ~5 nm(2) in models was supported by the experimental ΔA of ~4 nm(2) derived from the slope of open probability versus membrane tension in HEK-293T cells and their cytoskeleton-depleted blebs. In response to steps of suction, wild-type channel produced transient currents in cell-attached patches and mostly sustained currents upon patch excision. TREK-1 motifs not present in canonical K channels include divergent cytoplasmic N- and C-termini, and a characteristic 50-residue extracellular loop in the first homologous repeat. Deletion of the extracellular loop (Δ76-124) reduced the average current density in patches, increased spontaneous activity and generated a larger sub-population of high-conductance channels, while activation by tension augmented by arachidonic acid was fully retained. Further deletion of the C-terminal end (Δ76-124/Δ334-411) removed voltage dependency but otherwise produced no additional effect. In an attempt to generate a cysteine-free version of the channel, we mutated two remaining cysteines 159 and 219 in the transmembrane region. C219A did not compromise channel activity, whereas the C159A/S mutants were essentially inactive. Treatment with β-mercaptoethanol suggested that none of these cysteines form functionally-important disulfides.

A practical introduction to molecular dynamics simulations: applications to homology modeling

In this chapter, practical concepts and guidelines are provided for the use of molecular dynamics (MD) simulation for the refinement of homology models. First, an overview of the history and a theoretical background of MD are given. Literature examples of successful MD refinement of homology models are reviewed before selecting the Cytochrome P450 2J2 structure as a case study. We describe the setup of a system for classical MD simulation in a detailed stepwise fashion and how to perform the refinement described in the publication of Li et al. (Proteins 71:938-949, 2008). This tutorial is based on version 11 of the AMBER Molecular Dynamics software package (http://ambermd.org/). However, the approach discussed is equally applicable to any condensed phase MD simulation environment.

Structural models of TREK channels and their gating mechanism

Mechanosensitive TREK channels belong to the family of K2P channels, a family of widely distributed, well modulated channels that uniquely have two similar or identical subunits, each with two TM1-P-TM2 motifs. Our goal is to build viable structural models of TREK channels, as representatives of K2P channels family. The structures available to be used as templates belong to the 2TM channels superfamily. These have low sequence similarity and different structural features: four symmetrically arranged subunits, each having one TM1-P-TM2 motif. Our model building strategy used two subunits of the template (KcsA) to build one subunit of the target (TREK-1). Our models of the Closed channel were adjusted to differ substantially from those of the template, e.g., TM2 of the 2nd repeat is near the axis of the pore whereas TM2 of the 1st repeat is far from the axis. Segments linking the two repeats and immediately following the last TM segment were modeled ab initio as α-helices based on helical periodicities of hydrophobic and hydrophilic residues, highly conserved and poorly conserved residues, and statistically related positions from multiple sequence alignments. The models were further refined by two-fold symmetry-constrained MD simulations using a protocol we developed previously. We also built models of the Open state and suggest a possible tension-activated gating mechanism characterized by helical motion with two-fold symmetry. Our models are consistent with deletion/truncation mutagenesis and thermodynamic analysis of gating described in the accompanying paper.

Hydration properties of mechanosensitive channel pores define the energetics of gating

Opening of ion channels directly by tension in the surrounding membrane appears to be the most ancient and simple mechanism of gating. Bacterial mechanosensitive channels MscL and MscS are the best-studied tension-gated nanopores, yet the key physical factors that define their gating are still hotly debated. Here we present estimations, simulations and experimental results showing that hydration of the pore might be one of the major parameters defining the thermodynamics and kinetics of mechanosensitive channel gating. We associate closing of channel pores with complete dehydration of the hydrophobic gate (occlusion by 'vapor lock') and formation of two water-vapor interfaces above and below the constriction. The opening path is the expansion of these interfaces, ultimately leading to wetting of the hydrophobic pore, which does not appear to be the exact reverse of the closing path, thus producing hysteresis. We discuss specifically the role of polar groups (glycines) buried in narrow closed conformations but exposed in the open states that change the wetting characteristics of the pore lining and stabilize conductive states of the channels.