eBDIMS server: protein transition pathways with ensemble analysis in 2D-motion spaces

CryoEM structures of anion exchanger 1 capture multiple states of inward- and outward-facing conformations

Anion exchanger 1 (AE1, band 3) is a major membrane protein of red blood cells and plays a key role in acid-base homeostasis, urine acidification, red blood cell shape regulation, and removal of carbon dioxide during respiration. Though structures of the transmembrane domain (TMD) of three SLC4 transporters, including AE1, have been resolved previously in their outward-facing (OF) state, no mammalian SLC4 structure has been reported in the inward-facing (IF) conformation. Here we present the cryoEM structures of full-length bovine AE1 with its TMD captured in both IF and OF conformations. Remarkably, both IF-IF homodimers and IF-OF heterodimers were detected. The IF structures feature downward movement in the core domain with significant unexpected elongation of TM11. Molecular modeling and structure guided mutagenesis confirmed the functional significance of residues involved in TM11 elongation. Our data provide direct evidence for an elevator-like mechanism of ion transport by an SLC4 family member.

Exploring the Conformational Impact of Glycine Receptor TM1-2 Mutations Through Coarse-Grained Analysis and Atomistic Simulations

Pentameric ligand-gated ion channels (PLGICs) are a family of proteins that convert chemical signals into ion fluxes through cellular membranes. Their structures are highly conserved across all kingdoms from bacteria to eukaryotes. Beyond their classical roles in neurotransmission and neurological disorders, PLGICs have been recently related to cell proliferation and cancer. Here, we focus on the best characterized eukaryotic channel, the glycine receptor (GlyR), to investigate its mutational patterns in genomic-wide tumor screens and compare them with mutations linked to hyperekplexia (HPX), a Mendelian neuromotor disease that disrupts glycinergic currents. Our analysis highlights that cancer mutations significantly accumulate across TM1 and TM2, partially overlapping with HPX changes. Based on 3D-clustering, conservation, and phenotypic data, we select three mutations near the pore, expected to impact GlyR conformation, for further study by molecular dynamics (MD). Using principal components from experimental GlyR ensembles as framework, we explore the motions involved in transitions from the human closed and desensitized structures and how they are perturbed by mutations. Our MD simulations show that WT GlyR spontaneously explores opening and re-sensitization transitions that are significantly impaired by mutations, resulting in receptors with altered permeability and desensitization properties in agreement with HPX functional data.

Markov state models of proton- and pore-dependent activation in a pentameric ligand-gated ion channel

Ligand-gated ion channels conduct currents in response to chemical stimuli, mediating electrochemical signaling in neurons and other excitable cells. For many channels, the details of gating remain unclear, partly due to limited structural data and simulation timescales. Here, we used enhanced sampling to simulate the pH-gated channel GLIC, and construct Markov state models (MSMs) of gating. Consistent with new functional recordings, we report in oocytes, our analysis revealed differential effects of protonation and mutation on free-energy wells. Clustering of closed- versus open-like states enabled estimation of open probabilities and transition rates, while higher-order clustering affirmed conformational trends in gating. Furthermore, our models uncovered state- and protonation-dependent symmetrization. This demonstrates the applicability of MSMs to map energetic and conformational transitions between ion-channel functional states, and how they reproduce shifts upon activation or mutation, with implications for modeling neuronal function and developing state-selective drugs.

Probing solution structure of the pentameric ligand-gated ion channel GLIC by small-angle neutron scattering

Pentameric ligand-gated ion channels undergo subtle conformational cycling to control electrochemical signal transduction in many kingdoms of life. Several crystal structures have now been reported in this family, but the functional relevance of such models remains unclear. Here, we used small-angle neutron scattering (SANS) to probe ambient solution-phase properties of the pH-gated bacterial ion channel GLIC under resting and activating conditions. Data collection was optimized by inline paused-flow size-exclusion chromatography, and exchanging into deuterated detergent to hide the micelle contribution. Resting-state GLIC was the best-fit crystal structure to SANS curves, with no evidence for divergent mechanisms. Moreover, enhanced-sampling molecular-dynamics simulations enabled differential modeling in resting versus activating conditions, with the latter corresponding to an intermediate ensemble of both the extracellular and transmembrane domains. This work demonstrates state-dependent changes in a pentameric ion channel by SANS, an increasingly accessible method for macromolecular characterization with the coming generation of neutron sources.

ClustENMD: Efficient sampling of biomolecular conformational space at atomic resolution

Summary

Efficient sampling of conformational space is essential for elucidating functional/allosteric mechanisms of proteins and generating ensembles of conformers for docking applications. However, unbiased sampling is still a challenge especially for highly flexible and/or large systems. To address this challenge, we describe a new implementation of our computationally efficient algorithm ClustENMD that is integrated with ProDy and OpenMM softwares. This hybrid method performs iterative cycles of conformer generation using elastic network model for deformations along global modes, followed by clustering and short molecular dynamics simulations. ProDy framework enables full automation and analysis of generated conformers and visualization of their distributions in the essential subspace.

Availability and implementation

ClustENMD is open-source and freely available under MIT License from https://github.com/prody/ProDy.

Supplementary information

Supplementary data are available at Bioinformatics online.

Transcription factor allosteric regulation through substrate coordination to zinc

The development of new synthetic biology circuits for biotechnology and medicine requires deeper mechanistic insight into allosteric transcription factors (aTFs). Here we studied the aTF UxuR, a homodimer of two domains connected by a highly flexible linker region. To explore how ligand binding to UxuR affects protein dynamics we performed molecular dynamics simulations in the free protein, the aTF bound to the inducer D-fructuronate or the structural isomer D-glucuronate. We then validated our results by constructing a sensor plasmid for D-fructuronate in Escherichia coli and performed site-directed mutagenesis. Our results show that zinc coordination is necessary for UxuR function since mutation to alanines prevents expression de-repression by D-fructuronate. Analyzing the different complexes, we found that the disordered linker regions allow the N-terminal domains to display fast and large movements. When the inducer is bound, UxuR can sample an open conformation with a more pronounced negative charge at the surface of the N-terminal DNA binding domains. In opposition, in the free and D-glucuronate bond forms the protein samples closed conformations, with a more positive character at the surface of the DNA binding regions. These molecular insights provide a new basis to harness these systems for biological systems engineering.

HOPMA: Boosting Protein Functional Dynamics with Colored Contact Maps

In light of the recent very rapid progress in protein structure prediction, accessing the multitude of functional protein states is becoming more central than ever before. Indeed, proteins are flexible macromolecules, and they often perform their function by switching between different conformations. However, high-resolution experimental techniques such as X-ray crystallography and cryogenic electron microscopy can catch relatively few protein functional states. Many others are only accessible under physiological conditions in solution. Therefore, there is a pressing need to fill this gap with computational approaches. We present HOPMA, a novel method to predict protein functional states and transitions by using a modified elastic network model. The method exploits patterns in a protein contact map, taking its 3D structure as input, and excludes some disconnected patches from the elastic network. Combined with nonlinear normal mode analysis, this strategy boosts the protein conformational space exploration, especially when the input structure is highly constrained, as we demonstrate on a set of more than 400 transitions. Our results let us envision the discovery of new functional conformations, which were unreachable previously, starting from the experimentally known protein structures. The method is computationally efficient and available at https://github.com/elolaine/HOPMA and https://team.inria.fr/nano-d/software/nolb-normal-modes.

A Coarse-Grained Methodology Identifies Intrinsic Mechanisms That Dissociate Interacting Protein Pairs

We address the problem of triggering dissociation events between proteins that have formed a complex. We have collected a set of 25 non-redundant, functionally diverse protein complexes having high-resolution three-dimensional structures in both the unbound and bound forms. We unify elastic network models with perturbation response scanning (PRS) methodology as an efficient approach for predicting residues that have the propensity to trigger dissociation of an interacting protein pair, using the three-dimensional structures of the bound and unbound proteins as input. PRS reveals that while for a group of protein pairs, residues involved in the conformational shifts are confined to regions with large motions, there are others where they originate from parts of the protein unaffected structurally by binding. Strikingly, only a few of the complexes have interface residues responsible for dissociation. We find two main modes of response: In one mode, remote control of dissociation in which disruption of the electrostatic potential distribution along protein surfaces play the major role; in the alternative mode, mechanical control of dissociation by remote residues prevail. In the former, dissociation is triggered by changes in the local environment of the protein, e.g., pH or ionic strength, while in the latter, specific perturbations arriving at the controlling residues, e.g., via binding to a third interacting partner is required for decomplexation. We resolve the observations by relying on an electromechanical coupling model which reduces to the usual elastic network result in the limit of the lack of coupling. We validate the approach by illustrating the biological significance of top residues selected by PRS on select cases where we show that the residues whose perturbation leads to the observed conformational changes correspond to either functionally important or highly conserved residues in the complex.

Towards gaining sight of multiscale events: utilizing network models and normal modes in hybrid methods

With the explosion of normal mode analyses (NMAs) based on elastic network models (ENMs) in the last decade, and the proven precision of MD simulations for visualizing interactions at atomic scale, many hybrid methods have been proposed in recent years. These aim at exploiting the best of both worlds: the atomic precision of MD that often fall short of exploring time and length scales of biological interest, and the capability of ENM-NMA to predict the cooperative and often functional rearrangements of large structures and assemblies, albeit at low resolution. We present an overview of recent progress in the field with examples of successful applications highlighting the utility of such hybrid methods and pointing to emerging future directions guided by advances in experimental characterization of biomolecular systems structure and dynamics.

The molecular basis for sugar import in malaria parasites

Elucidating the mechanism of sugar import requires a molecular understanding of how transporters couple sugar binding and gating events. Whereas mammalian glucose transporters (GLUTs) are specialists¹, the hexose transporter from the malaria parasite Plasmodium falciparum PfHT1^2,3 has acquired the ability to transport both glucose and fructose sugars as efficiently as the dedicated glucose (GLUT3) and fructose (GLUT5) transporters. Here, to establish the molecular basis of sugar promiscuity in malaria parasites, we determined the crystal structure of PfHT1 in complex with D-glucose at a resolution of 3.6 Å. We found that the sugar-binding site in PfHT1 is very similar to those of the distantly related GLUT3 and GLUT5 structures^4,5. Nevertheless, engineered PfHT1 mutations made to match GLUT sugar-binding sites did not shift sugar preferences. The extracellular substrate-gating helix TM7b in PfHT1 was positioned in a fully occluded conformation, providing a unique glimpse into how sugar binding and gating are coupled. We determined that polar contacts between TM7b and TM1 (located about 15 Å from D-glucose) are just as critical for transport as the residues that directly coordinate D-glucose, which demonstrates a strong allosteric coupling between sugar binding and gating. We conclude that PfHT1 has achieved substrate promiscuity not by modifying its sugar-binding site, but instead by evolving substrate-gating dynamics.

Large-Scale Conformational Changes and Protein Function: Breaking the in silico Barrier

Large-scale conformational changes are essential to link protein structures with their function at the cell and organism scale, but have been elusive both experimentally and computationally. Over the past few years developments in cryo-electron microscopy and crystallography techniques have started to reveal multiple snapshots of increasingly large and flexible systems, deemed impossible only short time ago. As structural information accumulates, theoretical methods become central to understand how different conformers interconvert to mediate biological function. Here we briefly survey current in silico methods to tackle large conformational changes, reviewing recent examples of cross-validation of experiments and computational predictions, which show how the integration of different scale simulations with biological information is already starting to break the barriers between the in silico, in vitro, and in vivo worlds, shedding new light onto complex biological problems inaccessible so far.

Oncogenic mutations at the EGFR ectodomain structurally converge to remove a steric hindrance on a kinase-coupled cryptic epitope

Epidermal growth factor receptor (EGFR) signaling is initiated by a large ligand-favored conformational change of the extracellular domain (ECD) from a closed, self-inhibited tethered monomer, to an open untethered state, which exposes a loop required for strong dimerization and activation. In glioblastomas (GBMs), structurally heterogeneous missense and deletion mutations concentrate at the ECD for unclear reasons. We explore the conformational impact of GBM missense mutations, combining elastic network models (ENMs) with multiple molecular dynamics (MD) trajectories. Our simulations reveal that the main missense class, located at the I-II interface away from the self-inhibitory tether, can unexpectedly favor spontaneous untethering to a compact intermediate state, here validated by small-angle X-ray scattering (SAXS). Significantly, such intermediate is characterized by the rotation of a large ECD fragment (N-TR1), deleted in the most common GBM mutation, EGFRvIII, and that makes accessible a cryptic epitope characteristic of cancer cells. This observation suggested potential structural equivalence of missense and deletion ECD changes in GBMs. Corroborating this hypothesis, our FACS, in vitro, and in vivo data demonstrate that entirely different ECD variants all converge to remove N-TR1 steric hindrance from the 806-epitope, which we show is allosterically coupled to an intermediate kinase and hallmarks increased oncogenicity. Finally, the detected extraintracellular coupling allows for synergistic cotargeting of the intermediate with mAb806 and inhibitors, which is proved herein.