Frédéric Poitevin

One-microsecond molecular dynamics simulation of channel gating in a nicotinic receptor homologue

Recently discovered bacterial homologues of eukaryotic pentameric ligand-gated ion channels, such as the Gloeobacter violaceus receptor (GLIC), are increasingly used as structural and functional models of signal transduction in the nervous system. Here we present a one-microsecond-long molecular dynamics simulation of the GLIC channel pH stimulated gating mechanism. The crystal structure of GLIC obtained at acidic pH in an open-channel form is equilibrated in a membrane environment and then instantly set to neutral pH. The simulation shows a channel closure that rapidly takes place at the level of the hydrophobic furrow and a progressively increasing quaternary twist. Two major events are captured during the simulation. They are initiated by local but large fluctuations in the pore, taking place at the top of the M2 helix, followed by a global tertiary relaxation. The two-step transition of the first subunit starts within the first 50 ns of the simulation and is followed at 450 ns by its immediate neighbor in the pentamer, which proceeds with a similar scenario. This observation suggests a possible two-step domino-like tertiary mechanism that takes place between adjacent subunits. In addition, the dynamical properties of GLIC described here offer an interpretation of the paradoxical properties of a permeable A13′F mutant whose crystal structure determined at 3.15 Å shows a pore too narrow to conduct ions.

AquaSAXS: a web server for computation and fitting of SAXS profiles with non-uniformally hydrated atomic models

Small Angle X-ray Scattering (SAXS) techniques are becoming more and more useful for structural biologists and biochemists, thanks to better access to dedicated synchrotron beamlines, better detectors and the relative easiness of sample preparation. The ability to compute the theoretical SAXS profile of a given structural model, and to compare this profile with the measured scattering intensity, yields crucial structural informations about the macromolecule under study and/or its complexes in solution. An important contribution to the profile, besides the macromolecule itself and its solvent-excluded volume, is the excess density due to the hydration layer. AquaSAXS takes advantage of recently developed methods, such as AquaSol, that give the equilibrium solvent density map around macromolecules, to compute an accurate SAXS/WAXS profile of a given structure and to compare it to the experimental one. Here, we describe the interface architecture and capabilities of the AquaSAXS web server (http://lorentz.dynstr.pasteur.fr/aquasaxs.php).

String method solution of the gating pathways for a pentameric ligand-gated ion channel

Significance High-resolution structures of pentameric ligand-gated ion channels have created an opportunity to discover the mechanisms of rapid synaptic transduction in the brain. This study describes the mechanisms of allosteric channel gating using string method simulations, applied to a complete atomistic ion channel, combined with a transition analysis approach to extract free-energy surfaces from swarms of trajectories. We reproduce pH-modulated activity of the channel, identify the molecular interactions associated with interdomain communication, and quantify the energetics of the gating process. These results provide general mechanistic understanding of the function of pentameric ligand-gated channels, with potential applications in the design of improved anesthetics, neuromodulatory drugs, antiparasitics, and pesticides. Pentameric ligand-gated ion channels control synaptic neurotransmission by converting chemical signals into electrical signals. Agonist binding leads to rapid signal transduction via an allosteric mechanism, where global protein conformational changes open a pore across the nerve cell membrane. We use all-atom molecular dynamics with a swarm-based string method to solve for the minimum free-energy gating pathways of the proton-activated bacterial GLIC channel. We describe stable wetted/open and dewetted/closed states, and uncover conformational changes in the agonist-binding extracellular domain, ion-conducting transmembrane domain, and gating interface that control communication between these domains. Transition analysis is used to compute free-energy surfaces that suggest allosteric pathways; stabilization with pH; and intermediates, including states that facilitate channel closing in the presence of an agonist. We describe a switching mechanism that senses proton binding by marked reorganization of subunit interface, altering the packing of β-sheets to induce changes that lead to asynchronous pore-lining M2 helix movements. These results provide molecular details of GLIC gating and insight into the allosteric mechanisms for the superfamily of pentameric ligand-gated channels.

Comparative Normal Mode Analysis of the Dynamics of DENV and ZIKV Capsids

Key steps in the life cycle of a virus, such as the fusion event as the virus infects a host cell and its maturation process, relate to an intricate interplay between the structure and the dynamics of its constituent proteins, especially those that define its capsid, much akin to an envelope that protects its genomic material. We present a comprehensive, comparative analysis of such interplay for the capsids of two viruses from the flaviviridae family, Dengue (DENV) and Zika (ZIKV). We use for that purpose our own software suite, DD-NMA, which is based on normal mode analysis. We describe the elements of DD-NMA that are relevant to the analysis of large systems, such as virus capsids. In particular, we introduce our implementation of simplified elastic networks and justify their parametrization. Using DD-NMA, we illustrate the importance of packing interactions within the virus capsids on the dynamics of the E proteins of DENV and ZIKV. We identify differences between the computed atomic fluctuations of the E proteins in DENV and ZIKV and relate those differences to changes observed in their high resolution structures. We conclude with a discussion on additional analyses that are needed to fully characterize the dynamics of the two viruses.

Modified Poisson–Boltzmann equations for characterizing biomolecular solvation

Methods for computing electrostatic interactions often account implicitly for the solvent, due to the much smaller number of degrees of freedom involved. In the Poisson–Boltzmann (PB) approach the electrostatic potential is obtained by solving the Poisson–Boltzmann equation (PBE), where the solvent region is modeled as a homogeneous medium with a high dielectric constant. PB however is not exempt of problems. It does not take into account for example the sizes of the ions in the atmosphere surrounding the solute, nor does it take into account the inhomogeneous dielectric response of water due to the presence of a highly charged surface. In this paper we review two major modifications of PB that circumvent these problems, namely the size-modified PB (SMPB) equation and the Dipolar Poisson–Boltzmann Langevin (DPBL) model. In SMPB, steric effects between ions are accounted for with a lattice gas model. In DPBL, the solvent region is no longer modeled as a homogeneous dielectric media but rather as an assembly of self-orienting interacting dipoles of variable density. This model results in a dielectric profile that transits smoothly from the solute to the solvent region as well as in a variable solvent density that depends on the charges of the solute. We show successful applications of the DPBL formalism to computing the solvation free energies of isolated ions in water. Further developments of more accurately modified PB models are discussed.

Nutrient transport suggests an evolutionary basis for charged archaeal surface layer proteins

Surface layers (S-layers) are two-dimensional, proteinaceous, porous lattices that form the outermost cell envelope component of virtually all archaea and many bacteria. Despite exceptional sequence diversity, S-layer proteins (SLPs) share important characteristics such as their ability to form crystalline sheets punctuated with nano-scale pores, and their propensity for charged amino acids, leading to acidic or basic isoelectric points. However, the precise function of S-layers, or the role of charged SLPs and how they relate to cellular metabolism is unknown. Nano-scale lattices affect the diffusion behavior of low-concentration solutes, even if they are significantly smaller than the pore size. Here, we offer a rationale for charged S-layer proteins in the context of the structural evolution of S-layers. Using the ammonia-oxidizing archaea (AOA) as a model for S-layer geometry, and a 2D electrodiffusion reaction computational framework to simulate diffusion and consumption of the charged solute ammonium (NH4+), we find that the characteristic length scales of nanoporous S-layers elevate the concentration of NH4+ in the pseudo-periplasmic space. Our simulations suggest an evolutionary, mechanistic basis for S-layer charge and shed light on the unique ability of some AOA to oxidize ammonia in environments with nanomolar NH4+ availability, with broad implications for comparisons of ecologically distinct populations.

Intermolecular correlations are necessary to explain diffuse scattering from protein crystals

A comprehensive comparison of disorder models indicates that intermolecular correlations must be accounted for to explain the diffuse scattering observed from three protein crystals.

Structure of the 30S ribosomal decoding complex at ambient temperature.

The ribosome translates nucleotide sequences of messenger RNA to proteins through selection of cognate transfer RNA according to the genetic code. To date, structural studies of ribosomal decoding complexes yielding high-resolution data have predominantly relied on experiments performed at cryogenic temperatures. New light sources like the X-ray free electron laser (XFEL) have enabled data collection frommacromolecular crystals at ambient temperature. Here, we report an X-ray crystal structure of the Thermus thermophilus 30S ribosomal subunit decoding complex to 3.45 Å resolution using data obtained at ambient temperature at the Linac Coherent Light Source (LCLS). We find that this ambient-temperature structure is largely consistent with existing cryogenic-temperature crystal structures, with key residues of the decoding complex exhibiting similar conformations, including adenosine residues 1492 and 1493. Minor variations were observed, namely an alternate conformation of cytosine 1397 near the mRNA channel and the A-site. Our serial crystallography experiment illustrates the amenability of ribosomal microcrystals to routine structural studies at ambient temperature, thus overcoming a long-standing experimental limitation to structural studies of RNA and RNA–protein complexes at nearphysiological temperatures.

The Renormalization Group and Its Applications to Generating Coarse-Grained Models of Large Biological Molecular Systems

Understanding the dynamics of biomolecules is the key to understanding their biological activities. Computational methods ranging from all-atom molecular dynamics simulations to coarse-grained normal-mode analyses based on simplified elastic networks provide a general framework to studying these dynamics. Despite recent successes in studying very large systems with up to a 100,000,000 atoms, those methods are currently limited to studying small- to medium-sized molecular systems due to computational limitations. One solution to circumvent these limitations is to reduce the size of the system under study. In this paper, we argue that coarse-graining, the standard approach to such size reduction, must define a hierarchy of models of decreasing sizes that are consistent with each other, i.e., that each model contains the information of the dynamics of its predecessor. We propose a new method, Decimate, for generating such a hierarchy within the context of elastic networks for normal-mode analysis. This method is based on the concept of the renormalization group developed in statistical physics. We highlight the details of its implementation, with a special focus on its scalability to large systems of up to millions of atoms. We illustrate its application on two large systems, the capsid of a virus and the ribosome translation complex. We show that highly decimated representations of those systems, containing down to 1% of their original number of atoms, still capture qualitatively and quantitatively their dynamics. Decimate is available as an OpenSource resource.

Beyond Poisson–Boltzmann: Numerical Sampling of Charge Density Fluctuations

We present a method aimed at sampling charge density fluctuations in Coulomb systems. The derivation follows from a functional integral representation of the partition function in terms of charge density fluctuations. Starting from the mean-field solution given by the Poisson-Boltzmann equation, an original approach is proposed to numerically sample fluctuations around it, through the propagation of a Langevin-like stochastic partial differential equation (SPDE). The diffusion tensor of the SPDE can be chosen so as to avoid the numerical complexity linked to long-range Coulomb interactions, effectively rendering the theory completely local. A finite-volume implementation of the SPDE is described, and the approach is illustrated with preliminary results on the study of a system made of two like-charge ions immersed in a bath of counterions.