Primary structure and properties of the cathepsin G/chymotrypsin inhibitor from the larval hemolymph of Apis mellifera.

A member of the Ascaris inhibitor family exhibiting anti-cathepsin G and anti-chymotrypsin activity was purified from the larval hemolymph of the honey bee (Apis mellifera). Three forms of the inhibitor, designated AMCI 1-3, were isolated using gel filtration and anion-exchange chromatographies followed by reverse-phase HPLC. The amino-acid analyses indicated that AMCI-1 and AMCI-2 have an identical composition whereas AMCI-3 is shorter by two residues (Thr, Arg). All three forms contain as many as 10 cysteine residues and lack tryptophan, tyrosine, and histidine. The sequence of the isoinhibitors showed that the major form (AMCI-1) consisting of 56 amino-acid residues was a single-chain protein of molecular mass 5972 Da, whereas the other two forms were two-chain proteins with a very high residue identity. The AMCI-2 appeared to be derived from AMCI-1, as a result of the Lys24-Thr25 peptide bond splitting, while AMCI-3 was truncated at its N-terminus by the dipeptide Thr25-Arg26. The association constants for the binding of bovine alpha-chymotrypsin to all purified forms of the inhibitor were high and nearly identical, ranging from 4.8 x 10(10) M-1 for AMCI-1 to 2.7 x 10(9) M-1 for AMCI-3. The sensitivity of cathepsin G to inhibition by each inhibitor was different. Only the association constant for the interaction of this enzyme with AMCI-1 was high (2 x 10(8) M-1) whereas those for AMCI-2 and AMCI-3 were significantly lower, and appeared to be 3.7 x 10(7) M-1 and 4.5 x 10(6) M-1, respectively. The reactive site of the inhibitor, as identified by cathepsin G degradation and chemical modification, was found to be at Met30-Gln31. A search in the Protein Sequence Swiss-Prot databank revealed a significant degree of identity (44%) between the primary structure of AMCI and the trypsin isoinhibitor from Ascaris sp (ATI). On the basis of the cysteine residues alignment, the position of the reactive site as well as some sequence homology, the cathepsin G/chymotrypsin inhibitor from larval hemolymph of the honey bee may be considered to be a member of the Ascaris inhibitor family.     

Solution structure of BSTI: a new trypsin inhibitor from skin secretions of Bombina bombina

The three-dimensional solution structure of BSTI, a trypsin inhibitor from the European frog Bombina bombina, has been solved using (1)H NMR spectroscopy. The 60 amino acid protein contains five disulfide bonds, which were unambiguously determined to be Cys (4--38), Cys (13--34), Cys (17--30), Cys (21--60), and Cys (40--54) by experimental restraints and subsequent structure calculations. The main elements of secondary structure are four beta-strands, arranged as two small antiparallel beta-sheets. The overall fold of BSTI is disk shaped and is characterized by the lack of a hydrophobic core. The presumed active site is located on a loop comprising residues 21--34, which is a relatively disordered region similar to that seen in many other protease inhibitors. However, the overall fold is different to other known protease inhibitors with the exception of a small family of inhibitors isolated from nematodes of the family Ascaris and recently also from the haemolymph of Apis mellifera. BSTI may thus be classified as a new member of this recently discovered family of protease inhibitors.     

A structural and thermodynamic escape mechanism from a drug resistant mutation of the HIV-1 protease

The efficacy of HIV-1 protease inhibition therapies is often compromised by the appearance of mutations in the protease molecule that lower the binding affinity of inhibitors while maintaining viable catalytic activity and substrate affinity. The V82F/I84V double mutation is located within the binding site cavity and affects all protease inhibitors in clinical use. KNI-764, a second-generation inhibitor currently under development, maintains significant potency against this mutation by entropically compensating for enthalpic losses, thus minimizing the loss in binding affinity. KNI-577 differs from KNI-764 by a single functional group critical to the inhibitor response to the protease mutation. This single difference changes the response of the two inhibitors to the mutation by one order of magnitude. Accordingly, a structural understanding of the inhibitor response will provide important guidelines for the design of inhibitors that are less susceptible to mutations conveying drug resistance. The structures of the two compounds bound to the wild type and V82F/I84V HIV-1 protease have been determined by X-ray crystallography at 2.0 A resolution. The presence of two asymmetric functional groups, linked by rotatable bonds to the inhibitor scaffold, allows KNI-764 to adapt to the mutated binding site cavity more readily than KNI-577, with a single asymmetric group. Both inhibitors lose about 2.5 kcal/mol in binding enthalpy when facing the drug-resistant mutant protease; however KNI-764 gains binding entropy while KNI-577 loses binding entropy. The gain in binding entropy by KNI-764 accounts for its low susceptibility to the drug-resistant mutation. The heat capacity change associated with binding becomes more negative when KNI-764 binds to the mutant protease, consistent with increased desolvation. With KNI-577, the opposite effect is observed. Structurally, the crystallographic B factors increase for KNI-764 when it is bound to the drug-resistant mutant. The opposite is observed for KNI-577. Consistent with these observations, it appears that KNI-764 is able to gain binding entropy by a two-fold mechanism: it gains solvation entropy by burying itself deeper within the binding pocket and gains conformational entropy by losing interaction with the protease.     

Comparison of SARS and NL63 papain-like protease binding sites and binding site dynamics: inhibitor design implications

The human severe acute respiratory syndrome coronavirus (SARS-CoV) and the NL63 coronaviruses are human respiratory pathogens for which no effective antiviral treatment exists. The papain-like cysteine proteases encoded by the coronavirus (SARS-CoV: PLpro; NL63: PLP1 and PLP2) represent potential targets for antiviral drug development. Three recent inhibitor-bound PLpro structures highlight the role of an extremely flexible six-residue loop in inhibitor binding. The high binding site plasticity is a major challenge in computational drug discovery/design efforts. From conventional molecular dynamics and accelerated molecular dynamics (aMD) simulations, we find that with conventional molecular dynamics simulation, PLpro translationally samples the open and closed conformation of BL2 loop on a picosecond-nanosecond timescale but does not reproduce the peptide bond inversion between loop residues Tyr269 and Gln270 that is observed on inhibitor GRL0617 binding. Only aMD simulation, starting from the closed loop conformation, reproduced the 180° ?-ψ dihedral rotation back to the open loop state. The Tyr-Gln peptide bond inversion appears to involve a progressive conformational change of the full loop, starting at one side, and progressing to the other. We used the SARS-CoV apo X-ray structure to develop a model of the NL63-PLP2 catalytic site. Superimposition of the PLP2 model on the PLpro X-ray structure identifies binding site residues in PLP2 that contribute to the distinct substrate cleavage site specificities between the two proteases. The topological and electrostatic differences between the two protease binding sites also help explain the selectivity of non-covalent PLpro inhibitors.     

The effect of inhibitor binding on the structural stability and cooperativity of the HIV-1 protease.

The effects of the peptide inhibitor acetyl pepstatin on the structural stability of the HIV-1 protease have been measured by high sensitivity calorimetric techniques. At 25 degrees C and pH 3.6, acetyl pepstatin binds to HIV-1 protease with an affinity of 1.6 x 10(7 )M-1 and an enthalpy of 7.3 +/- 0.5 kcal/mol, indicating that binding is not favored enthalpically and that the favorable Gibbs energy originates from a large positive entropy. Since the binding of acetyl pepstatin is associated with a negative change in heat capacity (-450 cal/K*mol) the association reaction becomes enthalpically favored at temperatures higher than 40 degrees C. The presence of the inhibitor stabilizes the dimeric structure of the protease in a fashion that can be quantitatively described by a set of thermodynamic linkage equations. The combination of titration and differential scanning calorimetry provides an accurate way of determining binding constants for high affinity inhibitors that cannot be determined by titration calorimetry alone. A structure-based thermodynamic analysis of the binding process indicates that the stabilization effect is not distributed uniformly throughout the protease molecule. The binding of the inhibitor selectively stabilizes those conformational states in which the binding site is formed, triggering a redistribution of the state probabilities in the ensemble of conformations populated under native conditions. As a result, the stability constants for individual residues do not exhibit the same change in magnitude upon inhibitor binding. Residues in certain areas of the protein are affected significantly whereas residues in other areas are not affected at all. In particular, inhibitor binding has a significant effect on those regions that define the binding site, especially the flap region which becomes structurally stable as a result of the additional binding free energy. The induced stabilization propagates to regions not in direct contact with the inhibitor, particularly to the strand between residues Pro9 and Ala22 and the helix between Arg87 and Gly94. On the other hand, the stability of the strand between Asp60 and Leu76 is not significantly affected by inhibitor binding. The structural distribution of binding effects define cooperative pathways within the protease molecule. Proteins 1999;36:147-156.     

Binding to the open conformation of HIV-1 protease

A recent crystal structure of HIV-1 protease (HIVp) was the first to experimentally observe a ligand targeting an open-flap conformation. Researchers studying a symmetric pyrrolidine inhibitor found that two ligands cocrystallized with the protease, forcing an unusual configuration and unique crystallographic contacts. One molecule is centered in the traditional binding site (α pose) and the other binds between the flaps (β pose). The ligands stack against each other in a region termed the "eye" site. Ligands bound to the eye site should prevent flap closure, but it is unclear if the pyrrolidine inhibitors or the crystal packing are causing the open state. Molecular dynamics simulations were used to examine the solution-state behavior of three possible binding modes: the ternary complex of HIVp+αβ and the binary complexes, HIVp+α and HIVp+β. We show that HIVp+α is the most stable of the three states. During conformational sampling, α takes an asymmetric binding pose, with one naphthyl ring occupying the eye site and the other reoriented down to occupy positions seen with traditional inhibitors. This finding supports previous studies that reveal a requirement for asymmetric binding at the eye site. In fact, if the α pose is modified to splay both naphthyl rings across the binding site like traditional inhibitors, one ring consistently flips to occupy the eye site. Our simulations reveal that interactions to the eye site encourage a conformationally restrained state, and understanding those contacts may aid the design of ligands to specifically target alternate conformations of the protease.     

Inhibition of yeast mitochondrial F1-ATPase, F0F1-ATPase and submitochondrial particles by rhodamines and ethidium bromide

ATP hydrolysis by F1-ATPase is strongly inhibited by cationic rhodamines; neutral rhodamines are very poor inhibitors. Rhodamine 6G is a noncompetitive inhibitor of purified F0F1-ATPase and submitochondrial particles, however, an uncompetitive inhibitor of F1-ATPase (KI approximately equal to 2.4 microM for all three enzyme forms). Ethidium bromide is a noncompetitive inhibitor of F0F1-ATPase, submitochondrial particles and also F1-ATPase (KI approximately equal to 270 microM). Neither of the inhibitors affects the negative cooperativity (nH approximately equal to 0.7). The non-identical binding sites for rhodamine 6G and ethidium bromide are located on the F1-moiety and are topologically distinct from the catalytic site. Binding of the inhibitors prevents the conformational changes essential for energy transduction. It is concluded that the inhibitor binding sites are involved in proton translocation. In F1-ATPase, binding of MgATP at a catalytic site causes conformational changes, which allosterically induce the correct structure of the rhodamine 6G binding site. In F0F1-ATPase, this conformation of the F1-moiety exists a priori, due to allosteric interactions with F0-subunits. The binding site for ethidium bromide on F1-ATPase does not require substrate binding at the catalytic site and is not affected by F0F1-subunit interactions.     

Conformational flexibility and binding energy profile of c-Abl tyrosine kinase complexed with Imatinib: an insight from MD study

Structural biology of kinase and in particular of tyrosine kinase has given detailed insights into the intrinsic flexibility of the catalytic domain and has provided a rational basis for obtaining selective inhibitors. In this paper, we have studied the conformational flexibility of c-Abl tyrosine kinase complexed with Imatinib (STI), in the presence of TIP3P water in physiological conditions at neutral pH. The conformational studies suggest that the flexibility of activation loop is responsible to facilitate the nucleotide binding and release. Owing to the conformational adaptability, adenosine triphosphate (ATP) binds at a particular site in the loop region of the tyrosine kinase. The molecular mechanics Poisson–Boltzmann surface area methods are analysed, as is a free-energy pathways method, which shows the stable binding with free energy ? 6.04 kcal/mol for STI. The binding energy calculated by the Sietraj method is approximately the same as the experimental binding energy of STI with c-Abl kinase. It is suggested that the conserved glutamic acid and lysine residues are necessary for the stability and optimum activity of inhibitor. This study may be helpful in rational drug designing of new kinase inhibitors.     

Hydrophobic sliding: a possible mechanism for drug resistance in human immunodeficiency virus type 1 protease

Hydrophobic residues outside the active site of HIV-1 protease frequently mutate in patients undergoing protease inhibitor therapy; however, the mechanism by which these mutations confer drug resistance is not understood. From analysis of molecular dynamics simulations, 19 core hydrophobic residues appear to facilitate the conformational changes that occur in HIV-1 protease. The hydrophobic core residues slide by each other, exchanging one hydrophobic van der Waal contact for another, with little energy penalty, while maintaining many structurally important hydrogen bonds. Such hydrophobic sliding may represent a general mechanism by which proteins undergo conformational changes. Mutation of these residues in HIV-1 protease would alter the packing of the hydrophobic core, affecting the conformational flexibility of the protease. Therefore these residues impact the dynamic balance between processing substrates and binding inhibitors, and thus contribute to drug resistance.     

The mechanism of inhibition of antibody-based inhibitors of membrane-type serine protease 1 (MT-SP1)

The mechanisms of inhibition of two novel scFv antibody inhibitors of the serine protease MT-SP1/matriptase reveal the basis of their potency and specificity. Kinetic experiments characterize the inhibitors as extremely potent inhibitors with K(I) values in the low picomolar range that compete with substrate binding in the S1 site. Alanine scanning of the loops surrounding the protease active site provides a rationale for inhibitor specificity. Each antibody binds to a number of residues flanking the active site, forming a unique three-dimensional binding epitope. Interestingly, one inhibitor binds in the active site cleft in a substrate-like manner, can be processed by MT-SP1 at low pH, and is a standard mechanism inhibitor of the protease. The mechanisms of inhibition provide a rationale for the effectiveness of these inhibitors, and suggest that the development of specific antibody-based inhibitors against individual members of closely related enzyme families is feasible, and an effective way to develop tools to tease apart complex biological processes.     

Complexation of two proteic insect inhibitors to the active site of chymotrypsin suggests decoupled roles for binding and selectivity

The crystal structures of two homologous inhibitors (PMP-C and PMP-D2v) from the insect Locusta migratoria have been determined in complex with bovine alpha-chymotrypsin at 2.1- and 3.0-A resolution, respectively. PMP-C is a potent bovine alpha-chymotrypsin inhibitor whereas native PMP-D2 is a weak inhibitor of bovine trypsin. One unique mutation at the P1 position converts PMP-D2 into a potent bovine alpha-chymotrypsin inhibitor. The two peptides have a similar overall conformation, which consists of a triple-stranded antiparallel beta-sheet connected by three disulfide bridges, thus defining a novel family of serine protease inhibitors. They have in common the protease interaction site, which is composed of the classical protease binding loop (position P5 to P'4, corresponding to residues 26-34) and of an internal segment (residues 15-18), held together by two disulfide bridges. Structural divergences between the two inhibitors result in an additional interaction site between PMP-D2v (position P10 to P6, residues 21-25) and the residues 172-175 of alpha-chymotrypsin. This unusual interaction may be responsible for species selectivity. A careful comparison of data on bound and free inhibitors (from this study and previous NMR studies, respectively) suggests that complexation to the protease stabilizes the flexible binding loop (from P5 to P'4).     

Accessory mutations balance the marginal stability of the HIV-1 protease in drug resistance

The HIV-1 protease is a major target of inhibitor drugs in AIDS therapies. The therapies are impaired by mutations of the HIV-1 protease that can lead to resistance to protease inhibitors. These mutations are classified into major mutations, which usually occur first and clearly reduce the susceptibility to protease inhibitors, and minor, accessory mutations that occur later and individually do not substantially affect the susceptibility to inhibitors. Major mutations are predominantly located in the active site of the HIV-1 protease and can directly interfere with inhibitor binding. Minor mutations, in contrast, are typically located distal to the active site. A central question is how these distal mutations contribute to resistance development. In this article, we present a systematic computational investigation of stability changes caused by major and minor mutations of the HIV-1 protease. As most small single-domain proteins, the HIV-1 protease is only marginally stable. Mutations that destabilize the folded, active state of the protease therefore can shift the conformational equilibrium towards the unfolded, inactive state. We find that the most frequent major mutations destabilize the HIV-1 protease, whereas roughly half of the frequent minor mutations are stabilizing. An analysis of protease sequences from patients in treatment indicates that the stabilizing minor mutations are frequently correlated with destabilizing major mutations, and that highly resistant HIV-1 proteases exhibit significant fractions of stabilizing mutations. Our results thus indicate a central role of minor mutations in balancing the marginal stability of the protease against the destabilization induced by the most frequent major mutations.     

Requirement for hydrophobic Phe residues in Pleurotus ostreatus proteinase A inhibitor 1 for stable inhibition

Pleurotus ostreatus proteinase A inhibitor 1 (POIA1) has been shown to be unique among the various serine protease inhibitors in that its C-terminal region appears to be the reactive site responsible for its inhibitory action toward proteases. To investigate in more detail the mechanism of inhibition by POIA1, we have been studying its structural requirements for stable inhibition of proteases. In this study, we focused on hydrophobic Phe residues, which are generally located in the interior of protein molecules. A Phe-->Ala replacement at position 44 or 56 was introduced into a 'parent' mutant of POIA1 that had been converted into a strong and resistant inhibitor of subtilisin BPN' by replacement of its six C-terminal residues with those of the propeptide of subtilisin BPN' and the effects on inhibitory properties and structural stability were examined. Both of the mutated POIA1 molecules not only were found to exhibit decreased ability to bind to subtilisin BPN' (80-fold for the F44A mutant and 13-fold for the F56A mutant), but were also converted to temporary inhibitors that were degraded by the protease. The structural stability of the mutated POIA1 was also lowered, as shown by a 13 degrees C decrease in melting temperature for the F56A mutant. In particular, the F44A mutant was found to lose its tertiary structure, as judged from the circular dichroism spectrum, demonstrating that Phe44 is a strict requirement for structural formation by the POIA1 molecule. These results clearly indicate that stabilization of POIA1 by hydrophobic residues in its molecular interior is required for stable inhibition of the protease. This requirement for a stable tertiary structure is shared with other serine protease inhibitors, but other structural requirements seem to differ, in that strong binding with the protease is required for POIA1 whereas conformational rigidity around the reactive site is essential for many other protease inhibitors.     

Solution Structure of a Novel C<Subscript>2</Subscript>-Symmetrical Bifunctional Bicyclic Inhibitor Based on SFTI-1

A novel bifunctional bicyclic inhibitor has been created that combines features both from the Bowman–Birk inhibitor (BBI) proteins, which have two distinct inhibitory sites, and from sunflower trypsin inhibitor-1 (SFTI-1), which has a compact bicyclic structure. The inhibitor was designed by fusing together a pair of reactive loops based on a sequence derived from SFTI-1 to create a backbone-cyclized disulfide-bridged 16-mer peptide. This peptide has two symmetrically spaced trypsin binding sites. Its synthesis and biological activity have been reported in a previous communication [Jaulent and Leatherbarrow, 2004, PEDS 17, 681]. In the present study we have examined the three-dimensional structure of the molecule. We find that the new inhibitor, which has a symmetrical 8-mer half-cystine CTKSIPP′I′ motif repeated through a C₂ symmetry axis also shows a complete symmetry in its three-dimensional structure. Each of the two loops adopts the expected canonical conformation common to all BBIs as well as SFTI-1. We also find that the inhibitor displays a strong and unique structural identity, with a notable lack of minor conformational isomers that characterise most reactive site loop mimics examined to date as well as SFTI-1. This suggests that the presence of the additional cyclic loop acts to restrict conformational mobility and that the deliberate introduction of cyclic symmetry may offer a general route to locking the conformation of β-hairpin structures. Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Thermodynamic basis of resistance to HIV-1 protease inhibition: calorimetric analysis of the V82F/I84V active site resistant mutant

One of the most serious side effects associated with the therapy of HIV-1 infection is the appearance of viral strains that exhibit resistance to protease inhibitors. The active site mutant V82F/I84V has been shown to lower the binding affinity of protease inhibitors in clinical use. To identify the origin of this effect, we have investigated the binding thermodynamics of the protease inhibitors indinavir, ritonavir, saquinavir, and nelfinavir to the wild-type HIV-1 protease and to the V82F/I84V resistant mutant. The main driving force for the binding of all four inhibitors is a large positive entropy change originating from the burial of a significant hydrophobic surface upon binding. At 25 degrees C, the binding enthalpy is unfavorable for all inhibitors except ritonavir, for which it is slightly favorable (-2.3 kcal/mol). Since the inhibitors are preshaped to the geometry of the binding site, their conformational entropy loss upon binding is small, a property that contributes to their high binding affinity. The V82F/I84V active site mutation lowers the affinity of the inhibitors by making the binding enthalpy more positive and making the entropy change slightly less favorable. The effect on the enthalpy change is, however, the major one. The predominantly enthalpic effect of the V82F/I84V mutation is consistent with the idea that the introduction of the bulkier Phe side chain at position 82 and the Val side chain at position 84 distort the binding site and weaken van der Waals and other favorable interactions with inhibitors preshaped to the wild-type binding site. Another contribution of the V82F/I84V to binding affinity originates from an increase in the energy penalty associated with the conformational change of the protease upon binding. The V82F/I84V mutant is structurally more stable than the wild-type protease by about 1.4 kcal/mol. This effect, however, affects equally the binding affinity of substrate and inhibitors.     

Conformational change in rhomboid protease GlpG induced by inhibitor binding to its S' subsites

Rhomboid protease conducts proteolysis inside the hydrophobic environment of the membrane. The conformational flexibility of the protease is essential for the enzyme mechanism, but the nature of this flexibility is not completely understood. Here we describe the crystal structure of rhomboid protease GlpG in complex with a phosphonofluoridate inhibitor, which is covalently bonded to the catalytic serine and extends into the S' side of the substrate binding cleft. Inhibitor binding causes subtle but extensive changes in the membrane protease. Many transmembrane helices tilt and shift positions, and the gap between S2 and S5 is slightly widened so that the inhibitor can bind between them. The side chain of Phe-245 from a loop (L5) that acts as a cap rotates and uncovers the opening of the substrate binding cleft to the lipid bilayer. A concurrent turn of the polypeptide backbone at Phe-245 moves the rest of the cap and exposes the catalytic serine to the aqueous solution. This study, together with earlier crystallographic investigation of smaller inhibitors, suggests a simple model for explaining substrate binding to rhomboid protease.     

The effect of prime-site occupancy on the hepatitis C virus NS3 protease structure

We recently reported a new class of inhibitors of the chymotrypsin-like serine protease NS3 of the hepatitis C virus. These inhibitors exploit the binding potential of the S' site of the protease, which is not generally used by the natural substrates. The effect of prime-site occupancy was analyzed by circular dichroism spectroscopy and limited proteolysis-mass spectrometry. Generally, nonprime inhibitors cause a structural change in NS3. Binding in the S' site produces additional conformational changes with different binding modes, even in the case of the NS3/4A cofactor complex. Notably, inhibitor binding either in the S or S' site also has profound effects on the stabilization of the protease. In addition, the stabilization propagates to regions not in direct contact with the inhibitor. In particular, the N-terminal region, which according to structural studies is endowed with low structural stability and is not stabilized by nonprime inhibitors, was now fully protected from proteolytic degradation. From the perspective of drug design, P-P' inhibitors take advantage of binding pockets, which are not exploited by the natural HCV substrates; hence, they are an entry point for a novel class of NS3/4A inhibitors. Here we show that binding of each inhibitor is associated with a specific structural rearrangement. The development of a range of inhibitors belonging to different classes and an understanding of their interactions with the protease are required to address the issue of the most likely outcome of viral protease inhibitor therapy, that is, viral resistance.     

Multi-drug resistance profile of PR20 HIV-1 protease is attributed to distorted conformational and drug binding landscape: molecular dynamics insights

The PR20 HIV-1 protease, a variant with 20 mutations, exhibits high levels of multi-drug resistance; however, to date, there has been no report detailing the impact of these 20 mutations on the conformational and drug binding landscape at a molecular level. In this report, we demonstrate the first account of a comprehensive study designed to elaborate on the impact of these mutations on the dynamic features as well as drug binding and resistance profile, using extensive molecular dynamics analyses. Comparative MD simulations for the wild-type and PR20 HIV proteases, starting from bound and unbound conformations in each case, were performed. Results showed that the apo conformation of the PR20 variant of the HIV protease displayed a tendency to remain in the open conformation for a longer period of time when compared to the wild type. This led to a phenomena in which the inhibitor seated at the active site of PR20 tends to diffuse away from the binding site leading to a significant change in inhibitor–protein association. Calculating the per-residue fluctuation (RMSF) and radius of gyration, further validated these findings. MM/GBSA showed that the occurrence of 20 mutations led to a drop in the calculated binding free energies (ΔG_bind) by ~25.17 kcal/mol and ~5 kcal/mol for p2-NC, a natural peptide substrate, and darunavir, respectively, when compared to wild type. Furthermore, the residue interaction network showed a diminished inter-residue hydrogen bond network and changes in inter-residue connections as a result of these mutations. The increased conformational flexibility in PR20 as a result of loss of intra- and inter-molecular hydrogen bond interactions and other prominent binding forces led to a loss of protease grip on ligand. It is interesting to note that the difference in conformational flexibility between PR20 and WT conformations was much higher in the case of substrate-bound conformation as compared to DRV. Thus, developing analogues of DRV by retaining its key pharmacophore features will be the way forward in the search for novel protease inhibitors against multi-drug resistant strains.     

HIV-1 protease inhibitors: enthalpic versus entropic optimization of the binding affinity

Existing experimental as well as computational screening methods select potential ligands or drug candidates on the basis of binding affinity. Since the binding affinity is a function of the enthalpy (DeltaH) and entropy (DeltaS) changes, it is apparent that improved binding can be achieved in different ways: by optimizing DeltaH, DeltaS, or a combination of both. However, the behavior of enthalpically or entropically optimized inhibitors is fundamentally different, including their response to mutations that may elicit drug resistance. In the design of HIV-1 protease inhibitors, high binding affinity has usually been achieved by preshaping lead compounds to the geometry of the binding site and by incorporating a high degree of hydrophobicity. The thermodynamic consequence of that approach is that the binding affinity of the resulting inhibitors becomes entropically favorable but enthalpically unfavorable. Specifically, the resulting high binding affinity is due to an increased solvation entropy (hydrophobic effect) combined with a reduced loss of conformational entropy of the inhibitor upon binding (structural rigidity). Here we report that tripeptide inhibitors derived from the transframe region of Gag-Pol (Glu-Asp-Leu and Glu-Asp-Phe) bind to the HIV-1 protease with a favorable enthalpy change. This behavior is qualitatively different from that of known inhibitors and points to new strategies for inhibitor design. Since the binding affinities of enthalpically favorable and enthalpically unfavorable inhibitors have opposite temperature dependence, it is possible to design fast screening protocols that simultaneously select inhibitors on the basis of affinity and enthalpy.     

Mechanism of action of aspartic proteases. III. Conformational characteristics of HIV-1 protease inhibitor JG-365]

A set of conformations was shown to be characteristic of the free-state spatial structure of substrate-like inhibitor JG-365 for aspartic protease from HIV-1. Among them, the lowest-energy conformations have a folded form of the peptide backbone. The inhibitor has a noncleavable hydroxyethylamine group with an additional chiral center in its structure. Our calculations showed that only the S-isomer of the inhibitor displays conformational characteristics that practically coincide with those of the native substrate for HIV-1 protease. One of the calculated conformations with a completely extended main chain and a relative energy of 9.5 kcal/mol very closely mimics the experimentally observed structure of the inhibitor in the enzyme-inhibitor complex. The realization of this structure is unlikely for a free inhibitor, because it has only a small number of interresidual noncovalent interactions in the extended conformation; these are presumably compensated for by intermolecular interactions at the active site of the enzyme.