The Psychonauts’ Benzodiazepines; Quantitative Structure-Activity Relationship (QSAR) Analysis and Docking Prediction of Their Biological Activity

© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).Designer benzodiazepines (DBZDs) represent a serious health concern and are increasingly re-ported in polydrug consumption-related fatalities. When new DBZDs are identified, very limited information is available on their pharmacodynamics. Here, computational models (e.g., quantita-tive structure-activity relationship/QSAR and Molecular Docking) were used to analyse DBZDs identified online by an automated web crawler (NPSfinder®) and to predict their possible activi-ty/affinity on the gamma-aminobutyric acid receptors (GABA-ARs). The computational software MOE was used to calculate 2D QSAR models, perform docking studies on crystallised GABA-A receptors (6HUO, 6HUP) and generate pharmacophore queries from the docking conformational results. 101 DBZDs were identified online by NPSfinder®. The validated QSAR model predicted high biological activity values for 41% of these DBDZs. These predictions were supported by the docking studies (good binding affinity) and the pharmacophore modelling confirmed the im-portance of the presence and location of hydrophobic and polar functions identified by QSAR. This study confirms once again the importance of web-based analysis in the assessment of drug scenarios (DBZDs), and how computational models could be used to acquire fast and reliable in-formation on biological activity for index novel DBZDs, as preliminary data for further investiga-tions.Peer reviewe

Modellierung von Metalloenzymen: 3D-QSAR-Untersuchungen an Carboanhydrase-Isoenzymen und virtuelles Screening nach Peptiddeformylase-Inhibitoren

In der modernen Arzneistoffentwicklung unterscheidet man die Phasen der Leitstrukturfindung und der Leitstrukturoptimierung. Die vorliegende Dissertationsschrift beinhaltet Beiträge zu beiden Bereichen. Der erste Teil der Arbeit befasst sich mit der Entwicklung und Evaluierung von Computermodellen zur Vorhersage von Affinität und Selektivität und entstammt daher dem Bereich der Leitstrukturoptimierung. Selektivitätsaspekte spielen eine wichtige Rolle, da sie das Risiko von Nebenwirkungen maßgeblich beeinflussen. Zur Modellierung und Vorhersage von Affinitäts- und Selektivitätsparametern wurden QSAR-Methoden angewendet. Das Modellsystem stellten Carboanhydrasen (CAs) dar; diese zinkhaltigen Hydrolasen katalysieren die reversible Hydratisierung von Kohlendioxid zu Bicarbonat und einem Proton. Sie sind daher in eine Vielzahl (patho)physiologischer Prozesse involviert und stellen interessante therapeutische Targets dar. Die zahlreichen Isoenzyme der CAs besitzen im Bereich der Bindetasche hohe Ähnlichkeiten in Bezug auf physikochemische Eigenschaften, so dass die Entwicklung selektiver Inhibitoren kein triviales Problem darstellt. Im Mittelpunkt der Untersuchungen standen insbesondere 3D-QSAR-Verfahren. Es wurden statistisch hochsignifikante und robuste Modelle abgeleitet, um Affinität und Selektivität von Sulfonamidinhibitoren bezüglich der Isoenzyme CA I, II und IV vorherzusagen. Es zeigte sich, dass die geringen Unterschiede im strukturbasierten Alignment unter Verwendung der drei Bindetaschen nur geringen Einfluss auf die statistischen Parameter besitzen und dass bessere Ergebnisse erzielt werden, wenn für alle Isoenzyme das auf CA II basierende Alignment benutzt wird anstelle des Alignments in der jeweiligen Bindetasche. Ursache hierfür ist wahrscheinlich die Vielzahl an Kristallstrukturen, die für CA II existieren und damit das Alignment verlässlicher machen. Die erhaltenen Isokonturkarten erlaubten eine Interpretation der Modelle im Hinblick auf die Bedeutung physikochemischer Eigenschaften für die Affinität/Selektivität. Der Vergleich zu qualitativen proteinbasierten Isokonturkarten unterstreicht den komplementären Charakter beider Methoden: Während die ligandbasierten QSAR-Verfahren implizit teilweise die Struktur der Bindetasche widerspiegeln, aber auch von den Eigenheiten des Trainingsdatensatzes abhängen, vermögen die proteinbasierten Analysen auch Informationen über Bereiche der Bindetasche zu geben, die keine Interaktionen mit Liganden des Trainingsdatensatzes ausbilden. Ein weiteres Ziel bestand darin, QSAR-Methoden für das Screening größerer Datenbanken zu verwenden. Dies erlaubt die Identifizierung besonders interessanter (d.h. affiner/selektiver) Kandidaten zur Synthese im Sinne einer Leitstrukturoptimierung. Für 3D-QSAR-Methoden musste zunächst ein Protokoll zur Automatisierung des Alignments entwickelt und validiert werden. Es zeigte sich hierbei, dass ein ligandbasiertes Alignment vergleichbare Ergebnisse zu manuellen stukturbasierten Alignmentmethoden erzielt. Die 3D-Modelle erwiesen sich als überlegen im Vergleich zu fragmentbasierten 2D-Methoden oder insbesondere zu den eigenschaftsbasierten 1D-Methoden. Als praktisches Anwendungsbeispiel der entwickelten Modelle wurde eine mehrere tausend Einträge umfassende virtuelle Ligandbibliothek aufgebaut und mit den leistungsfähigsten Modellen bewertet. Der zweite Teil der Arbeit beinhaltet ein virtuelles Screening nach neuartigen Inhibitoren der Peptiddeformylasen (PDFs) und gehört somit in den Bereich der Leitstrukturfindung. PDFs sind (meist) eisenhaltige Enzyme, die die Deformylierung von in Mitochondrien, Plastiden oder Bakterien synthetisierten Proteinen katalysieren. Ausgehend von Kristallstrukturen potenter PDF-Inhibitoren wurden 3D-Pharmakophormodelle entwickelt und validiert. Diese waren in der Lage, strukturell diverse, aus der Literatur bekannte Inhibitoren zu identifizieren (hinreichende Sensitivität) und gleichzeitig die zu durchsuchenden Datenbanken stark zu reduzieren (hinreichende Spezifität). Die Pharmakophormodelle wurden zum Screening von Datenbanken kommerziell erhältlicher Moleküle mit wirkstoffartigen Eigenschaften benutzt. Durch Docking und Scoring wurden schließlich aus etwa zwei Millionen Verbindungen elf Substanzen identifiziert und erworben, die einer biologischen Testung unterzogen werden sollen. Erste vorliegende Messergebnisse zeigen, dass mindestens zwei der Substanzen mit einem IC50-Wert von 60 nM bzw. 190 nM potente Inhibitoren der PDF1B aus E. coli sind. Dies belegt die Güte der Modelle und des angewendeten Screening-Protokolls. Inhibitoren der PDF könnten Anwendung als Herbizide, Antibiotika und Antimalaria-Therapeutika finden

IN SILICO METHODS FOR DRUG DESIGN AND DISCOVERY

Computer-aided drug design (CADD) methodologies are playing an ever-increasing role in drug discovery that are critical in the cost-effective identification of promising drug candidates. These computational methods are relevant in limiting the use of animal models in pharmacological research, for aiding the rational design of novel and safe drug candidates, and for repositioning marketed drugs, supporting medicinal chemists and pharmacologists during the drug discovery trajectory.Within this field of research, we launched a Research Topic in Frontiers in Chemistry in March 2019 entitled “In silico Methods for Drug Design and Discovery,” which involved two sections of the journal: Medicinal and Pharmaceutical Chemistry and Theoretical and Computational Chemistry. For the reasons mentioned, this Research Topic attracted the attention of scientists and received a large number of submitted manuscripts. Among them 27 Original Research articles, five Review articles, and two Perspective articles have been published within the Research Topic. The Original Research articles cover most of the topics in CADD, reporting advanced in silico methods in drug discovery, while the Review articles offer a point of view of some computer-driven techniques applied to drug research. Finally, the Perspective articles provide a vision of specific computational approaches with an outlook in the modern era of CADD

Using fuzzy methods for rule extraction in the discrimination of class C GPCR subtypes from their subsequences

G-Protein-Coupled receptors (GPCR) are cell membrane proteins that regulate many of the cell functions and transduce signals between the intracellular and extracellular domains. This makes them relevant in pharmacology as therapeutic targets. As members of this superfamily, class C GPCRs in particular regulate a number of important physiological functions. Proteins of the class must be studied from their primary sequences, as only one of their 3-D structures has been fully determined, earlier this year. Protein function investigation requires the identification of motifs, or functional subsequences. In this thesis, we will describe the discrimination of class C GPCR subtypes through interpretable rules from a specific alignment free transformation of the sequences, namely amino acid composition. The Fuzzy Inductive Reasoning methodology was used as the basis to extract these linguistic rules

Antimicrobial Drug Repurposing Through Molecular Modelling: Acquisition, Analyis and Prediction

Antimicrobial resistance has sparked unprecedented medical crises around the world, not only increasing the mortality rate but also impacting nosocomial resources. Methicillin-resistant Staphylococcus aureus (MRSA) has consistently evaded the available range of antibiotics and is a typical case study for new generation drugs. Drug development has been conventionally suffering from exceedingly high costs and overdrawn timelines. Drug Repurposing can be a solution to alleviate those burdens. Put simply, DR is a mechanism to identify new usages of existing drugs, typically targeted to treat diseases different to the ones that these were initially intended for. This inherently interdisciplinary research targets to identify the best MRSA drug candidates analysing protein (BIG) data, in the process developing a combination of techniques from stochastic mathematics, statistics and data analytics that can generically identify drug targets from the databank. Structure-based virtual screening was used to repurpose an extensive range of marketed drugs and Phase I/II/III trials. Molecular docking methods were used for virtual screening against MRSA targets based on sequence alignment to match gene sequences against proteins in the Protein Data Bank (PDB). Ligands from the Database of Useful Decoys - Enhanced were docked against MRSA-oriented target proteins using 10 open-source docking programmes for benchmark. The novel consensus scoring methods prove superior to other reported consensus scores in terms of discrimination between decoys and active ligands concerning MRSA drug target identification. The consensus scoring predictions are then applied to docking data between MRSA targets and compounds from the Repurposing Hub to identify a list of potential drug candidates for anti-MRSA treatment. MRSA is currently an apocalypse across the world with limited prevention and medications. This study provided more potential candidates to help fight against MRSA. The consensus scoring developed in this study can be generically implemented to unlock other antimicrobial drug candidates