Search

Developments in structural biology mean that the number of different ion channel structures has increased significantly in recent years. Structures of ion channels enable us to rationalize how mutations may lead to channelopathies. However, determining the structures of ion channels is still not trivial, especially as they necessarily exist in many distinct functional states. Therefore, the use of computational modelling can provide complementary information that can refine working hypotheses of both wild type and mutant ion channels. The simplest but still powerful tool is homology modelling. Many structures are available now that can provide suitable templates for many different types of ion channels, allowing a full three-dimensional interpretation of mutational effects. These structural models, and indeed the structures themselves obtained by X-ray crystallography, and more recently cryo-electron microscopy, can be subjected to molecular dynamics simulations, either as a tool to help explore the conformational dynamics in detail or simply as a means to refine the models further. Here we review how these approaches have been used to improve our understanding of how diseases might be linked to specific mutations in ion channel proteins. This article is part of the Special Issue entitled 'Channelopathies.', (Copyright © 2017 The Authors. Published by Elsevier Ltd.. All rights reserved.)

γ-Aminobutyric acid receptors (GABA A Rs) are vital for controlling excitability in the brain. This is emphasized by the numerous neuropsychiatric disorders that result from receptor dysfunction. A critical component of most native GABA A Rs is the α subunit. Its transmembrane domain is the target for many modulators, including endogenous brain neurosteroids that impact anxiety, stress and depression, and for therapeutic drugs, such as general anesthetics. Understanding the basis for the modulation of GABA A R function requires high-resolution structures. Here we present the first atomic structures of a GABA A R chimera at 2.8-Å resolution, including those bound with potentiating and inhibitory neurosteroids. These structures define new allosteric binding sites for these modulators that are associated with the α-subunit transmembrane domain. Our findings will enable the exploitation of neurosteroids for therapeutic drug design to regulate GABA A Rs in neurological disorders.

Protease inhibition by metastable serine protease inhibitors (serpins) is mediated by one of the largest functional intradomain conformational changes known in biology. In this extensive structural rearrangement, protease-serpin complex formation triggers cleavage of the serpin reactive center loop (RCL), its subsequent insertion into central β-sheet A, and covalent trapping of the target protease. In this study, we present the first detailed accelerated molecular dynamics simulation of the insertion of the fully cleaved RCL in α-1-antitrypsin (α 1 AT), the archetypal member of the family of human serpins. Our results reveal internal water pathways that allow the initial incorporation of side chains of RCL residues into the protein interior. We observed structural plasticity of the helix F (hF) element that blocks the RCL path in the native state, which is in excellent agreement with previous experimental reports. Furthermore, the simulation suggested a novel role of hF and the connected turn (thFs3A) as chaperones that support the insertion process by reducing the conformational space available to the RCL. Transient electrostatic interactions of RCL residues potentially fine-tune the serpin inhibitory activity. On the basis of our simulation, we generated the α 1 AT mutants K168E, E346K, and K168E/E346K and analyzed their inhibitory activity along with their intrinsic stability and heat-induced polymerization. Remarkably, the E346K mutation exhibited enhanced inhibitory activity along with an increased rate of premature structural collapse (polymerization), suggesting a significant role of E346 in the gatekeeping of the strain in the metastable native state.

Human islet amyloid polypeptide (hIAPP) is a 37-residue peptide hormone, which upon misfolding changes from the physiologically active monomer into pathological amyloid fibril aggregates in the pancreas of type 2 diabetes mellitus patients. During this process, the insulin-producing pancreatic β-cells are damaged; however, the underlying mechanism of this mode of cytotoxicity remains elusive. It is known that anionic lipids accelerate amyloid fibril formation, implicating the importance of the cellular membrane in the process, and that a pH close to the level in the β-cell secretory granules (pH 5.5) inhibits amyloid fibril formation. Using all-atom molecular dynamics simulations, we have investigated the membrane-associated monomer state of α-helical hIAPP, analyzed specific interactions of hIAPP with a mixed anionic-zwitterionic lipid membrane and examined the influence of pH on the structure and dynamics of hIAPP and its interaction with the membrane. We find that hIAPP primarily interacts with the membrane by forming favorable interactions between anionic lipids and the positively charged residues in the N-terminal part of the peptide. Rationalizing experimental findings, the simulations show that the N-terminal part of the peptide interacts with the membrane in the lipid headgroup region. At neutral pH, the C-terminal part of the peptide, which contains the residues that initiate fibril formation, displays a highly dynamic, unfolded state, which interacts with the membrane significantly less than the N-terminal part. Such an unfolded form can be proposed to contribute to the acceleration of fibril formation. At low pH, protonation of His18 mediates a stronger interaction of the C-terminal part with the membrane, resulting in the immobilization of the C-terminal part on the membrane surface that might constitute a mechanism by which low pH inhibits fibril formation.

Current enhanced sampling molecular dynamics methods for studying large conformational changes in proteins suffer from certain limitations. These include, among others, the need for user defined collective variables, the prerequisite of both start and end point structures of the conformational change, and the need for a priori knowledge of the amount by which to boost specific parts of the potential. In this paper, a framework is proposed for a molecular dynamics method for studying ligand-induced conformational changes, in which the nonbonded interactions between the ligand and the protein are used to calculate a biasing force. The method requires only a single input structure, and does not entail the use of collective variables. We provide a proof-of-concept for accelerating conformational changes in three simple test molecules, as well as promising results for two proteins known to undergo domain closure upon ligand binding. For the ribose-binding protein, backbone root-mean-square deviations as low as 0.75 Å compared to the crystal structure of the closed conformation are obtained within 50 ns simulations, whereas no domain closures are observed in unbiased simulations. A skewed closed structure is obtained for the glutamine-binding protein at high bias values, indicating that specific protein-ligand interactions might suppress important protein-protein interactions.

The human transforming growth factor β-induced protein (TGFBIp) is involved in several types of corneal dystrophies where protein aggregation and amyloid fibril formation severely impair vision. Most disease-causing mutations are located in the last of four homologous fasciclin-1 (FAS1) domains of the protein, and it has been shown that when isolated, the fourth FAS1 domain (FAS1-4) mimics the behavior of full-length TGFBIp. In this study, we use molecular dynamics simulations and principal component analysis to study the wild-type FAS1-4 domain along with three disease-causing mutations (R555W, R555Q, and A546T) to decipher any internal difference in dynamical properties of the domains that may explain their varied stabilities and aggregation properties. In addition, we use a protein-protein docking method in combination with chemical cross-linking experiments and mass spectrometry of the cross-linked species to obtain information about interaction faces between identical FAS1-4 domains. The results show that the pathogenic mutations A546T and R555W affect the packing in the hydrophobic core of FAS1-4 in different directions. We further show that the FAS1-4 monomers associate using their β-rich regions, consistent with peptides observed to be part of the amyloid fibril core in lattice corneal dystrophy patients.

An epimer of the known glycosidase inhibitor noeurostegine, galacto-noeurostegine, was synthesised in 21 steps from levoglucosan and found to be a potent, competitive and highly selective galactosidase inhibitor of Aspergillus oryzae β-galactosidase. Galacto-noeurostegine was not found to be an inhibitor of green coffee bean α-galactosidase, yeast α-glucosidase and E. coli β-galactosidase, whereas potent but non-competitive inhibition against sweet almond β-glucosidase was established. The 2-deoxy-galacto-noeurostegine analogue was also prepared and found to be a less potent inhibitor of the same enzymes.