Probing of CD4 binding pocket of HIV-1 gp120 glycoprotein using unnatural phenylalanine analogues

CD4-gp120 interaction is the first step for HIV-1 entry into host cells. A highly conserved pocket in gp120 protein is an attractive target for developing gp120 inhibitors or novel HIV detection tools. Here we incorporate seven phenylalanine derivatives having different sizes and steric conformations into position 43 of domain 1 of CD4 (mD1.2) to explore the architecture of the ‘Phe43 cavity’ of HIV-1 gp120. The results show that the conserved hydrophobic pocket in gp120 tolerates a hydrophobic side chain of residue 43 of CD protein, which is 12.2 Å in length and 8.0 Å in width. This result provides useful information for developing novel gp120 inhibitors or new HIV detection tools.     

Probing binding requirements of NAD kinase with modified substrate (NAD) analogues

Synthesis of novel NAD(+) analogues that cannot be phosphorylated by NAD kinase is reported. In these analogues the C2' hydroxyl group of the adenosine moiety was replaced by fluorine in the ribo or arabino configuration (1 and 2, respectively) or was inverted into arabino configuration to give compound 3. Compounds 1 and 2 showed inhibition of human NAD kinase, whereas analogue 3 inhibited both the human and Mycobacterium tuberculosis NAD kinase. An uncharged benzamide adenine dinucleotide (BAD) was found to be the most potent competitive inhibitor (K(i)=90 microM) of the human enzyme reported so far.     

Acylguanidine inhibitors of beta-secretase: optimization of the pyrrole ring substituents extending into the S1' substrate binding pocket

The proteolytic enzyme beta-secretase (BACE-1) produces amyloid beta (Abeta) peptide, the primary constituent of neurofibrillary plaques, implicated in Alzheimer's disease, by cleavage of the amyloid precursor protein. A small molecule inhibitor of BACE-1, (diaminomethylene)-2,5-diphenyl-1H-pyrrole-1-acetamide (1, BACE-1 IC(50)=3.7 microM), was recently described, representing a new small molecule lead. Initial SAR investigation demonstrated the potential of accessing the nearby S(3) and S(1)(') substrate binding pockets of the BACE-1 enzyme by building substituents off one of the phenyl substituents and guanidinyl functional group. We report here the optimization of guanidinyl functional group substituents on 1, leading to potent submicromolar BACE-1 inhibitors.     

Computational design of new cyclic urea inhibitors for improved binding of HIV-1 aspartic protease

We report in this paper the design, by means of computational techniques, of new cyclic urea inhibitors of the HIV aspartic protease. The relationship between the complexation energies of the enzyme with known inhibitors and the experimentally determined log K(i) have been studied and used to predict inhibition constants for new inhibitors.     

HIV-1 protease substrate binding and product release pathways explored with coarse-grained molecular dynamics

We analyze the encounter of a peptide substrate with the native HIV-1 protease, the mechanism of substrate incorporation in the binding cleft, and the dissociation of products after substrate hydrolysis. To account for the substrate, we extend a coarse-grained model force field, which we previously developed to study the flap opening dynamics of HIV-1 protease on a microsecond timescale. Molecular and Langevin dynamics simulations show that the flaps need to open for the peptide to bind and that the protease interaction with the substrate influences the flap opening frequency and interval. On the other hand, release of the products does not require flap opening because they can slide out from the binding cleft to the sides of the enzyme. Our data show that in the protease-substrate complex the highest fluctuations correspond to the 17- and 39-turns and the substrate motion is anticorrelated with the 39-turn. Moreover, the active site residues and the flap tips move in phase with the peptide. We suggest some mechanistic principles for how the flexibility of the protein may be involved in ligand binding and release.     

Closing of the flaps of HIV-1 protease induced by substrate binding: a model of a flap closing mechanism in retroviral aspartic proteases

The active site of aspartic proteases is covered by one or more flaps, which control access to the active site and play a significant role in the binding of the substrate. An extensive conformational change of the flaps takes place upon binding of substrate to the active site. A long molecular dynamics simulation was performed on the complex consisting of a peptide (CA-p2) from a natural substrate cleavage site of the gag/pol polyprotein placed in the active site of HIV-1 protease (PR) with an open flap conformation. During the simulation, the substrate induced the closing of the flaps into the closed conformation in an asymmetrical way through a hydrophobic intermediate state cluster. The nature of the residues of HIV-1 PR identified to be important in the flap closing mechanism is conserved across known structures of retroviral aspartic proteases family. The flap closing mechanism described in HIV-1 PR is proposed to be a general model for flap closing in retroviral aspartic proteases.     

Selective alteration of substrate specificity by replacement of aspartic acid-189 with lysine in the binding pocket of trypsin

To test the role of Asp-189 which is located at the base of the substrate binding pocket in determining the specificity of trypsin toward basic substrates, this residue was replaced with a lysine residue by site-directed mutagenesis. Both rat trypsinogen and Lys-189 trypsinogen were expressed and secreted into the periplasmic space of Escherichia coli. The proteins were purified to homogeneity and activated by porcine enterokinase, and their catalytic activities were determined on natural and synthetic substrates. Lys-189 trypsin displayed no catalytic activity toward arginyl and lysyl substrates. Further, there was no compensatory change in specificity toward acidic substrates; no cleavage of aspartyl or glutamyl bonds was detected. Additional studies of substrate specificity involving gas-phase sequence analyses of digested natural substrates revealed an inherent but low chymotrypsin-like activity of trypsin. This activity was retained but modified by the Asp to Lys change at position 189. In addition to hydrolyzing phenylalanyl and tyrosyl peptide bonds, the mutant enzyme has the unique property of cleaving leucyl bonds. On the basis of computer graphic modeling studies of the Lys-189 side chain, it appears that the positively charged NH2 group is directed outside the substrate binding pocket. The resulting hydrophobic cavity may explain the altered substrate specificity of the mutant enzyme. The relatively low chymotrypsin-like activity of both recombinant enzymes may be due to distorted positioning of the scissile bond with respect to the catalytic triad rather than to the lack of sufficient interaction between the hydrophobic side chains and the substrate binding pocket of the enzyme.     

Comparison of inhibitor binding in HIV-1 protease and in non-viral aspartic proteases: the role of the flap

The crystal structure of HIV-1 protease with an inhibitor has been compared with the structures of non-viral aspartic proteases complexed with inhibitors. In the dimeric HIV-1 protease, two 4-stranded beta-sheets are formed by half of the inhibitor, residues 27-29, and the flap from each monomer. In the monomeric non-viral enzyme the single flap does not form a beta-sheet with an inhibitor. The HIV-1 protease shows more interactions with a longer peptide inhibitor than are observed in non-viral aspartic protease-inhibitor complexes. This, and the large movement of the flaps, restricts the conformation of the protease cleavage sites in the retroviral polyprotein precursor.     

P1' oxadiazole protease inhibitors with excellent activity against native and protease inhibitor-resistant HIV-1

HIV-1 protease inhibitors (PI's) bearing 1,3,4-oxadiazoles at the P1' position were prepared by a novel method involving the diastereoselective installation of a carboxylic acid and conversion to the P1' heterocycle. The compounds are picomolar inhibitors of native HIV-1 protease, with most of the compounds maintaining excellent antiviral activity against a panel of PI-resistant strains.     

HIV-1 protease inhibitors with picomolar potency against PI-resistant HIV-1 by modification of the P1' substituent

Transposition of the pyridyl nitrogen from the P(3) substituent to the P(1)' substituent in HIV-1 protease inhibitors (PI) affords compounds such as 3 with an improved inhibitory profile against multiple P450 isoforms. These compounds also displayed increased potency, with 3 inhibiting viral spread (CIC(95)) at <8 nM for every strain of PI-resistant HIV-1 tested. The poor to modest bioavailability of these compounds may correlate in part to their aqueous solubility.     

Cyclic sulfamide HIV-1 protease inhibitors, with sidechains spanning from P2/P2' to P1/P1'

Previous studies of HIV protease inhibitors have shown that it is possible to elongate the P1/P1' sidechains to reach the S3/S3' binding sites. By analogy, we expected that it would be possible to design inhibitors reaching between the S1/S1' and S2/S2' binding sites. Molecular modeling suggested that this could be achieved with appropriate ortho-substitution of the P2/P2' benzyl groups in our cyclic sulfamide inhibitors. Four different spacer groups were investigated. The compounds were smoothly prepared from tartaric acid in five steps and exhibit low to moderate activity, the most potent inhibitor possessing a Ki value of 0.53 microM.     

Combinatorial design of nonsymmetrical cyclic urea inhibitors of aspartic protease of HIV-1

The aspartic protease (PR) of the human immunodeficiency virus type 1 (HIV-1) is an important target for the design of specific antiviral agents dedicated to treatment of HIV-1 infection. We have employed computer-assisted combinatorial chemistry methods to design a small focused virtual library of nonsymmetrically substituted cyclic urea inhibitors of the PR. Nonsymmetrical compounds with decreased peptidic character were namely found to inhibit the PR with comparable inhibition potencies as their C2-pseudosymmetric counterparts and to possess superior pharmacokinetic properties. To generate the virtual library of fully nonsymmetrical cyclic urea analogs, diverse reagents were selected from databases of available chemicals with characteristics similar to those of the building blocks of known potent PR inhibitors. The X-ray structure of the protease-inhibitor complex PR-XV-638 was used as the receptor model in the structure-based focusing and in silico screening of the virtual library. A target-specific LUDI-type scoring function, parameterized for a QSAR training set of known cyclic urea inhibitors and validated on a set of compounds not included into the training set, was used to predict the inhibition constants (Ki) of the generated analogs toward the HIV-1 PR. The fragments most frequently occurring in the analogs with the highest predicted inhibition potencies (Ki*<10 pM) were then selected to constitute a highly focused library subset containing novel nonsymmetrical cyclic ureas with predicted Ki*s 1 order of magnitude lower than the most potent known cyclic urea inhibitors. ADME properties calculated for the most promising analogs suggested that the cyclic ureas are endowed with a wide range of favorable pharmacokinetic properties, which may favor the discovery of a potent orally administrable antiviral drug.     

Strength of hydrogen bond network takes crucial roles in the dissociation process of inhibitors from the HIV-1 protease binding pocket

To understand the underlying mechanisms of significant differences in dissociation rate constant among different inhibitors for HIV-1 protease, we performed steered molecular dynamics (SMD) simulations to analyze the entire dissociation processes of inhibitors from the binding pocket of protease at atomistic details. We found that the strength of hydrogen bond network between inhibitor and the protease takes crucial roles in the dissociation process. We showed that the hydrogen bond network in the cyclic urea inhibitors AHA001/XK263 is less stable than that of the approved inhibitor ABT538 because of their large differences in the structures of the networks. In the cyclic urea inhibitor bound complex, the hydrogen bonds often distribute at the flap tips and the active site. In contrast, there are additional accessorial hydrogen bonds formed at the lateral sides of the flaps and the active site in the ABT538 bound complex, which take crucial roles in stabilizing the hydrogen bond network. In addition, the water molecule W301 also plays important roles in stabilizing the hydrogen bond network through its flexible movement by acting as a collision buffer and helping the rebinding of hydrogen bonds at the flap tips. Because of its high stability, the hydrogen bond network of ABT538 complex can work together with the hydrophobic clusters to resist the dissociation, resulting in much lower dissociation rate constant than those of cyclic urea inhibitor complexes. This study may provide useful guidelines for design of novel potent inhibitors with optimized interactions.     

Probing the transducin nucleotide binding site with GDP analogues

An affinity study between the G protein of the visual photoreceptor, transducin, and eight different non-hydrolyzable GDP analogues is described. Imidodiphosphate derivatives have been shown to exhibit good affinities to transducin. This very important heterotrimeric G protein is shown to be highly restrictive with regard to structural modifications of the nucleotide at the pyrophosphate moiety, at the 3' position on ribose, as well as at the N1 position of the purine.     

Probing the S1 specificity pocket of the aminopeptidases that generate antigenic peptides

ERAP1 (endoplasmic reticulum aminopeptidase 1), ERAP2 and IRAP (insulin-regulated aminopeptidase) are three homologous enzymes that play critical roles in the generation of antigenic peptides. These aminopeptidases excise amino acids from N-terminally extended precursors of antigenic peptides in order to generate the correct length epitopes for binding on to MHC class I molecules. The specificity of these peptidases can affect antigenic peptide selection, but has not yet been investigated in detail. In the present study we utilized a collection of 82 fluorigenic substrates to define a detailed selectivity profile for each of the three enzymes and to probe structural and functional features of the S1 (primary specificity) pocket. Molecular modelling of the three S1 pockets reveals substrate-enzyme interactions that are critical determinants for specificity. The substrate selectivity profiles suggest that IRAP largely combines the S1 specificity of ERAP1 and ERAP2, consistent with its proposed biological function. IRAP, however, does not achieve this dual specificity by simply combining structural features of ERAP1 and ERAP2, but rather by an unique amino acid change at position 541. The results of the present study provide insights on antigenic peptide selection and may prove valuable in designing selective inhibitors or activity markers for this class of enzymes.     

Free energy calculations on binding and catalysis by alpha-lytic protease: the role of substrate size in the P1 pocket

We present free energy calculations using molecular dynamics on different substrates of alpha-lytic protease in the gas phase, in solution, while forming a noncovalent Michaelis complex with the enzyme, and in a tetrahedral structure representing a transition state/intermediate for acylation by the enzyme. Various P1 substrates were studied, with P1 = Gly, Ala, Val, and Leu. In qualitative agreement with experiment, the enzyme was calculated to bind and catalyze most effectively substrates with P1 = Ala over those with P1 = Gly, Val or Leu. Also, the calculated relative solvation free energies of Gly----Ala and Ala----Val were in qualitative agreement with experimental values in corresponding model systems. However, the level of quantitative agreement with experiment achieved in our earlier study of relative binding and catalysis of native subtilisin and an Asn-155----Ala mutant was not achieved. We surmise that this is due to the greater difficulty in quantitatively simulating effects that are predominantly van der Waals and hydrophobic compared to those that are hydrogen bonding/electrostatic.     

A poke in the eye: inhibiting HIV-1 protease through its flap-recognition pocket

A novel mechanism of inhibiting HIV-1 protease (HIVp) is presented. Using computational solvent mapping to identify complementary interactions and the Multiple Protein Structure method to incorporate protein flexibility, we generated a receptor-based pharmacophore model of the flexible flap region of the semiopen, apo state of HIVp. Complementary interactions were consistently observed at the base of the flap, only within a cleft with a specific structural role. In the closed, bound state of HIVp, each flap tip docks against the opposite monomer, occupying this cleft. This flap-recognition site is filled by the protein and cannot be identified using traditional approaches based on bound, closed structures. Virtual screening and dynamics simulations show how small molecules can be identified to complement this cleft. Subsequent experimental testing confirms inhibitory activity of this new class of inhibitor. This may be the first new inhibitor class for HIVp since dimerization inhibitors were introduced 17 years ago.     

Evaluation of the substrate envelope hypothesis for inhibitors of HIV-1 protease

Crystallographic data show that various substrates of HIV protease occupy a remarkably uniform region within the binding site; this region has been termed the substrate envelope. It has been suggested that an inhibitor that fits within the substrate envelope should tend to evade viral resistance because a protease mutation that reduces the affinity of the inhibitor will also tend to reduce the affinity of substrate, and will hence decrease the activity of the enzyme. Accordingly, inhibitors that fit the substrate envelope better should be less susceptible to clinically observed resistant mutations, since these must also allow substrates to bind. The present study describes a quantitative measure of the volume of a bound inhibitor falling outside the substrate envelope, and observes that this quantity correlates with the inhibitor's losses in affinity to clinically relevant mutants. This measure may thus be useful as a penalty function in the design of robust HIV protease inhibitors.     

Comprehensive mutagenesis of HIV-1 protease: a computational geometry approach

A computational geometry technique based on Delaunay tessellation of protein structure, represented by C(alpha) atoms, is used to study effects of single residue mutations on sequence-structure compatibility in HIV-1 protease. Profiles of residue scores derived from the four-body statistical potential are constructed for all 1881 mutants of the HIV-1 protease monomer and compared with the profile of the wild-type protein. The profiles for an isolated monomer of HIV-1 protease and the identical monomer in a dimeric state with an inhibitor are analyzed to elucidate changes to structural stability. Protease residues shown to undergo the greatest impact are those forming the dimer interface and flap region, as well as those known to be involved in inhibitor binding.     

Search for substrate-based inhibitors fitting the S2' space of malarial aspartic protease plasmepsin II.

Plasmepsin (Plm) has been identified as an important target for the development of new antimalarial drugs, since its inhibition leads to the starvation of Plasmodium falciparum. A series of substrate-based dipeptide-type Plm II inhibitors containing the hydroxymethylcarbonyl isostere as a transition-state mimic were synthesized. The general design principle was provision of a conformationally restrained hydroxyl group (corresponding to the set residue at the P2' position in native substrates) and a bulky unit to fit the S2' pocket.