A comparative study of warheads for design of cysteine protease inhibitors.

The effects on potency of cruzain inhibition of replacing a nitrile group with alternative warheads were explored. The oxime was almost an order of magnitude more potent than the corresponding nitrile and has the potential to provide access to the prime side of the catalytic site. Dipeptide aldehydes and azadipeptide nitriles were found to be two orders of magnitude more potent cruzain inhibitors than the corresponding dipeptide nitriles although potency differences were modulated by substitution at P1 and P3. Replacement of the α methylene of a dipeptide aldehyde with cyclopropane led to a loss of potency of almost three orders of magnitude. The vinyl esters and amides that were characterized as reversible inhibitors were less potent than the corresponding nitrile by between one and two orders of magnitude.

The influence of hydrogen bonding on partition coefficients.

This Perspective explores how consideration of hydrogen bonding can be used to both predict and better understand partition coefficients. It is shown how polarity of both compounds and substructures can be estimated from measured alkane/water partition coefficients. When polarity is defined in this manner, hydrogen bond donors are typically less polar than hydrogen bond acceptors. Analysis of alkane/water partition coefficients in conjunction with molecular electrostatic potential calculations suggests that aromatic chloro substituents may be less lipophilic than is generally believed and that some of the effect of chloro-substitution stems from making the aromatic π-cloud less available to hydrogen bond donors. Relationships between polarity and calculated hydrogen bond basicity are derived for aromatic nitrogen and carbonyl oxygen. Aligned hydrogen bond acceptors appear to present special challenges for prediction of alkane/water partition coefficients and this may reflect 'frustration' of solvation resulting from overlapping hydration spheres. It is also shown how calculated hydrogen bond basicity can be used to model the effect of aromatic aza-substitution on octanol/water partition coefficients.

Expression, purification and kinetic characterization of His-tagged glyceraldehyde-3-phosphate dehydrogenase from Trypanosoma cruzi.

Trypanosomes are flagellated protozoa responsible for serious parasitic diseases that have been classified by the World Health Organization as tropical sicknesses of major importance. One important drug target receiving considerable attention is the enzyme glyceraldehyde-3-phosphate dehydrogenase from the protozoan parasite Trypanosoma cruzi, the causative agent of Chagas disease (T. cruzi Glyceraldehyde-3-phosphate dehydrogenase (TcGAPDH); EC 1.2.1.12). TcGAPDH is a key enzyme in the glycolytic pathway of T. cruzi and catalyzes the oxidative phosphorylation of D-glyceraldehyde-3-phosphate (G3P) to 1,3-bisphosphoglycerate (1,3-BPG) coupled to the reduction of oxidized nicotinamide adenine dinucleotide, (NAD(+)) to NADH, the reduced form. Herein, we describe the cloning of the T. cruzi gene for TcGAPDH into the pET-28a(+) vector, its expression as a tagged protein in Escherichia coli, purification and kinetic characterization. The His(6)-tagged TcGAPDH was purified by affinity chromatography. Enzyme activity assays for the recombinant His(6)-TcGAPDH were carried out spectrophotometrically to determine the kinetic parameters. The apparent Michaelis-Menten constant (K(M)(app)) determined for D-glyceraldehyde-3-phosphate and NAD(+) were 352±21 and 272±25 µM, respectively, which were consistent with the values for the untagged enzyme reported in the literature. We have demonstrated by the use of Isothermal Titration Calorimetry (ITC) that this vector modification resulted in activity preserved for a higher period. We also report here the use of response surface methodology (RSM) to determine the region of optimal conditions for enzyme activity. A quadratic model was developed by RSM to describe the enzyme activity in terms of pH and temperature as independent variables. According to the RMS contour plots and variance analysis, the maximum enzyme activity was at 29.1°C and pH 8.6. Above 37°C, the enzyme activity starts to fall, which may be related to previous reports that the quaternary structure begins a process of disassembly.

Similarity search combined with docking and molecular dynamics for novel hAChE inhibitor scaffolds.

The main purpose of this study was to address the performance of virtual screening methods based on ligands and the protein structure of acetylcholinesterase (AChE) in order to retrieve novel human AChE (hAChE) inhibitors. In addition, a protocol was developed to identify novel hit compounds and propose new promising AChE inhibitors from the ZINC database with 10 million commercially available compounds. In this sense, 3D similarity searches using rapid overlay of chemical structures and similarity analysis through comparison of electrostatic overlay of docked hits were used to retrieve AChE inhibitors from collected databases. Molecular dynamics simulation of 100 ns was carried out to study the best docked compounds from similarity searches. Some key residues were identified as crucial for the dual binding mode of inhibitor with the interaction site. All results indicated the relevant use of EON and docking strategy for identifying novel hit compounds as promising potential anticholinesterase candidates, and seven new structures were selected as potential hAChE inhibitors. Graphical abstract Compound N01 in the 4M0E hAChE crystallography structure from docking results. Yellow dashed lines Hydrogen bonds, blue dashed lines π-stacking interactions, green dashed lines cation-π interactions.

Anti-trypanosomal activity of non-peptidic nitrile-based cysteine protease inhibitors.

The cysteine protease cruzipain is considered to be a validated target for therapeutic intervention in the treatment of Chagas disease. Anti-trypanosomal activity against the CL Brener strain of T. cruzi was observed in the 0.1 µM to 1 µM range for three nitrile-based cysteine protease inhibitors based on two scaffolds known to be associated with cathepsin K inhibition. The two compounds showing the greatest potency against the trypanosome were characterized by EC50 values (0.12 µM and 0.25 µM) that were an order of magnitude lower than the corresponding Ki values measured against cruzain, a recombinant form of cruzipain, in an enzyme inhibition assay. This implies that the anti-trypanosomal activity of these two compounds may not be explained only by the inhibition of the cruzain enzyme, thereby triggering a putative polypharmacological profile towards cysteine proteases.

Molecular Design, Synthesis and Trypanocidal Activity of Dipeptidyl Nitriles as Cruzain Inhibitors.

A series of compounds based on the dipeptidyl nitrile scaffold were synthesized and assayed for their inhibitory activity against the T. cruzi cysteine protease cruzain. Structure activity relationships (SARs) were established using three, eleven and twelve variations respectively at the P1, P2 and P3 positions. A Ki value of 16 nM was observed for the most potent of these inhibitors which reflects a degree of non-additivity in the SAR. An X-ray crystal structure was determined for the ligand-protein complex for the structural prototype for the series. Twenty three inhibitors were also evaluated for their anti-trypanosomal effects and an EC50 value of 28 µM was observed for the most potent of these. Although there remains scope for further optimization, the knowledge gained from this study is also transferable to the design of cruzain inhibitors based on warheads other than nitrile as well as alternative scaffolds.

Non-peptidic cruzain inhibitors with trypanocidal activity discovered by virtual screening and in vitro assay.

A multi-step cascade strategy using integrated ligand- and target-based virtual screening methods was developed to select a small number of compounds from the ZINC database to be evaluated for trypanocidal activity. Winnowing the database to 23 selected compounds, 12 non-covalent binding cruzain inhibitors with affinity values (K i) in the low micromolar range (3-60 µM) acting through a competitive inhibition mechanism were identified. This mechanism has been confirmed by determining the binding mode of the cruzain inhibitor Nequimed176 through X-ray crystallographic studies. Cruzain, a validated therapeutic target for new chemotherapy for Chagas disease, also shares high similarity with the mammalian homolog cathepsin L. Because increased activity of cathepsin L is related to invasive properties and has been linked to metastatic cancer cells, cruzain inhibitors from the same library were assayed against it. Affinity values were in a similar range (4-80 µM), yielding poor selectivity towards cruzain but raising the possibility of investigating such inhibitors for their effect on cell proliferation. In order to select the most promising enzyme inhibitors retaining trypanocidal activity for structure-activity relationship (SAR) studies, the most potent cruzain inhibitors were assayed against T. cruzi-infected cells. Two compounds were found to have trypanocidal activity. Using compound Nequimed42 as precursor, an SAR was established in which the 2-acetamidothiophene-3-carboxamide group was identified as essential for enzyme and parasite inhibition activities. The IC50 value for compound Nequimed42 acting against the trypomastigote form of the Tulahuen lacZ strain was found to be 10.6±0.1 µM, tenfold lower than that obtained for benznidazole, which was taken as positive control. In addition, by employing the strategy of molecular simplification, a smaller compound derived from Nequimed42 with a ligand efficiency (LE) of 0.33 kcal mol(-1) atom(-1) (compound Nequimed176) is highlighted as a novel non-peptidic, non-covalent cruzain inhibitor as a trypanocidal agent candidate for optimization.

A combined study using ligand-based design, synthesis, and pharmacological evaluation of analogues of the acetaminophen ortho-regioisomer with potent analgesic activity.

A ligand-based drug design study was performed to acetaminophen regioisomers as analgesic candidates employing quantum chemical calculations at the DFT/B3LYP level of theory and the 6-31G* basis set. To do so, many molecular descriptors were used such as highest occupied molecular orbital, ionization potential, H-O bond dissociation energies, and spin densities, which might be related to quench reactivity of the tyrosyl radical to give N-acetyl-p-benzosemiquinone-imine through an initial electron withdrawing or hydrogen atom abstraction. Based on this in silico work, the most promising molecule, orthobenzamol, was synthesized and tested. The results expected from the theoretical prediction were confirmed in vivo using mouse models of nociception such as writhing, paw licking, and hot plate tests. All biological results suggested an antinociceptive activity mediated by opioid receptors. Furthermore, at 90 and 120 min, this new compound had an effect that was comparable to morphine, the standard drug for this test. Finally, the pharmacophore model is discussed according to the electronic properties derived from quantum chemistry calculations.

The GRID/CPCA approach in drug discovery.

IMPORTANCE OF THE FIELD: Bimolecular recognition is the basis for almost all processes in biological systems. The geometrical and chemical complementarity of small molecule ligands and their macromolecular biological targets, matching paired interacting parts, can result in binding that will eventually yield a biological response. AREAS COVERED IN THIS REVIEW: The topics covered include the integration of molecular interaction fields and chemometrics, via the GRID/CPCA (consensus principal component analysis) method that is actively contributing to the optimization of potency and selectivity of ligands towards a chosen target. Key applications that hallmark the usefulness of the method are critically presented. By comparison of the GRID/CPCA and GRID/PCA, the breakthroughs and challenges are highlighted. WHAT THE READER WILL GAIN: Molecular recognition studies support the development of pharmacophore-based descriptors, which provide the means to identify new ligand templates ('scaffold-hopping'). The GRID/CPCA approach can simultaneously reveal common trends in more than one block of data for more than two target proteins, with several three-dimensional structures per protein. It offers the benefit of improving the weighting between different interaction energy probes within the GRID parameterization. An important consequence is that hydrophobic interactions can be assessed for selectivity between proteins. TAKE HOME MESSAGE: Molecular-field-based methods along with CPCA analysis are extremely powerful to understand bimolecular interactions. Because drug discovery and development is a costly, time consuming and high-risk activity and GRID/CPCA is at the forefront of the computer-aided design, it should be used as early as possible for discovering new drugs.

Novel insights for dihydroorotate dehydrogenase class 1A inhibitors discovery.

The enzyme dihydroorotate dehydrogenase (DHODH) has been suggested as a promising target for the design of trypanocidal agents. We report here the discovery of novel inhibitors of Trypanosoma cruzi DHODH identified by a combination of virtual screening and ITC methods. Monitoring of the enzymatic reaction in the presence of selected ligands together with structural information obtained from X-ray crystallography analysis have allowed the identification and validation of a novel site of interaction (S2 site). This has provided important structural insights for the rational design of T. cruzi and Leishmania major DHODH inhibitors. The most potent compound (1) in the investigated series inhibits TcDHODH enzyme with Kiapp value of 19.28 µM and possesses a ligand efficiency of 0.54 kcal mol(-1) per non-H atom. The compounds described in this work are promising hits for further development.