Molecular Modeling of Human Pentameric α7Neuronal Nicotinic Acetylcholine Receptor and Its Interaction with its Agonist and Competitive Antagonist

Prospective, observational, single-centre cohort study with an independent control group matched for age and sex aimed at investigating the significance of cholinergic activity in patients with schizophrenia: study protocol of the CLASH-study

INTRODUCTION

Alterations in the cholinergic metabolism may cause various clinical symptoms of schizophrenia. In addition to the 'monoamine hypothesis,' neuroinflammation is also discussed as a cause of schizophrenia. To date, there has been no evidence of alterations in the central cholinergic transmitter balance in patients with schizophrenia under clinical conditions. By contrast, studies in critically ill patients have established the measurement of acetylcholinesterase activity as a suitable surrogate parameter of central cholinergic transmitter balance/possible pathophysiological changes. Butyrylcholinesterase activity has been established as a parameter indicating possible (neuro)inflammatory processes. Both parameters can now be measured using a point-of-care approach. Therefore, the primary objective of this study is to investigate whether acetylcholinesterase and butyrylcholinesterase activity differs in patients with various forms of schizophrenia. Secondary objectives address the possible association between acetylcholinesterase and butyrylcholinesterase activity and (1) schizophrenic symptoms using the Positive and Negative Syndrome Scale, (2) the quantity of antipsychotics taken and (3) the duration of illness.

METHODS AND ANALYSIS

The study is designed as a prospective, observational cohort study with one independent control group. It is being carried out at the Department of Psychiatry and Psychotherapy III, Ulm University Hospital, Germany. Patient enrolment started in October 2020, and the anticipated end of the study is in January 2022. The enrolment period was set from October 2020 to December 2021 (extension required due to SARS-CoV-2 pandemic). The sample size is calculated at 50 patients in each group. Esterase activity is measured on hospital admission (acute symptomatology) and after referral to a postacute ward over a period of three consecutive days. The matched control group will be created after reaching 50 patients with schizophrenia. This will be followed by a comprehensive statistical analysis of the data set.

ETHICS AND DISSEMINATION

The study was registered prospectively in the German Clinical Trials Register (DRKS-ID: DRKS00023143,URL: https://www.drks.de/drks_web/navigate.do?navigationId=trial.HTML&TRIAL_ID=DRKS00023143) after approval by the ethics committee of the University of Ulm, Germany Trial Code No. 280/20.

TRIAL REGISTRATION NUMBER

DRKS00023143; Pre-results.

Critical Molecular Determinants of α7 Nicotinic Acetylcholine Receptor Allosteric Activation: SEPARATION OF DIRECT ALLOSTERIC ACTIVATION AND POSITIVE ALLOSTERIC MODULATION

The α7 nicotinic acetylcholine receptors (nAChRs) are uniquely sensitive to selective positive allosteric modulators (PAMs), which increase the efficiency of channel activation to a level greater than that of other nAChRs. Although PAMs must work in concert with "orthosteric" agonists, compounds such as GAT107 ((3aR,4S,9bS)-4-(4-bromophenyl)-3a,4,5,9b-tetrahydro-3H-cyclopenta[c]quinoline-8-sulfonamide) have the combined properties of agonists and PAMs (ago-PAM) and produce very effective channel activation (direct allosteric activation (DAA)) by operating at two distinct sites in the absence of added agonist. One site is likely to be the same transmembrane site where PAMs like PNU-120596 function. We show that the other site, required for direct activation, is likely to be solvent-accessible at the extracellular domain vestibule. We identify key attributes of molecules in this family that are able to act at the DAA site through variation at the aryl ring substituent of the tetrahydroquinoline ring system and with two different classes of competitive antagonists of DAA. Analyses of molecular features of effective allosteric agonists allow us to propose a binding model for the DAA site, featuring a largely non-polar pocket accessed from the extracellular vestibule with an important role for Asp-101. This hypothesis is supported with data from site-directed mutants. Future refinement of the model and the characterization of specific GAT107 analogs will allow us to define critical structural elements that can be mapped onto the receptor surface for an improved understanding of this novel way to target α7 nAChR therapeutically.

Looking below the surface of nicotinic acetylcholine receptors

The amino acid sequences of nicotinic acetylcholine receptors (nAChRs) from diverse species can be compared across extracellular, transmembrane, and intracellular domains. The intracellular domains are most divergent among subtypes, yet relatively consistent among species. The diversity indicates that each nAChR subtype has a unique language for communication with its host cell. The conservation across species also suggests that the intracellular domains have defining functional roles for each subtype. Secondary structure prediction indicates two relatively conserved alpha helices within the intracellular domains of all nAChRs. Among all subtypes, the intracellular domain of α7 nAChR is one of the most well conserved, and α7 nAChRs have effects in non-neuronal cells independent of generating ion currents, making it likely that the α7 intracellular domain directly mediates signal transduction. There are potential phosphorylation and protein-binding sites in the α7 intracellular domain, which are conserved and may be the basis for α7-mediated signal transduction.

Approaching the 5-HT₃ receptor heterogeneity by computational studies of the transmembrane and intracellular domains

5-hydroxytryptamine type-3 receptor (5-HT₃), an important target of many neuroactive drugs, is a cation selective transmembrane pentamer whose functional stoichiometries and subunit arrangements are still debated, due to the extreme complexity of the system. The three dimensional structure of the 5-HT₃R subunits has not been solved so far. Moreover, most of the available structural and functional data is related to the extracellular ligand-binding domain, whereas the transmembrane and the intracellular receptor domains are far less characterised, although they are crucial for receptor function. Here, for the first time, 3D homology models of the transmembrane and the intracellular receptor domains of all the known human 5-HT₃ subunits have been built and assembled into homopentameric (5-HT(3A)R, 5-HT(3B)R, 5-HT(3C)R, 5-HT(3D)R and 5-HT(3E)R) and heteropentameric receptors (5-HT(3AB), 5-HT(3AC), 5-HT(3AD) and 5-HT(3AE)), on the basis of the known three-dimensional structures of the nicotinic-acetylcholine receptor and of the ligand gated channel from Erwinia chrysanthemi. The comparative analyses of sequences, modelled structures, and computed electrostatic properties of the single subunits and of the assembled pentamers shed new light both on the stoichiometric composition and on the physicochemical requirements of the functional receptors. In particular, it emerges that a favourable environment for the crossing of the pore at the transmembrane and intracellular C terminus domain levels by Ca²⁺ ions is granted by the maximum presence of two B subunits in the 5-HT₃ pentamer.

Biological evaluation and interaction of a newly designed anti-cancer Pd(II) complex and human serum albumin

The pharmacokinetics and pharmacodynamics of any drug will depend, largely, on the interaction that has with human serum albumin (HSA), the most abundant plasma protein. The interaction between newly synthesized Pd(II) complexe, 2,2'-bipyridin Butylglycinato Pd(II) nitrate, an anti-tumor component, with HSA was studied at different temperatures by fluorescence, far UV circular dichroism (CD), UV-visible spectrophotometry and theoretical approaches. The Pd(II) complex has a strong ability to quench the intrinsic fluorescence of HSA through a dynamic quenching procedure. The binding parameters and thermodynamic parameters, including δH°, δS° and δG° were calculated by fluorescence quenching method, indicated that hydrophobic forces play a major role in the interaction of Pd(II) complex with HSA. Based on Autodock, FRET (fluorescence resonance energy transfer) and fluorescence quenching data, it may be concluded that one of the binding sites in the complex of HSA is near the only one Trp of HSA (Trp214) in sub domain IIA of the protein. Far-UV-CD results indicated that Pd(II)-complex induced increase in the α-helical content of the protein. The anti-tumor property of the synthesized Pd(II) complex was studied by testing it on human tumor cell line K562. The 50% cytotoxic concentration (Cc₅₀) of complex was determined using MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide) assay. Also, fluorescence staining with DAPI (4,6-diamidino-2-phenylindole) revealed some typical nuclear changes that are characteristic of apoptosis which is induced at Cc₅₀ concentration of Pd(II) complex in K562 cell line after 24 h incubation. Our results suggest that Pd(II) complex is a promising anti-proliferative agent and should execute its biological effects by inducing apoptosis.

Molecular dynamics and docking studies on cardiac troponin C

Cardiac troponin C (cTnC) is the Ca²⁺ dependent switch for contraction in heart muscle making it a potential target for drug research in the therapy of heart failure. Calcium binding on Troponin C (TnC) triggers a series of conformational changes exposing a hydrophobic pocket in the N-domain of TnC (cNTnC), which leads to force generation. Mutations and acidic pH have been related to altering the sensitivity of TnC affecting the efficiency of the heart. Bepridil, identified as a calcium sensitizer to TnC, has been experimentally found to bind to the N-domain pocket of TnC but with negative cooperativity. Screening and de novo design were carried out using LUDI and AUTOLUDI programs in this work to identify and design potential ligands that can bind to the hydrophobic pocket of TnC. Two docking centers and multiple searching radii including 5 Å, 5.5 Å, 6 Å, 6.5 Å, 7.0 Å and 7.5 Å were used in LUDI to screen the ZINC database. Based on the LUDI docking results, 8 molecules were identified from the database with good potential to bind into the binding pocket and they were used as template molecules to generate a series of new molecules by AUTOLUDI design. Out of all the newly-designed molecules, 14 new ligands were recognized to be potential ligands that can bind and fit well into the binding pocket. These molecules can be used as starting molecules to develop TnC ligands. The binding stability and binding affinity of these molecules to the protein was further analyzed by molecular dynamics simulations. The results show that the binding energies, interactions and complex stabilities of 6 ligands are comparable to or better than bepridil.

Biphasic association of T7 RNA polymerase and a nucleotide analogue, cibacron blue as a model to understand the role of initiating nucleotide in the mechanism of enzyme action

T7 RNA polymerase (T7 RNAP) is an enzyme that utilizes ribonucleotides to synthesize the nascent RNA chain in a template-dependent manner. Here we have studied the interaction of T7 RNAP with cibacron blue, an anthraquinone monochlorotriazine dye, its effect on the function of the enzyme and the probable mode of binding of the dye. We have used difference absorption spectroscopy and isothermal titration calorimetry to show that the dye binds T7 RNAP in a biphasic manner. The first phase of the binding is characterized by inactivation of the enzyme. The second binding site overlaps with the common substrate-binding site of the enzyme. We have carried out docking experiment to map the binding site of the dye in the promoter bound protein. Competitive displacement of the dye from the high affinity site by labeled GTP and isothermal titration calorimetry of high affinity GTP bound enzyme with the dye suggests a strong correlation between the high affinity dye binding and the high affinity GTP binding in T7 RNAP reported earlier from our laboratory.

In search of allosteric modulators of a7-nAChR by solvent density guided virtual screening

Nicotinic acetylcholine receptors (nAChR) are pentameric ligand gated ion channels whose activity can be modulated by endogenous neurotransmitters as well as by synthetic ligands that bind the same or distinct sites from the natural ligand. The subtype of α7 nAChR has been considered as a potenial therapeutic target for Alzheimer's disease, schizophrenia and other neurological and psychiatric disorders. Here we have developed a homology model of α7 nAChR based on two high resolution crystal structures with Brookhaven Protein Data Bank (PDB) codes 2QC1 and 2WN9 for threading on one monomer and then for building a pentamer, respectively. A number of small molecule binding sites are identified using Pocket Finder (J. An, M. Tortov, and R. Abagyan, Molecular & Cellular Proteomics, 4.6, 752-761 (2005)) of Internal Coordinate Mechanics (ICM). Remarkably, these computer-identified sites match perfectly with ordered solvent densities found in the high-resolution crystal structure of α1 nAChR, suggesting that the surface cavities in the α7 nAChR model are likely binding sites of small molecules. A high throughput virtual screening by flexible ligand docking of 5008 small molecule compounds was performed at three potential allosteric modulator (AM) binding sites of α7 nAChR using Molsoft ICM software (R. Abagyan, M. Tortov and D. Kuznetsov, J Comput Chem 15, 488-506, (1994)). Some experimentally verified allosteric modulators of α7 like CCMI comp-6, LY 7082101, 5-HI, TQS, PNU-120596, genistein, and NS-1738 ranked among top 100 compounds, while the rest of the compounds in the list could guide further search for new allosteric modulators.

Screening from the world's largest TCM database against H1N1 virus

The swine influenza virus (H1N1) 2009 pandemic highlights the importance of having effective anti-viral strategies. Recently, oseltamivir (Tamiflu) resistant influenza viruses are identified; which further emphasizes the urgency in developing new antiviral agents. In influenza virus replication cycle, viral surface glycoprotein, hemagglutinin, is responsible for viral entry into host cells. Hence, a potentially effective antiviral strategy is to inhibit viral entry mechanism. To develop novel antiviral agent that inhibits viral entry, we analyzed 20,000 traditional Chinese medicine (TCM) ingredients in hemagglutinin subtype H1 sialic acid binding site found on H1N1 virus. We then performed molecular dynamics simulations to investigate receptor-ligand interaction of the candidates obtained from docking. Here, we report three TCM derivatives that have high binding affinities to H1 sialic acid binding site residues based on structure-based calculations. The top three derivatives, xylopine_2, rosmaricine_14 and rosmaricine_15, all have an amine group that interact with Glu83 and a pyridinium group that interact with Asp103. Molecular dynamics simulations show that these derivatives form strong hydrogen bonding with Glu83 but interact transiently with Asp103. We therefore suggest that an enhanced hemagglutinin inhibitor, based on our scaffold, should be designed to bind both Glu83 and Asp103 with high affinity.

Protein-protein interactions: principles, techniques, and their potential role in new drug development

A vast network of genes is inter-linked through protein-protein interactions and is critical component of almost every biological process under physiological conditions. Any disruption of the biologically essential network leads to pathological conditions resulting into related diseases. Therefore, proper understanding of biological functions warrants a comprehensive knowledge of protein-protein interactions and the molecular mechanisms that govern such processes. The importance of protein-protein interaction process is highlighted by the fact that a number of powerful techniques/methods have been developed to understand how such interactions take place under various physiological and pathological conditions. Many of the key protein-protein interactions are known to participate in disease-associated signaling pathways, and represent novel targets for therapeutic intervention. Thus, controlling protein-protein interactions offers a rich dividend for the discovery of new drug targets. Availability of various tools to study and the knowledge of human genome have put us in a unique position to understand highly complex biological network, and the mechanisms involved therein. In this review article, we have summarized protein-protein interaction networks, techniques/methods of their binding/kinetic parameters, and the role of these interactions in the development of potential tools for drug designing.

The key residues of active sites on the catalytic fragment for paclitaxel interacting with poly (ADP-ribose) polymerase

Poly(ADP-ribose) polymerase (PARP) is regarded as a target protein for paclitaxel (PTX) to bind. An important issue is to identify the key residues as active sites for PTX interacting with PARP, which will help to understand the potential drug activity of PTX against cancer cells. Using docking method and MD simulation, we have constructed a refined structure of PTX docked on the catalytic function domain of PARP (PDB code: 1A26). The residues Glu327(988), Tyr246(907), Lys242(903), His165(826), Asp105(766), Gln102(763) and Gln98(759) in PARP are identified as potential sites involved in interaction with PTX according to binding energy (E(b)) between PTX and single residue calculated with B3LYP/6-31G(d,p). These residues form an active binding pocket located on the surface of the catalytic fragment, possibly interacting with the required groups of PTX leading to its activity against cancer cells. It is noted that most of the active sites make conatct with the "southern hemisphere" of PTX except for one residue, Tyr246(907), which interacts with the "northern hemisphere" of PTX. The conformation of PTX in complex with the catalytic fragment is observed as being T-shaped, similar to that complexed with β-tubulin. The total Eb of -269.9 kJ/mol represents the potent interaction between PTX and the catalytic fragment, implying that PTX can readily bind to the active pocket. The tight association of PTX with the catalytic fragment would inhibit PARP activation, suggesting a potential application of PTX as an effective antineoplastic agent.

Discrimination of agonists versus antagonists of nicotinic ligands based on docking onto AChBP structures

Numerous high-resolution crystallographic structures of the acetylcholine binding protein (AChBP), a molluscan cholinergic protein, homologous to the extracellular domain of nicotinic acetylcholine receptors, are available. This offers opportunities to model the interaction between various ligands and the acetylcholine binding site. Herein we present a study of the interplay between ligand binding and motions of the C-loop capping the binding site. Nicotinic agonists and antagonists were docked on AChBP X-ray structures. It is shown that the studied agonists and antagonists can be discriminated according to their higher affinities for structures respectively obtained in the presence of agonists or antagonists, highlighting the fact that AChBP structures retain a pharmacological footprint of the compound used in crystallography experiments. A detailed analysis of the binding site cavities suggests that this property is mainly related to the shape of the cavities.

A novel OmpY porin from Yersinia pseudotuberculosis: structure, channel-forming activity and trimer thermal stability

A novel OmpY porin was predicted based on the Yersinia pseudotuberculosis genome analysis. Whereas it has the different genomic annotation such as "outer membrane protein N" (ABS46310.1) in str. IP 31758 or "outer membrane protein C2, porin" (YP_070481.1) in str. IP32953, it might be warranted to rename the OmpN/OmpC2 to OmpY, "outer membrane protein Y", where letter "Y" pertained to Yersinia. Both phylogenetic analysis and genomic localization clearly support that the OmpY porin belongs to a new group of general bacterial porins. The recombinant OmpY protein with its signal sequence was overexpressed in porin-deficient Escherichia coli strain. The mature rOmpY was shown to insert into outer membrane as a trimer. The OmpY porin, isolated from the outer membrane, was studied employing spectroscopic, electrophoretic and bilayer lipid membranes techniques. The far UV CD spectrum of rOmpY was essentially identical to that of Y. pseudotuberculosis OmpF. The near UV CD spectrum of rOmpY was weaker and smoother than that of OmpF. The rOmpY single-channel conductance was 180 ± 20 pS in 0.1 M NaCl and was lower than that of the OmpF porin. As was shown by electrophoretic and bilayer lipid membrane experiments, the rOmpY trimers were less thermostable than the OmpF trimers. The porins differed in the trimer-monomer transition temperature by about 20°C. The three-dimensional structural models of the Y. pseudotuberculosis OmpY and OmpF trimers were generated and the intra- and intermonomeric interactions stabilizing the porins were investigated. The difference in the thermal stability of OmpY and OmpF trimers was established to correlate with the difference in intermonomeric polar contacts.

SMpred: a support vector machine approach to identify structural motifs in protein structure without using evolutionary information

Knowledge of three dimensional structure is essential to understand the function of a protein. Although the overall fold is made from the whole details of its sequence, a small group of residues, often called as structural motifs, play a crucial role in determining the protein fold and its stability. Identification of such structural motifs requires sufficient number of sequence and structural homologs to define conservation and evolutionary information. Unfortunately, there are many structures in the protein structure databases have no homologous structures or sequences. In this work, we report an SVM method, SMpred, to identify structural motifs from single protein structure without using sequence and structural homologs. SMpred method was trained and tested using 132 proteins domains containing 581 motifs. SMpred method achieved 78.79% accuracy with 79.06% sensitivity and 78.53% specificity. The performance of SMpred was evaluated with MegaMotifBase using 188 proteins containing 1161 motifs. Out of 1161 motifs, SMpred correctly identified 1503 structural motifs reported in MegaMotifBase. Further, we showed that SMpred is useful approach for the length deviant superfamilies and single member superfamilies. This result suggests the usefulness of our approach for facilitating the identification of structural motifs in protein structure in the absence of sequence and structural homologs. The dataset and executable for the SMpred algorithm is available at http://www3.ntu.edu.sg/home/EPNSugan/index_files/SMpred.htm.

Homology modeling of Cry10Aa toxin from B. thuringiensis israelensis and B. thuringiensis subsp. LDC-9

A three dimensional model was developed for Cry10Aa protein sequence of B. thuringiensis LDC-9 and B. thuringiensis israelensis that has not been solved empirically by X-ray crystallography or NMR. Homology modeling was employed for the structure prediction using Cry2Aa as template protein, a high-resolution X-ray crystallography structure. The model predicted for the B. thuringiensis LDC-9 Cry10Aa protein reveals a partial N-terminal domain only due to its partial sequence of 104 amino acids. B. thuringiensis israelensis Cry10Aa model contains three domains such as domain I, a bundle of eight alpha helices with the central relatively hydrophobic helix surrounded by amphipathic helices while domain II and III contain mostly beta-sheets. Significant structural differences within domain II in this model among all Cry protein structures indicates that it is involved in recognition and binding to cell surfaces. Comparison of B. thuringiensis israelensis predicted structure with available experimentally determined Cry structures reveals identical folds. The distribution of electrostatic potential on the surface of the molecules in the model is non-uniform and identifies one side of the alpha-helical domain as negatively charged indicating orientation of toxic molecules toward the cell membrane during the initial binding with a cell surface receptor. The collective knowledge of Cry toxin structures will lead to a more critical understanding of the structural basis for receptor binding and pore formation, as well as allowing the scope of diversity to be better appreciated. This model will serve as a starting point for the design of mutagenesis experiments aimed to improve the toxicity and to provide a new tool for the elucidation of the mechanism of action of these mosquitocidal proteins.

Anti-infectious bronchitis virus (IBV) activity of 1,8-cineole: effect on nucleocapsid (N) protein

In the present study, anti-IBV (infectious bronchitis virus) activity of 1,8-cineole was studied by MTT assay, as well as docking and molecular dynamic (MD) simulations. The CC50 of 1,8-cineole was above 10 mM. And the maximum noncytotoxic concentration (TD0) of 1,8-cineole was determined to be 3.90 ± 0.22 mM, which was much higher than that of ribavirin (0.78 ± 0.15 mM). 1,8-cineole could inhibit IBV with an IC(50) of 0.61 mM. MTT assay showed that the inhibition of IBV by 1, 8-cineole appears to occur moderately before entering the cell but much strongly after penetration of the virus into the cell. In silico simulations indicated that the binding site of 1,8-cineole was located at the N terminus of phosphorylated nucleocapsid (N) protein, with interaction energy equaling -40.33 kcal mol(-1). The residues TyrA92, ProA134, PheA137, AspA138 and TyrA140 had important roles during the binding process and are fully or partially conserved in various IBV strains. Based on spatial and energetic criteria, 1,8-cineole interfered with the binding between RNA and IBV N-protein. Results presented here may suggest that 1,8-cineole possesses anti-IBV properties, and therefore is a potential source of anti-IBV ingredients for the pharmaceutical industry.

In silico point mutation and evolutionary trace analysis applied to nicotinic acetylcholine receptors in deciphering ligand-binding surfaces

The nicotinic acetylcholine receptors (nAChRs) are members of the Cys-loop superfamily and contain ligand gated ion channels (LGIC). These receptors are located mostly in the central nervous system (CNS) and peripheral nervous system (PNS). nAChRs reside at pre-synaptic regions to mediate acetylcholine neurotransmission and in the post synaptic membrane to propagate nerve impulses through neurons via acetylcholine. Malfunction of this neurotransmitter receptor is believed to cause various neurodegenerative diseases, such as Alzheimer's disease, Parkinson's disease and schizophrenia, and nAChRs are thus important drug targets. In the present work, starting from an earlier model of pentameric alpha7nAChR, a considerable effort has been taken to investigate interaction with ligands by performing docking studies with a diverse array of agonists and antagonists. Analysis of these docking complexes reveals identification of possible ligand-interacting residues. Some of these residues, e.g. Ser34, Gln55, Ser146, and Tyr166, which are evolutionarily conserved, were specifically subjected to virtual mutations based on their amino acid properties and found to be highly sensitive in the presence of antagonists by docking. Further, the study was extended using evolutionary trace analysis, revealing conserved and class-specific residues close to the putative ligand-binding site, further supporting the results of docking experiments.