Dithiol proteins as guardians of the intracellular redox milieu in parasites: old and new drug targets in trypanosomes and malaria-causing plasmodia

Parasitic diseases such as sleeping sickness, Chagas' heart disease, and malaria are major health problems in poverty-stricken areas. Antiparasitic drugs that are not only active but also affordable and readily available are urgently required. One approach to finding new drugs and rediscovering old ones is based on enzyme inhibitors that paralyze antioxidant systems in the pathogens. These antioxidant ensembles are essential to the parasites as they are attacked in the human host by strong oxidants such as peroxynitrite, hypochlorite, and H2O2. The pathogen-protecting system consists of some 20 thiol and dithiol proteins, which buffer the intraparasitic redox milieu at a potential of -250 mV. In trypanosomes and leishmania the network is centered around the unique dithiol trypanothione (N1,N8-bis(glutathionyl)spermidine). In contrast, malaria parasites have a more conservative dual antioxidative system based on glutathione and thioredoxin. Inhibitors of antioxidant enzymes such as trypanothione reductase are, indeed, parasiticidal but they can also delay or prevent resistance against a number of other antiparasitic drugs.     

Design criteria for molecular mimics of fragments of the β-turn. 1. Cα atom analysis

Peptides represent an extensive class of biologically active molecules. They may be used as leads in the development of novel therapeutic agents provided the pharmacophoric information present within them can be translated into non-peptide analogs that lack the peptide backbone and are stable to proteolysis. This is the rationale for peptidomimetic drug design. Frequently, the -turn has been implicated as a conformation important for biological recognition of peptides. Empirical evidence from known peptidomimetics, coupled with a theoretical model of peptide binding and the observation that glycine and proline residues are common within the -turn, has suggested the design of molecules to mimic placement of between two and four of the side-chains. The moderate number of different -turn conformations, combined with the combinatoric nature of side-chain selection complicates the procedure. In this paper, cluster analysis has been used to classify the arrangement of C_{_} atoms about the various fragments of the -turn. Recombination of the observed patterns provides a general model for the -turn which may be used as an effective screen for potential peptidomimetic scaffolds in chemical databases.     

Factor Xa: Simulation studies with an eye to inhibitor design

Factor Xa is a serine protease which activates thrombin and plays a key regulatory role in the blood-coagulation cascade. Factor Xa is at the crossroads of the extrinsic and intrinsic pathways of coagulation and, hence, has become an important target for the design of anti-thrombotics (inhibitors). It is not known to be involved in other processes than hemostasis and its binding site is different to that of other serine proteases, thus facilitating selective inhibition. The design of high-affinity selective inhibitors of factor Xa requires knowledge of the structural and dynamical characteristics of its active site. The three-dimensional structure of factor Xa was resolved by X-ray crystallography and refined at 2.2 Å resolution by Padmanabhan and collaborators. In this article we present results from molecular dynamics simulations of the catalytic domain of factor Xa in aqueous solution. The simulations were performed to characterise the mobility and flexibility of the residues delimiting the unoccupied binding site of the enzyme, and to determine hydrogen bonding propensities (with protein and with solvent atoms) of those residues in the active site that could interact with a substrate or a potential inhibitor. The simulation data is aimed at facilitating the design of high-affinity selective inhibitors of factor Xa.     

Efficiency of connected binary block designs when a single observation is unavailable

In this paper the problem of finding the design efficiency is considered when a single observation is unavailable in a connected binary block design. The explicit expression of efficiency is found for the resulting design when the original design is a balanced incomplete block design or a group divisible, singular or semiregular or regular with ₁>0, design. The efficiency does not depend on the position of the unavailable observation. For a regular group divisible design with ₁>0, the efficiency depends on the position of the unavailable observation. The bounds, both lower and upper, on the efficiency are given in this situation. The efficiencies of designs resulting from a balanced incomplete block design and a group divisible design are in fact high when a single observation is unavailable.The work of the first author is sponsored by the Air Force Office of Scientific Research under Grant AFOSR-90-0092.On leave from Indian Statistical Institute, Calcutta, India. The work of the third author was supported by a grant from the CMDS, Indian Institute of Management, Calcutta.     

Modified AutoDock for accurate docking of protein kinase inhibitors

Protein kinases are an important class of enzymes controlling virtually all cellular signaling pathways. Consequently, selective inhibitors of protein kinases have attracted significant interest as potential new drugs for many diseases. Computational methods, including molecular docking, have increasingly been used in the inhibitor design process [1]. We have considered several docking packages in order to strengthen our kinase inhibitor work with computational capabilities. In our experience, AutoDock offered a reasonable combination of accuracy and speed, as opposed to methods that specialize either in fast database searches or detailed and computationally intensive calculations.However, AutoDock did not perform well in cases where extensive hydrophobic contacts were involved, such as docking of SB203580 to its target protein kinase p38. Another shortcoming was a hydrogen bonding energy function, which underestimated the attraction component and, thus, did not allow for sufficiently accurate modeling of the key hydrogen bonds in the kinase-inhibitor complexes.We have modified the parameter set used to model hydrogen bonds, which increased the accuracy of AutoDock and appeared to be generally applicable to many kinase-inhibitor pairs without customization. Binding to largely hydrophobic sites, such as the active site of p38, was significantly improved by introducing a correction factor selectively affecting only carbon and hydrogen energy grids, thus, providing an effective, although approximate, treatment of solvation.     

Modelling of GC and HPLC separations of simazin and atrazin by experimental design methodology

Summary Experimental Design methodology allows the modelling and optimization of the chromatographic separation of similar pesticides (triazine family) by GC and HPLC. The GC separation of simazin and atrazin is well modelled by a first degree equation, involving injected volume, carrier gas pressure and rising oven temperature. The LC is modelled by a second degree equation, depending on injected volume, eluent flow and composition. These calculated models allow easy optimization of the separations, using isoresponse curves.     

Computational combinatorial ligand design: Application to human α-thrombin

Summary A new method is presented for computer-aided ligand design by combinatorial selection of fragments that bind favorably to a macromolecular target of known three-dimensional structure. Firstly, the multiple-copy simultaneous-search procedure (MCSS) is used to exhaustively search for optimal positions and orientations of functional groups on the surface of the macromolecule (enzyme or receptor fragment). The MCSS minima are then sorted according to an approximated binding free energy, whose solvation component is expressed as a sum of separate electrostatic and nonpolar contributions. The electrostatic solvation energy is calculated by the numerical solution of the linearized Poisson-Boltzmann equation, while the nonpolar contribution to the binding free energy is assumed to be proportional to the loss in solvent-accessible surface area. The program developed for computational combinatorial ligand design (CCLD) allows the fast and automatic generation of a multitude of highly diverse compounds, by connecting in a combinatorial fashion the functional groups in their minimized positions. The fragments are linked as two atoms may be either fused, or connected by a covalent bond or a small linker unit. To avoid the combinatorial explosion problem, pruning of the growing ligand is performed according to the average value of the approximated binding free energy of its fragments. The method is illustrated here by constructing candidate ligands for the active site of human -thrombin. The MCSS minima with favorable binding free energy reproduce the interaction patterns of known inhibitors. Starting from these fragments, CCLD generates a set of compounds that are closely related to high-affinity thrombin inhibitors. In addition, putative ligands with novel binding motifs are suggested. Probable implications of the MCSS-CCLD approach for the evolving scenario of drug discovery are discussed.     

Functionality map analysis of the active site cleft of human thrombin

Summary The Multiple Copy Simultaneous Search methodology has been used to construct functionality maps for an extended region of human thrombin, including the active site. This method allows the determination of energetically favorable positions and orientations for functional groups defined by the user on the three-dimensional surface of a protein. The positions of 10 functional group sites are compared with those of corresponding groups of four thrombin-inhibitor complexes. Many, but not all features, of known thrombin inhibitors are reproduced by the method. The results indicate that certain aspects of the binding modes of these inhibitors are not optimal. In addition, suggestions are made for improving binding by interaction with functional group sites on the thrombin surface that are not used by the thrombin inhibitors. Abbreviations: MCSS, multiple copy simultaneous search; PPACK, d-phenylalanyl-l-propyl-l-arginine chloromethane; NAPAP, N -(2-naphthylsulfonylglycyl)-d-para-amidinophenylalanylpiperidine; argatroban, (2R,4R)-4-methyl-1-[N -(3-methyl-1,2,3,4-tetrahydro-8-quinolinylsulfonyl)-l-arginyl]-2-piperidine carboxylic acid; rms, root mean square. The thrombin residues are numbered according to the chymotrypsin-based numbering by Bode et al. [8]. P1, P2, P3, etc., denote the peptide inhibitor residues on the amino-terminal side of the scissile peptide bond, and S1, S2, S3, etc., the corresponding subsites of thrombin     

Thermodynamic Parameters for the Association of Fluorinated Benzenesulfonamides with Bovine Carbonic Anhydrase II

This paper describes a calorimetric study of the association of a series of seven fluorinated benzenesulfonamide ligands (C₆H_nF_5?nSO₂NH₂) with bovine carbonic anhydrase II (BCA). Quantitative structure–activity relationships between the free energy, enthalpy, and entropy of binding and pK_a and log P of the ligands allowed the evaluation of the thermodynamic parameters in terms of the two independent effects of fluorination on the ligand: its electrostatic potential and its hydrophobicity. The parameters were partitioned to the three different structural interactions between the ligand and BCA: the Zn^II cofactor–sulfonamide bond (≈65 % of the free energy of binding), the hydrogen bonds between the ligand and BCA (≈10 %), and the contacts between the phenyl ring of the ligand and BCA (≈25 %). Calorimetry revealed that all of the ligands studied bind in a 1:1 stoichiometry with BCA; this result was confirmed by ¹⁹F NMR spectroscopy and X‐ray crystallography (for complexes with human carbonic anhydrase II).     

Trifluoromethylation of Alkyl Radicals: Breakthrough and Challenges†

Trifluoromethylation of alkyl radicals is emerging as a powerful tool for C(sp³)–CF₃ bond formations. Based on the hypothesis of CF₃ group transfer from Cu(II)–CF₃ to alkyl radicals, a number of trifluoromethylation reactions have been developed, including trifluoromethylation of alkyl halides, decarboxylative trifluoromethylation of aliphatic carboxylic acids, C(sp³)–H trifluoromethylation, amino‐ and carbo‐trifluoromethylation of alkenes, etc. Challenges in this intriguing field are also discussed.