Un-gating and allosteric modulation of a pentameric ligand-gated ion channel captured by molecular dynamics.

Pentameric ligand-gated ion channels (pLGICs) mediate intercellular communication at synapses through the opening of an ion pore in response to the binding of a neurotransmitter. Despite the increasing availability of high-resolution structures of pLGICs, a detailed understanding of the functional isomerization from closed to open (gating) and back is currently missing. Here, we provide the first atomistic description of the transition from open to closed (un-gating) in the glutamate-gated chloride channel (GluCl) from Caenorhabditis Elegans. Starting with the active-state structure solved in complex with the neurotransmitter L-glutamate and the positive allosteric modulator (PAM) ivermectin, we analyze the spontaneous relaxation of the channel upon removal of ivermectin by explicit solvent/membrane Molecular Dynamics (MD) simulations. The µs-long trajectories support the conclusion that ion-channel deactivation is mediated by two distinct quaternary transitions, i.e. a global receptor twisting followed by the radial expansion (or blooming) of the extracellular domain. At variance with previous models, we show that pore closing is exclusively regulated by the global twisting, which controls the position of the ß1-ß2 loop relative to the M2-M3 loop at the EC/TM domain interface. Additional simulations with L-glutamate restrained to the crystallographic binding mode and ivermectin removed indicate that the same twisting isomerization is regulated by agonist binding at the orthosteric site. These results provide a structural model for gating in pLGICs and suggest a plausible mechanism for the pharmacological action of PAMs in this neurotransmitter receptor family. The simulated un-gating converges to the X-ray structure of GluCl resting state both globally and locally, demonstrating the predictive character of state-of-art MD simulations.

On the permeation of large organic cations through the pore of ATP-gated P2X receptors.

Pore dilation is thought to be a hallmark of purinergic P2X receptors. The most commonly held view of this unusual process posits that under prolonged ATP exposure the ion pore expands in a striking manner from an initial small-cation conductive state to a dilated state, which allows the passage of larger synthetic cations, such as N-methyl-d-glucamine (NMDG+). However, this mechanism is controversial, and the identity of the natural large permeating cations remains elusive. Here, we provide evidence that, contrary to the time-dependent pore dilation model, ATP binding opens an NMDG+-permeable channel within milliseconds, with a conductance that remains stable over time. We show that the time course of NMDG+ permeability superimposes that of Na+ and demonstrate that the molecular motions leading to the permeation of NMDG+ are very similar to those that drive Na+ flow. We found, however, that NMDG+ "percolates" 10 times slower than Na+ in the open state, likely due to a conformational and orientational selection of permeating molecules. We further uncover that several P2X receptors, including those able to desensitize, are permeable not only to NMDG+ but also to spermidine, a large natural cation involved in ion channel modulation, revealing a previously unrecognized P2X-mediated signaling. Altogether, our data do not support a time-dependent dilation of the pore on its own but rather reveal that the open pore of P2X receptors is wide enough to allow the permeation of large organic cations, including natural ones. This permeation mechanism has considerable physiological significance.

Space and Time Evolution of the Electrostatic Potential During the Activation of a Visual Pigment.

Animal and microbial retinal proteins employ the Schiff base of retinal as their chromophore. Here, the possible consequences of the charge translocation associated with the light-induced dynamics of the chromophore of a visual opsin are investigated along a representative semiclassical trajectory. We show that the evolution of the electrostatic potential projected by the chromophore onto the surrounding protein displays intense but topographically localized sudden variations in proximity of the decay region. pKa calculations carried out on selected snapshots used as probes, indicate that the only residue which may be sensitive to the electrostatic potential shift is Glu181. Accordingly, our results suggest that the frail Tyr191/268-Glu181-Wat2-Ser186 hydrogen bond network may be perturbed by the transient variations of the electrostatic potential.

Photo-switchable tweezers illuminate pore-opening motions of an ATP-gated P2X ion channel.

P2X receptors function by opening a transmembrane pore in response to extracellular ATP. Recent crystal structures solved in apo and ATP-bound states revealed molecular motions of the extracellular domain following agonist binding. However, the mechanism of pore opening still remains controversial. Here we use photo-switchable cross-linkers as 'molecular tweezers' to monitor a series of inter-residue distances in the transmembrane domain of the P2X2 receptor during activation. These experimentally based structural constraints combined with computational studies provide high-resolution models of the channel in the open and closed states. We show that the extent of the outer pore expansion is significantly reduced compared to the ATP-bound structure. Our data further reveal that the inner and outer ends of adjacent pore-lining helices come closer during opening, likely through a hinge-bending motion. These results provide new insight into the gating mechanism of P2X receptors and establish a versatile strategy applicable to other membrane proteins.

A gating mechanism of pentameric ligand-gated ion channels.

Pentameric ligand-gated ion channels (pLGICs) play a central role in intercellular communication in the nervous system and are involved in fundamental processes such as attention, learning, and memory. They are oligomeric protein assemblies that convert a chemical signal into an ion flux through the postsynaptic membrane, but the molecular mechanism of gating ions has remained elusive. Here, we present atomistic molecular dynamics simulations of the prokaryotic channels from Gloeobacter violaceus (GLIC) and Erwinia chrysanthemi (ELIC), whose crystal structures are thought to represent the active and the resting states of pLGICs, respectively, and of the eukaryotic glutamate-gated chloride channel from Caenorhabditis elegans (GluCl), whose open-channel structure was determined complexed with the positive allosteric modulator ivermectin. Structural observables extracted from the trajectories of GLIC and ELIC are used as progress variables to analyze the time evolution of GluCl, which was simulated in the absence of ivermectin starting from the structure with bound ivermectin. The trajectory of GluCl with ivermectin removed shows a sequence of structural events that couple agonist unbinding from the extracellular domain to ion-pore closing in the transmembrane domain. Based on these results, we propose a structural mechanism for the allosteric communication leading to deactivation/activation of the GluCl channel. This model of gating emphasizes the coupling between the quaternary twisting and the opening/closing of the ion pore and is likely to apply to other members of the pLGIC family.

Unwrapping of nucleosomal DNA ends: a multiscale molecular dynamics study.

To permit access to DNA-binding proteins involved in the control and expression of the genome, the nucleosome undergoes structural remodeling including unwrapping of nucleosomal DNA segments from the nucleosome core. Here we examine the mechanism of DNA dissociation from the nucleosome using microsecond timescale coarse-grained molecular dynamics simulations. The simulations exhibit short-lived, reversible DNA detachments from the nucleosome and long-lived DNA detachments not reversible on the timescale of the simulation. During the short-lived DNA detachments, 9 bp dissociate at one extremity of the nucleosome core and the H3 tail occupies the space freed by the detached DNA. The long-lived DNA detachments are characterized by structural rearrangements of the H3 tail including the formation of a turn-like structure at the base of the tail that sterically impedes the rewrapping of DNA on the nucleosome surface. Removal of the H3 tails causes the long-lived detachments to disappear. The physical consistency of the CG long-lived open state was verified by mapping a CG structure representative of this state back to atomic resolution and performing molecular dynamics as well as by comparing conformation-dependent free energies. Our results suggest that the H3 tail may stabilize the nucleosome in the open state during the initial stages of the nucleosome remodeling process.

Cys(x)His(y)-Zn2+ interactions: possibilities and limitations of a simple pairwise force field.

In zinc proteins, the Zn2+ cation frequently binds with a tetrahedral coordination to cysteine and histidine side chains. We examine the possibilities and limitations of a classical, pairwise force field for molecular dynamics of such systems. Hartree Fock and density functional calculations are used to obtain geometries, charge distributions, and association energies of side chain analogues bound to Zn2+. Both ionized and neutral cysteines are considered. Two parameterizations are obtained, then tested and compared through molecular dynamics simulations of two small, homologous proteins in explicit solvent: Protein Kinase C and the Cysteine Rich Domain (CRD) of Raf, which have two Cys3His-Zn2+ groups each. The lack of explicit polarizability and charge transfer in the force field leads to poor accuracy for the association energies, and to parameters--including the zinc charge, that depend on the number of bound cysteines and their protonation state. Nevertheless, the structures sampled with the best parameterization are in good overall agreement with experiment, and have zinc coordination geometries compatible with related structures in the Cambridge Structural Database and the Protein Data Bank. Non-optimized parameters lead to poorer structures. This suggests that while a simple force field is not appropriate for processes involving exchange between water and amino acids in the zinc coordination sphere (e.g. protein unfolding), it can be useful for equilibrium simulations of stable Cys3His zinc fingers.

The influence of a transmembrane pH gradient on protonation probabilities of bacteriorhodopsin: the structural basis of the back-pressure effect.

Bacteriorhodopsin pumps protons across a membrane using the energy of light. The proton pumping is inhibited when the transmembrane proton gradient that the protein generates becomes larger than four pH units. This phenomenon is known as the back-pressure effect. Here, we investigate the structural basis of this effect by predicting the influence of a transmembrane pH gradient on the titration behavior of bacteriorhodopsin. For this purpose we introduce a method that accounts for a pH gradient in protonation probability calculations. The method considers that in a transmembrane protein, which is exposed to two different aqueous phases, each titratable residue is accessible for protons from one side of the membrane depending on its hydrogen-bond pattern. This method is applied to several ground-state structures of bacteriorhodopsin, which residues already present complicated titration behaviors in the absence of a proton gradient. Our calculations show that a pH gradient across the membrane influences in a non-trivial manner the protonation probabilities of six titratable residues which are known to participate in the proton transfer: D85, D96, D115, E194, E204, and the Schiff base. The residues connected to one side of the membrane are influenced by the pH on the other side because of their long-range electrostatic interactions within the protein. In particular, D115 senses the pH at the cytoplasmic side of the membrane and transmits this information to D85 and the Schiff base. We propose that the strong electrostatic interactions found between D85, D115, and the Schiff base as well as the interplay of their respective protonation states under the influence of a transmembrane pH gradient are responsible for the back-pressure effect on bacteriorhodopsin.

Cys(x)His(y)-Zn2+ interactions: thiol vs. thiolate coordination.

In zinc proteins, the Zn2+ cation frequently binds with a tetrahedral coordination to cysteine and histidine side chains, for example, in many DNA-binding proteins, where it plays primarily a structural role. We examine the possibility of thiolate protonation in Cys(x)His(y)-Zn2+ groups, both in proteins and in solution, through a combination of theoretical calculations and database analysis. Seventy-five percent of the thiolate-coordinated zincs in the Cambridge Structural Database are tetrahedral, while di-alkanethiol coordination always involves five or more ligands. Ab initio quantum calculations are performed on (ethanethiol/thiolate)(3)imidazole-Zn2+ complexes in vacuum, yielding geometries and gas phase basicities. Protonating one (respectively two) thiolates increases the Zn-S(thiol) distance by 0.4 A (respectively 0.3 A), providing a structural marker for protonation. The stabilities of the complexes in solution are compared by combining the gas phase basicities with continuum dielectric solvation calculations. In a continuum solvent with permittivity epsilon = 4, 20, or 80, one of three thiolates is predicted to be protonated at neutral pH. By extension, Cys4-Zn2+ groups are expected to be protonated in the same conditions. In contrast, most Cys3His and Cys4 geometries in the Protein Data Bank (PDB) appear consistent with all-thiolate Zn2+ coordination. This apparent discrepancy is resolved by two recent surveys of zinc protein structures, which suggest that these all-thiolate sites are stabilized by charged and polar groups nearby in the protein, thus overcoming their intrinsic instability. However, the experimental resolution is not sufficient in all the PDB structures to rule out a thiol/thiolate mixture, and protonated thiolates may occur in some proteins not solved at high resolution or not represented in the PDB, as suggested by recent mass spectrometry experiments; this possibility should be allowed for in X-ray structure refinement.