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).

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.

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.

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.

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)P₂ or PIP₂), a minor phospholipid component and of cholesterol on the function-related conformational equilibria of the human A_2A adenosine receptor (A_2AR), a representative receptor in class A GPCR. Several GPCRs of class A interacted with PIP₂ 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 A_2A adenosine receptor (A_2AR) and the A₁ adenosine receptor (A₁R) 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 (PIP₂)-cholesterol, the plasma mimetic membrane best described the cholesterol binding sites previously detected for the inactive state of A_2AR 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 A_2AR, 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 A₁R, three binding sites, and for the inactive A₁R, 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 A_2AR and A₁R. The stability of the identified binding sites to A₁R or A_2AR 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 A_2AR have shallow binding free energy minima compared to the related binding sites in A₁R, suggesting a stronger binding for cholesterol to A₁R. 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 A₁R and A_2AR 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 A₁R (i.e., generally higher association rates) compared to A_2AR, 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 A_2AR and A₁R. Taken together, binding sites of active A₁R may be considered as promising allosteric targets.

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 A₁R/A₃R 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 A₁R 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 L250^6.51 to alanine increased the binding affinity of A17 at A₁R. We explored the structure-activity relationships against A₁R 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.

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 A₁R, 55 nM affinity and a residence time of ∼73 min, for the human A₃R and 1.7 μΜ affinity for the human A_2BR 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 N^6.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 Y271^7.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: [⁶⁸Ga]Ga-1, ((HBED-CC)-Ahx-Lys-NH-CO-NH Glu or PSMA-11), the GRPR specific: [⁶⁸Ga]Ga-2, ((HBED-CC)-4-amino-1-carboxymethyl piperidine-[D-Phe⁶, Sta¹³]BN(6-14), a bombesin (BN) analogue), and 3 (the BN analogue: 4-amino-1-carboxymethyl piperidine-[(R)-Phe⁶, Sta¹³]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 [⁶⁸Ga]Ga-1 (PSMA) for PCa and BCa cells and also of [⁶⁸Ga]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.

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.

Improved Synthesis of the Antitubercular Agent SQ109

We present here an improved procedure for the preparation of the promising antitubercular drug SQ109 that is currently in phase Ib/III of clinical trials against Mycobacterium tuberculosis. We investigated and tested the literature synthetic procedure that enables the development of structure–activity relationships and report the observed inconsistencies as well as presenting improvements or novelties for the more efficient preparation of SQ109. Most significantly we applied a novel reduction step of the aminoamide precursor using Me3SiCl­/LiAlH4 under mild conditions. These findings are important for research groups investigating the efficacy of this drug and analogues in academia and industry.

Rimantadine Binds to and Inhibits the Influenza A M2 Proton Channel without Enantiomeric Specificity

The influenza A M2 wild-type (WT) proton channel is the target of the anti-influenza drug rimantadine. Rimantadine has two enantiomers, though most investigations into drug binding and inhibition have used a racemic mixture. Solid-state NMR experiments using the full length-M2 WT have shown significant spectral differences that were interpreted to indicate tighter binding for (R)- vs (S)-rimantadine. However, it was unclear if this correlates with a functional difference in drug binding and inhibition. Using X-ray crystallography, we have determined that both (R)- and (S)-rimantadine bind to the M2 WT pore with slight differences in the hydration of each enantiomer. However, this does not result in a difference in potency or binding kinetics, as shown by similar values for kon, koff, and Kd in electrophysiological assays and for EC50 values in cellular assays. We concluded that the slight differences in hydration for the (R)- and (S)-rimantadine enantiomers are not relevant to drug binding or channel inhibition. To further explore the effect of the hydration of the M2 pore on binding affinity, the water structure was evaluated by grand canonical ensemble molecular dynamics simulations as a function of the chemical potential of the water. Initially, the two layers of ordered water molecules between the bound drug and the channel's gating His37 residues mask the drug's chirality. As the chemical potential becomes more unfavorable, the drug translocates down to the lower water layer, and the interaction becomes more sensitive to chirality. These studies suggest the feasibility of displacing the upper water layer and specifically recognizing the lower water layers in novel drugs.

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.

PET Diagnostic Molecules Utilizing Multimeric Cyclic RGD Peptide Analogs for Imaging Integrin αvβ3 Receptors

Multimeric ligands consisting of multiple pharmacophores connected to a single backbone have been widely investigated for diagnostic and therapeutic applications. In this review, we summarize recent developments regarding multimeric radioligands targeting integrin α_vβ₃ receptors on cancer cells for molecular imaging and diagnostic applications using positron emission tomography (PET). Integrin α_vβ₃ receptors are glycoproteins expressed on the cell surface, which have a significant role in tumor angiogenesis. They act as receptors for several extracellular matrix proteins exposing the tripeptide sequence arginine-glycine-aspartic (RGD). Cyclic RDG peptidic ligands c(RGD) have been developed for integrin α_vβ₃ tumor-targeting positron emission tomography (PET) diagnosis. Several c(RGD) pharmacophores, connected with the linker and conjugated to a chelator or precursor for radiolabeling with different PET radionuclides (¹⁸F, ⁶⁴Cu, and ⁶⁸Ga), have resulted in multimeric ligands superior to c(RGD) monomers. The binding avidity, pharmacodynamic, and PET imaging properties of these multimeric c(RGD) radioligands, in relation to their structural characteristics are analyzed and discussed. Furthermore, specific examples from preclinical studies and clinical investigations are included.

Structure and drug binding of the SARS-CoV-2 envelope protein transmembrane domain in lipid bilayers

An essential protein of the SARS-CoV-2 virus, the envelope protein E, forms a homopentameric cation channel that is important for virus pathogenicity. Here we report a 2.1-Å structure and the drug-binding site of E's transmembrane domain (ETM), determined using solid-state NMR spectroscopy. In lipid bilayers that mimic the endoplasmic reticulum-Golgi intermediate compartment (ERGIC) membrane, ETM forms a five-helix bundle surrounding a narrow pore. The protein deviates from the ideal α-helical geometry due to three phenylalanine residues, which stack within each helix and between helices. Together with valine and leucine interdigitation, these cause a dehydrated pore compared with the viroporins of influenza viruses and HIV. Hexamethylene amiloride binds the polar amino-terminal lumen, whereas acidic pH affects the carboxy-terminal conformation. Thus, the N- and C-terminal halves of this bipartite channel may interact with other viral and host proteins semi-independently. The structure sets the stage for designing E inhibitors as antiviral drugs.

Structure and inhibition of the SARS-CoV-2 main protease reveal strategy for developing dual inhibitors against Mpro and cathepsin L

The main protease (M^pro) of SARS-CoV-2 is a key antiviral drug target. While most M^pro inhibitors have a γ-lactam glutamine surrogate at the P1 position, we recently found that several M^pro inhibitors have hydrophobic moieties at the P1 site, including calpain inhibitors II and XII, which are also active against human cathepsin L, a host protease that is important for viral entry. In this study, we solved x-ray crystal structures of M^pro in complex with calpain inhibitors II and XII and three analogs of GC-376 The structure of M^pro with calpain inhibitor II confirmed that the S1 pocket can accommodate a hydrophobic methionine side chain, challenging the idea that a hydrophilic residue is necessary at this position. The structure of calpain inhibitor XII revealed an unexpected, inverted binding pose. Together, the biochemical, computational, structural, and cellular data presented herein provide new directions for the development of dual inhibitors as SARS-CoV-2 antivirals.

Ebselen, Disulfiram, Carmofur, PX-12, Tideglusib, and Shikonin Are Nonspecific Promiscuous SARS-CoV-2 Main Protease Inhibitors

Among the drug targets being investigated for SARS-CoV-2, the viral main protease (M^pro) is one of the most extensively studied. M^pro is a cysteine protease that hydrolyzes the viral polyprotein at more than 11 sites. It is highly conserved and has a unique substrate preference for glutamine in the P1 position. Therefore, M^pro inhibitors are expected to have broad-spectrum antiviral activity and a high selectivity index. Structurally diverse compounds have been reported as M^pro inhibitors. In this study, we investigated the mechanism of action of six previously reported M^pro inhibitors, ebselen, disulfiram, tideglusib, carmofur, shikonin, and PX-12, using a consortium of techniques including FRET-based enzymatic assay, thermal shift assay, native mass spectrometry, cellular antiviral assays, and molecular dynamics simulations. Collectively, the results showed that the inhibition of M^pro by these six compounds is nonspecific and that the inhibition is abolished or greatly reduced with the addition of reducing reagent 1,4-dithiothreitol (DTT). Without DTT, these six compounds inhibit not only M^pro but also a panel of viral cysteine proteases including SARS-CoV-2 papain-like protease and 2A^pro and 3C^pro from enterovirus A71 (EV-A71) and EV-D68. However, none of the compounds inhibits the viral replication of EV-A71 or EV-D68, suggesting that the enzymatic inhibition potency IC₅₀ values obtained in the absence of DTT cannot be used to faithfully predict their cellular antiviral activity. Overall, we provide compelling evidence suggesting that these six compounds are nonspecific SARS-CoV-2 M^pro inhibitors and urge the scientific community to be stringent with hit validation.

The balance between side-chain and backbone-driven association in folding of the α-helical influenza A transmembrane peptide

The correct balance between attractive, repulsive and peptide hydrogen bonding interactions must be attained for proteins to fold correctly. To investigate these important contributors, we sought a comparison of the folding between two 25-residues peptides, the influenza A M2 protein transmembrane domain (M2TM) and the 25-Ala (Ala₂₅ ). M2TM forms a stable α-helix as is shown by circular dichroism (CD) experiments. Molecular dynamics (MD) simulations with adaptive tempering show that M2TM monomer is more dynamic in nature and quickly interconverts between an ensemble of various α-helical structures, and less frequently turns and coils, compared to one α-helix for Ala₂₅ . DFT calculations suggest that folding from the extended structure to the α-helical structure is favored for M2TM compared with Ala₂₅ . This is due to CH⋯O attractive interactions which favor folding to the M2TM α-helix, and cannot be described accurately with a force field. Using natural bond orbital (NBO) analysis and quantum theory atoms in molecules (QTAIM) calculations, 26 CH⋯O interactions and 22 NH⋯O hydrogen bonds are calculated for M2TM. The calculations show that CH⋯O hydrogen bonds, although individually weaker, have a cumulative effect that cannot be ignored and may contribute as much as half of the total hydrogen bonding energy, when compared to NH⋯O, to the stabilization of the α-helix in M2TM. Further, a strengthening of NH⋯O hydrogen bonding interactions is calculated for M2TM compared to Ala₂₅ . Additionally, these weak CH⋯O interactions can dissociate and associate easily leading to the ensemble of folded structures for M2TM observed in folding MD simulations.

Structure and Drug Binding of the SARS-CoV-2 Envelope Protein in Phospholipid Bilayers

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the ongoing COVID-19 pandemic. Successful development of vaccines and antivirals against SARS-CoV-2 requires a comprehensive understanding of the essential proteins of the virus. The envelope (E) protein of SARS-CoV-2 assembles into a cation-selective channel that mediates virus budding, release, and host inflammation response. E blockage reduces virus pathogenicity while E deletion attenuates the virus. Here we report the 2.4 Å structure and drug-binding site of E's transmembrane (TM) domain, determined using solid-state nuclear magnetic resonance (NMR) spectroscopy. In lipid bilayers that mimic the endoplasmic reticulum Golgi intermediate compartment (ERGIC) membrane, ETM forms a five-helix bundle surrounding a narrow central pore. The middle of the TM segment is distorted from the ideal a-helical geometry due to three regularly spaced phenylalanine residues, which stack within each helix and between neighboring helices. These aromatic interactions, together with interhelical Val and Leu interdigitation, cause a dehydrated pore compared to the viroporins of influenza and HIV viruses. Hexamethylene amiloride and amantadine bind shallowly to polar residues at the N-terminal lumen, while acidic pH affects the C-terminal conformation. These results indicate that SARS-CoV-2 E forms a structurally robust but bipartite channel whose N- and C-terminal halves can interact with drugs, ions and other viral and host proteins semi-independently. This structure establishes the atomic basis for designing E inhibitors as antiviral drugs against SARS-CoV-2.

Chemical Probes for Blocking of Influenza A M2 Wild-type and S31N Channels

We report on using the synthetic aminoadamantane-CH₂-aryl derivatives 1-6 as sensitive probes for blocking M2 S31N and influenza A virus (IAV) M2 wild-type (WT) channels as well as virus replication in cell culture. The binding kinetics measured using electrophysiology (EP) for M2 S31N channel are very dependent on the length between the adamantane moiety and the first ring of the aryl headgroup realized in 2 and 3 and the girth and length of the adamantane adduct realized in 4 and 5. Study of 1-6 shows that, according to molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculations, all bind in the M2 S31N channel with the adamantyl group positioned between V27 and G34 and the aryl group projecting out of the channel with the phenyl (or isoxazole in 6) embedded in the V27 cluster. In this outward binding configuration, an elongation of the ligand by only one methylene in rimantadine 2 or using diamantane or triamantane instead of adamantane in 4 and 5, respectively, causes incomplete entry and facilitates exit, abolishing effective block compared to the amantadine derivatives 1 and 6. In the active M2 S31N blockers 1 and 6, the phenyl and isoxazolyl head groups achieve a deeper binding position and high k_on/low k_off and high k_on/high k_off rate constants, compared to inactive 2-5, which have much lower k_on and higher k_off. Compounds 1-5 block the M2 WT channel by binding in the longer area from V27-H37, in the inward orientation, with high k_on and low k_off rate constants. Infection of cell cultures by influenza virus containing M2 WT or M2 S31N is inhibited by 1-5 or 1-4 and 6, respectively. While 1 and 6 block infection through the M2 block mechanism in the S31N variant, 2-4 may block M2 S31N virus replication in cell culture through the lysosomotropic effect, just as chloroquine is thought to inhibit SARS-CoV-2 infection.

Ebselen, disulfiram, carmofur, PX-12, tideglusib, and shikonin are non-specific promiscuous SARS-CoV-2 main protease inhibitors

There is an urgent need for vaccines and antiviral drugs to combat the COVID-19 pandemic. Encouraging progress has been made in developing antivirals targeting SARS-CoV-2, the etiological agent of COVID-19. Among the drug targets being investigated, the viral main protease (M pro ) is one of the most extensively studied drug targets. M pro is a cysteine protease that hydrolyzes the viral polyprotein at more than 11 sites and it is highly conserved among coronaviruses. In addition, M pro has a unique substrate preference for glutamine in the P1 position. Taken together, it appears that M pro inhibitors can achieve both broad-spectrum antiviral activity and a high selectivity index. Structurally diverse compounds have been reported as M pro inhibitors, with several of which also showed antiviral activity in cell culture. In this study, we investigated the mechanism of action of six previously reported M pro inhibitors, ebselen, disulfiram, tideglusib, carmofur, shikonin, and PX-12 using a consortium of techniques including FRET-based enzymatic assay, thermal shift assay, native mass spectrometry, cellular antiviral assays, and molecular dynamics simulations. Collectively, the results showed that the inhibition of M pro by these six compounds is non-specific and the inhibition is abolished or greatly reduced with the addition of reducing reagent DTT. In the absence of DTT, these six compounds not only inhibit M pro , but also a panel of viral cysteine proteases including SARS-CoV-2 papain-like protease, the 2A pro and 3C pro from enterovirus A71 (EV-A71) and EV-D68. However, none of the compounds inhibits the viral replication of EV-A71 or EV-D68, suggesting that the enzymatic inhibition potency IC 50 values obtained in the absence of DTT cannot be used to faithfully predict their cellular antiviral activity. Overall, we provide compelling evidence suggesting that ebselen, disulfiram, tideglusib, carmofur, shikonin, and PX-12 are non-specific SARS-CoV-2 M pro inhibitors, and urge the scientific community to be stringent with hit validation.

Structure and inhibition of the SARS-CoV-2 main protease reveals strategy for developing dual inhibitors against Mpro and cathepsin L

The main protease (Mpro) of SARS-CoV-2, the pathogen responsible for the COVID-19 pandemic, is a key antiviral drug target. While most SARS-CoV-2 Mpro inhibitors have a γ-lactam glutamine surrogate at the P1 position, we recently discovered several Mpro inhibitors have hydrophobic moieties at the P1 site, including calpain inhibitors II/XII, which are also active against human cathepsin L, a host-protease that is important for viral entry. To determine the binding mode of these calpain inhibitors and establish a structure-activity relationship, we solved X-ray crystal structures of Mpro in complex with calpain inhibitors II and XII, and three analogues of GC-376, one of the most potent Mpro inhibitors in vitro. The structure of Mpro with calpain inhibitor II confirmed the S1 pocket of Mpro can accommodate a hydrophobic methionine side chain, challenging the idea that a hydrophilic residue is necessary at this position. Interestingly, the structure of calpain inhibitor XII revealed an unexpected, inverted binding pose where the P1' pyridine inserts in the S1 pocket and the P1 norvaline is positioned in the S1' pocket. The overall conformation is semi-helical, wrapping around the catalytic core, in contrast to the extended conformation of other peptidomimetic inhibitors. Additionally, the structures of three GC-376 analogues UAWJ246, UAWJ247, and UAWJ248 provide insight to the sidechain preference of the S1', S2, S3 and S4 pockets, and the superior cell-based activity of the aldehyde warhead compared with the α-ketoamide. Taken together, the biochemical, computational, structural, and cellular data presented herein provide new directions for the development of Mpro inhibitors as SARS-CoV-2 antivirals.

Investigation of the Drug Resistance Mechanism of M2-S31N Channel Blockers through Biomolecular Simulations and Viral Passage Experiments

Recent efforts in drug development against influenza A virus (IAV) M2 proton channel S31N mutant resulted in conjugates of amantadine linked with aryl head heterocycles. To understand the mechanism of drug resistance, we chose a representative M2-S31N inhibitor, compound 3, as a chemical probe to identify resistant mutants. To increase the possibility of identifying novel resistant mutants, serial viral passage experiments were performed with multiple strains of H1N1 and H3N2 viruses in different cell lines. This approach not only identified M2 mutations around the drug-binding site, including the pore-lining residues (V27A, V27F, N31S, and G34E) and an interhelical residue (I32N), but also a new allosteric mutation (R45H), in addition to L46P previously identified, located at the C-terminus of M2 that is more than 10 Å away from the drug-binding site. The effects of each mutation were next investigated using electrophysiology, recombinant viruses, and molecular dynamics (MD) simulations. The reduced sensitivity in channel blockage correlated with increased drug resistance in antiviral assays using recombinant viruses. The MD simulations show that the V27A, V27F, G34E, and R45H mutations increase the diameter and hydration state of the pore in complex with compound 3. The Molecular Mechanics Generalized Born (MM-GBSA) calculations result in more positive binding free energies for the complexes of resistant M2 (V27A, V27F, G34E, R45H) with compound 3 compared to the stable complexes (S31N and I32N). Overall, this is the first systematic study of the drug resistance mechanism of M2-S31N channel blockers using multiple viruses in different cell lines.

X-ray Crystal Structures of the Influenza M2 Proton Channel Drug-Resistant V27A Mutant Bound to a Spiro-Adamantyl Amine Inhibitor Reveal the Mechanism of Adamantane Resistance

The V27A mutation confers adamantane resistance on the influenza A matrix 2 (M2) proton channel and is becoming more prevalent in circulating populations of influenza A virus. We have used X-ray crystallography to determine structures of a spiro-adamantyl amine inhibitor bound to M2(22-46) V27A and also to M2(21-61) V27A in the Inwardclosed conformation. The spiro-adamantyl amine binding site is nearly identical for the two crystal structures. Compared to the M2 "wild type" (WT) with valine at position 27, we observe that the channel pore is wider at its N-terminus as a result of the V27A mutation and that this removes V27 side chain hydrophobic interactions that are important for binding of amantadine and rimantadine. The spiro-adamantyl amine inhibitor blocks proton conductance in the WT and V27A mutant channels by shifting its binding site in the pore depending on which residue is present at position 27. Additionally, in the structure of the M2(21-61) V27A construct, the C-terminus of the channel is tightly packed relative to that of the M2(22-46) construct. We observe that residues Asp44, Arg45, and Phe48 face the center of the channel pore and would be well-positioned to interact with protons exiting the M2 channel after passing through the His37 gate. A 300 ns molecular dynamics simulation of the M2(22-46) V27A-spiro-adamantyl amine complex predicts with accuracy the position of the ligands and waters inside the pore in the X-ray crystal structure of the M2(22-46) V27A complex.

The key role of protodeauration in the gold-catalyzed reaction of 1,3-diynes with pyrrole and indole to form complex heterocycles

The mechanism of indole and carbazole formation via a formal [4 + 2] cycloaddition strategy is dominated by the protodeauration step.

The boundary lipid around DMPC-spanning influenza A M2 transmembrane domain channels: Its structure and potential for drug accommodation

We have investigated the perturbation of influenza A M2TM in DMPC bilayers. We have shown that (a) DSC and SAXS detect changes in membrane organization caused by small changes (micromolar) in M2TM or aminoadamantane concentration and aminoadamantane structure, by comparison of amantadine and spiro[pyrrolidine-2,2'-adamantane] (AK13), (b) that WAXS and MD can suggest details of ligand topology. DSC and SAXS show that at a low M2TM micromolar concentration in DPMC bilayers, two lipid domains are observed, which likely correspond to M2TM boundary lipids and bulk-like lipids. At higher M2TM concentrations, one domain only is identified, which constitutes essentially all of the lipid molecules behaving as boundary lipids. According to SAXS, WAXS, and DSC in the absence of M2TM, both aminoadamantane drugs exert a similar perturbing effect on the bilayer at low concentrations. At the same concentrations of the drug when M2TM is present, amantadine and, to a lesser extent, AK13 cause, according to WAXS, a significant disordering of chain-stacking, which also leads to the formation of two lipid domains. This effect is likely due, according to MD simulations, to the preference of the more lipophilic AK13 to locate closer to the lateral surfaces of M2TM when compared to amantadine, which forms stronger ionic interactions with phosphate groups. The preference of AK13 to concentrate inside the lipid bilayer close to the exterior of the hydrophobic M2TM helices may contribute to its higher binding affinity compared to amantadine.

Insights to the Binding of a Selective Adenosine A3 Receptor Antagonist Using Molecular Dynamic Simulations, MM-PBSA and MM-GBSA Free Energy Calculations, and Mutagenesis

Adenosine A₃ receptor (A₃R) is a promising drug target cancer and for a number of other conditions like inflammatory diseases, including asthma and rheumatoid arthritis, glaucoma, chronic obstructive pulmonary disease, and ischemic injury. Currently, there is no experimentally determined structure of A₃R. We explored the binding profile of O4-{[3-(2,6-dichlorophenyl)-5-methylisoxazol-4-yl]carbonyl}-2-methyl-1,3-thiazole-4-carbohydroximamide (K18), which is a new specific and competitive antagonist at the orthosteric binding site of A₃R. MD simulations and MM-GBSA calculations of the WT A₃R in complex with K18 combined with in vitro mutagenic studies show that the most plausible binding conformation for the dichlorophenyl group of K18 is oriented toward trans-membrane helices (TM) 5, 6 and reveal important residues for binding. Further, MM-GBSA calculations distinguish mutations that reduce or maintain or increase antagonistic activity. Our studies show that selectivity of K18 toward A₃R is defined not only by direct interactions with residues within the orthosteric binding area but also by remote residues playing a significant role. Although V169^5.30 is considered to be a selectivity filter for A₃R binders, when it was mutated to glutamic acid, K18 maintained antagonistic potency, in agreement with our previous results obtained for agonists binding profile investigation. Mutation of the direct interacting residue L90^3.32 in the low region and the remote L264^7.35 in the middle/upper region to alanine increases antagonistic potency, suggesting an empty space in the orthosteric area available for increasing antagonist potency. These results approve the computational model for the description of K18 binding at A₃R, which we previously performed for agonists binding to A₃R, and the design of more effective antagonists based on K18.

Structural Characterization of Agonist Binding to an A3 Adenosine Receptor through Biomolecular Simulations and Mutagenesis Experiments

The adenosine A₃ receptor (A₃R) binds adenosine and is a drug target against cancer cell proliferation. Currently, there is no experimental structure of A₃R. Here, we have generated a molecular model of A₃R in complex with two agonists, the nonselective 1-(6-amino-9H-purin-9-yl)-1-deoxy-N-ethyl-β-d-ribofuranuronamide (NECA) and the selective 1-deoxy-1-[6-[[(3-iodophenyl)methyl]amino]-9H-purin-9-yl]-N-methyl-β-d-ribofuranuronamide (IB-MECA). Molecular dynamics simulations of the wild-type A₃R in complex with both agonists, combined with in vitro mutagenic studies revealed important residues for binding. Further, molecular mechanics-generalized Born surface area calculations were able to distinguish mutations that reduce or negate agonistic activity from those that maintained or increased the activity. Our studies reveal that selectivity of IB-MECA toward A₃R requires not only direct interactions with residues within the orthosteric binding area but also with remote residues. Although V169^5.30 is considered to be a selectivity filter for A₃R binders, when it was mutated to glutamic acid or alanine, the activity of IB-MECA increased by making new van der Waals contacts with TM5. This result may have implications in the design of new A₃R agonists.

The L46P mutant confers a novel allosteric mechanism of resistance toward the influenza A virus M2 S31N proton channel blockers

The Food and Drug Administration-approved influenza A antiviral amantadine inhibits the wild-type (WT) AM2 channel but not the S31N mutant predominantly found in circulating strains. In this study, serial viral passages were applied to select resistance against a newly developed isoxazole-conjugated adamantane inhibitor that targets the AM2 S31N channel. This led to the identification of the novel drug-resistant mutation L46P located outside the drug-binding site, which suggests an allosteric resistance mechanism. Intriguingly, when the L46P mutant was introduced to AM2 WT, the channel remained sensitive toward amantadine inhibition. To elucidate the molecular mechanism, molecular dynamics simulations and binding free energy molecular mechanics-generalized born surface area (MM-GBSA) calculations were performed on WT and mutant channels. It was found that the L46P mutation caused a conformational change in the N terminus of transmembrane residues 22-31 that ultimately broadened the drug-binding site of AM2 S31N inhibitor 4, which spans residues 26-34, but not of AM2 WT inhibitor amantadine, which spans residues 31-34. The MM-GBSA calculations showed stronger binding stability for 4 in complex with AM2 S31N compared with 4 in complex with AM2 S31N/L46P, and equal binding free energies of amantadine in complex with AM2 WT and AM2 L46P. Overall, these results demonstrate a unique allosteric resistance mechanism toward AM2 S31N channel blockers, and the L46P mutant represents the first experimentally confirmed drug-resistant AM2 mutant that is located outside of the pore where drug binds. SIGNIFICANCE STATEMENT: AM2 S31N is a high-profile antiviral drug target, as more than 95% of currently circulating influenza A viruses carry this mutation. Understanding the mechanism of drug resistance is critical in designing the next generation of AM2 S31N channel blockers. Using a previously developed AM2 S31N channel blocker as a chemical probe, this study was the first to identify a novel resistant mutant, L46P. The L46P mutant is located outside of the drug-binding site. Molecular dynamics simulations showed that L46P causes a dilation of drug-binding site between residues 22 and 31, which affects the binding of AM2 S31N channel blockers, but not the AM2 WT inhibitor amantadine.

Pharmacophore Generation and 3D-QSAR Model Development Using PHASE

Nowadays, the prediction of biological activity of novel compounds is one of the major challenges in drug design. Toward this aim a useful procedure is the development and application of predictive computational models using three-dimensional quantitative structure-activity relationship (3D-QSAR) methods, which can decrease the cost and time of biological experiments. In this chapter, the use of application PHASE is analyzed, which is a recent but already widespread method for pharmacophore- or atom-based 3D-QSAR model building. The main steps of procedure provided by PHASE are described, and a general workflow and important practical notes are referred. An attempt in order to design new chemotypes with enhanced cytotoxicity against K562 cells is also provided as an example for the 3D pharmacophore model generation on 33 novel (E)-α-benzylthiochalcones.

Host-Guest Interactions between Candesartan and Its Prodrug Candesartan Cilexetil in Complex with 2-Hydroxypropyl-β-cyclodextrin: On the Biological Potency for Angiotensin II Antagonism

Renin-angiotensin aldosterone system inhibitors are for a long time extensively used for the treatment of cardiovascular and renal diseases. AT1 receptor blockers (ARBs or sartans) act as antihypertensive drugs by blocking the octapeptide hormone Angiotensin II to stimulate AT1 receptors. The antihypertensive drug candesartan (CAN) is the active metabolite of candesartan cilexetil (Atacand, CC). Complexes of candesartan and candesartan cilexetil with 2-hydroxylpropyl-β-cyclodextrin (2-HP-β-CD) were characterized using high-resolution electrospray ionization mass spectrometry and solid state ¹³C cross-polarization/magic angle spinning nuclear magnetic resonance (CP/MAS NMR) spectroscopy. The ¹³C CP/MAS results showed broad peaks especially in the aromatic region, thus confirming the strong interactions between cyclodextrin and drugs. This experimental evidence was in accordance with molecular dynamics simulations and quantum mechanical calculations. The synthesized and characterized complexes were evaluated biologically in vitro. It was shown that as a result of CAN's complexation, CAN exerts higher antagonistic activity than CC. Therefore, a formulation of CC with 2-HP-β-CD is not indicated, while the formulation with CAN is promising and needs further investigation. This intriguing result is justified by the binding free energy calculations, which predicted efficient CC binding to 2-HP-β-CD, and thus, the molecule's availability for release and action on the target is diminished. In contrast, CAN binding was not favored, and this may allow easy release for the drug to exert its bioactivity.

1,2-Αnnulated Adamantane Heterocyclic Derivatives as Anti-Influenza Α Virus Agents

In this report we review our results on the development of 1,2-annulated adamantane heterocyclic derivatives and we discuss the structure-activity relationships obtained from their biological evaluation against influenza A virus. We have designed and synthesized numerous potent 1,2-annulated adamantane analogues of amantadine and rimantadine against influenza A targeting M2 protein the last 20 years. For their synthesis we utilized the key intermediates 2-(2-oxoadamantan-1-yl)acetic acid and 3-(2-oxoadamantan-1-yl)propanoic acid, which were obtained by a simple, fast and efficient synthetic protocol. The latter involved the treatment of protoadamantanone with different electrophiles and a carbon-skeleton rearrangement. These ketoesters offered a new pathway to the synthesis of 1,2-disubstituted adamantanes, which constitute starting materials for many molecules with pharmacological potential, such as the 1,2-annulated adamantane heterocyclic derivatives. To obtain additional insight for their binding to M2 protein three structurally similar 1,2-annulated adamantane piperidines, differing in nitrogen position, were studied using molecular dynamics (MD) simulations in palmitoyl-oleoyl-phosphatidyl-choline (POPC) hydrated bilayers.

Comparative Perturbation Effects Exerted by the Influenza A M2 WT Protein Inhibitors Amantadine and the Spiro[pyrrolidine-2,2'-adamantane] Variant AK13 to Membrane Bilayers Studied Using Biophysical Experiments and Molecular Dynamics Simulations

Aminoadamantane drugs are lipophilic amines that block the membrane-embedded influenza A M2 WT (wild type) ion channel protein. The comparative effects of amantadine ( Amt) and its synthetic spiro[pyrrolidine-2,2'-adamantane] (AK13) analogue in dimyristoylphosphatidylcholine (DMPC) bilayers were studied using a combination of experimental biophysical methods, differential scanning calorimetry (DSC), X-ray diffraction, solid-state NMR (ssNMR) spectroscopy, and molecular dynamics (MD) simulations. All three experimental methods pointed out that the two analogues perturbed drastically the DMPC bilayers with AK13 to be more effective at high concentrations. AK13 was tolerated in lipid bilayers at very high concentrations, while Amt was crystallized. This is an important consideration in the formulations of drugs as it designates a limitation of Amt incorporation. MD simulations verify provided details about the strong interactions of the drugs in the interface region between phosphoglycerol backbone and lipophilic segments. The two drugs form hydrogen bonding with both water and sn-2 carbonyls in their amine form or water and phosphate oxygens in their ammonium form. Such localization of the drugs explains the DMPC bilayers reorientation and their strong perturbing effect evidenced by all biophysical methodologies applied.

Inhibitors of the M2 Proton Channel Engage and Disrupt Transmembrane Networks of Hydrogen-Bonded Waters

Water-mediated interactions play key roles in drug binding. In protein sites with sparse polar functionality, a small-molecule approach is often viewed as insufficient to achieve high affinity and specificity. Here we show that small molecules can enable potent inhibition by targeting key waters. The M2 proton channel of influenza A is the target of the antiviral drugs amantadine and rimantadine. Structural studies of drug binding to the channel using X-ray crystallography have been limited because of the challenging nature of the target, with the one previously solved crystal structure limited to 3.5 Å resolution. Here we describe crystal structures of amantadine bound to M2 in the Inward_closed conformation (2.00 Å), rimantadine bound to M2 in both the Inward_closed (2.00 Å) and Inward_open (2.25 Å) conformations, and a spiro-adamantyl amine inhibitor bound to M2 in the Inward_closed conformation (2.63 Å). These X-ray crystal structures of the M2 proton channel with bound inhibitors reveal that ammonium groups bind to water-lined sites that are hypothesized to stabilize transient hydronium ions formed in the proton-conduction mechanism. Furthermore, the ammonium and adamantyl groups of the adamantyl-amine class of drugs are free to rotate in the channel, minimizing the entropic cost of binding. These drug-bound complexes provide the first high-resolution structures of drugs that interact with and disrupt networks of hydrogen-bonded waters that are widely utilized throughout nature to facilitate proton diffusion within proteins.

Discovery of Novel Adenosine Receptor Antagonists through a Combined Structure- and Ligand-Based Approach Followed by Molecular Dynamics Investigation of Ligand Binding Mode

An intense effort is made by pharmaceutical and academic research laboratories to identify and develop selective antagonists for each adenosine receptor (AR) subtype as potential clinical candidates for "soft" treatment of various diseases. Crystal structures of subtypes A_2A and A₁ARs offer exciting opportunities for structure-based drug design. In the first part of the present work, Maybridge HitFinder library of 14400 compounds was utilized to apply a combination of structure-based against the crystal structure of A_2AAR and ligand-based methodologies. The docking poses were rescored by CHARMM energy minimization and calculation of the desolvation energy using Poisson-Boltzmann equation electrostatics. Out of the eight selected and tested compounds, five were found positive hits (63% success). Although the project was initially focused on targeting A_2AAR, the identified antagonists exhibited low micromolar or micromolar affinity against A_2A/A₃, ARs, or A₃AR, respectively. Based on these results, 19 compounds characterized by novel chemotypes were purchased and tested. Sixteen of them were identified as AR antagonists with affinity toward combinations of the AR family isoforms (A_2A/A₃, A₁/A₃, A₁/A_2A/A₃, and A₃). The second part of this work involves the performance of hundreds of molecular dynamics (MD) simulations of complexes between the ARs and a total of 27 ligands to resolve the binding interactions of the active compounds, which were not achieved by docking calculations alone. This computational work allowed the prediction of stable and unstable complexes which agree with the experimental results of potent and inactive compounds, respectively. Of particular interest is that the 2-amino-thiophene-3-carboxamides, 3-acylamino-5-aryl-thiophene-2-carboxamides, and carbonyloxycarboximidamide derivatives were found to be selective and possess a micromolar to low micromolar affinity for the A₃ receptor.

Unraveling the Binding, Proton Blockage, and Inhibition of Influenza M2 WT and S31N by Rimantadine Variants

Recently, the binding kinetics of a ligand-target interaction, such as the residence time of a small molecule on its protein target, are seen as increasingly important for drug efficacy. Here, we investigate these concepts to explain binding and proton blockage of rimantadine variants bearing progressively larger alkyl groups to influenza A virus M2 wild type (WT) and M2 S31N protein proton channel. We showed that resistance of M2 S31N to rimantadine analogues compared to M2 WT resulted from their higher k_off rates compared to the k_on rates according to electrophysiology (EP) measurements. This is due to the fact that, in M2 S31N, the loss of the V27 pocket for the adamantyl cage resulted in low residence time inside the M2 pore. Both rimantadine enantiomers have similar channel blockage and binding k_on and k_off against M2 WT. To compare the potency between the rimantadine variants against M2, we applied approaches using different mimicry of M2, i.e., isothermal titration calorimetry and molecular dynamics simulation, EP, and antiviral assays. It was also shown that a small change in an amino acid at site 28 of M2 WT, which does not line the pore, seriously affects M2 WT blockage kinetics.

Molecular Dynamics Simulations on the Bioactive Molecule of hIAPP22-29 (NFGAILSS) and Rational Drug Design

This chapter includes information about the structure in equilibrium of the bioactive molecule hIAPP22-29 (NFGAILSS). The experimental structure was derived using X-ray and its 2D NOESY NMR experiments in d ₆-DMSO and d-HFIP solvents. This molecule contains eight of the ten amino acids of the 20-29 region of the human islet amyloid polypeptide (hIAPP) often referred as the "amyloidogenic core." Amyloid deposits are well-known to cause as many as 20 pathological neurodegenerative disorders such as Alzheimer, Parkinson, Huntington, and Creutzfeldt-Jakob. The experimental structure was relaxed using molecular dynamics (MD) in simulation boxes consisting in DMSO and HFIP; the latter not provided by the applied software. The calculations were performed in GPUs and supercomputers, and some basic scripting is described for reference. The simulations confirmed the inter- and intramolecular forces that led to an "amyloidogenic core" observed from NOE experiments. The results showed that in DMSO and HFIP environment, Phe is not in spatial proximity with Leu or Ile, and this is consistent with an amyloidogenic core. However, in an amphipathic environment such as the model lipid bilayers, this communication is possible and may influence peptide amyloidogenic properties. The knowledge gained through this study may contribute to the rational drug design of novel peptides or organic molecules acting by modifying preventing amyloidogenic properties of the hIAPP peptide.

Governing effects in the mechanism of the gold-catalyzed cycloisomerization of allenic hydroxylamine derivatives

The formation of chiral heterocycles via cycloisomerization reactions of allene derivatives has gained relevance due to their associated efficiency and atom-economy. The only drawback that keeps these reactions away from being routine synthetic strategies is the control in the regioselectivity (most often 5-endo vs. 6-endo). In this work, we computationally explore the experimental chemistry reported by Krause using N-hydroxy-α-aminoallenes and hydroxylamine ethers as substrates and provide a rationale for the different reactivity observed. The drastic effects observed experimentally when changing the nature of the gold catalyst have also been studied mechanistically. These results are expected to help in the design of improved regioselective protocols for the formation of medium sized chiral heterocycles from allene substrates.

Affinity of Rimantadine Enantiomers against Influenza A/M2 Protein Revisited

Recent findings from solid state NMR (ssNMR) studies suggested that the (R)-enantiomer of rimantadine binds to the full M2 protein with higher affinity than the (S)-enantiomer. Intrigued by these findings, we applied functional assays, such as antiviral assay and electrophysiology (EP), to evaluate the binding affinity of rimantadine enantiomers to the M2 protein channel. Unexpectedly, no significant difference was found between the two enantiomers. Our experimental data based on the full M2 protein function were further supported by alchemical free energy calculations and isothermal titration calorimetry (ITC) allowing an evaluation of the binding affinity of rimantadine enantiomers to the M2TM pore. Both enantiomers have similar channel blockage, affinity, and antiviral potency.