Alessandro Grottesi

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 Å) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (∼1 Å) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge‐bending of the pore‐lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.

Kv channel gating requires a compatible S4-S5 linker and bottom part of S6, constrained by non-interacting residues.

Voltage-dependent K+ channels transfer the voltage sensor movement into gate opening or closure through an electromechanical coupling. To test functionally whether an interaction between the S4-S5 linker (L45) and the cytoplasmic end of S6 (S6T) constitutes this coupling, the L45 in hKv1.5 was replaced by corresponding hKv2.1 sequence. This exchange was not tolerated but could be rescued by also swapping S6T. Exchanging both L45 and S6T transferred hKv2.1 kinetics to an hKv1.5 background while preserving the voltage dependence. A one-by-one residue substitution scan of L45 and S6T in hKv1.5 further shows that S6T needs to be α-helical and forms a “crevice” in which residues I422 and T426 of L45 reside. These residues transfer the mechanical energy onto the S6T crevice, whereas other residues in S6T and L45 that are not involved in the interaction maintain the correct structure of the coupling.

Insight into the mechanism of inactivation and pH sensitivity in potassium channels from molecular dynamics simulations

Potassium (K (+)) channels can regulate ionic conduction through their pore by a mechanism, involving the selectivity filter, known as C-type inactivation. This process is rapid in the hERG K (+) channel and is fundamental to its physiological role. Although mutations within hERG are known to remove this process, a structural basis for the inactivation mechanism has yet to be characterized. Using MD simulations based on homology modeling, we observe that the carbonyl of the filter aromatic, Phe627, forming the S 0 K (+) binding site, swiftly rotates away from the conduction axis in the wild-type channel. In contrast, in well-characterized non-inactivating mutant channels, this conformational change occurs less frequently. In the non-inactivating channels, interactions with a water molecule located behind the selectivity filter are critical to the enhanced stability of the conducting state. We observe comparable conformational changes in the acid sensitive TASK-1 channel and propose a common mechanism in these channels for regulating efflux of K (+) ions through the selectivity filter.

Molecular simulations and lipid–protein interactions: potassium channels and other membrane proteins

Molecular dynamics simulations may be used to probe the interactions of membrane proteins with lipids and with detergents at atomic resolution. Examples of such simulations for ion channels and for bacterial outer membrane proteins are described. Comparison of simulations of KcsA (an alpha-helical bundle) and OmpA (a beta-barrel) reveals the importance of two classes of side chains in stabilizing interactions with the head groups of lipid molecules: (i) tryptophan and tyrosine; and (ii) arginine and lysine. Arginine residues interacting with lipid phosphate groups play an important role in stabilizing the voltage-sensor domain of the KvAP channel within a bilayer. Simulations of the bacterial potassium channel KcsA reveal specific interactions of phosphatidylglycerol with an acidic lipid-binding site at the interface between adjacent protein monomers. A combination of molecular modelling and simulation reveals a potential phosphatidylinositol 4,5-bisphosphate-binding site on the surface of Kir6.2.

Voltage-gated ion channels.

Further experimental and computational studies are required before we reach a complete structural understanding of the mechanisms of voltage-sensing and voltage-gating of Kv channels. Progress is likely to be made using a range of indirect techniques to complement the structural data that have come from X-ray crystallography. In particular, spectroscopic [7xMolecular architecture of the KvAP voltage-dependent K+ channel in a lipid bilayer. Cuello, L.G., Cortes, D.M., and Perozo, E. Science. 2004; 306: 491–495Crossref | PubMed | Scopus (173)See all References[7] and chemical modification [8xMolecular mechanism of voltage sensor movements in a potassium channel. Elliott, D.J.S., Neale, E.J., Aziz, Q., Dunham, J.P., Munsey, T.S., Hunter, M., and Sivaprasadarao, A. EMBO J. 2005; in pressSee all References[8] studies of Kv channels in situ in lipid bilayers will help to resolve the structure of the resting (closed state) channel. These methods will then have to be deployed in combination with a transmembrane voltage difference to reveal the change in sensor conformation and orientation during the transition from the closed to the open state.It is likely that computational methods will be used to integrate data from these diverse sources in a molecular mechanism. Given the proposed changes in conformation of the voltage sensor upon activation of Kv channels, it is important to characterize the intrinsic flexibility of this domain. Molecular dynamics simulations offer one possibility for exploring the conformational dynamics of the sensor domain (our unpublished data).A number of other techniques may also yield information on the location of the voltage sensor relative to the membrane. For example, toxins that interact with the sensor, for example Vstx1 [9xA membrane-access mechanism of ion channel inhibition by voltage sensor toxins from spider venom. Lee, S.Y. and MacKinnon, R. Nature. 2004; 430: 232–235Crossref | PubMed | Scopus (169)See all References[9], provide valuable probes. Their location – and hence by extension the location of the voltage-sensor – relative to a lipid bilayer may be established via molecular simulations. Similarly, electron microscopy [10xElectron microscopic analysis of KvAP voltage-dependent K+ channels in an open conformation. Jiang, Q.X., Wang, D.N., and MacKinnon, R. Nature. 2004; 430: 806–810Crossref | PubMed | Scopus (82)See all References[10] may reveal the overall shape of Kv molecules trapped in different conformational states. Thus, by combining information from these disparate sources, a complete mechanism of Kv voltage-sensing and gating may emerge.

Early events in protein aggregation: molecular flexibility and hydrophobicity/charge interaction in amyloid peptides as studied by molecular dynamics simulations.

In a previous article (Zbilut et al., Biophys J 2003;85:3544–3557), we demonstrated how an aggregation versus folding choice could be approached considering hydrophobicity distribution and charge. In this work, our aim is highlighting the mutual interaction of charge and hydrophobicity distribution in the aggregation process. Use was made of two different peptides, both derived from a transmembrane protein (amyloid precursor protein; APP), namely, Aβ(1‐28) and Aβ(1‐40). Aβ(1‐28) has a much lower aggregation propensity than Aβ(1‐40). The results obtained by means of molecular dynamics simulations show that, when submitted to the most “aggregation‐prone” environment, corresponding to the isoelectric point and consequently to zero net charge, both peptides acquire their maximum flexibility, but Aβ(1‐40) has a definitely higher conformational mobility than Aβ(1‐28). The absence of a hydrophobic “tail,” which is the most mobile part of the molecule in Aβ(1‐40), is the element lacking in Aβ(1‐28) for obtaining a “fully aggregating” phenotype. Our results suggest that conformational flexibility, determined by both hydrophobicity and charge effect, is the main mechanistic determinant of aggregation propensity. Proteins 2005.

Episodic ataxia type 1 mutations affect fast inactivation of K+ channels by a reduction in either subunit surface expression or affinity for inactivation domain

Episodic ataxia type 1 (EA1) is an autosomal dominant disorder characterized by continuous myokymia and episodic attacks of ataxia. Mutations in the gene KCNA1 that encodes the voltage-gated potassium channel Kv1.1 are responsible for EA1. In several brain areas, Kv1.1 coassembles with Kv1.4, which confers N-type inactivating properties to heteromeric channels. It is therefore likely that the rate of inactivation will be determined by the number of Kv1.4 inactivation particles, as set by the precise subunit stoichiometry. We propose that EA1 mutations affect the rate of N-type inactivation either by reduced subunit surface expression, giving rise to a reduced number of Kv1.1 subunits in heterotetramer Kv1.1-Kv1.4 channels, or by reduced affinity for the Kv1.4 inactivation domain. To test this hypothesis, quantified amounts of mRNA for Kv1.4 or Kv1.1 containing selected EA1 mutations either in the inner vestibule of Kv1.1 on S6 or in the transmembrane regions were injected into Xenopus laevis oocytes and the relative rates of inactivation and stoichiometry were determined. The S6 mutations, V404I and V408A, which had normal surface expression, reduced the rate of inactivation by a decreased affinity for the inactivation domain while the mutations I177N in S1 and E325D in S5, which had reduced subunit surface expression, increased the rate of N-type inactivation due to a stoichiometric increase in the number of Kv1.4 subunits.

Molecular dynamics simulations of a K+ channel blocker: Tc1 toxin from Tityus cambridgei.

Toxins that block voltage‐gated potassium (Kv) channels provide a possible template for improved homology models of the Kv pore. In assessing the interactions of Kv channels and their toxins it is important to determine the dynamic flexibility of the toxins. Multiple 10 ns duration molecular dynamics simulations combined with essential dynamics analysis have been used to explore the flexibility of four different Kv channel‐blocking toxins. Three toxins (Tc1, AgTx and ChTx) share a common fold. They also share a common pattern of conformational dynamics, as revealed by essential dynamics analysis of the simulation results. This suggests that some aspects of dynamic behaviour are conserved across a single protein fold class. In each of these three toxins, the residue exhibiting minimum flexibility corresponds to a conserved lysine residue that is suggested to interact with the filter domain of the channel. Thus, comparative simulations reveal functionally important conservation of molecular dynamics as well as protein fold across a family of related toxins.

Gain-of-function defects of astrocytic Kir4.1 channels in children with autism spectrum disorders and epilepsy

Dysfunction of the inwardly-rectifying potassium channels Kir4.1 (KCNJ10) represents a pathogenic mechanism contributing to Autism-Epilepsy comorbidity. To define the role of Kir4.1 variants in the disorder, we sequenced KCNJ10 in a sample of affected individuals, and performed genotype-phenotype correlations. The effects of mutations on channel activity, protein trafficking, and astrocyte function were investigated in Xenopus laevis oocytes, and in human astrocytoma cell lines. An in vivo model of the disorder was also explored through generation of kcnj10a morphant zebrafish overexpressing the mutated human KCNJ10. We detected germline heterozygous KCNJ10 variants in 19/175 affected children. Epileptic spasms with dysregulated sensory processing represented the main disease phenotype. When investigated on astrocyte-like cells, the p.R18Q mutation exerted a gain-of-function effect by enhancing Kir4.1 membrane expression and current density. Similarly, the p.R348H variant led to gain of channel function through hindrance of pH-dependent current inhibition. The frequent polymorphism p.R271C seemed, instead, to have no obvious functional effects. Our results confirm that variants in KCNJ10 deserve attention in autism-epilepsy, and provide insight into the molecular mechanisms of autism and seizures. Similar to neurons, astrocyte dysfunction may result in abnormal synaptic transmission and electrical discharge, and should be regarded as a possible pharmacological target in autism-epilepsy.

Conformational switch of a flexible loop in human laminin receptor determines laminin-1 interaction

The 37/67-kDa human laminin receptor (LamR) is a cell surface protein that interacts with molecules located in the extra-cellular matrix. In particular, interactions between LamR and laminins play a major role in mediating changes in the cellular environment that affect cell adhesion, neurite outgrowth, tumor growth and metastasis. The exact interaction mode of laminin-1 and LamR is not fully understood. Laminin-1 is thought to bind to LamR through interaction with the so-called peptide G (residues 161–180) and the C-terminal helix (residues 205–229). Here we performed 100-ns atomistic force field-based molecular dynamics simulations to explore the structure and dynamics of LamR related to laminin-1 interactions. Our main finding is that loop 188–197 in the C-terminal region is highly flexible. It undergoes a major change resulting in a conformational switch that partially solvent exposes the R180 residue in the final part of the G peptide. So, R180 could contribute to laminin-1 binding. Projection of the simulations along the first two principal components also confirms the importance of this conformational switch in the LamR. This may be a basic prerequisite to clarify the key structural determinants of the interaction of LamR with laminin-1.