Evidence for CCl/CBr⋅⋅⋅π Interactions as an Important Contribution to Protein-Ligand Binding Affinity

Metacaspase (Pf MCA-1) as antimalarial drug target: An in silico approach and their biological validation

AIMS

Acquired drug resistance of Plasmodium is a global issue for the treatment of malaria. There are various proteases in the genome of Plasmodium falciparum (P. falciparum) including metacaspase-1 (PfMCA-1) that are essential and are being considered as an attractive drug target. It is aimed to identify novel therapeutics against malaria and their action on PfMCA-1 along with other apoptotic pathway events.

MAIN METHODS

High throughput virtual screening of 55,000 compounds derived from Maybridge library was performed against PfMCA-1. Based on the docking score, sixteen compounds were selected for in vitro antimalarial screening against drug sensitive and resistant strains of P. falciparum using SYBR green-based assay. Subsequently, three lead molecules were selected and subjected to the evaluation of cytotoxicity, caspase like protease activity, mitochondrial membrane potential, ROS generation and DNA fragmentation via TUNEL assay.

KEY FINDINGS

The in silico and in vitro approaches have brought forward some Maybridge library compounds with antiplasmodial activity most likely by enhancing the metacaspase activity. The compound CD11095 has shown better antimalarial efficacy, and KM06591 depicted higher caspase mediated killing, elevated TUNEL positive cells and moderate ROS generation. Mitochondrial membrane depolarization was augmented by RJC0069. Exposure of P. falciparum to CD11095, KM06591 and RJC0069 has ended up in parasite growth arrest via multiple mechanisms.

SIGNIFICANCE

It is proposed that the Maybridge molecules CD11095, KM06591 and RJC0069 have antimalarial activity. Their mechanism of action was found to be by enhancing the metacaspases-like protease activity, mitochondrial depolarization and DNA fragmentation which stipulates significant insights towards promising candidates for drug development.

Automated De Novo Design in Medicinal Chemistry: Which Types of Chemistry Does a Generative Neural Network Learn?

Artificial intelligence offers promising solutions for property prediction, compound design, and retrosynthetic planning, which are expected to significantly accelerate the search for pharmacologically relevant molecules. Here, we investigate aspects of artificial intelligence based de novo design pertaining to its integration into real-life workflows. First, different chemical spaces were used as training sets for reinforcement learning (RL) in combination with different reward functions. With the trained neuronal networks different biologically active molecules could be regenerated. Excluding molecules with substructures such as five-membered rings from training spaces nevertheless produced results containing these moieties. Furthermore, different scoring functions in RL were investigated and produced different design ensembles. In summary, some of these design proposals are close in chemical space to the query, thus supporting lead optimization, while 3D-shape or QSAR (quantitative structure-activity relationship) models produced significantly different proposals by sampling a broader region of the chemical space, thus supporting lead generation. Therefore, RL provides a good framework to tailored design approaches for different discovery phases.

Characterizing hydration sites in protein-ligand complexes towards the design of novel ligands

Water is an essential part of protein binding sites and mediates interactions to ligands. Its displacement by ligand parts affects the free binding energy of resulting protein-ligand complexes. Therefore the characterization of solvation properties is important for design. Of particular interest is the propensity of localized water to be favorably displaced by a ligand. This review discusses two popular computational approaches addressing these questions, namely WaterMap based on statistical mechanics analysis of MD simulations and 3D RISM based on integral equation theory of liquids. The theoretical background and recent applications in structure-based design will be presented.

Direct Experimental Evidence for Halogen-Aryl π Interactions in Solution from Molecular Torsion Balances

We dissected halogen-aryl π interactions experimentally using a bicyclic N-arylimide based molecular torsion balances system, which is based on the influence of the non-bonded interaction on the equilibria between folded and unfolded states. Through comparison of balances modulated by higher halogens with fluorine balances, we determined the magnitude of the halogen-aryl π interactions in our unimolecular systems to be larger than -5.0 kJ mol^-1 , which is comparable with the magnitude estimated in the biomolecular systems. Our study provides direct experimental evidence of halogen-aryl π interactions in solution, which until now have only been revealed in the solid state and evaluated theoretically by quantum-mechanical calculations.

The Important Role of Halogen Bond in Substrate Selectivity of Enzymatic Catalysis

The use of halogen bond is widespread in drug discovery, design, and clinical trials, but is overlooked in drug biosynthesis. Here, the role of halogen bond in the nitrilase-catalyzed synthesis of ortho-, meta-, and para-chlorophenylacetic acid was investigated. Different distributions of halogen bond induced changes of substrate binding conformation and affected substrate selectivity. By engineering the halogen interaction, the substrate selectivity of the enzyme changed, with the implication that halogen bond plays an important role in biosynthesis and should be used as an efficient and reliable tool in enzymatic drug synthesis.

Non-covalent C-Cl…π interaction in acetylene-carbon tetrachloride adducts: Matrix isolation infrared and ab initio computational studies

Non-covalent halogen-bonding interactions between π cloud of acetylene (C2H2) and chlorine atom of carbon tetrachloride (CCl4) have been investigated using matrix isolation infrared spectroscopy and quantum chemical computations. The structure and the energies of the 1:1 C2H2-CCl4 adducts were computed at the B3LYP, MP2 and M05-2X levels of theory using 6-311++G(d,p) basis set. The computations indicated two minima for the 1:1 C2H2-CCl4 adducts; with the C-Cl…π adduct being the global minimum, where π cloud of C2H2 is the electron donor. The second minimum corresponded to a C-H…Cl adduct, in which C2H2 is the proton donor. The interaction energies for the adducts A and B were found to be nearly identical. Experimentally, both C-Cl…π and C-H…Cl adducts were generated in Ar and N2 matrixes and characterized using infrared spectroscopy. This is the first report on halogen bonded adduct, stabilized through C-Cl…π interaction being identified at low temperatures using matrix isolation infrared spectroscopy. Atoms in Molecules (AIM) and Natural Bond Orbital (NBO) analyses were performed to support the experimental results. The structures of 2:1 ((C2H2)2-CCl4) and 1:2 (C2H2-(CCl4)2) multimers and their identification in the low temperature matrixes were also discussed.

Halogen and Hydrogen Bonding Benzothiophene Diol Derivatives: A Study Using ab initio Calculations and X-Ray Crystal Structure Measurements

The aim of this study is to describe and compare the supramolecular interactions, in the solid state, of chloro-, bromo-, and iodobenzothiophene diols. The compounds were obtained through organo-catalyzed reactions starting from 3-substituted halobenzothiophene carbaldehydes. Energies of the noncovalent interactions were obtained by density functional theory calculations. Bond distances and angles were found to be in accordance with those determined by X-ray structure analysis. anti-Bromobenzothiophene derivatives showed strong halogen⋅⋅⋅π interactions between bromine and the heterocyclic phenyl ring, corresponding to an energy of 7.5 kcal mol(-1). syn-Bromo and syn-iodo derivatives appeared to be isostructural, showing X⋅⋅⋅O (carbonyl) interactions, π stacking, and formation of extended hydrogen bonding networks. In contrast, the chloro derivatives displayed no halogen bonding interactions.

Molecular recognition in chemical and biological systems

Structure-based ligand design in medicinal chemistry and crop protection relies on the identification and quantification of weak noncovalent interactions and understanding the role of water. Small-molecule and protein structural database searches are important tools to retrieve existing knowledge. Thermodynamic profiling, combined with X-ray structural and computational studies, is the key to elucidate the energetics of the replacement of water by ligands. Biological receptor sites vary greatly in shape, conformational dynamics, and polarity, and require different ligand-design strategies, as shown for various case studies. Interactions between dipoles have become a central theme of molecular recognition. Orthogonal interactions, halogen bonding, and amide⋅⋅⋅π stacking provide new tools for innovative lead optimization. The combination of synthetic models and biological complexation studies is required to gather reliable information on weak noncovalent interactions and the role of water.

‘Honeycomb’ nanotube assembly based on thiacalix[4]arene derivatives by weak interactions

Crystallisation of six thiacalix[4]arene derivatives from hexane–chloroform leads to ‘honeycomb’ nanotube architectures and each tubular stack is surrounded by six close tubular neighbours via weak interactions, such as S⋯π interactions, C–H⋯π interactions, and so on.

Highly enantioselective synthesis of 1,2,3-substituted cyclopropanes by using α-Iodo- and α-chloromethylzinc carbenoids

Herein, we report the enantio- and diastereoselective formation of trans-iodo- and trans-chlorocyclopropanes from α-iodo- and α-chlorozinc carbenoids by using a dioxaborolane-derived chiral ligand. The synthetically useful iodocyclopropane building blocks were derivatized by an electrophilic trapping of the corresponding cyclopropyl lithium species or a Negishi coupling to give access to a variety of enantioenriched 1,2,3-substituted cyclopropanes. The synthetic utility of this method was demonstrated by the formal synthesis of an HIV-1 protease inhibitor. In addition, the related stereoselective bromocyclopropanation was also investigated. New insights about the relative electrophilicity of haloiodomethylzinc carbenoids are also presented.

High-resolution crystal structures of factor XIa coagulation factor in complex with nonbasic high-affinity synthetic inhibitors

Factor XI (FXI) is a key enzyme in the coagulation pathway and an attractive target for the development of anticoagulant drugs. A small number of high-resolution crystal structures of FXIa in complex with small synthetic inhibitors have been published to date. All of these ligands have a basic P1 group and bind exclusively in the nonprime side of the active site of FXIa. Here, two structures of FXIa in complex with nonbasic inhibitors that occupy both the prime and nonprime sides of the active site are presented. These new structures could be valuable in the design and optimization of new FXIa synthethic inhibitors.

Molecular recognition at the active site of factor Xa: cation-π interactions, stacking on planar peptide surfaces, and replacement of structural water

Factor Xa, a serine protease from the blood coagulation cascade, is an ideal enzyme for molecular recognition studies, as its active site is highly shape-persistent and features distinct, concave sub-pockets. We developed a family of non-peptidic, small-molecule inhibitors with a central tricyclic core orienting a neutral heterocyclic substituent into the S1 pocket and a quaternary ammonium ion into the aromatic box in the S4 pocket. The substituents were systematically varied to investigate cation-π interactions in the S4 pocket, optimal heterocyclic stacking on the flat peptide walls lining the S1 pocket, and potential water replacements in both the S1 and the S4 pockets. Structure-activity relationships were established to reveal and quantify contributions to the binding free enthalpy, resulting from single-atom replacements or positional changes in the ligands. A series of high-affinity ligands with inhibitory constants down to K(i)=2 nM were obtained and their proposed binding geometries confirmed by X-ray co-crystal structures of protein-ligand complexes.

Contribution of explicit solvent effects to the binding affinity of small-molecule inhibitors in blood coagulation factor serine proteases

The prevention of blood coagulation is important in treating thromboembolic disorders, and several serine proteases involved in the coagulation cascade have been classified as pharmaceutically relevant. Whereas structure-based drug design has contributed to the development of some serine protease inhibitors, traditional computational methods have not been able to fully describe structure-activity relationships (SAR). Here, we study the SAR for a number of serine proteases by using a method that calculates the thermodynamic properties (enthalpy and entropy) of the water that solvates the active site. We show that the displacement of water from specific subpockets (such as S1-4 and the ester binding pocket) of the active site by the ligand can govern potency, especially for cases in which small chemical changes (i.e., a methyl group or halogen) result in a substantial increase in potency. Furthermore, we describe how relative binding free energies can be estimated by combining the water displacement energy with complementary terms from an implicit solvent molecular mechanics description binding.