Copper-Catalyzed One-Pot Synthesis of Thiazolidin-2-imines.

The synthesis of thiazolines, thiazolidines, and thiazolidinones has been extensively studied, due to their biological activity related to neurodegenerative diseases, such as Parkinson's and Alzheimer's, as well as their antiparasitic and antihypertensive properties. The closely related thiazolidin-2-imines have been studied less, and efficient strategies for synthesizing them, mainly based on the reaction of propargylamines with isothiocyanates, have been explored less. The use of one-pot approaches, providing modular, straightforward, and sustainable access to these compounds, has also received very little attention. Herein, we report a novel, one-pot, multicomponent, copper-catalyzed reaction among primary amines, ketones, terminal alkynes, and isothiocyanates, toward thiazolidin-2-imines bearing quaternary carbon centers on the five-membered ring, in good to excellent yields. Density functional theory calculations, combined with experimental mechanistic findings, suggest that the copper(I)-catalyzed reaction between the in situ-formed propargylamines and isothiocyanates proceeds with a lower energy barrier in the pathway leading to the S-cyclized product, compared to that of the N-cyclized one, toward the chemo- and regioselective formation of 5-exo-dig S-cyclized thiazolidin-2-imines.

Computational Workflow for Refining AlphaFold Models in Drug Design Using Kinetic and Thermodynamic Binding Calculations: A Case Study for the Unresolved Inactive Human Adenosine A3 Receptor.

A structure-based drug design pipeline that considers both thermodynamic and kinetic binding data of ligands against a receptor will enable the computational design of improved drug molecules. For unresolved GPCR-ligand complexes, a workflow that can apply both thermodynamic and kinetic binding data in combination with alpha-fold (AF)-derived or other homology models and experimentally resolved binding modes of relevant ligands in GPCR-homologs needs to be tested. Here, as test case, we studied a congeneric set of ligands that bind to a structurally unresolved G protein-coupled receptor (GPCR), the inactive human adenosine A3 receptor (hA3R). We tested three available homology models from which two have been generated from experimental structures of hA1R or hA2AR and one model was a multistate alphafold 2 (AF2)-derived model. We applied alchemical calculations with thermodynamic integration coupled with molecular dynamics (TI/MD) simulations to calculate the experimental relative binding free energies and residence time (τ)-random accelerated MD (τ-RAMD) simulations to calculate the relative residence times (RTs) for antagonists. While the TI/MD calculations produced, for the three homology models, good Pearson correlation coefficients, correspondingly, r = 0.74, 0.62, and 0.67 and mean unsigned error (mue) values of 0.94, 1.31, and 0.81 kcal mol-1, the τ-RAMD method showed r = 0.92 and 0.52 for the first two models but failed to produce accurate results for the multistate AF2-derived model. With subsequent optimization of the AF2-derived model by reorientation of the side chain of R1735.34 located in the extracellular loop 2 (EL2) that blocked ligand's unbinding, the computational model showed r = 0.84 for kinetic data and improved performance for thermodynamic data (r = 0.81, mue = 0.56 kcal mol-1). Overall, after refining the multistate AF2 model with physics-based tools, we were able to show a strong correlation between predicted and experimental ligand relative residence times and affinities, achieving a level of accuracy comparable to an experimental structure. The computational workflow used can be applied to other receptors, helping to rank candidate drugs in a congeneric series and enabling the prioritization of leads with stronger binding affinities and longer residence times.

Accurate calculation of affinity changes to the close state of influenza A M2 transmembrane domain in response to subtle structural changes of adamantyl amines using free energy perturbation methods in different lipid bilayers.

Experimental binding free energies of 27 adamantyl amines against the influenza M2(22-46) WT tetramer, in its closed form at pH 8, were measured by ITC in DPC micelles. The measured Kd's range is ~44 while the antiviral potencies (IC50) range is ~750 with a good correlation between binding free energies computed with Kd and IC50 values (r = 0.76). We explored with MD simulations (ff19sb, CHARMM36m) the binding profile of complexes with strong, moderate and weak binders embedded in DMPC, DPPC, POPC or a viral mimetic membrane and using different experimental starting structures of M2. To predict accurately differences in binding free energy in response to subtle changes in the structure of the ligands, we performed 18 alchemical perturbative single topology FEP/MD NPT simulations (OPLS2005) using the BAR estimator (Desmond software) and 20 dual topology calculations TI/MD NVT simulations (ff19sb) using the MBAR estimator (Amber software) for adamantyl amines in complex with M2(22-46) WT in DMPC, DPPC, POPC. We observed that both methods with all lipids show a very good correlation between the experimental and calculated relative binding free energies (r = 0.77-0.87, mue = 0.36-0.92 kcal mol-1) with the highest performance achieved with TI/MBAR and lowest performance with FEP/BAR in DMPC bilayers. When antiviral potencies are used instead of the Kd values for computing the experimental binding free energies we obtained also good performance with both FEP/BAR (r = 0.83, mue = 0.75 kcal mol-1) and TI/MBAR (r = 0.69, mue = 0.77 kcal mol-1).

Structure-Based Lead Optimization of Enterovirus D68 2A Protease Inhibitors.

Enterovirus D68 (EV-D68) virus is a nonpolio enterovirus that typically causes respiratory illness and, in severe cases, can lead to paralysis and death in children. There is currently no vaccine or antiviral for EV-D68. We previously discovered the viral 2A protease (2Apro) as a viable antiviral drug target and identified telaprevir as a 2Apro inhibitor. 2Apro is a viral cysteine protease that cleaves the viral VP1-2A polyprotein junction. In this study, we report the X-ray crystal structures of EV-D68 2Apro, wild-type, and the C107A mutant and the structure-based lead optimization of telaprevir. Guided by the X-ray crystal structure, we predicted the binding pose of telaprevir in 2Apro using molecular dynamics simulations. We then utilized this model to inform structure-based optimization of the telaprevir's reactive warhead and P1-P4 substitutions. These efforts led to the discovery of 2Apro inhibitors with improved antiviral activity than telaprevir. These compounds represent promising lead compounds for further development as EV-D68 antivirals.

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory.

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.

A Study of the Activity of Adamantyl Amines against Mutant Influenza A M2 Channels Identified a Polycyclic Cage Amine Triple Blocker, Explored by Molecular Dynamics Simulations and Solid-State NMR.

We compared the anti-influenza potencies of 57 adamantyl amines and analogs against influenza A virus with serine-31 M2 proton channel, usually termed as WT M2 channel, which is amantadine sensitive. We also tested a subset of these compounds against viruses with the amantadine-resistant L26F, V27A, A30T, G34E M2 mutant channels. Four compounds inhibited WT M2 virus in vitro with mid-nanomolar potency, with 27 compounds showing sub-micromolar to low micromolar potency. Several compounds inhibited L26F M2 virus in vitro with sub-micromolar to low micromolar potency, but only three compounds blocked L26F M2-mediated proton current as determined by electrophysiology (EP). One compound was found to be a triple blocker of WT, L26F, V27A M2 channels by EP assays, but did not inhibit V27A M2 virus in vitro, and one compound inhibited WT, L26F, V27A M2 in vitro without blocking V27A M2 channel. One compound blocked only L26F M2 channel by EP, but did not inhibit virus replication. The triple blocker compound is as long as rimantadine, but could bind and block V27A M2 channel due to its larger girth as revealed by molecular dynamics simulations, while MAS NMR informed on the interaction of the compound with M2(18-60) WT or L26F or V27A.

Molecular Biophysics of Class A G Protein Coupled Receptors-Lipids Interactome at a Glance-Highlights from the A2A Adenosine Receptor.

G protein-coupled receptors (GPCRs) are embedded in phospholipid membrane bilayers with cholesterol representing 34% of the total lipid content in mammalian plasma membranes. Membrane lipids interact with GPCRs structures and modulate their function and drug-stimulated signaling through conformational selection. It has been shown that anionic phospholipids form strong interactions between positively charged residues in the G protein and the TM5-TM6-TM 7 cytoplasmic interface of class A GPCRs stabilizing the signaling GPCR-G complex. Cholesterol with a high content in plasma membranes can be identified in more specific sites in the transmembrane region of GPCRs, such as the Cholesterol Consensus Motif (CCM) and Cholesterol Recognition Amino Acid Consensus (CRAC) motifs and other receptor dependent and receptor state dependent sites. Experimental biophysical methods, atomistic (AA) MD simulations and coarse-grained (CG) molecular dynamics simulations have been applied to investigate these interactions. We emphasized here the impact of phosphatidyl inositol-4,5-bisphosphate (PtdIns(4,5)P2 or PIP2), a minor phospholipid component and of cholesterol on the function-related conformational equilibria of the human A2A adenosine receptor (A2AR), a representative receptor in class A GPCR. Several GPCRs of class A interacted with PIP2 and cholesterol and in many cases the mechanism of the modulation of their function remains unknown. This review provides a helpful comprehensive overview for biophysics that enter the field of GPCRs-lipid systems.

Study of SQ109 analogs binding to mycobacterium MmpL3 transporter using MD simulations and alchemical relative binding free energy calculations.

N-geranyl-N΄-(2-adamantyl)ethane-1,2-diamine (SQ109) is a tuberculosis drug that has high potency against Mycobacterium tuberculosis (Mtb) and may function by blocking cell wall biosynthesis. After the crystal structure of MmpL3 from Mycobacterium smegmatis in complex with SQ109 became available, it was suggested that SQ109 inhibits Mmpl3 mycolic acid transporter. Here, we showed using molecular dynamics (MD) simulations that the binding profile of nine SQ109 analogs with inhibitory potency against Mtb and alkyl or aryl adducts at C-2 or C-1 adamantyl carbon to MmpL3 was consistent with the X-ray structure of MmpL3 - SQ109 complex. We showed that rotation of SQ109 around carbon-carbon bond in the monoprotonated ethylenediamine unit favors two gauche conformations as minima in water and lipophilic solvent using DFT calculations as well as inside the transporter's binding area using MD simulations. The binding assays in micelles suggested that the binding affinity of the SQ109 analogs was increased for the larger, more hydrophobic adducts, which was consistent with our results from MD simulations of the SQ109 analogues suggesting that sizeable C-2 adamantyl adducts of SQ109 can fill a lipophilic region between Y257, Y646, F260 and F649 in MmpL3. This was confirmed quantitatively by our calculations of the relative binding free energies using the thermodynamic integration coupled with MD simulations method with a mean assigned error of 0.74 kcal mol-1 compared to the experimental values.

Comparative Study of Receptor-, Receptor State-, and Membrane-Dependent Cholesterol Binding Sites in A2A and A1 Adenosine Receptors Using Coarse-Grained Molecular Dynamics Simulations.

We used coarse-grained molecular dynamics (CG MD) simulations to study protein-cholesterol interactions for different activation states of the A2A adenosine receptor (A2AR) and the A1 adenosine receptor (A1R) and predict new cholesterol binding sites indicating amino acid residues with a high residence time in three biologically relevant membranes. Compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)-cholesterol and POPC-phosphatidylinositol-bisphosphate (PIP2)-cholesterol, the plasma mimetic membrane best described the cholesterol binding sites previously detected for the inactive state of A2AR and revealed the binding sites with long-lasting amino acid residues. We observed that using the plasma mimetic membrane and plotting residues with cholesterol residence time ≥2 µs, our CG MD simulations captured most obviously the cholesterol-protein interactions. For the inactive A2AR, we identified one more binding site in which cholesterol is bound to residues with a long residence time compared to the previously detected, for the active A1R, three binding sites, and for the inactive A1R, two binding sites. We calculated that for the active states, cholesterol binds to residues with a much longer residence time compared to the inactive state for both A2AR and A1R. The stability of the identified binding sites to A1R or A2AR with CG MD simulations was additionally investigated with potential of mean force calculations using umbrella sampling. We observed that the binding sites with residues to which cholesterol has a long residence time in A2AR have shallow binding free energy minima compared to the related binding sites in A1R, suggesting a stronger binding for cholesterol to A1R. The differences in binding sites in which cholesterol is stabilized and interacts with residues with a long residence time between active and inactive states of A1R and A2AR can be important for differences in functional activity and orthosteric agonist or antagonist affinity and can be used for the design of allosteric modulators, which can bind through lipid pathways. We observed a stronger binding for cholesterol to A1R (i.e., generally higher association rates) compared to A2AR, which remains to be demonstrated. For the active states, cholesterol binds to residues with much longer residence times compared to the inactive state for both A2AR and A1R. Taken together, binding sites of active A1R may be considered as promising allosteric targets.

Synthesis and Testing of Analogs of the Tuberculosis Drug Candidate SQ109 against Bacteria and Protozoa: Identification of Lead Compounds against Mycobacterium abscessus and Malaria Parasites.

SQ109 is a tuberculosis drug candidate that has high potency against Mycobacterium tuberculosis and is thought to function at least in part by blocking cell wall biosynthesis by inhibiting the MmpL3 transporter. It also has activity against bacteria and protozoan parasites that lack MmpL3, where it can act as an uncoupler, targeting lipid membranes and Ca2+ homeostasis. Here, we synthesized 18 analogs of SQ109 and tested them against M. smegmatis, M. tuberculosis, M. abscessus, Bacillus subtilis, and Escherichia coli, as well as against the protozoan parasites Trypanosoma brucei, T. cruzi, Leishmania donovani, L. mexicana, and Plasmodium falciparum. Activity against the mycobacteria was generally less than with SQ109 and was reduced by increasing the size of the alkyl adduct, but two analogs were ∼4-8-fold more active than SQ109 against M. abscessus, including a highly drug-resistant strain harboring an A309P mutation in MmpL3. There was also better activity than found with SQ109 with other bacteria and protozoa. Of particular interest, we found that the adamantyl C-2 ethyl, butyl, phenyl, and benzyl analogs had 4-10× increased activity against P. falciparum asexual blood stages, together with low toxicity to a human HepG2 cell line, making them of interest as new antimalarial drug leads. We also used surface plasmon resonance to investigate the binding of inhibitors to MmpL3 and differential scanning calorimetry to investigate binding to lipid membranes. There was no correlation between MmpL3 binding and M. tuberculosis or M. smegmatis cell activity, suggesting that MmpL3 is not a major target in mycobacteria. However, some of the more active species decreased lipid phase transition temperatures, indicating increased accumulation in membranes, which is expected to lead to enhanced uncoupler activity.

A proof-of-concept study of the secondary structure of influenza A, B M2 and MERS- and SARS-CoV E transmembrane peptides using folding molecular dynamics simulations in a membrane mimetic solvent.

Here, we have carried out a proof-of-concept molecular dynamics (MD) simulation with adaptive tempering in a membrane mimetic environment to study the folding of single-pass membrane peptides. We tested the influenza A M2 viroporin, influenza B M2 viroporin, and protein E from coronaviruses MERS-Cov-2 and SARS-CoV-2 peptides with known experimental secondary structures in membrane bilayers. The two influenza-derived peptides are significantly different in the peptide sequence and secondary structure and more polar than the two coronavirus-derived peptides. Through a total of more than 50 µs of simulation time that could be accomplished in trifluoroethanol (TFE), as a membrane model, we characterized comparatively the folding behavior, helical stability, and helical propensity of these transmembrane peptides that match perfectly their experimental secondary structures, and we identified common motifs that reflect their quaternary organization and known (or not) biochemical function. We showed that BM2 is organized into two structurally distinct parts: a significantly more stable N-terminal half, and a fast-converting C-terminal half that continuously folds and unfolds between α-helical structures and non-canonical structures, which are mostly turns. In AM2, both the N-terminal half and C-terminal half are very flexible. In contrast, the two coronavirus-derived transmembrane peptides are much more stable and fast helix-formers when compared with the influenza ones. In particular, the SARS-derived peptide E appears to be the fastest and most stable helix-former of all the four viral peptides studied, with a helical structure that persists almost without disruption for the whole of its 10 µs simulation. By comparing the results with experimental observations, we benchmarked TFE in studying the conformation of membrane and hydrophobic peptides. This work provided accurate results suggesting a methodology to run long MD simulations and predict structural properties of biologically important membrane peptides.

Dual A1/A3 Adenosine Receptor Antagonists: Binding Kinetics and Structure-Activity Relationship Studies Using Mutagenesis and Alchemical Binding Free Energy Calculations.

Drugs targeting adenosine receptors (AR) can provide treatment for diseases. We report the identification of 7-(phenylamino)-pyrazolo[3,4-c]pyridines L2-L10, A15, and A17 as low-micromolar to low-nanomolar A1R/A3R dual antagonists, with 3-phenyl-5-cyano-7-(trimethoxyphenylamino)-pyrazolo[3,4-c]pyridine (A17) displaying the highest affinity at both receptors with a long residence time of binding, as determined using a NanoBRET-based assay. Two binding orientations of A17 produce stable complexes inside the orthosteric binding area of A1R in molecular dynamics (MD) simulations, and we selected the most plausible orientation based on the agreement with alanine mutagenesis supported by affinity experiments. Interestingly, for drug design purposes, the mutation of L2506.51 to alanine increased the binding affinity of A17 at A1R. We explored the structure-activity relationships against A1R using alchemical binding free energy calculations with the thermodynamic integration coupled with the MD simulation (TI/MD) method, applied on the whole G-protein-coupled receptor-membrane system, which showed a good agreement (r = 0.73) between calculated and experimental relative binding free energies.

Synthesis and in vitro proof-of-concept studies on bispecific iron oxide magnetic nanoparticles targeting PSMA and GRP receptors for PET/MR imaging of prostate cancer.

Prostate cancer (PCa) is the most common malignancy worldwide in men. This is a proof-of-concept study describing the development of 68Ga-magnetic iron oxide nanoparticles (mNP) targeting prostate specific membrane antigen (PSMA) and gastrin releasing peptide (GRPR) receptors as potential tools for diagnosis of PCa with PET/MRI. Two pharmacophores targeting PSMA, 1, and GRPR, 2, were coupled to mNPs carrying -SH (mNP-S1/2) or -NH2 (mNP-N1/2) groups. The mNP-S1/2 and mNP-N1/2 were characterized for their size, zeta potential, structure, and efficiency of functionalization using dynamic light scattering (DLS), FT-IR and RP-HPLC. A direct 68Ga-labelling procedure was followed, where 68Ga-mNP-N1/2 proved superior to 68Ga-mNP-S1/2 regarding radiolabelling efficiency, and thus were further evaluated in vitro. Toxicity studies in PCa cells (LNCaP, PC-3) showed low toxicity, and minimal hemolysis of red blood cells. In vitro assays in cells expressing PSMA (LNCaP), and GRPR (PC-3), showed specific time-dependent binding (40 min to plateau), high avidity (PC-3: Kd = 28.27 nM, LNCaP: Kd = 11.49 nM) and high internalization rates for 68Ga-mNP-N1/2 in both cell lines.

Discovery of a High Affinity Adenosine A1/A3 Receptor Antagonist with a Novel 7-Amino-pyrazolo[3,4-d]pyridazine Scaffold.

Here we describe the design and synthesis of pyrazolo[3,4-d]pyridazines as adenosine receptor (AR) ligands. We demonstrate that the introduction of a 3-phenyl group, together with a 7-benzylamino and 1-methyl group at the pyrazolopyridazine scaffold, generated the antagonist compound 10b, which displayed 21 nM affinity and a residence time of ∼60 min, for the human A1R, 55 nM affinity and a residence time of ∼73 min, for the human A3R and 1.7 µΜ affinity for the human A2BR while not being toxic. Strikingly, the 2-methyl analog of 10b, 15b, had no significant affinity. Docking calculations and molecular dynamics simulations of the ligands inside the orthosteric binding area suggested that the 2-methyl group in 15b hinders the formation of hydrogen bonding interactions with N6.55 which are considered critical for the stabilization inside the orthosteric binding cavity. We, therefore, demonstrate that 10a is a novel scaffold for the development of high affinity AR ligands. From the mutagenesis experiments the biggest effect was observed for the Y2717.46A mutation which caused an ∼10-fold reduction in the binding affinity of 10b.

Investigation of Tumor Cells and Receptor-Ligand Simulation Models for the Development of PET Imaging Probes Targeting PSMA and GRPR and a Possible Crosstalk between the Two Receptors.

Prostate-specific membrane antigen (PSMA) and gastrin-releasing peptide receptor (GRPR) have both been used in nuclear medicine as targets for molecular imaging and therapy of prostate (PCa) and breast cancer (BCa). Three bioconjugate probes, the PSMA specific: [68Ga]Ga-1, ((HBED-CC)-Ahx-Lys-NH-CO-NH Glu or PSMA-11), the GRPR specific: [68Ga]Ga-2, ((HBED-CC)-4-amino-1-carboxymethyl piperidine-[D-Phe6, Sta13]BN(6-14), a bombesin (BN) analogue), and 3 (the BN analogue: 4-amino-1-carboxymethyl piperidine-[(R)-Phe6, Sta13]BN(6-14) connected with the fluorescent dye, BDP-FL), were synthesized and tested in vitro with PCa and BCa cell lines, more specifically, with PCa cells, PC-3 and LNCaP, with BCa cells, T47D, MDA-MB-231, and with the in-house created PSMA-overexpressing PC-3(PSMA), T47D(PSMA), and MDA-MB-231(PSMA). In addition, biomolecular simulations were conducted on the association of 1 and 2 with PSMA and GRPR. The PSMA overexpression resulted in an increase of cell-bound radioligand [68Ga]Ga-1 (PSMA) for PCa and BCa cells and also of [68Ga]Ga-2 (GRPR), especially in those cell lines already expressing GRPR. The results were confirmed by fluorescence-activated cell sorting with a PE-labeled PSMA-specific antibody and the fluorescence tracer 3. The docking calculations and molecular dynamics simulations showed how 1 enters the PSMA funnel region and how pharmacophore Glu-urea-Lys interacts with the arginine patch, the S1', and S1 subpockets by forming hydrogen and van der Waals bonds. The chelating moiety of 1, that is, HBED-CC, forms additional stabilizing hydrogen bonding and van der Waals interactions in the arene-binding site. Ligand 2 is diving into the GRPR transmembrane (TM) helical cavity, thereby forming hydrogen bonds through its amidated end, water-mediated hydrogen bonds, and π-π interactions. Our results provide valuable information regarding the molecular mechanisms involved in the interactions of 1 and 2 with PSMA and GRPR, which might be useful for the diagnostic imaging and therapy of PCa and BCa.

Amantadine has potential for the treatment of COVID-19 because it inhibits known and novel ion channels encoded by SARS-CoV-2.

The dire need for COVID-19 treatments has inspired strategies of repurposing approved drugs. Amantadine has been suggested as a candidate, and cellular as well as clinical studies have indicated beneficial effects of this drug. We demonstrate that amantadine and hexamethylene-amiloride (HMA), but not rimantadine, block the ion channel activity of Protein E from SARS-CoV-2, a conserved viroporin among coronaviruses. These findings agree with their binding to Protein E as evaluated by solution NMR and molecular dynamics simulations. Moreover, we identify two novel viroporins of SARS-CoV-2; ORF7b and ORF10, by showing ion channel activity in a X. laevis oocyte expression system. Notably, amantadine also blocks the ion channel activity of ORF10, thereby providing two ion channel targets in SARS-CoV-2 for amantadine treatment in COVID-19 patients. A screen of known viroporin inhibitors on Protein E, ORF7b, ORF10 and Protein 3a from SARS-CoV-2 revealed inhibition of Protein E and ORF7b by emodin and xanthene, the latter also blocking Protein 3a. This illustrates a general potential of well-known ion channel blockers against SARS-CoV-2 and specifically a dual molecular basis for the promising effects of amantadine in COVID-19 treatment. We therefore propose amantadine as a novel, cheap, readily available and effective way to treat COVID-19.

Inside and Out of the Pore: Comparing Interactions and Molecular Dynamics of Influenza A M2 Viroporin Complexes in Standard Lipid Bilayers.

Ion channels located at viral envelopes (viroporins) have a critical function for the replication of infectious viruses and are important drug targets. Over the last decade, the number and duration of molecular dynamics (MD) simulations of the influenza A M2 ion channel owing to the increased computational efficiency. Here, we aimed to define the system setup and simulation conditions for the correct description of the protein-pore and the protein-lipid interactions for influenza A M2 in comparison with experimental data. We performed numerous MD simulations of the influenza A M2 protein in complex with adamantane blockers in standard lipid bilayers using OPLS2005 and CHARMM36 (C36) force fields. We explored the effect of varying the M2 construct (M2(22-46) and M2(22-62)), the lipid buffer size and type (stiffer DMPC or softer POPC with or without 20% cholesterol), the simulation time, the H37 protonation site (Nδ or Νε), the conformational state of the W41 channel gate, and M2's cholesterol binding sites (BSs). We report that the 200 ns MD with M2(22-62) (having Nε Η37) in the 20 Å lipid buffer with the C36 force field accurately describe: (a) the M2 pore structure and interactions inside the pore, that is, adamantane channel blocker location, water clathrate structure, and water or chloride anion blockage/passage from the M2 pore in the presence of a channel blocker and (b) interactions between M2 and the membrane environment as reflected by the calculation of the M2 bundle tilt, folding of amphipathic helices, and cholesterol BSs. Strikingly, we also observed that the C36 1 µs MD simulations using M2(22-62) embedded in a 20 Å POPC:cholesterol (5:1) scrambled membrane produced frequent interactions with cholesterol, which when combined with computational kinetic analysis, revealed the experimentally observed BSs of cholesterol and suggested three similarly long-interacting positions in the top leaflet that have previously not been observed experimentally. These findings promise to be useful for other viroporin systems.

Discovery of SARS-CoV-2 Papain-like Protease Inhibitors through a Combination of High-Throughput Screening and a FlipGFP-Based Reporter Assay.

The papain-like protease (PLpro) of SARS-CoV-2 is a validated antiviral drug target. Through a fluorescence resonance energy transfer-based high-throughput screening and subsequent lead optimization, we identified several PLpro inhibitors including Jun9-72-2 and Jun9-75-4 with improved enzymatic inhibition and antiviral activity compared to GRL0617, which was reported as a SARS-CoV PLpro inhibitor. Significantly, we developed a cell-based FlipGFP assay that can be applied to predict the cellular antiviral activity of PLpro inhibitors in the BSL-2 setting. X-ray crystal structure of PLpro in complex with GRL0617 showed that binding of GRL0617 to SARS-CoV-2 induced a conformational change in the BL2 loop to a more closed conformation. Molecular dynamics simulations showed that Jun9-72-2 and Jun9-75-4 engaged in more extensive interactions than GRL0617. Overall, the PLpro inhibitors identified in this study represent promising candidates for further development as SARS-CoV-2 antivirals, and the FlipGFP-PLpro assay is a suitable surrogate for screening PLpro inhibitors in the BSL-2 setting.