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.

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.

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.

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.

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.

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.

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.

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.

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ß3 receptors on cancer cells for molecular imaging and diagnostic applications using positron emission tomography (PET). Integrin αvß3 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ß3 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 (18F, 64Cu, and 68Ga), 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.

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 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 (Ala25 ). 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 Ala25 . DFT calculations suggest that folding from the extended structure to the α-helical structure is favored for M2TM compared with Ala25 . 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 Ala25 . 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.

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.

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 13C cross-polarization/magic angle spinning nuclear magnetic resonance (CP/MAS NMR) spectroscopy. The 13C 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.

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 A3 receptor (A3R) 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 A3R. 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 A3R. MD simulations and MM-GBSA calculations of the WT A3R 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 A3R is defined not only by direct interactions with residues within the orthosteric binding area but also by remote residues playing a significant role. Although V1695.30 is considered to be a selectivity filter for A3R 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 L903.32 in the low region and the remote L2647.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 A3R, which we previously performed for agonists binding to A3R, and the design of more effective antagonists based on K18.

Effects of Cholesterol on GPCR Function: Insights from Computational and Experimental Studies.

The extensive experimental and computational evidences revealed that cholesterol is involved in the drug binding to G protein-coupled receptor (GPCR) targets that is influenced by the membrane environment and external functions. These multifunctional factors make the understanding of the molecular mechanism of action in greater detail an entirely difficult task. Significant efforts have been made for better understanding the role of multi-directional specific, receptor-dependent interactions of cholesterol, and its effects on drug design and development. Additional efforts must be made in this complex system in order to shed more light on cholesterol molecular basis of action. The results of molecular simulations that complemented experimental data may reveal new aspects of GPCR-cholesterol interactions and may provide a comprehensive understanding of receptor function.

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 Inwardclosed conformation (2.00 Å), rimantadine bound to M2 in both the Inwardclosed (2.00 Å) and Inwardopen (2.25 Å) conformations, and a spiro-adamantyl amine inhibitor bound to M2 in the Inwardclosed 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 A2A and A1ARs 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 A2AAR 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 A2AAR, the identified antagonists exhibited low micromolar or micromolar affinity against A2A/A3, ARs, or A3AR, 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 (A2A/A3, A1/A3, A1/A2A/A3, and A3). 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 A3 receptor.

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.

Assessing the attractive/repulsive force balance in axial cyclohexane C-Hax ···Yax contacts: A combined computational analysis in monosubstituted cyclohexanes.

The interactions of axial substituents in monosubstituted cyclohexane rings are studied in this work using an array of different computational techniques. Additionally, the anomalous axial preference for some bulky substituents is related to stabilizing dispersion interactions. We find that the C-Hax ···Yax contacts for various substituents with distances ranging from 2 to ∼5 Å may include attractive dispersion forces that can affect the conformational equilibrium; these forces co-exist with Pauli repulsive forces effected by Yax group due to van der Waals sphere penetration. At distances between 2 and 3 Å stabilizing electron transfer interactions were calculated and the combination of natural bond orbital and QTAIM analysis showed that, in certain cases, Yax = t Bu, Cax -O or Cax = O or Sax = O or Cax = S this interaction can be characterized as an improper H-bond. DFT-D3 and non-covalent interactions calculations (NCIs) in cyclohexane derivatives with Yax = SiOR3 including HYax ···Hcy surfaces at distances ranging between 4 and 6 Å suggest that dispersion has a clear effect on the experimentally observed stabilization of the axial conformer. NCIs computed from the reduced density gradient help to visually identify and analyze these interactions. © 2016 Wiley Periodicals, Inc.