Interactions of different inhibitors with active-site aspartyl residues of HIV-1 protease and possible relevance to pepsin

Non-active site mutations in the HIV protease: Diminished drug binding affinity is achieved through modulating the hydrophobic sliding mechanism

The global HIV/AIDS epidemic still currently affects approximately 38 million individuals globally. The protease enzyme of the human immunodeficiency virus is a major drug target in antiviral therapy, however, under the influence of reverse transcriptase and in the context of drug pressure, the rapid PR mutation rate contributes significantly to clinical failure. The set of cooperative non-active site mutations, I13V/I62V/V77I, have been associated with reduced inhibitor susceptibility and are the focus of the current study. When compared to the wild-type protease the mutant protease exhibited decreased binding affinities towards ATV and DRV by 64- and 12-fold, respectively, and decreased the overall favourable Gibbs free energy for ATV, DRV, RTV and SQV. Moreover, these mutations decreased the thermal stability of the protease when in complex with ATV and DRV by approximately 6.4 and 4.2 °C, respectively. The crystal structure of the mutant protease revealed that the location of these mutations and their effect on the hydrophobic sliding mechanism may be crucial in their role in resistance.

Exploring the effect of ritonavir and TMC-310911 on SARS-CoV-2 and SARS-CoV main proteases: potential from a molecular perspective

AIM

As coronavirus (CoV) disease 2019-associated pneumonia spreads globally, there has been an urgent need to combat the spread and develop vaccines.

MATERIALS & METHODS

We used an integrated computational algorithm to explore the binding mechanism of TMC-310911/ritonavir (RVT) with SARS-CoV-2 and SARS-CoV main proteases.

RESULTS

RVT and TMC-310911 had favorable interactions with the proteases, and these high interactions are facilitated by some significant residues such as Asn133, Gly195 and Gln192. Our study further implicated two important rings in the structure of RVT as a possible chemical culprit in its therapeutic activity.

CONCLUSION

Although there are conflicting clinical results on the therapeutic potency of RVT in the treatment of coronavirus disease 2019, our findings provided molecular insight into the binding mechanism of TMC-310911 and RVT with SARS-CoV-2 and SARS-CoV main proteases.

Sensing HIV Protease and Its Inhibitor Using "Helical Epitope"-Imprinted Polymers

A helical epitope-peptide (lle⁸⁵-Gly⁹⁴) was selected from the α-helix structure of the HIV protease (PR) as the template, which represents an intricate interplay between structure conformation and dimerization. The peptide template was mixed with water, trifluoroethanol (TFE), and acetonitrile (ACN) at a certain ratio to enlarge the helical conformation in the solution for the fabrication of helical epitope-mediated molecularly imprinted polymers (HEMIPs) on a quartz crystal microbalance (QCM) chip. The template molecules were then removed under equilibrium batch rebinding conditions involving 5% acetic acid/water. The resulting HEMIPs chip exhibited a high affinity toward template peptide HIV PR_85-94, His-tagged HIV PR, and HIV PR, with dissociation constants (K_d) as 160, 43.3, and 78.5 pM, respectively. The detection limit of the developed HIV PR_85-94 QCM sensor is 0.1 ng/mL. The HEMIPs chip exhibited a high affinity and selectivity to bind HIV PR and subsequently to an inhibitor of HIV PR (nelfinavir). The HIV PR binding site was properly oriented on the HEMIPs-chip to develop a HIV PR/HEMIPs chip, which can effectively bind nelfinavir to establish a sandwich assay. The nelfinavir then attached to the HIV PR/HEMIPs chip, which can be easily removed involving 0.8% acetic acid/water. Therefore, HIV PR/HEMIPs chip can be useful to screen for other HIV PR inhibitors. This technique may improve drug targeting for HIV therapy and also strengthen investigations into other virus assays.

Class A G Protein-Coupled Receptor Antagonist Famotidine as a Therapeutic Alternative Against SARS-CoV2: An In Silico Analysis

The pandemic associated with Severe Acute Respiratory Syndrome Coronavirus type 2 (SARS-CoV2) and its disease named COVID-19 challenged the scientific community to discover effective therapeutic solutions in a short period. Repurposing existing drugs is one viable approach that emphasizes speed during these urgent times. Famotidine, a class A G protein-coupled receptor antagonist used for the treatment of gastroesophageal reflux was recently identified in an in silico screening. Additionally, a recent retrospective clinical report showed that the treatment with famotidine provided a good outcome in patients infected with SARS-CoV2. A clinical trial testing effectiveness of famotidine in combination with hydroxychloroquine is currently ongoing in the United States (US). In the 1990s, famotidine was described as an antiviral agent against human immunodeficiency virus (HIV). Interestingly, some HIV protease inhibitors are presently being used against SARS-CoV2. However, it is not clear if famotidine could be effective against SARS-CoV2. Thus, by using a computational analysis, we aimed to examine if the antiviral effect of famotidine could be related to the inhibition of proteases involved in the virus replication. Our results showed that famotidine could interact within the catalytic site of the three proteases associated with SARS-CoV2 replication. However, weak binding affinity of famotidine to these proteases suggests that a successful famotidine therapy could likely be achieved only in combination with other antiviral drugs. Finally, analysis of famotidine's pharmacokinetic parameters indicated that its effect against SARS-CoV2 infection could be reached only upon intravenous administration. This work will contribute to the pharmacological knowledge of famotidine as an antiviral agent against SARS-CoV2.

A dual target of Plasmepsin IX and X: Unveiling the atomistic superiority of a core chemical scaffold in malaria therapy

Plasmepsin IX and X, members of the prominent aspartic family of proteases whose function were hitherto unknown have only recently been established as key mediators of erythrocyte invasion and egress of the virulent malarial parasite. Inhibitor 49c, a potent antimalarial peptidomimetic inhibitor initially developed to target Plasmepsin II has lately been proven to exhibit potent inhibitory activity against Plasmepsin IX and X. However, the molecular and structural dynamics supporting its inhibitory activity remain inconclusive. Hindering the motion of the flap and hinge region of an aspartic protease remains essential for disabling the catalytic activity of the enzyme. Integrating molecular dynamic simulations coupled with other advanced biocomputational tools, we reveal the enhanced structural mechanistic competence of 49c in complex with Plasmepsin IX and X relative to Pepstatin. Pepstatin, a known aspartic protease inhibitor which actively hinders the opening and closing of the flap tip and flexible loop and consequently limits access to the catalytic aspartic residues, however, its administration has been related to elevated levels of toxicity. Thermodynamic calculations reveal a higher relative binding free energy associated with Plasmepsin IX and X in complex with 49c as opposed to Pepstatin. A relatively compact and structurally rigid 49c bound complexes sequel into the restriction of the flap and hinge residues by restraining cohesive movement, consequently hindering their "twisting motion" from transpiring. Findings unveil an atomistic perspective into the structural superiority of 49c in complex with Plasmepsin IX and X.

Exploring the flap dynamics of the South African HIV subtype C protease in presence of FDA-approved inhibitors: MD study

HIV-1 protease (HIV PR) is considered as one of the most attractive targets for the treatment of HIV and the impact of flap dynamics of HIV PR on the binding affinities of protease inhibitors (PIs) is a crucial ongoing research field. Recently, our research group evaluated the binding affinities of different FDA approved PIs against the South African HIV-1 subtype C (C-SA) protease (PR). The CSA-HIV PR displayed weaker binding affinity for most of the clinical PIs compared to HIV-1 B subtype for West and Central Europe, the Americas. In the current work, the flap dynamics of four different systems of HIV-1 C-SA PR complexed to FDA approved second generation PIs and its impact on binding was explored over the molecular dynamic trajectories. It was observed that the interactions of the selected drugs with the binding site residues of the protease may not be the major contributor for affinity towards PIs. Various post-MD analyses were performed, also entropic contributions, solvation free energies and hydrophobic core formation interactions were studied to assess how the flap dynamics of C-SA PR which is affected by such factors. From these contributions, large van der Waals interactions and low solvation free energies were found to be major factors for the higher activity of ATV against C-SA HIV PR. Furthermore, a comparatively stable hydrophobic core may be responsible for higher stability of the PR flaps of the ATV complex. The outcome of this study provides significant guidance to how the flap dynamics of C-SA PR is affected by various factors as a result of the binding affinity of various protease inhibitors. It will also assist with the design of potent inhibitors against C-SA HIV PR that apart from binding in the active site of PR can interacts with the flaps to prevent opening of the flaps resulting in inactivation of the protease.

Evolution under Drug Pressure Remodels the Folding Free-Energy Landscape of Mature HIV-1 Protease

Using high-pressure NMR spectroscopy and differential scanning calorimetry, we investigate the folding landscape of the mature HIV-1 protease homodimer. The cooperativity of unfolding was measured in the absence or presence of a symmetric active site inhibitor for the optimized wild type protease (PR), its inactive variant PRD25N, and an extremely multidrug-resistant mutant, PR20. The individual fit of the pressure denaturation profiles gives rise to first order, ∆GNMR, and second order, ∆VNMR (the derivative of ∆GNMR with pressure); apparent thermodynamic parameters for each amide proton considered. Heterogeneity in the apparent ∆VNMR values reflects departure from an ideal cooperative unfolding transition. The narrow to broad distribution of ∆VNMR spanning the extremes from inhibitor-free PR20D25N to PR-DMP323 complex, and distinctively for PRD25N-DMP323 complex, indicated large variations in folding cooperativity. Consistent with this data, the shape of thermal unfolding transitions varies from asymmetric for PR to nearly symmetric for PR20, as dimer-inhibitor ternary complexes. Lack of structural cooperativity was observed between regions located close to the active site, including the hinge and tip of the glycine-rich flaps, and the rest of the protein. These results strongly suggest that inhibitor binding drastically decreases the cooperativity of unfolding by trapping the closed flap conformation in a deep energy minimum. To evade this conformational trap, PR20 evolves exhibiting a smoother folding landscape with nearly an ideal two-state (cooperative) unfolding transition. This study highlights the malleability of retroviral protease folding pathways by illustrating how the selection of mutations under drug pressure remodels the free-energy landscape as a primary mechanism.

Mutations Proximal to Sites of Autoproteolysis and the α-Helix That Co-evolve under Drug Pressure Modulate the Autoprocessing and Vitality of HIV-1 Protease

N-Terminal self-cleavage (autoprocessing) of the HIV-1 protease precursor is crucial for liberating the active dimer. Under drug pressure, evolving mutations are predicted to modulate autoprocessing, and the reduced catalytic activity of the mature protease (PR) is likely compensated by enhanced conformational/dimer stability and reduced susceptibility to self-degradation (autoproteolysis). One such highly evolved, multidrug resistant protease, PR20, bears 19 mutations contiguous to sites of autoproteolysis in retroviral proteases, namely clusters 1-3 comprising residues 30-37, 60-67, and 88-95, respectively, accounting for 11 of the 19 mutations. By systematically replacing corresponding clusters in PR with those of PR20, and vice versa, we assess their influence on the properties mentioned above and observe no strict correlation. A 10-35-fold decrease in the cleavage efficiency of peptide substrates by PR20, relative to PR, is reflected by an only ∼4-fold decrease in the rate of Gag processing with no change in cleavage order. Importantly, optimal N-terminal autoprocessing requires all 19 PR20 mutations as evaluated in vitro using the model precursor TFR-PR20 in which PR is flanked by the transframe region. Substituting PR20 cluster 3 into TFR-PR (TFR-PR(PR20-3)) requires the presence of PR20 cluster 1 and/or 2 for autoprocessing. In accordance, substituting PR clusters 1 and 2 into TFR-PR20 affects the rate of autoprocessing more drastically (>300-fold) compared to that of TFR-PR(PR20-3) because of the cumulative effect of eight noncluster mutations present in TFR-PR20(PR-12). Overall, these studies imply that drug resistance involves a complex synchronized selection of mutations modulating all of the properties mentioned above governing PR regulation and function.

Dimerization of HIV-1 protease occurs through two steps relating to the mechanism of protease dimerization inhibition by darunavir

Dimerization of HIV-1 protease (PR) subunits is an essential process for PR's acquisition of proteolytic activity, which plays a critical role in the maturation of HIV-1. Recombinant wild-type PR (PR(WT)) proved to dimerize, as examined with electrospray ionization mass spectrometry; however, two active site interface PR mutants (PR(T26A) and PR(R87K)) remained monomeric. On the other hand, two termini interface PR mutants (PR(1-C95A) and PR(97/99)) took both monomeric and dimeric forms. Differential scanning fluorimetry indicated that PR(1-C95A) and PR(97/99) dimers were substantially less stable than PR(WT) dimers. These data indicate that intermolecular interactions of two monomers occur first at the active site interface, generating unstable or transient dimers, and interactions at the termini interface subsequently occur, generating stable dimers. Darunavir (DRV), an HIV-1 protease inhibitor, inhibits not only proteolytic activity but also PR dimerization. DRV bound to protease monomers in a one-to-one molar ratio, inhibiting the first step of PR dimerization, whereas conventional protease inhibitors (such as saquinavir) that inhibit enzymatic activity but not dimerization failed to bind to monomers. DRV also bound to mutant PRs containing the transframe region-added PR (TFR-PR(D25N) and TFR-PR(D25N-7AA)), whereas saquinavir did not bind to TFR-PR(D25N) or TFR-PR(D25N-7AA). Notably, DRV failed to bind to mutant PR containing four amino acid substitutions (V32I, L33F, I54M, and I84V) that confer resistance to DRV on HIV-1. To our knowledge, the present report represents the first demonstration of the two-step PR dimerization dynamics and the mechanism of dimerization inhibition by DRV, which should help design further, more potent novel PIs.

The role of select subtype polymorphisms on HIV-1 protease conformational sampling and dynamics

HIV-1 protease is an essential enzyme for viral particle maturation and is a target in the fight against HIV-1 infection worldwide. Several natural polymorphisms are also associated with drug resistance. Here, we utilized both pulsed electron double resonance, also called double electron-electron resonance, and NMR (15)N relaxation measurements to characterize equilibrium conformational sampling and backbone dynamics of an HIV-1 protease construct containing four specific natural polymorphisms commonly found in subtypes A, F, and CRF_01 A/E. Results show enhanced backbone dynamics, particularly in the flap region, and the persistence of a novel conformational ensemble that we hypothesize is an alternative flap orientation of a curled open state or an asymmetric configuration when interacting with inhibitors.

HIV-1 protease with 20 mutations exhibits extreme resistance to clinical inhibitors through coordinated structural rearrangements

The escape mutant of HIV-1 protease (PR) containing 20 mutations (PR20) undergoes efficient polyprotein processing even in the presence of clinical protease inhibitors (PIs). PR20 shows >3 orders of magnitude decreased affinity for PIs darunavir (DRV) and saquinavir (SQV) relative to PR. Crystal structures of PR20 crystallized with yttrium, substrate analogue p2-NC, DRV, and SQV reveal three distinct conformations of the flexible flaps and diminished interactions with inhibitors through the combination of multiple mutations. PR20 with yttrium at the active site exhibits widely separated flaps lacking the usual intersubunit contacts seen in other inhibitor-free dimers. Mutations of residues 35-37 in the hinge loop eliminate interactions and perturb the flap conformation. Crystals of PR20/p2-NC contain one uninhibited dimer with one very open flap and one closed flap and a second inhibitor-bound dimer in the closed form showing six fewer hydrogen bonds with the substrate analogue relative to wild-type PR. PR20 complexes with PIs exhibit expanded S2/S2' pockets and fewer PI interactions arising from coordinated effects of mutations throughout the structure, in agreement with the strikingly reduced affinity. In particular, insertion of the large aromatic side chains of L10F and L33F alters intersubunit interactions and widens the PI binding site through a network of hydrophobic contacts. The two very open conformations of PR20 as well as the expanded binding site of the inhibitor-bound closed form suggest possible approaches for modifying inhibitors to target extreme drug-resistant HIV.

Dynamics of thermodynamically stable, kinetically trapped, and inhibitor-bound states of pepsin

The pepsin folding mechanism involves a prosegment (PS) domain that catalyzes folding, which is then removed, resulting in a kinetically trapped native state. Although native pepsin (Np) is kinetically stable, it is irreversibly denatured due to a large folding barrier, and in the absence of the PS it folds to a more thermodynamically stable denatured state, termed refolded pepsin (Rp). This system serves as a model to understand the nature of kinetic barriers and folding transitions between compact states. Quasielastic neutron scattering (QENS) was used to characterize and compare the flexibility of Np, as a kinetically trapped state, with that of Rp, as a thermodynamically stable fold. Additionally, the dynamics of Np were compared with those of a partially unfolded form and a thermally stabilized, inhibitor-bound form. QENS revealed length-scale-dependent differences between Np and Rp on a picosecond timescale and indicated greater flexibility in Np, leading to the conclusion that kinetic stabilization likely does not correspond to reduced internal dynamics. Furthermore, large differences were observed upon inhibition, indicating that QENS of proteins in solution may prove useful for examining the role of conformational entropy changes in ligand binding.

Inhibition of autoprocessing of natural variants and multidrug resistant mutant precursors of HIV-1 protease by clinical inhibitors

Self-cleavage at the N terminus of HIV-1 protease from the Gag-Pol precursor (autoprocessing) is crucial for stabilizing the protease dimer required for onset of mature-like catalytic activity, viral maturation, and propagation. Among nine clinical protease inhibitors (PIs), darunavir and saquinavir were the most effective in inhibiting wild-type HIV-1 group M precursor autoprocessing, with an IC(50) value of 1-2 μM, 3-5 orders of magnitude higher than their binding affinities to the corresponding mature protease. Accordingly, both group M and N precursor-PI complexes exhibit T(m)s 17-21 °C lower than those of the corresponding mature protease-PI complexes suggestive of markedly reduced stabilities of the precursor dimer-PI ensembles. Autoprocessing of group N (natural variant) and three group M precursors bearing 11-20 mutations associated with multidrug resistance was either weakly responsive or fully unresponsive to inhibitors at concentrations up to a practical limit of approximately 150 μM PI. This observation parallels decreases of up to 8 × 10(3)-fold (e.g., 5 pM to 40 nM) in the binding affinity of darunavir and saquinavir to mature multidrug resistant proteases relative to wild type, suggesting that inhibition of some of these mutant precursors will occur only in the high μM to mM range in extreme PI-resistance, which is an effect arising from coordinated multiple mutations. An extremely darunavir-resistant mutant precursor is more responsive to inhibition by saquinavir. These findings raise the questions whether clinical failure of PI therapy is related to lack of inhibition of autoprocessing and whether specific inhibitors can be designed with low-nM affinity to target autoprocessing.

The L76V drug resistance mutation decreases the dimer stability and rate of autoprocessing of HIV-1 protease by reducing internal hydrophobic contacts

The mature HIV-1 protease (PR) bearing the L76V drug resistance mutation (PR(L76V)) is significantly less stable, with a >7-fold higher dimer dissociation constant (K(d)) of 71 ± 24 nM and twice the sensitivity to urea denaturation (UC(50) = 0.85 M) relative to those of PR. Differential scanning calorimetry showed decreases in T(m) of 12 °C for PR(L76V) in the absence of inhibitors and 5-7 °C in the presence of inhibitors darunavir (DRV), saquinavir (SQV), and lopinavir (LPV), relative to that of PR. Isothermal titration calorimetry gave a ligand dissociation constant of 0.8 nM for DRV, ∼160-fold higher than that of PR, consistent with DRV resistance. Crystal structures of PR(L76V) in complexes with DRV and SQV were determined at resolutions of 1.45-1.46 Å. Compared to the corresponding PR complexes, the mutated Val76 lacks hydrophobic interactions with Asp30, Lys45, Ile47, and Thr74 and exhibits closer interactions with Val32 and Val56. The bound DRV lacks one hydrogen bond with the main chain of Asp30 in PR(L76V) relative to PR, possibly accounting for the resistance to DRV. SQV shows slightly improved polar interactions with PR(L76V) compared to those with PR. Although the L76V mutation significantly slows the N-terminal autoprocessing of the precursor TFR-PR(L76V) to give rise to the mature PR(L76V), the coselected M46I mutation counteracts the effect by enhancing this rate but renders the TFR-PR(M46I/L76V) precursor less responsive to inhibition by 6 μM LPV while preserving inhibition by SQV and DRV. The correlation of lowered stability, higher K(d), and impaired autoprocessing with reduced internal hydrophobic contacts suggests a novel molecular mechanism for drug resistance.

Autocatalytic maturation, physical/chemical properties, and crystal structure of group N HIV-1 protease: relevance to drug resistance

The mature protease from Group N human immunodeficiency virus Type 1 (HIV-1) (PR1(N)) differs in 20 amino acids from the extensively studied Group M protease (PR1(M)) at positions corresponding to minor drug-resistance mutations (DRMs). The first crystal structure (1.09 Å resolution) of PR1(N) with the clinical inhibitor darunavir (DRV) reveals the same overall structure as PR1(M), but with a slightly larger inhibitor-binding cavity. Changes in the 10s loop and the flap hinge propagate to shift one flap away from the inhibitor, whereas L89F and substitutions in the 60s loop perturb inhibitor-binding residues 29-32. However, kinetic parameters of PR1(N) closely resemble those of PR1(M), and calorimetric results are consistent with similar binding affinities for DRV and two other clinical PIs, suggesting that minor DRMs coevolve to compensate for the detrimental effects of drug-specific major DRMs. A miniprecursor (TFR 1-61-PR1(N)) comprising the transframe region (TFR) fused to the N-terminus of PR1(N) undergoes autocatalytic cleavage at the TFR/PR1(N) site concomitant with the appearance of catalytic activity characteristic of the dimeric, mature enzyme. This cleavage is inhibited at an equimolar ratio of precursor to DRV (∼6 μM), which partially stabilizes the precursor dimer from a monomer. However, cleavage at L34/W35 within the TFR, which precedes the TFR 1-61/PR1(N) cleavage at pH ≤ 5, is only partially inhibited. Favorable properties of PR1(N) relative to PR1(M) include its suitability for column fractionation by size under native conditions and >10-fold higher dimer dissociation constant (150 nM). Exploiting these properties may facilitate testing of potential dimerization inhibitors that perturb early precursor processing steps.

Revealing the dimer dissociation and existence of a folded monomer of the mature HIV-2 protease

Purification and in vitro protein-folding schemes were developed to produce monodisperse samples of the mature wild-type HIV-2 protease (PR2), enabling a comprehensive set of biochemical and biophysical studies to assess the dissociation of the dimeric protease. An E37K substitution in PR2 significantly retards autoproteolytic cleavage during expression. Furthermore, it permits convenient measurement of the dimer dissociation of PR2(E37K) (elevated K(d) approximately 20 nM) by enzyme kinetics. Differential scanning calorimetry reveals a T(m) of 60.5 for PR2 as compared with 65.7 degrees C for HIV-1 protease (PR1). Consistent with weaker binding of the clinical inhibitor darunavir (DRV) to PR2, the T(m) of PR2 increases by 14.8 degrees C in the presence of DRV as compared with 22.4 degrees C for PR1. Dimer interface mutations, such as a T26A substitution in the active site (PR2(T26A)) or a deletion of the C-terminal residues 96-99 (PR2(1-95)), drastically increase the K(d) (>10(5)-fold). PR2(T26A) and PR2(1-95) consist predominantly of folded monomers, as determined by nuclear magnetic resonance (NMR) and size-exclusion chromatography coupled with multiangle light scattering and refractive index measurements (SMR), whereas wild-type PR2 and its active-site mutant PR2(D25N) are folded dimers. Addition of twofold excess active-site inhibitor promotes dimerization of PR2(T26A) but not of PR2(1-95), indicating that subunit interactions involving the C-terminal residues are crucial for dimer formation. Use of SMR and NMR with PR2 facilitates probing for potential inhibitors that restrict protein folding and/or dimerization and, thus, may provide insights for the future design of inhibitors to circumvent drug resistance.