Protein conformational changes are at the heart of cell functions, from signalling to ion transport. However, the transient nature of the intermediates along transition pathways hampers their experimental detection, making the underlying mechanisms elusive. Here we retrieve dynamic information on the actual transition routes from principal component analysis (PCA) of structurally-rich ensembles and, in combination with coarse-grained simulations, explore the conformational landscapes of five well-studied proteins. Modelling them as elastic networks in a hybrid elastic-network Brownian dynamics simulation (eBDIMS), we generate trajectories connecting stable end-states that spontaneously sample the crystallographic motions, predicting the structures of known intermediates along the paths. We also show that the explored non-linear routes can delimit the lowest energy passages between end-states sampled by atomistic molecular dynamics. The integrative methodology presented here provides a powerful framework to extract and expand dynamic pathway information from the Protein Data Bank, as well as to validate sampling methods in general.

Accurate force‐field (FF) parameters are key to reliable prediction of properties obtained from molecular modeling (MM) and molecular dynamics (MD) simulations. With ever‐widening applicability of MD simulations, robust parameters need to be generated for a wider range of chemical species. The CHARMM General Force Field program (CGenFF, https://cgenff.umaryland.edu/) is a tool for obtaining initial parameters for a given small molecule based on analogy with the available CGenFF parameters. However, improvement of these parameters is often required and performing their optimization remains tedious and time consuming. In addition, tools for optimization of small molecule parameters in the context of the Drude polarizable FF are not yet available. To overcome these issues, the FFParam package has been designed to facilitate the parametrization process. The package includes a graphical user interface (GUI) created using Qt libraries. FFParam supports Gaussian and Psi4 for performing quantum mechanical calculations and CHARMM and OpenMM for MM calculations. A Monte Carlo simulated annealing (MCSA) algorithm has been implemented for automated fitting of partial atomic charge, atomic polarizabilities and Thole scale parameters. The LSFITPAR program is called for automated fitting of bonded parameters. Accordingly, FFParam provides all the features required for generation and analysis of CHARMM and Drude FF parameters for small molecules. FFParam‐GUI includes a text editor, graph plotter, molecular visualization, and text to table converter to meet various requirements of the parametrization process. It is anticipated that FFParam will facilitate wider use of CGenFF as well as promote future use of the Drude polarizable FF.

Improving our understanding of the mechanisms and effects of anesthetics is a critically important part of neuroscience. The currently dominant theory is that anesthetics and similar molecules act by binding to Cys-loop receptors in the postsynaptic terminal of nerve cells and potentiate or inhibit their function. Although structures for some of the most important mammalian channels have still not been determined, a number of important results have been derived from work on homologous cationic channels in bacteria. However, partly due to the lack of a nervous system in bacteria, there are a number of questions about how these results relate to higher organisms. The recent determination of a structure of the eukaryotic chloride channel, GluCl, is an important step toward accurate modeling of mammalian channels, because it is more similar in function to human Cys-loop receptors such as GABAAR or GlyR. One potential issue with using GluCl to model other receptors is the presence of the large ligand ivermectin (IVM) positioned between all five subunits. Here, we have performed a series of microsecond molecular simulations to study how the dynamics and structure of GluCl change in the presence versus absence of IVM. When the ligand is removed, subunits move at least 2 Å closer to each other compared to simulations with IVM bound. In addition, the pore radius shrinks to 1.2 Å, all of which appears to support a model where IVM binding between subunits stabilizes an open state, and that the relaxed nonIVM conformations might be suitable for modeling other channels. Interestingly, the presence of IVM also has an effect on the structure of the important loop C located at the neurotransmitter-binding pocket, which might help shed light on its partial agonist behavior.

Pentameric ligand-gated ion channels or Cys-loop receptors are responsible for fast inhibitory or excitatory synaptic transmission. The antipsychotic compound chlorpromazine is a widely used tool to probe the ion channel pore of the nicotinic acetylcholine receptor, which is a prototypical Cys-loop receptor. In this study, we determine the molecular determinants of chlorpromazine binding in the Erwinia ligand-gated ion channel (ELIC). We report the X-ray crystal structures of ELIC in complex with chlorpromazine or its brominated derivative bromopromazine. Unexpectedly, we do not find a chlorpromazine molecule in the channel pore of ELIC, but behind the β8-β9 loop in the extracellular ligand-binding domain. The β8-β9 loop is localized downstream from the neurotransmitter binding site and plays an important role in coupling of ligand binding to channel opening. In combination with electrophysiological recordings from ELIC cysteine mutants and a thiol-reactive derivative of chlorpromazine, we demonstrate that chlorpromazine binding at the β8-β9 loop is responsible for receptor inhibition. We further use molecular-dynamics simulations to support the X-ray data and mutagenesis experiments. Together, these data unveil an allosteric binding site in the extracellular ligand-binding domain of ELIC. Our results extend on previous observations and further substantiate our understanding of a multisite model for allosteric modulation of Cys-loop receptors.ligand-gated ion channel | X-ray crystallography | allosteric modulation | Cys-loop receptor | nicotinic acetylcholine receptor C hlorpromazine (CPZ) (Fig. 1), a phenothiazine-derived antipsychotic drug, was introduced in psychiatry in the early 1950s, revolutionizing the treatment of psychotic disorders (1, 2). The main mechanism of action of CPZ consists in the blockage of dopamine receptors (2-4), but the numerous side effects associated with this drug indicate that it interacts with other physiologically relevant targets. CPZ was indeed shown to interfere with several voltage-and ligand-gated channels: it inhibits neuronal voltage-gated K + channels (5-7), BK Ca channels (8), and the human α 1E subunit-mediated Ca 2+ channels (9); CPZ was also shown to inhibit GABAergic currents (10, 11), specifically through GABA A receptors (GABA A Rs) (12), and to inhibit serotonin type-3 receptors (5-HT 3 Rs) (13, 14) and nicotinic acetylcholine receptors (nAChRs) (15, 16), members of the Cys-loop receptor family.The Cys-loop receptor family is composed of membranespanning ligand-gated ion channels that are responsible for fast excitatory or inhibitory synaptic neurotransmission. They are composed of five identical or nonidentical subunits, each of them comprising an N-terminal extracellular domain, which contains the neurotransmitter binding site, four transmembrane helices, that when assembled allow ions to pass through the membrane, and an intracellular domain, responsible for channel conductance, receptor modulation, and trafficking (17, 18). Initial structural insight into the ...