Use of steady state kinetic methods to elucidate the kinetic and chemical mechanisms of retroviral proteases

Transition states of native and drug-resistant HIV-1 protease are the same

HIV-1 protease is an important target for the treatment of HIV/AIDS. However, drug resistance is a persistent problem and new inhibitors are needed. An approach toward understanding enzyme chemistry, the basis of drug resistance, and the design of powerful inhibitors is to establish the structure of enzymatic transition states. Enzymatic transition structures can be established by matching experimental kinetic isotope effects (KIEs) with theoretical predictions. However, the HIV-1 protease transition state has not been previously resolved using these methods. We have measured primary (14)C and (15)N KIEs and secondary (3)H and (18)O KIEs for native and multidrug-resistant HIV-1 protease (I84V). We observed (14)C KIEs ((14)V/K) of 1.029 ± 0.003 and 1.025 ± 0.005, (15)N KIEs ((15)V/K) of 0.987 ± 0.004 and 0.989 ± 0.003, (18)O KIEs ((18)V/K) of 0.999 ± 0.003 and 0.993 ± 0.003, and (3)H KIEs ((3)V/K) KIEs of 0.968 ± 0.001 and 0.976 ± 0.001 for the native and I84V enzyme, respectively. The chemical reaction involves nucleophilic water attack at the carbonyl carbon, proton transfer to the amide nitrogen leaving group, and C-N bond cleavage. A transition structure consistent with the KIE values involves proton transfer from the active site Asp-125 (1.32 Å) with partial hydrogen bond formation to the accepting nitrogen (1.20 Å) and partial bond loss from the carbonyl carbon to the amide leaving group (1.52 Å). The KIEs measured for the native and I84V enzyme indicate nearly identical transition states, implying that a true transition-state analogue should be effective against both enzymes.

Validating the vitality strategy for fighting drug resistance

The current challenge in designing effective drugs against HIV-1 is to find novel candidates with high potency, but with a lower susceptibility to mutations associated with drug resistance. Trying to address this challenge, we developed in our previous study (Ishikita and Warshel, Angew Chem Int Ed Engl 2008; 47:697-700) a novel computational strategy for fighting drug resistance by predicting the likely moves of the virus through constraints on binding and catalysis. This has been based on calculating the ratio between the vitality values ((K(i) k(cat)/K(M))(mutant)/(K(i) k(cat)/K(M))(wild-type)) and using it as a guide for predicting the moves of the virus. The corresponding calculations of the binding affinity, K(i), were carried out using the semi-macroscopic version of the protein dipole Langevin dipole (PDLD/S) in its linear response approximation (LRA) in its β version (PDLD/S-LRA/β). We also calculate the proteolytic efficiency, k(cat)/K(M), by evaluating the transition state (TS) binding free energies using the PDLD/S-LRA/β method. Here we provide an extensive validation of our strategy by calculating the vitality of six existing clinical and experimental drug candidates. It is found that the computationally determined vitalities correlate reasonably well with those derived from the corresponding experimental data. This indicates that the calculated vitality may be used to identify mutations that would be most effective for the survival of the virus. Thus, it should be possible to use our approach in screening for mutations that would provide the most effective resistance to any proposed antiviral drug. This ability should be very useful in guiding the design of drug molecules that will lead to the slowest resistance.

Computational insights into the development of novel therapeutic strategies for Alzheimer's disease

BACKGROUND

β-amyloidosis and oxidative stress have been implicated as root causes of Alzheimer's disease (AD). Current potential therapeutic strategies for the treatment of AD include inhibition of amyloid β (Aβ) production, stimulation of Aβ degradation and prevention of Aβ oligomerization. However, efforts in this direction are hindered by the lack of understanding of the biochemical processes occurring at the atomic level in AD.

DISCUSSION

A radically different approach to achieve this goal would be the application of comprehensive theoretical and computational techniques such as molecular dynamics, quantum mechanics, hybrid quantum mechanics/molecular mechanics, bioinformatics and rotational spectroscopy to investigate complex chemical and physical processes in β-amyloidosis and the oxidative stress mechanism.

CONCLUSION

Results obtained from these studies will provide an atomic level understanding of biochemical processes occurring in AD and advance efforts to develop effective therapeutic strategies for this disease.

Substrate shape determines specificity of recognition for HIV-1 protease: analysis of crystal structures of six substrate complexes

The homodimeric HIV-1 protease is the target of some of the most effective antiviral AIDS therapy, as it facilitates viral maturation by cleaving ten asymmetric and nonhomologous sequences in the Gag and Pol polyproteins. Since the specificity of this enzyme is not easily determined from the sequences of these cleavage sites alone, we solved the crystal structures of complexes of an inactive variant (D25N) of HIV-1 protease with six peptides that correspond to the natural substrate cleavage sites. When the protease binds to its substrate and buries nearly 1000 A2 of surface area, the symmetry of the protease is broken, yet most internal hydrogen bonds and waters are conserved. However, no substrate side chain hydrogen bond is conserved. Specificity of HIV-1 protease appears to be determined by an asymmetric shape rather than a particular amino acid sequence.

The herpesvirus proteases as targets for antiviral chemotherapy

Viruses of the family Herpesviridae are responsible for a diverse set of human diseases. The available treatments are largely ineffective, with the exception of a few drugs for treatment of herpes simplex virus (HSV) infections. For several members of this DNA virus family, advances have been made recently in the biochemistry and structural biology of the essential viral protease, revealing common features that may be possible to exploit in the development of a new class of anti-herpesvirus agents. The herpesvirus proteases have been identified as belonging to a unique class of serine protease, with a Ser-His-His catalytic triad. A new, single domain protein fold has been determined by X-ray crystallography for the proteases of at least three different herpesviruses. Also unique for serine proteases, dimerization has been shown to be required for activity of the cytomegalovirus and HSV proteases. The dimerization requirement seriously impacts methods needed for productive, functional analysis and inhibitor discovery. The conserved functional and catalytic properties of the herpesvirus proteases lead to common considerations for this group of proteases in the early phases of inhibitor discovery. In general, classical serine protease inhibitors that react with active site residues do not readily inactivate the herpesvirus proteases. There has been progress however, with activated carbonyls that exploit the selective nucleophilicity of the active site serine. In addition, screening of chemical libraries has yielded novel structures as starting points for drug development. Recent crystal structures of the herpesvirus proteases now allow more direct interpretation of ligand structure-activity relationships. This review first describes basic functional aspects of herpesvirus protease biology and enzymology. Then we discuss inhibitors identified to date and the prospects for their future development.

Real-time measurements of dark substrate catalysis

We have developed a novel procedure to monitor the real-time cleavage of natural unmodified peptides (dark substrates). In the competition-based assay, the initial cleavage rate of a fluorogenic peptide substrate is measured in the presence of a second substrate that is not required to exhibit any optical property change upon cleavage. Using a unique experimental design and steady-state enzyme kinetics for a two-substrate system, we were able to determine both Km and k(cat) values for cleavage of the dark substrate. The method was applied to HIV-1 protease and to the V82F/I84V drug resistant mutant enzyme. Using two different substrates, we showed that the kinetic parameters derived from the competition assay are in good agreement with those determined independently using standard direct assay. This method can be applied to other enzyme systems as long as they have one substrate for which catalysis can be conveniently monitored in real time.

Structural role of the 30's loop in determining the ligand specificity of the human immunodeficiency virus protease

The structural basis of ligand specificity in human immunodeficiency virus (HIV) protease has been investigated by determining the crystal structures of three chimeric HIV proteases complexed with SB203386, a tripeptide analogue inhibitor. The chimeras are constructed by substituting amino acid residues in the HIV type 1 (HIV-1) protease sequence with the corresponding residues from HIV type 2 (HIV-2) in the region spanning residues 31-37 and in the active site cavity. SB203386 is a potent inhibitor of HIV-1 protease (Ki = 18 nM) but has a decreased affinity for HIV-2 protease (Ki = 1280 nM). Crystallographic analysis reveals that substitution of residues 31-37 (30's loop) with those of HIV-2 protease renders the chimera similar to HIV-2 protease in both the inhibitor binding affinity and mode of binding (two inhibitor molecules per protease dimer). However, further substitution of active site residues 47 and 82 has a compensatory effect which restores the HIV-1-like inhibitor binding mode (one inhibitor molecule in the center of the protease active site) and partially restores the affinity. Comparison of the three chimeric protease structures with those of HIV-1 and SIV proteases complexed with the same inhibitor reveals structural changes in the flap regions and the 80's loops, as well as changes in the dimensions of the active site cavity. The study provides structural evidence of the role of the 30's loop in conferring inhibitor specificity in HIV proteases.

Efficient expression and rapid purification of human T-cell leukemia virus type 1 protease

Human T-cell leukemia virus type 1 (HTLV-1) is an oncovirus that is clinically associated with adult T-cell leukemia. We report here the construction of a pET19-based expression clone containing HTLV-1 protease fused to a decahistidine-containing leader peptide. The recombinant protein is efficiently expressed in Escherichia coli, and the fusion protein can be easily purified by affinity chromatography. Active mature protease in yields in excess of 3 mg/liter of culture can then be obtained by a novel two-step refolding and autoprocessing procedure. The purified enzyme exhibited Km and Kcat, values of 0.3 mM and 0.143 sec(-1) at pH 5.3 and was inhibited by pepstatin A.

Rate-determining steps in HIV-1 protease catalysis. The hydrolysis of the most specific substrate

The human immunodeficiency virus type-1 (HIV-1) encodes a protease which is essential for the production of infectious virus. The protease prefers substrates that contain glutamic acid or glutamine at the P2' position. The catalytic role of these residues has been studied by using a highly specific fluorogen substrate, 2-aminobenzoyl-Thr-Ile-Nle-Phe(NO2)-Gln-Arg (substrate QR), and its counterpart (substrate ER) containing Glu in place of Gln. The newly designed substrate ER that contains a pair of charged residues at P2' and P3' sites is the most specific substrate described so far for HIV-1 protease. The specificity rate constant (kcat/Km = 2.1 x 10(7) M-1 s-1) approaches, but does not reach, the diffusion limit. This follows from the appreciable solvent kinetic deuterium isotope effects on the rate constants, indicating that, independent of the salt concentration, the rate-limiting step of the catalysis is a chemical process rather than a physical one. The reaction also has positive entropy of activation. On the other hand, the rate-limiting step for substrate QR changes with increasing salt concentration from a physical to chemical step, while the negative activation entropy becomes positive. The rate increase with substrate ER is 50-fold with respect to substrate QR in the presence of 0.1 M NaCl and diminishes to 3.5-fold at 2.0 M NaCl concentration, as a consequence of a considerable rate increase at high salt concentration with substrate QR but not with substrate ER. The Km value is much lower for the substrate ER (0.8 microM) than for substrate QR (15 microM), indicating a more effective binding for substrate ER at 0.1 M NaCl. Unexpectedly, the strong binding appears to be achieved by the unionized form of Glu in P2', as follows from the remarkably different pH-rate profiles for substrates QR and ER. The effective binding elicited by the glutamic acid may be utilized in designing inhibitors for therapeutic purposes.

Conformational stability and catalytic activity of HIV-1 protease are both enhanced at high salt concentration

The activity of human immunodeficiency virus protease is markedly increased at elevated salt concentration. The structural basis of this effect has been explored by several independent methods by using both the wild-type enzyme and its triple mutant (Q7K/L33I/L63I) (Mildner, A. M., Rothrock, D. J., Leone, J. W., Bannow, C. A., Lull, J. M., Reardon, I. M., Sarcich, J. L., Howe, W. J., Tomich, C.-S. C., Smith, C. W., Heinrikson, R. L., and Tomasselli, A. G. (1994) Biochemistry 33, 9405-9413), designed to better resist autolysis. Monitoring the intrinsic fluorescence of the two enzymes during urea-mediated denaturation has shown that at high NaCl concentration, both the conformational stability ( DeltaG0) and the transition midpoint (D1/2) between the folded and unfolded states increase, indicating that the salt stabilizes the enzyme structure. These equilibrium data are supported by kinetic studies on the urea-mediated unfolding by measuring fluorescence change, red shifting in the maximum of the emission spectrum, and far- and near-UV CD. The salt effects observed in urea-mediated unfolding reactions prevail upon heat denaturation. All these findings support the existence of a two-state equilibrium between the folded and unfolded proteins. The pH dependence of fluorescence intensity indicated that the conformation of human immunodeficiency virus type 1 protease should change in the catalytically competent pH region. It is concluded that preferential hydration stabilizes the protease structure in the presence of salt, providing entropic contribution to enhance the catalytic activity.