Multiple routes and milestones in the folding of HIV-1 protease monomer

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.

Repurposing antiviral drugs against HTLV-1 protease by molecular docking and molecular dynamics simulation

Human T-cell leukemia virus type I (HTLV-1) belongs to the delta retrovirus family and the etiological agent of adult T-cell leukemia (ATL(. While the current HTLV-1 therapy, relies on using Zidovudine plus IFN-γ, there is no FDA approved drugs against it. In silico drug repurposing is a fast and accurate way for screening US-FDA approved drugs to find a therapeutic option for the HTLV-1 infection. So that, this research aims to analyze a dataset of approved antiviral drugs as a potential prospect for an anti-viral drug against HTLV-1 infection. Molecular docking simulation was performed to identify interactions of the antiviral drugs with the key residues in the HTLV-1 protease binding site. Then, molecular dynamics simulation was also performed for the potential protein-ligand complexes to confirm the stable behavior of the ligands inside the binding pocket. The best docking scores with the target was found to be Simeprevir, Atazanavir, and Saquinavir compounds which indicate that these drugs can firmly bind to the HTLV-1 protease. The MD simulation confirmed the stability of Simeprevir-protease, Atazanavir-Protease, and Saquinavir-Protease interactions. Clearly, these compounds should be further evaluated in experimental assays and clinical trials to confirm their actual activity against HTLV-1 infection.Communicated by Ramaswamy H. Sarma.

Mechanism of darunavir binding to monomeric HIV-1 protease: a step forward in the rational design of dimerization inhibitors

HIV protease (HIVPR) is a key target in AIDS therapeutics. All ten FDA-approved drugs that compete with substrates in binding to this dimeric enzyme's active site have become ineffective due to the emergence of drug resistant mutants. Blocking the dimerization interface of HIVPR is thus being explored as an alternate strategy. The latest drug, darunavir (DRV), which exhibited a high genetic barrier to viral resistance, is said to have a dual mode of action - (i) binding to the dimeric active site, and (ii) preventing the dimerization by binding to the HIVPR monomer. Despite several reports on DRV complexation with dimeric HIVPR, the mode and mechanism of the binding of DRV to the HIVPR monomer are poorly understood. In this study, we utilized all-atomic MD simulations and umbrella sampling techniques to identify the best possible binding mode of DRV to the monomeric HIVPR and its mechanism of association. The results suggest that DRV binds between the active site and the flap of the monomer, and the flap plays a crucial role in directing the drug to bind and driving the other protein domains to undergo induced fit changes for stronger complexation. The obtained binding mode of DRV was validated by comparing with various mutational data from clinical isolates to reported in vitro mutations. The identified binding pose was also able to successfully reproduce the experimental K_i value in the picomolar range. The residue-level information extracted from this study could accelerate the structure-based drug designing approaches targeting HIVPR dimerization.

HIV-1 Protease Dimerization Dynamics Reveals a Transient Druggable Binding Pocket at the Interface

The binding mechanism of HIV-1 protease monomers leading to the catalytically competent dimeric enzyme has been investigated by means of state-of-the-art atomistic simulations. The emerging picture allows a deeper understanding of experimental observations and reveals that water molecules trapped at the interface have an important role in slowing down the kinetics of the association process. Unexpectedly, a cryptic binding pocket is identified at the interface of the complex, corresponding to a partially bound dimer that lacks enzymatic function. The pocket has a transient nature with a lifetime longer than 1 μs, and it displays very favorable druggability features. Docking as well as MM-GBSA free-energy calculations further support the possibility to target the new binding site by means of inhibitors able to prevent the complete dimerization by capturing the inactive conformation. This discovery could open the way to the rational design of a new class of anti-HIV drugs.

Probing the structural dynamics of proteins and nucleic acids with optical tweezers

Conformational changes are an essential feature of most molecular processes in biology. Optical tweezers have emerged as a powerful tool for probing conformational dynamics at the single-molecule level because of their high resolution and sensitivity, opening new windows on phenomena ranging from folding and ligand binding to enzyme function, molecular machines, and protein aggregation. By measuring conformational changes induced in a molecule by forces applied by optical tweezers, new insight has been gained into the relationship between dynamics and function. We discuss recent advances from studies of how structure forms in proteins and RNA, including non-native structures, fluctuations in disordered proteins, and interactions with chaperones assisting native folding. We also review the development of assays probing the dynamics of complex protein-nucleic acid and protein-protein assemblies that reveal the dynamic interactions between biomolecular machines and their substrates.

The Parkinson disease-linked LRRK2 protein mutation I2020T stabilizes an active state conformation leading to increased kinase activity

The effect of leucine-rich repeat kinase 2 (LRRK2) mutation I2020T on its kinase activity has been controversial, with both increased and decreased effects being reported. We conducted steady-state and pre-steady-state kinetic studies on LRRKtide and its analog LRRKtide(S). Their phosphorylation differs by the rate-limiting steps: product release is rate-limiting for LRRKtide and phosphoryl transfer is rate-limiting for LRRKtide(S). As a result, we observed that the I2020T mutant is more active than wild type (WT) LRRK2 for LRRKtide(S) phosphorylation, whereas it is less active than WT for LRRKtide phosphorylation. Our pre-steady-state kinetic data suggest that (i) the I2020T mutant accelerates the rates of phosphoryl transfer of both reactions by 3-7-fold; (ii) this increase is masked by a rate-limiting product release step for LRRKtide phosphorylation; and (iii) the observed lower activity of the mutant for LRRKtide phosphorylation is a consequence of its instability: the concentration of the active form of the mutant is 3-fold lower than WT. The I2020T mutant has a dramatically low KATP and therefore leads to resistance to ATP competitive inhibitors. Two well known DFG-out or type II inhibitors are also weaker toward the mutant because they inhibit the mutant in an unexpected ATP competitive mechanism. The I2020 residue lies next to the DYG motif of the activation loop of the LRRK2 kinase domain. Our modeling and metadynamic simulations suggest that the I2020T mutant stabilizes the DYG-in active conformation and creates an unusual allosteric pocket that can bind type II inhibitors but in an ATP competitive fashion.

Natural polymorphisms and unusual mutations in HIV-1 protease with potential antiretroviral resistance: a bioinformatic analysis

Background

The correlations of genotypic and phenotypic tests with treatment, clinical history and the significance of mutations in viruses of HIV-infected patients are used to establish resistance mutations to protease inhibitors (PIs). Emerging mutations in human immunodeficiency virus type 1 (HIV-1) protease confer resistance to PIs by inducing structural changes at the ligand interaction site. The aim of this study was to establish an in silico structural relationship between natural HIV-1 polymorphisms and unusual HIV-1 mutations that confer resistance to PIs.

Results

Protease sequences isolated from 151 Mexican HIV-1 patients that were naïve to, or subjected to antiretroviral therapy, were examined. We identified 41 unrelated resistance mutations with a prevalence greater than 1%. Among these mutations, nine exhibited positive selection, three were natural polymorphisms (L63S/V/H) in a codon associated with drug resistance, and six were unusual mutations (L5F, D29V, L63R/G, P79L and T91V). The D29V mutation, with a prevalence of 1.32% in the studied population, was only found in patients treated with antiretroviral drugs. Using in silico modelling, we observed that D29V formed unstable protease complexes when were docked with lopinavir, saquinavir, darunavir, tipranavir, indinavir and atazanavir.

Conclusions

The structural correlation of natural polymorphisms and unusual mutations with drug resistance is useful for the identification of HIV-1 variants with potential resistance to PIs. The D29V mutation likely confers a selection advantage in viruses; however, in silico, presence of this mutation results in unstable enzyme/PI complexes, that possibly induce resistance to PIs.

The maturation of HIV-1 protease precursor studied by discrete molecular dynamics

The equilibrium properties of a HIV-1-protease precursor are studied by means of an efficient molecular dynamics scheme, which allows for the simulation of the folding of the protein monomers and their dimerization into an active form and compare them with those of the mature protein. The results of the model provide, with atomic detail, an overall account of several experimental findings, including the NMR conformation of the mature dimer, the calorimetric properties of the system, the effects of the precursor tail on the dimerization constant, the secondary chemical shifts of the monomer, and the paramagnetic relaxation enhancement data associated with the conformations of the precursor. It is found that although the mature protein can dimerize in a unique, single way, the precursor populates several dimeric conformations in which monomers are always native-like, but their binding can be non-native.

Type II kinase inhibitors show an unexpected inhibition mode against Parkinson's disease-linked LRRK2 mutant G2019S

A number of well-known type II inhibitors (ATP-noncompetitive) that bind kinases in their DFG-out conformation were tested against wild-type LRRK2 and the most common Parkinson's disease-linked mutation, G2019S. We found that traditional type II inhibitors exhibit surprising variability in their inhibition mechanism between the wild type (WT) and the G2019S mutant of LRRK2. The type II kinase inhibitors were found to work in an ATP-competitive fashion against the G2019S mutant, whereas they appear to follow the expected noncompetitive mechanism against WT. Because the G2019S mutation lies in the DXG motif (DYG in LRRK2 but DFG in most other kinases) of the activation loop, we explored the structural consequence of the mutation on loop dynamics using an enhanced sampling method called metadynamics. The simulations suggest that the G2019S mutation stabilizes the DYG-in state of LRRK2 through a series of hydrogen bonds, leading to an increase in the conformational barrier between the active and inactive forms of the enzyme and a relative stabilization of the active form. The conformational bias toward the active form of LRRK2 mutants has two primary consequences. (1) The mutant enzyme becomes hyperactive, a known contributor to the Parkinsonian phenotype, as a consequence of being "locked" into the activated state, and (2) the mutation creates an unusual allosteric pocket that can bind type II inhibitors but in an ATP-competitive fashion. Our results suggest that developing type II inhibitors, which are generally considered superior to type I inhibitors because of desirable selectivity profiles, might be especially challenging for the G2019S LRRK2 mutant.

How Does Darunavir Prevent HIV-1 Protease Dimerization?

The drug Darunavir (DRV) is a potent inhibitor of HIV-1 protease (PR), a homodimeric essential enzyme of the AIDS virus. Recent experimental data suggest that DRV is able to prevent dimerization of HIV-1 PR, which, together with its high affinity for the mature enzyme, has been linked to the high genetic barrier to the development of viral resistance. The mechanism of dimerization inhibition and the binding mode(s) of DRV to monomeric HIV-1 PR are unknown. Here, multiple molecular dynamics simulations with explicit solvent (for a total of 11 μs with the CHARMM force field and 1 μs with the Amber force field) show that the monomer of HIV-1 PR is structurally stable and reveal a major binding mode of DRV. This binding mode is stabilized by favorable interactions between the apolar groups of DRV and the hydrophobic residues Ile32, Ile47, Ile50, Ile54, Pro79, Val82, and Ile84. The binding mode to monomeric HIV-1 PR identified by molecular dynamics is different from the two binding modes observed in the crystal structure of the complex with dimeric HIV-1 PR. As an example, there are no interactions between DRV and the catalytic Asp25 in the binding mode to monomeric HIV-1 PR revelead by the simulations. In contrast, the simulations show extensive and stable interactions between DRV and the flap (residues 46-55), which are likely to sterically hinder the formation of the flap interface as observed in the dimeric structure. Which of the two mechanisms of inhibition (dimerization inhibition by association with the flap or binding to the active site of the mature enzyme) dominates might depend on the HIV-1 PR mutations, and it is likely that dimerization inhibition is predominant for multiple mutations at the active site in multidrug resistant strains.