Dynamical Network of HIV-1 Protease Mutants Reveals the Mechanism of Drug Resistance and Unhindered Activity

How do the mutations in PfK13 protein promote anti-malarial drug resistance?

Plasmodium falciparum develops resistance to artemisinin upon exposure to the anti-malarial drug. Various mutations in the Plasmodium falciparum Kelch13 (PfK13) protein such as Y493H, R539T, I543T and C580Y have been associated with anti-malarial drug resistance. These mutations impede the regular ubiquitination process that eventually invokes drug resistance. However, the relationship between the mutation and the mechanism of drug resistance has not yet been fully elucidated. The comparative protein dynamics are studied by performing the classical molecular dynamics (MD) simulations and subsequent analysis of the trajectories adopting root-mean-square fluctuations, the secondary-structure predictions and the dynamical cross-correlation matrix analysis tools. Here, we observed that the mutations in the Kelch-domain do not have any structural impact on the mutated site; however, it significantly alters the overall dynamics of the protein. The loop-region of the BTB-domain especially for Y493H and C580Y mutants is found to have the enhanced dynamical fluctuations. The enhanced fluctuations in the BTB-domain could affect the protein-protein (PfK13-Cullin) binding interactions in the ubiquitination process and eventually lead to anti-malarial drug resistance.Communicated by Ramaswamy H. Sarma.

Turning up the heat mimics allosteric signaling in imidazole-glycerol phosphate synthase

Allosteric drugs have the potential to revolutionize biomedicine due to their enhanced selectivity and protection against overdosage. However, we need to better understand allosteric mechanisms in order to fully harness their potential in drug discovery. In this study, molecular dynamics simulations and nuclear magnetic resonance spectroscopy are used to investigate how increases in temperature affect allostery in imidazole glycerol phosphate synthase. Results demonstrate that temperature increase triggers a cascade of local amino acid-to-amino acid dynamics that remarkably resembles the allosteric activation that takes place upon effector binding. The differences in the allosteric response elicited by temperature increase as opposed to effector binding are conditional to the alterations of collective motions induced by either mode of activation. This work provides an atomistic picture of temperature-dependent allostery, which could be harnessed to more precisely control enzyme function.

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.

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.

Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease

Computational analysis of protein–ligand interactions is of crucial importance for drug discovery. Assessment of ligand binding energy allows us to have a glimpse of the potential of a small organic molecule to be a ligand to the binding site of a protein target. Available scoring functions, such as in docking programs, all rely on equations that sum each type of protein–ligand interactions in order to predict the binding affinity. Most of the scoring functions consider electrostatic interactions involving the protein and the ligand. Electrostatic interactions constitute one of the most important part of total interactions between macromolecules. Unlike dispersion forces, they are highly directional and therefore dominate the nature of molecular packing in crystals and in biological complexes and contribute significantly to differences in inhibition strength among related enzyme inhibitors. In this study, complexes of HIV-1 protease with inhibitor molecules (JE-2147 and darunavir) were analyzed by using charge densities from the transferable aspherical-atom University at Buffalo Databank (UBDB). Moreover, we analyzed the electrostatic interaction energy for an ensemble of structures, using molecular dynamic simulations to highlight the main features of electrostatic interactions important for binding affinity.

Understanding the mechanism of HIV-1 protease inhibition by monoclonal antibodies

HIV-1 protease is an essential enzyme in the life cycle of human immunodeficiency virus (HIV) and hence is one of the most important targets for antiviral drug design. Although there are ten FDA approved drugs against HIV protease (PR), their long term usage elicits mutations leading to drug resistance. As a result, novel therapeutic approaches are being explored including synthetic antibodies. Recently, a murine monoclonal antibody, mAB1696 (mAB) was reported to inhibit PR by preventing dimerization. Crystallographic data could reveal only six protease residues that interact with mAB. The present study employs a range of computational techniques, starting from protein-protein docking to all-atomic molecular dynamics simulations to generate plausible 3D structures of PR-mAB complex. Results show that mAB interacts very strongly with several PR dimer interface residues, such as Gln7, Arg8 (N-terminal), Cys95, Leu97 (C-terminal), Thr26, Gly27 (active site), Gly49, Ile50 (flap), apart from its interactions with the PR epitope region, Pro1-Trp6 (N-terminal). These observations support the hypothesis that binding of mAB prevents the dimerization of PR. The interactions and binding conformations identified in this study could form the basis for designing allosteric inhibitors preventing the dimerization of HIV-1 Protease.

Disulfide Reduction Allosterically Destabilizes the β-Ladder Subdomain Assembly within the NS1 Dimer of ZIKV

The Zika virus (ZIKV) was responsible for a recent debilitating epidemic that till date has no cure. A potential way to reduce ZIKV virulence is to limit the action of the nonstructural proteins involved in its viral replication. One such protein, NS1, encoded as a monomer by the viral genome, plays a major role via symmetric oligomerization. We examine the homodimeric structure of the dominant β-ladder segment of NS1 with extensive all atom molecular dynamics. We find it stably bounded by two spatially separated interaction clusters (C1 and C2) with significant differences in the nature of their interactions. Four pairs of distal, intramonomeric disulfide bonds are found to be coupled to the stability, local structure, and wettability of the interfacial region. Symmetric reduction of the intramonomeric disulfides triggers marked dynamical heterogeneity, interfacial wettability, and asymmetric salt-bridging propensity. Harnessing the model-free Lipari-Szabo based formalism for estimation of conformational entropy (S_conf), we find clear signatures of heterogeneity in the monomeric conformational entropies. The observed asymmetry, very small in the unperturbed state, expands significantly in the reduced states. This allosteric effect is most noticeable in the electrostatically bound C2 cluster that underlies the greatest stability in the unperturbed state. Allosteric induction of conformational and thermodynamic asymmetry is expected to affect the pathways leading to symmetric higher-ordered oligomerization, and thereby affect crucial replication pathways.

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.

Molecular Basis of Differential Stability and Temperature Sensitivity of ZIKA versus Dengue Virus Protein Shells

Rapid spread of ZIKA virus (ZIKV) and its association with severe birth defects have raised worldwide concern. Recent studies have shown that ZIKV retains its infectivity and remains structurally stable at temperatures up to 40 °C, unlike dengue and other flaviviruses. In spite of recent cryo-EM structures that showed similar architecture of ZIKA and dengue virus (DENV) E protein shells, little is known that makes ZIKV so temperature insensitive. Here, we attempt to unravel the molecular basis of greater thermal stability of ZIKV over DENV2 by executing atomistic molecular dynamics (MD) simulations on the viral E protein shells at 37 °C. Our results suggest that ZIKA E protein shell retains its structural integrity through stronger inter-raft communications facilitated by a series of electrostatic and H-bonding interactions among multiple inter-raft residues. In comparison, the DENV2 E protein shell surface was loosly packed that exhibited holes at all 3-fold vertices, in close agreement with another EM structure solved at 37 °C. The residue-level information obtained from our study could pave way for designing small molecule inhibitors and specific antibodies to inhibit ZIKV E protein assembly and membrane fusion.

The Determination of HIV-1 RT Mutation Rate, Its Possible Allosteric Effects, and Its Implications on Drug Resistance

The high mutation rate of the human immunodeficiency virus type 1 (HIV-1) plays a major role in treatment resistance, from the development of vaccines to therapeutic drugs. In addressing the crux of the issue, various attempts to estimate the mutation rate of HIV-1 resulted in a large range of 10⁻⁵–10⁻³ errors/bp/cycle due to the use of different types of investigation methods. In this review, we discuss the different assay methods, their findings on the mutation rates of HIV-1 and how the locations of mutations can be further analyzed for their allosteric effects to allow for new inhibitor designs. Given that HIV is one of the fastest mutating viruses, it serves as a good model for the comprehensive study of viral mutations that can give rise to a more horizontal understanding towards overall viral drug resistance as well as emerging viral diseases.

Exploring the potency of currently used drugs against HIV-1 protease of subtype D variant by using multiscale simulations

Acquired immune deficiency syndrome (AIDS) is caused by the human immunodeficiency virus (HIV), type 1 and 2. Further, the diversity in HIV-1 has given rise to many serotypes and recombinant strains. The currently used protease inhibitors have been developed for subtype B, although non-B subtype strains account for ∼ 90% of the global HIV infections. Subtype D is spreading rapidly and infecting a large population in North Africa and the Middle East. In the current study, molecular dynamics simulations in conjunction with the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) scheme was used to investigate the potency of four drugs, namely atazanavir (ATV), darunavir (DRV), lopinavir (LPV) and tipranavir (TPV) against the subtype D variant. Our calculations predicted that the potency of the inhibitors decreased in the order TPV > ATV > DRV > LPV. TPV was found to be the most potent against subtype D due to an increase in van der Waals and electrostatic interactions and reduction in the desolvation energy compared to other inhibitors. This result is further supported by the hydrogen bond interactions between inhibitors and protease. Furthermore, our analyses suggested that the binding of TPV induced a more closed conformation of the flap compared to apo or other complexes. It was observed that TPV/PR^D has a lower cavity volume relative to the other three complexes leading to a tighter binding. The open conformation of the flap was observed for LPV/PR^D. We expect that this study might be useful for designing more potent inhibitors against HIV-1 subtype D. Communicated by Ramaswamy H. Sarma.

Reviewing HIV-1 Gag Mutations in Protease Inhibitors Resistance: Insights for Possible Novel Gag Inhibitor Designs

HIV protease inhibitors against the viral protease are often hampered by drug resistance mutations in protease and in the viral substrate Gag. To overcome this drug resistance and inhibit viral maturation, targeting Gag alongside protease rather than targeting protease alone may be more efficient. In order to successfully inhibit Gag, understanding of its drug resistance mutations and the elicited structural changes on protease binding needs to be investigated. While mutations on Gag have already been mapped to protease inhibitor resistance, there remain many mutations, particularly the non-cleavage mutations, that are not characterized. Through structural studies to unravel how Gag mutations contributes to protease drug resistance synergistically, it is thus possible to glean insights to design novel Gag inhibitors. In this review, we discuss the structural role of both novel and previously reported Gag mutations in PI resistance, and how new Gag inhibitors can be designed.

Eigenvector centrality for characterization of protein allosteric pathways

Determining the principal energy-transfer pathways responsible for allosteric communication in biomolecules remains challenging, partially due to the intrinsic complexity of the systems and the lack of effective characterization methods. In this work, we introduce the eigenvector centrality metric based on mutual information to elucidate allosteric mechanisms that regulate enzymatic activity. Moreover, we propose a strategy to characterize the range of correlations that underlie the allosteric processes. We use the V-type allosteric enzyme imidazole glycerol phosphate synthase (IGPS) to test the proposed methodology. The eigenvector centrality method identifies key amino acid residues of IGPS with high susceptibility to effector binding. The findings are validated by solution NMR measurements yielding important biological insights, including direct experimental evidence for interdomain motion, the central role played by helix h[Formula: see text], and the short-range nature of correlations responsible for the allosteric mechanism. Beyond insights on IGPS allosteric pathways and the nature of residues that could be targeted by therapeutic drugs or site-directed mutagenesis, the reported findings demonstrate the eigenvector centrality analysis as a general cost-effective methodology to gain fundamental understanding of allosteric mechanisms at the molecular level.

Elucidating the Interdependence of Drug Resistance from Combinations of Mutations

HIV-1 protease is responsible for the cleavage of 12 nonhomologous sites within the Gag and Gag-Pro-Pol polyproteins in the viral genome. Under the selective pressure of protease inhibition, the virus evolves mutations within (primary) and outside of (secondary) the active site, allowing the protease to process substrates while simultaneously countering inhibition. The primary protease mutations impede inhibitor binding directly, while the secondary mutations are considered accessory mutations that compensate for a loss in fitness. However, the role of secondary mutations in conferring drug resistance remains a largely unresolved topic. We have shown previously that mutations distal to the active site are able to perturb binding of darunavir (DRV) via the protein's internal hydrogen-bonding network. In this study, we show that mutations distal to the active site, regardless of context, can play an interdependent role in drug resistance. Applying eigenvalue decomposition to collections of hydrogen bonding and van der Waals interactions from a series of molecular dynamics simulations of 15 diverse HIV-1 protease variants, we identify sites in the protease where amino acid substitutions lead to perturbations in nonbonded interactions with DRV and/or the hydrogen-bonding network of the protease itself. While primary mutations are known to drive resistance in HIV-1 protease, these findings delineate the significant contributions of accessory mutations to resistance. Identifying the variable positions in the protease that have the greatest impact on drug resistance may aid in future structure-based design of inhibitors.

Evolution of inhibitor-resistant natural mutant forms of HIV-1 protease probed by pre-steady state kinetic analysis

Pre-steady state kinetic analysis of mechanistic features of substrate binding and processing is crucial for insight into the evolution of inhibitor-resistant forms of HIV-1 protease. These data may provide a correct vector for rational drug design assuming possible intrinsic dynamic effects. These data should also give some clues to the molecular mechanism of protease action and resistance to inhibitors. Here we report pre-steady state kinetics of the interaction of wild type or mutant forms of HIV-1 protease with a FRET-labeled peptide. The three-stage "minimal" kinetic scheme with first and second reversible steps of substrate binding and with following irreversible peptide cleavage step adequately described experimental data. For the first time, a set of "elementary" kinetic parameters of wild type HIV-1 protease and its natural mutant inhibitor-resistant forms MDR-HM, ANAM-11 and prDRV4 were compared. Inhibitors of the first and second generation were used to estimate the inhibitory effects on HIV-1 protease activity. The resulting set of kinetic data supported that the mutant forms are kinetically unaffected by inhibitors of the first generation, proving their functional resistance to these compounds. The second generation inhibitor darunavir inhibited mutant forms MDR-HM and ANAM-11, but was ineffective against prDRV4. Our kinetic data revealed that these inhibitors induced different conformational changes in the enzyme and, thereby they have different mode of binding in the enzyme active site. These data confirmed hypothesis that the driving force of the inhibitor-resistance evolution is disruption of enzyme-inhibitor complex by changing of the contact network in the inhibitor binding site.

The effect of V155M mutation on the complex of hSTING and 2′3′-cGAMP: an in silico study case

All atoms MD simulations for the complexes of WT/V155M hSTING and 2′3′-cGAMP were performed to study the effect of the V155M mutation on SAVI disease at the molecular level.

Structural analyses of 2015-updated drug-resistant mutations in HIV-1 protease: an implication of protease inhibitor cross-resistance

Background

Strategies to control HIV for improving the quality of patient lives have been aided by the Highly Active Anti-Retroviral Therapy (HAART), which consists of a cocktail of inhibitors targeting key viral enzymes. Numerous new drugs have been developed over the past few decades but viral resistances to these drugs in the targeted viral enzymes are increasingly reported. Nonetheless the acquired mutations often reduce viral fitness and infectivity. Viral compensatory secondary-line mutations mitigate this loss of fitness, equipping the virus with a broad spectrum of resistance against these drugs. While structural understanding of the viral protease and its drug resistance mutations have been well established, the interconnectivity and development of structural cross-resistance remain unclear. This paper reports the structural analyses of recent clinical mutations on the drug cross-resistance effects from various protease and protease inhibitors (PIs) complexes.

Methods

Using the 2015 updated clinical HIV protease mutations, we constructed a structure-based correlation network and a minimum-spanning tree (MST) based on the following features: (i) topology of the PI-binding pocket, (ii) allosteric effects of the mutations, and (iii) protease structural stability.

Results and conclusion

Analyis of the network and the MST of dominant mutations conferring resistance to the seven PIs (Atazanavir-ATV, Darunavir-DRV, Indinavir-IDV, Lopinavir-LPV, Nelfinavir-NFV, Saquinavir-SQV, and Tipranavir-TPV) showed that cross-resistance can develop easily across NFV, SQV, LPV, IDV, and DRV, but not for ATV or TPV. Through estimation of the changes in vibrational entropies caused by each reported mutation, some secondary mutations were found to destabilize protease structure. Our findings provide an insight into the mechanism of PI cross-resistance and may also be useful in guiding the selection of PI in clinical treatment to delay the onset of cross drug resistance.

Electronic supplementary material

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

Large-Scale Structure-Based Prediction and Identification of Novel Protease Substrates Using Computational Protein Design

Characterizing the substrate specificity of protease enzymes is critical for illuminating the molecular basis of their diverse and complex roles in a wide array of biological processes. Rapid and accurate prediction of their extended substrate specificity would also aid in the design of custom proteases capable of selectively and controllably cleaving biotechnologically or therapeutically relevant targets. However, current in silico approaches for protease specificity prediction, rely on, and are therefore limited by, machine learning of sequence patterns in known experimental data. Here, we describe a general approach for predicting peptidase substrates de novo using protein structure modeling and biophysical evaluation of enzyme-substrate complexes. We construct atomic resolution models of thousands of candidate substrate-enzyme complexes for each of five model proteases belonging to the four major protease mechanistic classes-serine, cysteine, aspartyl, and metallo-proteases-and develop a discriminatory scoring function using enzyme design modules from Rosetta and AMBER's MMPBSA. We rank putative substrates based on calculated interaction energy with a modeled near-attack conformation of the enzyme active site. We show that the energetic patterns obtained from these simulations can be used to robustly rank and classify known cleaved and uncleaved peptides and that these structural-energetic patterns have greater discriminatory power compared to purely sequence-based statistical inference. Combining sequence and energetic patterns using machine-learning algorithms further improves classification performance, and analysis of structural models provides physical insight into the structural basis for the observed specificities. We further tested the predictive capability of the model by designing and experimentally characterizing the cleavage of four novel substrate motifs for the hepatitis C virus NS3/4 protease using an in vivo assay. The presented structure-based approach is generalizable to other protease enzymes with known or modeled structures, and complements existing experimental methods for specificity determination.

Biophysical and computational methods to analyze amino acid interaction networks in proteins

Globular proteins are held together by interacting networks of amino acid residues. A number of different structural and computational methods have been developed to interrogate these amino acid networks. In this review, we describe some of these methods, including analyses of X-ray crystallographic data and structures, computer simulations, NMR data, and covariation among protein sequences, and indicate the critical insights that such methods provide into protein function. This information can be leveraged towards the design of new allosteric drugs, and the engineering of new protein function and protein regulation strategies.