Molecular dynamics simulations of dihydro-β-erythroidine bound to the human α4β2 nicotinic acetylcholine receptor

Mitophagy in perioperative neurocognitive disorder: mechanisms and therapeutic strategies

Perioperative neurocognitive disorder (PND) is a common neurological complication after surgery/anesthesia in elderly patients that affect postoperative outcome and long-term quality of life, which increases the cost of family and social resources. The pathological mechanism of PND is complex and not fully understood, and the methods of prevention and treatment of PND are very limited, so it is particularly important to analyze the mechanism of PND. Research indicates that mitochondrial dysfunction is pivotal in the initiation and progression of PND, although the precise mechanisms remain elusive and could involve disrupted mitophagy. We reviewed recent studies on the link between mitophagy and PND, highlighting the role of key proteins in abnormal mitophagy and discussing therapeutic strategies aimed at mitophagy regulation. This provides insights into the mechanisms underlying PND and potential therapeutic targets.

Dissociation of Nicotine from Acetylcholine-Binding Protein under Terahertz Waves Radiation

The binding of nicotine (NCT) to acetylcholine-binding protein (AChBP) plays an important role in synaptic transmission and neurotransmitter regulation. However, effectively regulating their binding or dissociation processes remains a challenging problem. In this study, we employed all-atom molecular dynamics (MD) simulations to systematically investigate the impact of external terahertz (THz) waves on the binding kinetics between AChBP and NCT. We first identified the key residues (i.e., W143) and the key interactions (i.e., hydrogen bonding and cation-π interaction) in AChBP-NCT binding without THz waves. We then investigated the binding and dissociation of charged NCT with AChBP at three different frequencies (i.e., 13.02, 21.44, 42.55 THz). Importantly, the predominant vibrational modes at 13.02 THz can drive the rotation of the pentagonal ring on NCT. This leads to the disruption of hydrogen bonds between NCT and W143 and a reduced likelihood of forming cation-π interactions, resulting in the dissociation of NCT from AChBP. Additionally, we further investigated the influence of electric field intensities on the dissociation kinetics and found that when the electric field intensity exceeds a critical value (∼0.60 V/nm), the probability of ligand dissociation gradually rises as the intensity increases. In general, this study contributes to a better understanding of the effects of THz waves on protein-ligand interactions, which might also shed some light on potential applications in nicotine addiction treatment and therapeutic strategies for neurodegenerative diseases.

Molecular Dynamics Refinement of Open State Serotonin 5-HT3A Receptor Structures

Pentameric ligand-gated ion channels play an important role in mediating fast neurotransmissions. As a member of this receptor family, cation-selective 5-HT₃ receptors are a clinical target for treating nausea and vomiting associated with chemotherapy and radiation therapy (Thompson and Lummis, 2006). Multiple cryo-electron microscopy (cryo-EM) structures of 5-HT₃ receptors have been determined in distinct functional states (e.g., open, closed, etc.) (Basak et al., 2018; Basak et al., 2018; Polovinkin et al., 2018; Zhang et al., 2015). However, recent work has shown that the transmembrane pores of the open 5-HT₃ receptor structures rapidly collapse and become artificially asymmetric in molecular dynamics (MD) simulations. To avoid this hydrophobic collapse, Dämgen and Biggin developed an equilibration protocol that led to a stable open state structure of the glycine receptor in MD simulations (Dämgen and Biggin, 2020). However, the protocol failed to yield open-like structures of the 5-HT₃ receptor in our simulations. Here, we present a refined equilibration protocol that involves the rearrangement of the transmembrane helices to achieve stable open state structures of the 5-HT₃ receptor that allow both water and ion permeation through the channel. Notably, channel gating is mediated through collective movement of the transmembrane helices, involving not only pore lining M2 helices but also their cross-talk with the adjacent M1 and M3 helices. Thus, the successful application of our refined equilibration protocol underscores the importance of the conformational coupling between the transmembrane helices in stabilizing open-like structures of the 5-HT₃ receptor.

Design, synthesis, and mechanism of action of novel μ-conotoxin KIIIA analogues for inhibition of the voltage-gated sodium channel Nav1.7

μ-Conotoxin KIIIA, a selective blocker of sodium channels, has strong inhibitory activity against several Na_v isoforms, including Na_v1.7, and has potent analgesic effects, but it contains three pairs of disulfide bonds, making structural modification difficult and synthesis complex. To circumvent these difficulties, we designed and synthesized three KIIIA analogues with one disulfide bond deleted. The most active analogue, KIIIA-1, was further analyzed, and its binding pattern to hNa_v1.7 was determined by molecular dynamics simulations. Guided by the molecular dynamics computational model, we designed and tested 32 second-generation and 6 third-generation analogues of KIIIA-1 on hNa_v1.7 expressed in HEK293 cells. Several analogues showed significantly improved inhibitory activity on hNa_v1.7, and the most potent peptide, 37, was approximately 4-fold more potent than the KIIIA Isomer I and 8-fold more potent than the wildtype (WT) KIIIA in inhibiting hNa_v1.7 current. Intraperitoneally injected 37 exhibited potent in vivo analgesic activity in a formalin-induced inflammatory pain model, with activity reaching ∼350-fold of the positive control drug morphine. Overall, peptide 37 has a simplified disulfide-bond framework and exhibits potent in vivo analgesic effects and has promising potential for development as a pain therapy in the future.

Structural Insights into the Role of β3 nAChR Subunit in the Activation of Nicotinic Receptors

The β3 subunit of nicotinic acetylcholine receptors (nAChRs) participates in heteropentameric assemblies with some α and other β neuronal subunits forming a plethora of various subtypes, differing in their electrophysiological and pharmacological properties. While β3 has for several years been considered an accessory subunit without direct participation in the formation of functional binding sites, recent electrophysiology data have disputed this notion and indicated the presence of a functional (+) side on the extracellular domain (ECD) of β3. In this study, we present the 2.4 Å resolution crystal structure of the monomeric β3 ECD, which revealed rather distinctive loop C features as compared to those of α nAChR subunits, leading to intramolecular stereochemical hindrance of the binding site cavity. Vigorous molecular dynamics simulations in the context of full length pentameric β3-containing nAChRs, while not excluding the possibility of a β3 (+) binding site, demonstrate that this site cannot efficiently accommodate the agonist nicotine. From the structural perspective, our results endorse the accessory rather than functional role of the β3 nAChR subunit, in accordance with earlier functional studies on β3-containing nAChRs.

Mechanism of interactions between α-conotoxin RegIIA and carbohydrates at the human α3β4 nicotinic acetylcholine receptor

Conotoxins are marine peptide toxins from marine cone snails. The α-conotoxin RegIIA can selectively act on human (h) α3β4 nicotinic acetylcholine receptor (nAChR), and is an important lead for drug development. The high-resolution cryo-electron microscopy structure of the α3β4 nAChR demonstrates several carbohydrates are located near the orthosteric binding sites, which may affect α-conotoxin binding. Oligosaccharide chains can modify the physical and chemical properties of proteins by changing the conformation, hydrophobicity, quality and size of the protein. The purpose of this study is to explore the effect of oligosaccharide chains on the binding modes and activities of RegIIA and its derivatives at hα3β4 nAChRs. Through computational simulations, we designed and synthesized RegIIA mutants at position 14 to explore the importance of residue H14 to the activity of the peptide. Molecular dynamics simulations suggest that the oligosaccharide chains affect the binding of RegIIA at the hα3β4 nAChR through direct interactions with H14 and by affecting the C-loop conformation of the binding sites. Electrophysiology studies on H14 analogues suggest that in addition to forming direct interactions with the carbohydrates, the residue might play an important role in maintaining the conformation of the peptide. Overall, this study further clarifies the structure-activity relationship of α-conotoxin RegIIA at the hα3β4 nAChR and, also provides important experimental and theoretical basis for the development of new peptide drugs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-021-00108-9.

Antagonistic Mechanism of α-Conotoxin BuIA toward the Human α3β2 Nicotinic Acetylcholine Receptor

Nicotinic acetylcholine receptors (nAChRs) are pentameric ligand-gated ion channels that are abundantly expressed in the central and peripheral nervous systems, playing an important role in mediating neurotransmitter release and inter-synaptic signaling. Dysfunctional nAChRs are associated with neurological disorders, and studying the structure and function of nAChRs is essential for development of drugs or strategies for treatment of related diseases. α-Conotoxins are selective antagonists of the nAChR and are an important class of drug leads. So far, the antagonistic mechanism of α-conotoxins toward the nAChRs is still unclear. In this study, we built an α3β2 nAChR homology model and investigated its conformational transition mechanism upon binding with a highly potent inhibitor, α-conotoxin BuIA, through μs molecular dynamic simulations and site-directed mutagenesis studies. The results suggested that the α3β2 nAChR underwent global conformational transitions and was stabilized into a closed state with three hydrophobic gates present in the transmembrane domain by BuIA. Finally, the probable antagonistic mechanism of BuIA was proposed. Overall, the closed-state model of the α3β2 nAChR bound with BuIA is not only essential for understanding the antagonistic mechanism of α-conotoxins but also particularly valuable for development of therapeutic inhibitors in future.

Simulations of Promising Indolizidine-α6-β2 Nicotinic Acetylcholine Receptor Complexes

Smoking-cessation drugs bind many off-target nicotinic acetylcholine receptors (nAChRs) and cause severe side effects if they are based on nicotine. New drugs that bind only those receptors, such as α6β2* nAChR, implicated in nicotine addiction would avoid the off-target binding. Indolizidine (-)-237D (IND (-)-237D), a bicyclic alkaloid, has been shown to block α6β2* containing nAChRs and functionally inhibit the nicotine-evoked dopamine release. To improve the affinity of indolizidine (-)-237D for α6β2*, we built a library of 2226 analogs. We screened virtually the library against a homology model of α6β2 nAChR that we derived from the recent crystal structure of α4β2 nAChR. We also screened the crystal structure of α4β2 nAChR as a control on specificity. We ranked the compounds based on their predicted free energy of binding. We selected the top eight compounds bound in their best pose and subjected the complexes to 100 ns molecular dynamics simulations to assess the stability of the complexes. All eight analogs formed stable complexes for the duration of the simulations. The results from this work highlight nine distinct analogs of IND (-)-237D with high affinity towards α6β2* nAChR. These leads can be synthesized and tested in in vitro and in vivo studies as lead candidates for drugs to treat nicotine addiction.

Erythrina Alkaloid Analogues as nAChR Antagonists-A Flexible Platform for Leads in Drug Discovery

Erythrina alkaloids and their central nervous system effects have been studied for over a century, mainly due to their potent antagonistic actions at β2-containing nicotinic acetylcholine receptors (nAChRs). In the present work, we report a synthetic approach giving access to a diverse set of Erythrina natural product analogues and present the enantioselective total synthesis of (+)-Cocculine and (+)-Cocculidine, both found to be potent antagonists of the β2-containing nAChRs.

Structure and gating mechanism of the α7 nicotinic acetylcholine receptor

The α7 nicotinic acetylcholine receptor plays critical roles in the central nervous system and in the cholinergic inflammatory pathway. This ligand-gated ion channel assembles as a homopentamer, is exceptionally permeable to Ca²⁺, and desensitizes faster than any other Cys-loop receptor. The α7 receptor has served as a prototype for the Cys-loop superfamily yet has proven refractory to structural analysis. We present cryo-EM structures of the human α7 nicotinic receptor in a lipidic environment in resting, activated, and desensitized states, illuminating the principal steps in the gating cycle. The structures also reveal elements that contribute to its function, including a C-terminal latch that is permissive for channel opening, and an anionic ring in the extracellular vestibule that contributes to its high conductance and calcium permeability. Comparisons among the α7 structures provide a foundation for mapping the gating cycle and reveal divergence in gating mechanisms in the Cys-loop receptor superfamily.

A Computational Analysis of the Factors Governing the Dynamics of α7 nAChR and Its Homologs

The α7 nicotinic acetylcholine receptor is a homopentameric ion channel from the Cys-loop receptor superfamily targeted for psychiatric indications and inflammatory pain. Molecular dynamics studies of the receptor have focused on residue mobility and global conformational changes to address receptor function. However, a comparative analysis of α7 with its homologs that cannot trigger channel opening has not been made so far. To identify the residues involved in α7 activation, we ran triplicate 500-ns molecular dynamics simulations with an α7 extracellular domain homology model and two acetylcholine-binding protein homologs. We tested the effect of ligand binding and amino acid sequence on the structure and dynamics of the three proteins. We found that mobile regions identified based on root mean-square deviation and root mean-square fluctuation values are not always consistent among the individual α7 extracellular domain simulations. Comparison of the replica-average properties of the three proteins based on dynamic cross-correlation maps showed that ligand binding affects the coupling between the C-loop and the Cys-loop, vestibular loop, and β1-β2 loops. In addition, the main-immunogenic-region-like domain of α7 went through correlated motions with multiple domains of the receptor. These correlated motions were absent or diminished in α7 homologs, suggesting a unique role in α7 activation.

Nicotinic receptor pharmacology in silico: Insights and challenges

Nicotinic acetylcholine receptors (nAChR) are homo- or hetero-pentameric ligand-gated ion channels of the Cys-loop superfamily and play important roles in the nervous system and muscles. Studies on nAChR benefit from in silico modeling due to the lack of high-resolution structures for most receptor subtypes and challenges in experiments addressing the complex mechanism of activation involving allosteric sites. Although there is myriad of computational modeling studies on nAChR, the multitude of the methods and parameters used in these studies makes modeling nAChR a daunting task, particularly for the non-experts in the field. To address this problem, the modeling literature on Torpedo nAChR and α7 nAChR were focused on as examples of heteromeric and homomeric nAChR, and the key in silico modeling studies between the years 1995-2019 were concisely reviewed. This was followed by a critical analysis of these studies by comparing the findings with each other and with the emerging experimental and computational data on nAChR. Based on these critical analyses, suggestions were made to guide the future researchers in the field of in silico modeling of nAChR. This article is part of the special issue on 'Contemporary Advances in Nicotine Neuropharmacology'.

Molecular dynamics simulations of dihydro-β-erythroidine bound to the human α4β2 nicotinic acetylcholine receptor

BACKGROUND AND PURPOSE

The heteromeric α4β2 nicotinic acetylcholine receptor (nAChR) is abundant in the human brain and is associated with a range of CNS disorders. This nAChR subtype has been recently crystallised in a conformation that was proposed to represent a desensitised state. Here, we investigated the conformational transition mechanism of this nAChR from a desensitised to a closed/resting state.

EXPERIMENTAL APPROACH

The competitive antagonist dihydro-β-erythroidine (DHβE) was modelled by replacement of the agonist nicotine in the α4β2 nAChR experimental structure. DHβE is used both in vitro and in vivo for its ability to block α4β2 nAChRs. This system was studied by three molecular dynamics simulations with a combined simulation time of 2.6 μs. Electrophysiological studies of mutated receptors were performed to validate the simulation results.

KEY RESULTS

The relative positions of the extracellular and transmembrane domains in the models are distinct from those of the desensitised state structure and are compatible with experimental structures of Cys-loop receptors captured in a closed/resting state.

CONCLUSIONS AND IMPLICATIONS

Our model suggests that the side chains of α4 L257 (9') and α4 L264 (16') are the main constrictions in the transmembrane pore. The involvement of position 9' in channel gating is well established, but position 16' was only previously identified as a gate for the bacterial channels, ELIC and GLIC. L257 but not L264 was found to influence the slow component of desensitisation. The structure of the antagonist-bound state proposed here should be valuable for the development of therapeutic or insecticide compounds.