Strengthening an Intramolecular Non-Classical Hydrogen Bond to Get in Shape for Binding

In this research article, we report on the strengthening of a non-classical hydrogen bond (C-H⋅⋅⋅O) by introducing electron withdrawing groups at the carbon atom. The approach is demonstrated on the example of derivatives of the physiological E-selectin ligand sialyl Lewisx (1, sLex). Its affinity is mainly due to a beneficial entropy term, which is predominantly caused by the pre-organization of sLex in its binding conformation. We have shown, that among the elements responsible for the pre-organization, the stabilization by a non-classical hydrogen bond between the H-C5 of l-fucose and the ring oxygen O5 of the neighboring d-galactose moiety is essential and yields 7.4 kJ mol-1. This effect could be further strengthened by replacing l-fucose by 6,6,6-trifluoro-l-fucose leading to an improved non-classical H-bond of 14.9 kJ mol-1, i.e., an improved pre-organization in the bioactive conformation. For a series of glycomimetics of sLex (1), this outcome could be confirmed by high field NMR-shifts of the H-C5Fuc, by X-ray diffraction analysis of glycomimetics co-crystallized with E-selectin as well as by isothermal titration calorimetry. Furthermore, the electron-withdrawing character of the CF3-group beneficially influences the pharmacokinetic properties of sLex mimetics. Thus, acid-stability, a prerequisite for gastrointestinal stability, could be substantially improved.

Evaluation of Anticholinesterase Activity of the Fungicides Mefentrifluconazole and Pyraclostrobin

Triazoles are compounds with various biological activities, including fungicidal action. They became popular through cholinesterase studies after the successful synthesis of the dual binding femtomolar triazole inhibitor of acetylcholinesterase (AChE, EC 3.1.1.7) by Sharpless et al. via in situ click chemistry. Here, we evaluate the anticholinesterase effect of the first isopropanol triazole fungicide mefentrifluconazole (Ravystar®), developed to overcome fungus resistance in plant disease management. Mefentrifluconazole is commercially available individually or in a binary fungicidal mixture, i.e., with pyraclostrobin (Ravycare®). Pyraclostrobin is a carbamate that contains a pyrazole ring. Carbamates are known inhibitors of cholinesterases and the carbamate rivastigmine is already in use for the treatment of Alzheimer's disease. We tested the type and potency of anticholinesterase activity of mefentrifluconazole and pyraclostrobin. Mefentrifluconazole reversibly inhibited human AChE and BChE with a seven-fold higher potency toward AChE (Ki = 101 ± 19 μM). Pyraclostrobin (50 μM) inhibited AChE and BChE progressively with rate constants of (t1/2 = 2.1 min; ki = 6.6 × 103 M-1 min-1) and (t1/2 = 1.5 min; ki = 9.2 × 103 M-1 min-1), respectively. A molecular docking study indicated key interactions between the tested fungicides and residues of the lipophilic active site of AChE and BChE. Additionally, the physicochemical properties of the tested fungicides were compared to values for CNS-active drugs to estimate the blood-brain barrier permeability. Our results can be applied in the design of new molecules with a lesser impact on humans and the environment.

Trialkylphosphonium Oxoborate as C(sp3 )-H Oxyanion Hole Catalyst for Diels-Alder Reaction

We have developed a catalytic protocol for Diels-Alder reaction using trialkylphosphonium oxoborates as oxyanion hole catalysts. The reaction can be operated under ambient conditions. Dienes could easily polymerize under acidic condition. Nonetheless, these acid-sensitive substrates are compatible with the catalytic protocol and the reaction scope covers a wide range of substrates.

Trialkylphosphonium oxoborates as C(sp3)-H oxyanion holes and their application in catalytic chemoselective acetalization

The use of trialkylphosphonium oxoborates (TOB) as catalysts is reported. The site-isolated borate counter anion in a TOB catalyst increases the availability of C(sp3)-H to interact with electron donor substrates. The catalytic protocol is applicable to a wide range of substrates in the acetalization reaction and provides excellent chemoselectivity in the acetalization over thioacetalization in the presence of alcohols and thiols, which is otherwise hard to achieve using typical acid catalysts. Experimental and computational studies revealed that the TOB catalysts have multiple preorganized C(sp3)-Hs that serve as a mimic of oxyanion holes, which can stabilize the oxyanion intermediates via multiple C(sp3)-H non-classical hydrogen bond interactions.

Intra- and Intermolecular Cation-π Interactions between Onium Salts and Alkynes/Acetylene: Experimental and Theoretical Insights

Cation-π interactions between various onium salts, alkynes, and acetylene were studied, taking into account the substituents of the triple bond, the nature of the anions, and the polarity of the solvent, through a combination of MP2 calculations and experiments. In an intramolecular setting, these data (including single-crystal X-ray crystallography) concurred with the stability of folded conformers of alkynyl onium salts, even substituted with electron-withdrawing groups. To examine the contribution of these interactions on the alkyne electronic population, a thorough in silico study was carried out using natural bonding orbital analysis of the conformers. Intramolecular interactions from sulfonium salt tethered to phenylalkyne were highlighted, as illustrated above by the computed folded conformation (MP2) along with noncovalent interaction (NCI) analysis. Furthermore, investigations of intermolecular interactions, involving acetylene or phenylacetylene with various onium ions, revealed the high energy interactions of their complexes with phenyldimethylsulfonium chloride, as illustrated above with the complex PhC≡CH/PhMe2SCl (MP2 calculations and NCI analysis).

A Catalyst-Controlled Enantiodivergent Bromolactonization

A catalyst-controlled enantiodivergent bromolactonization of olefinic acids has been developed. Quinine-derived amino-amides bearing the same chiral core but different achiral aryl substituents were used as the catalysts. Switching the methoxy substituent in the aryl amide system from meta- to ortho-position results in a complete switch in asymmetric induction to afford the desired lactone in good enantioselectivity and yield. Mechanistic studies, including chemical experiments and density functional theory calculations, reveal that the differences in steric and electronic effects of the catalyst substituent alter the reaction mechanism.

Zwitterion-Catalyzed Isomerization of Maleic to Fumaric Acid Diesters

Fumaric acid diesters are important building blocks for organic synthesis. A class of zwitterionic organocatalysts based on an amide anion/iminium cation charge pair were found to be effective in catalyzing the isomerization of maleic acid diesters to give fumaric acid diesters. Comparison of the performance of different zwitterionic organocatalysts toward the reaction revealed that nonclassical hydrogen bonding was involved in the stabilization of the Michael adduct intermediate.

Physical Mechanisms Governing Substituent Effects on Arene-Arene Interactions in a Protein Milieu

Arene-arene interactions play important roles in protein-ligand complex formation. Here, we investigate the characteristics of arene-arene interactions between small organic molecules and aromatic amino acids in protein interiors. The study is based on X-ray crystallographic data and quantum mechanical calculations using the enzyme acetylcholinesterase and selected inhibitory ligands as a model system. It is shown that the arene substituents of the inhibitors dictate the strength of the interaction and the geometry of the resulting complexes. Importantly, the calculated interaction energies correlate well with the measured inhibitor potency. Non-hydrogen substituents strengthened all interaction types in the protein milieu, in keeping with results for benzene dimer model systems. The interaction energies were dispersion-dominated, but substituents that induced local dipole moments increased the electrostatic contribution and thus yielded more strongly bound complexes. These findings provide fundamental insights into the physical mechanisms governing arene-arene interactions in the protein milieu and thus into molecular recognition between proteins and small molecules.

Zwitterion-Catalyzed Intermolecular Bromoesterifications

Intermolecular haloesterification is an important class of transformations. The resulting products are valuable building blocks. However, it is often necessary to use super-stoichiometric amount of acid in order to compensate the low reactivity. Herein, we report a zwitterion-catalyzed intermolecular bromoesterification using acid and olefin in an equimolar ratio. Mechanistic study revealed that the charge pair in the zwitterion works synergistically in activating both NBS and carboxylic acid.

Sagrillo F, R. Freitas MC, A. Alves M, Thota S, Tinoco LW, Magalhães A, Sant’Anna CMR, Fraga CAM. Structure–property relationship studies of 3-acyl-substituted furans: the serendipitous identification and characterization of a new non-classical hydrogen bond donor moiety

A serendipitous identification and characterization of a new non-classical hydrogen bond donor moiety found in N-acylhydrazones containing 3-acyl-substituted furan subunit is presented.

Structural insights into the anti-cancer activity of quercetin on G-tetrad, mixed G-tetrad, and G-quadruplex DNA using quantum chemical and molecular dynamics simulations

Human telomerase referred as 'terminal transferase' is a nucleoprotein enzyme which inhibits the disintegration of telomere length and act as a drug target for the anticancer therapy. The tandem repeating structure of telomere sequence forms the guanine-rich quadruplex structures that stabilize stacked tetrads. In our present work, we have investigated the interaction of quercetin with DNA tetrads using DFT. Geometrical analysis revealed that the influence of quercetin drug induces the structural changes into the DNA tetrads. Among DNA tetrads, the quercetin stacked with GCGC tetrad has the highest interaction energy of -88.08 kcal/mol. The binding mode and the structural stability are verified by the absorption spectroscopy method. The longer wavelength was found at 380 nm and it exhibits bathochromic shift. The findings help us to understand the binding nature of quercetin drug with DNA tetrads and it also inhibits the telomerase activity. Further, the quercetin drug interacted with G-quadruplex DNA by using molecular dynamics (MD) simulation studies for 100 ns simulation at different temperatures and different pH levels (T = 298 K, 320 K and pH = 7.4, 5.4). The structural stability of the quercetin with G-quadruplex structure is confirmed by RMSD. For the acidic condition (pH = 5.4), the binding affinity is higher toward G-quadruplex DNA, this result resembles that the quercetin drug is well interacted with G-quadruplex DNA at acidic condition (pH = 7.4) than the neutral condition. The obtained results show that quercetin drug stabilizes the G-quadruplex DNA, which regulates telomerase enzyme and it potentially acts as a novel anti-cancer agent.Communicated by Ramaswamy H. Sarma.

N+-C-H···O Hydrogen bonds in protein-ligand complexes

In the context of drug design, C-H···O hydrogen bonds have received little attention so far, mostly because they are considered weak relative to other noncovalent interactions such as O-H···O hydrogen bonds, π/π interactions, and van der Waals interactions. Herein, we demonstrate the significance of hydrogen bonds between C-H groups adjacent to an ammonium cation and an oxygen atom (N⁺-C-H···O hydrogen bonds) in protein-ligand complexes. Quantum chemical calculations revealed details on the strength and geometrical requirements of these N⁺-C-H···O hydrogen bonds, and a subsequent survey of the Protein Data Bank (PDB) based on these criteria suggested that numerous protein-ligand complexes contain such N⁺-C-H···O hydrogen bonds. An ensuing experimental investigation into the G9a-like protein (GLP)-inhibitor complex demonstrated that N⁺-C-H···O hydrogen bonds affect the activity of the inhibitors against the target enzyme. These results should provide the basis for the use of N⁺-C-H···O hydrogen bonds in drug discovery.

Noncovalent Inhibitors of Mosquito Acetylcholinesterase 1 with Resistance-Breaking Potency

Resistance development in insects significantly threatens the important benefits obtained by insecticide usage in vector control of disease-transmitting insects. Discovery of new chemical entities with insecticidal activity is highly desired in order to develop new insecticide candidates. Here, we present the design, synthesis, and biological evaluation of phenoxyacetamide-based inhibitors of the essential enzyme acetylcholinesterase 1 (AChE1). AChE1 is a validated insecticide target to control mosquito vectors of, e.g., malaria, dengue, and Zika virus infections. The inhibitors combine a mosquito versus human AChE selectivity with a high potency also for the resistance-conferring mutation G122S; two properties that have proven challenging to combine in a single compound. Structure-activity relationship analyses and molecular dynamics simulations of inhibitor-protein complexes have provided insights that elucidate the molecular basis for these properties. We also show that the inhibitors demonstrate in vivo insecticidal activity on disease-transmitting mosquitoes. Our findings support the concept of noncovalent, selective, and resistance-breaking inhibitors of AChE1 as a promising approach for future insecticide development.

Structure-Based Design of Inhibitors Targeting PrfA, the Master Virulence Regulator of Listeria monocytogenes

Listeria monocytogenes is a bacterial pathogen that controls much of its virulence through the transcriptional regulator PrfA. In this study, we describe structure-guided design and synthesis of a set of PrfA inhibitors based on ring-fused 2-pyridone heterocycles. Our most effective compound decreased virulence factor expression, reduced bacterial uptake into eukaryotic cells, and improved survival of chicken embryos infected with L. monocytogenes compared to previously identified compounds. Crystal structures identified an intraprotein "tunnel" as the main inhibitor binding site (A_I), where the compounds participate in an extensive hydrophobic network that restricts the protein's ability to form functional DNA-binding helix-turn-helix (HTH) motifs. Our studies also revealed a hitherto unsuspected structural plasticity of the HTH motif. In conclusion, we have designed 2-pyridone analogues that function as site-A_I selective PrfA inhibitors with potent antivirulence properties.

Ranking Power of the SQM/COSMO Scoring Function on Carbonic Anhydrase II-Inhibitor Complexes

Accurate prediction of protein-ligand binding affinities is essential for hit-to-lead optimization and virtual screening. The reliability of scoring functions can be improved by including quantum effects. Here, we demonstrate the ranking power of the semiempirical quantum mechanics (SQM)/implicit solvent (COSMO) scoring function by using a challenging set of 10 inhibitors binding to carbonic anhydrase II through Zn²⁺ in the active site. This new dataset consists of the high-resolution (1.1-1.4 Å) crystal structures and experimentally determined inhibitory constant (K_i ) values. It allows for evaluation of the common approximations, such as representing the solvent implicitly or by using a single target conformation combined with a set of ligand docking poses. SQM/COSMO attained a good correlation of R² of 0.56-0.77 with the experimental inhibitory activities, benefiting from careful handling of both noncovalent interactions (e.g. charge transfer) and solvation. This proof-of-concept study of SQM/COSMO ranking for metalloprotein-ligand systems demonstrates its potential for hit-to-lead applications.

Comparison of the Binding of Reversible Inhibitors to Human Butyrylcholinesterase and Acetylcholinesterase: A Crystallographic, Kinetic and Calorimetric Study

Acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) hydrolyze the neurotransmitter acetylcholine and, thereby, function as coregulators of cholinergic neurotransmission. Although closely related, these enzymes display very different substrate specificities that only partially overlap. This disparity is largely due to differences in the number of aromatic residues lining the active site gorge, which leads to large differences in the shape of the gorge and potentially to distinct interactions with an individual ligand. Considerable structural information is available for the binding of a wide diversity of ligands to AChE. In contrast, structural data on the binding of reversible ligands to BChE are lacking. In a recent effort, an inhibitor competition approach was used to probe the overlap of ligand binding sites in BChE. Here, we extend this study by solving the crystal structures of human BChE in complex with five reversible ligands, namely, decamethonium, thioflavin T, propidium, huprine, and ethopropazine. We compare these structures to equivalent AChE complexes when available in the protein data bank and supplement this comparison with kinetic data and observations from isothermal titration calorimetry. This new information now allows us to define the binding mode of various ligand families and will be of importance in designing specific reversible ligands of BChE that behave as inhibitors or reactivators.

New pyridine derivatives as inhibitors of acetylcholinesterase and amyloid aggregation

A new series of pyridine derivatives with carbamic or amidic function has been designed and synthesized to act as cholinesterase inhibitors. The synthesized compounds were tested toward EeAChE and hAChE and toward eqBChE and hBChE. The carbamate 8 was the most potent hAChE inhibitor (IC₅₀ = 0.153 ± 0.016 μM) while the carbamate 11 was the most potent inhibitor of hBChE (IC₅₀ = 0.828 ± 0.067 μM). A molecular docking study indicated that the carbamate 8 was able to bind AChE by interacting with both CAS and PAS, in agreement with the mixed inhibition mechanism. Furthermore, the carbamates 8, 9 and 11 were able to inhibit Aβ₄₂ self-aggregation and possessed quite low toxicity against human astrocytoma T67 and HeLa cell lines, being the carbamate 8 the less toxic compound on both cell lines.

An Unusual Dimeric Inhibitor of Acetylcholinesterase: Cooperative Binding of Crystal Violet

Acetylcholinesterase (AChE) is an essential enzyme that terminates cholinergic transmission by a rapid hydrolysis of the neurotransmitter acetylcholine. AChE is an important target for treatment of various cholinergic deficiencies, including Alzheimer’s disease and myasthenia gravis. In a previous high throughput screening campaign, we identified the dye crystal violet (CV) as an inhibitor of AChE. Herein, we show that CV displays a significant cooperativity for binding to AChE, and the molecular basis for this observation has been investigated by X-ray crystallography. Two monomers of CV bind to residues at the entrance of the active site gorge of the enzyme. Notably, the two CV molecules have extensive intermolecular contacts with each other and with AChE. Computational analyses show that the observed CV dimer is not stable in solution, suggesting the sequential binding of two monomers. Guided by the structural analysis, we designed a set of single site substitutions, and investigated their effect on the binding of CV. Only moderate effects on the binding and the cooperativity were observed, suggesting a robustness in the interaction between CV and AChE. Taken together, we propose that the dimeric cooperative binding is due to a rare combination of chemical and structural properties of both CV and the AChE molecule itself.

Conformational energy range of ligands in protein crystal structures: The difficult quest for accurate understanding

In this review, we address a fundamental question: What is the range of conformational energies seen in ligands in protein-ligand crystal structures? This value is important biophysically, for better understanding the protein-ligand binding process; and practically, for providing a parameter to be used in many computational drug design methods such as docking and pharmacophore searches. We synthesize a selection of previously reported conflicting results from computational studies of this issue and conclude that high ligand conformational energies really are present in some crystal structures. The main source of disagreement between different analyses appears to be due to divergent treatments of electrostatics and solvation. At the same time, however, for many ligands, a high conformational energy is in error, due to either crystal structure inaccuracies or incorrect determination of the reference state. Aside from simple chemistry mistakes, we argue that crystal structure error may mainly be because of the heuristic weighting of ligand stereochemical restraints relative to the fit of the structure to the electron density. This problem cannot be fixed with improvements to electron density fitting or with simple ligand geometry checks, though better metrics are needed for evaluating ligand and binding site chemistry in addition to geometry during structure refinement. The ultimate solution for accurately determining ligand conformational energies lies in ultrahigh-resolution crystal structures that can be refined without restraints.

Superior Performance of the SQM/COSMO Scoring Functions in Native Pose Recognition of Diverse Protein-Ligand Complexes in Cognate Docking

General and reliable description of structures and energetics in protein-ligand (PL) binding using the docking/scoring methodology has until now been elusive. We address this urgent deficiency of scoring functions (SFs) by the systematic development of corrected semiempirical quantum mechanical (SQM) methods, which correctly describe all types of noncovalent interactions and are fast enough to treat systems of thousands of atoms. Two most accurate SQM methods, PM6-D3H4X and SCC-DFTB3-D3H4X, are coupled with the conductor-like screening model (COSMO) implicit solvation model in so-called "SQM/COSMO" SFs and have shown unique recognition of native ligand poses in cognate docking in four challenging PL systems, including metalloprotein. Here, we apply the two SQM/COSMO SFs to 17 diverse PL complexes and compare their performance with four widely used classical SFs (Glide XP, AutoDock4, AutoDock Vina, and UCSF Dock). We observe superior performance of the SQM/COSMO SFs and identify challenging systems. This method, due to its generality, comparability across the chemical space, and lack of need for any system-specific parameters, gives promise of becoming, after comprehensive large-scale testing in the near future, a useful computational tool in structure-based drug design and serving as a reference method for the development of other SFs.

Structure of a prereaction complex between the nerve agent sarin, its biological target acetylcholinesterase, and the antidote HI-6

Organophosphorus nerve agents interfere with cholinergic signaling by covalently binding to the active site of the enzyme acetylcholinesterase (AChE). This inhibition causes an accumulation of the neurotransmitter acetylcholine, potentially leading to overstimulation of the nervous system and death. Current treatments include the use of antidotes that promote the release of functional AChE by an unknown reactivation mechanism. We have used diffusion trap cryocrystallography and density functional theory (DFT) calculations to determine and analyze prereaction conformers of the nerve agent antidote HI-6 in complex with Mus musculus AChE covalently inhibited by the nerve agent sarin. These analyses reveal previously unknown conformations of the system and suggest that the cleavage of the covalent enzyme-sarin bond is preceded by a conformational change in the sarin adduct itself. Together with data from the reactivation kinetics, this alternate conformation suggests a key interaction between Glu202 and the O-isopropyl moiety of sarin. Moreover, solvent kinetic isotope effect experiments using deuterium oxide reveal that the reactivation mechanism features an isotope-sensitive step. These findings provide insights into the reactivation mechanism and provide a starting point for the development of improved antidotes. The work also illustrates how DFT calculations can guide the interpretation, analysis, and validation of crystallographic data for challenging reactive systems with complex conformational dynamics.