Energetics of Glutamate Binding to an Ionotropic Glutamate Receptor

AMPA receptor neurotransmission and therapeutic applications: A comprehensive review of their multifaceted modulation

The neuropharmacological community has shown a strong interest in AMPA receptors as critical components of excitatory synaptic transmission during the last fifteen years. AMPA receptors, members of the ionotropic glutamate receptor family, allow rapid excitatory neurotransmission in the brain. AMPA receptors, which are permeable to sodium and potassium ions, manage the bulk of the brain's rapid synaptic communications. This study thoroughly examines the recent developments in AMPA receptor regulation, focusing on a shift from single chemical illustrations to a more extensive investigation of underlying processes. The complex interplay of these modulators in modifying the function and structure of AMPA receptors is the main focus, providing insight into their influence on the speed of excitatory neurotransmission. This research emphasizes the potential of AMPA receptor modulation as a therapy for various neurological disorders such as epilepsy and Alzheimer's disease. Analyzing these regulators' sophisticated molecular details enhances our comprehension of neuropharmacology, representing a significant advancement in using AMPA receptors for treating intricate neurological conditions.

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.

Sustained-release behavior and the antitumor effect of charge-convertible poly(amino acid)s drug-loaded nanoparticles

Enhancing tissue permeability and achieving drug aggregation is the key to targeted tumor therapy. A series triblock copolymers of poly(ethylene glycol)-poly(L-lysine)-poly(L-glutamine) were synthesized by ring-opening polymerization, and charge-convertible nano-delivery system was constructed by loading doxorubicin (DOX) with 2-(hexaethylimide) ethanol on side chain. In normal environment (pH = 7.4), the zeta potential of the drug-loaded nanoparticle solution is negative, which is conducive to avoiding the identification and clearance of nanoparticles by the reticulo-endothelial system, while potential-reversal can be achieved in the tumor microenvironment, which effectively promotes cellular uptake. Nanoparticles could effectively reduce the distribution of DOX in normal tissues and achieve targeted aggregation at tumor sites, which can effectively improve the antitumor effect, while would not causing toxicity and damage to normal body.

Determination of Multidirectional Pathways for Ligand Release from the Receptor: A New Approach Based on Differential Evolution

Steered molecular dynamics (SMD) simulation is a powerful method in computer-aided drug design as it can be used to access the relative binding affinity with high precision but with low computational cost. The success of SMD depends on the choice of the direction along which the ligand is pulled from the receptor-binding site. In most simulations, the unidirectional pathway was used, but in some cases, this choice resulted in the ligand colliding with the complex surface of the exit tunnel. To overcome this difficulty, several variants of SMD with multidirectional pulling have been proposed, but they are not completely devoid of disadvantages. Here, we have proposed to determine the direction of pulling with a simple scoring function that minimizes the receptor–ligand interaction, and an optimization algorithm called differential evolution is used for energy minimization. The effectiveness of our protocol was demonstrated by finding expulsion pathways of Huperzine A and camphor from the binding site of Torpedo California acetylcholinesterase and P450cam proteins, respectively, and comparing them with the previous results obtained using memetic sampling and random acceleration molecular dynamics. In addition, by applying this protocol to a set of ligands bound with LSD1 (lysine specific demethylase 1), we obtained a much higher correlation between the work of pulling force and experimental data on the inhibition constant IC50 compared to that obtained using the unidirectional approach based on minimal steric hindrance.

Structure, Function, and Pharmacology of Glutamate Receptor Ion Channels

Many physiologic effects of l-glutamate, the major excitatory neurotransmitter in the mammalian central nervous system, are mediated via signaling by ionotropic glutamate receptors (iGluRs). These ligand-gated ion channels are critical to brain function and are centrally implicated in numerous psychiatric and neurologic disorders. There are different classes of iGluRs with a variety of receptor subtypes in each class that play distinct roles in neuronal functions. The diversity in iGluR subtypes, with their unique functional properties and physiologic roles, has motivated a large number of studies. Our understanding of receptor subtypes has advanced considerably since the first iGluR subunit gene was cloned in 1989, and the research focus has expanded to encompass facets of biology that have been recently discovered and to exploit experimental paradigms made possible by technological advances. Here, we review insights from more than 3 decades of iGluR studies with an emphasis on the progress that has occurred in the past decade. We cover structure, function, pharmacology, roles in neurophysiology, and therapeutic implications for all classes of receptors assembled from the subunits encoded by the 18 ionotropic glutamate receptor genes. SIGNIFICANCE STATEMENT: Glutamate receptors play important roles in virtually all aspects of brain function and are either involved in mediating some clinical features of neurological disease or represent a therapeutic target for treatment. Therefore, understanding the structure, function, and pharmacology of this class of receptors will advance our understanding of many aspects of brain function at molecular, cellular, and system levels and provide new opportunities to treat patients.

Neural Network Sampling of the Free Energy Landscape for Nitrogen Dissociation on Ruthenium

In heterogeneous catalysis, free energy profiles of reactions govern the mechanisms, rates, and equilibria. Energetics are conventionally computed using the harmonic approximation (HA), which requires determination of critical states a priori. Here, we use neural networks to efficiently sample and directly calculate the free energy surface (FES) of a prototypical heterogeneous catalysis reaction-the dissociation of molecular nitrogen on ruthenium-at density-functional-theory-level accuracy. We find that the vibrational entropy of surface atoms, often neglected in HA for transition metal catalysts, contributes significantly to the reaction barrier. The minimum free energy path for dissociation reveals an "on-top" adsorbed molecular state prior to the transition state. While a previously reported flat-lying molecular metastable state can be identified in the potential energy surface, it is absent in the FES at relevant reaction temperatures. These findings demonstrate the importance of identifying critical points self-consistently on the FES for reactions that involve considerable entropic effects.

A multiscale coarse-grained model of the SARS-CoV-2 virion

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and ongoing development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data become publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion.

A Multiscale Coarse-grained Model of the SARS-CoV-2 Virion

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations, however, are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and on-going development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Structural data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data becomes publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion.

Significance Statement

This study reports the construction of a molecular model for the SARS-CoV-2 virion and details our multiscale approach towards model refinement. The resulting model and methods can be applied to and enable the simulation of SARS-CoV-2 virions.

Atomic-scale characterization of mature HIV-1 capsid stabilization by inositol hexakisphosphate (IP6)

Inositol hexakisphosphates (IP6) are cellular cofactors that promote the assembly of mature capsids of HIV. These negatively charged molecules coordinate an electropositive ring of arginines at the center of pores distributed throughout the capsid surface. Kinetic studies indicate that the binding of IP6 increases the stable lifetimes of the capsid by several orders of magnitude from minutes to hours. Using all-atom molecular dynamics simulations, we uncover the mechanisms that underlie the unusually high stability of mature capsids in complex with IP6 We find that capsid hexamers and pentamers have differential binding modes for IP6 Ligand density calculations show three sites of interaction with IP6 including at a known capsid inhibitor binding pocket. Free energy calculations demonstrate that IP6 preferentially stabilizes pentamers over hexamers to enhance fullerene modes of assembly. These results elucidate the molecular role of IP6 in stabilizing and assembling the retroviral capsid.

Druggability Simulations and X-Ray Crystallography Reveal a Ligand-Binding Site in the GluA3 AMPA Receptor N-Terminal Domain

Ionotropic glutamate receptors (iGluRs) mediate the majority of excitatory neurotransmission in the brain. Their dysfunction is implicated in many neurological disorders, rendering iGluRs potential drug targets. Here, we performed a systematic analysis of the druggability of two major iGluR subfamilies, using molecular dynamics simulations in the presence of drug-like molecules. We demonstrate the applicability of druggability simulations by faithfully identifying known agonist and modulator sites on AMPA receptors (AMPARs) and NMDA receptors. Simulations produced the expected allosteric changes of the AMPAR ligand-binding domain in response to agonist. We also identified a novel ligand-binding site specific to the GluA3 AMPAR N-terminal domain (NTD), resulting from its unique conformational flexibility that we explored further with crystal structures trapped in vastly different states. In addition to providing an in-depth analysis into iGluR NTD dynamics, our approach identifies druggable sites and permits the determination of pharmacophoric features toward novel iGluR modulators.

Computational conformational analysis of α-thrombin inhibitors possessing distinct scaffolds in aqueous solution and on Ala-sheet

Recent computational simulations on protein-ligand binding/unbinding have precisely been uncovering the ligand-binding process at the atomic level. In the process, the non-specific binding of ligands to the target site is suggested to occur before binding to the target. We in this study analyzed the conformations of ligands under the non-specific binding on a protein surface to figure out the differences in the conformational characteristics in aqueous solution using the 55-ns molecular dynamic simulation. As for the protein surface, we constructed an artificial β-sheet, composed of poly-alanine residues (Ala-sheet). For the ligands, the four α-thrombin inhibitors possessing two scaffolds with distinct hydrophobicity profiles were used. During the simulation, all the inhibitors kept interaction with Ala-sheet and had the limited conformational fluctuations compared with in aqueous solution. The representative conformations obtained from the cluster analysis showed that two of hydrophobic inhibitors adopted the extended conformations in aqueous solution and also on Ala-sheet. For the other two hydrophilic inhibitors, the conformations in aqueous solution adopted the bent conformation with two terminal hydrophobic rings closely packed. On Ala-sheet, contrarily, the two hydrophobic rings were open and took the extended conformations, which were placed on the sheet as a foothold. The charged moieties in the hydrophilic inhibitors were protruded into aqueous environment with the extended conformation. The conformational characteristics of the inhibitors in aqueous solution and Ala-sheet varied likely by chemical features or structures of the inhibitors, but each was considered to be physicochemically reasonable.

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.

Activation and desensitization of ionotropic glutamate receptors by selectively triggering pre-existing motions

Ionotropic glutamate receptors (iGluRs) are ligand-gated ion channels that are key players in synaptic transmission and plasticity. They are composed of four subunits, each containing four functional domains, the quaternary packing and collective structural dynamics of which are important determinants of their molecular mechanism of function. With the explosion of structural studies on different members of the family, including the structures of activated open channels, the mechanisms of action of these central signaling machines are now being elucidated. We review the current state of computational studies on two major members of the family, AMPA and NMDA receptors, with focus on molecular simulations and elastic network model analyses that have provided insights into the coupled movements of extracellular and transmembrane domains. We describe the newly emerging mechanisms of activation, allosteric signaling and desensitization, as mainly a selective triggering of pre-existing soft motions, as deduced from computational models and analyses that leverage structural data on intact AMPA and NMDA receptors in different states.