Inhibitor binding induces active site stabilization of the HCV NS3 protein serine protease domain

Hepatitis C: The Story of a Long Journey through First, Second, and Third Generation NS3/4A Peptidomimetic Inhibitors. What Did We Learn?

Hepatitis C viral (HCV) infection is the leading cause of liver failure and still represents a global health burden. Over the past decade, great advancements made HCV curable, and sustained viral remission significantly improved to more than 98%. Historical treatment with pegylated interferon alpha and ribavirin has been displaced by combinations of direct-acting antivirals. These regimens include drugs targeting different stages of the HCV life cycle. However, the emergence of viral resistance remains a big concern. The design of peptidomimetic inhibitors (PIs) able to fit and fill the conserved substrate envelope region within the active site helped avoid contact with the vulnerable sites of the most common resistance-associated substitutions Arg155, Ala156, and Asp168. Herein, we give an overview of HCV NS3 PIs discovered during the past decade, and we deeply discuss the rationale behind the structural optimization efforts essential to achieve pangenotypic activity.

Piperazine-derived small molecules as potential Flaviviridae NS3 protease inhibitors. In vitro antiviral activity evaluation against Zika and Dengue viruses

Since 2011 Direct Acting antivirals (DAAs) drugs targeting different non-structural (NS) viral proteins (NS3, NS5A or NS5B inhibitors) have been approved for clinical use in HCV therapies. However, currently there are not licensed therapeutics to treat Flavivirus infections and the only licensed DENV vaccine, Dengvaxia, is restricted to patients with preexisting DENV immunity. Similarly to NS5 polymerase, the NS3 catalytic region is evolutionarily conserved among the Flaviviridae family sharing strong structural similarity with other proteases belonging to this family and therefore is an attractive target for the development of pan-flavivirus therapeutics. In this work we present a library of 34 piperazine-derived small molecules as potential Flaviviridae NS3 protease inhibitors. The library was developed through a privileged structures-based design and then biologically screened using a live virus phenotypic assay to determine the half-maximal inhibitor concentration (IC₅₀) of each compound against ZIKV and DENV. Two lead compounds, 42 and 44, with promising broad-spectrum activity against ZIKV (IC₅₀ 6.6 µM and 1.9 µM respectively) and DENV (IC₅₀ 6.7 µM and 1.4 µM respectively) and a good security profile were identified. Besides, molecular docking calculations were performed to provide insights about key interactions with residues in NS3 proteases' active sites.

Discovery of small-molecule protein stabilizers toward exogenous alpha-l-iduronidase to reduce the accumulated heparan sulfate in mucopolysaccharidosis type I cells

Synthesis of a series of l-iduronic acid (IdoA)- and imino-IdoA-typed C-glycosides for modulating α-l-iduronidase (IDUA) activity is described. In an enzyme inhibition study, IdoA-typed C-glycosides were more potent than imino-IdoA analogs, with the most potent IdoA-typed C-glycoside 27c showing an IC₅₀ value of 1 μM. On the other hand, co-treatment of 12 with rh-α-IDUA in mucopolysaccharidosis type I (MPS I) fibroblasts exhibited a nearly 3-fold increase of the IDUA activity, resulting in a clear reduction of the accumulated heparan sulfate (HS) compared to the exogenous enzyme treatment alone. This is the first report of small molecules facilitating IDUA stabilization, enhancing enzyme activity, and reducing accumulated HS in MPS I cell-based assays, which reveals that small molecules as rh-α-IDUA stabilizers to improve enzyme replacement therapy (ERT) efficacy toward MPS I is feasible and promising.

Probing the Dynamic Structure-Function and Structure-Free Energy Relationships of the Coronavirus Main Protease with Biodynamics Theory

The SARS-CoV-2 main protease (M^pro) is of major interest as an antiviral drug target. Structure-based virtual screening efforts, fueled by a growing list of apo and inhibitor-bound SARS-CoV/CoV-2 M^pro crystal structures, are underway in many laboratories. However, little is known about the dynamic enzyme mechanism, which is needed to inform both assay development and structure-based inhibitor design. Here, we apply biodynamics theory to characterize the structural dynamics of substrate-induced M^pro activation under nonequilibrium conditions. The catalytic cycle is governed by concerted dynamic structural rearrangements of domain 3 and the m-shaped loop (residues 132–147) on which Cys145 (comprising the thiolate nucleophile and half of the oxyanion hole) and Gly143 (comprising the second half of the oxyanion hole) reside. In particular, we observed the following: (1) Domain 3 undergoes dynamic rigid-body rotation about the domain 2–3 linker, alternately visiting two primary conformational states (denoted as M₁^pro ↔ M₂^pro); (2) The Gly143-containing crest of the m-shaped loop undergoes up and down translations caused by conformational changes within the rising stem of the loop (Lys137–Asn142) in response to domain 3 rotation and dimerization (denoted as M_1/down^pro ↔ 2·M_2/up^pro) (noting that the Cys145-containing crest is fixed in the up position). We propose that substrates associate to the M_1/down^pro state, which promotes the M_2/down^pro state, dimerization (denoted as 2·M_2/up^pro–substrate), and catalysis. Here, we explore the state transitions of M^pro under nonequilibrium conditions, the mechanisms by which they are powered, and the implications thereof for efficacious inhibition under in vivo conditions.

Probing contacts of inhibitor locked in transition states in the catalytic triad of DENV2 type serine protease and its mutants by 1H, 19F and 15 N NMR spectroscopy

Background

Detailed structural knowledge of enzyme-inhibitor complexes trapped in intermediate state is the key for a fundamental understanding of reaction mechanisms taking place in enzymes and is indispensable as a structure-guided drug design tool. Solution state NMR uniquely allows the study of active sites of enzymes in equilibrium between different tautomeric forms. In this study 1H, 19F and 15 N NMR spectroscopy has been used to probe the interaction contacts of inhibitors locked in transition states of the catalytic triad of a serine protease. It was demonstrated on the serotype II Dengue virus NS2B:NS3pro serine protease and its mutants, H51N and S135A, in complex with high-affinity ligands containing trifluoromethyl ketone (tfk) and boronic groups in the C-terminal of tetra-peptides.

Results

Monitoring 19F resonances, shows that only one of the two isomers of the tfk tetra-peptide binds with NS2B:NS3pro and that access to the bulk of the active site is limited. Moreover, there were no bound water found in proximity of the active site for any of the ligands manifesting in a favorable condition for formation of low barrier hydrogen bonds (LBHB) in the catalytic triad. Based on this data we were able to identify a locked conformation of the protein active site. The data also indicates that the different parts of the binding site most likely act independently of each other.

Conclusions

Our reported findings increases the knowledge of the detailed function of the catalytic triad in serine proteases and could facilitate the development of rational structure based inhibitors that can selectively target the NS3 protease of Dengue type II (DENV2) virus. In addition the results shows the usefulness of probing active sites using ¹⁹F NMR spectroscopy.

Structure-thermodynamics-relationships of hepatitis C viral NS3 protease inhibitors

Thermodynamic parameters were determined for structurally-related inhibitors of HCV NS3 protease to assess how binding entropies and enthalpies vary with incremental changes at the P2 and P3 inhibitor subsites. Changing the heterocyclic substituent at P2 from a pyridyl to a 7-methoxy-2-phenyl-4-quinolyl group leads to a 710-fold increase in affinity. Annelating a benzene ring onto a pyridine ring leads to quinoline-derived inhibitors having higher affinities, but the individual enthalpy and entropy contributions are markedly different for each ligand pair. Introducing a phenyl group at C2 of the heterocyclic ring at P2 uniformly leads to higher affinity analogs with more favorable binding entropies, while adding a methoxy group at C7 of the quinoline ring at P2 provides derivatives with more favorable binding enthalpies. Significant enthalpy/entropy compensation is observed for structural changes made to inhibitors lacking a 2-phenyl substituent, whereas favorable changes in both binding enthalpies and entropies accompany structural modifications when a 2-phenyl group is present. Overall, binding energetics of inhibitors having a 2-phenyl-4-quinolyl group at P2 are dominated by entropic effects, whereas binding of the corresponding norphenyl analogs are primarily enthalpy driven. Notably, the reversal from an entropy driven association to an enthalpy driven one for this set of inhibitors also correlates with alternate binding modes. When the steric bulk of the side chain at P3 is increased from a hydrogen atom to a tert-butyl group, there is a 770-fold improvement in affinity. The 30-fold increase resulting from the first methyl group is solely the consequence of a more favorable change in entropy, whereas subsequent additions of methyl groups leads to modest increases in affinity that arise primarily from incremental improvements in binding enthalpies accompanied with smaller favorable entropic contributions.

Synthesis of (3S,4S,5S)-trihydroxylpiperidine derivatives as enzyme stabilizers to improve therapeutic enzyme activity in Fabry patient cell lines

A series of 3S,4S,5S-trihydroxylated piperidines bearing structural diversity at C-2 or C-6 positions has been synthesized and tested to determine their ability to stabilize the activity of recombinant human α-Galactosidase A (rh-α-Gal A). Hit molecules were identified by rapid inhibitory activity screening, and then further investigated for their ability to protect this enzyme from thermo-induced denaturation and enhance its activity in Fabry patient cell lines. Our study resulted in the identification of a new class of small molecules as enzyme stabilizers for the potential treatment of Fabry disease. Of these, stabilizer 21 was the most effective, showing a 12-fold increase in rh-α-Gal A activity in Fabry disease cell lines.

A combinatorial approach towards the synthesis of non-hydrolysable triazole–iduronic acid hybrid inhibitors of human α-l-iduronidase: discovery of enzyme stabilizers for the potential treatment of MPSI

Synthesis and bioevaluation of substituent-diverse triazole–iduronic acid hybrid molecules are highlighted.

The NS4A Cofactor Dependent Enhancement of HCV NS3 Protease Activity Correlates with a 4D Geometrical Measure of the Catalytic Triad Region

We are developing a 4D computational methodology, based on 3D structure modeling and molecular dynamics simulation, to analyze the active site of HCV NS3 proteases, in relation to their catalytic activity. In our previous work, the 4D analyses of the interactions between the catalytic triad residues (His57, Asp81, and Ser139) yielded divergent, gradual and genotype-dependent, 4D conformational instability measures, which strongly correlate with the known disparate catalytic activities among genotypes. Here, the correlation of our 4D geometrical measure is extended to intra-genotypic alterations in NS3 protease activity, due to sequence variations in the NS4A activating cofactor. The correlation between the 4D measure and the enzymatic activity is qualitatively evident, which further validates our methodology, leading to the development of an accurate quantitative metric to predict protease activity in silico. The results suggest plausible “communication” pathways for conformational propagation from the activation subunit (the NS4A cofactor binding site) to the catalytic subunit (the catalytic triad). The results also strongly suggest that the well-sampled (via convergence quantification) structural dynamics are more connected to the divergent catalytic activity observed in HCV NS3 proteases than to rigid structures. The method could also be applicable to predict patients’ responses to interferon therapy and better understand the innate interferon activation pathway.

Structural modeling of HCV NS3/4A serine protease drug-resistance mutations using end-point continuum solvation and side-chain flexibility calculations

Computational methods of modeling protein-ligand interactions have gained widespread application in modern drug discovery. In continuum solvation-based methods of binding affinity estimation, limited description of solvent environment and protein flexibility is traded for a time scale that fits medicinal chemistry test cycles. The results of this speed-accuracy trade-off have been promising in terms of modeling structure-activity relationships of ligand series against protein targets. The potential of these approaches in recapitulating structural and energetic effects of resistance mutations, which involve large changes in binding affinity, remains relatively unexplored. We used continuum solvation binding affinity predictions and graph theory-based flexibility calculations to model thirteen drug resistance mutations in HCV NS3/4A serine protease, against three small-molecule inhibitors, with a 2-fold objective: quantitative assessment of binding energy predictions against experimental data and elucidation of structural/energetic determinants of resistance. The results show statistically significant correlation between predicted and experimental binding affinities, with R(2) and predictive index of up to 0.83 and 0.91, respectively. The level of accuracy was consistent with what has been reported for the inverse problem of binding affinity estimation of congeneric ligands against the same target. The quality of predictions was poor for mutations involving induced-fit effects, primarily, because of the lack of entropy terms. Flexibility analysis explained this discrepancy by indicating characteristic changes in side-chain mobility of a key binding site residue. The combined results from two approaches provide novel insights regarding the molecular mechanism of resistance. NS3/4A inhibitors, with large P2 substituents, derive high affinity with optimal van der Waals interactions in the S2 subsite, in order to overcome unfavorable desolvation and entropic cost of induced-fit effects. High-level resistance mutations tend to increase the desolvation and/or entropic barrier to ligand binding. The lead optimization strategies should, therefore, address the balance of these opposing energetic contributions in both the wild-type and mutant target.

Immunoreactivity assessment of hepatitis C virus NS3 protease and NS5A proteins expressed in TOPO cloning system

BACKGROUND

Hepatitis C virus (HCV) is a major cause of acute and chronic liver disease. Numerous screening assays based on the detection of immunoresponses to HCV structural and nonstructural proteins have been designed. Various studies have demonstrated genotype-specific differences in anti-HCV antibody responses to different HCV proteins.

METHODS

Full-length NS3 protease and N-terminally truncated NS5A were expressed using pET TOPO 102/D system. Antigenicity of the purified recombinant proteins was assessed by immunoblotting and indirect enzyme-linked immunosorbent assay (ELISA). Furthermore, anti-HCV antibody responses to the recombinant proteins were evaluated in three prevalent genotypes in Iran.

RESULTS

We were able to express and purify NS5A and NS3 protease using TOPO cloning system. The HCV NS3 protease and NS5A produced in BL21 Star (DE3) was immunoreactive. Our results demonstrate that NS3 protease and NS5A have good immunoreactivity, but they are not sufficient for detecting all HCV-positive sera. No significant genotype-specific differences were detected in immunoresponses to the recombinant proteins.

CONCLUSION

In conclusion, we successfully isolated, expressed, and purified substantial amount of HCV NS3 protease and N-terminally truncated NS5A, and used them as capturing antigens in a screening ELISA assay with high sensitivity, reproducibility, and specificity. Accordingly, it is well confirmed that TOPO cloning system can be used as a dynamic system in order to express higher amount of immunoreactive viral proteins.

A new method for induced fit docking (GENIUS) and its application to virtual screening of novel HCV NS3-4A protease inhibitors

Hepatitis C virus (HCV) is an etiologic agent of chronic liver disease, and approximately 170 million people worldwide are infected with the virus. HCV NS3-4A serine protease is essential for the replication of this virus, and thus has been investigated as an attractive target for anti-HCV drugs. In this study, we developed our new induced-fit docking program (genius), and applied it to the discovery of a new class of NS3-4A protease inhibitors (IC(50)=1-10 μM including high selectivity index). The new inhibitors thus identified were modified, based on the docking models, and revealed preliminary structure-activity relationships. Moreover, the genius in silico screening performance was validated by using an enrichment factor. We believe our designed scaffold could contribute to the improvement of HCV chemotherapy.

Combined X-ray, NMR, and kinetic analyses reveal uncommon binding characteristics of the hepatitis C virus NS3-NS4A protease inhibitor BI 201335

Hepatitis C virus infection, a major cause of liver disease worldwide, is curable, but currently approved therapies have suboptimal efficacy. Supplementing these therapies with direct-acting antiviral agents has the potential to considerably improve treatment prospects for hepatitis C virus-infected patients. The critical role played by the viral NS3 protease makes it an attractive target, and despite its shallow, solvent-exposed active site, several potent NS3 protease inhibitors are currently in the clinic. BI 201335, which is progressing through Phase IIb trials, contains a unique C-terminal carboxylic acid that binds noncovalently to the active site and a bromo-quinoline substitution on its proline residue that provides significant potency. In this work we have used stopped flow kinetics, x-ray crystallography, and NMR to characterize these distinctive features. Key findings include: slow association and dissociation rates within a single-step binding mechanism; the critical involvement of water molecules in acid binding; and protein side chain rearrangements, a bromine-oxygen halogen bond, and profound pK(a) changes within the catalytic triad associated with binding of the bromo-quinoline moiety.

A comparative QM/MM study of the reaction mechanism of the Hepatitis C virus NS3/NS4A protease with the three main natural substrates NS5A/5B, NS4B/5A and NS4A/4B

The reaction mechanism of the NS3/NS4A protease with the NS4B/5A and NS4A/4B natural substrates has been investigated using the QM/MM (quantum mechanics/molecular mechanics) approach, and some calculations have been performed on the reaction with the NS5A/5B natural substrate. This study widely extends a previous contribution of our group on the reaction mechanism with the NS5A/5B substrate, the main goal here being to understand the differences found between the reaction mechanism of each natural substrate and the role played by the enzymatic residues in the catalytic cycle. This knowledge will ultimately help in developing new NS3/NS4A protease inhibitors. The two first steps of the mechanism have been considered: Acylation and breaking of the peptide bond, with emphasis on the former one (rate limiting process). Energy and free energy profiles for both steps have been calculated at the AM1/MM level and corrected by means of MP2 ab initio calculations, being evident the importance of correlation energy. Acylation is the rate limiting step in all cases and occurs through a tetracoordinated intermediate, as previously suggested for other serine proteases. Specificities in the NS4B/5A mechanism can be attributed to the presence of a Proline residue in the substrate P2 position. The analysis of structures and energies confirm the importance of the oxyanion hole in the electrostatic stabilization of the tetracoordinated intermediate. Finally, the role of other residues, e.g., Arg-155 and Asp-79, has been explained, and the viability of Arg-155 mutants and its resistance to some protease inhibitors has been understood thanks to virtual mutation studies.

Quantitative expression of protein heterogeneity: Response of amino acid side chains to their local environment

A general method has been developed to characterize the hydrophobicity or hydrophilicity of the microenvironment (MENV), in which a given amino acid side chain is immersed, by calculating a quantitative property descriptor (QPD) based on the relative (to water) hydrophobicity of the MENV. Values of the QPD were calculated for a test set of 733 proteins to analyze the modulating effects on amino acid residue properties by the MENV in which they are imbedded. The QPD values and solvent accessibility were used to derive a partitioning of residues based on the MENV hydrophobicities. From this partitioning, a new hydrophobicity scale was developed, entirely in the context of protein structure, where amino acid residues are immersed in one or more "MENVpockets." Thus, the partitioning is based on the residues "sampling" a large number of "solvents" (MENVs) that represent a very large range of hydrophobicity values. It was found that the hydrophobicity of around 80% of amino acid side chains and their MENV are complementary to each other, but for about 20%, the MENV and their imbedded residue can be considered as mismatched. Many of these mismatches could be rationalized in terms of the structural stability of the protein and/or the involvement of the imbedded residue in function. The analysis also indicated a remarkable conservation of local environments around highly conserved active site residues that have similar functions across protein families, but where members have relatively low sequence homology. Thus, quantitative evaluation of this QPD is suggested, here, as a tool for structure-function prediction, analysis, and parameter development for the calculation of properties in proteins.

Novel potent macrocyclic inhibitors of the hepatitis C virus NS3 protease: Use of cyclopentane and cyclopentene P2-motifs

Several highly potent novel HCV NS3 protease inhibitors have been developed from two inhibitor series containing either a P2 trisubstituted macrocyclic cyclopentane- or a P2 cyclopentene dicarboxylic acid moiety as surrogates for the widely used N-acyl-(4R)-hydroxyproline in the P2 position. These inhibitors were optimized for anti HCV activities through examination of different ring sizes in the macrocyclic systems and further by exploring the effect of P4 substituent removal on potency. The target molecules were synthesized from readily available starting materials, furnishing the inhibitor compounds in good overall yields. It was found that the 14-membered ring system was the most potent in these two series and that the corresponding 13-, 15-, and 16-membered macrocyclic rings delivered less potent inhibitors. Moreover, the corresponding P1 acylsulfonamides had superior potencies over the corresponding P1 carboxylic acids. It is noteworthy that it has been possible to develop highly potent HCV protease inhibitors that altogether lack the P4 substituent. Thus the most potent inhibitor described in this work, inhibitor 20, displays a K(i) value of 0.41 nM and an EC(50) value of 9 nM in the subgenomic HCV replicon cell model on genotype 1b. To the best of our knowledge this is the first example described in the literature of a HCV protease inhibitor displaying high potency in the replicon assay and lacking the P4 substituent, a finding which should facilitate the development of orally active small molecule inhibitors against the HCV protease.

A quantum mechanics/molecular mechanics study of the reaction mechanism of the hepatitis C virus NS3 protease with the NS5A/5B substrate

Combined quantum mechanics and molecular mechanics (QM/MM) calculations were carried out to characterize the reaction mechanism of the NS3 protease with its preferred substrate (NS5A/5B). The main purpose of this study was to locate the barrier states and intermediates along the distinguished coordinate path (DCP) involved in this process. These structures, and in particular the one corresponding to the first barrier state and intermediate (B1 and I1), could be a starting point for the synthesis of inhibitors of this protease, which could be used to treat hepatitis C. The two first steps of the reaction mechanism were studied, i.e., the acylation step and the breaking of the peptide bond. The first step takes place through a tetracoordinated intermediate, as suggested from previous works on other Serine proteases. The importance of the different amino acid residues was also considered (perturbation study where the MM charges of each residue were set to zero independently). The residues of the oxyanion hole were confirmed as the most important for the electrostatic stabilization of the tetracoordinate intermediate. Moreover, the role of other residues, e.g., Arg-155 and Asp-79, was also explained.

Synthesis of novel potent hepatitis C virus NS3 protease inhibitors: Discovery of 4-hydroxy-cyclopent-2-ene-1,2-dicarboxylic acid as a N-acyl-l-hydroxyproline bioisostere

Potent tetrapeptidic inhibitors of the HCV NS3 protease have been developed incorporating 4-hydroxy-cyclopent-2-ene-1,2-dicarboxylic acid as a new N-acyl-l-hydroxyproline mimic. The hydroxycyclopentene template was synthesized in eight steps from commercially available (syn)-tetrahydrophthalic anhydride. Three different amino acids were explored in the P1-position and in the P2-position the hydroxyl group of the cyclopentene template was substituted with 7-methoxy-2-phenyl-quinolin-4-ol. The P3/P4-positions were then optimized from a set of six amino acid derivatives. All inhibitors were evaluated in an in vitro assay using the full-length NS3 protease. Several potent inhibitors were identified, the most promising exhibiting a K(i) value of 1.1nM.

Homology modeling and molecular dynamics study of West Nile virus NS3 protease: A molecular basis for the catalytic activity increased by the NS2B cofactor

The West Nile virus (WNV) NS3 serine protease, which plays an important role in assembly of infective virion, is an attractive target for anti-WNV drug development. Cofactors NS2B and NS4A increase the catalytic activity of NS3 in dengue virus and Hepatitis C virus, respectively. Recent studies on the WNV-NS3 characterize the catalytically active form of NS3 by tethering the 40-residue cofactor NS2B. It is suggested that NS2B is essential for the NS3 activity in WNV, while there is no information of the WNV-NS3-related crystal structure. To understand the role of NS2B/substrate in the NS3 catalytic activity, we built a series of models: WNV-NS3 and WNV-NS3-NS2B and WNV-NS3-NS2B-substrate using homology modeling and molecular modeling techniques. Molecular dynamics (MD) simulations were performed for 2.75 ns on each model, to investigate the structural stabilization and catalytic triad motion of the WNV NS3 protease with and without NS2B/substrate. The simulations show that the NS3 rearrangement occurs upon the NS2B binding, resulting in the stable D75-OD1...H51-NH hydrogen bonding. After the substrate binds to the NS3-NS2B active site, the NS3 protease becomes more stable, and the catalytic triad is formed. These results provide a structural basis for the activation and stabilization of the enzyme by its cofactor and substrate.

Potent inhibitors of the hepatitis C virus NS3 protease: Use of a novel P2 cyclopentane-derived template

The HCV NS3 protease is essential for replication of the hepatitis C virus (HCV) and therefore constitutes a promising new drug target for anti-HCV therapy. Several potent and promising HCV NS3 protease inhibitors, some of which display low nanomolar activities, were identified from a series of novel inhibitors incorporating a trisubstituted cyclopentane dicarboxylic acid moiety as a surrogate for the widely used N-acyl-(4R)-hydroxyproline in the P2 position.

Calpain Inhibition by α-Ketoamide and Cyclic Hemiacetal Inhibitors Revealed by X-ray Crystallography,

Calpains are intracellular calcium-activated cysteine proteases whose unregulated proteolysis following the loss of calcium homeostasis can lead to acute degeneration during ischemic episodes and trauma, as well as Alzheimer's disease and cataract formation. The determination of the crystal structure of the proteolytic core of mu-calpain (muI-II) in a calcium-bound active conformation has made structure-guided design of active site inhibitors feasible. We present here high-resolution crystal structures of rat muI-II complexed with two reversible calpain-specific inhibitors employing cyclic hemiacetal (SNJ-1715) and alpha-ketoamide (SNJ-1945) chemistries that reveal new details about the interactions of inhibitors with this enzyme. The SNJ-1715 complex confirms that the free aldehyde is the reactive species of the cornea-permeable cyclic hemiacetal. The alpha-ketoamide warhead of SNJ-1945 binds with the hydroxyl group of the tetrahedral adduct pointing toward the catalytic histidine rather than the oxyanion hole. The muI-II-SNJ-1945 complex shows residue Glu261 displaced from the S1' site by the inhibitor, resulting in an extended "open" conformation of the domain II gating loop and an unobstructed S1' site. This conformation offers an additional template for structure-based drug design extending to the primed subsites. An important role for the highly conserved Glu261 is proposed.

Structure and dynamics of coxsackievirus B4 2A proteinase, an enyzme involved in the etiology of heart disease

The 2A proteinases (2A(pro)) from the picornavirus family are multifunctional cysteine proteinases that perform essential roles during viral replication, involving viral polyprotein self-processing and shutting down host cell protein synthesis through cleavage of the eukaryotic initiation factor 4G (eIF4G) proteins. Coxsackievirus B4 (CVB4) 2A(pro) also cleaves heart muscle dystrophin, leading to cytoskeletal dysfunction and the symptoms of human acquired dilated cardiomyopathy. We have determined the solution structure of CVB4 2A(pro) (extending in an N-terminal direction to include the C-terminal eight residues of CVB4 VP1, which completes the VP1-2A(pro) substrate region). In terms of overall fold, it is similar to the crystal structure of the mature human rhinovirus serotype 2 (HRV2) 2A(pro), but the relatively low level (40%) of sequence identity leads to a substantially different surface. We show that differences in the cI-to-eI2 loop between HRV2 and CVB4 2A(pro) translate to differences in the mechanism of eIF4GI recognition. Additionally, the nuclear magnetic resonance relaxation properties of CVB4 2A(pro), particularly of residues G1 to S7, F64 to S67, and P107 to G111, reveal that the substrate region is exchanging in and out of a conformation in which it occupies the active site with association and dissociation rates in the range of 100 to 1,000 s(-1). This exchange influences the conformation of the active site and points to a mechanism for how self-processing can occur efficiently while product inhibition is avoided.

Identification of residues in the dengue virus type 2 NS2B cofactor that are critical for NS3 protease activation

Proteolytic processing of the dengue virus polyprotein is mediated by host cell proteases and the virus-encoded NS2B-NS3 two-component protease. The NS3 protease represents an attractive target for the development of antiviral inhibitors. The three-dimensional structure of the NS3 protease domain has been determined, but the structural determinants necessary for activation of the enzyme by the NS2B cofactor have been characterized only to a limited extent. To test a possible functional role of the recently proposed Phix(3)Phi motif in NS3 protease activation, we targeted six residues within the NS2B cofactor by site-specific mutagenesis. Residues Trp62, Ser71, Leu75, Ile77, Thr78, and Ile79 in NS2B were replaced with alanine, and in addition, an L75A/I79A double mutant was generated. The effects of these mutations on the activity of the NS2B(H)-NS3pro protease were analyzed in vitro by sodium dodecyl sulfate-polyacrylamide gel electrophoresis analysis of autoproteolytic cleavage at the NS2B/NS3 site and by assay of the enzyme with the fluorogenic peptide substrate GRR-AMC. Compared to the wild type, the L75A, I77A, and I79A mutants demonstrated inefficient autoproteolysis, whereas in the W62A and the L75A/I79A mutants self-cleavage appeared to be almost completely abolished. With exception of the S71A mutant, which had a k(cat)/K(m) value for the GRR-AMC peptide similar to that of the wild type, all other mutants exhibited drastically reduced k(cat) values. These results indicate a pivotal function of conserved residues Trp62, Leu75, and Ile79 in the NS2B cofactor in the structural activation of the dengue virus NS3 serine protease.

Potent Inhibitors of the Hepatitis C Virus NS3 Protease: Design and Synthesis of Macrocyclic Substrate-Based β-Strand Mimics

The virally encoded NS3 protease is essential to the life cycle of the hepatitis C virus (HCV), an important human pathogen causing chronic hepatitis, cirrhosis of the liver, and hepatocellular carcinoma. The design and synthesis of 15-membered ring beta-strand mimics which are capable of inhibiting the interactions between the HCV NS3 protease enzyme and its polyprotein substrate will be described. The binding interactions between a macrocyclic ligand and the enzyme were explored by NMR and molecular dynamics, and a model of the ligand/enzyme complex was developed.

Structure-based design of inhibitors of NS3 serine protease of hepatitis C virus

We have designed small focused combinatorial library of hexapeptide inhibitors of NS3 serine protease of the hepatitis C virus (HCV) by structure-based molecular design complemented by combinatorial optimisation of the individual residues. Rational residue substitutions were guided by the structure and properties of the binding pockets of the enzyme's active site. The inhibitors were derived from peptides known to inhibit the NS3 serine protease by using unusual amino acids and alpha-ketocysteine or difluoroaminobutyric acid, which are known to bind to the S1 pocket of the catalytic site. Inhibition constants (Ki) of the designed library of inhibitors were predicted from a QSAR model that correlated experimental Ki of known peptidic inhibitors of NS3 with the enthalpies of enzyme-inhibitor interaction computed via molecular mechanics and the solvent effect contribution to the binding affinity derived from the continuum model of solvation. The library of the optimised inhibitors contains promising drug candidates-water-soluble anionic hexapeptides with predicted Ki* in the picomolar range.

Identification of protons position in acid-base enzyme catalyzed reactions: The hepatitis C viral NS3 protease

General acid-base catalysis is a key element of the catalytic activity of most enzymes. Therefore, any explicit molecular modeling of enzyme-catalyzed chemical reactions requires correct identification of protons location on the catalytic groups. In this work, we apply our quantum mechanical/self-consistent reaction field in virtual solvent [QM/SCRF(VS)] method for identification of the position of protons shared by the enzyme catalytic groups and the polar groups of the inhibitor in a covalent tetrahedral complex (TC) of the hepatitis C virus NS3 protease with a peptidyl alpha-ketoacid inhibitor. To identify the relevant protonation states, we have analyzed relative stabilities of R and S configurations of the TC that depend on the specific proton distribution over the polar groups and correlated it with experimental NMR and X-ray crystallography data, both at low and neutral pH ranges. The tentative assignment of the single resonance in the (13)C NMR spectrum of the hemiketal carbon at physiological pH to the S configuration of TC is confirmed. Both R and S configurations are equally stable at acidic pH in our modeling, in good agreement with the (13)C NMR observation.

Application of the lambda-dynamics method to evaluate the relative binding free energies of inhibitors to HCV protease

The lambda-dynamics method was used to calculate the relative binding free energies of inhibitors to the hepatitis C virus (HCV) protease. A total of seven HCV protease p-side product inhibitors were used in this study. The inhibitors are 6-mer peptides spanning P6-P1 (Ac-Asp-d-Glu-Leu-Ile-Cha-P1-CO(2)H). For this protein, S1 is a major hydrophobic pocket for binding. Binding of various residues to this pocket was investigated through free energy simulations and experimental inhibition constants. Several 300 ps lambda-dynamics simulations in explicit solvent were performed. The relative binding free energy was estimated from these simulations. From a single simulation, the inhibitors can be correctly classified into highly potent and weakly potent groups. The multiple simulations give an accurate rank ordering of inhibitor potency; computed and experimental binding free energies agree with 0.6 kcal/mol for five of the seven inhibitors. In addition, free energy perturbation (FEP) calculations were carried out to validate the results from lambda-dynamics. A total of 6 ligand pairs were compared. For each pair, 5-11 windows were used to map one ligand to the other. The cumulative simulation time was over 2 ns for each ligand pair. For four of the six ligand pairs, the lambda-dynamics free energy difference fits better than the FEP difference to the experimental value. The fact that the lambda-dynamics method achieved similar results in only a fraction of the total simulation time for FEP further demonstrates the robustness of the lambda-dynamics method.

The design of a potent inhibitor of the hepatitis C virus NS3 protease:BILN 2061?From the NMR tube to the clinic

The virally encoded serine protease NS3/NS4A is essential to the life cycle of the hepatitis C virus (HCV), an important human pathogen causing chronic hepatitis, cirrhosis of the liver, and hepatocellular carcinoma. Until very recently, the design of inhibitors for the HCV NS3 protease was limited to large peptidomimetic compounds with poor pharmacokinetic properties, making drug discovery an extremely challenging endeavor. In our quest for the discovery of a small-molecule lead that could block replication of the hepatitis C virus by binding to the HCV NS3 protease, the critical protein-polypeptide interactions between the virally encoded NS3 serine protease and its polyprotein substrate were investigated. Lead optimization of a substrate-based hexapeptide, guided by structural data, led to the understanding of the molecular dynamics and electronic effects that modulate the affinity of peptidomimetic ligands for the active site of this enzyme. Macrocyclic beta-strand scaffolds were designed that allowed the discovery of potent, highly selective, and orally bioavailable compounds. These molecules were the first HCV NS3 protease inhibitors reported that inhibit replication of HCV subgenomic RNA in a cell-based replicon assay at low nanomolar concentrations. Optimization of their biopharmaceutical properties led to the discovery of the clinical candidate BILN 2061. Oral administration of BILN 2061 to patients infected with the hepatitis C genotype 1 virus resulted in an impressive reduction of viral RNA levels, establishing proof-of-concept for HCV NS3 protease inhibitors as therapeutic agents in humans.

Serpin mechanism of hepatitis C virus nonstructural 3 (NS3) protease inhibition: induced fit as a mechanism for narrow specificity

Hepatitis C virus (HCV) nonstructural 3 (NS3) serine protease disrupts important cellular antiviral signaling pathways and plays a pivotal role in the proteolytic maturation of the HCV polyprotein precursor. This recent discovery has fostered the search for NS3 protease inhibitors. However, the enzyme's unusual induced fit behavior and peculiar molecular architecture have imposed considerable obstacles to the development of small molecule inhibitors. In this article, we demonstrate that such unique induced fit behavior and the chymotrypsin-like catalytic domain can provide the structural plasticity necessary to generate protein-based inhibitors of the NS3 protease. We took advantage of the macromolecular scaffold of a Drosophila serpin, SP6, which intrinsically supports chymotrypsin-like enzyme inhibition, to design a novel class of potent and selective inhibitors. We show that altering the SP6 reactive site loop (RSL) resulted in the development of the first effective (K(i) of 34 nm) and selective serpin, SP6(EVC/S), directed at the NS3 protease. SP6(EVC/S) operates as a suicide substrate inhibitor, and its partitioning between the complex-forming and proteolytic pathways for the NS3 protease is HCV NS4A cofactor-dependent and -specific. Once bound to the protease active site, SP6(EVC/S) partitions with equal probability to undergo proteolysis by NS3 at the C-terminal site of the engineered RSL, (P(6))Glu-Ile-(P(4))Val-Met-Thr-(P(1))Cys- downward arrow -(P(1)')Ser, or to form a covalent acyl-enzyme complex characteristic of cognate protease-serpin pairs. Our results also reveal a novel cofactor-induced serpin mechanism of enzyme inhibition that could be explored for developing effective and selective inhibitors of other important induced fit viral proteases of the Flaviviridae family such as the West Nile virus NS3 endoprotease.

Structure-based design of hepatitis C virus inhibitors

The limitations of current treatment for chronic hepatitis C virus (HCV) infection have prompted the development of novel therapeutic strategies targeting events specific to viral replication. Over the past decade, advances in the study of HCV molecular biology have led to the identification of cis-acting RNA sequences and viral enzymatic activities which present attractive targets for inhibition. High-resolution, three-dimensional structures of the HCV serine protease, helicase and RNA-dependent RNA polymerase have been determined through X-ray crystallographic studies. More recently, solution structures of these proteins and the HCV internal ribosome entry site have been evaluated by nuclear magnetic resonance spectroscopy and electron microscopy. Mutational analysis and structural characterization of these macromolecules in complex with bound substrates, cofactors and inhibitors has further defined the various electrochemical interactions which mediate protein-protein, protein-RNA and other intermolecular contacts. This review will discuss the available structural data with respect to the rational design of HCV enzyme inhibitors and the development of antisense-based therapeutic strategies, such as RNA interference.

Mimicry between the hepatitis C virus polyprotein and antigenic targets of nuclear and smooth muscle antibodies in chronic hepatitis C virus infection

Autoantibodies to smooth muscle (SMA) and nuclear components (ANA) arise in the natural course of chronic infection with hepatitis C virus. In view of the growing evidence for 'molecular mimicry' as a mechanism of autoimmunity we investigated whether cross-reactive immune reactions between host smooth muscle/nuclear components and HCV antigens may contribute to the formation of SMA and ANA in chronic HCV infection. Computer-assisted protein database search methods were used to identify three smooth muscle (smoothelin698-717, myosin1035-1054, vimentin69-88) and three nuclear (matrin722-741, histone H2A11-30, replication protein A133-152) host antigens with the highest local sequence similarity to the HCV polyprotein and 20-mer peptides corresponding to these regions were constructed. Sera from 51 children with chronic HCV infection [median age: 8 (2-16); 27 boys], 26 SMA positive and five ANA positive, were tested for reactivity to the synthesized HCV peptides and their human homologues by enzyme linked immunosorbent assay (ELISA). Sera from patients with HBV infection and chronic liver disease of different aetiologies were used as controls. 'Double reactivity' to HCV peptides and smooth muscle/nuclear homologues was associated strongly with HCV infection (P < 0.001 for both). Humoral cross-reactivity was established as the basis for double recognition by competition ELISA. Double-reactivity to smooth muscle and HCV peptide antigens correlated with SMA positivity by indirect immunofluouresence (P = 0.05). Of 15 patients double-reactive to myosin1035-1054 and its HCV homologue, 13 recognized whole myosin by immunoblot. These results suggest that ANA and SMA in chronic HCV infection may arise, at least in part, as a consequence of cross-reactive immune responses to HCV and host smooth muscle/nuclear antigens.

Novel, potent phenethylamide inhibitors of the hepatitis C virus (HCV) NS3 protease: probing the role of P2 aryloxyprolines with hybrid structures

Synthesis of hybrid HCV NS3 protease/NS4A inhibitors having the 4,4-difluoroaminobutyric acid (difluoroAbu) phenethylamides as P1-P1' and quinolyloxyprolines as P2 fragments led to 7 (IC(50) 54 nM). Molecular modelling suggests that this potent tripeptide inhibitor utilizes interactions in the S1', S1, S2, S3 and S4 sites of the protease.

Inhibiting viral proteases: challenges and opportunities

Inhibitor design against viral targets must take into account the peculiar characteristics of viral biology-in particular, the plasticity of their replicative machinery. This includes maturational cleavage of the polyprotein, which is mediated by virally encoded proteases. Designing against a movable target is particularly challenging, but at the same time it offers new opportunities. Here we describe our experience with the NS3/4A (NS: nonstructural) serine protease of human hepatitis C virus (HCV). By extensive use of combinatorial peptide libraries, various inhibitor types were generated, including product inhibitors, serine traps, P-P' inhibitors, and prime side inhibitors. The latter represent a first case for a serine protease. A key finding, derived from structural studies utilizing these inhibitors, was that NS3 is an induced-fit protease, requiring both the NS4A cofactor protein and the substrate to fully activate its catalytic machinery. In the absence of cofactor and/or substrate, NS3 exists in solution as a large conformational ensemble, which can be matched by a correspondingly large set of peptide inhibitors, each one stabilizing a given conformer. In the perspective of inhibiting viral proteases in general, we suggest that combinatorial ligand ensembles may be a powerful tool, to contrast the adaptive potential of the viral quasispecies.

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.

Inhibition of hepatitis C virus NS3 protease by peptides derived from complementarity-determining regions (CDRs) of the monoclonal antibody 8D4: tolerance of a CDR peptide to conformational changes of a target

We have synthesized and characterized peptides derived from complementarity-determining regions (CDRs) of 8D4, a mouse monoclonal antibody against NS3 protease domain of hepatitis C virus. 8D4 inhibits enzymatic activity without its cofactor, NS4A peptide. One of the synthetic peptides derived from CDRs, CDR1 of the heavy-chain (CDR-H1) peptide strongly inhibited NS3 protease activity competitively in the absence of NS4A and non-competitively in the presence of NS4A. Moreover, cyclic CDR-H1 peptides bridged by disulfide inhibited NS3 protease more potently. The chain length of the CDR-H1 peptide is critical for strong inhibition, even when the peptide is circularized. This finding suggests the importance of peptide conformation. In contrast to a cognate antibody molecule, CDR-derived peptides may provide good ligands for target molecules by having a tolerance to conformational changes of the targets caused by cofactor binding or mutation.

Evolution, synthesis and SAR of tripeptide alpha-ketoacid inhibitors of the hepatitis C virus NS3/NS4A serine protease

N-terminal truncation of the hexapeptide ketoacid 1 gave rise to potent tripeptide inhibitors of the hepatitis C virus NS3 protease/NS4A cofactor complex. Optimization of these tripeptides led to ketoacid 30 with an IC50 of 0.38 microM. The SAR of these tripeptides is discussed in the light of the recently published crystal structures of a ternary tripetide/NS3/NS4A complexes.

Hepatitis C virus NS3 protease requires its NS4A cofactor peptide for optimal binding of a boronic acid inhibitor as shown by NMR

NMR spectroscopy was used to characterize the hepatitis C virus (HCV) NS3 protease in a complex with the 24 residue peptide cofactor from NS4A and a boronic acid inhibitor, Ac-Asp-Glu-Val-Val-Pro-boroAlg-OH. Secondary-structure information, NOE constraints between protease and cofactor, and hydrogen-deuterium exchange rates revealed that the cofactor was an integral strand in the N-terminal beta-sheet of the complex as observed in X-ray crystal structures. Based upon chemical-shift perturbations, inhibitor-protein NOEs, and the protonation state of the catalytic histidine, the boronic acid inhibitor was bound in the substrate binding site as a transition state mimic. In the absence of cofactor, the inhibitor had a lower affinity for the protease. Although the inhibitor binds in the same location, differences were observed at the catalytic site of the protease.

Recent advances in the molecular biology of hepatitis C virus

The Hepatitis C virus is a positive-stranded RNA virus which is the causal agent for a chronic liver infection afflicting more than 170,000,000 people world-wide. The HCV genome is approximately 9.6 kb in length and the proteome encoded is a polyprotein of a little more than 3000 amino acid residues. This polyprotein is processed by a combination of host and viral proteases into structural and non-structural proteins. The functions of most of these proteins have been established by analogy to other viruses and by sequence homology to known proteins, as well as subsequent biochemical analysis. Two of the non-structural proteins, NS4b and NS5a, are still of unknown function. The development of antivirals for this infectious agent has been hampered by the lack of robust and economical cell culture and animal infection systems. Recent progress in the molecular virology of HCV has come about due to the definition of molecular clones, which are infectious in the chimpanzee, the development of a subgenomic replicon system in Huh7 cells, and the description of a transgenic mouse model for HCV infection. Recent progress in the structural biology of the virus has led to the determination of high resolution three-dimensional structures of a number of the key virally encoded enzymes, including the NS3 protease, NS3 helicase, and NS5b RNA-dependent RNA polymerase. In some cases these structures have been determined in complex with substrates, co-factors (NS4a), and inhibitors. Finally, a variety of techniques have been used to define host factors, which may be required for HCV replication, although this work is just beginning.

Assessment of substrate specificity of hepatitis G virus NS3 protease by a genetic method

The RNA genome of hepatitis G virus (HGV) encodes a large polyprotein that is processed to mature proteins by viral-encoded proteases. The HGV NS3 protease is responsible for the cleavage of the HGV polyprotein at four different locations. No conserved sequence motif has been identified for the cleavage sites of the NS3 protease. To determine the substrate specificity of the NS3 protease, amino acid sequences cleaved by the NS3 protease were obtained from randomized sequence libraries by using a screening method referred to as GASP (Genetic Assay for Site-specific Proteolysis). Based on statistical analyses of the obtained cleavable sequences, a consensus substrate sequence was deduced: Gln-Glu-Thr-Leu-Val downward arrow Ser, with the scissile bond located between Val and Ser. The relevance of this peptide as a cleavable substrate was further supported by molecular modeling of the NS3 protease. Our result would provide an insight on the molecular activity of the NS3 protease and may be useful for the design of substrate-based inhibitors.