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.     

Serine protease mechanism: structure of an inhibitory complex of alpha-lytic protease and a tightly bound peptide boronic acid

The structure of the complex formed between alpha-lytic protease, a serine protease secreted by Lysobacter enzymogenes, and N-tert-butyloxycarbonylalanylprolylvaline boronic acid (Ki = 0.35 nM) has been studied by X-ray crystallography to a resolution of 2.0 A. The active-site serine forms a covalent, nearly tetrahedral adduct with the boronic acid moiety of the inhibitor. The complex is stabilized by seven hydrogen bonds between the enzyme and inhibitor with additional stabilization arising from van der Waals interactions between enzyme and inhibitor side chains and the burying of 330 A2 of hydrophobic surface area. Hydrogen bonding between Asp-102 and His-57 remains intact in the enzyme-inhibitor complex, and His N epsilon 2 is well positioned to donate its hydrogen to the leaving group. Little change in the positions of protease residues was observed on complex formation (root mean square main chain deviation = 0.13 A), suggesting that in its native state the enzyme is complementary to tetrahedral reaction intermediates or to the nearly tetrahedral transition state for the reaction.     

Crystal structure of the anticarcinogenic Bowman-Birk inhibitor from snail medic (Medicago scutellata) seeds complexed with bovine trypsin

The structure of the ternary complex of the anticarcinogenic Bowman-Birk protease inhibitor purified from snail medic (Medicago scutellata) seeds (MSTI) and two molecules of bovine trypsin has been solved by X-ray diffraction analysis of single crystals to a resolution of 2.0 A. This is the highest resolution model of a ternary complex of this type currently available. The two binding loops of the MSTI differ in only one amino acid and have in both cases an arginine in position P1. The distances between the residues of the inhibitor at the binding interface and the trypsin side chains that recognize them are almost identical in the two sites. When compared to the NMR model of the uncomplexed MSTI, the inhibitor in the functional assembly with trypsin shows the largest differences in the two P2' residues. Compared with the similar ternary complex of the soybean trypsin inhibitor, this model shows very small differences in the polypeptide chain of the trypsin binding sites and its largest difference in the area between Asp 26 and His 32 of the MSTI which in the soybean inhibitor has an extra Leu inserted in position 29.     

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.     

Anisotropic dynamics of the JE-2147-HIV protease complex: drug resistance and thermodynamic binding mode examined in a 1.09 A structure

The structure of HIV protease (HIV Pr) bound to JE-2147 (also named AG1776 or KNI-764) is determined here to 1.09 A resolution. This highest-resolution structure for HIV Pr allows refinement of anisotropic displacement parameters (ADPs) for all atoms. Clustering based on the directional information in ADPs defines two sets of subdomains such that within each set, subdomains undergo similar anisotropic motion. These sets are (a) the core of monomer A grouped with both substrate-binding flaps and (b) the core of monomer B coupled to both catalytic aspartates (25A/B). The four-stranded beta-sheet (1-4 A/B and 95-99 A/B) that forms a significant part of the dimer interface exhibits large anisotropic amplitudes that differ from those of the other sets of subdomains. JE-2147 is shown here to be a picomolar inhibitor (K(i) = 41 +/- 18 pM). The structure is used to interpret the mechanism of association of JE-2147, a second-generation inhibitor for which binding is enthalpically driven, with respect to first-generation inhibitors for which binding is predominantly entropically driven [Velazquez-Campoy, A., et al. (2001) Arch. Biochem. Biophys. 390, 169-175]. Relative to the entropically driven inhibitor complexes, the JE-2147-HIV Pr complex exhibits an approximately 0.5 A movement of the substrate flaps in toward the substrate, suggesting a more compatible enthalpically driven association. Domains of the protease identified by clustering of ADPs also suggest a model of enthalpy-entropy compensation for all HIV Pr inhibitors in which dynamic coupling of the flaps is offset by an increased level of motion of the beta-sheet domain of the dimer interface (1-4 A/B and 95-99 A/B).     

Three-dimensional solution structure of the HIV-1 protease complexed with DMP323, a novel cyclic urea-type inhibitor, determined by nuclear magnetic resonance spectroscopy.

The three-dimensional solution structure of the HIV-1 protease homodimer, MW 22.2 kDa, complexed to a potent, cyclic urea-based inhibitor, DMP323, is reported. This is the first solution structure of an HIV protease/inhibitor complex that has been elucidated. Multidimensional heteronuclear NMR spectra were used to assemble more than 4,200 distance and angle constraints. Using the constraints, together with a hybrid distance geometry/simulated annealing protocol, an ensemble of 28 NMR structures was calculated having no distance or angle violations greater than 0.3 A or 5 degrees, respectively. Neglecting residues in disordered loops, the RMS deviation (RMSD) for backbone atoms in the family of structures was 0.60 A relative to the average structure. The individual NMR structures had excellent covalent geometry and stereochemistry, as did the restrained minimized average structure. The latter structure is similar to the 1.8-A X-ray structure of the protease/DMP323 complex (Chang CH et al., 1995, Protein Science, submitted); the pairwise backbone RMSD calculated for the two structures is 1.22 A. As expected, the mismatch between the structures is greatest in the loops that are disordered and/or flexible. The flexibility of residues 37-42 and 50-51 may be important in facilitating substrate binding and product release, because these residues make up the respective hinges and tips of the protease flaps. Flexibility of residues 4-8 may play a role in protease regulation by facilitating autolysis.     

Structural and conformational changes concomitant with the E1-E2 transition in H(+)K(+)-ATPase: a comparative protein modeling study

Comparative modeling studies on conserved regions of the gastric H(+)K(+)-ATPase reveal that the E1-E2 conformational transition induces significant tertiary structural changes while conserving the secondary structure. The residues 516-530 of the cytoplasmic domain and TM10 within the transmembrane (TM) regions undergo maximum tertiary structural changes. The luminal regions exhibit comparatively lesser tertiary structural deviations. Residues 249-304 show maximum secondary structural deviation in the conformational transition. The Cys-815 and Cys-323 residues involved in inhibitor binding are found to have smaller buried side chain areas in the E1 conformation compared to E2. Retention of activity correlates well with the buried side chain area when selected amino acid residues in TM6 are mutated using modeling techniques with bulkier amino acid residues. Conformational specificity for ion binding is corroborated with the fraction of side chains exposed to polar atoms of the residues E345, D826, V340, A341, V343, and E822.     

Determination of local conformation of proteins from 1H-NMR spectroscopy data

A method is proposed to determine conformations of amino acid residues of the protein and effective correlation time tau c from cross-peak intensities in two-dimensional nuclear Overhauser enhancement (NOESY) spectra. The method consists in fitting complete relaxation matrix of dipeptide unit protons to experimental cross-peak intensities by varying phi, psi, chi torsional angles and tau c. To verify the method, NOESY spectra of basic pancreatic trypsin inhibitor (BPTI) were theoretically generated at mixing times tau m = 25-300 ms and tau c = 4 ns and used for local structure determination. The method works well with optimum for measurement of NOE intensities tau m 100-200 ms. As a result, the backbone phi, psi torsion angles were unambiguously determined at tau m = 100 ms for all but Gly residues of BPTI, and chi 1 angles were determined for the majority of side chains. The obtained dipeptide unit conformations are very close to the BPTI crystallographic structure: root mean square deviation (RMSD) of interproton distances within dipeptide units, on the average, is 0.08 A (maximal deviation 0.44 A), and RMSD of phi and psi angles are 18 and 9 degrees, respectively (maximal deviations are 44 and 22 degrees).     

Crystal structures of a multidrug-resistant human immunodeficiency virus type 1 protease reveal an expanded active-site cavity

The goal of this study was to use X-ray crystallography to investigate the structural basis of resistance to human immunodeficiency virus type 1 (HIV-1) protease inhibitors. We overexpressed, purified, and crystallized a multidrug-resistant (MDR) HIV-1 protease enzyme derived from a patient failing on several protease inhibitor-containing regimens. This HIV-1 variant contained codon mutations at positions 10, 36, 46, 54, 63, 71, 82, 84, and 90 that confer drug resistance to protease inhibitors. The 1.8-angstrom (A) crystal structure of this MDR patient isolate reveals an expanded active-site cavity. The active-site expansion includes position 82 and 84 mutations due to the alterations in the amino acid side chains from longer to shorter (e.g., V82A and I84V). The MDR isolate 769 protease "flaps" stay open wider, and the difference in the flap tip distances in the MDR 769 variant is 12 A. The MDR 769 protease crystal complexes with lopinavir and DMP450 reveal completely different binding modes. The network of interactions between the ligands and the MDR 769 protease is completely different from that seen with the wild-type protease-ligand complexes. The water molecule-forming hydrogen bonds bridging between the two flaps and either the substrate or the peptide-based inhibitor are lacking in the MDR 769 clinical isolate. The S1, S1', S3, and S3' pockets show expansion and conformational change. Surface plasmon resonance measurements with the MDR 769 protease indicate higher k(off) rates, resulting in a change of binding affinity. Surface plasmon resonance measurements provide k(on) and k(off) data (K(d) = k(off)/k(on)) to measure binding of the multidrug-resistant protease to various ligands. This MDR 769 protease represents a new antiviral target, presenting the possibility of designing novel inhibitors with activity against the open and expanded protease forms.     

Inhibition of the hepatitis C virus NS3/4A protease. The crystal structures of two protease-inhibitor complexes

The hepatitis C virus NS3 protein contains a serine protease domain with a chymotrypsin-like fold, which is a target for development of therapeutics. We report the crystal structures of this domain complexed with NS4A cofactor and with two potent, reversible covalent inhibitors spanning the P1-P4 residues. Both inhibitors bind in an extended backbone conformation, forming an anti-parallel beta-sheet with one enzyme beta-strand. The P1 residue contributes most to the binding energy, whereas P2-P4 side chains are partially solvent exposed. The structures do not show notable rearrangements of the active site upon inhibitor binding. These results are significant for the development of antivirals.     

Structural studies of FIV and HIV-1 proteases complexed with an efficient inhibitor of FIV protease

Three forms of feline immunodeficiency virus protease (FIV PR), the wild type (wt) and two single point mutants, V59I and Q99V, as well as human immunodeficiency virus type 1 protease (HIV-1 PR), were cocrystallized with the C2-symmetric inhibitor, TL-3. The mutants of FIV PR were designed to replace residues involved in enzyme-ligand interactions by the corresponding HIV-1 PR residues at the structurally equivalent position. TL-3 shows decreased (improved) inhibition constants with these FIV PR mutants relative to wt FIV PR. Despite similar modes of binding of the inhibitor to all PRs (from P3 to P3'), small differences are evident in the conformation of the Phe side chains of TL-3 at the P1 and P1' positions in the complexes with the mutated FIV PRs. The differences mimick the observed binding of TL-3 in HIV-1 PR and correlate with a significant improvement in the inhibition constants of TL-3 with the two mutant FIV PRs. Large differences between the HIV-1 and FIV PR complexes are evident in the binding modes of the carboxybenzyl groups of TL-3 at P4 and P4'. In HIV-1 PR:TL-3, these groups bind over the flap region, whereas in the FIV PR complexes, the rings are located along the major axis of the active site. A significant difference in the location of the flaps in this region of the HIV-1 and FIV PRs correlates with the observed conformational changes in the binding mode of the peptidomimetic inhibitor at the P4 and P4' positions. These findings provide a structural explanation of the observed Ki values for TL-3 with the different PRs and will further assist in the development of improved inhibitors.     

Solution structure of the mature HIV-1 protease monomer: insight into the tertiary fold and stability of a precursor

We present the first solution structure of the HIV-1 protease monomer spanning the region Phe1-Ala95 (PR1-95). Except for the terminal regions (residues 1-10 and 91-95) that are disordered, the tertiary fold of the remainder of the protease is essentially identical to that of the individual subunit of the dimer. In the monomer, the side chains of buried residues stabilizing the active site interface in the dimer, such as Asp25, Asp29, and Arg87, are now exposed to solvent. The flap dynamics in the monomer are similar to that of the free protease dimer. We also show that the protease domain of an optimized precursor flanked by 56 amino acids of the N-terminal transframe region is predominantly monomeric, exhibiting a tertiary fold that is quite similar to that of PR1-95 structure. This explains the very low catalytic activity observed for the protease prior to its maturation at its N terminus as compared with the mature protease, which is an active stable dimer under identical conditions. Adding as few as 2 amino acids to the N terminus of the mature protease significantly increases its dissociation into monomers. Knowledge of the protease monomer structure and critical features of its dimerization may aid in the screening and design of compounds that target the protease prior to its maturation from the Gag-Pol precursor.     

Efficiency of a second-generation HIV-1 protease inhibitor studied by molecular dynamics and absolute binding free energy calculations

A subnanomolar inhibitor of human immunodeficiency virus type 1 (HIV-1) protease, designated QF34, potently inhibits the wild-type and drug-resistant enzyme. To explain its broad activity, the binding of QF34 to the wild-type HIV-1 protease is investigated by molecular dynamics simulations and compared to the binding of two inhibitors that are used clinically, saquinavir (SQV) and indinavir (IDV). Analysis of the flexibility of protease residues and inhibitor segments in the complex reveals that segments of QF34 were more mobile during the dynamics studies than the segments of SQV and IDV. The dynamics of hydrogen bonding show that QF34 forms a larger number of stable hydrogen bonds than the two inhibitors that are used clinically. Absolute binding free energies were calculated with molecular mechanics-generalized Born surface area (MM-GBSA) methodology using three protocols. The most consistent results were obtained using the single-trajectory approach, due to cancellation of errors and inadequate sampling in the separate-trajectory protocols. For all three inhibitors, energy components in favor of binding include van der Waals and electrostatic terms, whereas polar solvation and entropy terms oppose binding. Decomposition of binding energies reveals that more protease residues contribute significantly to the binding of QF34 than to the binding of SQV and IDV. Moreover, contributions from protease main chains and side chains are balanced in the case of QF34 (52:48 ratio, respectively), whereas side chain contributions prevail in both SQV and IDV (main-chain:side-chain ratios of 41:59 and 45:55, respectively). The presented results help explain the ability of QF34 to inhibit multiple resistant mutants and should be considered in the design of broad-specificity second-generation HIV-1 protease inhibitors.     

Role of hydroxyl group and R/S configuration of isostere in binding properties of HIV-1 protease inhibitors.

The crystal structure of the complex between human immunodeficiency virus type 1 (HIV-1) protease and a peptidomimetic inhibitor of ethyleneamine type has been refined to R factor of 0.178 with diffraction limit 2.5 A. The peptidomimetic inhibitor Boc-Phe-Psi[CH2CH2NH]-Phe-Glu-Phe-NH2 (denoted here as OE) contains the ethyleneamine replacement of the scissile peptide bond. The inhibitor lacks the hydroxyl group which is believed to mimic tetrahedral transition state of proteolytic reaction and thus is suspected to be necessary for good properties of peptidomimetic HIV-1 protease inhibitors. Despite the missing hydroxyl group the inhibition constant of OE is 1.53 nm and it remains in the nanomolar range also towards several available mutants of HIV-1 protease. The inhibitor was found in the active site of protease in an extended conformation with a unique hydrogen bond pattern different from hydroxyethylene and hydroxyethylamine inhibitors. The isostere nitrogen forms a hydrogen bond to one catalytic aspartate only. The other aspartate forms two weak hydrogen bridges to the ethylene group of the isostere. A comparison with other inhibitors of this series containing isostere hydroxyl group in R or S configuration shows different ways of accommodation of inhibitor in the active site. Special attention is devoted to intermolecular contacts between neighbouring dimers responsible for mutual protein adhesion and for a special conformation of Met46 and Phe53 side chains not expected for free protein in water solution.     

X-ray crystal structure of the protease inhibitor domain of Alzheimer's amyloid beta-protein precursor

Alzheimer's amyloid beta-protein precursor contains a Kunitz protease inhibitor domain (APPI) potentially involved in proteolytic events leading to cerebral amyloid deposition. To facilitate the identification of the physiological target of the inhibitor, the crystal structure of APPI has been determined and refined to 1.5-A resolution. Sequences in the inhibitor-protease interface of the correct protease target will reflect the molecular details of the APPI structure. While the overall tertiary fold of APPI is very similar to that of the Kunitz inhibitor BPTI, a significant rearrangement occurs in the backbone conformation of one of the two protease binding loops. A number of Kunitz inhibitors have similar loop sequences, indicating the structural alteration is conserved and potentially an important determinant of inhibitor specificity. In a separate region of the protease binding loops, APPI side chains Met-17 and Phe-34 create an exposed hydrophobic surface in place of Arg-17 and Val-34 in BPTI. The restriction this change places on protease target sequences is seen when the structure of APPI is superimposed on BPTI complexed to serine proteases, where the hydrophobic surface of APPI faces a complementary group of nonpolar side chains on kallikrein A versus polar side chains on trypsin.     

Folded monomer of HIV-1 protease

The mature human immunodeficiency virus type 1 protease rapidly folds into an enzymatically active stable dimer, exhibiting an intricate interplay between structure formation and dimerization. We now show by NMR and sedimentation equilibrium studies that a mutant protease containing the R87K substitution (PR(R87K)) within the highly conserved Gly(86)-Arg(87)-Asn(88) sequence forms a monomer with a fold similar to a single subunit of the dimer. However, binding of the inhibitor DMP323 to PR(R87K) produces a stable dimer complex. Based on the crystal structure and our NMR results, we postulate that loss of specific interactions involving the side chain of Arg(87) destabilizes PR(R87K) by perturbing the inner C-terminal beta-sheet (residues 96-99 from each monomer), a region that is sandwiched between the two beta-strands formed by the N-terminal residues (residues 1-4) in the mature protease. We systematically examined the folding, dimerization, and catalytic activities of mutant proteases comprising deletions of either one of the terminal regions (residues 1-4 or 96-99) or both. Although both N- and C-terminal beta-strands were found to contribute to dimer stability, our results indicate that the inner C-terminal strands are absolutely essential for dimer formation. Knowledge of the monomer fold and regions critical for dimerization may aid in the rational design of novel inhibitors of the protease to overcome the problem of drug resistance.     

Structure refinement and diffuse streak scattering of silk (Bombyx mori).

Reexamination of the crystal structure of silk (Bombyx mori) was carried out by X-ray diffraction method. Four molecular chains are contained in the rectangular unit cell with parameters, a = 9.38 A, b = 9.49 A, and c (fiber axis) = 6.98 A, and the space group P2(1)-C(2)2. Silk assumes the statistical crystal structure, in which two antipolar-antiparallel sheet structures with different orientations statistically occupy a crystal site with the ratio 2:1. The molecular conformation is essentially the same pleated sheet structure as Marsh, Corey and Pauling proposed. However, the sheet structure formed by hydrogen bonds assumes the antipolar antiparallel structure different from that proposed by Marsh, Corey and Pauling, in which the methyl groups of alanine residues alternately point to both sides of the sheet structure along the hydrogen bonding direction. The crystalline region of silk is composed of stacking of two antipolar antiparallel sheet structures with different orientations.     

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.