Multi-drug resistance profile of PR20 HIV-1 protease is attributed to distorted conformational and drug binding landscape: molecular dynamics insights

Metabolites profiling and cheminformatics bioprospection of selected medicinal plants against the main protease and RNA-dependent RNA polymerase of SARS-CoV-2

Despite the existence of some vaccines, SARS-CoV-2 (S-2) infections persist for various reasons relating to vaccine reluctance, rapid mutation rate, and an absence of specific treatments targeted to the infection. Due to their availability, low cost and low toxicity, research into potentially repurposing phytometabolites as therapeutic alternatives has gained attention. Therefore, this study explored the antiviral potential of metabolites of some medicinal plants [Spondias mombin, Macaranga barteri and Dicerocaryum eriocarpum (Sesame plant)] identified using liquid chromatography-mass spectrometry (LCMS) as possible inhibitory agents against the S-2 main protease (S-2 MP) and RNA-dependent RNA polymerase (RP) using computational approaches. Molecular docking was used to identify the compounds with the best affinities for the selected therapeutics targets. Afterwards, compounds with poor physicochemical characteristics, pharmacokinetics, and drug-likeness were screened out. The top-ranked compounds were further subjected to a 120-ns molecular dynamics (MD) simulation. Only quercetin 3-O-rhamnoside (-48.77 kcal/mol) had higher binding free energy than the reference standard (zafirlukast) (-44.99 kcal/mol) against S-2 MP. Conversely, all the top-ranked compounds (ellagic acid hexoside, spiraeoside, apigenin-4'-glucoside and chrysoeriol 7-glucuronide) except gnetin L (-24.24 kcal/mol) had higher binding free energy (-55.19 kcal/mol, -52.75 kcal/mol, -47.22 kcal/mol and -43.35 kcal/mol) respectively, against S-2 RP relative to the reference standard (-34.79 kcal/mol). The MD simulations study further revealed that the investigated inhibitors are thermodynamically stable and form structurally compatible complexes that impede the regular operation of the respective S-2 therapeutic targets. Although, these S-2 therapeutic candidates are promising, further in vitro and in vivo evaluation is required and highly recommended.Communicated by Ramaswamy H. Sarma.

Coevolution-based computational approach to detect resistance mechanism of epidermal growth factor receptor

Tyrosine kinase epidermal growth factor receptor (EGFR) correlates the neoplastic cell metastasis, angiogenesis, neoplastic incursion, and apoptosis. Due to the involvement of EGFR in these biological processes, it becomes a most potent target for treating non-small cell lung cancer (NSCLC). The tyrosine kinase inhibitors (TKI) have endorsed high efficacy and anticipation to patients but unfortunately, within a year of treatment, drug targets develop resistance due to mutations. The present study detected the compensatory mutations in EGFR to know the evolutionary mechanism of drug resistance. The results of this study demonstrate that compensatory mutations enlarge the drug-binding pocket which may lead to the altered orientation of the ligand (gefitinib and erlotinib) causing drug resistance. This indicates that coevolutionary forces play a significant role in fine-tuning the structure of EGFR protein against the drugs. The analysis provides insight into the evolution-induced structural aspects of drug resistance changes in EGFR which in turn be useful in designing drugs with better efficacy.

Influence of Clerodendrum volubile leaf extract on doxorubicin-induced toxicity and inhibition of carbonyl reductase mediated metabolism

OBJECTIVES

Doxorubicin (DOX) is a commonly used chemotherapeutic drug. However, its non-target organ toxicities pose a serious problem. This study is to assess the protective role of Clerodendrum volubile leaf extract (CVE) against DOX-induced toxicities in rats. In addition, the inhibitory activities of three phytochemical compounds (Rutin, Gallic acid and Rosmarinic acid) from CVE against Carbonyl reductase 1 (CBR1) were examined.

METHODS

Rats were randomly divided into 5 groups: (a) Control group rats were given 0.9% NaCl as vehicle, (b) DOX group: A single dose of DOX (25 mg/kg; i.p.) was administered and rats were sacrificed 4 days after DOX injection, while groups (c-e) CVE-treated DOX rat groups were given 125, 250 and 500 mg/kg body weight of extracts orally for 12 consecutive days; 8 days before, and 4 days after the DOX administration. Computational techniques were used to determine the inhibitory activities of the compounds against CBR1.

RESULTS

DOX intoxication caused a significant increase (p<0.05) in serum marker enzymes: ALT, AST, ALP, LDH, CK activities. The levels of liver and heart tissues antioxidant parameters: GPx, SOD, CAT, and GSH were significantly (p<0.05) decreased in DOX-intoxicated rats with concomitant elevation of malondialdehyde levels. Pretreatment with CVE reversed the above trends. From the structural analysis, Rutin and RSA exhibited the highest binding free energies against CBR1, and also exhibited structural stability when bound with CBR1.

CONCLUSIONS

Our study indicates the protective effect of CVE when used in combination with doxorubicin thus improving its chemotherapeutic application via inhibition of CBR-mediated metabolism.

HIV-1 protease with 10 lopinavir and darunavir resistance mutations exhibits altered inhibition, structural rearrangements and extreme dynamics

Antiretroviral drug resistance is a therapeutic obstacle for people with HIV. HIV protease inhibitors darunavir and lopinavir are recommended for resistant infections. We characterized a protease mutant (PR10x) derived from a highly resistant clinical isolate including 10 mutations associated with resistance to lopinavir and darunavir. Compared to the wild-type protease, PR10x exhibits ∼3-fold decrease in catalytic efficiency and Ki values of 2-3 orders of magnitude worse for darunavir, lopinavir, and potent investigational inhibitor GRL-519. Crystal structures of the mutant were solved in a ligand-free form and in complex with GRL-519. The structures show altered interactions in the active site, flap-core interface, hydrophobic core, hinge region, and 80s loop compared to the corresponding wild-type protease structures. The ligand-free crystal structure exhibits a highly curled flap conformation which may amplify drug resistance. Molecular dynamics simulations performed for 1 μs on ligand-free dimers showed extremely large fluctuations in the flaps for PR10x compared to equivalent simulations on PR with a single L76V mutation or wild-type protease. This analysis offers insight about the synergistic effects of mutations in highly resistant variants.

Repurposing drugs and identification of inhibitors of integral proteins (spike protein and main protease) of SARS-CoV-2

The outbreak of Coronavirus infection (COVID-19) has prompted the World Health Organisation (WHO) to declare the outbreak, a Public Health Emergency of International concern. As part of the efforts to discover lead compounds for clinical use, 53 molecules were screened using molecular docking and dynamic simulations (MDS) techniques to identify potential inhibitors of SARS-CoV-2 spike protein (COVID-19 Sgp) and main protease (COVID-19 Mpro) or both. Lopinavir (LPV), nelfinavir (NEF), hydroxychloroquine (HCQ), remdesivir (RDV) and an irreversible inhibitor of SARS-CoV (N3) were used as standard drugs for COVID-19 Mpro, while zafirlukast (ZFK) and cefoperazone (CSP)) as standard drugs for COVID-19 Sgp. After 100 ns of MDS, with reference to standard drugs (N3, -52.463 Kcal/mol, NEF, -51.618 Kcal/mol, RDV, -48.780 Kcal/mol, LPV, -46.788 Kcal/mol, DRV, -33.655 Kcal/mol and HCQ, -21.065 Kcal/mol), five molecules, HCR, GRN, C3G, EGCG, and K7G were predicted to be promising inhibitors of COVID-19 Mpro with binding energies of -53.877 kcal/mol, -50.653 Kcal/mol, -48.600 kcal/mol, -47.798 kcal/mol and -46.902 kcal/mol, respectively. These lead molecules were then docked at receptor-binding domain (RBD) of COVID-19 Sgp to examine their inhibitory effects. C3G, GRN and K7G exhibited higher binding energies of -42.310 kcal/mol, -32.210 kcal/mol, -26.922 kcal/mol than the recorded values for the reference drugs (CSP, -35.509 kcal/mol, ZFK, -24.242 kcal/mol), respectively. The results of the binding energy and structural analyses from this study revealed that C3G, GRN and K7G could serve as potential dual inhibitors of COVID-19 Sgp and COVID-19 Mpro, while HCR and EGCG would be inhibitors of COVID-19 Mpro.Communicated by Ramaswamy H. Sarma.

Drug Resistance Mechanism of M46I-Mutation-Induced Saquinavir Resistance in HIV-1 Protease Using Molecular Dynamics Simulation and Binding Energy Calculation

Drug-resistance-associated mutation in essential proteins of the viral life cycle is a major concern in anti-retroviral therapy. M46I, a non-active site mutation in HIV-1 protease has been clinically associated with saquinavir resistance in HIV patients. A 100 ns molecular dynamics (MD) simulation and MM-PBSA calculations were performed to study the molecular mechanism of M46I-mutation-based saquinavir resistance. In order to acquire deeper insight into the drug-resistance mechanism, the flap curling, closed/semi-open/open conformations, and active site compactness were studied. The M46I mutation significantly affects the energetics and conformational stability of HIV-1 protease in terms of RMSD, RMSF, Rg, SASA, and hydrogen formation potential. This mutation significantly decreased van der Waals interaction and binding free energy (∆G) in the M46I–saquinavir complex and induced inward flap curling and a wider opening of the flaps for most of the MD simulation period. The predominant open conformation was reduced, but inward flap curling/active site compactness was increased in the presence of saquinavir in M46I HIV-1 protease. In conclusion, the M46I mutation induced structural dynamics changes that weaken the protease grip on saquinavir without distorting the active site of the protein. The produced information may be utilized for the discovery of inhibitor(s) against drug-resistant HIV-1 protease.

Investigating potency of TMC-126 against wild-type and mutant variants of HIV-1 protease: a molecular dynamics and free energy study

A detailed computational study was performed to investigate the conformational changes of flap region and the mechanism underlying the binding of the inhibitor TMC-126 to HIV-1 protease (PR1) and its mutant variants through molecular dynamics simulations in conjunction with the molecular mechanics Poisson-Boltzmann (MM-PBSA) free energy calculation. Further, we have studied the effectiveness of the inhibitor against HIV-2 protease (PR2). The MM-PBSA calculation suggests that TMC-126 loses its potency against mutant variants and PR2 compared to wild-type PR1 mainly due to the loss in intermolecular electrostatic interactions. The potency of the inhibitor decreases in the order: wild type PR1 > M46L > MDR20 > I50V > PR2 > V32I > A28S. Our study reveals that the flap of PR1 adopts a semi-open conformation due to the mutation I50V or MDR20. The dissimilar nature of the movement of the flap tip of both monomers is evident from the dynamic cross-correlation map. The protein structural network analysis displays that mutation causes structural rearrangements and changes the communication path between residues. Overall, we believe our study may help explore and accelerate the development of novel HIV-1/HIV-2 protease inhibitors with better potency.

Revertant mutation V48G alters conformational dynamics of highly drug resistant HIV protease PRS17

Drug resistance is a serious problem for controlling the HIV/AIDS pandemic. Current antiviral drugs show several orders of magnitude worse inhibition of highly resistant clinical variant PRS17 of HIV-1 protease compared with wild-type protease. We have analyzed the effects of a common resistance mutation G48V in the flexible flaps of the protease by assessing the revertant PRS17_V48G for changes in enzyme kinetics, inhibition, structure, and dynamics. Both PRS17 and the revertant showed about 10-fold poorer catalytic efficiency than wild-type enzyme (0.55 and 0.39 μM^-1min^-1 compared to 6.3 μM^-1min^-1). Clinical inhibitors, amprenavir and darunavir, showed 2-fold and 8-fold better inhibition, respectively, of the revertant than of PRS17, although the inhibition constants for PRS17_V48G were still 25 to 1,200-fold worse than for wild-type protease. Crystal structures of inhibitor-free revertant and amprenavir complexes with revertant and PRS17 were solved at 1.3-1.5 Å resolution. The amprenavir complexes of PRS17_V48G and PRS17 showed no significant differences in the interactions with inhibitor, although changes were observed in the conformation of Phe53 and the interactions of the flaps. The inhibitor-free structure of the revertant showed flaps in an open conformation, however, the flap tips do not have the unusual curled conformation seen in inhibitor-free PRS17. Molecular dynamics simulations were run for 1 μs on the two inhibitor-free mutants and wild-type protease. PRS17 exhibited higher conformational fluctuations than the revertant, while the wild-type protease adopted the closed conformation and showed the least variation. The second half of the simulations captured the transition of the flaps of PRS17 from a closed to a semi-open state, whereas the flaps of PRS17_V48G tucked into the active site and the wild-type protease retained the closed conformation. These results suggest that mutation G48V contributes to drug resistance by altering the conformational dynamics of the flaps.

Single-Agent and Fixed-Dose Combination HIV-1 Protease Inhibitor Drugs in Fission Yeast (Schizosaccharomyces pombe)

Successful combination antiretroviral therapies (cART) eliminate active replicating HIV-1, slow down disease progression, and prolong lives. However, cART effectiveness could be compromised by the emergence of viral multidrug resistance, suggesting the need for new drug discoveries. The objective of this study was to further demonstrate the utility of the fission yeast cell-based systems that we developed previously for the discovery and testing of HIV protease (PR) inhibitors (PIs) against wild-type or multi-PI drug resistant _M11PR that we isolated from an infected individual. All thirteen FDA-approved single-agent and fixed-dose combination HIV PI drugs were tested. The effect of these drugs on HIV PR activities was tested in pure compounds or formulation drugs. All FDA-approved PI drugs, except for a prodrug FPV, were able to suppress the wild-type PR-induced cellular and enzymatic activities. Relative drug potencies measured by EC₅₀ in fission yeast were discussed in comparison with those measured in human cells. In contrast, none of the FDA-approved drugs suppressed the multi-PI drug resistant _M11PR activities. Results of this study show that fission yeast is a reliable cell-based system for the discovery and testing of HIV PIs and further demonstrate the need for new PI drugs against viral multi-PI resistance.

An insight on the nature of biochemical interactions between glycyrrhizin, myricetin and CYP3A4 isoform

Drug interaction studies are imperative to gain insights into the beneficial or harmful effects of therapeutic and dietary agents. This study investigated the mechanism of modulatory roles of glycyrrhizin (GLH) and myricetin (MYC) on the human CYP3A4 isoform using in silico and in vitro methods. While MYC had concentration-dependent inhibitory effect on CYP3A4 (IC₅₀ : 10.5 ± 0.55 μM) with characteristic K_m and V_max values of 1.13 μM and 1.54 nM/min, respectively, GLH exhibited no inhibitory effect on CYP3A4 activity in vitro. These observations are consistent with the results of in silico evaluations where the effect of MYC compared well with that of ketoconazole (a known CYP3A4 inhibitor) against CYP3A4. Overall, the established interactions between the study compounds and CYP3A4 could potentiate clinically vital drug-drug interactions and has lent credence to the mechanism of modulatory effect of MYC and GLH on CYP3A4 that could guide their safe use as therapeutic agents. PRACTICAL IMPLICATIONS: Myricetin (MYR) and glycyrrhizin (GLH) occur freely in commonly ingested foods and their supplements are recommended for the treatment of several debilitating diseases such as diabetes, cancer, and cardiovascular complications. This study provided an insight on the possible interactions that could be established between these compounds (MYR and GLH) and CYP3A4 when ingested and metabolized by the liver. The results suggested possibilities of potential clinical drug-drug interactions and advocates for their cautious use within the therapeutic dose in food supplements or medications to avoid probable liver damage.

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.

Acquired HIV-1 Protease Conformational Flexibility Associated with Lopinavir Failure May Shape the Outcome of Darunavir Therapy after Antiretroviral Therapy Switch

Understanding the underlying molecular interaction during a therapy switch from lopinavir (LPV) to darunavir (DRV) is essential to achieve long-term virological suppression. We investigated the kinetic and structural characteristics of multidrug-resistant South African HIV-1 subtype C protease (HIV-1 PR) during therapy switch from LPV to DRV using enzyme activity and inhibition assay, fluorescence spectroscopy, and molecular dynamic simulation. The HIV-1 protease variants were from clinical isolates with a combination of drug resistance mutations; MUT-1 (M46I, I54V, V82A, and L10F), MUT-2 (M46I, I54V, L76V, V82A, L10F, and L33F), and MUT-3 (M46I, I54V, L76V, V82A, L90M, and F53L). Enzyme kinetics analysis shows an association between increased relative resistance to LPV and DRV with the progressive decrease in the mutant HIV-1 PR variants' catalytic efficiency. A direct relationship between high-level resistance to LPV and intermediate resistance to DRV with intrinsic changes in the three-dimensional structure of the mutant HIV-1 PR as a function of the multidrug-resistance mutation was observed. In silico analysis attributed these structural adjustments to the multidrug-resistance mutations affecting the LPV and DRV binding landscape. Though DRV showed superiority to LPV, as a lower concentration was needed to inhibit the HIV-1 PR variants, the inherent structural changes resulting from mutations selected during LPV therapy may dynamically shape the DRV treatment outcome after the therapy switch.

Development of A Fission Yeast Cell-Based Platform for High Throughput Screening of HIV-1 Protease Inhibitors

BACKGROUND

HIV-1 protease inhibitor (PI) is one of the most potent classes of drugs in combinational antiretroviral therapies (cART). When a PI is used in combination with other anti- HIV drugs, cART can often suppress HIV-1 below detection thus prolonging the patient's lives. However, the challenge often faced by patients is the emergence of HIV-1 drug resistance. Thus, PIs with high genetic-barrier to drug-resistance are needed.

OBJECTIVE

The objective of this study was to develop a novel and simple fission yeast (Schizosaccharomyces pombe) cell-based system that is suitable for high throughput screening (HTS) of small molecules against HIV-1 protease (PR).

METHODS

A fission yeast RE294-GFP strain that stably expresses HIV-1 PR and green fluorescence protein (GFP) under the control of an inducible nmt1 promoter was used. Production of HIV-1 PR induces cellular growth arrest, which was used as the primary endpoint for the search of PIs and was quantified by an absorbance-based method. Levels of GFP production were used as a counter-screen control to eliminate potential transcriptional nmt1 inhibitors.

RESULTS

Both the absorbance-based HIV-1 PR assay and the GFP-based fluorescence assay were miniaturized and optimized for HTS. A pilot study was performed using a small drug library mixed with known PI drugs and nmt1 inhibitors. With empirically adjusted and clearly defined double-selection criteria, we were able to correctly identify the PIs and to exclude all hidden nmt1 inhibitors.

CONCLUSION

We have successfully developed and validated a fission yeast cell-based HTS platform for the future screening and testing of HIV-1 PR inhibitors.

Molecular dynamic mechanism(s) of inhibition of bioactive antiviral phytochemical compounds targeting cytochrome P450 3A4 and P-glycoprotein

P-glycoprotein (ABCB1) and cytochrome P450 3A4 (CYP3A4) metabolize almost all known human immunodeficiency virus' protease inhibitor drugs (PIs). Over induction of these proteins' activities has been linked to rapid metabolism of PIs which are then pumped out of the circulatory system, eventually leading to drug-resistance in HIV-positive patients. This study aims to determine, with the use of computational tools, the inhibitory potential of four phytochemical compounds (PCs) (epigallocatechin gallate (EGCG), kaempferol-7-glucoside (K7G), luteolin (LUT) and ellagic acid (EGA)) in inhibiting the activities of these drug-metabolizing proteins. The comparative analysis of the MM/GBSA results revealed that the binding affinity (ΔG_bind) of EGCG and K7G for CYP3A4 and ABCB1 are higher than LUT and EGA and fall between the ΔG_bind of the inhibitors of CYP3A4 and ABCB1 (Ritonavir (strong inhibitor) and Lopinavir (moderate inhibitor)). The structural analysis (RMSD, RMSF, RoG and protein-ligand interaction plots) also confirmed that EGCG and K7G showed similar inhibitory activities with the inhibitors. The study has shown that EGCG and K7G have inhibitory activities against the two proteins and assumes they could decrease intracellular efflux of PIs, consequently increasing the optimal concentration of PIs in the systemic circulation.Communicated by Ramaswamy H. Sarma.

Highly drug-resistant HIV-1 protease reveals decreased intra-subunit interactions due to clusters of mutations

Drug-resistance is a serious problem for treatment of the HIV/AIDS pandemic. Potent clinical inhibitors of HIV-1 protease show several orders of magnitude worse inhibition of highly drug-resistant variants. Hence, the structure and enzyme activities were analyzed for HIV protease mutant HIV-1 protease (EC 3.4.23.16) (PR) with 22 mutations (PRS5B) from a clinical isolate that was selected by machine learning to represent high-level drug-resistance. PRS5B has 22 mutations including only one (I84V) in the inhibitor binding site; however, clinical inhibitors had poor inhibition of PRS5B activity with kinetic inhibition value (K_i ) values of 4-1000 nm or 18- to 8000-fold worse than for wild-type PR. High-resolution crystal structures of PRS5B complexes with the best inhibitors, amprenavir (APV) and darunavir (DRV) (K_i ~ 4 nm), revealed only minor changes in protease-inhibitor interactions. Instead, two distinct clusters of mutations in distal regions induce coordinated conformational changes that decrease favorable internal interactions across the entire protein subunit. The largest structural rearrangements are described and compared to other characterized resistant mutants. In the protease hinge region, the N83D mutation eliminates a hydrogen bond connecting the hinge and core of the protease and increases disorder compared to highly resistant mutants PR with 17 mutations and PR with 20 mutations with similar hinge mutations. In a distal β-sheet, mutations G73T and A71V coordinate with accessory mutations to bring about shifts that propagate throughout the subunit. Molecular dynamics simulations of ligand-free dimers show differences consistent with loss of interactions in mutant compared to wild-type PR. Clusters of mutations exhibit both coordinated and antagonistic effects, suggesting PRS5B may represent an intermediate stage in the evolution of more highly resistant variants. DATABASES: Structural data are available in Protein Data Bank under the accession codes 6P9A and 6P9B for PRS5B/DRV and PRS5B/APV, respectively.

Potent antiviral HIV-1 protease inhibitor combats highly drug resistant mutant PR20

Drug-resistance threatens effective treatment of HIV/AIDS. Clinical inhibitors, including darunavir (1), are ineffective for highly resistant protease mutant PR20, however, antiviral compound 2 derived from 1 with fused tricyclic group at P2, extended amino-benzothiazole P2' ligand and two fluorine atoms on P1 shows 16-fold better inhibition of PR20 enzyme activity. Crystal structures of PR20 and wild-type PR complexes reveal how the extra groups of 2 counteract the expanded ligand-binding pocket, dynamic flaps, and faster dimer dissociation of PR20.

Decoding HIV resistance: from genotype to therapy

Genetic variation in HIV poses a major challenge for prevention and treatment of the AIDS pandemic. Resistance occurs by mutations in the target proteins that lower affinity for the drug or alter the protein dynamics, thereby enabling viral replication in the presence of the drug. Due to the prevalence of drug-resistant strains, monitoring the genotype of the infecting virus is recommended. Computational approaches for predicting resistance from genotype data and guiding therapy are discussed. Many prediction methods rely on rules derived from known resistance-associated mutations, however, statistical or machine learning can improve the classification accuracy and assess unknown mutations. Adding classifiers such as information on the atomic structure of the protein can further enhance the predictions.

A fission yeast cell-based system for multidrug resistant HIV-1 proteases

Background

HIV-1 protease (PR) is an essential enzyme for viral production. Thus, PR inhibitors (PIs) are the most effective class of anti-HIV drugs. However, the main challenge to the successful use of PI drugs in patient treatment is the emergence of multidrug resistant PRs (_mdrPRs). This study aimed to develop a fission yeast cell-based system for rapid testing of new PIs that combat _mdrPRs.

Results

Three _mdrPRs were isolated from HIV-infected patients that carried seven (_M7PR), ten (_M10PR) and eleven (_M11PR) PR gene mutations, respectively. They were cloned and expressed in fission yeast under an inducible promoter to allow the measurement of PR-specific proteolysis and drug resistance. The results showed that all three _mdrPRs maintained their abilities to proteolyze HIV viral substrates (MA↓CA and p6) and to confer drug resistance. Production of these proteins in the fission yeast caused cell growth inhibition, oxidative stress and altered mitochondrial morphologies that led to cell death. Five investigational PIs were used to test the utility of the established yeast system with an FDA-approved PI drug Darunavir (DRV) as control. All six compounds suppressed the wildtype PR (_wtPR) and the _M7PR-mediated activities. However, none of them were able to suppress the _M10PR or the _M11PR.

Conclusions

The three clinically isolated _mdrPRs maintained their viral proteolytic activities and drug resistance in the fission yeast. Furthermore, those viral _mdrPR activities were coupled with the induction of growth inhibition and cell death, which could be used to test the PI activities. Indeed, the five investigational PIs and DRV suppressed the _wtPR in fission yeast as they did in mammalian cells. Significantly, two of the high level _mdrPRs (_M10PR and _M11PR) were resistant to all of the existing PI drugs including DRV. This observation underscores the importance of continued searching for new PIs against _mdrPRs.

Electronic supplementary material

The online version of this article (doi:10.1186/s13578-016-0131-5) contains supplementary material, which is available to authorized users.

Structural Studies of a Rationally Selected Multi-Drug Resistant HIV-1 Protease Reveal Synergistic Effect of Distal Mutations on Flap Dynamics

We report structural analysis of HIV protease variant PR^S17 which was rationally selected by machine learning to represent wide classes of highly drug-resistant variants. Crystal structures were solved of PR^S17 in the inhibitor-free form and in complex with antiviral inhibitor, darunavir. Despite its 17 mutations, PR^S17 has only one mutation (V82S) in the inhibitor/substrate binding cavity, yet exhibits high resistance to all clinical inhibitors. PR^S17 has none of the major mutations (I47V, I50V, I54ML, L76V and I84V) associated with darunavir resistance, but has 10,000-fold weaker binding affinity relative to the wild type PR. Comparable binding affinity of 8000-fold weaker than PR is seen for drug resistant mutant PR20, which bears 3 mutations associated with major resistance to darunavir (I47V, I54L and I84V). Inhibitor-free PR^S17 shows an open flap conformation with a curled tip correlating with G48V flap mutation. NMR studies on inactive PR^S17_D25N unambiguously confirm that the flaps adopt mainly an open conformation in solution very similar to that in the inhibitor-free crystal structure. In PR^S17, the hinge loop cluster of mutations, E35D, M36I and S37D, contributes to the altered flap dynamics by a mechanism similar to that of PR20. An additional K20R mutation anchors an altered conformation of the hinge loop. Flap mutations M46L and G48V in PR^S17/DRV complex alter the Phe53 conformation by steric hindrance between the side chains. Unlike the L10F mutation in PR20, L10I in PR^S17 does not break the inter-subunit ion pair or diminish the dimer stability, consistent with a very low dimer dissociation constant comparable to that of wild type PR. Distal mutations A71V, L90M and I93L propagate alterations to the catalytic site of PR^S17. PR^S17 exhibits a molecular mechanism whereby mutations act synergistically to alter the flap dynamics resulting in significantly weaker binding yet maintaining active site contacts with darunavir.

A review on architecture of the gag-pol ribosomal frameshifting RNA in human immunodeficiency virus: a variability survey of virus genotypes

Programmed '-1' ribosomal frameshifting is necessary for expressing the pol gene overlapped from a gag of human immunodeficiency virus. A viral RNA structure that requires base pairing across the overlapping sequence region suggests a mechanism of regulating ribosome and helicase traffic during expression. To get precise roles of an element around the frameshift site, a review on architecture of the frameshifting RNA is performed in combination of reported information with augments of a representative set of 19 viral samples. In spite of a different length for the viral RNAs, a canonical comparison on the element sequence allocation is performed for viewing variability associations between virus genotypes. Additionally, recent and historical insights recognized in frameshifting regulation are looked back as for indel and single nucleotide polymorphism of RNA. As specially noted, structural changes at a frameshift site, the spacer sequence, and a three-helix junction element, as well as two Watson-Crick base pairs near a bulge and a C-G pair close a loop, are the most vital strategies for the virus frameshifting regulations. All of structural changes, which are dependent upon specific sequence variations, facilitate an elucidation about the RNA element conformation-dependent mechanism for frameshifting. These facts on disrupting base pair interactions also allow solving the problem of competition between ribosome and helicase on a same RNA template, common to single-stranded RNA viruses. In a broad perspective, each new insight of frameshifting regulation in the competition systems introduced by the RNA element construct changes will offer a compelling target for antiviral therapy.

Conformational variation of an extreme drug resistant mutant of HIV protease

Molecular mechanisms leading to high level drug resistance have been analyzed for the clinical variant of HIV-1 protease bearing 20 mutations (PR20); which has several orders of magnitude worse affinity for tested drugs. Two crystal structures of ligand-free PR20 with the D25N mutation of the catalytic aspartate (PR20D25N) revealed three dimers with different flap conformations. The diverse conformations of PR20D25N included a dimer with one flap in a unique "tucked" conformation; directed into the active site. Analysis of molecular dynamics (MD) simulations of the ligand-free PR20 and wild-type enzymes showed that the mutations in PR20 alter the correlated interactions between two monomers in the dimer. The two flaps tend to fluctuate more independently in PR20 than in the wild type enzyme. Combining the results of structural analysis by X-ray crystallography and MD simulations; unusual flap conformations and weakly correlated inter-subunit motions may contribute to the high level resistance of PR20.

Insights on the structural perturbations in human MTHFR Ala222Val mutant by protein modeling and molecular dynamics

Methylenetetrahydrofolate reductase (MTHFR) protein catalyzes the only biochemical reaction which produces methyltetrahydrofolate, the active form of folic acid essential for several molecular functions. The Ala222Val polymorphism of human MTHFR encodes a thermolabile protein associated with increased risk of neural tube defects and cardiovascular disease. Experimental studies have shown that the mutation does not affect the kinetic properties of MTHFR, but inactivates the protein by increasing flavin adenine dinucleotide (FAD) loss. The lack of completely solved crystal structure of MTHFR is an impediment in understanding the structural perturbations caused by the Ala222Val mutation; computational modeling provides a suitable alternative. The three-dimensional structure of human MTHFR protein was obtained through homology modeling, by taking the MTHFR structures from Escherichia coli and Thermus thermophilus as templates. Subsequently, the modeled structure was docked with FAD using Glide, which revealed a very good binding affinity, authenticated by a Glide XP score of -10.3983 (kcal mol(-1)). The MTHFR was mutated by changing Alanine 222 to Valine. The wild-type MTHFR-FAD complex and the Ala222Val mutant MTHFR-FAD complex were subjected to molecular dynamics simulation over 50 ns period. The average difference in backbone root mean square deviation (RMSD) between wild and mutant variant was found to be ~.11 Å. The greater degree of fluctuations in the mutant protein translates to increased conformational stability as a result of mutation. The FAD-binding ability of the mutant MTHFR was also found to be significantly lowered as a result of decreased protein grip caused by increased conformational flexibility. The study provides insights into the Ala222Val mutation of human MTHFR that induces major conformational changes in the tertiary structure, causing a significant reduction in the FAD-binding affinity.