An integrated molecular dynamics, principal component analysis and residue interaction network approach reveals the impact of M184V mutation on HIV reverse transcriptase resistance to lamivudine

Repurposing of SARS-CoV-2 compounds against Marburg Virus using MD simulation, mm/GBSA, PCA analysis, and free energy landscape

The significant mortality rate associated with Marburg virus infection made it the greatest hazard among infectious diseases. Drug repurposing using in silico methods has been crucial in identifying potential compounds that could prevent viral replication by targeting the virus's primary proteins. This study aimed at repurposing the drugs of SARS-CoV-2 for identifying potential candidates against the matrix protein VP40 of the Marburg virus. Virtual screening was performed where the control compound, Nilotinib, showed a binding score of -9.99 kcal/mol. Based on binding scores, hit compounds 9549298, 11960895, 44545852, 51039094, and 89670174 were selected that had a lower binding score than the control. Subsequent molecular dynamics (MD) simulation revealed that compound 9549298 consistently formed a hydrogen bond with the residue Gln290. This was observed both in molecular docking and MD simulation poses, indicating a strong and significant interaction with the protein. 11960895 had the most stable and consistent RMSD pattern exhibited in 100 ns simulation, while 9549298 had the most identical RMSD plot compared to the control molecule. MM/PBSA analysis showed that the binding free energy (ΔG) of 9549298 and 11960895 was lower than the control, with -30.84 and -38.86 kcal/mol, respectively. It was observed by the PCA (principal component analysis) and FEL (free energy landscape) analysis that compounds 9549298 and 11960895 had lesser conformational variation. Overall, this study proposed 9549298 and 11960895 as potential binders of VP40 MARV that can cause its inhibition, however it inherently lacks experimental validation. Furthermore, the study proposes in-vitro experiments as the next step to validate these computational findings, offering a practical approach to further explore these compounds' potential as antiviral agents.Communicated by Ramaswamy H. Sarma.

From Data to Knowledge: Systematic Review of Tools for Automatic Analysis of Molecular Dynamics Output

An increasing number of crystal structures available on one side, and the boost of computational power available for computer-aided drug design tasks on the other, have caused that the structure-based drug design tools are intensively used in the drug development pipelines. Docking and molecular dynamics simulations, key representatives of the structure-based approaches, provide detailed information about the potential interaction of a ligand with a target receptor. However, at the same time, they require a three-dimensional structure of a protein and a relatively high amount of computational resources. Nowadays, as both docking and molecular dynamics are much more extensively used, the amount of data output from these procedures is also growing. Therefore, there are also more and more approaches that facilitate the analysis and interpretation of the results of structure-based tools. In this review, we will comprehensively summarize approaches for handling molecular dynamics simulations output. It will cover both statistical and machine-learning-based tools, as well as various forms of depiction of molecular dynamics output.

Leveraging on Active Site Similarities; Identification of Potential Inhibitors of Zinc-Finger and UFSP domain Protein (ZUFSP)

BACKGROUND

ZUFSP (Zinc-finger and UFSP domain protein) is a novel representative member of the recently characterized seventh class of deubiquitinating enzymes (DUBs). Due to the roles DUBs play in genetic instability, they have become a major drug target in cancer and neurodegenerative diseases. ZUFSP, being a DUB enzyme has also been implicated in genetic stability. However, no lead compound has been developed to target ZUFSP.

OBJECTIVE/METHODS

Therefore, in this study, we used a combined drug repurposing, virtual screening and per-Residue Energy Decomposition (PRED) to identify ZUFSP inhibitors with therapeutic potential. 3-bromo-6-{[4-hydroxy-1-3(3-phenylbutanoyl)piperidin-4-yl]methyl}-4H,5H,6H,7H-thieno[2,3- C]pyridine-7-one (BHPTP) which is an inhibitor of USP7 was repurposed to target ZUFSP. The rationale behind this is based on the similarity of the active between USP7 and ZUFSP.

RESULTS

PRED of the binding between BHPTP and ZUFSP revealed Cys223, Arg408, Met410, Asn460, and Tyr465 as the crucial residues responsible for this interaction. The pharmacophoric moieties of BHPTP responsible for this binding along with other physiochemical properties were used as a filter to retrieve potential ligands. 799 compounds were retrieved, ZINC083241427, ZINC063648749, and ZINC063648753 were selected due to the binding energy they exhibited. Cheminformatics analysis revealed that the compounds possess high membrane permeability, however, BHPTP had a low membrane permeability. Furthermore, the compounds are drug like, having obeyed Lipinski's rule of five.

CONCLUSION

Taken together, findings from this study put ZINC083241427, ZINC063648749, and ZINC063648753 as potential ZUFSP inhibitor, however, more experimental validation is required to unravel the mechanism of actions of these compounds.

Dynamic residue interaction network analysis of the oseltamivir binding site of N1 neuraminidase and its H274Y mutation site conferring drug resistance in influenza A virus

Background

Oseltamivir (OTV)-resistant influenza virus exhibits His-to-Tyr mutation at residue 274 (H274Y) in N1 neuraminidase (NA). However, the molecular mechanisms by which the H274Y mutation in NA reduces its binding affinity to OTV have not been fully elucidated.

Methods

In this study, we used dynamic residue interaction network (dRIN) analysis based on molecular dynamics simulation to investigate the correlation between the OTV binding site of NA and its H274Y mutation site.

Results

dRIN analysis revealed that the OTV binding site and H274Y mutation site of NA interact via the three interface residues connecting them. H274Y mutation significantly enhanced the interaction between residue 274 and the three interface residues in NA, thereby significantly decreasing the interaction between OTV and its surrounding loop 150 residues. Thus, we concluded that such changes in residue interactions could reduce the binding affinity of OTV to NA, resulting in drug resistant influenza viruses. Using dRIN analysis, we succeeded in understanding the characteristic changes in residue interactions due to H274Y mutation, which can elucidate the molecular mechanism of reduction in OTV binding affinity to influenza NA. Finally, the dRIN analysis used in this study can be widely applied to various systems such as individual proteins, protein-ligand complexes, and protein-protein complexes, to characterize the dynamic aspects of the interactions.

Interplay among Structural Stability, Plasticity, and Energetics Determined by Conformational Attuning of Flexible Loops in PD-1

The dynamics and plasticity of the PD-1/PD-L1 axis are the bottlenecks for the discovery of small-molecule antagonists to perturb this interaction interface significantly. Understanding the process of this protein-protein interaction (PPI) is of fundamental biological interest in structure-based drug designing. Food and Drug Administration (FDA)-approved anti-PD-1 monoclonal antibodies (mAbs) are the first-in-class with distinct binding modes to access this axis clinically; however, their mechanistic aspects remain elusive. Here, we have unveiled the interactive interfaces with PD-L1 and mAbs to investigate the native plasticity of PD-1 at global (structural and dynamical) and local (residue side-chain orientations) levels. We found that the structural stability and coordinated C_α movements are increased in the presence of PD-1's binding partners. The rigorous analysis of these PPIs using computational biophysical approaches revealed PD-1's intrinsic plasticity, its concerted loops' movement (BC, FG, and CC'), distal side-chain motions, and the thermodynamic landscape, which are perturbed remarkably from its unbound to bound states. Based on intra-/inter-residues' contact networks and energetics, the hot-spots have been identified that were found to be essential to arrest the dynamical motions of PD-1 significantly for the rational design of therapeutic agents by mimicking the mAbs mechanism.

Discovery of human coronaviruses pan-papain-like protease inhibitors using computational approaches

The papain-like protease (PLpro) is vital for the replication of coronaviruses (CoVs), as well as for escaping innate-immune responses of the host. Hence, it has emerged as an attractive antiviral drug-target. In this study, computational approaches were employed, mainly the structure-based virtual screening coupled with all-atom molecular dynamics (MD) simulations to computationally identify specific inhibitors of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) PLpro, which can be further developed as potential pan-PLpro based broad-spectrum antiviral drugs. The sequence, structure, and functional conserveness of most deadly human CoVs PLpro were explored, and it was revealed that functionally important catalytic triad residues are well conserved among SARS-CoV, SARS-CoV-2, and middle east respiratory syndrome coronavirus (MERS-CoV). The subsequent screening of a focused protease inhibitors database composed of ∼7,000 compounds resulted in the identification of three candidate compounds, ADM_13083841, LMG_15521745, and SYN_15517940. These three compounds established conserved interactions which were further explored through MD simulations, free energy calculations, and residual energy contribution estimated by MM-PB(GB)SA method. All these compounds showed stable conformation and interacted well with the active residues of SARS-CoV-2 PLpro, and showed consistent interaction profile with SARS-CoV PLpro and MERS-CoV PLpro as well. Conclusively, the reported SARS-CoV-2 PLpro specific compounds could serve as seeds for developing potent pan-PLpro based broad-spectrum antiviral drugs against deadly human coronaviruses. Moreover, the presented information related to binding site residual energy contribution could lead to further optimization of these compounds.

Identification of novel human USP2 inhibitor and its putative role in treatment of COVID-19 by inhibiting SARS-CoV-2 papain-like (PLpro) protease

Human ubiquitin carboxyl-terminal hydrolase-2 (USP2) inhibitors, such as thiopurine analogs, have been reported to inhibit SARS-CoV papain-like proteases (PLpro). The PLpro have significant functional implications in the innate immune response during SARS-CoV-2 infection and considered an important antiviral target. Both proteases share strikingly similar USP fold with right-handed thumb-palm-fingers structural scaffold and conserved catalytic triad Cys-His-Asp/Asn. In this urgency situation of COVID-19 outbreak, there is a lack of in-vitro facilities readily available to test SARS-CoV-2 inhibitors in whole-cell assays. Therefore, we adopted an alternate route to identify potential USP2 inhibitor through integrated in-silico efforts. After an extensive virtual screening protocol, the best compounds were selected and tested. The compound Z93 showed significant IC50 value against Jurkat (9.67 μM) and MOTL-4 cells (11.8 μM). The binding mode of Z93 was extensively analyzed through molecular docking, followed by MD simulations, and molecular interactions were compared with SARS-CoV-2. The relative binding poses of Z93 fitted well in the binding site of both proteases and showed consensus π-π stacking and H-bond interactions with histidine and aspartate/asparagine residues of the catalytic triad. These results led us to speculate that compound Z93 might be the first potential chemical lead against SARS-CoV-2 PLpro, which warrants in-vitro evaluations.

Sequence Characterization and Molecular Modeling of Clinically Relevant Variants of the SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (Mpro) is essential to viral replication and cleaves highly specific substrate sequences, making it an obvious target for inhibitor design. However, as for any virus, SARS-CoV-2 is subject to constant neutral drift and selection pressure, with new Mpro mutations arising over time. Identification and structural characterization of Mpro variants is thus critical for robust inhibitor design. Here we report sequence analysis, structure predictions, and molecular modeling for seventy-nine Mpro variants, constituting all clinically observed mutations in this protein as of April 29, 2020. Residue substitution is widely distributed, with some tendency toward larger and more hydrophobic residues. Modeling and protein structure network analysis suggest differences in cohesion and active site flexibility, revealing patterns in viral evolution that have relevance for drug discovery.

Discovery of novel Hepatitis C virus inhibitor targeting multiple allosteric sites of NS5B polymerase

HCV is a viral infection posing a severe global threat when left untreated progress to end-stage liver disease, including cirrhosis and HCC. The NS5B polymerase of HCV is the most potent target that harbors four allosteric binding sites that could interfere with the HCV infection. We present the discovery of a novel synthetic compound that harbors the potential of NS5B polymerase inhibition. All eight compounds belonging to the benzothiazine family of heterocycles displayed no cellular cytotoxicity in HepG2 cells at nontoxic dose concentration (200 μM). Subsequently, among eight compounds of the series, merely compound 5b exhibited significant inhibition of the expression of the HCV NS5B gene as compared to DMSO control in semi-quantitative PCR. Based on our western blot result, 5b at the range of 50, 100 and 200 μM induced 20, 40, and 70% inhibition of NS5B protein respectively. To estimate the binding potential, 5b was docked at respective allosteric sites followed by molecular dynamics (MD) simulations for a period of 20 ns. In addition, binding free energy calculation by MM-GB/PBSA method revealed a conserved interaction profile of residues lining the allosteric sites in agreement with the reported NS5B co-crystallized inhibitors. The presented results provide important information about a novel compound 5b which may facilitate the the discovery of novel inhibitors that tends to target multiple sites on NS5B polymerase.

Understanding the Hsp90 N-terminal Dynamics: Structural and Molecular Insights into the Therapeutic Activities of Anticancer Inhibitors Radicicol (RD) and Radicicol Derivative (NVP-YUA922)

Heat shock protein 90 (Hsp90) is a crucial component in carcinogenesis and serves as a molecular chaperone that facilitates protein maturation whilst protecting cells against temperature-induced stress. The function of Hsp90 is highly dependent on adenosine triphosphate (ATP) binding to the N-terminal domain of the protein. Thus, inhibition through displacement of ATP by means of competitive binding with a suitable organic molecule is considered an attractive topic in cancer research. Radicicol (RD) and its derivative, resorcinylic isoxazole amine NVP-AUY922 (NVP), have shown promising pharmacodynamics against Hsp90 activity. To date, the underlying binding mechanism of RD and NVP has not yet been investigated. In this study, we provide a comprehensive understanding of the binding mechanism of RD and NVP, from an atomistic perspective. Density functional theory (DFT) calculations enabled the analyses of the compounds' electronic properties and results obtained proved to be significant in which NVP was predicted to be more favorable with solvation free energy value of -23.3 kcal/mol and highest stability energy of 75.5 kcal/mol for a major atomic delocalization. Molecular dynamic (MD) analysis revealed NVP bound to Hsp90 (NT-NVP) is more stable in comparison to RD (NT-RD). The Hsp90 protein exhibited a greater binding affinity for NT-NVP (-49.4 ± 3.9 kcal/mol) relative to NT-RD (-28.9 ± 4.5 kcal/mol). The key residues influential in this interaction are Gly 97, Asp 93 and Thr 184. These findings provide valuable insights into the Hsp90 dynamics and will serve as a guide for the design of potent novel inhibitors for cancer treatment.

Gibbs Free Energy Calculation of Mutation in PncA and RpsA Associated With Pyrazinamide Resistance

A central approach for better understanding the forces involved in maintaining protein structures is to investigate the protein folding and thermodynamic properties. The effect of the folding process is often disturbed in mutated states. To explore the dynamic properties behind mutations, molecular dynamic (MD) simulations have been widely performed, especially in unveiling the mechanism of drug failure behind mutation. When comparing wild type (WT) and mutants (MTs), the structural changes along with solvation free energy (SFE), and Gibbs free energy (GFE) are calculated after the MD simulation, to measure the effect of mutations on protein structure. Pyrazinamide (PZA) is one of the first-line drugs, effective against latent Mycobacterium tuberculosis isolates, affecting the global TB control program 2030. Resistance to this drug emerges due to mutations in pncA and rpsA genes, encoding pyrazinamidase (PZase) and ribosomal protein S1 (RpsA) respectively. The question of how the GFE may be a measure of PZase and RpsA stabilities, has been addressed in the current review. The GFE and SFE of MTs have been compared with WT, which were already found to be PZA-resistant. WT structures attained a more stable state in comparison with MTs. The physiological effect of a mutation in PZase and RpsA may be due to the difference in energies. This difference between WT and MTs, depicted through GFE plots, might be useful in predicting the stability and PZA-resistance behind mutation. This study provides useful information for better management of drug resistance, to control the global TB problem.

Structural features in common of HBV and HIV-1 resistance against chirally-distinct nucleoside analogues entecavir and lamivudine

Chronic hepatitis B virus (HBV) infection is a major public health problem that affects millions of people worldwide. Nucleoside analogue reverse transcriptase (RT) inhibitors, such as entecavir (ETV) and lamivudine (3TC), serve as crucial anti-HBV drugs. However, structural studies of HBV RT have been hampered due to its unexpectedly poor solubility. Here, we show that human immunodeficiency virus type-1 (HIV-1) with HBV-associated amino acid substitutions Y115F/F116Y/Q151M in its RT (HIV^{Y115F/F116Y/Q151M}) is highly susceptible to ETV and 3TC. Additionally, we experimentally simulated previously reported ETV/3TC resistance for HBV using HIV^{Y115F/F116Y/Q151M} with F160M/M184V (L180M/M204V in HBV RT) substituted. We determined crystal structures for HIV-1 RT^{Y115F/F116Y/Q151M}:DNA complexed with 3TC-triphosphate (3TC-TP)/ETV-triphosphate (ETV-TP)/dCTP/dGTP. These structures revealed an atypically tight binding conformation of 3TC-TP, where the Met184 side-chain is pushed away by the oxathiolane of 3TC-TP and exocyclic methylene of ETV-TP. Structural analysis of RT^{Y115F/F116Y/Q151M/F160M/M184V}:DNA:3TC-TP also demonstrated that the loosely bound 3TC-TP is misaligned at the active site to prevent a steric clash with the side chain γ-methyl of Val184. These findings shed light on the common structural mechanism of HBV and HIV-1 resistance to 3TC and ETV and should aid in the design of new agents to overcome drug resistance to 3TC and ETV.

Understanding the Pyrimethamine Drug Resistance Mechanism via Combined Molecular Dynamics and Dynamic Residue Network Analysis

In this era of precision medicine, insights into the resistance mechanism of drugs are integral for the development of potent therapeutics. Here, we sought to understand the contribution of four point mutations (N51I, C59R, S108N, and I164L) within the active site of the malaria parasite enzyme dihydrofolate reductase (DHFR) towards the resistance of the antimalarial drug pyrimethamine. Homology modeling was used to obtain full-length models of wild type (WT) and mutant DHFR. Molecular docking was employed to dock pyrimethamine onto the generated structures. Subsequent all-atom molecular dynamics (MD) simulations and binding free-energy computations highlighted that pyrimethamine's stability and affinity inversely relates to the number of mutations within its binding site and, hence, resistance severity. Generally, mutations led to reduced binding affinity to pyrimethamine and increased conformational plasticity of DHFR. Next, dynamic residue network analysis (DRN) was applied to determine the impact of mutations and pyrimethamine binding on communication dispositions of DHFR residues. DRN revealed residues with distinctive communication profiles, distinguishing WT from drug-resistant mutants as well as pyrimethamine-bound from pyrimethamine-free models. Our results provide a new perspective on the understanding of mutation-induced drug resistance.

Structure and dynamics of ionic liquid tolerant hyperthermophilic endoglucanase Cel12A from Rhodothermus marinus

Economic deconstruction of lignocellulose remains a challenge due to the complex architecture of cellulose, hemicellulose, and lignin. Advancements in pretreatment processes have introduced ionic liquids (ILs) as promising non-derivatizing solvents for reducing biomass recalcitrance and for promoting enzymatic hydrolysis. However, available commercial cellulases are destabilized or inactivated even in low concentration of residual ILs. Thus, a molecular understanding of IL-enzyme interactions is crucial for developing IL-tolerant enzymes with high catalytic activity. In this study, molecular insight behind the IL tolerance of hyperthermophilic endoglucanase Cel12A from Rhodothermus marinus (RmCel12A) has been investigated in 20%, 40%, and 60% 1-ethyl-3-methylimidazolium acetate (EmimAc) through molecular dynamic simulations at 368 K. Though the enzyme retained its stability in all EmimAc concentrations, the activity was affected due to the loss of essential dynamic motions. A protein structure network was constructed using the snapshots of protein structures from the simulation trajectories and the hub properties of residues R20, Y59, W68, W197, E203, and F220 were found to be lost in 60% EmimAc. Emim cations were observed to intrude the active site tunnel and interact with more number of catalytic residues with higher cumulative fractional occupancy in 60% EmimAc than in 20% or 40% EmimAc. Some non-catalytic residues have also been identified at the active site, which can be probable mutation targets for improving the IL tolerance. Our findings reveal the molecular understanding behind the origin of activity loss of RmCel12A and proposed insights for the further improvement of IL sensitivity.

Understanding the behavior of ionic liquid tolerant hyperthermophilic endoglucanase Cel12A from Rhodothermus marinus in different concentrations of EmimAc.

Perspectives towards antiviral drug discovery against Ebola virus

Ebola virus disease (EVD), caused by Ebola viruses, resulted in more than 11 500 deaths according to a recent 2018 WHO report. With mortality rates up to 90%, it is nowadays one of the most deadly infectious diseases. However, no Food and Drug Administration‐approved Ebola drugs or vaccines are available yet with the mainstay of therapy being supportive care. The high fatality rate and absence of effective treatment or vaccination make Ebola virus a category‐A biothreat pathogen. Fortunately, a series of investigational countermeasures have been developed to control and prevent this global threat. This review summarizes the recent therapeutic advances and ongoing research progress from research and development to clinical trials in the development of small‐molecule antiviral drugs, small‐interference RNA molecules, phosphorodiamidate morpholino oligomers, full‐length monoclonal antibodies, and vaccines. Moreover, difficulties are highlighted in the search for effective countermeasures against EVD with additional focus on the interplay between available in silico prediction methods and their evidenced potential in antiviral drug discovery.

Exploring the Lapse in Druggability: Sequence Analysis, Structural Dynamics and Binding Site Characterization of K-RasG12C Variant, a Feasible Oncotherapeutics Target

BACKGROUND

The difficulty in druggability of K-Ras variant has presented a challenge in the treatment of cancer diseases associated with its dysfunctionality. Despite the identification of different binding sites, limited information exists in the literature about their characteristics. Therefore, identification, crossvalidation and characterization of its druggable sites would aid the design of chemical compounds that will arrest its dysfunctionality related oncogenesis.

OBJECTIVE

This study entails the identification, cross-validation and characterization of K-Ras G12C variant's binding sites for potential druggability, coupled with the elucidation of alterations in 3D conformations and dynamics.

METHOD

Molecular dynamics simulation was carried out on the inactive, the active and the hyperactive K-RasG12Cvariant using the amber software package. The SiteMap software was employed in identifying and characterizing the druggable binding sites while the validation of the binding sites was carried out with the SiteHound and MetaPocket servers.

RESULTS

Four druggable binding sites were identified, validated and characterized based on physicochemical attributes such as size, volume, degree of enclosure or exposure, degree of contact, hydrophobic/hydrophilic character, hydrophobic/hydrophilic balance and hydrogen-bonding features. Conformational studies also revealed that the K-Ras variant exhibited notable structural instability, increased flexibility and a strongly anticorrelated movement compared to the inactive and active wildtype forms.

CONCLUSION

The attributes of the characterized druggable sites will be useful in designing site-specific K-Ras inhibitors for the treatment of K-Ras variant associated cancer diseases.

Broadening the horizon: Integrative pharmacophore-based and cheminformatics screening of novel chemical modulators of mitochondria ATP synthase towards interventive Alzheimer's disease therapy

The proven efficacy of J147 in the treatment of Alzheimer's disease (AD) has been emphatic, particularly since its selective modulatory roles towards mitochondrial ATP synthase (mATPase) were defined. This prospect, if methodically probed, could further pave way for the discovery of novel anti-AD drugs with improved pharmacokinetics and therapeutic potential. To this effect, for the first time, we employed a four-step paradigm that integrated our in-house per-residue energy decomposition (PRED) protocol coupled with molecular dynamics, cheminformatics and analytical binding free energy methods. This was geared towards the screening and identification of new leads that exhibit modulatory potentials towards mATPase in a J147-similar pattern. Interestingly, from a large-scale library of compounds, we funnelled down on three potential hits that demonstrated selective and high-affinity binding activities towards α-F1-ATP synthase (ATP5A) relative to J147. Moreover, these compounds exhibited higher binding propensity with a differential ΔGs greater than -1 kcal/mol comparative to J147, and also elicited distinct modulatory effects on ATP5A domain structures. More interestingly, per-residue pharmacophore modeling of these lead compounds revealed similar interactive patterns with crucial residues at the α-site region of ATP5A characterized by high energy contributions based on binding complementarity. Recurrent target residues involved in high-affinity interactions with the hit molecules relative to J147 include Arg1112 and Gln426. Furthermore, assessments of pharmacokinetics revealed that the lead compounds were highly drug-like with minimal violations of the Lipinski's rule of five. As developed in this study, the most extrapolative pharmacophore model of the selected hits encompassed three electron donors and one electron acceptor. We speculate that these findings will be fundamental to the reformation of anti-AD drug discovery procedures.

In silico structural elucidation of RNA-dependent RNA polymerase towards the identification of potential Crimean-Congo Hemorrhagic Fever Virus inhibitors

The Crimean-Congo Hemorrhagic Fever virus (CCHFV) is a segmented negative single-stranded RNA virus (-ssRNA) which causes severe hemorrhagic fever in humans with a mortality rate of ~50%. To date, no vaccine has been approved. Treatment is limited to supportive care with few investigational drugs in practice. Previous studies have identified viral RNA dependent RNA Polymerase (RdRp) as a potential drug target due to its significant role in viral replication and transcription. Since no crystal structure is available yet, we report the structural elucidation of CCHFV-RdRp by in-depth homology modeling. Even with low sequence identity, the generated model suggests a similar overall structure as previously reported RdRps. More specifically, the model suggests the presence of structural/functional conserved RdRp motifs for polymerase function, the configuration of uniform spatial arrangement of core RdRp sub-domains, and predicted positively charged entry/exit tunnels, as seen in sNSV polymerases. Extensive pharmacophore modeling based on per-residue energy contribution with investigational drugs allowed the concise mapping of pharmacophoric features and identified potential hits. The combination of pharmacophoric features with interaction energy analysis revealed functionally important residues in the conserved motifs together with in silico predicted common inhibitory binding modes with highly potent reference compounds.

Dynamic perspectives into the mechanisms of mutation-induced p53-DNA binding loss and inactivation using active perturbation theory: Structural and molecular insights toward the design of potent reactivators in cancer therapy

The DNA-binding ability of p53 represents the crux of its tumor suppressive activities, which involves transcriptional activation of target genes responsible for apoptosis and cell-cycle arrest. Mutational occurrences within or in close proximity to the DNA-binding surface of p53 have accounted for the loss of direct DNA-binding ability and inactivation implicated in many cases of cancer. Moreover, the design of therapeutic compounds that can restore DNA-binding ability in p53 mutants has been identified as a way forward in curtailing their oncogenic activities. However, there is still the need for more insights into evaluate the perturbations that occur at the DNA-binding interface of mp53 relative to DNA-binding loss, inactivation, and design of potent reactivators, hence the purpose of this study. Therefore, we evaluated p53-structural (R175H) and contact (R273C) mutational effects using tunnel perturbation analysis and other computational tools. We identified significant perturbations in the active tunnels of p53, which resulted in altered geometry and loss, unlike in the wild-type p53. This corroborated with structural, DNA-binding, and interaction network analysis, which showed that loss of flexibility, repulsion of DNA-interactive residues, and instability occurred at the binding interface of both mutants. Also, these mutations altered bonding interactions and network topology at the DNA-binding interface, resulting in the reduction of p53-DNA binding proximity and affinity. Therefore, these findings would aid the structure-based design of novel chemical entities capable of restoring p53-DNA binding and activation.

Induced Mutation Proves a Potential Target for TB Therapy: A Molecular Dynamics Study on LprG

Molecular dynamics (MD) simulations of wild-type and V91W mutant Mycobacterium tuberculosis-LprG (Mtb-LprG) were performed with the goal to provide a comprehensive understanding of the Mtb-LprG as a potential antimycobacterial target. A long-range MD simulations and post-MD analyzes led us to various results that plainly explained the impact of V91W mutation on Mtb-LprG. Herein, the results revealed that the wild-type is less stable compared to V91W mutant. This was further supported by root mean square fluctuation, where the V91W mutant showed a higher degree of flexibility compared to the wild-type. Dynamic cross-correlation analysis revealed that induced mutation leads to higher residual flexibility in the mutant structure as compared to the wild-type structure thus resulting in the existence of negatively correlated motions. The difference in principal component analysis scatter plot across the first two normal modes suggests a greater mobility of the V91W mutant conformation compared to the wild-type. Thermodynamic calculations revealed that the van der Waals (E_vdw) forces contribute the most towards binding free energy in a case of the V91W mutant as compared to the wild-type LprG complex. In addition, the residue interaction networks revealed more of E_vdw interaction existence among residues in case of the V91W mutant. This study supports the Mtb-LprG as a potential antimycobacterial target and also serves as a cornerstone to identifying new potential targets that have no inhibitors.

Across the blood-brain barrier: Neurotherapeutic screening and characterization of naringenin as a novel CRMP-2 inhibitor in the treatment of Alzheimer's disease using bioinformatics and computational tools

The discovery and developmental processes of CNS drugs have been limited by the inability of potential drug molecules to pass through the blood-brain barrier (BBB). This presents a significant setback in the treatment of neurodegenerative disorders such as Alzheimer's disease (AD), hence the need for compounds that can adhere strictly to the selective criteria of suitable CNS drugs. Collapsin response mediator protein-2 (CRMP-2) has been recently identified as a viable target in neurotherapeutics due to its involvement in the etiology of AD. As shown in previous studies, Naringenin (NAR), a small molecule derivative of Drynaria rhizome (DR) extract, specifically binds CRMP-2 and reduces its phosphorylation. This was shown to facilitate axonal regrowth, with improvement in cognition and learning. Herein, we report the first account of the use of cheminformatics techniques to define the CNS drug-suitability of NAR using selective criteria, coupled with the prediction of possible biological activities and toxicities. Also, we evaluated the mechanistic activity of NAR by modeling its molecular interaction with human CRMP-2 (hCRMP-2). Physicochemical analyses revealed the suitability of NAR as a CNS drug and its ability to transverse the BBB. Possible neurogenic, anti-carcinogenic and cardioprotective activities were also predicted. NAR exhibited favorable binding to CRMP-2 and formed strong bonds with active site residues, which accounts for its stabilization and affinity. Moreover, NAR induced notable conformational changes in CRMP-2, an occurrence that could possibly disrupt kinase-mediated phosphorylation. These findings will aid in the optimization of NAR and improve its neurotherapeutic activities in the treatment of AD.

The impact of Thr91 mutation on c-Src resistance to UM-164: molecular dynamics study revealed a new opportunity for drug design

The emergence of a drug resistant non-receptor tyrosine kinase (c-Src) in triple-negative breast cancer (TNBC) remains a prime concern in relation to the burden of TNBC among people living with breast cancer and drug development. Thr91 mutation was found to induce a complete loss of protein conformation required for drug fitness. Herein, we provide the first account of the molecular impact of the Thr91 mutation on c-Src resistance to experimental drug UM-164 using various computational approaches, namely molecular dynamics simulation, principal component analysis (PCA), dynamic cross-correlation matrices (DCCM) analysis, hydrogen bond occupancy, thermodynamics calculation, ligand-residue interaction and residue interaction networks (RINs). Findings from this study revealed that Thr91 mutation leads to a steric conflict between UM-164 and the side chain of methionine (Met91); this mutation distorts the UM-164 optimum orientation on the conformational space of mutant c-Src compared to the wild-type; decreases hydrogen bond formation between the residues in the mutant protein structure; decreases the UM-164 binding energy in the mutant by -13.416 kcal mol^-1; reduces the residue correlation in the mutant protein structure; induces a change in the overall protein structure conformation from an inactive to active conformation; and distorts the ligand atomic interaction network and the residue interaction network. This report provides important insights that will assist in the further design of novel dual kinase inhibitors to minimise the chances of drug resistance in triple negative breast cancer.

MD-TASK: a software suite for analyzing molecular dynamics trajectories

SUMMARY

Molecular dynamics (MD) determines the physical motions of atoms of a biological macromolecule in a cell-like environment and is an important method in structural bioinformatics. Traditionally, measurements such as root mean square deviation, root mean square fluctuation, radius of gyration, and various energy measures have been used to analyze MD simulations. Here, we present MD-TASK, a novel software suite that employs graph theory techniques, perturbation response scanning, and dynamic cross-correlation to provide unique ways for analyzing MD trajectories.

AVAILABILITY AND IMPLEMENTATION

MD-TASK has been open-sourced and is available for download from https://github.com/RUBi-ZA/MD-TASK , implemented in Python and supported on Linux/Unix.

CONTACT

o.tastanbishop@ru.ac.za.

Tailored-pharmacophore model to enhance virtual screening and drug discovery: a case study on the identification of potential inhibitors against drug-resistant Mycobacterium tuberculosis (3R)-hydroxyacyl-ACP dehydratases

Aim: Virtual screening (VS) is powerful tool in discovering molecular inhibitors that are most likely to bind to drug targets of interest. Herein, we introduce a novel VS approach, so-called ‘tailored-pharmacophore’, in order to explore inhibitors that overcome drug resistance. Methodology & results: The emergence and spread of drug resistance strains of tuberculosis is one of the most critical issues in healthcare. A tailored-pharmacophore approach was found promising to identify in silico predicted hit with better binding affinities in case of the resistance mutations in MtbHadAB as compared with thiacetazone, a prodrug used in the clinical treatment of tuberculosis. Conclusion: This approach can potentially be enforced for the discovery and design of drugs against a wide range of resistance targets.

VP40 of the Ebola Virus as a Target for EboV Therapy: Comprehensive Conformational and Inhibitor Binding Landscape from Accelerated Molecular Dynamics

The first account of the dynamic features of the loop region of VP40 of the Ebola virus was studied using accelerated molecular dynamics simulations and reported herein. Among the proteins of the Ebola virus, the matrix protein (VP40) plays a significant role in the virus lifecycle thereby making it a promising therapeutic target. Of interest is the newly elucidated N-terminal domain loop region of VP40 comprising residues K127, T129, and N130 which when mutated to alanine have demonstrated an unrecognized role for N-terminal domain-plasma membrane interaction for efficient VP40-plasma membrane localization, oligomerization, matrix assembly, and egress. The molecular understanding of the conformational features of VP40 in complex with a known inhibitor still remains elusive. Using accelerated molecular dynamics approaches, we conducted a comparative study on VP40 apo and bound systems to understand the conformational features of VP40 at the molecular level and to determine the effect of inhibitor binding with the aid of a number of post-dynamic analytical tools. Significant features were seen in the presence of an inhibitor as per molecular mechanics/generalized born surface area binding free energy calculations. Results revealed that inhibitor binding to VP40 reduces the flexibility and mobility of the protein as supported by root mean square fluctuation and root mean square deviation calculations. The study revealed a characteristic "twisting" motion and coiling of the loop region of VP40 accompanied by conformational changes in the dimer interface upon inhibitor binding. We believe that results presented in this study will ultimately provide useful insight into the binding landscape of VP40 which could assist researchers in the discovery of potent Ebola virus inhibitors for anti-Ebola therapies.

Sliding Clamp of DNA Polymerase III as a Drug Target for TB Therapy: Comprehensive Conformational and Binding Analysis from Molecular Dynamic Simulations

Tuberculosis, caused by Mycobacterium tuberculosis, is one of the most common causes of death in the world. Mycobacterium tuberculosis -sliding clamp is a protein essential for many important DNA transactions including replication and DNA repair proteins, thus, a potential drug target for tuberculosis. Further investigation is needed in understanding DNA polymerase sliding clamp structure, especially from a computational perspective. In this study, we employ a wide-range of comparative molecular dynamic analyses on two systems: Mycobacterium tuberculosis - sliding clamp enzyme in its apo and bound form. The results reported in this study shows apo conformation to be less stable, as compared to bound conformation with an average radius of gyration of 25.812 and 25.459 Å, respectively. This was further supported by root mean square fluctuation, where an apo enzyme showed a higher degree of flexibility. However, the presence of the ligand lowers radius of gyration and root mean square fluctuation and also leads to an existence of negative correlated motions. Principal component analysis further justifies the same findings, whereby the apo enzyme exhibits a higher fluctuation compared to the bound complex. In addition, a stable 3₁₀ helix located at the binding site appears to be unstable in the presence of the ligand. Hence, it is possible that the binding of the ligand may have caused a rearrangement of the structure, leading to a change in the unwinding of 3₁₀ helix. Findings reported in this study further enhance the understanding of Mycobacterium tuberculosis -DnaN and also give a lead to the development of potent tuberculosis drugs.

CD4-gp120 interaction interface - a gateway for HIV-1 infection in human: molecular network, modeling and docking studies

The major causative agent for Acquired Immune Deficiency Syndrome (AIDS) is Human Immunodeficiency Virus-1 (HIV-1). HIV-1 is a predominant subtype of HIV which counts on human cellular mechanism virtually in every aspect of its life cycle. Binding of viral envelope glycoprotein-gp120 with human cell surface CD4 receptor triggers the early infection stage of HIV-1. This study focuses on the interaction interface between these two proteins that play a crucial role for viral infectivity. The CD4-gp120 interaction interface has been studied through a comprehensive protein-protein interaction network (PPIN) analysis and highlighted as a useful step towards identifying potential therapeutic drug targets against HIV-1 infection. We prioritized gp41, Nef and Tat proteins of HIV-1 as valuable drug targets at early stage of viral infection. Lack of crystal structure has made it difficult to understand the biological implication of these proteins during disease progression. Here, computational protein modeling techniques and molecular dynamics simulations were performed to generate three-dimensional models of these targets. Besides, molecular docking was initiated to determine the desirability of these target proteins for already available HIV-1 specific drugs which indicates the usefulness of these protein structures to identify an effective drug combination therapy against AIDS.

Mechanism of Inhibition of Hsp90 Dimerization by Gyrase B Inhibitor Coumermycin A1 (C-A1) Revealed by Molecular Dynamics Simulations and Thermodynamic Calculations

Heat shock protein (Hsp) 90 an emerging and attracting target in the anti-HIV drug discovery process due to the key role it plays in the pathogenicity of HIV-1 virus. In this research study, long-range all-atom molecular dynamics simulations were engaged for the bound and the unbound proteins to enhance the understanding of the molecular mechanisms of the Hsp90 dimerization and inhibition. Results evidently showed that coumermycin A1 (C-A1), a recently discovered Hsp90 inhibitor, binds at the dimer's active site of the Hsp90 protein and leads to a substantial parting between dimeric opposed residues, which include Arg591.B, Lys594.A, Ser663.A, Thr653.B, Ala665.A, Thr649.B, Leu646.B and Asn669.A. Significant differences in magnitudes were observed in radius of gyration, root-mean-square deviation and root-mean-square fluctuation, which confirms a reasonably more flexible state in the apo conformation associated with it dimerization. In contrast, the bound conformer of Hsp90 showed less flexibility. This visibly highpoints the inhibition process resulting from the binding of the ligand. These findings were further validated by principal component analysis. We believe that the detailed dynamic analyses of Hsp90 presented in this study, would give an imperative insight and better understanding to the function and mechanisms of inhibition. Furthermore, information obtained from the binding mode of the inhibitor would be of great assistance in the design of more potent inhibitors against the HIV target Hsp90.

Perturbations in inter-domain associations may trigger the onset of pathogenic transformations in PrP(C): insights from atomistic simulations

Conversion of the predominantly α-helical cellular prion protein (PrP(C)) to the misfolded β-sheet enriched Scrapie form (PrP(Sc)) is a critical event in prion pathogenesis. However, the conformational triggers that lead to the isoform conversion (PrP(C) to PrP(Sc)) remain obscure, and conjectures about the role of unusually hydrophilic, short helix H1 of the C-terminal globular domain in the transition are varied. Helix H1 is anchored to helix H3 via a few stabilizing polar interactions. We have employed fully atomistic molecular dynamics simulations to study the effects triggered by a minor perturbation in the network of these non-bonded interactions in PrP(C). The elimination of just one of the key H1-H3 hydrogen bonds led to a cascade of conformational changes that are consistent with those observed in partially unfolded intermediates of PrP(C), with pathogenic mutations and in low pH environments. Our analyses reveal that the perturbation results in the enhanced conformational flexibility of the protein. The resultant enhancement in the dynamics leads to overall increased solvent exposure of the hydrophobic core residues and concomitant disruption of the H1-H3 inter-domain salt bridge network. This study lends credence to the hypothesis that perturbing the cooperativity of the stabilizing interactions in the PrP(C) globular domain can critically affect its dynamics and may lead to structural transitions of pathological relevance.

E2 superfamily of ubiquitin-conjugating enzymes: constitutively active or activated through phosphorylation in the catalytic cleft

Protein phosphorylation is a modification that offers a dynamic and reversible mechanism to regulate the majority of cellular processes. Numerous diseases are associated with aberrant regulation of phosphorylation-induced switches. Phosphorylation is emerging as a mechanism to modulate ubiquitination by regulating key enzymes in this pathway. The molecular mechanisms underpinning how phosphorylation regulates ubiquitinating enzymes, however, are elusive. Here, we show the high conservation of a functional site in E2 ubiquitin-conjugating enzymes. In catalytically active E2s, this site contains aspartate or a phosphorylatable serine and we refer to it as the conserved E2 serine/aspartate (CES/D) site. Molecular simulations of substrate-bound and -unbound forms of wild type, mutant and phosphorylated E2s, provide atomistic insight into the role of the CES/D residue for optimal E2 activity. Both the size and charge of the side group at the site play a central role in aligning the substrate lysine toward E2 catalytic cysteine to control ubiquitination efficiency. The CES/D site contributes to the fingerprint of the E2 superfamily. We propose that E2 enzymes can be divided into constitutively active or regulated families. E2s characterized by an aspartate at the CES/D site signify constitutively active E2s, whereas those containing a serine can be regulated by phosphorylation.

Insights on the structural perturbations in human MTHFR Ala222Val mutant by protein modeling and molecular dynamics

Methylenetetrahydrofolate reductase (MTHFR) protein catalyzes the only biochemical reaction which produces methyltetrahydrofolate, the active form of folic acid essential for several molecular functions. The Ala222Val polymorphism of human MTHFR encodes a thermolabile protein associated with increased risk of neural tube defects and cardiovascular disease. Experimental studies have shown that the mutation does not affect the kinetic properties of MTHFR, but inactivates the protein by increasing flavin adenine dinucleotide (FAD) loss. The lack of completely solved crystal structure of MTHFR is an impediment in understanding the structural perturbations caused by the Ala222Val mutation; computational modeling provides a suitable alternative. The three-dimensional structure of human MTHFR protein was obtained through homology modeling, by taking the MTHFR structures from Escherichia coli and Thermus thermophilus as templates. Subsequently, the modeled structure was docked with FAD using Glide, which revealed a very good binding affinity, authenticated by a Glide XP score of -10.3983 (kcal mol(-1)). The MTHFR was mutated by changing Alanine 222 to Valine. The wild-type MTHFR-FAD complex and the Ala222Val mutant MTHFR-FAD complex were subjected to molecular dynamics simulation over 50 ns period. The average difference in backbone root mean square deviation (RMSD) between wild and mutant variant was found to be ~.11 Å. The greater degree of fluctuations in the mutant protein translates to increased conformational stability as a result of mutation. The FAD-binding ability of the mutant MTHFR was also found to be significantly lowered as a result of decreased protein grip caused by increased conformational flexibility. The study provides insights into the Ala222Val mutation of human MTHFR that induces major conformational changes in the tertiary structure, causing a significant reduction in the FAD-binding affinity.

Flap flexibility amongst plasmepsins I, II, III, IV, and V: Sequence, structural, and molecular dynamics analyses

Herein, for the first time, we comparatively report the opening and closing of apo plasmepsin I - V. Plasmepsins belong the aspartic protease family of enzymes, and are expressed during the various stages of the P. falciparum lifecycle, the species responsible for the most lethal and virulent malaria to infect humans. Plasmepsin I, II, IV and HAP degrade hemoglobin from infected red blood cells, whereas plasmepsin V transport proteins crucial to the survival of the malaria parasite across the endoplasmic reticulum. Flap-structures covering the active site of aspartic proteases (such as HIV protease) are crucial to the conformational flexibility and dynamics of the protein, and ultimately control the binding landscape. The flap-structure in plasmepsins is made up of a flip tip in the N-terminal lying perpendicular to the active site, adjacent to the flexible loop region in the C-terminal. Using molecular dynamics, we propose three parameters to better describe the opening and closing of the flap-structure in apo plasmepsins. Namely, the distance, d1, between the flap tip and the flexible region; the dihedral angle, ϕ, to account for the twisting motion; and the TriCα angle, θ1. Simulations have shown that as the flap-structure twists, the flap and flexible region move apart opening the active site, or move toward each other closing the active site. The data from our study indicate that of all the plasmepsins investigated in the present study, Plm IV and V display the highest conformational flexibility and are more dynamic structures versus Plm I, II, and HAP.

Effect of T68A/N126Y mutations on the conformational and ligand binding landscape of Coxsackievirus B3 3C protease

3C protease of Coxsackievirus B3 (CVB3) plays an essential role in the viral replication cycle, and therefore, emerged as an attractive therapeutic target for the treatment of human diseases caused by CVB3 infection. In this study, we report the first account of the molecular impact of the T68A/N126Y double mutant (Mutant(Bound)) using an integrated computational approach. Molecular dynamics simulation and post-dynamics binding free energy, principal component analysis (PCA), hydrogen bond occupancy, SASA, R(g) and RMSF confirm that T68A/N126Y instigated an increased conformational flexibility due to the loss of intra- and inter-molecular hydrogen bond interactions and other prominent binding forces, which led to a decreased protease grip on the ligand (3CPI). The double mutations triggered a distortion orientation of 3CPI in the active site and decreases the binding energy, ΔG(bind) (∼3 kcal mol(-1)), compared to the wild type (Wild(Bound)). The van der Waals and electrostatic energy contributions coming from residues 68 and 126 are lower for Mutant(Bound) when compared with Wild(Bound). In addition, variation in the overall enzyme motion as evident from the PCA, distorted hydrogen bonding network and loss of protein-ligand interactions resulted in a loss of inhibitor efficiency. The comprehensive molecular insight gained from this study should be of great importance in understanding the drug resistance against CVB3 3C protease; also, it will assist in the designing of novel Coxsackievirus B3 inhibitors with high ligand efficacy on resistant strains.

Proline substitutions in a Mip-like peptidyl-prolyl cis-trans isomerase severely affect its structure, stability, shape and activity

FKBP22, an Escherichia coli-specific peptidyl-prolyl cis-trans isomerase, shows substantial homology with the Mip-like virulence factors. Mip-like proteins are homodimeric and possess a V-shaped conformation. Their N-terminal domains form dimers, whereas their C-terminal domains bind protein/peptide substrates and distinct inhibitors such as rapamycin and FK506. Interestingly, the two domains of the Mip-like proteins are separated by a lengthy, protease-susceptible α-helix. To delineate the structural requirement of this domain-connecting region in Mip-like proteins, we have investigated a recombinant FKBP22 (rFKBP22) and its three point mutants I65P, V72P and A82P using different probes. Each mutant harbors a Pro substitution mutation at a distinct location in the hinge region. We report that the three mutants are not only different from each other but also different from rFKBP22 in structure and activity. Unlike rFKBP22, the three mutants were unfolded by a non-two state mechanism in the presence of urea. In addition, the stabilities of the mutants, particularly I65P and V72P, differed considerably from that of rFKBP22. Conversely, the rapamycin binding affinity of no mutant was different from that of rFKBP22. Of the mutants, I65P showed the highest levels of structural/functional loss and dissociated partly in solution. Our computational study indicated a severe collapse of the V-shape in I65P due to the anomalous movement of its C-terminal domains. The α-helical nature of the domain-connecting region is, therefore, critical for the Mip-like proteins.

Heat-shock protein 90 (Hsp90) as anticancer target for drug discovery: an ample computational perspective

There are over 100 different types of cancer, and each is classified based on the type of cell that is initially affected. If left untreated, cancer can result in serious health problems and eventually death. Recently, the paradigm of cancer chemotherapy has evolved to use a combination approach, which involves the use of multiple drugs each of which targets an individual protein. Inhibition of heat-shock protein 90 (Hsp90) is one of the novel key cancer targets. Because of its ability to target several signaling pathways, Hsp90 inhibition emerged as a useful strategy to treat a wide variety of cancers. Molecular modeling approaches and methodologies have become 'close counterparts' to experiments in drug design and discovery workflows. A wide range of molecular modeling approaches have been developed, each of which has different objectives and outcomes. In this review, we provide an up-to-date systematic overview on the different computational models implemented toward the design of Hsp90 inhibitors as anticancer agents. Although this is the main emphasis of this review, different topics such as background and current statistics of cancer, different anticancer targets including Hsp90, and the structure and function of Hsp90 from an experimental perspective, for example, X-ray and NMR, are also addressed in this report. To the best of our knowledge, this review is the first account, which comprehensively outlines various molecular modeling efforts directed toward identification of anticancer drugs targeting Hsp90. We believe that the information, methods, and perspectives highlighted in this report would assist researchers in the discovery of potential anticancer agents.

Multi-drug resistance profile of PR20 HIV-1 protease is attributed to distorted conformational and drug binding landscape: molecular dynamics insights

The PR20 HIV-1 protease, a variant with 20 mutations, exhibits high levels of multi-drug resistance; however, to date, there has been no report detailing the impact of these 20 mutations on the conformational and drug binding landscape at a molecular level. In this report, we demonstrate the first account of a comprehensive study designed to elaborate on the impact of these mutations on the dynamic features as well as drug binding and resistance profile, using extensive molecular dynamics analyses. Comparative MD simulations for the wild-type and PR20 HIV proteases, starting from bound and unbound conformations in each case, were performed. Results showed that the apo conformation of the PR20 variant of the HIV protease displayed a tendency to remain in the open conformation for a longer period of time when compared to the wild type. This led to a phenomena in which the inhibitor seated at the active site of PR20 tends to diffuse away from the binding site leading to a significant change in inhibitor-protein association. Calculating the per-residue fluctuation (RMSF) and radius of gyration, further validated these findings. MM/GBSA showed that the occurrence of 20 mutations led to a drop in the calculated binding free energies (ΔGbind) by ~25.17 kcal/mol and ~5 kcal/mol for p2-NC, a natural peptide substrate, and darunavir, respectively, when compared to wild type. Furthermore, the residue interaction network showed a diminished inter-residue hydrogen bond network and changes in inter-residue connections as a result of these mutations. The increased conformational flexibility in PR20 as a result of loss of intra- and inter-molecular hydrogen bond interactions and other prominent binding forces led to a loss of protease grip on ligand. It is interesting to note that the difference in conformational flexibility between PR20 and WT conformations was much higher in the case of substrate-bound conformation as compared to DRV. Thus, developing analogues of DRV by retaining its key pharmacophore features will be the way forward in the search for novel protease inhibitors against multi-drug resistant strains.

Single H5N1 influenza A neuraminidase mutation develops resistance to oseltamivir due to distorted conformational and drug binding landscape: multiple molecular dynamics analyses

Clinical studies showed that a single mutation, I117V, develops severe resistance to oseltamivir, the first orally active influenza A neuraminidase inhibitor, in highly pathogenic H5N1 influenza A viruses.

Reaching beyond HIV/HCV: nelfinavir as a potential starting point for broad-spectrum protease inhibitors against dengue and chikungunya virus

Re-purposing HIV/HCV inhibitors against DENV and CHIKV using computer aided drug design.

Flap dynamics of plasmepsin proteases: insight into proposed parameters and molecular dynamics

Herein, for the first time, we report the flap opening and closing in Plasmepsin proteases – plasmepsin II (PlmII) was used as a prototype model.