Exploring the drug resistance of V32I and M46L mutant HIV-1 protease to inhibitor TMC114: flap dynamics and binding mechanism

Roles of the gate loop in β-arrestin-1 conformational dynamics and phosphorylated receptor interaction

Arrestins interact with phosphorylated G protein-coupled receptors (GPCRs) and regulate the homologous desensitization and internalization of GPCRs. The gate loop in arrestins is a critical region for both stabilization of the basal state and interaction with phosphorylated receptors. We investigated the roles of specific residues in the gate loop (K292, K294, and H295) using β-arrestin-1 and phosphorylated C-tail peptide of vasopressin receptor type 2 (V2Rpp) as a model system. We measured the binding affinity of V2Rpp and analyzed conformational dynamics of β-arrestin-1. Our results suggest that K294 plays a critical role in the interaction with V2Rpp without influencing the overall conformation of the V2Rpp-bound state. The residues K292 and H295 contribute to the stability of the polar core in the basal state and form a specific conformation of the finger loop in the V2Rpp-bound state.

Evaluation of the inhibitory potency of anti-dengue phytocompounds against DENV-2 NS2B-NS3 protease: virtual screening, ADMET profiling and molecular dynamics simulation investigations

Dengue fever has been a worldwide concern, with 50-100 million new infections each year mainly due to five different serotypes of the Dengue virus (DENV). Designing a perfect anti-dengue agent that can inhibit all the serotypes by distinguishing antigenic differences is quite difficult. Previous anti-dengue researches have included chemical compounds screening against DENV enzymes. The ongoing analysis is meant for investigation of the plant-based compounds as antagonistic to DENV-2 focusing on the specific NS2B-NS3Pro target, a trypsin like serine protease that cuts the DENV polyprotein into separate proteins crucial for viral reproduction. Initially, a virtual library of more than 130 phytocompounds was prepared from previously published reports of plants with anti-dengue properties, which were then virtually screened and shortlisted against the WT, H51N and S135A mutant of DENV-2 NS2B-NS3Pro. The three top-most compounds were viewed as Gallocatechin (GAL), Flavokawain-C (FLV), and Isorhamnetin (ISO) showing docking scores of -5.8, -5.7, -5.7 kcal/mol for WT, -7.5, -6.8, -7.6 kcal/mol for the H51N, and -6.9, -6.5, -6.1 kcal/mol for the S135A mutant protease, respectively. 100 ns long MD simulations and MM-GBSA based free energy calculations were performed on the NS2B-NS3Pro complexes to witness the relative binding affinity of the compounds and favourable molecular interactions network. A comprehensive analysis of the study reveals some promising outcomes with ISO as the topmost compound with favourable pharmacokinetic properties for the WT and mutants (H51N and S135A) as well, suggesting as a novel anti-NS2B-NS3Pro agent with better adapting characters in both the mutants.Communicated by Ramaswamy H. Sarma.

Ebola and Marburg virus VP35 coiled-coil validated as antiviral target by tripartite split-GFP complementation

Ebola virus (EBOV) and Marburg virus (MARV) are highly pathogenic viruses in humans, against which approved antivirals are lacking. During EBOV and MARV infection, coiled-coil mediated oligomerization is essential for the virion protein 35 (VP35) polymerase co-factor function and type I interferon antagonism, making VP35 coiled-coil an elective drug target. We established a tripartite split-green fluorescent protein (GFP) fluorescence complementation (FC) system based on recombinant GFP-tagged EBOV and MARV VP35, which probes VP35 coiled-coil assembly by monitoring fluorescence on E. coli colonies, or in vitro in 96/384-multiwell. Oligomerization-defective VP35 mutants showed that correct coiled-coil knobs-into-holes pairing within VP35 oligomer is pre-requisite for GFP tags and GFP detector to reconstitute fluorescing full-length GFP. The method was validated by screening a small compound library, which identified Myricetin and 4,5,6,7-Tetrabromobenzotriazole as inhibitors of EBOV and MARV VP35 oligomerization-dependent FC with low-micromolar IC₅₀ values. These findings substantiate the VP35 coiled-coil value as antiviral target.

•

Ebola and Marburg virus VP35 oligomerize via trimeric and tetrameric coiled-coil

•

VP35 coiled-coil assembly triggers fluorescence of a tripartite split-GFP system

•

Mutations perturbing VP35 coiled-coil hamper split-GFP complementation

•

Myricetin and TBBT inhibit split-GFP complementation mediated by VP35 coiled-coil

Targeting the DENV NS2B-NS3 protease with active antiviral phytocompounds: structure-based virtual screening, molecular docking and molecular dynamics simulation studies

Dengue fever has been a global health concern. Mitigation is a challenging problem due to non-availability of workable treatments. The most difficult objective is to design a perfect anti-dengue agent capable of inhibiting infections caused by all four serotypes. Various tactics have been employed in the past to discover dengue antivirals, including screening of chemical compounds against dengue virus enzymes. The objective of the current study is to investigate phytocompounds as anti-dengue remedies that target the non-structural 2B and non-structural 3 protease (NS2B-NS3^pro), a possible therapeutic target for dengue fever. Initially, 300 + antiviral phytocompounds were collected from Duke’s phytochemical and ethnobotanical database and 30 phytocompounds with anti-dengue properties were identified from previously reported studies, which were virtually screened against NS2B-NS3^pro using molecular docking and toxicity evaluation. The top five most screened ligands were naringin, hesperidin, gossypol, maslinic acid and rhodiolin with binding affinities of − 8.7 kcal/mol, − 8.5 kcal/mol, − 8.5 kcal/mol, − 8.5 kcal/mol and − 8.1 kcal/mol, respectively. The finest docked compounds complexed with NS2B-NS3^pro were subjected for molecular dynamics (MD) simulations and binding free energy estimations through molecular mechanics generalized born surface area–based calculations. The results of the study are intriguing in the context of computer-aided screening and the binding affinities of the phytocompounds, proposing maslinic acid (MAS) as a potent bioactive antiviral for the development of phytocompound-based anti-dengue agent.

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.

Molecular insights to the binding interactions of APNS containing HIV-protease inhibitors against SARS-CoV-2 Mpro: an in silico approach towards drug repurposing

SARS-CoV-2 M^pro is one of the most vital enzymes of the new coronavirus-2 (SARS-CoV-2) and is a crucial target for drug discovery. Unfortunately, there is not any potential drugs available to combat the action of SARS-CoV-2 M^pro. Based on the reports HIV-protease inhibitors can be applied against the SARS by targeting the SARS-CoV-1 M^pro, we have chosen few clinically trialed experimental and allophenylnorstatine (APNS) containing HIV-protease inhibitors (JE-2147, JE-533, KNI-227, KNI-272 & KNI-1931), to examine their binding affinities with SARS-CoV-2 M^pro and to assess their potential to check for a possible drug candidate against the protease. Here, we have chosen a methodology to understand the binding mechanism of these five inhibitors to SARS-CoV-2 M^pro by merging molecular docking, molecular dynamics (MD) simulation and MM-PBSA based free energy calculations. Our estimations disclose that JE-2147 is highly effective (ΔG_Bind = -28.31 kcal/mol) due to an increased favorable van der Waals (ΔE_vdw) interactions and decreased solvation (ΔG_solv) energies between the inhibitor and viral protease. JE-2147 shows a higher level of interactions as compared to JE-533 (-6.85 kcal/mol), KNI-227 (-18.36 kcal/mol), KNI-272 (-15.69 kcal/mol) and KNI-1931 (-21.59 kcal/mol) against SARS-CoV-2 M^pro. Binding contributions of important residues (His41, Met49, Cys145, His164, Met165, Glu166, Pro168, Gln189, etc.) from the active site or near the active site regions with ≥1.0 kcal/mol suggest a potent binding of the inhibitors. It is anticipated that the current study of binding interactions of these APNS containing inhibitors can pitch some valuable insights to design the significantly effective anti-SARS-CoV-2 M^pro drugs.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.

Interactions of the α3β2 Nicotinic Acetylcholine Receptor Interfaces with α-Conotoxin LsIA and its Carboxylated C-terminus Analogue: Molecular Dynamics Simulations

Notably, α-conotoxins with carboxy-terminal (C-terminal) amidation are inhibitors of the pentameric nicotinic acetylcholine receptors (nAChRs), which are therapeutic targets for neurological diseases and disorders. The (α3)₂(β2)₃ nAChR subunit arrangement comprises a pair of α3(+)β2(−) and β2(+)α3(−) interfaces, and a β2(+)β2(−) interface. The β2(+)β2(−) interface has been suggested to have higher agonist affinity relative to the α3(+)β2(−) and β2(+)α3(−) interfaces. Nevertheless, the interactions formed by these subunit interfaces with α-conotoxins are not well understood. Therefore, in order to address this, we modelled the interactions between α-conotoxin LsIA and the α3β2 subtype. The results suggest that the C-terminal carboxylation of LsIA predominantly influenced the enhanced contacts of the conotoxin via residues P7, P14 and C17 on LsIA at the α3(+)β2(−) and β2(+)α3(−) interfaces. However, this enhancement is subtle at the β2(+)β2(−) site, which can compensate the augmented interactions by LsIA at α3(+)β2(−) and β2(+)α3(−) binding sites. Therefore, the divergent interactions at the individual binding interface may account for the minor changes in binding affinity to α3β2 subtype by C-terminal carboxylation of LsIA versus its wild type, as shown in previous experimental results. Overall, these findings may facilitate the development of new drug leads or subtype-selective probes.

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.

Recent Developments and Applications of the MMPBSA Method

The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.

Unraveling HIV protease flaps dynamics by Constant pH Molecular Dynamics simulations

The active site of HIV protease (HIV-PR) is covered by two flaps. These flaps are known to be essential for the catalytic activity of the HIV-PR, but their exact conformations at the different stages of the enzymatic pathway remain subject to debate. Understanding the correct functional dynamics of the flaps might aid the development of new HIV-PR inhibitors. It is known that, the HIV-PR catalytic efficiency is pH-dependent, likely due to the influence of processes such as charge transfer and protonation/deprotonation of ionizable residues. Several Molecular Dynamics (MD) simulations have reported information about the HIV-PR flaps. However, in MD simulations the protonation of a residue is fixed and thus it is not possible to study the correlation between conformation and protonation state. To address this shortcoming, this work attempts to capture, through Constant pH Molecular Dynamics (CpHMD), the conformations of the apo, substrate-bound and inhibitor-bound HIV-PR, which differ drastically in their flap arrangements. The results show that the HIV-PR flaps conformations are defined by the protonation of the catalytic residues Asp25/Asp25' and that these residues are sensitive to pH changes. This study suggests that the catalytic aspartates can modulate the opening of the active site and substrate binding.

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.