A quantitative structure-activity relationship study on some HIV-1 protease inhibitors using molecular connectivity index

A quantitative structure-activity relationship (QSAR) study has been made on two different series of tetrahydropyrimidinones acting as HIV-1 protease inhibitors. A structural parameter, the first order valence molecular connectivity index ((1)chi(v)), has been used to account for the variation in the activity. The protease inhibition activity as well as the antiviral potency of the compounds are found to be significantly correlated with (1)chi(v) of P(2)/P(2') substituents attached to the two nitrogens N1 and N3, suggesting that substituents containing less electronegative and more saturated atoms, meaning thereby the less polar or more hydrophobic substituents, will be more advantageous. Further, if P(2) and P(2') are dissimilar, the former is found to be more effective than the latter. This difference is attributed to a conformational change in the enzyme that may be more favorable to P(2) binding than to P(2') binding.     

Trypsin activation. Effect of the Ile-Val dipeptide concentration on Kazal inhibitor binding to bovine trypsinogen

The effect of Ile-Val concentration (up to 2.0 M) on the thermodynamic parameters for the binding of the porcine pancreatic secretory trypsin inhibitor (Kazal inhibitor) to trypsinogen has been investigated at pH 5.5 between 7 degrees C and 42 degrees C. Thermodynamic parameters for Kazal inhibitor binding to the Ile-Val:zymogen adduct are more favorable than those observed for inhibitor association to the free proenzyme, but less so than those reported for beta-trypsin:Kazal inhibitor adduct formation (even under saturating dipeptide concentrations), suggesting that the effector dipeptide does not induce a complete rigidification of the proenzyme's activation domain. Considering the dependence of the association equilibrium constant for Kazal inhibitor binding to trypsinogen from Ile-Val concentration, thermodynamic parameters for the effector dipeptide binding to the free proenzyme and to its binary complex with Kazal inhibitor have been obtained. Differences in affinity for Ile-Val binding to the free zymogen and its binary complexes with inhibitors and substrates are indicative of the presence of different activation levels of the proenzyme, none of them exactly coincident with that of beta-trypsin. Such different discrete states should correspond to those involved in the zymogen-to-active-enzyme transition which should not be considered as an all-or-nothing process, but as a multistep event.     

A calorimetric study of the binding of S-alkylglutathiones to glutathione S-transferase.

The binding of three competitive glutathione analogue inhibitors (S-alkylglutathione derivatives) to glutathione S-transferase from Schistosoma japonicum, SjGST, has been investigated by isothermal titration microcalorimetry at pH 6.5 over a temperature range of 15--30 degrees C. Calorimetric measurements in various buffer systems with different ionization heats suggest that no protons are exchanged during the binding of S-alkylglutathione derivatives. Thus, at pH 6.5, the protons released during the binding of substrate may be from its thiol group. Calorimetric analyses show that S-methyl-, S-butyl-, and S-octylglutathione bind to two equal and independent sites in the dimer of SjGST. The affinity of these inhibitors to SjGST is greater as the number of methylene groups in the hydrocarbon side chain increases. In all cases studied, Delta G(0) remains invariant as a function of temperature, while Delta H(b) and Delta S(0) both decrease as the temperature increases. The binding of three S-alkylglutathione derivatives to the enzyme is enthalpically favourable at all temperatures studied. The temperature dependence of the enthalpy change yields negative heat capacity changes, which become less negative as the length of the side chain increases.     

Urethanyl-3-amidinophenylalanine derivatives as inhibitors of factor Xa. X-ray crystal structure of a trypsin/inhibitor complex and modeling studies

Hydrophobic urethanyl derivatives of 3-amidinophenylalanine methyl ester were found to be relatively potent and selective factor Xa inhibitors. These compounds consist of the arginine-mimetic 3-benzamidino group as P1 residue and of hydrophobic residues as potential interaction partners for the S3/S4 aryl binding site of the enzyme. Attempts to possibly identify their binding mode to factor Xa via the X-ray crystal structure of a trypsin/inhibitor complex and analogy modeling on the crystal structure of factor Xa failed. However, synthesis of enantiomerically pure (R)- and (S)-derivatives, combined with modeling experiments, led to an hypothetical non-substrate like binding mode, which was fully confirmed by the remarkably enhanced inhibitory potency of derivatives in which the methyl ester was replaced by arylamides for interactions with the S3/S4 enzyme binding subsites. With adamantyloxycarbonyl-(R)-3-amidinophenylalanine-phenethylamide+ ++ a nanomolar inhibiton was obtained, thus indicating this new class of factor Xa inhibitors as a highly promising lead structure.     

Thermodynamic determination of the binding constants of angiotensin-converting enzyme inhibitors by a displacement method

Somatic angiotensin I-converting enzyme (s-ACE) plays a central role in blood pressure regulation and has been the target of most antihypertensive drugs. A displacement isothermal titration calorimetry method has been used to accurately determine the binding constant of three strong s-ACE inhibitors. Under the experimental conditions studied in this work, the relative potency of the inhibitors was determined to be enalaprilat>lisinopril>captopril. We analyze the thermodynamic behaviour of the binding process using the new structural information provided by the ACE structures, as well as the conformational changes that occur upon binding.     

Quantitative assessment of the noncovalent inhibition of sickle hemoglobin gelation by phenyl derivatives and other known agents.

The ability of a variety of phenyl derivatives to inhibit sickle cell hemoglobin gelation was placed on a quantitative scale by parallel equilibrium and kinetic assays. Modifications of the phenyl ring studied include polar, nonpolar, and charged substituents, added aromatic rings, and loss of aromaticity. Other noncovalent inhibitors previously reported to have high potency were measured and placed on the same quantitative scale. Some phenyl derivatives were found to be as effective an any other known noncovalent antigelling agent. The phenyl compounds penetrate easily into red cells, and their potency is tolerant to chemical modification, which holds out the possibility of designing low-toxicity derivatives. On the negative side, the level of potency obtainable appears to be inadequate for clinical use. The best phenyl inhibitors display a functionally defined inhibitory constant (K1) of 75 mM, and it can be estimated that inhibitor concentrations over 20 mM would be necessary to obtain minimal clinically significant benefit. Furthermore, with the variety of modifications tested here, no impressive increase in activity could be achieved over that found in the simplest phenyl compounds.     

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.     

Trypsin from Greenland cod, Gadus ogac. Isolation and comparative properties

Trypsin(ogen) was isolated from the pyloric ceca of Greenland cod. Greenland cod trypsin catalyzed hydrolysis of N alpha-benzoyl-DL-arginine p-nitroanilide, tosyl arginine methyl ester and protein and was inhibited by the serine protease inhibitor PMSF and other well-known trypsin inhibitors. Greenland cod trypsin was more stable at alkaline pH than at acid pH; and was inactivated by relatively low thermal treatment. Like other trypsins, the enzyme was rich in potential acidic amino acid residues but poor in basic amino acid residues and had a molecular weight of 23,500; but it had less potential disulfide pairs, less alpha-helix and a lower H phi ave than other trypsins previously characterized. Reactions catalyzed by Greenland cod trypsin were not very responsive to temperature change, such that specific activity was relatively high at low reaction temperature.     

A novel mode of protein kinase inhibition exploiting hydrophobic motifs of autoinhibited kinases: discovery of ATP-independent inhibitors of fibroblast growth factor receptor

Protein kinase inhibitors with enhanced selectivity can be designed by optimizing binding interactions with less conserved inactive conformations because such inhibitors will be less likely to compete with ATP for binding and therefore may be less impacted by high intracellular concentrations of ATP. Analysis of the ATP-binding cleft in a number of inactive protein kinases, particularly in the autoinhibited conformation, led to the identification of a previously undisclosed non-polar region in this cleft. This ATP-incompatible hydrophobic region is distinct from the previously characterized hydrophobic allosteric back pocket, as well as the main pocket. Generalized hypothetical models of inactive kinases were constructed and, for the work described here, we selected the fibroblast growth factor receptor (FGFR) tyrosine kinase family as a case study. Initial optimization of a FGFR2 inhibitor identified from a library of commercial compounds was guided using structural information from the model. We describe the inhibitory characteristics of this compound in biophysical, biochemical, and cell-based assays, and have characterized the binding mode using x-ray crystallographic studies. The results demonstrate, as expected, that these inhibitors prevent activation of the autoinhibited conformation, retain full inhibitory potency in the presence of physiological concentrations of ATP, and have favorable inhibitory activity in cancer cells. Given the widespread regulation of kinases by autoinhibitory mechanisms, the approach described herein provides a new paradigm for the discovery of inhibitors by targeting inactive conformations of protein kinases.     

Thermodynamic analysis of the binding of galactose and poly-N-acetyllactosamine derivatives to human galectin-3.

Galectin-3, with a wide tissue distribution and marked developmental regulation, provides significant insights into the progression of various disease and developmental stages. Recognized by its specificity for galactose, a detailed characterization of its sugar binding ability has been investigated by isothermal titration calorimetry. The results presented here complement well with the earlier studies utilizing hapten inhibition assays. Among the various lactose derivatives studied, A-tetrasaccharide emerged with the highest affinity for binding to galectin-3 combining site. This blood group saccharide exhibited a binding affinity 37-fold higher and a 102 kJ/mol more favorable change in enthalpy over lactose at 280 K indicating the existence of additional subsites for both the alpha1-3-linked N-acetylgalactosamine at the non-reducing end and the alpha1-2-linked L-fucosyl residue. The thermodynamic parameters evaluated for other ligands substantiate further the carbohydrate recognition domain to be part of an extended binding site. Binding thermodynamics of galectin-3 with the galactose derivatives are essentially enthalpically driven and exhibit compensatory changes in DeltaH degrees and TDeltaS owing to solvent reorganization.     

Protease inhibitors: Part 4. Synthesis of weakly basic thrombin inhibitors incorporating pyridinium-sulfanilylaminoguanidine moieties

Three series of derivatives have been prepared by reaction of sulfanilylaminoguanidine with pyrylium salts, with the pyridinium derivatives of glycine and with the pyridinium derivatives of beta-alanine, respectively. The new compounds were assayed as inhibitors of two serine proteases, thrombin and trypsin. The study showed that in contrast to the leads, possessing KI's around 100-300 nM against thrombin, and 450-1420 nM against trypsin, respectively, the new derivatives showed inhibition constants in the range of 15-50 nM against thrombin, whereas their affinity for trypsin remained relatively low. Derivatives of beta-alanine were more active than the corresponding glycine derivatives, which in turn were more inhibitory than the pyridinium derivatives of sulfanilylaminoguanidine possessing the same substitution pattern at the pyridinium ring. Thus, the present study proposes two novel approaches for the preparation of high affinity, specific thrombin inhibitors: a novel S1 anchoring moiety in the already large family of arginine/amidine-based inhibitors, i.e., the SO2NHNHC(=NH)NH2 group, and novel non-peptidomimetic scaffolds obtained by incorporating alkyl-/aryl-substituted-pyridinium moieties in the hydrophobic binding site(s). The first one is important for obtaining bioavailable thrombin inhibitors, devoid of the high basicity of the commonly used arginine/amidine-based inhibitors, whereas the second one may lead to improved water solubility of such compounds.     

Epithelial chloride channel. Development of inhibitory ligands

Chloride channels are present in the majority of epithelial cells, where they mediate absorption or secretion of NaCl. Although the absorptive and secretory channels are well characterized in terms of their electrophysiological behavior, there is a lack of pharmacological ligands that can aid us in further functional and eventually molecular characterization. To obtain such ligands, we prepared membrane vesicles from bovine kidney cortex and apical membrane vesicles from trachea and found that they contain a chloride transport process that is electrically conductive. This conductance was reduced by preincubating the vesicles in media containing ATP or ATP-gamma-S, but not beta-methylene ATP, which suggests that the membranes contain a kinase that can close the channels. We then screened compounds derived from three classes: indanyloxyacetic acid (IAA), anthranilic acid (AA), and ethacrynic acid. We identified potent inhibitors from the IAA and the AA series. We tritiated IAA-94 and measured binding of this ligand to the kidney cortex membrane vesicles and found a high-affinity binding site whose dissociation constant (0.6 microM) was similar to the inhibition constant (1 microM). There was a good correlation between the inhibitory potency of several IAA derivatives and their efficacy in displacing [3H]IAA-94 from its binding site. Further, other chloride channel inhibitors, including AA derivatives, ethacrynic acid, bumetanide, and DIDS, also displaced the ligand from its binding site. A similar conductance was found in apical membrane vesicles from bovine trachea that was also inhibited by IAA-94 and AA-130B, but the inhibitory effects of these compounds were weaker than their effects on the renal cortex channel. The two drugs were also less potent in displacing [3H]IAA-94 from the tracheal binding site.     

Naphthylisopropylamine and N-benzylamphetamine derivatives as monoamine oxidase inhibitors

A series of naphthylisopropylamine and N-benzyl-4-methylthioamphetamine derivatives were evaluated as monoamine oxidase inhibitors. Their potencies were compared with those of a series of amphetamine derivatives, to test if the increase of electron richness of the aromatic ring and overall size of the molecule might improve their potency as enzyme inhibitors. Molecular dockings were performed to gain insight regarding the binding mode of these inhibitors and rationalize their different potencies. In the case of naphthylisopropylamine derivatives, the increased electron-donating capacity and size of the aromatic moiety resulting from replacement of the phenyl ring of amphetamine derivatives by a naphthalene system resulted in more potent compounds. In the other case, extension of the arylisopropylamine molecule by N-benzylation of the amino group led to a decrease in potency as monoamine oxidase inhibitors.     

Protease Inhibitors: Part 4. Synthesis of Weakly Basic Thrombin Inhibitors Incorporating Pyridinium-Sulfanilylaminoguanidine Moieties

Three series of derivatives have been prepared by reaction of sulfanilylaminoguanidine with pyrylium salts, with the pyridinium derivatives of glycine and with the pyridinium derivatives of β-alanine, respectively. The new compounds were assayed as inhibitors of two serine proteases, thrombin and trypsin. The study showed that in contrast to the leads, possessing K₁'s around 100–300nM against thrombin, and 450–1420nM against trypsin, respectively, the new derivdtives showed inhibition constants in the range of 15–50nM against thrombin, whereas their affinity for trypsin remained relatively low. Derivatives of β-alanine were more active than the corresponding glycine derivatives, which in turn were more inhibitory than the pyridinium derivatives of sulfanilylaminoguanidine possessing the same substitution pattern at the pyridinium ring. Thus, the present study proposes two novel approaches for the preparation of high affinity, specific thrombin inhibitors: a novel S1 anchoring moiety in the already large family of arginine/amidine-based inhibitors, i.e., the SO₂NHNHC(=NH)NH₂ group, and novel non-peptidomimetic scaffolds obtained by incorporating alkyl-/aryl-substituted-pyridinium moieties in the hydrophobic binding site(s). The first one is important for obtaining bioavailable thrombin inhibitors, devoid of the high basicity of the commonly used arginine/amidine-based inhibitors, whereas the second one may lead to improved water solubility of such compounds.     

Calcitonin gene-related peptide (CGRP) receptor antagonists: pyridine as a replacement for a core amide group

In our continuing efforts to identify CGRP receptor antagonists that can be dosed orally for the treatment of migraine headache, we have investigated a pyridine bioisosteric replacement of a polar amide portion of a previous lead compound, BMS-694153. Pyridine derivatives were discovered and their SAR was studied. Some of them showed excellent binding potency. However, oral bioavailability was low, even for compounds with good Caco-2 cell permeability.     

Think Twice: Understanding the High Potency of Bis(phenyl)methane Inhibitors of Thrombin

Successful design of potent and selective protein inhibitors, in terms of structure-based drug design, strongly relies on the correct understanding of the molecular features determining the ligand binding to the target protein. We present a case study of serine protease inhibitors with a bis(phenyl)methane moiety binding into the S3 pocket. These inhibitors bind with remarkable potency to the active site of thrombin, the blood coagulation factor IIa. A combination of X-ray crystallography and isothermal titration calorimetry provides conclusive insights into the driving forces responsible for the surprisingly high potency of these inhibitors. Analysis of six well-resolved crystal structures (resolution 1.58-2.25 Å) along with the thermodynamic data allows an explanation of the tight binding of the bis(phenyl)methane inhibitors. Interestingly, the two phenyl rings contribute to binding affinity for very different reasons — a fact that can only be elucidated by a structure-based approach. The first phenyl moiety occupies the hydrophobic S3 pocket, resulting in a mainly entropic advantage of binding. This observation is based on the displacement of structural water molecules from the S3 pocket that are observed in complexes with inhibitors that do not bind in the S3 pocket. The same classic hydrophobic effect cannot explain the enhanced binding affinity resulting from the attachment of the second, more solvent-exposed phenyl ring. For the bis(phenyl)methane inhibitors, an observed adaptive rotation of a glutamate residue adjacent to the S3 binding pocket attracted our attention. The rotation of this glutamate into salt-bridging distance with a lysine moiety correlates with an enhanced enthalpic contribution to binding for these highly potent thrombin binders. This explanation for the magnitude of the attractive force is confirmed by data retrieved by a Relibase search of several thrombin-inhibitor complexes deposited in the Protein Data Bank exhibiting similar molecular features.Special attention was attributed to putative changes in the protonation states of the interaction partners. For this purpose, two analogous inhibitors differing mainly in their potential to change the protonation state of a hydrogen-bond donor functionality were compared. Buffer dependencies of the binding enthalpy associated with complex formation could be traced by isothermal titration calorimetry, which revealed, along with analysis of the crystal structures (resolution 1.60 and 1.75 Å), that a virtually compensating proton interchange between enzyme, inhibitor and buffer is responsible for the observed buffer-independent thermodynamic signatures.     

Limited proteolysis reveals conformational changes in uncoupling protein-1 from brown adipose tissue mitochondria

Limited proteolytic digestion of uncoupling protein-1 (UCP1) from hamster brown adipose tissue mitochondria was studied. Under optimal conditions, trypsin and chymotrypsin cleave at Lys-292 and at Phe-102, yielding major products 31-kDa T1 and 22-kDa Ch1. Both T1 and Ch1 remained dimers, as in UCP1. Using fluorescent nucleotide derivative 2'-O-dansyl GTP, it is shown that T1 retains the nucleotide binding affinity (K(D)=1 microM for dansyl GTP) while Ch1 does not bind nucleotide. Previously kinetic binding and H(+) transport studies [Biochemistry 35 (1996) 7846] have shown that UCP1 forms tight complexes to varying degrees with nucleotides and their derivatives. Nucleotides strongly protect against tryptic digestion but less against chymotryptic digestion, because the chymotryptic product Ch1 does not bind nucleotide. The nucleotides and derivatives show the same potency profile in protecting against both trypsinolysis and chymotryptic digestion, suggesting that UCP1 undergoes a major conformational change upon nucleotide binding from an initial loose complex into a tight complex, in which the cleavage sites become masked from proteolysis.