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

Exploring HIV-1 Maturation: A New Frontier in Antiviral Development

HIV-1 virion maturation is an essential step in the viral replication cycle to produce infectious virus particles. Gag and Gag-Pol polyproteins are assembled at the plasma membrane of the virus-producer cells and bud from it to the extracellular compartment. The newly released progeny virions are initially immature and noninfectious. However, once the Gag polyprotein is cleaved by the viral protease in progeny virions, the mature capsid proteins assemble to form the fullerene core. This core, harboring two copies of viral genomic RNA, transforms the virion morphology into infectious virus particles. This morphological transformation is referred to as maturation. Virion maturation influences the distribution of the Env glycoprotein on the virion surface and induces conformational changes necessary for the subsequent interaction with the CD4 receptor. Several host factors, including proteins like cyclophilin A, metabolites such as IP6, and lipid rafts containing sphingomyelins, have been demonstrated to have an influence on virion maturation. This review article delves into the processes of virus maturation and Env glycoprotein recruitment, with an emphasis on the role of host cell factors and environmental conditions. Additionally, we discuss microscopic technologies for assessing virion maturation and the development of current antivirals specifically targeting this critical step in viral replication, offering long-acting therapeutic options.

Elucidating the Substrate Envelope of Enterovirus 68-3C Protease: Structural Basis of Specificity and Potential Resistance

Enterovirus-D68 (EV68) has emerged as a global health concern over the last decade with severe symptomatic infections resulting in long-lasting neurological deficits and death. Unfortunately, there are currently no FDA-approved antiviral drugs for EV68 or any other non-polio enterovirus. One particularly attractive class of potential drugs are small molecules inhibitors, which can target the conserved active site of EV68-3C protease. For other viral proteases, we have demonstrated that the emergence of drug resistance can be minimized by designing inhibitors that leverage the evolutionary constraints of substrate specificity. However, the structural characterization of EV68-3C protease bound to its substrates has been lacking. Here, we have determined the substrate specificity of EV68-3C protease through molecular modeling, molecular dynamics (MD) simulations, and co-crystal structures. Molecular models enabled us to successfully characterize the conserved hydrogen-bond networks between EV68-3C protease and the peptides corresponding to the viral cleavage sites. In addition, co-crystal structures we determined have revealed substrate-induced conformational changes of the protease which involved new interactions, primarily surrounding the S1 pocket. We calculated the substrate envelope, the three-dimensional consensus volume occupied by the substrates within the active site. With the elucidation of the EV68-3C protease substrate envelope, we evaluated how 3C protease inhibitors, AG7088 and SG-85, fit within the active site to predict potential resistance mutations.

The human retroviral-like aspartic protease 1 (ASPRV1): From in vitro studies to clinical correlations

The human retroviral-like aspartic protease 1 (ASPRV1) is a retroviral-like protein that was first identified in the skin due to its expression in the stratum granulosum layer of the epidermis. Accordingly, it is also referred to as skin-specific aspartic protease. Similar to the retroviral polyproteins, the full-length ASPRV1 also undergoes self-proteolysis, the processing of the precursor is necessary for the autoactivation of the protease domain. ASPRV1's functions are well-established at the level of the skin: it is part of the epidermal proteolytic network and has a significant contribution to skin moisturization via the limited proteolysis of filaggrin; it is only natural protein substrate identified so far. Filaggrin and ASPRV1 are also specific for mammalians, these proteins provide unique features for the skins of these species, and the importance of filaggrin processing in hydration is proved by the fact that some ASPRV1 mutations are associated with skin diseases such as ichthyosis. ASPRV1 was also found to be expressed in macrophage-like neutrophil cells, indicating that its functions are not limited to the skin. In addition, differential expression of ASPRV1 was detected in many diseases, with yet unknown significance. The currently known enzymatic characteristics-that had been revealed mainly by in vitro studies-and correlations with pathogenic phenotypes imply potentially important functions in multiple cell types, which makes the protein a promising target of functional studies. In this review we describe the currently available knowledge and future perspective in regard to ASPRV1.

Premature Activation of the HIV-1 Protease Is Influenced by Polymorphisms in the Hinge Region

HIV-1 protease inhibitors are an essential component of antiretroviral therapy. However, drug resistance is a pervasive issue motivating a persistent search for novel therapies. Recent reports found that when protease activates within the host cell's cytosol, it facilitates the pyroptotic killing of infected cells. This has led to speculation that promoting protease activation, rather than inhibiting it, could help to eradicate infected cells and potentially cure HIV-1 infection. Here, we used a nanoscale flow cytometry-based assay to characterize protease resistance mutations and polymorphisms. We quantified protease activity, viral concentration, and premature protease activation and confirmed previous findings that major resistance mutations generally destabilize the protease structure. Intriguingly, we found evidence that common polymorphisms in the hinge domain of protease can influence its susceptibility to premature activation. This suggests that viral heterogeneity could pose a considerable challenge for therapeutic strategies aimed at inducing premature protease activation in the future.

Target-based drug design strategies to overcome resistance to antiviral agents: opportunities and challenges

Viral infections have a major impact in human health. Ongoing viral transmission and escalating selective pressure have the potential to favor the emergence of vaccine- and antiviral drug-resistant viruses. Target-based approaches for the design of antiviral drugs can play a pivotal role in combating drug-resistant challenges. Drug design computational tools facilitate the discovery of novel drugs. This review provides a comprehensive overview of current drug design strategies employed in the field of antiviral drug resistance, illustrated through the description of a series of successful applications. These strategies include technologies that enhance compound-target affinity while minimizing interactions with mutated binding pockets. Furthermore, emerging approaches such as virtual screening, targeted protein/RNA degradation, and resistance analysis during drug design have been harnessed to curtail the emergence of drug resistance. Additionally, host targeting antiviral drugs offer a promising avenue for circumventing viral mutation. The widespread adoption of these refined drug design strategies will effectively address the prevailing challenge posed by antiviral drug resistance.

Substitution Models of Protein Evolution with Selection on Enzymatic Activity

Substitution models of evolution are necessary for diverse evolutionary analyses including phylogenetic tree and ancestral sequence reconstructions. At the protein level, empirical substitution models are traditionally used due to their simplicity, but they ignore the variability of substitution patterns among protein sites. Next, in order to improve the realism of the modeling of protein evolution, a series of structurally constrained substitution models were presented, but still they usually ignore constraints on the protein activity. Here, we present a substitution model of protein evolution with selection on both protein structure and enzymatic activity, and that can be applied to phylogenetics. In particular, the model considers the binding affinity of the enzyme-substrate complex as well as structural constraints that include the flexibility of structural flaps, hydrogen bonds, amino acids backbone radius of gyration, and solvent-accessible surface area that are quantified through molecular dynamics simulations. We applied the model to the HIV-1 protease and evaluated it by phylogenetic likelihood in comparison with the best-fitting empirical substitution model and a structurally constrained substitution model that ignores the enzymatic activity. We found that accounting for selection on the protein activity improves the fitting of the modeled functional regions with the real observations, especially in data with high molecular identity, which recommends considering constraints on the protein activity in the development of substitution models of evolution.

ER-trafficking triggers NRF1 ubiquitination to promote its proteolytic activation

The transcription factor NRF1 resides in the endoplasmic reticulum (ER) and is constantly transported to the cytosol for proteasomal degradation. However, when the proteasome is defective, NRF1 escapes degradation and undergoes proteolytic cleavage by the protease DDI2, generating a transcriptionally active form that restores proteostasis, including proteasome function. The mechanisms that regulate NRF1 proteolytic activation and transcriptional potential remain poorly understood. This study demonstrates that the ER is a crucial regulator of NRF1 function by orchestrating its ubiquitination through the E3 ubiquitin ligase HRD1. We show that HRD1-mediated NRF1 ubiquitination is necessary for DDI2-mediated processing in cells. Furthermore, we found that deficiency in both RAD23A and RAD23B impaired DDI2-mediated NRF1 processing, indicating that these genes are essential components of the DDI2 proteolytic machinery. Our findings highlight the intricate mechanism by which the ER activates NRF1 to coordinate the transcriptional activity of an adaptation response in cells.

Systematic Analyses of the Resistance Potential of Drugs Targeting SARS-CoV-2 Main Protease

Drugs that target the main protease (Mpro) of SARS-CoV-2 are effective therapeutics that have entered clinical use. Wide-scale use of these drugs will apply selection pressure for the evolution of resistance mutations. To understand resistance potential in Mpro, we performed comprehensive surveys of amino acid changes that can cause resistance to nirmatrelvir (Pfizer), and ensitrelvir (Xocova) in a yeast screen. We identified 142 resistance mutations for nirmatrelvir and 177 for ensitrelvir, many of which have not been previously reported. Ninety-nine mutations caused apparent resistance to both inhibitors, suggesting likelihood for the evolution of cross-resistance. The mutation with the strongest drug resistance score against nirmatrelvir in our study (E166V) was the most impactful resistance mutation recently reported in multiple viral passaging studies. Many mutations that exhibited inhibitor-specific resistance were consistent with the distinct interactions of each inhibitor in the substrate binding site. In addition, mutants with strong drug resistance scores tended to have reduced function. Our results indicate that strong pressure from nirmatrelvir or ensitrelvir will select for multiple distinct-resistant lineages that will include both primary resistance mutations that weaken interactions with drug while decreasing enzyme function and compensatory mutations that increase enzyme activity. The comprehensive identification of resistance mutations enables the design of inhibitors with reduced potential of developing resistance and aids in the surveillance of drug resistance in circulating viral populations.

Integrative modeling of diverse protein-peptide systems using CABS-dock

The CABS model can be applied to a wide range of protein-protein and protein-peptide molecular modeling tasks, such as simulating folding pathways, predicting structures, docking, and analyzing the structural dynamics of molecular complexes. In this work, we use the CABS-dock tool in two diverse modeling tasks: 1) predicting the structures of amyloid protofilaments and 2) identifying cleavage sites in the peptide substrates of proteolytic enzymes. In the first case, simulations of the simultaneous docking of amyloidogenic peptides indicated that the CABS model can accurately predict the structures of amyloid protofilaments which have an in-register parallel architecture. Scoring based on a combination of symmetry criteria and estimated interaction energy values for bound monomers enables the identification of protofilament models that closely match their experimental structures for 5 out of 6 analyzed systems. For the second task, it has been shown that CABS-dock coarse-grained docking simulations can be used to identify the positions of cleavage sites in the peptide substrates of proteolytic enzymes. The cleavage site position was correctly identified for 12 out of 15 analyzed peptides. When combined with sequence-based methods, these docking simulations may lead to an efficient way of predicting cleavage sites in degraded proteins. The method also provides the atomic structures of enzyme-substrate complexes, which can give insights into enzyme-substrate interactions that are crucial for the design of new potent inhibitors.

Accelerating antiviral drug discovery: lessons from COVID-19

During the coronavirus disease 2019 (COVID-19) pandemic, a wave of rapid and collaborative drug discovery efforts took place in academia and industry, culminating in several therapeutics being discovered, approved and deployed in a 2-year time frame. This article summarizes the collective experience of several pharmaceutical companies and academic collaborations that were active in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) antiviral discovery. We outline our opinions and experiences on key stages in the small-molecule drug discovery process: target selection, medicinal chemistry, antiviral assays, animal efficacy and attempts to pre-empt resistance. We propose strategies that could accelerate future efforts and argue that a key bottleneck is the lack of quality chemical probes around understudied viral targets, which would serve as a starting point for drug discovery. Considering the small size of the viral proteome, comprehensively building an arsenal of probes for proteins in viruses of pandemic concern is a worthwhile and tractable challenge for the community.

SARS-CoV-2 Mpro Protease Variants of Concern Display Altered Viral Substrate and Cell Host Target Galectin-8 Processing but Retain Sensitivity toward Antivirals

The main protease of SARS-CoV-2 (M^pro) is the most promising drug target against coronaviruses due to its essential role in virus replication. With newly emerging variants there is a concern that mutations in M^pro may alter the structural and functional properties of protease and subsequently the potency of existing and potential antivirals. We explored the effect of 31 mutations belonging to 5 variants of concern (VOCs) on catalytic parameters and substrate specificity, which revealed changes in substrate binding and the rate of cleavage of a viral peptide. Crystal structures of 11 M^pro mutants provided structural insight into their altered functionality. Additionally, we show M^pro mutations influence proteolysis of an immunomodulatory host protein Galectin-8 (Gal-8) and a subsequent significant decrease in cytokine secretion, providing evidence for alterations in the escape of host-antiviral mechanisms. Accordingly, mutations associated with the Gamma VOC and highly virulent Delta VOC resulted in a significant increase in Gal-8 cleavage. Importantly, IC50s of nirmatrelvir (Pfizer) and our irreversible inhibitor AVI-8053 demonstrated no changes in potency for both drugs for all mutants, suggesting M^pro will remain a high-priority antiviral drug candidate as SARS-CoV-2 evolves.

Comprehending the Structure, Dynamics, and Mechanism of Action of Drug-Resistant HIV Protease

Since the emergence of the Human Immunodeficiency Virus (HIV) in the 1980s, strategies to combat HIV-AIDS are continuously evolving. Among the many tested targets to tackle this virus, its protease enzyme (PR) was proven to be an attractive option that brought about numerous research publications and ten FDA-approved drugs to inhibit the PR activity. However, the drug-induced mutations in the enzyme made these small molecule inhibitors ineffective with prolonged usage. The research on HIV PR, therefore, remains a thrust area even today. Through this review, we reiterate the importance of understanding the various structural and functional components of HIV PR in redesigning the structure-based small molecule inhibitors. We also discuss at length the currently available FDA-approved drugs and how these drug molecules induced mutations in the enzyme structure. We then recapitulate the reported mechanisms on how these drug-resistant variants remain sufficiently active to cleave the natural substrates. We end with the future scope covering the recently proposed strategies that show promise to deal with the mutations.

The CARD8 inflammasome in HIV infection

The biggest challenge to immune control of HIV infection is the rapid within-host viral evolution, which allows selection of viral variants that escape from T cell and antibody recognition. Thus, it is impossible to clear HIV infection without targeting "immutable" components of the virus. Unlike the adaptive immune system that recognizes cognate epitopes, the CARD8 inflammasome senses the essential enzymatic activity of the HIV-1 protease, which is immutable for the virus. Hence, all subtypes of HIV clinical isolates can be recognized by CARD8. In HIV-infected cells, the viral protease is expressed as a subunit of the viral Gag-Pol polyprotein and remains functionally inactive prior to viral budding. A class of anti-HIV drugs, the non-nucleoside reverse transcriptase inhibitors (NNRTIs), can promote Gag-pol dimerization and subsequent premature intracellular activation of the viral protease. NNRTI treatment triggers CARD8 inflammasome activation, which leads to pyroptosis of HIV-infected CD4⁺ T cells and macrophages. Targeting the CARD8 inflammasome can be a potent and broadly effective strategy for HIV eradication.

The Effect of Treatment-Associated Mutations on HIV Replication and Transmission Cycles

HIV/AIDS mortality has been decreasing over the last decade. While promising, this decrease correlated directly with increased use of antiretroviral drugs. As a natural consequence of its high mutation rate, treatments provide selection pressure that promotes the natural selection of escape mutants. Individuals may acquire drug-naive strains, or those that have already mutated due to treatment. Even within a host, mutation affects HIV tropism, where initial infection begins with R5-tropic virus, but the clinical transition to AIDS correlates with mutations that lead to an X4-tropic switch. Furthermore, the high mutation rate of HIV has spelled failure for all attempts at an effective vaccine. Pre-exposure drugs are currently the most effective drug-based preventatives, but their effectiveness is also threatened by viral mutation. From attachment and entry to assembly and release, the steps in the replication cycle are also discussed to describe the drug mechanisms and mutations that arise due to those drugs. Revealing the patterns of HIV-1 mutations, their effects, and the coordinated attempt to understand and control them will lead to effective use of current preventative measures and treatment options, as well as the development of new ones.

Molecular basis for reduced cleavage activity and drug resistance in D30N HIV-1 protease

Nelfinavir is one of the FDA-approved HIV-1 protease inhibitors and a part of highly active anti-retroviral therapy (HAART) for the treatment of HIV-AIDS. Nelfinavir was the first HIV-1 protease inhibitor to be approved as a paediatric formulation. The application of HAART had resulted in significant improvement in the lives of AIDS patients. However, the emergence of drug resistance in HIV-1 protease has limited the use of many of these drugs including nelfinavir. A unique mutation observed frequently in patients treated with nelfinavir is D30N as it is selected exclusively by nelfinavir. The D30N mutation imparts very high resistance to nelfinavir but unlike other primary mutations does not give cross-resistance to the majority of other drugs. D30N mutation also significantly reduces cleavage activity of HIV-1 protease and affects viral fitness. Here, we have determined crystal structures of D30N HIV-1 protease in unliganded form and in complex with nelfinavir. These structures provide the rationale for reduced cleavage activity and the molecular basis of drug resistance induced by D30N mutation. The loss of coulombic interaction part of a crucial hydrogen bond between the drug and the protease is likely to play a major role in reduced affinity and resistance towards nelfinavir. The decreased catalytic activity of D30N HIV-1 protease due to altered interaction with the substrates and reduced stability of folding core may be the reason for the reduced replicative capacity of the virus harboring mutant HIV-1 protease.Communicated by Ramaswamy H. Sarma.

Resistor: An algorithm for predicting resistance mutations via Pareto optimization over multistate protein design and mutational signatures

Resistance to pharmacological treatments is a major public health challenge. Here, we introduce Resistor-a structure- and sequence-based algorithm that prospectively predicts resistance mutations for drug design. Resistor computes the Pareto frontier of four resistance-causing criteria: the change in binding affinity (ΔKa) of the (1) drug and (2) endogenous ligand upon a protein's mutation; (3) the probability a mutation will occur based on empirically derived mutational signatures; and (4) the cardinality of mutations comprising a hotspot. For validation, we applied Resistor to EGFR and BRAF kinase inhibitors treating lung adenocarcinoma and melanoma. Resistor correctly identified eight clinically significant EGFR resistance mutations, including the erlotinib and gefitinib "gatekeeper" T790M mutation and five known osimertinib resistance mutations. Furthermore, Resistor predictions are consistent with BRAF inhibitor sensitivity data from both retrospective and prospective experiments using KinCon biosensors. Resistor is available in the open-source protein design software OSPREY.

Sequence dependencies and biophysical features both govern cleavage of diverse cut-sites by HIV protease

The infectivity of HIV-1 requires its protease (PR) cleave multiple cut-sites with low sequence similarity. The diversity of cleavage sites has made it challenging to investigate the underlying sequence properties that determine binding and turnover of substrates by PR. We engineered a mutational scanning approach utilizing yeast display, flow cytometry, and deep sequencing to systematically measure the impacts of all individual amino acid changes at 12 positions in three different cut-sites (MA/CA, NC/p1, and p1/p6). The resulting fitness landscapes revealed common physical features that underlie cutting of all three cut-sites at the amino acid positions closest to the scissile bond. In contrast, positions more than two amino acids away from the scissile bond exhibited a strong dependence on the sequence background of the rest of the cut-site. We observed multiple amino acid changes in cut-sites that led to faster cleavage rates, including a preference for negative charge five and six amino acids away from the scissile bond at locations where the surface of protease is positively charged. Analysis of individual cut sites using full-length matrix-capsid proteins indicate that long-distance sequence context can contribute to cutting efficiency such that analyses of peptides or shorter engineered constructs including those in this work should be considered carefully. This work provides a framework for understanding how diverse substrates interact with HIV-1 PR and can be extended to investigate other viral PRs with similar properties.

Defining the substrate envelope of SARS-CoV-2 main protease to predict and avoid drug resistance

Coronaviruses can evolve and spread rapidly to cause severe disease morbidity and mortality, as exemplified by SARS-CoV-2 variants of the COVID-19 pandemic. Although currently available vaccines remain mostly effective against SARS-CoV-2 variants, additional treatment strategies are needed. Inhibitors that target essential viral enzymes, such as proteases and polymerases, represent key classes of antivirals. However, clinical use of antiviral therapies inevitably leads to emergence of drug resistance. In this study we implemented a strategy to pre-emptively address drug resistance to protease inhibitors targeting the main protease (Mpro) of SARS-CoV-2, an essential enzyme that promotes viral maturation. We solved nine high-resolution cocrystal structures of SARS-CoV-2 Mpro bound to substrate peptides and six structures with cleavage products. These structures enabled us to define the substrate envelope of Mpro, map the critical recognition elements, and identify evolutionarily vulnerable sites that may be susceptible to resistance mutations that would compromise binding of the newly developed Mpro inhibitors. Our results suggest strategies for developing robust inhibitors against SARS-CoV-2 that will retain longer-lasting efficacy against this evolving viral pathogen.

Beyond Inhibition: A Novel Strategy of Targeting HIV-1 Protease to Eliminate Viral Reservoirs

HIV-1 protease (PR) is a viral enzyme that cleaves the Gag and Gag-Pol polyprotein precursors to convert them into their functional forms, a process which is essential to generate infectious viral particles. Due to its broad substrate specificity, HIV-1 PR can also cleave certain host cell proteins. Several studies have identified host cell substrates of HIV-1 PR and described the potential impact of their cleavage on HIV-1-infected cells. Of particular interest is the interaction between PR and the caspase recruitment domain-containing protein 8 (CARD8) inflammasome. A recent study demonstrated that CARD8 can sense HIV-1 PR activity and induce cell death. While PR typically has low levels of intracellular activity prior to viral budding, premature PR activation can be achieved using certain non-nucleoside reverse transcriptase inhibitors (NNRTIs), resulting in CARD8 cleavage and downstream pyroptosis. Used together with latency reversal agents, the induction of premature PR activation to trigger CARD8-mediated cell killing may help eliminate latent reservoirs in people living with HIV. This represents a novel strategy of utilizing PR as an antiviral target through premature activation rather than inhibition. In this review, we discuss the viral and host substrates of HIV-1 protease and highlight potential applications and advantages of targeting CARD8 sensing of HIV-1 PR.

Antiretroviral Drug-Resistance Mutations on the Gag Gene: Mutation Dynamics during Analytic Treatment Interruption among Individuals Experiencing Virologic Failure

We describe drug-resistance mutation dynamics of the gag gene among individuals under antiretroviral virologic failure who underwent analytical treatment interruption (ATI). These mutations occur in and around the cleavage sites that form the particles that become the mature HIV-1 virus. The study involved a 12-week interruption in antiretroviral therapy (ART) and sequencing of the gag gene in 38 individuals experiencing virologic failure and harboring triple-class resistant HIV strains. Regions of the gag gene surrounding the NC-p2 and p1-p6 cleavage sites were sequenced at baseline before ATI and after 12 weeks from plasma HIV RNA using population-based Sanger sequencing. Fourteen of the sixteen patients sequenced presented at least one mutation in the gag gene at baseline, with an average of 4.93 mutations per patient. All the mutations had reverted to the wild type by the end of the study. Mutations in the gag gene complement mutations in the pol gene to restore HIV fitness. Those mutations around cleavage sites and within substrates contribute to protease inhibitor resistance and difficulty in re-establishing effective virologic suppression. ART interruption in the presence of antiretroviral resistant HIV strains was used here as a practical measure for more adapted HIV profiles in the absence of ART selective pressure.

Potential Resistance of SARS-CoV-2 Main Protease (Mpro) against Protease Inhibitors: Lessons Learned from HIV-1 Protease

Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome 2 (SARS-CoV-2), has been one of the most devastating pandemics of recent times. The lack of potent novel antivirals had led to global health crises; however, emergence and approval of potent inhibitors of the viral main protease (Mpro), such as Pfizer's newly approved nirmatrelvir, offers hope not only in the therapeutic front but also in the context of prophylaxis against the infection. By their nature, RNA viruses including human immunodeficiency virus (HIV) have inherently high mutation rates, and lessons learnt from previous and currently ongoing pandemics have taught us that these viruses can easily escape selection pressure through mutation of vital target amino acid residues in monotherapeutic settings. In this paper, we review nirmatrelvir and its binding to SARS-CoV-2 Mpro and draw a comparison to inhibitors of HIV protease that were rendered obsolete by emergence of resistance mutations, emphasizing potential pitfalls in the design of inhibitors that may be of important relevance to the long-term use of novel inhibitors against SARS-CoV-2.

Co-evolution of drug resistance and broadened substrate recognition in HIV protease variants isolated from an Escherichia coli genetic selection system

A genetic selection system for activity of HIV protease is described that is based on a synthetic substrate constructed as a modified AraC regulatory protein that when cleaved stimulate L-arabinose metabolism in an Escherichia coli araC strain. Growth stimulation on selective plates was shown to depend on active HIV protease and the scissile bond in the substrate. In addition, the growth of cells correlated well with the established cleavage efficiency of the sites in the viral polyprotein, Gag, when these sites were individually introduced into the synthetic substate of the selection system. Plasmids encoding protease variants selected based on stimulation of cell growth in the presence of saquinavir or cleavage of a site not cleaved by wild-type protease, were indistinguishable with respect to both phenotypes. Also, both groups of selected plasmids encoded side chain substitutions known from clinical isolates or displayed different side chain substitutions but at identical positions. One highly frequent side chain substitution, E34V, not regarded as a major drug resistance substitution was found in variants obtained under both selective conditions and is suggested to improve protease processing of the synthetic substrate. This substitution is away from the substrate-binding cavity and together with other substitutions in the selected reading frames supports the previous suggestion of a substrate-binding site extended from the active site binding pocket itself.

Analyses of HIV proteases variants at the threshold of viability reveals relationships between processing efficiency and fitness

Investigating the relationships between protein function and fitness provides keys for understanding biochemical mechanisms that underly evolution. Mutations with partial fitness defects can delineate the threshold of biochemical function required for viability. We utilized a previous deep mutational scan of HIV-1 protease (PR) to identify variants with 15–45 per cent defects in replication and analysed the biochemical function of eight variants (L10M, L10S, V32C, V32I, A71V, A71S, Q92I, Q92N). We purified each variant and assessed the efficiency of peptide cleavage for three cut sites (MA-CA, TF-PR, and PR-RT) as well as gel-based analyses of processing of purified Gag. The cutting activity of at least one site was perturbed relative to WT protease for all variants, consistent with cutting activity being a primary determinant of fitness effects. We examined the correlation of fitness defects with cutting activity of different sites. MA-CA showed the weakest correlation (R² = 0.02) with fitness, suggesting relatively weak coupling with viral replication. In contrast, cutting of the TF-PR site showed the strongest correlation with fitness (R² = 0.53). Cutting at the TF-PR site creates a new PR protein with a free N-terminus that is critical for activity. Our findings indicate that increasing the pool of active PR is rate limiting for viral replication, making this an ideal step to target with inhibitors.

Viral proteases: Structure, mechanism and inhibition

Viral proteases are diverse in structure, oligomeric state, catalytic mechanism, and substrate specificity. This chapter focuses on proteases from viruses that are relevant to human health: human immunodeficiency virus subtype 1 (HIV-1), hepatitis C (HCV), human T-cell leukemia virus type 1 (HTLV-1), flaviviruses, enteroviruses, and coronaviruses. The proteases of HIV-1 and HCV have been successfully targeted for therapeutics, with picomolar FDA-approved drugs currently used in the clinic. The proteases of HTLV-1 and the other virus families remain emerging therapeutic targets at different stages of the drug development process. This chapter provides an overview of the current knowledge on viral protease structure, mechanism, substrate recognition, and inhibition. Particular focus is placed on recent advances in understanding the molecular basis of diverse substrate recognition and resistance, which is essential toward designing novel protease inhibitors as antivirals.

Multiscale Approach for Computing Gated Ligand Binding from Molecular Dynamics and Brownian Dynamics Simulations

We develop an approach to characterize the effects of gating by a multiconformation protein consisting of macrostate conformations that are either accessible or inaccessible to ligand binding. We first construct a Markov state model of the apo-protein from atomistic molecular dynamics simulations from which we identify macrostates and their conformations, compute their relative macrostate populations and interchange kinetics, and structurally characterize them in terms of ligand accessibility. We insert the calculated first-order rate constants for conformational transitions into a multistate gating theory from which we derive a gating factor γ that quantifies the degree of conformational gating. Applied to HIV-1 protease, our approach yields a kinetic network of three accessible (semi-open, open, and wide-open) and two inaccessible (closed and a newly identified, "parted") macrostate conformations. The parted conformation sterically partitions the active site, suggesting a possible role in product release. We find that the binding kinetics of drugs and drug-like inhibitors to HIV-1 protease falls in the slow gating regime. However, because γ = 0.75, conformational gating only modestly slows ligand binding. Brownian dynamics simulations of the diffusional association of eight inhibitors to the protease─having a wide range of experimental association constants (∼10⁴-10¹⁰ M^-1 s^-1)─yields gated rate constants in the range of ∼0.5-5.7 × 10⁸ M^-1 s^-1. This indicates that, whereas the association rate of some inhibitors could be described by the model, for many inhibitors either subsequent conformational transitions or alternate binding mechanisms may be rate-limiting. For systems known to be modulated by conformational gating, the approach could be scaled computationally efficiently to screen association kinetics for a large number of ligands.

Stephan Oroszlan and the Proteolytic Processing of Retroviral Proteins: Following A Pro

Steve Oroszlan determined the sequences at the ends of virion proteins for a number of different retroviruses. This work led to the insight that the amino-terminal amino acid of the mature viral CA protein is always proline. In this remembrance, we review Steve’s work that led to this insight and show how that insight was a necessary precursor to the work we have done in the subsequent years exploring the cleavage rate determinants of viral protease processing sites and the multiple roles the amino-terminal proline of CA plays after protease cleavage liberates it from its position in a protease processing site.

Structures of Substrate Complexes of Foamy Viral Protease-Reverse Transcriptase

Reverse transcriptases (RTs) use their DNA polymerase and RNase H activities to catalyze the conversion of single-stranded RNA to double-stranded DNA (dsDNA), a crucial process for the replication of retroviruses. Foamy viruses (FVs) possess a unique RT, which is a fusion with the protease (PR) domain. The mechanism of substrate binding by this enzyme has been unknown. Here, we report a crystal structure of monomeric full-length marmoset FV (MFV) PR-RT in complex with an RNA/DNA hybrid substrate. We also describe a structure of MFV PR-RT with an RNase H deletion in complex with a dsDNA substrate in which the enzyme forms an asymmetric homodimer. Cryo-electron microscopy reconstruction of the full-length MFV PR-RT–dsDNA complex confirmed the dimeric architecture. These findings represent the first structural description of nucleic acid binding by a foamy viral RT and demonstrate its ability to change its oligomeric state depending on the type of bound nucleic acid.

IMPORTANCE Reverse transcriptases (RTs) are intriguing enzymes converting single-stranded RNA to dsDNA. Their activity is essential for retroviruses, which are divided into two subfamilies differing significantly in their life cycles: Orthoretrovirinae and Spumaretrovirinae. The latter family is much more ancient and comprises five genera. A unique feature of foamy viral RTs is that they contain N-terminal protease (PR) domains, which are not present in orthoretroviral enzymes. So far, no structural information for full-length foamy viral PR-RT interacting with nucleic substrates has been reported. Here, we present crystal and cryo-electron microscopy structures of marmoset foamy virus (MFV) PR-RT. These structures revealed the mode of binding of RNA/DNA and dsDNA substrates. Moreover, unexpectedly, the structures and biochemical data showed that foamy viral PR-RT can adopt both a monomeric configuration, which is observed in our structures in the presence of an RNA/DNA hybrid, and an asymmetric dimer arrangement, which we observed in the presence of dsDNA.

Elasticity-Associated Functionality and Inhibition of the HIV Protease

HIV protease plays a critical role in the life cycle of the virus through the generation of mature and infectious virions. Detailed knowledge of the structure of the enzyme and its substrate has led to the development of protease inhibitors. However, the development of resistance to all currently available protease inhibitors has contributed greatly to the decreased success of antiretroviral therapy. When therapy failure occurs, multiple mutations are found within the protease sequence starting with primary mutations, which directly impact inhibitor binding, which can also negatively impact viral fitness and replicative capacity by decreasing the binding affinity of the natural substrates to the protease. As such, secondary mutations which are located outside of the active site region accumulate to compensate for the recurrently deleterious effects of primary mutations. However, the resistance mechanism of these secondary mutations is not well understood, but what is known is that these secondary mutations contribute to resistance in one of two ways, either through increasing the energetic penalty associated with bringing the protease into the closed conformation, or, through decreasing the stability of the protein/drug complex in a manner that increases the dissociation rate of the drug, leading to diminished inhibition. As a result, the elasticity of the enzyme-substrate complex has been implicated in the successful recognition and catalysis of the substrates which may be inferred to suggest that the elasticity of the enzyme/drug complex plays a role in resistance. A realistic representation of the dynamic nature of the protease may provide a more powerful tool in structure-based drug design algorithms.

Deciphering Antifungal Drug Resistance in Pneumocystis jirovecii DHFR with Molecular Dynamics and Machine Learning

Drug resistance impacts the effectiveness of many new therapeutics. Mutations in the therapeutic target confer resistance; however, deciphering which mutations, often remote from the enzyme active site, drive resistance is challenging. In a series of Pneumocystis jirovecii dihydrofolate reductase variants, we elucidate which interactions are key bellwethers to confer resistance to trimethoprim using homology modeling, molecular dynamics, and machine learning. Six molecular features involving mainly residues that did not vary were the best indicators of resistance.

Immature HIV-1 assembles from Gag dimers leaving partial hexamers at lattice edges as potential substrates for proteolytic maturation

HIV-1 particle assembly is driven by the viral Gag protein, which oligomerizes into a hexameric array on the inner surface of the viral envelope, forming a truncated spherical lattice containing large and small gaps. Gag is then cut by the viral protease, disassembles, and rearranges to form the mature, infectious virus. Here, we present structures and molecular dynamics simulations of the edges of the immature Gag lattice. Our analysis shows that Gag dimers are the basic assembly unit of the HIV-1 particle, lattice edges are partial hexamers, and partial hexamers are prone to structural changes allowing protease to cut Gag. These findings provide insights into assembly of the immature virus, its structure, and how it disassembles during maturation.

The CA (capsid) domain of immature HIV-1 Gag and the adjacent spacer peptide 1 (SP1) play a key role in viral assembly by forming a lattice of CA hexamers, which adapts to viral envelope curvature by incorporating small lattice defects and a large gap at the site of budding. This lattice is stabilized by intrahexameric and interhexameric CA-CA interactions, which are important in regulating viral assembly and maturation. We applied subtomogram averaging and classification to determine the oligomerization state of CA at lattice edges and found that CA forms partial hexamers. These structures reveal the network of interactions formed by CA-SP1 at the lattice edge. We also performed atomistic molecular dynamics simulations of CA-CA interactions stabilizing the immature lattice and partial CA-SP1 helical bundles. Free energy calculations reveal increased propensity for helix-to-coil transitions in partial hexamers compared to complete six-helix bundles. Taken together, these results suggest that the CA dimer is the basic unit of lattice assembly, partial hexamers exist at lattice edges, these are in a helix-coil dynamic equilibrium, and partial helical bundles are more likely to unfold, representing potential sites for HIV-1 maturation initiation.

Cantilever-centric mechanism of cooperative non-active site mutations in HIV protease: Implications for flap dynamics

The HIV-1 protease is an important drug target in antiretroviral therapy due to the crucial role it plays in viral maturation. A greater understanding of the dynamics of the protease as a result of drug-induced mutations has been successfully elucidated using computational models in the past. We performed induced-fit docking studies and molecular dynamics simulations on the wild-type South African HIV-1 subtype C protease and two non-active site mutation-containing protease variants; HP3 PR and HP4 PR. The HP3 PR contained the I13V, I62V, and V77I mutations while HP4 PR contained the same mutations with the addition of the L33F mutation. The simulations were initiated in a cubic cell universe containing explicit solvent, with the protease variants beginning in the fully closed conformation. The trajectory for each simulation totalled 50 ns. The results indicate that the mutations increase the dynamics of the flap, hinge, fulcrum and cantilever regions when compared to the wild-type protease while in complex with protease inhibitors. Specifically, these mutations result in the protease favouring the semi-open conformation when in complex with inhibitors. Moreover, the HP4 PR adopted curled flap tip conformers which coordinated several water molecules into the active site in a manner that may reduce inhibitor binding affinity. The mutations affected the thermodynamic landscape of inhibitor binding as there were fewer observable chemical contacts between the mutated variants and saquinavir, atazanavir and darunavir. These data help to elucidate the biophysical basis for the selection of cooperative non-active site mutations by the HI virus.

Deciphering Complex Mechanisms of Resistance and Loss of Potency through Coupled Molecular Dynamics and Machine Learning

Drug resistance threatens many critical therapeutics through mutations in the drug target. The molecular mechanisms by which combinations of mutations, especially those remote from the active site, alter drug binding to confer resistance are poorly understood and thus difficult to counteract. A machine learning strategy was developed that coupled parallel molecular dynamics simulations with experimental potency to identify specific conserved mechanisms underlying resistance. Physical features were extracted from the simulations, analyzed, and integrated into one consistent and interpretable elastic network model. To rigorously test this strategy, HIV-1 protease variants with diverse mutations were used, with potencies ranging from picomolar to micromolar to the drug darunavir. Feature reduction resulted in a model with four specific features that predicts for both the training and test sets inhibitor binding free energy within 1 kcal/mol of the experimental value over this entire range of potency. These predictive features are physically interpretable, as they vary specifically with affinity and diagonally transverse across the protease homodimer. This physics-based strategy of parallel molecular dynamics and machine learning captures mechanisms by which complex combinations of mutations confer resistance and identify critical features that serve as bellwethers of affinity, which will be critical in future drug design.

Predicting HIV-1 Protease Cleavage Sites With Positive-Unlabeled Learning

Understanding the substrate specificity of HIV-1 protease plays an essential role in the prevention of HIV infection. A variety of computational models have thus been developed to predict substrate sites that are cleaved by HIV-1 protease, but most of them normally follow a supervised learning scheme to build classifiers by considering experimentally verified cleavable sites as positive samples and unknown sites as negative samples. However, certain noisy can be contained in the negative set, as false negative samples are possibly existed. Hence, the performance of the classifiers is not as accurate as they could be due to the biased prediction results. In this work, unknown substrate sites are regarded as unlabeled samples instead of negative ones. We propose a novel positive-unlabeled learning algorithm, namely PU-HIV, for an effective prediction of HIV-1 protease cleavage sites. Features used by PU-HIV are encoded from different perspectives of substrate sequences, including amino acid identities, coevolutionary patterns and chemical properties. By adjusting the weights of errors generated by positive and unlabeled samples, a biased support vector machine classifier can be built to complete the prediction task. In comparison with state-of-the-art prediction models, benchmarking experiments using cross-validation and independent tests demonstrated the superior performance of PU-HIV in terms of AUC, PR-AUC, and F-measure. Thus, with PU-HIV, it is possible to identify previously unknown, but physiologically existed substrate sites that are able to be cleaved by HIV-1 protease, thus providing valuable insights into designing novel HIV-1 protease inhibitors for HIV treatment.

Drug Design Strategies to Avoid Resistance in Direct-Acting Antivirals and Beyond

Drug resistance is prevalent across many diseases, rendering therapies ineffective with severe financial and health consequences. Rather than accepting resistance after the fact, proactive strategies need to be incorporated into the drug design and development process to minimize the impact of drug resistance. These strategies can be derived from our experience with viral disease targets where multiple generations of drugs had to be developed to combat resistance and avoid antiviral failure. Significant efforts including experimental and computational structural biology, medicinal chemistry, and machine learning have focused on understanding the mechanisms and structural basis of resistance against direct-acting antiviral (DAA) drugs. Integrated methods show promise for being predictive of resistance and potency. In this review, we give an overview of this research for human immunodeficiency virus type 1, hepatitis C virus, and influenza virus and the lessons learned from resistance mechanisms of DAAs. These lessons translate into rational strategies to avoid resistance in drug design, which can be generalized and applied beyond viral targets. While resistance may not be completely avoidable, rational drug design can and should incorporate strategies at the outset of drug development to decrease the prevalence of drug resistance.

A synergy of activity, stability, and inhibitor-interaction of HIV-1 protease mutants evolved under drug-pressure

A clinically-relevant, drug-resistant mutant of HIV-1 protease (PR), termed Flap+_(I54V) and containing L10I, G48V, I54V and V82A mutations, is known to produce significant changes in the entropy and enthalpy balance of drug-PR interactions, compared to wild-type PR. A similar mutant, Flap+_(I54A) , which evolves from Flap+_(I54V) and contains the single change at residue 54 relative to Flap+_(I54V) , does not. Yet, how Flap+_(I54A) behaves in solution is not known. To understand the molecular basis of V54A evolution, we compared nuclear magnetic resonance (NMR) spectroscopy, fluorescence spectroscopy, isothermal titration calorimetry, and enzymatic assay data from four PR proteins: PR (pWT), Flap+_(I54V) , Flap+_(I54A) , and Flap+_(I54) , a control mutant that contains only L10I, G48V and V82A mutations. Our data consistently show that selection to the smaller side chain at residue 54, not only decreases inhibitor affinity, but also restores the catalytic activity.

Probing the Interaction between HIV-1 Protease and the Homodimeric p66/p66' Reverse Transcriptase Precursor by Double Electron-Electron Resonance EPR Spectroscopy

Following excision from the Gag-Pol polyprotein, HIV-1 reverse transcriptase is released as an asymmetric homodimer comprising two p66 subunits that are structurally dissimilar but identical in amino acid sequence. Subsequent cleavage of the RNase H domain from only one of the subunits, denoted p66', results in the formation of the mature p66/p51 enzyme in which catalytic activity resides in the p66 subunit, and the p51 subunit (derived from p66') provides a supporting structural scaffold. Here, we probe the interaction of the p66/p66' asymmetric reverse transcriptase precursor with HIV-1 protease by pulsed Q-band double electron-electron resonance EPR spectroscopy to measure distances between nitroxide labels introduced at surface-engineered cysteine residues. The data suggest that the flexible, exposed linker between the RNaseH and connection domains in the open state of the p66' subunit binds to the active site of protease in a configuration that is similar to that of extended peptide substrates.

Incorporating the Coevolving Information of Substrates in Predicting HIV-1 Protease Cleavage Sites

Human immunodeficiency virus 1 (HIV-1) protease (PR) plays a crucial role in the maturation of the virus. The study of substrate specificity of HIV-1 PR as a new endeavor strives to increase our ability to understand how HIV-1 PR recognizes its various cleavage sites. To predict HIV-1 PR cleavage sites, most of the existing approaches have been developed solely based on the homogeneity of substrate sequence information with supervised classification techniques. Although efficient, these approaches are found to be restricted to the ability of explaining their results and probably provide few insights into the mechanisms by which HIV-1 PR cleaves the substrates in a site-specific manner. In this work, a coevolutionary pattern-based prediction model for HIV-1 PR cleavage sites, namely EvoCleave, is proposed by integrating the coevolving information obtained from substrate sequences with a linear SVM classifier. The experiment results showed that EvoCleave yielded a very promising performance in terms of ROC analysis and f-measure. We also prospectively assessed the biological significance of coevolutionary patterns by applying them to study three fundamental issues of HIV-1 PR cleavage site. The analysis results demonstrated that the coevolutionary patterns offered valuable insights into the understanding of substrate specificity of HIV-1 PR.

Exploring the concerted mechanistic pathway for HIV-1 PR-substrate revealed by umbrella sampling simulation

HIV-1 protease (HIV-1 PR) is an essential enzyme for the replication process of its virus, and therefore considered an important target for the development of drugs against the acquired immunodeficiency syndrome (AIDS). Our previous study shows that the catalytic mechanism of subtype B/C-SA HIV-1 PR follows a one-step concerted acyclic hydrolysis reaction process using a two-layered ONIOM B3LYP/6-31++G(d,p) method. This present work is aimed at exploring the proposed mechanism of the proteolysis catalyzed by HIV-1 PR and to ensure our proposed mechanism is not an artefact of a single theoretical technique. Hence, we present umbrella sampling method that is suitable for calculating potential mean force (PMF) for non-covalent ligand/substrate-enzyme association/dissociation interactions which provide thermodynamic details for molecular recognition. The free activation energy results were computed in terms of PMF analysis within the hybrid QM(DFTB)/MM approach. The theoretical findings suggest that the proposed mechanism corresponds in principle with experimental data. Given our observations, we suggest that the QM/MM MD method can be used as a reliable computational technique to rationalize lead compounds against specific targets such as the HIV-1 protease.

Peptidyl Fluoromethyl Ketones and Their Applications in Medicinal Chemistry

Peptidyl fluoromethyl ketones occupy a pivotal role in the current scenario of synthetic chemistry, thanks to their numerous applications as inhibitors of hydrolytic enzymes. The insertion of one or more fluorine atoms adjacent to a C-terminal ketone moiety greatly modifies the physicochemical properties of the overall substrate, especially by increasing the reactivity of this functionalized carbonyl group toward nucleophiles. The main application of these peptidyl α-fluorinated ketones in medicinal chemistry relies in their ability to strongly and selectively inhibit serine and cysteine proteases. These compounds can be used as probes to study the proteolytic activity of the aforementioned proteases and to elucidate their role in the insurgence and progress on several diseases. Likewise, if the fluorinated methyl ketone moiety is suitably connected to a peptidic backbone, it may confer to the resulting structure an excellent substrate peculiarity and the possibility of being recognized by a specific subclass of human or pathogenic proteases. Therefore, peptidyl fluoromethyl ketones are also currently highly exploited for the target-based design of compounds for the treatment of topical diseases such as various types of cancer and viral infections.

Electrostatics Plays a Crucial Role in HIV-1 Protease Substrate Binding, Drugs Fail to Take Advantage

HIV-1 protease (HIVPR) is an important drug target for combating AIDS. This enzyme is an aspartyl protease that is functionally active in its dimeric form. Nuclear magnetic resonance reports have convincingly shown that a pseudosymmetry exists at the HIVPR active site, where only one of the two aspartates remains protonated over the pH range of 2.5-7.0. To date, all HIVPR-targeted drug design strategies focused on maximizing the size-shape complementarity and van der Waals interactions of the small molecule drugs with the deprotonated, symmetric active site envelope of crystallized HIVPR. However, these strategies were ineffective with the emergence of drug resistant protease variants, primarily due to the steric clashes at the active site. In this study, we traced a specificity in the substrate binding motif that emerges primarily from the asymmetrical electrostatic potential present in the protease active site due to the uneven protonation. Our detailed results from atomistic molecular dynamics simulations show that while such a specific mode of substrate binding involves significant electrostatic interactions, none of the existing drugs or inhibitors could utilize this electrostatic hot spot. As the electrostatic is long-range interaction, it can provide sufficient binding strength without the necessity of increasing the bulkiness of the inhibitors. We propose that introducing the electrostatic component along with optimal fitting at the binding pocket could pave the way for promising designs that might be more effective against both wild type and HIVPR resistant variants.

Viral structural proteins as targets for antivirals

Viral structural proteins are emerging as effective targets for new antivirals. In a viral lifecycle, the capsid must assemble, disassemble, and respond to host proteins, all at the right time and place. These reactions work within a narrow range of conditions, making them susceptible to small molecule interference. In at least three specific viruses, this approach has had met with preliminary success. In rhinovirus and poliovirus, compounds like pleconaril bind capsid and block RNA release. Bevirimat binds to Gag protein in HIV, inhibiting maturation. In Hepatitis B virus, core protein allosteric modulators (CpAMs) promote spontaneous assembly of capsid protein leading to empty and aberrant particles. Despite the biological diversity between viruses and the chemical diversity between antiviral molecules, we observe common features in these antivirals' mechanisms of action. These approaches work by stabilizing protein-protein interactions.

Drug susceptibility and replication capacity of a rare HIV-1 subtype C protease hinge region variant

BACKGROUND

Protease inhibitors form the main component of second-line antiretroviral treatment in South Africa. Despite their efficacy, mutations arising within the HIV-1 gag and protease coding regions contribute to the development of resistance against this class of drug. In this paper we investigate a South African HIV-1 subtype C Gag-protease that contains a hinge region mutation and insertion (N37T↑V).

METHODS

In vitro single-cycle drug susceptibility and viral replication capacity assays were performed on W1201i, a wild-type reference isolate (MJ4) and a chimeric construct (MJ4GagN37T↑VPR). Additionally, enzyme assays were performed on the N37T↑V protease and a wild-type reference protease.

RESULTS

W1201i showed a small (threefold), but significant (P<0.0001) reduction in drug susceptibility to darunavir compared with MJ4. Substitution of W1201i-Gag with MJ4-Gag resulted in an additional small (twofold), but significant (P<0.01) reduction in susceptibility to lopinavir and atazanavir. The W1201i pseudovirus had a significantly (P<0.01) reduced replication capacity (16.4%) compared with the MJ4. However, this was dramatically increased to 164% (P<0.05) when W1201i-Gag was substituted with MJ4-Gag. Furthermore, the N37T↑V protease displayed reduced catalytic processing compared with the SK154 protease.

CONCLUSIONS

Collectively, these data suggest that the N37T↑V mutation and insertion increases viral infectivity and decreases drug susceptibility. These variations are classified as secondary mutations, and indirectly impact inhibitor binding, enzyme fitness and enzyme stability. Additionally, polymorphisms arising in Gag can modify the impact of protease with regards to viral replication and susceptibility to protease inhibitors.

Fine-Tuning of Sequence Specificity by Near Attack Conformations in Enzyme-Catalyzed Peptide Hydrolysis

The catalytic role of near attack conformations (NACs), molecular states that lie on the pathway between the ground state (GS) and transition state (TS) of a chemical reaction, is not understood completely. Using a computational approach that combines Bürgi–Dunitz theory with all-atom molecular dynamics simulations, the role of NACs in catalyzing the first stages of HIV-1 protease peptide hydrolysis was previously investigated using a substrate that represents the recognized SP1-NC cleavage site of the HIV-1 Gag polyprotein. NACs were found to confer no catalytic effect over the uncatalyzed reaction there ( Δ Δ G N ‡ ∼ 0 kcal/mol). Here, using the same approach, the role of NACs across multiple substrates that each represent a further recognized cleavage site is investigated. Overall rate enhancement varies by | Δ Δ G ‡ | ∼ 12–15 kcal/mol across this set, and although NACs contribute a small and approximately constant barrier to the uncatalyzed reaction (< Δ G N ‡ u > = 4.3 ± 0.3 kcal/mol), they are found to contribute little significant catalytic effect ( | Δ Δ G N ‡ | ∼ 0–2 kcal/mol). Furthermore, no correlation is exhibited between NAC contributions and the overall energy barrier ( R 2 = 0.01). However, these small differences in catalyzed NAC contributions enable rates to match those required for the kinetic order of processing. Therefore, NACs may offer an alternative and subtle mode compared to non-NAC contributions for fine-tuning reaction rates during complex evolutionary sequence selection processes—in this case across cleavable polyproteins whose constituents exhibit multiple functions during the virus life-cycle.

Gag-protease coevolution shapes the outcome of lopinavir-inclusive treatment regimens in chronically infected HIV-1 subtype C patients

MOTIVATION

Commonly, protease inhibitor failure is characterized by the development of multiple protease resistance mutations (PRMs). While the impact of PRMs on therapy failure are understood, the introduction of Gag mutations with protease remains largely unclear.

RESULTS

Here, we utilized phylogenetic analyses and Bayesian network learning as tools to understand Gag-protease coevolution and elucidate the pathways leading to Lopinavir failure in HIV-1 subtype C infected patients. Our analyses indicate that while PRMs coevolve in response to drug selection pressure within protease, the Gag mutations added to the existing network while specifically interacting with known Lopinavir failure PRMs. Additionally, the selection of mutations at specific positions in Gag-protease suggests that these coevolving mutational changes occurs to maintain structural integrity during Gag cleavage.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Concerted hydrolysis mechanism of HIV-1 natural substrate against subtypes B and C-SA PR: insight through molecular dynamics and hybrid QM/MM studies

It is well known that understanding the catalytic mechanism of HIV-1 PR is the rationale on which its inhibitors were developed; therefore, a better understanding of the mechanism of natural substrate hydrolysis is important. Herein, the reaction mechanism of HIV-1 natural substrates with subtypes B and common mutant in South Africa (subtype C-SA) protease were studied through transition state modelling, using a general acid-general base (GA-GB) one-step concerted process. The activation free energies of enzyme-substrate complexes were compared based on their rate of hydrolysis using a two-layered ONIOM (B3LYP/6-31++G(d,p):AMBER) method. We expanded our computational model to obtain a better understanding of the mechanism of hydrolysis as well as how the enzyme recognises or chooses the cleavage site of the scissile bonds. Using this model, a potential substrate-based inhibitor could be developed with better potency. The calculated activation energies of natural substrates in our previous study correlated well with experimental data. A similar trend was observed for the Gag and Gag-Pol natural substrates in the present work for both enzyme complexes except for the PR-RT substrate. Natural bond orbital (NBO) analysis was also applied to determine the extent of charge transfer within the QM part of both enzymes considered and the PR-RT natural substrate. The result of this study shows that the method can be utilized as a dependable computational technique to rationalize lead compounds against specific targets.

Characterizing Protein-Ligand Binding Using Atomistic Simulation and Machine Learning: Application to Drug Resistance in HIV-1 Protease

Over the past several decades, atomistic simulations of biomolecules, whether carried out using molecular dynamics or Monte Carlo techniques, have provided detailed insights into their function. Comparing the results of such simulations for a few closely related systems has guided our understanding of the mechanisms by which changes such as ligand binding or mutation can alter the function. The general problem of detecting and interpreting such mechanisms from simulations of many related systems, however, remains a challenge. This problem is addressed here by applying supervised and unsupervised machine learning techniques to a variety of thermodynamic observables extracted from molecular dynamics simulations of different systems. As an important test case, these methods are applied to understand the evasion by human immunodeficiency virus type-1 (HIV-1) protease of darunavir, a potent inhibitor to which resistance can develop via the simultaneous mutation of multiple amino acids. Complex mutational patterns have been observed among resistant strains, presenting a challenge to developing a mechanistic picture of resistance in the protease. In order to dissect these patterns and gain mechanistic insight into the role of specific mutations, molecular dynamics simulations were carried out on a collection of HIV-1 protease variants, chosen to include highly resistant strains and susceptible controls, in complex with darunavir. Using a machine learning approach that takes advantage of the hierarchical nature in the relationships among the sequence, structure, and function, an integrative analysis of these trajectories reveals key details of the resistance mechanism, including changes in the protein structure, hydrogen bonding, and protein-ligand contacts.

Picomolar to Micromolar: Elucidating the Role of Distal Mutations in HIV-1 Protease in Conferring Drug Resistance

Drug resistance continues to be a growing global problem. The efficacy of small molecule inhibitors is threatened by pools of genetic diversity in all systems, including antibacterials, antifungals, cancer therapeutics, and antivirals. Resistant variants often include combinations of active site mutations and distal "secondary" mutations, which are thought to compensate for losses in enzymatic activity. HIV-1 protease is the ideal model system to investigate these combinations and underlying molecular mechanisms of resistance. Darunavir (DRV) binds wild-type (WT) HIV-1 protease with a potency of <5 pM, but we have identified a protease variant that loses potency to DRV 150 000-fold, with 11 mutations in and outside the active site. To elucidate the roles of these mutations in DRV resistance, we used a multidisciplinary approach, combining enzymatic assays, crystallography, and molecular dynamics simulations. Analysis of protease variants with 1, 2, 4, 8, 9, 10, and 11 mutations showed that the primary active site mutations caused ∼50-fold loss in potency (2 mutations), while distal mutations outside the active site further decreased DRV potency from 13 nM (8 mutations) to 0.76 μM (11 mutations). Crystal structures and simulations revealed that distal mutations induce subtle changes that are dynamically propagated through the protease. Our results reveal that changes remote from the active site directly and dramatically impact the potency of the inhibitor. Moreover, we find interdependent effects of mutations in conferring high levels of resistance. These mechanisms of resistance are likely applicable to many other quickly evolving drug targets, and the insights may have implications for the design of more robust inhibitors.

Protease Inhibitors for the Treatment of HIV/AIDS: Recent Advances and Future Challenges

Acquired Immunodeficiency Syndrome (AIDS) is a chronic disease characterized by multiple life-threatening illnesses caused by a retro-virus, Human Immunodeficiency Virus (HIV). HIV infection slowly destroys the immune system and increases the risk of various other infections and diseases. Although, there is no immediate cure for HIV infection/AIDS, several drugs targeting various cruxes of HIV infection are used to slow down the progress of the disease and to boost the immune system. One of the key therapeutic strategies is Highly Active Antiretroviral Therapy (HAART) or ' AIDS cocktail' in a general sense, which is a customized combination of anti-retroviral drugs designed to combat the HIV infection. Since HAART's inception in 1995, this treatment was found to be effective in improving the life expectancy of HIV patients over two decades. Among various classes of HAART treatment regimen, Protease Inhibitors (PIs) are known to be widely used as a major component and found to be effective in treating HIV infection/AIDS. For the past several years, a variety of protease inhibitors have been reported. This review outlines the drug design strategies of PIs, chemical and pharmacological characteristics of some mechanism-based inhibitors, summarizes the recent developments in small molecule based drug discovery with HIV protease as a drug target. Further discussed are the pharmacology, PI drug resistance on HIV PR, adverse effects of HIV PIs and challenges/impediments in the successful application of HIV PIs as an important class of drugs in HAART regimen for the effective treatment of AIDS.