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

Molecularly imprinted polymer-based sensors for identification volatile compounds in pharmaceutical products: in silico rational design

The present study aimed to strategically design a Molecularly Imprinted Polymer (MIP) with selective extraction capabilities for volatile compounds found in pork. These specific volatile compounds, such as 3-methyl-1-butanol, 1-nonanal, octanal, hexanal, 2-pentyl-furan, 1-penten-3-one, N-morpholinomethyl-isopropyl-sulfide, methyl butyrate, and (E,E)-2,4-decadienal, are primarily responsible for the distinctive aroma and flavor characteristics associated with pork. Molecular dynamics simulations were employed to investigate the stability of the pre-polymerization system, simulating the interactions between the volatile compounds as templates, 4-hydroxyethyl methacrylate (HEMA) as monomers, and ethylene glycol dimethacrylate (EGDMA) as crosslinkers. Computational simulations revealed that the optimal mole ratio of 1:4:20 for templates, monomers, and crosslinkers resulted in the most favorable functional radial distribution and exhibited the strongest interactions. To validate the computational findings, additional analyses were performed utilizing Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA), radial distribution function (RDF), and hydrogen bond (HBond) occupancy. The calculated binding free energy demonstrated that all template molecules were capable to bind with both the monomers and crosslinkers, including 1-penten-3-one and N-morpholinomethyl-isopropyl-sulfide displaying the strongest interactions, with values of -12,674 kJ/mol and -11,646 kJ/mol, respectively. The congruence between the results obtained from the molecular simulation analyses highlights the crucial role of molecular dynamics simulations in the study and development of MIP for the analysis of marker compounds present in pork.Communicated by Ramaswamy H. Sarma.

Computational studies indicated the effectiveness of human metabolites against SARS-Cov-2 main protease

To fight against the devastating coronavirus disease 2019 (COVID-19), identifying robust anti-SARS-CoV-2 therapeutics from all possible directions is necessary. To contribute to this effort, we selected a human metabolites database containing waters and lipid-soluble metabolites to screen against the 3-chymotrypsin-like proteases (3CLpro) protein of SARS-CoV-2. The top 8 hits from virtual screening displayed a docking score varying between ~ - 11 and ~ - 14 kcal/mol. Molecular dynamics simulations complement the virtual screening study in conjunction with the molecular mechanics generalized Born surface area (MM/GBSA) scheme. Our analyses revealed that (HMDB0132640) has the best glide docking score, - 14.06 kcal/mol, and MM-GBSA binding free energy, - 18.08 kcal/mol. The other three lead molecules are also selected along with the top molecule through a critical inspection of their pharmacokinetic properties. HMDB0132640 displayed a better binding affinity than the other three compounds (HMDB0127868, HMDB0134119, and HMDB0125821) due to increased favorable contributions from the intermolecular electrostatic and van der Waals interactions. Further, we have investigated the ligand-induced structural dynamics of the main protease. Overall, we have identified new compounds that can serve as potential leads for developing novel antiviral drugs against SARS-CoV-2 and elucidated molecular mechanisms of their binding to the main protease. Identification of probable hits from human metabolites against SARS-CoV-2 using integrated computational approaches-Missed against MS.

Finding potential inhibitors against RNA-dependent RNA polymerase (RdRp) of bovine ephemeral fever virus (BEFV): an in-silico study

The bovine ephemeral fever virus (BEFV) is an enzootic agent that affects millions of bovines and causes major economic losses. Though the virus is seasonally reported with a very high morbidity rate (80-100%) from African, Australian, and Asiatic continents, it remains a neglected pathogen in many of its endemic areas, with no proper therapeutic drugs or vaccines presently available for treatment. The RNA-dependent RNA polymerase (RdRp) catalyzes the viral RNA synthesis and is an appropriate candidate for antiviral drug developments. We utilized integrated computational tools to build the 3D model of BEFV-RdRp and then predicted its probable active binding sites. The virtual screening and optimization against these active sites, using several small-molecule inhibitors from a different category of Life Chemical database and FDA-approved drugs from the ZINC database, was performed. We found nine molecules that have docking scores varying between -6.84 to -10.43 kcal/mol. Furthermore, these complexes were analyzed for their conformational dynamics and thermodynamic stability using molecular dynamics simulations in conjunction with the molecular mechanics generalized Born surface area (MM-GBSA) scheme. The binding free energy calculations depict that the electrostatic interactions play a dominant role in the RdRp-inhibitor binding. The hot spot residues, such as Arg565, Asp631, Glu633, Asp740, and Glu707, were found to control the RdRp-inhibitor interaction. The ADMET analysis strongly suggests favorable pharmacokinetics of these compounds that may prove useful for treating the BEFV ailment. Overall, we anticipate that these findings would help explore and develop a wide range of anti-BEFV therapy.Communicated by Ramaswamy H. Sarma.

Finding inhibitors and deciphering inhibitor-induced conformational plasticity in the Janus kinase via multiscale simulations

The Janus kinase (JAK) is a master regulator of the JAK/STAT pathway. Dysregulation of this signalling cascade causes neuroinflammation and autoimmune disorders. Therefore, JAKs have been characterized as an attractive target for developing anti-inflammatory drugs. Nowadays, designing efficient, effective, and specific targeted therapeutics without being cytotoxic has gained interest. We performed the virtual screening of natural products in combination with pharmacological analyses. Subsequently, we performed molecular dynamics simulations to study the stability of the ligand-bound complexes and ligand-induced inactive conformations. Notably, inactive kinases display remarkable conformational plasticity; however, ligand-induced molecular mechanisms of these conformations are still poorly understood. Herein, we performed a free energy landscape analysis to explore the conformational plasticity of the JAK1 kinase. Leonurine, STOCK1N-68642, STOCK1N-82656, and STOCK1N-85809 bound JAK1 exhibited a smooth transition from an active (αC-in) to a completely inactive conformation (αC-out). Ligand binding induces disorders in the αC-helix. Molecular mechanics Poisson Boltzmann surface area (MM/PBSA) calculation suggested three phytochemicals, namely STOCK1N-68642, Epicatechin, and STOCK1N-98615, have higher binding affinity compared to other ligand molecules. The ligand-induced conformational plasticity revealed by our simulations differs significantly from the available crystal structures, which might help in designing allosteric drugs.

Molecular docking used as an advanced tool to determine novel compounds on emerging infectious diseases: A systematic review

Emerging infectious diseases (EID) as well as reappearing irresistible infections are expanding worldwide. Utmost of similar cases, it was seen that the EIDs have long been perceived as a predominant conclusion of host-pathogen adaption. Here, one should get to analyze their host-pathogen interlink and their by needs to look ways, as an example, by exploitation process methodology particularly molecular docking and molecular dynamics simulation, have been utilized in recent time as the most outstanding tools. Hence, we have overviewed some of important factors that influences on EIDs especially HIV/AIDs, H1N1 and coronavirus. Moreover, here we specified the importance of molecular docking applications especially molecular dynamics simulations approach to determine novel compounds on the emerging infectious diseases. Additionally, in vivo and in vitro studies approach to determine novel compounds on the emerging infectious diseases that has implemented to evaluate the limiting affinities between small particles as well as macromolecule that can further, used as a target of HIV/AIDs, H1N1, and coronavirus were also discussed. These novel drug molecules approved in vivo and in vitro studies with reaffirm results and hence, it is clear that the computational methods (mainly molecular docking and molecular dynamics) are found to be more effective technique for drug discovery and medical practitioners.

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.

A plausible contributor to multiple sclerosis; presentation of antigenic myelin protein epitopes by major histocompatibility complexes

BACKGROUND

Multiple sclerosis (MS) can be induced upon successful presentation of myelin antigens by MHC I/II. Antigenic similarity between the myelin and viral proteins may worsen the immunological responses.

METHODOLOGY

Antigenic regions within myelin proteins; PLP1, MBP, MOG, and MAG were analyzed using SVMTrip and EMBOSS. Homology search identified sequence similarity between the predicted host epitopes and viral proteins. NetMHCpan predicted MHC I/II binding followed by peptide-protein docking through the HPEPDOCK server. Thereafter we analyzed conformational flexibility and stability of 15 protein-peptide complexes based on high docking scores. The binding free energy was calculated using conventional (MD) and Gaussian accelerated molecular dynamics simulation.

RESULTS

PLP1, MBP, MAG and MOG contained numerous antigenic epitopes. MBP and MOG epitopes had sequence similarity to HHV-6 BALF5; EBNA1 and CMV glycoprotein M (gM), and EBV LMP2B, gp350/220; HHV-8 ORFs respectively. Many herpes virus proteins like tegument, envelope glycoproteins, and ORFs of EBV, CMV, HHV-6, and HHV-8 demonstrated sequence similarity with MAG and PLP1. Some antigenic peptides were also linear B-cell epitopes and influenced cytokine production by T-cell. MHC I allele HLA-B*57:01 bound to PLP1 peptide and HLA-A*68:02 bound to a MAG peptide strongly. MHC II alleles HLA-DRB1*04:05 and HLA-DR1*01:01 associated with MAG- and MOG-derived peptides, respectively, demonstrating high HPEPDOCK scores. MD simulations established stable binding of certain peptides with the MHC namely HLA-B*51:01-MBP(DYKSAHKGFKGVDAQGTLSKIFKL), HLA-B*57:01-PLP1(PDKFVGITYALTVVWLLVFACSAVPVYIYF), HLA-DR1*01:01-MOG(VEDPFYWVSPGVLVLLAVLPVLLLQITVGLVFLCLQYR) and HLA-DRB1*04:05-MAG(TWVQVSLLHFVPTREA).

CONCLUSIONS

Cross-reactivity between self-antigens and pathogen derived immunodominant epitopes may induce MS. Our study supported the role of specific MHC alleles as a contributing MS risk factor.

Role of Doxorubicin on the Loading Efficiency of ICG within Silk Fibroin Nanoparticles

The effective loading or encapsulation of multimodal theranostic agents within a nanocarrier system plays an important role in the clinical development of cancer therapy. In recent years, the silk fibroin protein-based delivery system has been drawing significant attention to be used in nanomedicines due to its biocompatible and biodegradable nature. In this study, silk fibroin nanoparticles (SNPs) have been synthesized by a novel and cost-effective ultrasonic atomizer-based technique for the first time. The fabricated SNPs were coencapsulated by the FDA-approved indocyanine green (ICG) dye and the chemotherapeutic drug doxorubicin (DOX). The synthesized SNPs are spherical, with an average diameter of ∼37 ± 4 nm, and the ICG-DOX-coencapsulated SNPs (ID-SNPs) have a diameter size of ∼47 ± 6 nm. For the first time, here we demonstrate that DOX helps in the higher loading of ICG within the ID-SNPs, which enhances the encapsulation efficiency of ICG by ∼99%. This could be attributed to the interaction of ICG and DOX molecules with the silk fibroin protein, which helps ICG to get loaded more efficiently within these nanoparticles. The overall finding of this study suggests that the ID-SNPs could be utilized for enhanced ICG-complemented multimodal deep-tissue bioimaging and synergistic chemo-photothermal therapy.

Plant-derived active compounds as a potential nucleocapsid protein inhibitor of SARS-CoV-2: an in-silico study

The coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2. This virus has a high mismatch repair proofreading ability due to its unique exonuclease activity, making it knotty to treat. The nucleocapsid protein can serve as a potential antiviral drug target, as this protein is responsible for multiple captious functions during the viral life cycle. Herein, we have investigated the potential to repurpose active antiviral compounds of plant origins for treating the SARS-CoV-2 infection. In the present study, we followed the molecular docking methodology to screen druggable natural plants' active compounds against the nucleocapsid protein of SARS-CoV-2. The virtual screening of all 68 compounds revealed that the top seven active compounds, such as withanolide D, hypericin, silymarin, oxyacanthine, withaferin A, Acetyl aleuritolic acid, and rhein, exhibit good binding affinity with druggable ADME properties, toxicity, and Pass prediction. The stability of the docked complexes was studied by conducting molecular simulations of 100 ns. MM-GBSA calculated the binding free energy uncovered that withanolide D, hypericin, and silymarin result in highly stable binding conformations in three different sites of the nucleocapsid protein. However, further investigation is needed in order to validate the candidacy of these inhibitors for clinical trials. HighlightsNatural plants' active compounds may aid in the inhibition of SARS-CoV-2 replication and COVID-19 therapeutics.Hypericin, silymarin, withanolide D, oxyacanthine, withaferin A, Acetyl aleuritolic acid, and rhein are effective against SARS-CoV-2 N protein.Studied natural plants' active compounds could be useful against COVID-19 and its associated organs comorbidities.ADMET properties of selected compounds favor these compounds as druggable candidates.Communicated by Ramaswamy H. Sarma.

Unraveling the Molecular Mechanism of Recognition of Selected Next-Generation Antirheumatoid Arthritis Inhibitors by Janus Kinase 1

Rheumatoid arthritis (RA) is a chronic immune-related condition, primarily of joints, and is highly disabling and painful. The inhibition of Janus kinase (JAK)-related cytokine signaling pathways using small molecules is prevalent nowadays. The JAK family belongs to nonreceptor cytoplasmic protein tyrosine kinases (PTKs), including JAK1, JAK2, JAK3, and TYK2 (tyrosine kinase 2). JAK1 has received significant attention after being identified as a promising target for developing anti-RA therapeutics. Currently, no crystal structure is available for JAK1 in complex with the next-generation anti-RA drugs. In the current study, we investigated the mechanism of binding of baricitinib, filgotinib, itacitinib, and upadacitinib to JAK1 using a combined method of molecular docking, molecular dynamics simulation, and binding free energy calculation via the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) scheme. We found that the calculated binding affinity decreases in the order upadacitinib > itacitinib > filgotinib > baricitinib. Due to the increased favorable intermolecular electrostatic contribution, upadacitinib is more selective to JAK1 compared to the other three inhibitors. The cross-correlation and principal component analyses showed that different inhibitor bindings significantly affect the binding site dynamics of JAK1. Furthermore, our studies indicated that the hydrophobic residues and hydrogen bonds from the hinge region (Glu⁹⁵⁷ and Leu⁹⁵⁹) of JAK1 played an essential role in stabilizing the inhibitors. Protein structural network analysis reveals that the total number of links and hubs in JAK1/baricitinib (354, 48) is more significant than those in apo (328, 40) and the other three complexes. The JAK1/baricitinib complex shows the highest probability of the highest-ranked community, indicating a compact network of the JAK1/baricitinib complex, consistent with its higher stability than the rest of the four systems. Overall, our study may be crucial for the rational design of JAK1-selective inhibitors with better affinity.

Discovery of potential competitive inhibitors against With-No-Lysine kinase 1 for treating hypertension by virtual screening, inverse pharmacophore-based lead optimization, and molecular dynamics simulations

The With-No-Lysine (WNK) has received attention because of its involvement in hypertension. Genetic mutation in the genes of WNK, leading to its overexpression, has been reported in Familial Hyperkalaemic Hypertension, and thus WNK is considered a potential drug target. Herein, we have performed a high-throughput virtual screening of ~11,000 compounds, mainly the natural phytochemical compounds and kinase inhibitory libraries, to find potential competitive inhibitors against WNK1. Initially, candidates with a docking score of ~ -10.0 kcal/mol or less were selected to further screen their good pharmacological properties by applying absorption, distribution, metabolism, excretion, and toxicity (ADMET). Finally, six docked compounds bearing appreciable binding affinities and WNK1 selectivity were complimented with 500 ns long all-atom molecular dynamic simulations. Subsequently, the MMPBSA scheme (Molecular Mechanics Poisson Boltzmann Surface Area) suggested three phytochemical compounds, C00000947, C00020451, and C00005031, with favourable binding affinity against WNK1. Among them, C00000947 acts as the most potent competitive inhibitor of WNK1. Further, inverse pharmacophore-based lead optimization of the C00000947 leads to one potential compound, meciadanol, which shows better binding affinity and specificity than C00000947 towards WNK1, which may be further exploited to develop effective therapeutics against WNK1-associated hypertension after in vitro and in vivo validation.

Decoding the Host-Parasite Protein Interactions Involved in Cerebral Malaria Through Glares of Molecular Dynamics Simulations

Malaria causes millions of deaths every year. The malaria parasite spends a substantial part of its life cycle inside human erythrocytes. Inside erythrocytes, it synthesizes and displays various proteins onto the erythrocyte surface, such as Plasmodium falciparum erythrocytic membrane protein-1 (PfEMP1). This protein contains cysteine-rich interdomain region (CIDR) domains which have many subtypes based on sequence diversity and can cross-talk with host molecules. The CIDRα1.4 subtype can attach host endothelial protein C receptor (EPCR). This interaction facilitates infected erythrocyte adherence to brain endothelium and subsequent development of cerebral malaria. Through molecular dynamics simulations in conjunction with the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method, we explored the mechanism of interaction in the CIDRα1-EPCR complex. We examined the structural behavior of two CIDRα1 molecules (encoded by HB3-isolate var03-gene and IT4-isolate var07-gene) with EPCR unbound and bound (complex) forms. ^HB3var03CIDRα1 in apo and complexed with EPCR was comparatively more stable than ^IT4var07CIDRα1. Both of the complexes adopted two distinct conformational energy states. The hydrophobic residues played a crucial role in the binding of both complexes. For ^HB3var03CIDRα1-EPCR, the dominant energetic components were total polar interactions, while in ^IT4var07CIDRα1-EPCR, the primary interaction was van der Waals and nonpolar solvation energy. The study also revealed details such as correlated conformational motions and secondary structure evolution. Further, it elucidated various hotspot residues involved in protein-protein recognition. Overall, our study provides additional information on the structural behavior of CIDR molecules in unbound and receptor-bound states, which will help to design potent inhibitors.

Unraveling the Molecular Mechanism of Recognition of Human Interferon-Stimulated Gene Product 15 by Coronavirus Papain-Like Proteases: A Multiscale Simulation Study

The papain-like protease (PLpro) of the coronavirus (CoV) family plays an essential role in processing the viral polyprotein and immune evasion. Additional proteolytic activities of PLpro include deubiquitination and deISGylation, which can reverse the post-translational modification of cellular proteins conjugated with ubiquitin or (Ub) or Ub-like interferon-stimulated gene product 15 (ISG15). These activities regulate innate immune responses against viral infection. Thus, PLpro is a potential antiviral target. Here, we have described the structural and energetic basis of recognition of PLpro by the human ISG15 protein (hISG15) using atomistic molecular dynamics simulation across the CoV family, i.e., MERS-CoV (MCoV), SARS-CoV (SCoV), and SARS-CoV-2 (SCoV2). The cumulative simulation length for all trajectories was 32.0 μs. In the absence of the complete crystal structure of complexes, protein-protein docking was used. A mutation (R167E) was introduced across all three PLpro to study the effect of mutation on the protein-protein binding. Our study reveals that the apo-ISG15 protein remains closed while it adopts an open conformation when bound to PLpro, although the degree of openness varies across the CoV family. The binding free energy analysis suggests that hISG15 binds more strongly with SCoV2-PLpro compared to SCoV or MCoV. The intermolecular electrostatic interaction drives the hISG15-PLpro complexation. Our study showed that SCoV or MCoV-PLpro binds more strongly with the C-domain of hISG15, while SCoV2-PLpro binds more favorably the N-domain of hISG15. Overall, our study explains the molecular basis of differential deISGylating activities of PLpro among the CoV family and the specificity of SCoV2-PLpro toward hISG15.

Investigating potency of TMC-126 against wild-type and mutant variants of HIV-1 protease: a molecular dynamics and free energy study

A detailed computational study was performed to investigate the conformational changes of flap region and the mechanism underlying the binding of the inhibitor TMC-126 to HIV-1 protease (PR1) and its mutant variants through molecular dynamics simulations in conjunction with the molecular mechanics Poisson-Boltzmann (MM-PBSA) free energy calculation. Further, we have studied the effectiveness of the inhibitor against HIV-2 protease (PR2). The MM-PBSA calculation suggests that TMC-126 loses its potency against mutant variants and PR2 compared to wild-type PR1 mainly due to the loss in intermolecular electrostatic interactions. The potency of the inhibitor decreases in the order: wild type PR1 > M46L > MDR20 > I50V > PR2 > V32I > A28S. Our study reveals that the flap of PR1 adopts a semi-open conformation due to the mutation I50V or MDR20. The dissimilar nature of the movement of the flap tip of both monomers is evident from the dynamic cross-correlation map. The protein structural network analysis displays that mutation causes structural rearrangements and changes the communication path between residues. Overall, we believe our study may help explore and accelerate the development of novel HIV-1/HIV-2 protease inhibitors with better potency.

Molecular Mechanism of Inhibiting WNK Binding to OSR1 by Targeting the Allosteric Pocket of the OSR1-CCT Domain with Potential Antihypertensive Inhibitors: An In Silico Study

The oxidative-stress-responsive kinase 1 (OSR1) and the STE20/SPS1-related proline-alanine-rich kinase (SPAK) are physiological substrates of the with-no-lysine (WNK) kinase. They are the master regulators of cation Cl^- cotransporters that could be targeted for discovering novel antihypertensive agents. Both kinases have a conserved carboxy-terminal (CCT) domain that recognizes a unique peptide motif (Arg-Phe-Xaa-Val) present in their upstream kinases and downstream substrates. Here, we have combined molecular docking with molecular dynamics simulations and free-energy calculations to identify potential inhibitors that can bind to the allosteric pocket of the OSR1-CCT domain and impede its interaction with the WNK peptide. Our study revealed that STOCK1S-14279 and Closantel bound strongly to the allosteric pocket of OSR1 and displaced the WNK peptide from the primary pocket of OSR1. We showed that primarily Arg1004 and Gln1006 of the WNK4-peptide motif were involved in strong H-bond interactions with Glu453 and Arg451 of OSR1. Besides, our study revealed that atoms of Arg1004 were solvent-exposed in cases of STOCK1S-14279 and Closantel, implying that the WNK4 peptide was moved out of the pocket. Overall, the predicted potential inhibitors altogether abolish the OSR1-WNK4-peptide interaction, suggesting their potency as a prospective allosteric inhibitor against OSR1.

Computational Investigation of Structural Dynamics of SARS-CoV-2 Methyltransferase-Stimulatory Factor Heterodimer nsp16/nsp10 Bound to the Cofactor SAM

Recently, a highly contagious novel coronavirus disease 2019 (COVID-19), caused by SARS-CoV-2, has emerged, posing a global threat to public health. Identifying a potential target and developing vaccines or antiviral drugs is an urgent demand in the absence of approved therapeutic agents. The 5'-capping mechanism of eukaryotic mRNA and some viruses such as coronaviruses (CoVs) are essential for maintaining the RNA stability and protein translation in the virus. SARS-CoV-2 encodes S-adenosyl-L-methionine (SAM) dependent methyltransferase (MTase) enzyme characterized by nsp16 (2'-O-MTase) for generating the capped structure. The present study highlights the binding mechanism of nsp16 and nsp10 to identify the role of nsp10 in MTase activity. Furthermore, we investigated the conformational dynamics and energetics behind the binding of SAM to nsp16 and nsp16/nsp10 heterodimer by employing molecular dynamics simulations in conjunction with the Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) method. We observed from our simulations that the presence of nsp10 increases the favorable van der Waals and electrostatic interactions between SAM and nsp16. Thus, nsp10 acts as a stimulator for the strong binding of SAM to nsp16. The hydrophobic interactions were predominately identified for the nsp16-nsp10 interactions. Also, the stable hydrogen bonds between Ala83 (nsp16) and Tyr96 (nsp10), and between Gln87 (nsp16) and Leu45 (nsp10) play a vital role in the dimerization of nsp16 and nsp10. Besides, Computational Alanine Scanning (CAS) mutagenesis was performed, which revealed hotspot mutants, namely I40A, V104A, and R86A for the dimer association. Hence, the dimer interface of nsp16/nsp10 could also be a potential target in retarding the 2'-O-MTase activity in SARS-CoV-2. Overall, our study provides a comprehensive understanding of the dynamic and thermodynamic process of binding nsp16 and nsp10 that will contribute to the novel design of peptide inhibitors based on nsp16.

Molecular dynamics simulation of docking structures of SARS-CoV-2 main protease and HIV protease inhibitors

•

We consider possible repurposed-drugs candidates against SARS-CoV-2.

•

10 different HIV protease inhibitors were investigated.

•

In silico simulations were used to study protease inhibitors for SARS-CoV-2.

In this paper we investigate 10 different HIV protease inhibitors (HPIs) as possible repurposed-drugs candidates against SARS-CoV-2. To this end, we execute molecular docking and molecular dynamics simulations. The in silico data demonstrated that, despite their molecular differences, all HPIs presented a similar behavior for the parameters analyzed, with the exception of Nelfinavir that showed better results for most of the molecular dynamics parameters in comparison with the N3 inhibitor.

Structural dynamics and susceptibility of anti-HIV drugs against HBV reverse transcriptase

Hepatitis B virus (HBV), a small enveloped DNA virus, attacks the human liver causing both acute and chronic diseases. Current therapeutic drugs use the nucleos(t)ide analogues (NAs) as a competitive inhibitor against HBV reverse transcriptase (HBV-RT), an essential enzyme pivotally involved in viral replication. Unfortunately, this treatment still causes the development of resistant variants of HBV against NAs. As HBV-RT is homologous to the human immunodeficiency virus reverse transcriptase (HIV-RT), it is reasonable to treat HBV-RT with anti-HIV drugs. In the present study, we aimed to investigate the structural dynamics and susceptibility of the known anti-HIV drugs (stavudine [d4T], didanosine [DDI], and zidovudine [ZDV]) against HBV-RT enzyme in comparison to the anti-HBV drug lamivudine (3TC) and deoxythymidine triphosphate (dTTP) substrate using several computational approaches. The ΔGbindresidue calculations revealed that seven polar residues (K32, R41, D83, S85, D205, N236, and K239) and three hydrophobic residues (A86, A87, and F88) of HBV-RT as well as the adjacent DNA strands play an important role in the ligand binding. In addition, the H-bond pattern of d4T is similar to that of 3TC, especially at the residues A86 and A87. Such interactions promote the favorable conformation of ligand in the HBV-RT binding pocket, while the several different conformations of ligand are found in the unbound state. The predicted binding free energy results based on QM/MM-GBSA and MM/GB(PB)SA methods suggested that the susceptibility towards HBV-RT of d4T and ZDV is higher than that of 3TC and dTTP. Altogether, this work sheds light on the potentiality of d4T and ZDV as a promising drug for HBV-infected patients harboring 3TC resistance.Communicated by Ramaswamy H. Sarma.

Investigating Conformational Dynamics of Lewis Y Oligosaccharides and Elucidating Blood Group Dependency of Cholera Using Molecular Dynamics

Cholera is caused by Vibrio cholerae and is an example of a blood-group-dependent disease. Recent studies suggest that the receptor-binding B subunit of the cholera toxin (CT) binds histo-blood group antigens at a secondary binding site. Herein, we studied the conformational dynamics of Lewis Y (Le^Y) oligosaccharides, H-tetrasaccharides and A-pentasaccharides, in aqueous solution by conducting accelerated molecular dynamics (aMD) simulations. The flexible nature of both oligosaccharides was displayed in aMD simulations. Furthermore, aMD simulations revealed that for both oligosaccharides in the free form, ⁴C₁ and ¹C₄ puckers were sampled for all but GalNAc monosaccharides, while either the ⁴C₁ (GlcNAc, Gal, GalNAc) or ¹C₄ (Fuc2, Fuc3) pucker was sampled in the CT-bound forms. In aMD, the complete transition from the ⁴C₁ to ¹C₄ pucker was sampled for GlcNAc and Gal in both oligosaccharides. Further, we have observed a transition from the open to closed conformer in the case of A-pentasaccharide, while H-tetrasaccharide remains in the open conformation throughout the simulation. Both oligosaccharides adopted an open conformation in the CT binding site. Moreover, we have investigated the molecular basis of recognition of Le^Y oligosaccharides by the B subunit of the cholera toxin of classical and El Tor biotypes using the molecular mechanics generalized Born surface area (MM/GBSA) scheme. The O blood group determinant, H-tetrasaccharide, exhibits a stronger affinity to both biotypes compared to the A blood group determinant, A-pentasaccharide, which agrees with the experimental data. The difference in binding free energy between O and A blood group determinants mainly arises due to the increased entropic cost and desolvation energy in the case of A-pentasaccharide compared to that of H-tetrasaccharide. Our study also reveals that the terminal Fuc3 contributes most to the binding free energy compared to other carbohydrate residues as it forms multiple hydrogen bonds with CT. Overall, our study might help in designing glycomimetic drugs targeting the cholera toxin.

Computational investigation on the effect of Oleuropein aglycone on the α-synuclein aggregation

Parkinson's disease (PD) is considered to be the second most common progressive neurodegenerative brain disorder after Alzheimer's disease, which is caused by misfolding and aggregation of Alpha-synuclein (α-synuclein). It is characterized by distinct aggregated fibrillary form of α-synuclein known as the Lewy bodies and Lewy neurites. The most promising approach to combat PD is to prevent the misfolding and subsequent aggregation of α-synuclein. Recently, Oleuropein aglycone (OleA) has been reported to stabilize the monomeric structure of α-synuclein, subsequently favoring the growth of nontoxic aggregates. Therefore, understanding the conformational dynamics of α-synuclein monomer in the presence of OleA is significant. Here, we have investigated the effect of OleA on the conformational dynamics and the aggregation propensity of α-synuclein using molecular dynamics simulation. From molecular dynamics trajectory analysis, we noticed that when OleA is bound to α-synuclein, the intramolecular distance between non-amyloid-β component domain and C-terminal domain of α-synuclein was increased, whereas long-range hydrophobic interactions between the two region were reduced. Oleuropein aglycone was found to interact with the N-terminal domain of α-synuclein, making this region unavailable for interaction with membranes and lipids for the formation of cellular toxic aggregates. From the binding-free energy analysis, we found binding affinity between α-synuclein and OleA to be indeed high (ΔG_bind = -12.56 kcal mol^-1 from MM-PBSA and ΔG_bind = -27.41 kcal mol^-1from MM-GBSA). Our findings in this study thus substantiate the effect of OleA on the structure and stabilization of α-synuclein monomer that subsequently favors the growth of stable and nontoxic aggregates.Communicated by Ramaswamy H. Sarma.

Effect of CTerm of human albumin on the aggregation propensity of Aβ1-42 peptide: a potential of mean force study

Alzheimer's disease (AD) is the most common progressive neurodegenerative brain disorder. It is characterized by the presence of extracellular aggregated fibrillary form of amyloid beta (Aβ) peptide and intraneuronal neurofibrillary tangles caused by the hyperphosphorylation of tau protein. Monomeric form of Aβ peptide in α-conformation is not toxic but it can undergo self-aggregation to form β-conformation which is neurotoxic. The most promising approach to combat AD is to prevent the self-aggregation of Aβ peptide. Recently, it has been reported that C-terminal (CTerm) of human albumin (HA) binds to the Aβ_1-42 peptide and impairs the Aβ_1-42 aggregation and promotes disassembly of Aβ_1-42 aggregates. In this work, using potential of mean force (PMF) and binding free energy (BFE) calculations, we have demonstrated the effect of CTerm of HA on the dimerization of Aβ_1-42 peptide. From the PMF profile, we noticed Aβ_1-42-CTerm Heterodimer (10.99 kcal mol - 1) complex to have higher disassociation energy than Aβ_1-42-Aβ_1-42 homodimer (2.23 kcal mol - 1) complex. And also from the BFE calculations, we found that the binding affinity between Aβ_1-42 peptide and CTerm (ΔG_bind = -32.27 kcal mol - 1 from MM-GBSA and ΔG_bind = -2.83 kcal mol - 1 from MM-PBSA (molecular mechanics-Poisson - Boltzmann surface area)) to be stronger than the Aβ_1-42 peptide and another Aβ_1-42 peptide (ΔG_bind = -16.20 kcal mol - 1 from MM-GBSA and ΔG_bind = -1.95 kcal mol - 1 from MM-PBSA). In this study, our findings from PMF and BFE analysis of the two complexes provide salient structural, binding and unbinding features and thermodynamics that support the ability of CTerm of HA in affecting the dimerization of Aβ_1-42. Communicated by Ramaswamy H. Sarma.