Steered Molecular Dynamics Simulations Predict Conformational Stability of Glutamate Receptors

Unlocking the Reverse Targeting Mechanisms of Cannabidiol: Unveiling New Therapeutic Avenues

Cannabidiol (CBD) and Δ9-tetrahydrocannabinol (THC), the main components of Cannabis sativa plants, have attracted a significant amount of attention due to their biological activities. This study identified GPR18 as the target of partial agonist CBD activating the p42/p44 MAPK pathway leading to migration of endometrial epithelial cells. Induced fit docking (IFD) showed that the affinity of THC for GPR18 is higher than that of CBD, and molecular dynamics (MD) simulations showed that CBD-GPR18 complexes at 130/200 ns might have stable conformations, potentially activating GPR18 by changing the distances of key residues in its active pocket. In contrast, THC maintains "metastable" conformations, generating a "shrinking space" leading to full agonism of THC by adding mechanical constraints in GPR18's active pocket. Steered molecular dynamics (SMD) revealed GPR18's active pocket was influenced more by CBD's partial agonism compared with THC. This combined IFD-MD-SMD method may be used to explain the mechanism of activation of partial or full agonists of GPR18.

Allosteric competition and inhibition in AMPA receptors

Excitatory neurotransmission is principally mediated by α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA)-subtype ionotropic glutamate receptors (AMPARs). Negative allosteric modulators are therapeutic candidates that inhibit AMPAR activation and can compete with positive modulators to control AMPAR function through unresolved mechanisms. Here we show that allosteric inhibition pushes AMPARs into a distinct state that prevents both activation and positive allosteric modulation. We used cryo-electron microscopy to capture AMPARs bound to glutamate, while a negative allosteric modulator, GYKI-52466, and positive allosteric modulator, cyclothiazide, compete for control of the AMPARs. GYKI-52466 binds in the ion channel collar and inhibits AMPARs by decoupling the ligand-binding domains from the ion channel. The rearrangement of the ligand-binding domains ruptures the cyclothiazide site, preventing positive modulation. Our data provide a framework for understanding allostery of AMPARs and for rational design of therapeutics targeting AMPARs in neurological diseases.

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.

The α-Synuclein Monomer May Have Different Misfolding Mechanisms in the Induction of α-Synuclein Fibrils with Different Polymorphs

The aggregation of alpha-synuclein (α-Syn) is closely related to the occurrence of some neurodegenerative diseases such as Parkinson's disease. The misfolding of α-Syn monomer plays a key role in the formation of aggregates and extension of fibril. However, the misfolding mechanism of α-Syn remains elusive. Here, three different α-Syn fibrils (isolated from a diseased human brain, generated by in vitro cofactor-tau induction, and obtained by in vitro cofactor-free induction) were selected for the study. The misfolding mechanisms of α-Syn were uncovered by studying the dissociation of the boundary chains based on the conventional molecular dynamics (MD) and Steered MD simulations. The results showed that the dissociation paths of the boundary chains in the three systems were different. According to the reverse process of dissociation, we concluded that in the human brain system, the binding of the monomer and template starts from the C-terminal and gradually misfolds toward the N-terminal. In the cofactor-tau system, the monomer binding starts from residues 58-66 (contain β3), followed by the C-terminal coil (residues 67-79). Then, the N-terminal coil (residues 36-41) and residues 50-57 (contain β2) bind to the template, followed by residues 42-49 (contain β1). In the cofactor-free system, two misfolding paths were found. One is that the monomer binds to the N/C-terminal (β1/β6) and then binds to the remaining residues. The other one is that the monomer binds sequentially from the C- to N-terminal, similar to the human brain system. Furthermore, in the human brain and cofactor-tau systems, electrostatic interactions (especially from residues 58-66) are the main driving force during the misfolding process, whereas in the cofactor-free system, the contributions of electrostatic and van der Waals interactions are comparable. These results may provide a deeper understanding for the misfolding and aggregation mechanism of α-Syn.

Molecular Dynamics Simulations of the Cardiac Ryanodine Receptor Type 2 (RyR2) Gating Mechanism

Mutations in the cardiac ryanodine receptor type 2 (RyR2) have been linked to fatal cardiac arrhythmias such as catecholaminergic polymorphic ventricular tachycardia (CPVT). While many CPVT mutations are associated with an increase in Ca²⁺ leak from the sarcoplasmic reticulum, the mechanistic details of RyR2 channel gating are not well understood, and this poses a barrier in the development of new pharmacological treatments. To address this, we explore the gating mechanism of the RyR2 using molecular dynamics (MD) simulations. We test the effect of changing the conformation of certain structural elements by constructing chimera RyR2 structures that are derived from the currently available closed and open cryo-electron microscopy (cryo-EM) structures, and we then use MD simulations to relax the system. Our key finding is that the position of the S4-S5 linker (S4S5L) on a single subunit can determine whether the channel as a whole is open or closed. Our analysis reveals that the position of the S4S5L is regulated by interactions with the U-motif on the same subunit and with the S6 helix on an adjacent subunit. We find that, in general, channel gating is crucially dependent on high percent occupancy interactions between adjacent subunits. We compare our interaction analysis to 49 CPVT1 mutations in the literature and find that 73% appear near a high percent occupancy interaction between adjacent subunits. This suggests that disruption of cooperative, high percent occupancy interactions between adjacent subunits is a primary cause of channel leak and CPVT in mutant RyR2 channels.

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.

Fluorine Substituted Proline Enhances the Tubulin Binding Potential of a Tetrapeptide at the GTP Binding Pocket Causing the Inhibition of Microtubule Motility and an Antimitotic Effect

The microtubule is regarded as the key target for designing anticancer and neurotherapeutic drugs due to its functional importance in eukaryotic cells including neurons. The microtubule is a dynamic hollow polymer tube consisting of α,β-tubulin heterodimer. Polymerization of α,β-tubulin heterodimer resulted in microtubule formation. GTP plays a crucial role in microtubule polymerization. It binds at the exchangeable binding site of the β-tubulin heterodimer, and it is one of the most crucial therapeutic hot spots for designing anticancer therapeutics. In this manuscript, we have shown using an in silico strategy and various in vitro and cellular experiments that the binding affinity to the tubulin and cancer therapeutic potential of an exchangeable GTP/GDP binding antimitotic tetrapeptide (SP: Ser-Leu-Arg-Pro) is increased through changing proline with the multifluorine substituted proline. This study showcases the importance of the proline amino acid and its pyrrolidine ring in the regulation of binding with tubulin at the GTP binding pocket.

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.

Thermodynamic Integration in 3n Dimensions without Biases or Alchemy for Protein Interactions

Thermodynamic integration (TI), a powerful formalism for computing Gibbs free energy, has been implemented for many biophysical processes with alchemical schemes that require delicate human efforts to choose/design biasing potentials for sampling the desired biophysical events and to remove their artifactitious consequences afterwards. Theoretically, an alchemical scheme is exact but practically, an unsophisticated implementation of this exact formula can cause error amplifications. Small relative errors in the input parameters can be amplified many times in their propagation into the computed free energy [due to subtraction of similar numbers such as (105 ± 5)‒(100 ± 5) = 5 ± 7]. In this paper, we present an unsophisticated implementation of TI in 3n dimensions (3nD) (n=1,2,3…) for the potential of mean force along a 3nD path connecting one state in the bound state ensemble to one state in the unbound state ensemble. Fluctuations in these 3nD are integrated in the bound and unbound state ensembles but not along the 3nD path. Using TI3nD, we computed the standard binding free energies of three protein complexes: trometamol in Salmonella effector SpvD (n=1), biotin-avidin (n=2), and Colicin E9 endonuclease with cognate immunity protein Im9 (n=3). We employed three different protocols in three independent computations of E9-Im9 to show TI3nD's robustness. We also computed the hydration energies of ten biologically relevant compounds (n=1 for water, acetamide, urea, glycerol, trometamol, ammonium and n=2 for erythritol, 1,3-propanediol, xylitol, biotin). Each of the 15 computations is accomplishable within one (for hydration) to ten (for E9-Im9) days on an inexpensive GPU workstation. The computed results all agree with the available experimental data.

The structural arrangement at intersubunit interfaces in homomeric kainate receptors

Kainate receptors are glutamate-gated cation-selective channels involved in excitatory synaptic signaling and are known to be modulated by ions. Prior functional and structural studies suggest that the dimer interface at the agonist-binding domain plays a key role in activation, desensitization, and ion modulation in kainate receptors. Here we have used fluorescence-based methods to investigate the changes and conformational heterogeneity at these interfaces associated with the resting, antagonist-bound, active, desensitized, and ion-modulated states of the receptor. These studies show that in the presence of Na+ ions the interfaces exist primarily in the coupled state in the apo, antagonist-bound and activated (open channel) states. Under desensitizing conditions, the largely decoupled dimer interface at the agonist-binding domain as seen in the cryo-EM structure is one of the states observed. However, in addition to this state there are several additional states with lower levels of decoupling. Replacing Na+ with Cs+ does not alter the FRET efficiencies of the states significantly, but shifts the population to the more decoupled states in both resting and desensitized states, which can be correlated with the lower activation seen in the presence of Cs+.

Disclosing the Template-Induced Misfolding Mechanism of Tau Protein by Studying the Dissociation of the Boundary Chain from the Formed Tau Fibril Based on a Steered Molecular Dynamics Simulation

The level of tau aggregation into neurofibrillary tangles, including paired helical filament (PHF) and straight filament (SF), is closely associated with Alzheimer's disease. Despite the pathological importance of misfolding and aggregation of tau, the corresponding mechanism remains unclear. Therefore, to uncover the misfolding mechanism of the tau monomer upon induction of formed PHF and SF, in this study, a conventional molecular dynamics simulation combined with a steered molecular dynamics simulation was performed to study the dissociation of the boundary chain. Interestingly, our results show that the dissociation mechanisms of the boundary chain in PHF and SF are different. In PHF, the boundary chain begins to dissociate from regions β2 and β3 and ends at β8. However, in SF, it is simultaneously dissociated from β1 and β8 and ends at β5. The dissociation of the boundary chain is the reverse of template-induced misfolding of the monomer. Therefore, we can deduce the misfolding mechanism of the monomer upon induction of the template. For PHF, β8 first interacts with the template by hydrophobic interaction. Then β7, β6, β5, β4, and β1 sequentially bind to the template by electrostatic and hydrophobic interactions. After β1 binds to the template, β2 and β3 very quickly bind to the template through hydrophobic interaction. For SF, β5 of the monomer first interacts with the template by electrostatic attraction. Then β4 and β6, β3 and β7, and β2 and β8 bind to the template in turn. Finally, β1 and β8 are fully bound to the template by hydrophobic interaction. The obtained results will be vital for understanding the earlier events during misfolding and aggregation of tau.

Polyamine-mediated channel block of ionotropic glutamate receptors and its regulation by auxiliary proteins

Most excitatory neurotransmission in the mammalian brain is mediated by a family of plasma membrane-bound signaling proteins called ionotropic glutamate receptors (iGluRs). iGluRs assemble at central synapses as tetramers, forming a central ion-channel pore whose primary function is to rapidly transport Na⁺ and Ca²⁺ in response to binding the neurotransmitter l-glutamic acid. The pore of iGluRs is also accessible to bulkier cytoplasmic cations, such as the polyamines spermine, spermidine, and putrescine, which are drawn into the permeation pathway, but get stuck and block the movement of other ions. The degree of this polyamine-mediated channel block is highly regulated by processes that control the free cytoplasmic polyamine concentration, the membrane potential, or the iGluR subunit composition. Recently, an additional regulation by auxiliary proteins, most notably transmembrane AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid) receptor regulatory proteins (TARPs), cornichons, and neuropilin and tolloid-like proteins (NETOs), has been identified. Here, I review what we have learned of polyamine block of iGluRs and its regulation by auxiliary subunits. TARPs, cornichons, and NETOs attenuate the channel block by enabling polyamines to exit the pore. As a result, polyamine permeation occurs at more negative and physiologically relevant membrane potentials. The structural basis for enhanced polyamine transport remains unresolved, although alterations in both channel architecture and charge-screening mechanisms have been proposed. That auxiliary subunits can attenuate the polyamine block reveals an unappreciated impact of polyamine permeation in shaping the signaling properties of neuronal AMPA- and kainate-type iGluRs. Moreover, enhanced polyamine transport through iGluRs may have a role in regulating cellular polyamine levels.

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.

Functional Validation of Heteromeric Kainate Receptor Models

Kainate receptors require the presence of external ions for gating. Most work thus far has been performed on homomeric GluK2 but, in vivo, kainate receptors are likely heterotetramers. Agonists bind to the ligand-binding domain (LBD) which is arranged as a dimer of dimers as exemplified in homomeric structures, but no high-resolution structure currently exists of heteromeric kainate receptors. In a full-length heterotetramer, the LBDs could potentially be arranged either as a GluK2 homomer alongside a GluK5 homomer or as two GluK2/K5 heterodimers. We have constructed models of the LBD dimers based on the GluK2 LBD crystal structures and investigated their stability with molecular dynamics simulations. We have then used the models to make predictions about the functional behavior of the full-length GluK2/K5 receptor, which we confirmed via electrophysiological recordings. A key prediction and observation is that lithium ions bind to the dimer interface of GluK2/K5 heteromers and slow their desensitization.

Unitary Properties of AMPA Receptors with Reduced Desensitization

Wild-type AMPA receptors display a characteristic rapidly desensitizing phenotype. Many studies point to the dimer interface between pairs of extracellular ligand binding domains as the key region controlling the rate at which the receptors desensitize. However, mutations at the extracellular end of the pore-forming regions (near the putative ion channel gate) have also been shown to alter desensitization. Here we report the behavior of single GluA4 receptors carrying one of two mutations that greatly reduce desensitization at the level of ensemble currents: the dimer interface mutation L484Y and the Lurcher mutation (A623T, GluA4-Lc) in the extracellular end of M3 (the second true transmembrane helix). Analysis of unitary currents in patches with just one active receptor showed that each mutation greatly prolongs bursts of openings without prolonging the apparent duration of individual openings. Each mutation decreases the frequency with which individual receptors visit desensitized states, but both mutant receptors still desensitize multiple times per second. Cyclothiazide (CTZ) reduced desensitization of wild-type receptors and both types of mutant receptor. Analysis of shut-time distributions revealed a form of short-lived desensitization that was resistant to CTZ and was especially prominent for GluA4-Lc receptors. Despite reducing desensitization of GluA4 L484Y receptors, CTZ decreased the amplitude of ensemble currents through GluA2 and GluA4 LY receptor mutants. Single-channel analysis and comparison of the GluA2 L483Y ligand binding domain dimer in complex with glutamate with and without CTZ is consistent with the conclusion that CTZ binding to the dimer interface prevents effects of the LY mutation to modulate receptor activation, resulting in a reduction in the prevalence of large-conductance substates that accounts for the decrease in ensemble current amplitudes. Together, the results show that similar nondesensitizing AMPA-receptor phenotypes of population currents can arise from distinct underlying molecular mechanisms that produce different types of unitary activity.

Configurational Preference of the Glutamate Receptor Ligand Binding Domain Dimers

Ionotropic glutamate receptors are a family of tetrameric ion channels with functional states consisting of nonconducting, conducting, and desensitized states that are starting to become well characterized by electrophysiological and biophysical studies. However, the structure and relative energetics of these states beyond the general structure of the receptor are still not well understood. It is known that the interface between monomeric subunits of the tetramer plays a major role in distinguishing these functional states. We have used umbrella sampling and multimicrosecond molecular dynamics simulations of the GluA2 AMPA subtype glutamate receptor ligand-binding domain (LBD) dimers to characterize a natural propensity of the LBD dimers for various configurational states. Our results show a proposed desensitized conformation of the LBD dimer is a highly preferable conformation of the LBD dimer without the influence of other receptor domains or crystallographic conditions. This has been demonstrated by both free protein simulations of 5 μs duration, as well as by computed free energy difference between the active and desensitized states. At the same time, the simulations performed using the same protocols revealed that for the LBD mutant L483Y, known to lack desensitization, the postulated active state of the LBD dimer is indeed the preferred configurational state, which remained stable in the simulations. Our findings pave the path for developing more detailed hypotheses of the full receptor activation mechanism. Combined with the energetics of glutamate binding to the LBD and the energy required to open the transmembrane pore helices, our results strongly support a hypothesis that the low absolute free-energy state is the desensitized state of the intact AMPA receptor.

Role of the Ion Channel Extracellular Collar in AMPA Receptor Gating

AMPA subtype ionotropic glutamate receptors mediate fast excitatory neurotransmission and are implicated in numerous neurological diseases. Ionic currents through AMPA receptor channels can be allosterically regulated via different sites on the receptor protein. We used site-directed mutagenesis and patch-clamp recordings to probe the ion channel extracellular collar, the binding region for noncompetitive allosteric inhibitors. We found position and substitution-dependent effects for introduced mutations at this region on AMPA receptor gating. The results of mutagenesis suggested that the transmembrane domains M1, M3 and M4, which contribute to the ion channel extracellular collar, undergo significant relative displacement during gating. We used molecular dynamics simulations to predict an AMPA receptor open state structure and rationalize the results of mutagenesis. We conclude that the ion channel extracellular collar plays a distinct role in gating and represents a hub for powerful allosteric modulation of AMPA receptor function that can be used for developing novel therapeutics.

Computing the binding affinity of a ligand buried deep inside a protein with the hybrid steered molecular dynamics

Computing the ligand-protein binding affinity (or the Gibbs free energy) with chemical accuracy has long been a challenge for which many methods/approaches have been developed and refined with various successful applications. False positives and, even more harmful, false negatives have been and still are a common occurrence in practical applications. Inevitable in all approaches are the errors in the force field parameters we obtain from quantum mechanical computation and/or empirical fittings for the intra- and inter-molecular interactions. These errors propagate to the final results of the computed binding affinities even if we were able to perfectly implement the statistical mechanics of all the processes relevant to a given problem. And they are actually amplified to various degrees even in the mature, sophisticated computational approaches. In particular, the free energy perturbation (alchemical) approaches amplify the errors in the force field parameters because they rely on extracting the small differences between similarly large numbers. In this paper, we develop a hybrid steered molecular dynamics (hSMD) approach to the difficult binding problems of a ligand buried deep inside a protein. Sampling the transition along a physical (not alchemical) dissociation path of opening up the binding cavity---pulling out the ligand---closing back the cavity, we can avoid the problem of error amplifications by not relying on small differences between similar numbers. We tested this new form of hSMD on retinol inside cellular retinol-binding protein 1 and three cases of a ligand (a benzylacetate, a 2-nitrothiophene, and a benzene) inside a T4 lysozyme L99A/M102Q(H) double mutant. In all cases, we obtained binding free energies in close agreement with the experimentally measured values. This indicates that the force field parameters we employed are accurate and that hSMD (a brute force, unsophisticated approach) is free from the problem of error amplification suffered by many sophisticated approaches in the literature.