Allosteric Competition and Inhibition in AMPA Receptors.

Excitatory neurotransmission is principally mediated by AMPA-subtype ionotropic glutamate receptors (AMPARs). Dysregulation of AMPARs is the cause of many neurological disorders and how therapeutic candidates such as negative allosteric modulators inhibit AMPARs is unclear. Here, we show that non-competitive inhibition desensitizes AMPARs to activation and prevents positive allosteric modulation. We dissected the noncompetitive inhibition mechanism of action by capturing AMPARs bound to glutamate and the prototypical negative allosteric modulator, GYKI-52466, with cryo-electron microscopy. Noncompetitive inhibition by GYKI-52466, which binds in the transmembrane collar region surrounding the ion channel, negatively modulates AMPARs by decoupling glutamate binding in the ligand binding domain from the ion channel. Furthermore, during allosteric competition between negative and positive modulators, negative allosteric modulation by GKYI-52466 outcompetes positive allosteric modulators to control AMPAR function. Our data provide a new framework for understanding allostery of AMPARs and foundations for rational design of therapeutics targeting AMPARs in neurological diseases.

Naturally mutagenic sequence diversity in a human type II topoisomerase.

Type II topoisomerases transiently cleave duplex DNA as part of a strand passage mechanism that helps control chromosomal organization and superstructure. Aberrant DNA cleavage can result in genomic instability, and how topoisomerase activity is controlled to prevent unwanted breaks is poorly understood. Using a genetic screen, we identified mutations in the beta isoform of human topoisomerase II (hTOP2ß) that render the enzyme hypersensitive to the chemotherapeutic agent etoposide. Several of these variants were unexpectedly found to display hypercleavage behavior in vitro and to be capable of inducing cell lethality in a DNA repair-deficient background; surprisingly, a subset of these mutations were also observed in TOP2B sequences from cancer genome databases. Using molecular dynamics simulations and computational network analyses, we found that many of the mutations obtained from the screen map to interfacial points between structurally coupled elements, and that dynamical modeling could be used to identify other damage-inducing TOP2B alleles present in cancer genome databases. This work establishes that there is an innate link between DNA cleavage predisposition and sensitivity to topoisomerase II poisons, and that certain sequence variants of human type II topoisomerases found in cancer cells can act as DNA-damaging agents. Our findings underscore the potential for hTOP2ß to function as a clastogen capable of generating DNA damage that may promote or support cellular transformation.

Conformational changes in the essential E. coli septal cell wall synthesis complex suggest an activation mechanism.

The bacterial divisome is a macromolecular machine composed of more than 30 proteins that controls cell wall constriction during division. Here, we present a model of the structure and dynamics of the core complex of the E. coli divisome, supported by a combination of structure prediction, molecular dynamics simulation, single-molecule imaging, and mutagenesis. We focus on the septal cell wall synthase complex formed by FtsW and FtsI, and its regulators FtsQ, FtsL, FtsB, and FtsN. The results indicate extensive interactions in four regions in the periplasmic domains of the complex. FtsQ, FtsL, and FtsB support FtsI in an extended conformation, with the FtsI transpeptidase domain lifted away from the membrane through interactions among the C-terminal domains. FtsN binds between FtsI and FtsL in a region rich in residues with superfission (activating) and dominant negative (inhibitory) mutations. Mutagenesis experiments and simulations suggest that the essential domain of FtsN links FtsI and FtsL together, potentially modulating interactions between the anchor-loop of FtsI and the putative catalytic cavity of FtsW, thus suggesting a mechanism of how FtsN activates the cell wall synthesis activities of FtsW and FtsI.

Phosphorylation at Ser65 modulates ubiquitin conformational dynamics.

Ubiquitin phosphorylation at Ser65 increases the population of a rare C-terminally retracted (CR) conformation. Transition between the Major and CR ubiquitin conformations is critical for promoting mitochondrial degradation. The mechanisms by which the Major and CR conformations of Ser65-phosphorylated (pSer65) ubiquitin interconvert, however, remain unresolved. Here, we perform all-atom molecular dynamics simulations using the string method with swarms of trajectories to calculate the lowest free-energy path between these two conformers. Our analysis reveals the existence of a Bent intermediate in which the C-terminal residues of the ß5 strand shift to resemble the CR conformation, while pSer65 retains contacts resembling the Major conformation. This stable intermediate was reproduced in well-tempered metadynamics calculations but was less stable for a Gln2Ala mutant that disrupts contacts with pSer65. Lastly, dynamical network modeling reveals that the transition from the Major to CR conformations involves a decoupling of residues near pSer65 from the adjacent ß1 strand.

Sodium/Iodide Symporter Metastable Intermediates Provide Insights into Conformational Transition between Principal Thermodynamic States.

The Sodium/Iodide Symporter (NIS), a 13-helix transmembrane protein found in the thyroid and other tissues, transports iodide, a required constituent of thyroid hormones T3 and T4. Despite extensive experimental information and clinical data, structural details of the intermediate microstates comprising the conformational transition of NIS between its inwardly and outwardly open states remain unresolved. We present data from a combination of enhanced sampling and transition path molecular dynamics (MD) simulations that elucidate the principal intermediate states comprising the inwardly to outwardly open transition of fully bound and apo NIS under an enforced ionic gradient. Our findings suggest that in both the absence and presence of bound physiological ions, NIS principally occupies a proximally inward to inwardly open state. When fully bound, NIS is also found to occupy a rare "inwardly occluded" state. The results of this work provide novel insight into the populations of NIS intermediates and the free energy landscape comprising the conformational transition, adding to a mechanistic understanding of NIS ion transport. Moreover, the knowledge gained from this approach can serve as a basis for studies of NIS mutants to target therapeutic interventions.

Excitatory and inhibitory D-serine binding to the NMDA receptor.

N-methyl-D-aspartate receptors (NMDARs) uniquely require binding of two different neurotransmitter agonists for synaptic transmission. D-serine and glycine bind to one subunit, GluN1, while glutamate binds to the other, GluN2. These agonists bind to the receptor's bi-lobed ligand-binding domains (LBDs), which close around the agonist during receptor activation. To better understand the unexplored mechanisms by which D-serine contributes to receptor activation, we performed multi-microsecond molecular dynamics simulations of the GluN1/GluN2A LBD dimer with free D-serine and glutamate agonists. Surprisingly, we observed D-serine binding to both GluN1 and GluN2A LBDs, suggesting that D-serine competes with glutamate for binding to GluN2A. This mechanism is confirmed by our electrophysiology experiments, which show that D-serine is indeed inhibitory at high concentrations. Although free energy calculations indicate that D-serine stabilizes the closed GluN2A LBD, its inhibitory behavior suggests that it either does not remain bound long enough or does not generate sufficient force for ion channel gating. We developed a workflow using pathway similarity analysis to identify groups of residues working together to promote binding. These conformation-dependent pathways were not significantly impacted by the presence of N-linked glycans, which act primarily by interacting with the LBD bottom lobe to stabilize the closed LBD.

Development and characterization of functional antibodies targeting NMDA receptors.

N-methyl-D-aspartate receptors (NMDARs) are critically involved in basic brain functions and neurodegeneration as well as tumor invasiveness. Targeting specific subtypes of NMDARs with distinct activities has been considered an effective therapeutic strategy for neurological disorders and diseases. However, complete elimination of off-target effects of small chemical compounds has been challenging and thus, there is a need to explore alternative strategies for targeting NMDAR subtypes. Here we report identification of a functional antibody that specifically targets the GluN1-GluN2B NMDAR subtype and allosterically down-regulates ion channel activity as assessed by electrophysiology. Through biochemical analysis, x-ray crystallography, single-particle electron cryomicroscopy, and molecular dynamics simulations, we show that this inhibitory antibody recognizes the amino terminal domain of the GluN2B subunit and increases the population of the non-active conformational state. The current study demonstrates that antibodies may serve as specific reagents to regulate NMDAR functions for basic research and therapeutic objectives.

Finding Druggable Sites in Proteins Using TACTICS.

Structure-based drug discovery efforts require knowledge of where drug-binding sites are located on target proteins. To address the challenge of finding druggable sites, we developed a machine-learning algorithm called TACTICS (trajectory-based analysis of conformations to identify cryptic sites), which uses an ensemble of molecular structures (such as molecular dynamics simulation data) as input. First, TACTICS uses k-means clustering to select a small number of conformations that represent the overall conformational heterogeneity of the data. Then, TACTICS uses a random forest model to identify potentially bindable residues in each selected conformation, based on protein motion and geometry. Lastly, residues in possible binding pockets are scored using fragment docking. As proof-of-principle, TACTICS was applied to the analysis of simulations of the SARS-CoV-2 main protease and methyltransferase and the Yersinia pestis aryl carrier protein. Our approach recapitulates known small-molecule binding sites and predicts the locations of sites not previously observed in experimentally determined structures. The TACTICS code is available at https://github.com/Albert-Lau-Lab/tactics_protein_analysis.

Distinct axial and lateral interactions within homologous filaments dictate the signaling specificity and order of the AIM2-ASC inflammasome.

Inflammasomes are filamentous signaling platforms integral to innate immunity. Currently, little is known about how these structurally similar filaments recognize and distinguish one another. A cryo-EM structure of the AIM2PYD filament reveals that the architecture of the upstream filament is essentially identical to that of the adaptor ASCPYD filament. In silico simulations using Rosetta and molecular dynamics followed by biochemical and cellular experiments consistently demonstrate that individual filaments assemble bidirectionally. By contrast, the recognition between AIM2 and ASC requires at least one to be oligomeric and occurs in a head-to-tail manner. Using in silico mutagenesis as a guide, we also identify specific axial and lateral interfaces that dictate the recognition and distinction between AIM2 and ASC filaments. Together, the results here provide a robust framework for delineating the signaling specificity and order of inflammasomes.

Structural biology and thermodynamics of GluD receptors.

Glutamate delta (GluD) receptors are a functionally enigmatic subfamily of ionotropic glutamate receptors. Despite sharing similar sequences and structures with AMPA, NMDA, and kainate receptors, GluD receptors do not bind glutamate nor function as ligand-gated ion channels. Binding d-serine and engaging in transsynaptic protein-protein interactions, GluD receptors are thought to undergo complex conformational rearrangements for non-ionotropic signaling that regulates synaptic plasticity. Recent structural, biochemical, and computational studies have elucidated multiple conformational and thermodynamic factors governing the unique properties of GluD receptors. Here, we review advances in biophysical insights into GluD receptors and discuss the structural thermodynamic relationships that underpin their neurobiological functions.

D-Serine Potently Drives Ligand-Binding Domain Closure in the Ionotropic Glutamate Receptor GluD2.

Despite their classification as ionotropic glutamate receptors, GluD receptors are not functional ligand-gated ion channels and do not bind glutamate. GluD2 receptors bind D-serine and coordinate transsynaptic complexes that regulate synaptic plasticity. Instead of opening the ion channel pore, mechanical tension produced from closure of GluD2 ligand-binding domains (LBDs) drives conformational rearrangements for non-ionotropic signaling. We report computed conformational free energy landscapes for the GluD2 LBD in apo and D-serine-bound states. Unexpectedly, the conformational free energy associated with GluD2 LBD closure upon D-serine binding is greater than that for AMPA, NMDA, and kainate receptor LBDs upon agonist binding. This excludes insufficient force generation as an explanation for lack of ion channel activity in GluD2 receptors and suggests that non-ionotropic conformational rearrangements do more work than pore opening. We also report free energy landscapes for GluD2 LBD harboring a neurodegenerative mutation and demonstrate selective stabilization of closed conformations in the apo state.

Enhanced sampling of glutamate receptor ligand-binding domains.

The majority of excitatory synaptic transmission in the central nervous system is mediated by ionotropic glutamate receptors (iGluRs). These membrane-bound protein assemblies consist of modular domains that can be genetically isolated and expressed, which has resulted in a plethora of crystal structures of individual domains in different conformations bound to different ligands. These structures have presented opportunities for molecular dynamics (MD) simulation studies. To examine the free energies that govern molecular behavior, simulation strategies and algorithms have been developed, collectively called enhanced sampling methods This review focuses on the use of enhanced sampling MD simulations of isolated iGluR ligand-binding domains to characterize thermodynamic properties important to receptor function.

High Conformational Variability in the GluK2 Kainate Receptor Ligand-Binding Domain.

The kainate family of ionotropic glutamate receptors (iGluRs) mediates pre- and postsynaptic neurotransmission. Previously computed conformational potentials of mean force (PMFs) for iGluR ligand-binding domains (LBDs) revealed subtype-dependent conformational differences between α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) and N-methyl-D-aspartic acid (NMDA) iGluR subfamilies. Here we report PMFs for the kainate receptor GluK2 in apo and glutamate-bound states. Apo and glutamate-bound GluK2 LBDs preferentially access closed-cleft conformations. Apo GluK2 exhibits a surprisingly high degree of conformational flexibility, accessing open and closed states. Comparing across iGluR subtypes, these results are similar to glycine-binding GluN1 and GluN3A NMDA subunits and differ from glutamate-binding GluA2 and GluN2A subunits. To test the contribution of cross-lobe interactions on closed-cleft LBD stability, we computed PMFs for two GluK2 mutants, D462A and D656S. D462A, but not D656S, weakens closed-cleft conformations of the glutamate-bound LBD. Theoretical Boltzmann-weighted small-angle X-ray scattering profiles improve agreement with experimental results compared with calculations from the LBD crystal structure alone.

Glutamate and Glycine Binding to the NMDA Receptor.

At central nervous system synapses, agonist binding to postsynaptic ionotropic glutamate receptors (iGluRs) results in signaling between neurons. N-Methyl-D-aspartic acid (NMDA) receptors are a unique family of iGluRs that activate in response to the concurrent binding of glutamate and glycine. Here, we investigate the process of agonist binding to the GluN2A (glutamate binding) and GluN1 (glycine binding) NMDA receptor subtypes using long-timescale unbiased molecular dynamics simulations. We find that positively charged residues on the surface of the GluN2A ligand-binding domain (LBD) assist glutamate binding via a "guided-diffusion" mechanism, similar in fashion to glutamate binding to the GluA2 LBD of AMPA receptors. Glutamate can also bind in an inverted orientation. Glycine, on the other hand, binds to the GluN1 LBD via an "unguided-diffusion" mechanism, whereby glycine finds its binding site primarily by random thermal fluctuations. Free energy calculations quantify the glutamate- and glycine-binding processes.

Neurotransmitter Funneling Optimizes Glutamate Receptor Kinetics.

Ionotropic glutamate receptors (iGluRs) mediate neurotransmission at the majority of excitatory synapses in the brain. Little is known, however, about how glutamate reaches the recessed binding pocket in iGluR ligand-binding domains (LBDs). Here we report the process of glutamate binding to a prototypical iGluR, GluA2, in atomistic detail using unbiased molecular simulations. Charged residues on the LBD surface form pathways that facilitate glutamate binding by effectively reducing a three-dimensional diffusion process to a spatially constrained, two-dimensional one. Free energy calculations identify residues that metastably bind glutamate and help guide it into the binding pocket. These simulations also reveal that glutamate can bind in an inverted conformation and also reorient while in its pocket. Electrophysiological recordings demonstrate that eliminating these transient binding sites slows activation and deactivation, consistent with slower glutamate binding and unbinding. These results suggest that binding pathways have evolved to optimize rapid responses of AMPA-type iGluRs at synapses.

Energetics of Glutamate Binding to an Ionotropic Glutamate Receptor.

Ionotropic glutamate receptors (iGluRs) are ligand-gated ion channels that are responsible for the majority of excitatory transmission at the synaptic cleft. Mechanically speaking, agonist binding to the ligand binding domain (LBD) activates the receptor by triggering a conformational change that is transmitted to the transmembrane region, opening the ion channel pore. We use fully atomistic molecular dynamics simulations to investigate the binding process in the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor, an iGluR subtype. The string method with swarms of trajectories was applied to calculate the possible pathways glutamate traverses during ligand binding. Residues peripheral to the binding cleft are found to metastably bind the ligand prior to ligand entry into the binding pocket. Umbrella sampling simulations were performed to compute the free energy barriers along the binding pathways. The calculated free energy profiles demonstrate that metastable interactions contribute substantially to the energetics of ligand binding and form local minima in the overall free energy landscape. Protein-ligand interactions at sites outside of the orthosteric agonist-binding site may serve to lower the transition barriers of the binding process.

Molecular lock regulates binding of glycine to a primitive NMDA receptor.

The earliest metazoan ancestors of humans include the ctenophore Mnemiopsis leidyi The genome of this comb jelly encodes homologs of vertebrate ionotropic glutamate receptors (iGluRs) that are distantly related to glycine-activated NMDA receptors and that bind glycine with unusually high affinity. Using ligand-binding domain (LBD) mutants for electrophysiological analysis, we demonstrate that perturbing a ctenophore-specific interdomain Arg-Glu salt bridge that is notably absent from vertebrate AMPA, kainate, and NMDA iGluRs greatly increases the rate of recovery from desensitization, while biochemical analysis reveals a large decrease in affinity for glycine. X-ray crystallographic analysis details rearrangements in the binding pocket stemming from the mutations, and molecular dynamics simulations suggest that the interdomain salt bridge acts as a steric barrier regulating ligand binding and that the free energy required to access open conformations in the glycine-bound LBD is largely responsible for differences in ligand affinity among the LBD variants.

Dynamics of the Ligand Binding Domain Layer during AMPA Receptor Activation.

Ionotropic glutamate receptors are postsynaptic tetrameric ligand-gated channels whose activity mediates fast excitatory transmission. Glutamate binding to clamshell-shaped ligand binding domains (LBDs) triggers opening of the integral ion channel, but how the four LBDs orchestrate receptor activation is unknown. Here, we present a high-resolution x-ray crystal structure displaying two tetrameric LBD arrangements fully bound to glutamate. Using a series of engineered metal ion trapping mutants, we showed that the more compact of the two assemblies corresponds to an arrangement populated during activation of full-length receptors. State-dependent cross-linking of the mutants identified zinc bridges between the canonical active LBD dimers that formed when the tetramer was either fully or partially bound by glutamate. These bridges also stabilized the resting state, consistent with the recently published full-length apo structure. Our results provide insight into the activation mechanism of glutamate receptors and the complex conformational space that the LBD layer can sample.

Conformational analysis of NMDA receptor GluN1, GluN2, and GluN3 ligand-binding domains reveals subtype-specific characteristics.

The NMDA receptor family of glutamate receptor ion channels is formed by obligate heteromeric assemblies of GluN1, GluN2, and GluN3 subunits. GluN1 and GluN3 bind glycine, whereas GluN2 binds glutamate. Crystal structures of the GluN1 and GluN3A ligand-binding domains (LBDs) in their apo states unexpectedly reveal open- and closed-cleft conformations, respectively, with water molecules filling the binding pockets. Computed conformational free energy landscapes for GluN1, GluN2A, and GluN3A LBDs reveal that the apo-state LBDs sample closed-cleft conformations, suggesting that their agonists bind via a conformational selection mechanism. By contrast, free energy landscapes for the AMPA receptor GluA2 LBD suggest binding of glutamate via an induced-fit mechanism. Principal component analysis reveals a rich spectrum of hinge bending, rocking, twisting, and sweeping motions that are different for the GluN1, GluN2A, GluN3A, and GluA2 LBDs. This variation highlights the structural complexity of signaling by glutamate receptor ion channels.

A conformational intermediate in glutamate receptor activation.

Ionotropic glutamate receptors (iGluRs) transduce the chemical signal of neurotransmitter release into membrane depolarization at excitatory synapses in the brain. The opening of the transmembrane ion channel of these ligand-gated receptors is driven by conformational transitions that are induced by the association of glutamate molecules to the ligand-binding domains (LBDs). Here, we describe the crystal structure of a GluA2 LBD tetramer in a configuration that involves an ∼30° rotation of the LBD dimers relative to the crystal structure of the full-length receptor. The configuration is stabilized by an engineered disulfide crosslink. Biochemical and electrophysiological studies on full-length receptors incorporating either this crosslink or an engineered metal bridge show that this LBD configuration corresponds to an intermediate state of receptor activation. GluA2 activation therefore involves a combination of both intra-LBD (cleft closure) and inter-LBD dimer conformational transitions. Overall, these results provide a comprehensive structural characterization of an iGluR intermediate state.