Advances in the Field

Highlights and overview of selected recent Consortium-sponsored publications

Biomolecular DNP- supported NMR spectroscopy using site directed spin labeling

Authors: Elwin A.W. van der Cruijsen, Eline J. Koers, Claire Sauvee, Raymond E. Hulse, Markus Weingarth, Olivier Ouari, Eduardo Perozo, Paul Tordo, and Marc Baldus.

In the press, Journal of the American Chemical Society. (Consortium project: Dynamics of Ion Permeation and Conformation Coupling in KcsA)

Molecular probes that combine the benefits of enhanced spectroscopic sensitivity with site-specific localization have significantly expanded our ability to track molecular structure and function in applications ranging from cell biology to material science . In the field of magnetic resonance, dynamic nuclear polarization (DNP) has become a widely usable method to significantly enhance overall spectroscopic sensitivity in NMR and MRI. Here, the authors show that DNP can be established by creating local spin clusters via site-directed spin labeling using mono- or biradicals. Applied to a membrane-embedded potassium channel, we show that this approach can significantly enhance NMR sensitivity while preserving the intrinsic spectroscopic properties of (bi)radicals as paramagnetic relaxation enhancers.

The results suggest that the creation of local spin clusters can generate sizable DNP enhancements while preserving the intrinsic benefits of PRE-based NMR approaches. Our results are consistent with the idea that the magnitude in DNP enhancement are highly dependent on the nearest neighbor electron-electron distances.

13C NMR detects conformational change in the 100-kD membrane transporter Cl-C-ec1

Authors: Sherwin J. Abraham, Ricky C. Cheng, Thomas A. Chew, Chandra M. Khantwal, Corey W. Liu, Shimei Gong, Robert K. Nakamoto, and Merritt Maduke.

J Biomol NMR. 2015 Jan 29. (Consortium project: Conformational Dynamics in the CLC Channel Transporter Family. Consortium core: Membrane Protein Production)

Members of the Cl-C (‘‘Chloride-Channel’’) family play central roles in cardiovascular, neuronal, bone, and epithelial function. Cl-C transporters catalyze the exchange of Cl--for H+ across cellular membranes. To do so, they must couple Cl-- and H+ binding and unbinding to protein conformational change. However, the sole conformational changes distinguished crystallographically are small movements of a glutamate side chain that locally gates the ion transport pathways. Therefore, our understanding of whether and how global protein dynamics contribute to the exchange mechanism has been severely limited. To overcome the limitations of crystallography, the authors used solution-state 13Cmethyl NMR with labels on methionine, lysine, and engineered cysteine residues to investigate substrate (H+) dependent conformational change outside the restraints of crystallization. They show that methyl labels in several regions report H+-dependent spectral changes. They  identify one of these regions as Helix R, a helix that extends from the center of the protein, where it forms the part of the inner gate to the Cl-–permeation pathway, to the extracellular solution. The H+-dependent spectral change does not occur when a label is positioned just beyond Helix R, on the unstructured C-terminus of the protein.

Together, the results suggest that H+ binding is mechanistically coupled to Cl-osing of the intracellular access-pathway for Cl--. These studies set the stage for investigating the structural details and dynamics of this change.

Room-Temperature Distance Measurements of Immobilized Spin-Labeled Protein by DEER/PELDOR

Authors: Virginia Meyer, Michael A. Swanson, Laura J. Clouston, Przemyslaw J. Boratynski, Richard A. Stein, Hassane S. Mchaourab, Andrzej Rajca, Sandra S. Eaton, and Gareth R. Eaton.

Biophys J. 2015 Mar 10;108(5):1213-9. (Consortium Core: Membrane Protein Production)

Nitroxide spin labels are used for double electron-electron resonance (DEER) measurements of distances between sites in biomolecules. Rotation of gem-dimethyls in commonly used nitroxides causes spin echo dephasing times (Tm) to be too short to perform DEER measurements at temperatures between ~80 and 295 K, even in immobilized samples. A spirocyclohexyl spin label has been prepared that has longer Tm between 80 and 295 K in immobilized samples than conventional labels. Two of the spirocyclohexyl labels were attached to sites on T4 lysozyme introduced by site-directed spin labeling. Interspin distances up to ~4 nm were measured by DEER at temperatures up to 160 K in water/glycerol glasses. In a glassy trehalose matrix the Tm for the doubly labeled T4 lysozyme was long enough to measure an interspin distance of 3.2 nm at 295 K, which could not be measured for the same protein labeled with the conventional 1-oxyl-2,2,5,5-tetramethyl-3-pyrroline-3-(methyl) methanethio-sulfonate label.

The leap from 80 K DEER measurements to room temperature is an important step toward DEER measurement in physiological environments. The current fundamental requirement for DEER of protein immobilization provides additional avenues toward improvement of the technique.

Architecture of the nuclear pore complex coat

Authors: Tobias Stuwe, Ana R. Correia, Daniel H. Lin, Marcin Paduch, Vincent T. Lu, Anthony A. Kossiakoff, and André Hoelz.

Science. 2015 Mar 6;347(6226):1148-52. (Consortium core: Synthetic Antigen Binder (SAB) Generation and Crystallography)

The nuclear pore complex (NPC) constitutes the sole gateway for bidirectional nucleocytoplasmic transport. Despite half a century of structural characterization, the architecture of the NPC remains unknown. In this research report, the authors present the crystal structure of a reconstituted ~400-kilodalton coat nucleoporin complex (CNC) from Saccharomyces cerevisiae at a 7.4 angstrom resolution. The crystal structure revealed a curved Y-shaped architecture and the molecular details of the coat nucleoporin interactions forming the central “triskelion” of the Y. A structural comparison of the yeast CNC with an electron microscopy reconstruction of its human counterpart suggested the evolutionary conservation of the elucidated architecture. Moreover, 32 copies of the CNC crystal structure docked readily into a cryoelectron tomographic reconstruction of the fully assembled human NPC, thereby accounting for ~16 megadalton of its mass.

Multilevel Summation Method for Electrostatic Force Evaluation

Authors: David J. Hardy, Zhe Wu, James C. Phillips, John E. Stone, Robert D. Skeel, and Klaus Schulten.

J Chem Theory Comput. 2015 Feb 10;11(2):766-779. PMCID: PMC4325600. (Consortium core: Computational Modeling Core)

Significant long-range electrostatic interactions arise in many biomolecular systems, such as negatively charged DNA and RNA, polar or charged membranes, ion channels, and electrostatic steering of protein−protein and enzyme−substrate association. Accordingly, electrostatic interactions need to be accurately represented in molecular modeling calculations. The computational cost increases in principle as N2, where N is the number of charged partiCl-es in the system.

The multilevel summation method (MSM) offers an efficient algorithm utilizing convolution for evaluating long-range forces arising in molecular dynamics simulations. Shifting the balance of computation and communication, MSM provides key advantages over the ubiquitous partiCl-e−mesh Ewald (PME) method, offering better scaling on parallel computers and permitting more modeling flexibility, with support for periodic systems as does PME but also for semiperiodic and nonperiodic systems. The version of MSM available in the simulation program NAMD is described, and its performance and accuracy are compared with the PME method. The accuracy feasible for MSM in practical applications reproduces PME results for water property calculations of density, diffusion constant, dielectric constant, surface tension, radial distribution function, and distance dependent Kirkwood factor, even though the numerical accuracy of PME is higher than that of MSM. Excellent agreement between MSM and PME is found also for interface potentials of air−water and membrane−water interfaces, where long-range Coulombic interactions are crucial. Applications demonstrate also the suitability of MSM for systems with semiperiodic and nonperiodic boundaries. For this purpose, simulations have been performed with periodic boundaries along directions parallel to a membrane surface but not along the surface normal, yielding membrane pore formation induced by an imbalance of charge across the membrane. Using a similar semiperiodic boundary condition, ion conduction through a graphene nanopore driven by an ion gradient has been simulated. Furthermore, proteins have been simulated inside a single spherical water droplet. Finally, parallel scalability results show the ability of MSM to outperform PME when scaling a system of modest size (less than 100 K atoms) to over a thousand processors, demonstrating the suitability of MSM for large-scale parallel simulation.

Ongoing is the development of improved interpolation for MSM to provide higher accuracy for a given polynomial degree p without increasing the computational cost. Future work inCl-udes also the calculation of dispersion forces without truncation with MSM-based NAMD; these forces, in particular, their long-range contribution, are considered to be important for membrane properties. With support in NAMD also for long-range dispersion forces, the present CHARMM-prescribed 12 Å cutoff/splitting distance can be used as a true control for MSM accuracy. High performance simulations will then be able to achieve practical accuracy with a reduced splitting distance, where a splitting distance of between 8 and 9 Å is expected to double the overall simulation performance.

Potential Application of Alchemical Free Energy Simulations to Discriminate GPCR Ligand Efficacy

Authors: Hui Sun Lee, Chaok Seok, and Wonpil Im.

J. Chem. Theory Comput. 2015 Feb 10;11:1255-1266. (Consortium core: Computational Modeling Core)

G protein-coupled receptors (GPCRs) constitute the largest protein superfamily in the human genome with almost 1,000  members. They play key functional roles as major contributors of information flow from the outside to the inside of the cell, making them one of the most important protein families. As a result of their broad influence on human physiology and behavior, GPCRs are arguably the most promising targets for development of new and more effective therapeutic agents.

Based on the fact that GPCR-mediated signaling is modulated in a ligand-specific manner such as agonist, inverse agonist, and neutral antagonist (termed ligand efficacy), quantitative characterization of the ligand efficacy is essential for rational design of selective modulators for GPCR targets. As experimental approaches for this purpose are time-, cost-, and labor-intensive, computational tools that can systematically predict GPCR ligand efficacy can have a big impact on GPCR drug design. Here, the authors have performed free energy perturbation molecular dynamics simulations to calculate absolute binding free energy of an inverse agonist, a neutral antagonist, and an agonist to β2-adrenergic receptor (β2-AR) active and inactive states, respectively, in explicit lipid bilayers. Relatively short alchemical free energy calculations reveal that both the time series of the total binding free energy and decomposed energy contributions can be used as relevant physical properties to discriminate β2-AR ligand efficacy. This study illustrates a merit of the current approach over simple, fast docking calculations or highly expensive millisecond-time scale simulations.

It is the authors’ hope that their computational approach improves research and development efficiency in designing novel lead compounds targeting various GPCRs for the treatment of various human diseases.

Recent Consortium advances in the study of neurotransmitter transporters point the way to future research initiatives

Cartoon model of LeuT transport derived from EPR data, click to enlarge (Publication source)

Using the recently determined crystal structure of a prokaryotic leucine transporter (LeuT), the Transport Cycle in Neurotransmitter Uptake Project explores conformational changes and dynamic properties relevant to function in Neurotransmitter transporters translocation cycle using a combination of computational, functional, and spectroscopic approaches.

The conceptual design, scope and integration of this Project exemplifies the consortium approach to discovery of mechanistic principles of secondary active transport, including the conformational dynamics that govern alternating access in transporters and the allosteric interplay of substrate sites that transduce the energy stored in the electrochemical Na+ gradient into transport work. One of the most recent collaborative publications from members of this Project (lab of Harel Weinstein) and the Consortium at large (labs of Consortium member Olga Boudker and associate member Scott Blanchard) describes such conformational dynamics for a complex transporter system, illustrating the power of the combined experimental / computational approaches we have implemented for the study of these membrane proteins in the possible next phase of the grant:

Akyuz N, Georgieva ER, Zhou Z, Stolzenberg S, Cuendet MA, Khelashvili G, Altman RB, Terry DS, Freed JH, Weinstein H, Boudker O, Blanchard SC.
Transport domain unlocking sets the uptake rate of an aspartate transporter
Nature. 2015 Feb 5;518(7537):68-73. PMID: 25652997. PMCID: PMC4351760.

Crystal structure of the H276,395-GltPh.

Publication abstract: Glutamate transporters terminate neurotransmission by clearing synaptically released glutamate from the extracellular space, allowing repeated rounds of signalling and preventing glutamate-mediated excitotoxicity. Crystallographic studies of a glutamate transporter homologue from the archaeon Pyrococcus horikoshii, GltPh, showed that distinct transport domains translocate substrates into the cytoplasm by moving across the membrane within a central trimerization scaffold. Here we report direct observations of these ‘elevator-like’ transport domain motions in the context of reconstituted proteoliposomes and physiological ion gradients using single-molecule fluorescence resonance energy transfer (smFRET) imaging. We show that GltPh bearing two mutations introduced to impart characteristics of the human transporter exhibits markedly increased transport domain dynamics, which parallels an increased rate of substrate transport, thereby establishing a direct temporal relationship between transport domain motion and substrate uptake. Crystallographic and computational investigations corroborated these findings by revealing that the ‘humanizing’ mutations favour structurally ‘unlocked’ intermediate states in the transport cycle exhibiting increased solvent occupancy at the interface between the transport domain and the trimeric scaffold.

In a second publication from this Project, Matthias Quick (who is working with Project PI Jonathan Javitch), and Computational Modeling Core participant Lei Shi (from Harel Weinstein’s group), published findings that reveal the existence of two substrate sites in vSGLT, a member of the solute-sodium symporter (SSS) family of secondary transporters, as well as in PutP, another member of the SSS family. This publication follows a number of Consortium-sponsored interactions between Quick and Shi as well as Javitch and Weinstein. An earlier 2013 research initiative examined the aggregation dynamics and detergent binding to the second substrate site binding site of LeuT, and in the same year, their team published a key paper on the chloride binding site of neurotransmitter sodium transporters. The very first paper published in the name of the Project, was also the product of collaboration between this team of scientists, and showed how experimental conditions can block the second high-affinity site of LeuT. The important discovery of the present publication suggests that our findings in LeuT are generalizable to a different family of secondary transporters. These findings are the basis for forthcoming scientific collaborations between Quick and Shi, further demonstrating the broadening of the impact of the consortium and the synergy that the “glue” of the consortium has engendered.

Li Z, Lee AS, Bracher S, Jung H, Paz A, Kumar JP, Abramson J, Quick M, Shi L.
Identification of a second substrate-binding site in solute-sodium symporters
J Biol Chem., 2015 Jan 2;290(1):127-41. PMID: 25398883. PMCID: PMC4281715.

Crystallographically identified substrate-binding site of vSGLT is located more extracellularly than that of LeuT.

Publication abstract: The structure of the sodium/galactose transporter (vSGLT), a solute-sodium symporter (SSS) from Vibrio parahaemolyticus, shares a common structural fold with LeuT of the neurotransmitter-sodium symporter family. Structural alignments between LeuT and vSGLT reveal that the crystallographically identified galactose-binding site in vSGLT is located in a more extracellular location relative to the central substrate-binding site (S1) in LeuT. Our computational analyses suggest the existence of an additional galactose-binding site in vSGLT that aligns to the S1 site of LeuT. Radiolabeled galactose saturation binding experiments indicate that, like LeuT, vSGLT can simultaneously bind two substrate molecules under equilibrium conditions. Mutating key residues in the individual substrate-binding sites reduced the molar substrate-to-protein binding stoichiometry to ∼1. In addition, the related and more experimentally tractable SSS member PutP (the Na+/proline transporter) also exhibits a binding stoichiometry of 2. Targeting residues in the proposed sites with mutations results in the reduction of the binding stoichiometry and is accompanied by severely impaired translocation of proline. Our data suggest that substrate transport by SSS members requires both substrate-binding sites, thereby implying that SSSs and neurotransmitter-sodium symporters share common mechanistic elements in substrate transport.

CHARMM-GUI Membrane Builder publication by Im et al. featured as cover image for Journal of Computational Chemistry

Dr. Wonpil Im, University of Kansas

One of the ongoing projects in the laboratory of Wonpil Im (University of Kansas and MPSDC team member) has been to contribute to the development of CHARMM-GUI and CHARMM-GUI modules, such as the Ligand Binder and the Membrane Builder. These tools help users generate a series of CHARMM inputs towards a specific purpose, like calculating a protein/ligand binding free energy or building a protein/membrane complex for molecular dynamics simulations. The Consortium’s Computational Modeling Core has provided partial support for the completion of such modules.

Im and his team recently published a paper about the Membrane Builder module in Journal of Computational Chemistry, titled “CHARMM-GUI Membrane Builder toward realistic biological membrane simulations” (reference below). This publication was featured as the cover image for the issue, which can be viewed below. J Comput Chem attached the following note about the article’s relevance:

On page 1997 (DOI: 10.1002/jcc.23712), Emilia L. Wu, Xi Cheng, Sunhwan Jo, Huan Rui, Kevin C. Song, Eder M. Dávila‐Contreras, Yifei Qi, Jumin Lee, Viviana Monje‐Galvan, Richard M. Venable, Jeffery B. Klauda, and Wonpil Im report on the CHARMM‐GUI Membrane Builder, a web‐based user interface designed to interactively build all‐atom protein/membrane or membrane‐only systems for molecular dynamics simulation through an automated optimized process. The cover image illustrates the building process of a protein/membrane system by reading the protein structure, generating pore water/bilayer/bulk water/ions, and combining all the components.

(Click to enlarge)


Wu EL, Cheng X, Jo S, Rui H, Song KC, Dávila-Contreras EM, Qi Y, Lee J, Monje-Galvan V, Venable RM, Klauda JB, Im W. CHARMM-GUI Membrane Builder toward realistic biological membrane simulations. J Comput Chem. 2014 Oct 15;35(27):1997-2004. PMID: 25130509. PMCID: PMC4165794 [Available on 2015/10/15].

Klaus Schulten’s keynote lecture movie “Photosynthetic Membrane of Purple Bacteria – A Clockwork of Proteins and Processes” made available with audio

Dr. Klaus Schulten

In the coming days, we will be releasing media footage from the exciting Frontiers in Membrane Protein Structural Dynamics 2014 meeting held at the Hilton Hotel in Chicago from May 7th-9th.

As a first piece of footage, we can think of no better than keynote speaker Klaus Schulten (UIUC)’s fascinating atom-by-atom movie titled “Photosynthetic Membrane of Purple Bacteria – A Clockwork of Proteins and Processes“, now made available with Schulten’s narrative of the movie during the keynote lecture.

We initially posted this as a silent movie to accompany our interview with professor Schulten in which he addresses his scientific research interests both past and present, his perspective on some of the key challenges for the field of membrane protein biophysics in the coming 5-10 years, his keynote lecture, and the Membrane Protein Structural Dynamics Consortium.

Now you can watch the movie together with Schulten’s own narration:

Collaborative LeuT study provides new view of neurotransmitter transporters

Model of LeuT alternating access inferred from the crystal structures.Model of LeuT alternating access inferred from the crystal structures.

This week, the Transport Cycle in Neurotransmitter Uptake Systems bridging project of the Membrane Protein Structural Dynamics Consortium (MPSDC) published an important article in Nature Structural & Molecular Biology on the bacterial leucine transporter (LeuT), a transporter which is structurally and functionally similar to neurotransmitter transporter proteins that direct neurotransmitters from synapse and terminal nerve signaling. The publication, titled “Conformational dynamics of ligand-dependent alternating access in LeuT,” was spearheaded by Vanderbilt graduate student Kelli Kazmier and Professor of Molecular Physiciology & Biophysics Hassane Mchaourab, and also featured collaboration by Consortium colleagues Jonathan Javitch, Harel Weinstein, and Benoît Roux.

The Transport Cycle in Neurotransmitter Uptake Systems project explores the conformational changes and dynamic properties relevant to function in Neurotransmitter transporters translocation cycle using a combination of computational, functional, and spectroscopic approaches. Using the recently determined crystal structure of a prokaryotic leucine transporter (LeuT), the scientists collaborating in this project are modeling the transport mechanisms of these proteins.

In this study, Mchaourab and colleagues used spectroscopic tools to make dynamic measurements in LeuT, in order to elucidate sodium- and leucine-dependent conformational This work highlights the importance of assessing the mechanistic identity of crystal structures, demonstrates the importance of dynamics in understanding function and realizes the vision of the consortium in integrating teams of scientists towards defining mechanistic principles of membrane proteins.changes in the transporter. The results identify the structural motifs that underlie the shift of LeuT between its various states – outward-facing, inward-facing and occluded. The conformational changes reported present a dynamic picture of the alternating-access mechanism of LeuT and NSSs that is different from the inferences reached from currently available structural models.

The publication marks a significant advance for the project’s research objectives, and is demonstrative of the cutting-edge collaborations between experimentalists and computationalists within the Consortium. According to Mchaourab, this work “highlights the importance of assessing the mechanistic identity of crystal structures, demonstrates the importance of dynamics in understanding function and realizes the vision of the consortium in integrating teams of scientists towards defining mechanistic principles of membrane proteins.”

The publication was also featured in Research News @ Vanderbilt. Click to read »

Interview with keynote speaker Klaus Schulten

ABC Transporter paper by Tajkhorshid laboratory recommended on F1000Prime

One of our long-term projects on Structural Dynamics of ABC Transporters integrates computational, biochemical and spectroscopic approaches to understand the structural dynamics of the ATP binding cassette (ABC) transporters, which are associated with a number of human pathologies and play critical roles in the removal of cytotoxic agents.

Among the MPSDC participants in this project is Dr. Emad Tajkhorshid, whose contribution is applying molecular dynamics (MD) simulations integrating experimental constraints to develop structural models for key conformational states and characterize their inter-conversion during the transport cycle.

Tajkhorshid, together with postdoctoral researcher Mahmoud Moradi, recently published a paper on conformational transitions of ATP exporters which was recommended on the F1000Prime or “Faculty of 1000″ website. F1000 is a team of 5,000 Faculty Members – senior scientists and leading experts in all areas of biology and medicine — plus their associates who provide recommendations of important scientific articles, rating them and providing short explanations for their selections.

The publication by Moradi and Tajkhorshid, titled “Mechanistic picture for conformational transition of a membrane transporter at atomic resolution“, was published on November 19, 2013 in PNAS vol. 110, no. 47. The paper describes a nonequilibrium approach which they developed to characterize the conformational transition of MsbA, a member of the ATP-binding cassette exporter family, which is involved in transport of diverse substrates across the membrane. The F1000Prime recommendation, written by Qian Cui, tagged the publication as Good for Teaching, as having an Interesting Hypothesis, and for Technical Advance. Read the F1000 recommendation »

Dancing Proteins: Cell Membrane Transporter Motion May Revolutionize Drug Therapies (video)

The Beckman Institute at the University of Illinois at Urbana-Champaign produced the following video. In the video, Dr. Tajkhorshid describes how his laboratory has successfully simulated the molecular dance moves that a multidrug resistance membrane transporter undertakes as it pumps compounds out of a cell. This is the first time researchers have been able to simulate the motion of a complex membrane transporter in its native environment in full atomic detail and gives drug developers vital new targets to help combat drug-resistant cancers and other diseases.

Other plaudits

The popular NIH Biomedical Beat blog, which covers research news from NIGMS, featured the video on their website. As per the blog page, “In this video, Emad Tajkhorshid of the University of Illinois at Urbana-Champaign explains the molecular dance of ABC transporters, a family of molecular machines that utilize ATP to move substances across the cell membrane. Tajkhorshid and his team recently used computational methods to map the movements between two known structural models of MsbA, a bacterial version of a transporter in human cells that helps to export anti-cancer drugs. They then described the individual steps of the molecular motions during the transport cycle. Understanding the process at such a detailed level could suggest new targets for treating a range of diseases, including some drug-resistant cancers that often make more transporter proteins to kick out medications meant to kill them.”

Additionally, Tajkhorshid and Moradi were also featured in the University of Illinois News Bureau, in an article titled “Difficult dance steps: Team learns how membrane transporter moves.” The article helpfully describes the nature of the research and points to its innovation.

Photo taken from the University of Illinois News Bureau website. Photo by L. Brian Stauffer.

According to the article, “the new findings, reported in the Proceedings of the National Academy of Sciences, will help scientists figure out how other transporters work. The work also offers new insights into multi-drug-resistant (MDR) cancers, some of which use these transporters to export cancer-killing drugs.” Previously, it has been difficult to research large, membrane-bound proteins like MsbA because they are not easy to crystallize, and each crystal structure reflects only one of the many conformations of these shape-shifting proteins. This study marks “the first time that we are characterizing a very complex structural transition at atomic-level resolution for a large protein,” Dr. Tajkhorshid is quoted as saying.

Emad Tajkhorshid is Professor of Biochemistry, Biophysics, and Pharmacology and an affiliate of the Beckman Institute for Advanced Science and Technology at the University of Illinois at Urbana-Champaign.

CLC Channel/Transporter Family Project team combines computational and experimental methodologies in PNAS publication

The Conformational Dynamics in the CLC Channel/Transporter Family Project of the Membrane Protein Structural Dynamics Consortium (MPSDC) has published its first publication titled “Water access points and hydration pathways in CLC H+/Cl- transporters” in Proceedings of the National Academy of Sciences of the United States of America (PNAS). This Consortium Project is spearheaded by Principal Investigator Merritt Maduke. The laboratories of Marc Baldus, Emad Tajkhorshid, and Hassane Mchaourab are also collaborators in the Project’s ongoing research.

Water access points and hydration pathways in CLC H+/Cl- transporters

Figure 3. Entryways of water into the central hydrophobic region. (Han et al. 2013)

Abstract: CLC transporters catalyze transmembrane exchange of chloride for protons. Although a putative pathway for Cl has been established, the pathway of H+ translocation remains obscure. Through a highly concerted computational and experimental approach, we characterize microscopic details essential to understanding H+-translocation. An extended (0.4 µs) equilibrium molecular dynamics simulation of membrane-embedded, dimeric ClC-ec1, a CLC from Escherichia coli, reveals transient but frequent hydration of the central hydrophobic region by water molecules from the intracellular bulk phase via the interface between the two subunits. We characterize a portal region lined by E202, E203, and A404 as the main gateway for hydration. Supporting this mechanism, site-specific mutagenesis experiments show that ClC-ec1 ion transport rates decrease as the size of the portal residue at position 404 is increased. Beyond the portal, water wires form spontaneously and repeatedly to span the 15-Å hydrophobic region between the two known H+ transport sites [E148 (Gluex) and E203 (Gluin)]. Our finding that the formation of these water wires requires the presence of Cl explains the previously mystifying fact that Cl occupancy correlates with the ability to transport protons. To further validate the idea that these water wires are central to the H+ transport mechanism, we identified I109 as the residue that exhibits the greatest conformational coupling to water wire formation and experimentally tested the effects of mutating this residue. The results, by providing a detailed microscopic view of the dynamics of water wire formation and confirming the involvement of specific protein residues, offer a mechanism for the coupled transport of H+ and Cl ions in CLC transporters.

The results were published in Proceedings of the National Academy of Sciences of the United States of America, and are currently made available as an e-publication ahead of print. The citation is as follows:

Wei Han, Ricky C. Cheng, Merritt C. Maduke, and Emad Tajkhorshid.
Water access points and hydration pathways in CLC H+/Cl− transporters
PNAS 2013; published ahead of print December 30, 2013, doi:10.1073/pnas.1317890111

Learn more about this publication »


CLC transporters are biologically essential proteins that catalyze the transmembrane exchange of chloride for protons. The permeation pathway for chloride through the transporters has been well characterized. In this publication, Han et al. study the more elusive permeation pathway for protons. Through computational modeling, they show that water molecules can permeate deep inside the protein and form continuous wires. To test the hypothesis that these water wires mediate proton transport, they mutated residues predicted to impede water wire This research article reports results from the tightly coordinated efforts of a computational and an experimental lab brought together by the Consortium. The study addresses a critical question about the CLC transporter mechanism: how does H+ traverse the hydrophobic expanse of the CLC protein?formation and experimentally evaluated the effects of the mutations. The results from their concerted computational and experimental approach strongly support the role of water in proton transport by CLCs and provide a valuable framework for investigating their overall transport mechanism.

In a commentary piece published by PNAS this month, Mounir Tarek (National Center of Scientific Research at the University of Lorraine, France) describes the significance of this paper for future research of chloride channels, and highlights the fruitful combination of simulation and experimentation: “In PNAS, Han et al. used molecular dynamics (MD) simulations of the CLC-ecl, a CLC exchanger from Escherichia coli to specifically address this issue. The predictions of their calculations were tested by additional experiments, providing a robust description of the molecular prerequisites to proton transport in CLC-ecl and a framework for refining models of the Cl-/H+-coupled transport in CLCs.”

Indeed, this research article reports results from the tightly coordinated efforts of a computational and an experimental lab brought together by the Consortium. The study addresses a critical question about the CLC transporter mechanism: how does H+ traverse the hydrophobic expanse of the CLC protein? The MD simulations performed reveal water dynamics, water-wire formation, and side-chain conformational change not observed in any of the static crystal structures. The functional analyses validated predictions of the simulations and confirm the importance of water dynamics in the transport mechanism. The simulations further reveal that Cl- binding is critical for water-wire formation, thus providing a satisfying explanation for the puzzling experimental observation that Cl- occupancy correlates with the ability of CLCs to transport H+. These studies provide a crucial framework for understanding how H+ and Cl- binding and translocation steps are coordinated in the CLC transporters to control stoichiometric transport.

About the project

The CLC family of chloride channels and transporters is necessary for proper neuronal, cardiovascular, and epithelial function. One of the important aspects of this family of transport proteins is that minute changes in their amino acid sequence can result in a shift in their operation mode from a channel to a transporter. Studying the structural dynamics of CLCs can therefore provide fundamental information on the nature of structural and dynamical differences between passive channels and active transporters.

The Conformational Dynamics in the CLC Channel/Transporter Family project addresses the multiple structural conformations that underlie the dual function of ClC- proteins as both channels and coupled transporters. Using a combination of NMR (solution and solid-state) and molecular dynamics simulations, the multiple conformations that support closely-coupled, stoichiometric ion transport will be accessed by binding and unbinding its two ligands, (Cl- and H+). Additional efforts are made made to use conformation-specific ligands to “lock” CLC proteins in order to study these conformations by crystallography, EPR, and NMR.

Learn more about the project »

Overview: Structural Refinement Based on EPR Data from Restrained-Ensemble Simulation

Shahidul M. Islam,† Richard A. Stein,‡ Hassane S. Mchaourab,‡ and Benoît Roux

† Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois
‡ Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, Tennessee

Membrane proteins perform numerous physiological functions that are critical for human health. They account for 60% of drug targets and mutations in their primary sequence hamper their normal function, which can lead to various diseases. Membrane proteins act like “molecular machines”, changing their shape and visiting many conformational states to perform their function. Knowledge of all the important states is critical to understand these proteins; however, the process of obtaining such structural information is not at all straightforward. While X-ray crystallography is commonly the best technique to obtain high-resolution structural information of proteins, it is important to achieve a more complete picture of the accessible conformational states of a protein in its native environment, free from the constraints of the crystal lattice. Moreover, membrane proteins are underrepresented in the protein structure data bank.

Two important spectroscopic approaches occupy a central role in the efforts to understand structure and function of membrane proteins: nuclear magnetic resonance (NMR) and electron spin resonance (ESR). However, NMR investigations are limited by the size of the protein system, while ESR requires the introduction of spectroscopic probes into the system via site-directed spin-labeling (SDSL) techniques. Recently, a novel computational simulation technique was developed to exploit the information from distance distribution data obtained from ESR/DEER spectroscopy for the refinement of membrane protein structures. This simulation technique, called the Restrained-Ensemble Molecular Dynamics (REMD) simulation method, uses a global ensemble-based energy restraint to force the spin-spin distance distribution histograms calculated from a multiple-copy molecular dynamics simulation to match those obtained from ESR/DEER experiments.An important advantage of ESR technique is that it provides strong signal from the spin labels in the case of extremely large macromolecular complexes, such as in the case of large membrane proteins. Of particular interest, DEER (Double Electron-Electron Resonance) is a powerful ESR technique that makes it possible to measure the distance histogram between a pair of spin-labels inserted in a macromolecule.

The most commonly used nitroxide spin-label is MTSSL (1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl-methanethiosulfonate), which is typically linked to a cysteine residue in the protein through a disulphide bond (Figure 1A). The MTSSL moiety has five dihedral bonds, resulting in a highly flexible side-chain once linked to a protein. Because of this flexibility, spin label incorporation does not lead to significant changes of the original protein conformation; however, it introduces complexity in “translating” the ESR data on the spin labels to obtain structural information on the protein itself. A reliable characterization of the dynamical properties of the spin-label is therefore essential to fully exploit the spin-pair distance histograms for the purpose of structural refinement of membrane proteins. Computational modeling methods, such as the Multiscale Modeling of Macromolecular systems (MMM) software package of Yevhen Polyhach and Gunnar Jeschke1,2, and the PRONOX algorithm of Hatmal et al,3 have been developed to determine the inter-label distances distributions based on the analysis of spin-label rotamers. In spite of these efforts, there remains a need to develop a robust and effective computational method for making best use of ESR/DEER data in the context of structural refinement. All the previous computational simulation studies4-7 and modeling methods1,3 use the ESR/DEER distance histogram in post-analysis, to assess the correctness of models that were generated independently from the experimental data. In other words, none of the existing methods “drive” the structural model toward a 3D conformation that satisfies the ESR/DEER data.

(click to enlarge)

Figure 1. (A) Cartoon representation of T4 lysozyme (T4L) with 37 spin-labeled sites (colored in blue). As an Enrolling in Medicare Hospital (Medicare Part A) is a no brainer. example, the structural formula of two spin label side chains at positions 62 and 109 in T4L is shown here. The Restrained-Ensemble (RE) simulation method has been developed to restrain the calculated spin-spin distance distributions obtain from MD simulation to match with those obtained from the ESR/DEER, (B) From RE simulation, the dynamics of nitroxide oxygen of spin-labels around the Cα atom was obtained, (C) A simplified dummy nitroxide spin-label was designed and parameterized, (D) The simplified dummy spin-labels along with RE simulation refined distorted structures of T4L which is illustrated with root mean square deviation of the backbone atoms with respect to the X-ray structure (colored in blue). NOE-like distance restraints were unable to refine many of the distorted structures (colored in black).

Recently,8-10 a novel computational simulation technique was developed to exploit the information from distance distribution data obtained from ESR/DEER spectroscopy for the refinement of membrane protein structures. This simulation technique, called the Restrained-Ensemble Molecular Dynamics (REMD) simulation method, uses a global ensemble-based energy restraint to force the spin-spin distance distribution histograms calculated from a multiple-copy molecular dynamics simulation to match those obtained from ESR/DEER experiments. The REMD simulation method was applied to 51 ESR/DEER distance histogram data from spin-labels inserted at 37 different positions in a membrane protein, as well as the T4 lysozyme (T4L) (as shown in Figure 1). The T4 lysozyme is typically used as a model system by the ESR community due to the availability of extensive amount of structural data from X-ray crystallography, NMR and multifrequency ESR experiments. Millions of data were collected and stored in a rotameric library for the spin label side chains at various positions in T4L. The rotamer population distributions are shown to be consistent with available information from X-ray crystallography. From the all atom RE simulations, the authors designed and parameterized a simplified nitroxide dummy spin-label (Figure 1B and 1C), which was finally used for the purpose of structural refinement. The authors finally demonstrated that RE simulations with the dummy nitroxide spin-labels imposing the ESR/DEER experimental distance distribution data are able to systematically correct and refine a series of distorted T4L structures (Figure 1D). This computationally efficient approach allows experimental restraints from DEER experiments to be incorporated into RE simulations for efficient structural refinement. It is expected that this novel simulation method will revolutionize our perspective on the refinement of membrane proteins. The authors are currently investigating several important membrane proteins such as the ion channels and transmembrane proteins with this Restrained-Ensemble simulation method.


(1) Polyhach, Y.; Bordignon, E.; Jeschke, G. Rotamer libraries of spin labelled cysteines for protein studies, Phys. Chem. Chem. Phys., 2011, 13, 2356.

(2) Jeschke, G. DEER distance measurements on proteins, Annu Rev Phys Chem, 2012, 63, 419.

(3) Hatmal, M. M.; Li, Y.; Hegde, B. G.; Hegde, P. B.; Jao, C. C.; Langen, R.; Haworth, I. S. Computer modeling of nitroxide spin labels on proteins, Biopolymers, 2012, 97, 35.

(4) Polyhach, Y.; Godt, A.; Bauer, C.; Jeschke, G. Spin pair geometry revealed by high-field DEER in the presence of conformational distributions, J. Magn. Reson., 2007, 185, 118.

(5) Tikhonova, I. G.; Best, R. B.; Engel, S.; Gershengorn, M. C.; Hummer, G.; Costanzi, S. Atomistic insights into rhodopsin activation from a dynamic model, J. Am. Chem. Soc., 2008, 130, 10141.

(6) Ding, F.; Layten, M.; Simmerling, C. Solution structure of HIV-1 protease flaps probed by comparison of molecular dynamics simulation ensembles and EPR experiments, Journal of the American Chemical Society, 2008, 130, 7184.

(7) Boura, E.; Rozycki, B.; Herrick, D. Z.; Chung, H. S.; Vecer, J.; Eaton, W. A.; Cafiso, D. S.; Hummer, G.; Hurley, J. H. Solution structure of the ESCRT-I complex by small-angle X-ray scattering, EPR, and FRET spectroscopy, Proc. Natl. Acad. Sci. USA, 2011, 108, 9437.

(8) Roux, B.; Islam, S. M. Restrained-Ensemble Molecular Dynamics Simulations Based on Histograms from Double Electron-Electron Resonance Spectroscopy, J. Phys. Chem. B, 2013, In Press. (link)

(9) Roux, B.; Weare, J. On the statistical equivalence of restrained-ensemble simulations with the maximum entropy method, J. Chem. Phys., 2013, 138. (link)

(10) Islam, S. M.; Stein, R.; Mchaourab, H.; Roux, B. Structural Refinement from Restrained-Ensemble Simulations Based on EPR/DEER Data: Application to T4 Lysozyme, J. Phys. Chem. B, 2013, In Press. (link)

Flexibility of the Flexizyme

By Robert K. Nakamoto
Department of Molecular Physiology and Biological Physics, University of Virginia

There is almost nothing more useful to the biophysicist (except perhaps a talented grad student or post-doc) than to have probes in known positions within the structure of a protein. Spectroscopic and biochemical approaches use such probes to monitor the dynamics, environment and conformational shifts that occur during function. Probe placement is most commonly done by site-directed mutagenesis to introduce a unique cysteine and chemically modifying the purified protein with a sulfhydryl-reactive label. In recent years, investigators have developed new methods to incorporate custom-made amino acids directly into polypeptide chains by expanding the genetic code. In principle, amino acids with side chains of any chemical design can be incorporated into the growing polypeptide chain by the ribosome if the acyl-tRNA charged with the unnatural amino acid can be made. Usually a tRNA is used that recognizes one of the non-sense codons (UAG, AAG or UGA) or a codon not otherwise used in the coding sequence of the open reading frame. Of course, there are many limitations to this process but the main one is acylation of the 3’ base of the tRNA with the amino acid. This reaction is done in the cell by the amino acid-specific aminoacyl-tRNA synthetase. To broaden the genetic code, investigators have developed methods to “evolve” orthogonal synthetases to shift the specificities to their custom amino acid of choice (see Liu & Shultz for a review). Here, I describe a new and highly flexible method for synthesizing the acyl-tRNA with a wide range of unnatural amino acids using the “Flexizyme”.

Developed and optimized in the laboratory of Professor Hiroaki Suga at the University of Tokyo, the “Flexizyme” approach is based on the use of a ribozyme to incorporate unnatural amino acids in proteins. A Nature Protocols paper (Goto et al.) describes the entire system and provides a detailed list of reagents and instructions. Though highly flexible in many ways, it is important to highlight that the system only works in cell-free translation systems.

(click to enlarge)

From Flexizymes for genetic code reprogramming
Yuki Goto, Takayuki Katoh & Hiroaki Suga
Nature Protocols 6, 779–790 (2011)

Like in most systems designed to incorporate unnatural amino acids (see Liu & Shultz) a charged suppressor tRNA is used to integrate a designed amino acid into the nascent polypeptide chain. However, rather than acylating the tRNAsup with an aminoacyl-tRNA synthetase specifically evolved to accept only a given unnatural amino acid, the Flexizyme system acylates the tRNAsup with a ribozyme. This avoids altogether the difficult and somewhat daunting process of evolving an orthogonal tRNA synthetase to recognize each unnatural amino acid. The Suga laboratory has developed a set of Flexizymes capable of charging a tRNA with a diverse array of unnatural amino acids. Because the system is used only in vitro, it also avoids the additional worry of importing the unnatural amino acids into a cell, which can have toxic effects.

A key part of the system is the acyl-donor substrate, which is used by the Flexizyme to acylate the tRNAsup. Because the reaction is primarily driven by the nature of the leaving group, the amino acid part of the substrate can vary. Suga and his colleagues have developed a set of three Flexizymes optimized to work with four different leaving groups. The Flexizyme system is particularly attractive for incorporating spectroscopic probes such as fluorescence or nitroxide spin labels at user-defined positions. Normally, investigators try to create a background devoid of cysteines and unique sulfhydryls are introduced to provide specific reactive positions. The Flexizyme system abrogates the potential problems that replacing native Cys can have on protein stability and function, or that a given Cys may be inaccessible and unreactive.The investigator must determine which donor-leaving group/Flexizyme combination works best with their specific unnatural amino acid. For example, aromatic side chains generally work best with the cyanomethyl ester-leaving group and the “eFx” ribozyme, while non-aromatic side chains would use the 3,5-dinitrobenzyl ester with the “dFx” ribozyme (see Goto et al. for sequences and structures). The Flexizymes and tRNAs are easily made in the laboratory and the chemistries to synthesize the acyl-donor substrates are relatively straightforward. Assays to determine the efficiencies of the transfer reactions are described by Murakami et al.

In a very nice example, Öjemalm et al. used the Flexizyme system to create substrates to study the influence of polypeptide hydrophobicity on insertion into the endoplasmic reticulum membrane via the Sec61 translocon. The Flexizyme was used to charge the tRNAsup with a series of unnatural amino acids with linear aliphatic side chains up to 10 carbons in length.

The Flexizyme system is particularly attractive for incorporating spectroscopic probes such as fluorescence or nitroxide spin labels at user-defined positions. Normally, investigators try to create a background devoid of cysteines and unique sulfhydryls are introduced to provide specific reactive positions. The Flexizyme system abrogates the potential problems that replacing native Cys can have on protein stability and function, or that a given Cys may be inaccessible and unreactive. Furthermore, there is a wider choice of labels with the Flexizyme and the investigator can use the proper probe to monitor specific properties such as backbone dynamics or local environments, or to impart desired modifications into the structure.

In combination with robust cell free synthesis systems such as those developed by Volker Dõtsch and colleagues (see Klammt et al. for a review, and the Membrane Protein Expression and Purification Core page of this website for detailed protocols), the Flexizyme system can be a powerful tool for the efficient production of the modified protein in biochemical amounts. Obviously, the system is potentially limited by the ease with which the in vitro synthesized polypeptide can be folded into its native conformation. This can be particularly difficult with integral membrane proteins, since folding must be done in detergents, lipids or a complex combination of both. Recent successes are indeed encouraging and the advantages provided by combining the two approaches represent an untapped resource in the analysis of the structure and dynamics of membrane proteins.


Goto, Y., Katoh, T. and Suga, H. (2011) Flexizymes for genetic code reprogramming. Nature Protoc. 6, 779-790.

Klammt, C., Löhr, F., Schäfer, B., Haase, W., Dötsch, V., Rüterjans, H., Glaubitz, C. and Bernhard, F. (2004) High level cell-free expression and specific labeling of integral membrane proteins. Eur. J. Biochem. 271, 568–580.

Liu, C. C. and Schultz, P. G. (2010) Adding new chemistries to the genetic code. Annu. Rev. Biochem. 79, 413-444.

Murakami, H., Ohta, A., Ashigai, H. and Suga, H. (2006) A highly flexible tRNA acylation method for non-natural polypeptide synthesis. Nature Meth. 3, 357-359.

Öjemalm, K., Higuchi, T., Jiang, Y., Langel, Ü., Nilsson, I., White, S. H., Suga, H. and von Heijne, G. (2011) Apolar surface area determines the efficiency of translocon-mediated membrane-protein integration into the endoplasmic reticulum. Proc. Natl. Acad. Sci. U.S.A. 108, 359-364.

Recent Posts
Recent Comments
Recent Photos