Non-active site mutants of HIV-1 protease influence resistance and sensitisation towards protease inhibitors

Growth of infants delivered by mothers with HIV in Guangxi, China: An 18-month longitudinal follow-up study, 2015-2021

OBJECTIVES

The prevention of mother-to-child transmission of HIV has been a global success. But little is known about the growth parameters of infants delivered by mothers with HIV or the drug resistance of infants with HIV in China. The study aimed to assess growth parameters and drug resistance in Chinese infants exposed to HIV.

METHODS

We conducted an 18-month longitudinal follow-up study of 3283 infants (3222 without HIV; 61 with HIV) born to mothers with HIV in the Guangxi Zhuang Autonomous Region between January 2015 and December 2021. The weight and length of all participants was recorded. In addition, genetic subtypes and drug resistance analysis were performed for infants with HIV.

RESULTS

Compared with infants without HIV, those with HIV had significantly lower weight/length Z-scores, except at 18 months of age. The length/age Z-scores of infants with HIV was significantly reduced, except at 1 month of age. The weight/age Z-scores of infants with HIV were significantly lower at all follow-up time points. The weight/length Z-scores of male infants without HIV were significantly lower than for female infants without HIV at all follow-up time points. Male infants without HIV had lower length/age and weight/age Z-scores than female infants at the remaining follow-up points, except at 1 month of age. Of a total of 61 infants with HIV, subtype and drug-resistance data were obtained from 37 (60.66%) samples. Infants with HIV were dominated by the CRF01_AE genotype and showed a diversity of mutation sites dominated by non-nucleoside reverse transcriptase inhibitor resistance.

CONCLUSION

Our study demonstrates the growth of infants exposed to HIV in southwest China and provides detailed information on subtype distribution and drug resistance of those with HIV. Nutritional support and drug-resistance surveillance for infants exposed to HIV need to be strengthened.

Rational design, synthesis and biological evaluation of novel HIV-1 protease inhibitors containing 2-phenylacetamide derivatives as P2 ligands with potent activity against DRV-Resistant HIV-1 variants

A novel kind of potent HIV-1 protease inhibitors, containing diverse hydroxyphenylacetic acids as the P2-ligands and 4-substituted phenyl sulfonamides as the P2' ligands, were designed, synthesized and evaluated in this work. Majority of the target compounds exhibited good to excellent activity against HIV-1 protease with IC50 values below 200 nM. In particular, compound 18d with a 2-(3,4-dihydroxyphenyl) acetamide as the P2 ligand and a 4- methoxybenzene sulfonamide P2' ligand exhibited inhibitory activity IC50 value of 0.54 nM, which was better than that of the positive control darunavir (DRV). More importantly, no significant decline of the potency against HIV-1DRVRS (DRV-resistant mutation) and HIV-1NL4_3 variant (wild type) for 18d was detected. The molecular docking study of 18d with HIV-1 protease (PDB-ID: 1T3R, www.rcsb.org) revealed possible binding mode with the HIV-1 protease. These results suggested the validity of introducing phenol-derived moieties into the P2 ligand and deserve further optimization which was of great value for future discovery of novel HIV-1 protease.

Mechanism of drug resistance in HIV-1 protease subtype C in the presence of Atazanavir

AIDS is one of the deadliest diseases in the history of humankind caused by HIV. Despite the technological development, curtailing the viral infection inside human host still remains a challenge. Therapies such as HAART uses a combination of drugs to inhibit the viral activity. One of the important targets includes HIV protease and inhibiting its activity will minimize the production of mature structural proteins. However, the genetic diversity and the occurrence of drug resistant mutations adds complexity to effective drug design. In this study, we aimed at understanding the drug binding mechanism of one such subtype, namely subtype C and its insertion variant L38HL. We performed multiple molecular dynamics simulations along with binding free energy analysis of wild-type and L38HL bound to Atazanavir (ATV). From the analysis, we revealed that the insertion alters the hydrogen bond and hydrophobic interaction networks. The alterations in the interaction networks increase flexibility at the hinge-fulcrum interface. Further, the effects of these changes affect flap tip curling. Moreover, the changes in the hinge-fulcrum-cantilever interface alters the concerted motion of the functional regions leading to change in the direction of flap movement thus causing a subtle change in the active site volume. Additionally, formation of intramolecular hydrogen bonds in the ATV docked to L38HL restricted the movement of R1 and R2 groups thereby altering the interactions. Overall, the changes in the flexibility of flap together with the changes in the active site volume and compactness of the ligand provide insights for increased binding affinity of ATV with L38HL.

Viral Evolution Shaped by Host Proteostasis Networks

Understanding the factors that shape viral evolution is critical for developing effective antiviral strategies, accurately predicting viral evolution, and preventing pandemics. One fundamental determinant of viral evolution is the interplay between viral protein biophysics and the host machineries that regulate protein folding and quality control. Most adaptive mutations in viruses are biophysically deleterious, resulting in a viral protein product with folding defects. In cells, protein folding is assisted by a dynamic system of chaperones and quality control processes known as the proteostasis network. Host proteostasis networks can determine the fates of viral proteins with biophysical defects, either by assisting with folding or by targeting them for degradation. In this review, we discuss and analyze new discoveries revealing that host proteostasis factors can profoundly shape the sequence space accessible to evolving viral proteins. We also discuss the many opportunities for research progress proffered by the proteostasis perspective on viral evolution and adaptation.

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.

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.

Chemoinformatic Design and Profiling of Derivatives of Dasabuvir, Efavirenz, and Tipranavir as Potential Inhibitors of Zika Virus RNA-Dependent RNA Polymerase and Methyltransferase

Zika virus (ZIKV) infection is one of the mosquito-borne flaviviruses of human importance with more than 2 million suspected cases and more than 1 million people infected in about 30 countries. There are reported inhibitors of the zika virus replication machinery, but no approved effective antiviral therapy including vaccines directed against the virus for treatment or prevention is currently available. The study investigated the chemoinformatic design and profiling of derivatives of dasabuvir, efavirenz, and tipranavir as potential inhibitors of the zika virus RNA-dependent RNA polymerase (RdRP) and/or methyltransferase (MTase). The three-dimensional (3D) coordinates of dasabuvir, efavirenz, and tipranavir were obtained from the PubChem database, and their respective derivatives were designed with DataWarrior-5.2.1 using an evolutionary algorithm. Derivatives that were not mutagenic, tumorigenic, or irritant were selected; docked into RdRP and MTase; and further subjected to absorption, distribution, metabolism, excretion, and toxicity (ADMET) evaluation with Swiss-ADME and pkCSM web tools. Some of the designed compounds are Lipinski's rule-of-five compliant, with good synthetic accessibilities. Compounds 20d, 21d, 22d, and 1e are nontoxic with the only limitation of CYP1A2, CYP2C19, and/or CYP2C9 inhibition. Replacements of -CH₃ and -NH- in the methanesulfonamide moiety of dasabuvir with -OH and -CH₂- or -CH₂CH₂-, respectively, improved the safety/toxicity profile. Hepatotoxicity in 5d, 4d, and 18d is likely due to -NH- in their methanesulfonamide/sulfamic acid moieties. These compounds are potent inhibitors of N-7 and 2'-methylation activities of ZIKV methyltransferase and/or RNA synthesis through interactions with amino acid residues in the priming loop/"N-pocket" in the virus RdRP. Synthesis of these compounds and wet laboratory validation against ZIKV are recommended.

Identification of Aloe-derived natural products as prospective lead scaffolds for SARS-CoV-2 main protease (Mpro) inhibitors

In the past two years, the COVID-19 pandemic has caused over 5 million deaths and 250 million infections worldwide. Despite successful vaccination efforts and emergency approval of small molecule therapies, a diverse range of antivirals is still needed to combat the inevitable resistance that will arise from new SARS-CoV-2 variants. The main protease of SARS-CoV-2 (M^pro) is an attractive drug target due to the clinical success of protease inhibitors against other viruses, such as HIV and HCV. However, in order to combat resistance, various chemical scaffolds need to be identified that have the potential to be developed into potent inhibitors. To this end, we screened a high-content protease inhibitor library against M^proin vitro, in order to identify structurally diverse compounds that could be further developed into antiviral leads. Our high-content screening efforts retrieved 27 hits each with > 50% inhibition in our M^pro FRET assay. Of these, four of the top inhibitor compounds were chosen for follow-up due to their potency and drugability (Lipinski's rules of five criteria): anacardic acid, aloesin, aloeresin D, and TCID. Further analysis via dose response curves revealed IC₅₀ values of 6.8 μM, 38.9 μM, 125.3 μM, and 138.0 μM for each compound, respectively. Molecular docking studies demonstrated that the four inhibitors bound at the catalytic active site of M^pro with varying binding energies (-7.5 to -5.6 kcal/mol). Furthermore, M^pro FRET assay kinetic studies demonstrated that M^pro catalysis is better represented by a sigmoidal Hill model than the standard Michaelis-Menten hyperbola, indicating substantial cooperativity of the active enzyme dimer. This result suggests that the dimerization interface could be an attractive target for allosteric inhibitors. In conclusion, we identified two closely-related natural product compounds from the Aloe plant (aloesin and aloeresin D) that may serve as novel scaffolds for M^pro inhibitor design and additionally confirmed the strongly cooperative kinetics of M^pro proteolysis. These results further advance our knowledge of structure-function relationships in M^pro and offer new molecular scaffolds for inhibitor design.

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.

Pre-Exascale Computing of Protein-Ligand Binding Free Energies with Open Source Software for Drug Design

Nowadays, drug design projects benefit from highly accurate protein-ligand binding free energy predictions based on molecular dynamics simulations. While such calculations have been computationally expensive in the past, we now demonstrate that workflows built on open source software packages can efficiently leverage pre-exascale computing resources to screen hundreds of compounds in a matter of days. We report our results of free energy calculations on a large set of pharmaceutically relevant targets assembled to reflect industrial drug discovery projects.

One Plus One Makes Three: Triangular Coupling of Correlated Amino Acid Mutations

Correlated mutations have played a pivotal role in the recent success in protein fold prediction. Understanding nonadditive effects of mutations is crucial for altering protein structure, as mutations of multiple residues may change protein stability or binding affinity in a manner unforeseen by the investigation of single mutants. While the couplings between amino acids can be inferred from homologous protein sequences, the physical mechanisms underlying these correlations remain elusive. In this work we demonstrate that calculations based on the first-principles of statistical mechanics are capable of capturing the effects of nonadditivities in protein mutations. The identified thermodynamic couplings cover the short-range as well as previously unknown long-range correlations. We further explore a set of mutations in staphyloccocal nuclease to unravel an intricate interaction pathway underlying the correlations between amino acid mutations.

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.