Exploring molecular mechanism of allosteric inhibitor to relieve drug resistance of multiple mutations in HIV-1 protease by enhanced conformational sampling

Molecular basis for differential recognition of an allosteric inhibitor by receptor tyrosine kinases

Understanding kinase-inhibitor selectivity continues to be a major objective in kinase drug discovery. We probe the molecular basis of selectivity of an allosteric inhibitor (MSC1609119A-1) of the insulin-like growth factor-I receptor kinase (IGF1RK), which has been shown to be ineffective for the homologous insulin receptor kinase (IRK). Specifically, we investigated the structural and energetic basis of the allosteric binding of this inhibitor to each kinase by combining molecular modeling, molecular dynamics (MD) simulations, and thermodynamic calculations. We predict the inhibitor conformation in the binding pocket of IRK and highlight that the charged residues in the histidine-arginine-aspartic acid (HRD) and aspartic acid-phenylalanine-glycine (DFG) motifs and the nonpolar residues in the binding pocket govern inhibitor interactions in the allosteric pocket of each kinase. We suggest that the conformational changes in the IGF1RK residues M1054 and M1079, movement of the ⍺C-helix, and the conformational stabilization of the DFG motif favor the selectivity of the inhibitor toward IGF1RK. Our thermodynamic calculations reveal that the observed selectivity can be rationalized through differences observed in the electrostatic interaction energy of the inhibitor in each inhibitor/kinase complex and the hydrogen bonding interactions of the inhibitor with the residue V1063 in IGF1RK that are not attained with the corresponding residue V1060 in IRK. Overall, our study provides a rationale for the molecular basis of recognition of this allosteric inhibitor by IGF1RK and IRK, which is potentially useful in developing novel inhibitors with improved affinity and selectivity.

An Exploration of the Inhibitory Mechanism of Rationally Screened Benzofuran-1,3,4-Oxadiazoles and-1,2,4-Triazoles as Inhibitors of NS5B RdRp Hepatitis C Virus through Pharmacoinformatic Approaches

Benzofuran, 1,3,4-oxadiazole, and 1,2,4-triazole are privileged heterocyclic moieties that display the most promising and wide spectrum of biological activities against a wide variety of diseases. In the current study, benzofuran-1,3,4-oxadiazole BF1-BF7 and benzofuran-1,2,4-triazole compounds BF8-BF15 were tested against HCV NS5B RNA-dependent RNA polymerase (RdRp) utilizing structure-based screening via a computer-aided drug design (CADD) approach. A molecular docking approach was applied to evaluate the binding potential of benzofuran-appended 1,3,4-oxadiazole and 1,2,4-triazole BF1-BF15 molecules. Benzofuran-1,3,4-oxadiazole scaffolds BF1-BF7 showed lesser binding affinities (-12.63 to -14.04 Kcal/mol) than benzofuran-1,2,4-triazole scaffolds BF8-BF15 (-14.11 to -16.09 Kcal/mol) against the HCV NS5B enzyme. Molecular docking studies revealed the excellent binding affinity scores exhibited by benzofuran-1,2,4-triazole structural motifs BF-9 (-16.09 Kcal/mol), BF-12 (-15.75 Kcal/mol), and BF-13 (-15.82 Kcal/mol), respectively, which were comparatively better than benzofuran-based HCV NS5B inhibitors' standard reference drug Nesbuvir (-15.42 Kcal/mol). A molecular dynamics simulation assay was also conducted to obtain valuable insights about the enzyme-compounds interaction profile and structural stability, which indicated the strong intermolecular energies of the BF-9+NS5B complex and the BF-12+NS5B complex as per the MM-PBSA method, while the BF-12+NS5B complex was the most stable system as per the MM-GBSA calculation. The drug-likeness and ADMET studies of all the benzofuran-1,2,4-triazole derivatives BF8-BF15 revealed that these compounds possessed good medicinal chemistry profiles in agreement with all the evaluated parameters for being drugs. The molecular docking affinity scores, MM-PBSA/MM-GBSA and MD-simulation stability analysis, drug-likeness profiling, and ADMET study assessment indicated that N-4-fluorophenyl-S-linked benzofuran-1,2,4-triazole BF-12 could be a future promising anti-HCV NS5B RdRp inhibitor therapeutic drug candidate that has a structural agreement with the Nesbuvir standard reference drug.

The design of compounds with desirable properties - The anti-HIV case study

Efficacy and safety are among the most desirable characteristics of an ideal drug. The tremendous increase in computing power and the entry of artificial intelligence into the field of computational drug design are accelerating the process of identifying, developing, and optimizing potential drugs. Here, we present novel approach to design new molecules with desired properties. We combined various neural networks and linear regression algorithms to build models for cytotoxicity and anti-HIV activity based on Continual Molecular Interior analysis (CoMIn) and Cinderella's Shoe (CiS) derived molecular descriptors. After validating the reliability of the models, a genetic algorithm was coupled with the Des-Pot Grid algorithm to generate new molecules from a predefined pool of molecular fragments and predict their bioactivity and cytotoxicity. This combination led to the proposal of 16 hit molecules with high anti-HIV activity and low cytotoxicity. The anti-SARS-CoV-2 activity of the hits was predicted.

Decoding drug resistant mechanism of V32I, I50V and I84V mutations of HIV-1 protease on amprenavir binding by using molecular dynamics simulations and MM-GBSA calculations

Mutations V32I, I50V and I84V in the HIV-1 protease (PR) induce drug resistance towards drug amprenavir (APV). Multiple short molecular dynamics (MSMD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) method were utilized to investigate drug-resistant mechanism of V32I, I50V and I84V towards APV. Dynamic information arising from MSMD simulations suggest that V32I, I50V and I84V highly affect structural flexibility, motion modes and conformational behaviours of two flaps in the PR. Binding free energies calculated by MM-GBSA method suggest that the decrease in binding enthalpy and the increase in binding entropy induced by mutations V32I, I50V and I84V are responsible for drug resistance of the mutated PRs on APV. The energetic contributions of separate residues on binding of APV to the PR show that V32I, I50V and I84V highly disturb the interactions of two flaps with APV and mostly drive the decrease in binding ability of APV to the PR. Thus, the conformational changes of two flaps in the PR caused by V32I, I50V and I84V play key roles in drug resistance of three mutated PR towards APV. This study can provide useful dynamics information for the design of potent inhibitors relieving drug resistance.

The structural, dynamic, and thermodynamic basis of darunavir resistance of a heavily mutated HIV-1 protease using molecular dynamics simulation

The human immunodeficiency virus type 1 protease (HIV-1 PR) is an important enzyme in the life cycle of the HIV virus. It cleaves inactive pre-proteins of the virus and changes them into active proteins. Darunavir (DRV) suppresses the wild-type HIV-1 PR (WT-Pr) activity but cannot inhibit some mutant resistant forms (MUT-Pr). Increasing knowledge about the resistance mechanism can be helpful for designing more effective inhibitors. In this study, the mechanism of resistance of a highly MUT-Pr strain against DRV was investigated. For this purpose, complexes of DRV with WT-Pr (WT-Pr-D) and MUT-Pr (MUT-Pr-D) were studied by all-atom molecular dynamics simulation in order to extract the dynamic and energetic properties. Our data revealed that mutations increased the flap-tip flexibility due to the reduction of the flap-flap hydrophobic interactions. So, the protease’s conformation changed from a closed state to a semi-open state that can facilitate the disjunction of DRV from the active site. On the other hand, energy analysis limited to the final basins of the energy landscape indicated that the entropy of binding of DRV to MUT-Pr was more favorable than that of WT-Pr. However, the enthalpy penalty overcomes it and makes binding more unfavorable relative to the WT-Pr. The unfavorable interaction of DRV with R8, I50, I84, D25′, and A28′ residues in MUT-Pr-D relative to WT-Pr-D is the reason for this enthalpy penalty. Thus, mutations drive resistance to DRV. The hydrogen bond analysis showed that compared with WT-Pr, the hydrogen bonds between DRV and the active-site residues of MUT-Pr were disrupted.

Q61 mutant-mediated dynamics changes of the GTP-KRAS complex probed by Gaussian accelerated molecular dynamics and free energy landscapes

Understanding the molecular mechanism of the GTP-KRAS binding is significant for improving the target roles of KRAS in cancer treatment. In this work, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations were applied to decode the effect of Q61A, Q61H and Q61L on the activity of KRAS. Dynamics analyses based on MR-GaMD trajectory reveal that motion modes and dynamics behavior of the switch domain in KRAS are heavily affected by the three Q61 mutants. Information of free energy landscapes (FELs) shows that Q61A, Q61H and Q61L induce structural disorder of the switch domain and disturb the activity of KRAS. Analysis of the interaction network uncovers that the decrease in the stability of hydrogen bonding interactions (HBIs) of GTP with residues V29 and D30 induced by Q61A, Q61H and Q61L is responsible for the structural disorder of the switch-I and that in the occupancy of the hydrogen bond between GTP and residue G60 leads to the structural disorder of the switch-II. Thus, the high disorder of the switch domain caused by three current Q61 mutants produces a significant effect on binding of KRAS to its effectors. This work is expected to provide useful information for further understanding function and target roles of KRAS in anti-cancer drug development.

Multiple Molecular Dynamics Simulations and Free-Energy Predictions Uncover the Susceptibility of Variants of HIV-1 Protease against Inhibitors Darunavir and KNI-1657

HIV-1 protease (PR) is considered to be the main targets of anti-AIDS drug design because of its role in the proteolytic processing of viral polyproteins. However, the emergence of drug-resistant HIV has become a major problem in the therapy of HIV-1-infected patients. Focused on the complexes of wild type (WT) PR and two mutant PRs (V32I/L33F/I54M/V82I and V32I/L33F/I54M/I84 V) with inhibitors Darunavir (DRV) and KNI-1657 (KNI), respectively, we have conducted research on the conformational dynamics and the resistance mechanism caused by residue mutations through multiple molecular dynamics (MD) simulations combined with an energy (MM-PBSA and solvated interaction energy (SIE)) prediction. The results indicate that mutated residues of PR alter the distance between flap regions and catalytic sites, the volume of the inner catalytic site, and the curling degree of the flap tips, thereby affecting DRV and KNI inhibitor binding to PR. These mutated residues reduced the binding affinity of the two mutant PRs to DRV, resulting in drug resistance, whereas the two mutant PRs increase the binding affinity with KNI, indicating they enhance the sensitivity to KNI. Compared with the WT PR, the changes in van der Waals interaction and electrostatic interaction in the two variant PRs play a vital part in the binding of PR with DRV and KNI. These results may supply valuable guidance for the design of anti-AIDS drugs targeting PR.

Decoding molecular mechanism underlying binding of drugs to HIV-1 protease with molecular dynamics simulations and MM-GBSA calculations

HIV-1 protease (PR) is thought to be efficient targets of anti-AIDS drug design. Molecular dynamics (MD) simulations and multiple post-processing analysis technologies were applied to decipher molecular mechanism underlying binding of three drugs Lopinavir (LPV), Nelfinavir (NFV) and Atazanavir (ATV) to the PR. Binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) suggest that compensation between binding enthalpy and entropy plays a vital role in binding of drugs to PR. Dynamics analyses show that binding of LPV, NFV and ATV highly affects structural flexibility, motion modes and dynamics behaviour of the PR, especially for two flaps. Computational alanine scanning and interaction network analysis verify that although three drugs have structural difference, they share similar binding modes to the PR and common interaction clusters with the PR. The current findings also confirm that residues located interaction clusters, such as Asp25/Asp25', Gly27/Gly27', Ala28/Ala28', Asp29, Ile47/Ile47', Gly49/Gly49', Ile50/Ile50', Val82/Val82' and Ile84/Ile84, can be used as efficient targets of clinically available inhibitors towards the PR.

Insights into effect of the Asp25/Asp25' protonation states on binding of inhibitors Amprenavir and MKP97 to HIV-1 protease using molecular dynamics simulations and MM-GBSA calculations

The protonation states of two aspartic acids in the catalytic strands of HIV-1 protease (PR) remarkably affect bindings of inhibitors to PR. It is requisite for the design of potent inhibitors towards PR to investigate the influences of Asp25/Asp25' protonated states on dynamics behaviour of PR and binding mechanism of inhibitors to PR. In this work, molecular dynamics (MD) simulations, MM-GBSA method and principal component (PC) analysis were coupled to explore the effect of Asp25/Asp25' protonation states on conformational changes of PR and bindings of Amprenavir and MKP97 to PR. The results show that the Asp25/Asp25' protonation states exert different impacts on structural fluctuations, flexibility and motion modes of PR. Dynamics analysis verifies that Asp25/Asp25' protonated states highly affect conformational dynamics of two flaps in PR. The binding free energy calculations results suggest that the Asp25/Asp25' protonated states obviously strengthen bindings of inhibitors to PR compared to the non-protonation state. Calculations of residue-based free energy decomposition indicate that the Asp25/Asp25' protonation not only disturbs the interaction network of inhibitors with PR but also stabilizes bindings of inhibitors to PR by cancelling the electrostatic repulsive interaction. Therefore, special attentions should be paid to the Asp25/Asp25' protonation in the design of potent inhibitors towards PR.

Genotypic Methods for HIV Drug Resistance Monitoring: The Opportunities and Challenges Faced by China

AIDS is a globalized infectious disease. In 2014, UNAIDS launched a global project of "90-90-90" to end the HIV epidemic by 2030. The second and third 90 require 90% of HIV-1 infected individuals receiving antiretroviral therapy (ART) and durable virological suppression. However, wide use of ART will greatly increase the emergence and spreading of HIV drug resistance and current HIV drug resistance test (DRT) assays in China are seriously lagging behind, hindering to achieve virological suppression. Therefore, recommending an appropriate HIV DRT method is critical for HIV routine surveillance and prevention in China. In this review, we summarized the current existing HIV drug resistance genotypic testing methods around the world and discussed the advantages and disadvantages of these methods.

Novel radial distribution function approach in the study of point mutations: the HIV-1 protease case study

Background: Mutations are one of the engines of evolution. Under constant stress pressure, mutations can lead to the emergence of unwanted, drug-resistant entities. Methodology: The radial distribution function weighted by the number of valence shell electrons is used to design quantitative structure-activity relationship (QSAR) model relating descriptors with the inhibition constant for a series of wild-type HIV-1 protease inhibitor complexes. The residuals of complexes with mutant HIV-1 protease were correlated with the energy of the highest occupied molecular orbitals of the residues introduced to enzyme via point mutations. Conclusion: Successful identification of residues Ile3, Asp25, Val32 and Ile50 as the one whose substitution influences the inhibition constant the most, demonstrates the potential of the proposed methodology for the study of the effects of point mutations.

Effect of mutations on binding of ligands to guanine riboswitch probed by free energy perturbation and molecular dynamics simulations

Riboswitches can regulate gene expression by direct and specific interactions with ligands and have recently attracted interest as potential drug targets for antibacterial. In this work, molecular dynamics (MD) simulations, free energy perturbation (FEP) and molecular mechanics generalized Born surface area (MM-GBSA) methods were integrated to probe the effect of mutations on the binding of ligands to guanine riboswitch (GR). The results not only show that binding free energies predicted by FEP and MM-GBSA obtain an excellent correlation, but also indicate that mutations involved in the current study can strengthen the binding affinity of ligands GR. Residue-based free energy decomposition was applied to compute ligand-nucleotide interactions and the results suggest that mutations highly affect interactions of ligands with key nucleotides U22, U51 and C74. Dynamics analyses based on MD trajectories indicate that mutations not only regulate the structural flexibility but also change the internal motion modes of GR, especially for the structures J12, J23 and J31, which implies that the aptamer domain activity of GR is extremely plastic and thus readily tunable by nucleotide mutations. This study is expected to provide useful molecular basis and dynamics information for the understanding of the function of GR and possibility as potential drug targets for antibacterial.

Structural and binding insights into HIV-1 protease and P2-ligand interactions through molecular dynamics simulations, binding free energy and principal component analysis

HIV-1 protease (HIV-1-pr) plays an important role in viral replication and maturation, making it one of the most attractive targets for anti-retroviral therapy. To design new effective inhibitors able to combat drug resistance in mutant HIV-1-pr variants, it is essential to gain further understanding about the mechanisms by which the recently proposed inhibitors deactivate the mutant HIV-1-pr variants. In the present work, we explored the interactions between two P2-ligands (DRV, and one new derivative, 4UY) with wild type (WT) and two multiple mutant HIV-1-pr variants (p20 and p51) with all atom molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA). The trajectories of MD simulations show that both 4UY and DRV primarily bind with the active sites, flap and 80s loop regions of HIV-1-pr variants through either hydrogen bonds or hydrophobic interactions. More hydrogen bonds and hydrophobic interactions were located for 4UY/HIV-1-pr complexes than for DRV/HIV-1-pr counterparts. More importantly, 4UY was found to have an extra hydrogen bond with the backbone of Gly48' in the flap region of the HIV-1-prs. The flap tip-tip distance (I50-I50') and flap tip-active site distance (I50-D25 and I50'-D25') indicate that the flaps turn more closed in 4UY bound HIV-1-prs than DRV bound ones, and the former also have more compact hydrophobic cavities than the latter. Further, the vector projections from PCA indicate that 4UY/DRV inhibitor binding projects the closing of flap in HIV-1-pr variants. In line with the above trajectory analysis, the thermodynamics calculation with MM-PBSA method suggests much stronger binding affinity for 4UY/HIV-1-pr than DRV/HIV-1-pr by 4.3-6.4 kcal/mol. Although p20 and p51 also induce weaker binding due to multiple mutants for 4UY inhibitor by 1.9-1.8 kcal/mol, their bindings to the new P2 ligand (4UY) are indeed significantly enhanced as compared to DRV. The thermodynamic components responsible for the binding differences and the contribution from key residues to the binding were also discussed in detail.

Decoding molecular mechanism of inhibitor bindings to CDK2 using molecular dynamics simulations and binding free energy calculations

CDK2 can be used as an attractive target for development of efficient inhibitors curing multiple disease relating with CDK2. In this work, molecular dynamics (MD) simulations and binding free energy calculations were coupled to probe conformational changes of CDK2 due to inhibitor associations and binding mechanisms of inhibitors PM1, FMD and X64 to CDK2. The results suggest that the binding strength of FMD and X64 to CDK2 is stronger than that of PM1. Principal component (PC) analysis and cross-correlation map calculations based on the equilibrated MD trajectories demonstrate that the structural difference in inhibitors exerts important impact on motion modes and dynamics behavior of CDK2. Residue-based free energy decomposition method was adopted to estimate the inhibitor-residue spectrum. The results not only efficiently identify the hot interaction spot of inhibitors with CDK2 but also show that the hydrophobic rings R1, R2 and R3 as well as polar groups of three inhibitors play key roles in favorably binding of inhibitors to CDK2. This work is expected to contribute energetic basis and dynamics information to development of promising inhibitors toward CDK2.Communicated by Ramaswamy H. Sarma.