Investigating the mechanism of recognition and structural dynamics of nucleoprotein-RNA complex from Peste des petits ruminants virus via Gaussian accelerated molecular dynamics simulations

Identification of inhibitors for neurodegenerative diseases targeting dual leucine zipper kinase through virtual screening and molecular dynamics simulations

Neurodegenerative diseases lead to a gradual decline in cognitive and motor functions due to the progressive loss of neurons in the central nervous system. The role of dual leucine zipper kinase (DLK) in regulating stress responses and neuronal death pathways highlights its significance as a target against neurodegenerative diseases. The non-availability of FDA-approved drugs emphasizes a need to identify novel DLK-inhibitors. We screened NPAtlas (Natural products) and MedChemExpress (FDA-approved) libraries to identify potent ATP-competitive DLK inhibitors. ADMET analyses identified four compounds (two natural products and two FDA-approved) with favourable features. Subsequently, we performed molecular dynamics simulations to examine the binding-stability and ligand-induced conformational dynamics. Molecular mechanics Poisson Boltzmann surface area (MM-PBSA) calculations demonstrated CID139591660, dithranol, and danthron having greater affinity, while CID156581477 showed lower affinity than control sunitinib. PCA and network analysis results indicated structural and network alteration post-ligand binding. Furthermore, we identified an analogue of CID156581477 using the deep learning-based web server DeLA Drug which demonstrated a higher affinity than its parent compound and the control and identified several crucial interacting residues. Overall, our study provides significant theoretical guidance for designing potent novel DLK inhibitors and compounds that could emerge as promising drug candidates against DLK following laboratory validation.

Protein S-palmitoylation is markedly inhibited by 4″-alkyl ether lipophilic derivatives of EGCG, the major green tea polyphenol: In vitro and in silico studies

S-palmitoylation is a dynamic lipid-based protein post-translational modification facilitated by a family of protein acyltransferases (PATs) commonly known as DHHC-PATs or DHHCs. It is the only lipid modification that is reversible, and this very fact uniquely qualifies it for therapeutic interventions through the development of DHHC inhibitors. Herein, we report that 4″-alkyl ether lipophilic derivatives of EGCG can effectively inhibit protein S-palmitoylation in vitro. With the help of metabolic labeling followed by copper(I)-catalyzed azide-alkyne cycloaddition Click reaction, we demonstrate that 4″-C14 EGCG and 4″-C16 EGCG markedly inhibited S-palmitoylation in various mammalian cells including HEK 293T, HeLa, and MCF-7 using both in gel fluorescence as well as confocal microscopy. Further, these EGCG derivatives were able to attenuate the S-palmitoylation to the basal level in DHHC3-overexpressed cells, suggesting that they are plausibly targeting DHHCs. Confocal microscopy data qualitatively reflected spatial and temporal distribution of S-palmitoylated proteins in different sub-cellular compartments and the inhibitory effects of 4″-C14 EGCG and 4″-C16 EGCG were clearly observed in the native cellular environment. Our findings were further substantiated by in silico analysis which revealed promising binding affinity and interactions of 4″-C14 EGCG and 4″-C16 EGCG with key amino acid residues present in the hydrophobic cleft of the DHHC20 enzyme. We also demonstrated the successful inhibition of S-palmitoylation of GAPDH by 4″-C16 EGCG. Taken together, our in vitro and in silico data strongly suggest that 4″-C14 EGCG and 4″-C16 EGCG can act as potent inhibitors for S-palmitoylation and can be employed as a complementary tool to investigate S-palmitoylation.

Deciphering the molecular choreography of Janus kinase 2 inhibition via Gaussian accelerated molecular dynamics simulations: a dynamic odyssey

The Janus kinases (JAK) are crucial targets in drug development for several diseases. However, accounting for the impact of possible structural rearrangements on the binding of different kinase inhibitors is complicated by the extensive conformational variability of their catalytic kinase domain (KD). The dynamic KD contains mainly four prominent mobile structural motifs: the phosphate-binding loop (P-loop), the αC-helix within the N-lobe, the Asp-Phe-Gly (DFG) motif, and the activation loop (A-loop) within the C-lobe. These distinct structural orientations imply a complex signal transmission path for regulating the A-loop's flexibility and conformational preference for optimal JAK function. Nevertheless, the precise dynamical features of the JAK induced by different types of inhibitors still remain elusive. We performed comparative, microsecond-long, Gaussian accelerated molecular dynamics simulations in triplicate of three phosphorylated JAK2 systems: the KD alone, type-I ATP-competitive inhibitor (CI) bound KD in the catalytically active DFG-in conformation, and the type-II inhibitor (AI) bound KD in the catalytically inactive DFG-out conformation. Our results indicate significant conformational variations observed in the A-loop and αC helix motions upon inhibitor binding. Our studies also reveal that the DFG-out inactive conformation is characterized by the closed A-loop rearrangement, open catalytic cleft of N and C-lobe, the outward movement of the αC helix, and open P-loop states. Moreover, the outward positioning of the αC helix impacts the hallmark salt bridge formation between Lys882 and Glu898 in an inactive conformation. Finally, we compared their ligand binding poses and free energy by the MM/PBSA approach. The free energy calculations suggested that the AI's binding affinity is higher than CI against JAK2 due to an increased favorable contribution from the total non-polar interactions and the involvement of the αC helix. Overall, our study provides the structural and energetic insights crucial for developing more promising type I/II JAK2 inhibitors for treating JAK-related diseases.

4″-Alkyl EGCG Derivatives Induce Cytoprotective Autophagy Response by Inhibiting EGFR in Glioblastoma Cells

Epidermal growth factor receptor (EGFR)-targeted therapy has been proven vital in the last two decades for the treatment of multiple cancer types, including nonsmall cell lung cancer, glioblastoma, breast cancer and head and neck squamous cell carcinoma. Unfortunately, the majority of approved EGFR inhibitors fall into the drug resistance category because of continuous mutations and acquired resistance. Recently, autophagy has surfaced as one of the emerging underlying mechanisms behind resistance to EGFR-tyrosine kinase inhibitors (TKIs). Previously, we developed a series of 4″-alkyl EGCG (4″-Cn EGCG, n = 6, 8, 10, 12, 14, 16, and 18) derivatives with enhanced anticancer effects and stability. Therefore, the current study hypothesized that 4″-alkyl EGCG might induce cytoprotective autophagy upon EGFR inhibition, and inhibition of autophagy may lead to improved cytotoxicity. In this study, we have observed growth inhibition and caspase-3-dependent apoptosis in 4″-alkyl EGCG derivative-treated glioblastoma cells (U87-MG). We also confirmed that 4″-alkyl EGCG could inhibit EGFR in the cells, as well as mutant L858R/T790M EGFR, through an in vitro kinase assay. Furthermore, we have found that EGFR inhibition with 4″-alkyl EGCG induces cytoprotective autophagic responses, accompanied by the blockage of the AKT/mTOR signaling pathway. In addition, cytotoxicity caused by 4″-C10 EGCG, 4″-C12 EGCG, and 4″-C14 EGCG was significantly increased after the inhibition of autophagy by the pharmacological inhibitor chloroquine. These findings enhance our understanding of the autophagic response toward EGFR inhibitors in glioblastoma cells and suggest a potent combinatorial strategy to increase the therapeutic effectiveness of EGFR-TKIs.

Dynamics and Function of sRNA/mRNAs Under the Scrutiny of Computational Simulation Methods

Molecular dynamics simulations have proved extremely useful in investigating the functioning of proteins with atomic-scale resolution. Many applications to the study of RNA also exist, and their number increases by the day. However, implementing MD simulations for RNA molecules in solution faces challenges that the MD practitioner must be aware of for the appropriate use of this tool. In this chapter, we present the fundamentals of MD simulations, in general, and the peculiarities of RNA simulations, in particular. We discuss the strengths and limitations of the technique and provide examples of its application to elucidate small RNA's performance.

Roles of Accelerated Molecular Dynamics Simulations in Predictions of Binding Kinetic Parameters

Rational predictions on binding kinetics parameters of drugs to targets play significant roles in future drug designs. Full conformational samplings of targets are requisite for accurate predictions of binding kinetic parameters. In this review, we mainly focus on the applications of enhanced sampling technologies in calculations of binding kinetics parameters and residence time of drugs. The methods involved in molecular dynamics simulations are applied to not only probe conformational changes of targets but also reveal calculations of residence time that is significant for drug efficiency. For this review, special attention are paid to accelerated molecular dynamics (aMD) and Gaussian aMD (GaMD) simulations that have been adopted to predict the association or disassociation rate constant. We also expect that this review can provide useful information for future drug design.

Development of chalcone-like derivatives and their biological and mechanistic investigations as novel influenza nuclear export inhibitors

Concerning the emergence of resistance to current anti-influenza drugs, our previous phenotypic-based screening study identified the compound A9 as a promising lead compound. This chalcone analog, containing a 2,6-dimethoxyphenyl moiety, exhibited significant inhibitory activity against oseltamivir-resistant strains (H1N1 pdm09), with an EC50 value of 1.34 μM. However, it also displayed notable cytotoxicity, with a CC50 value of 41.46 μM. Therefore, compound A9 was selected as a prototype structure for further structural optimization in this study. Initially, it was confirmed that the substituting the α,β-unsaturated ketone with pent-1,4-diene-3-one as a linker group significantly reduced the cytotoxicity of the final compounds. Subsequently, the penta-1,4-dien-3-one group was utilized as a privileged fragment for further structural optimization. Following two subsequent rounds of optimizations, we identified compound IIB-2, which contains a 2,6-dimethoxyphenyl- and 1,4-pentadiene-3-one moieties. This compound exhibited inhibitory effects on oseltamivir-resistant strains comparable to its precursor (compound A9), while demonstrating reduced toxicity (CC50 > 100 μM). Furthermore, we investigated its mechanism of action against anti-influenza virus through immunofluorescence, Western blot, and surface plasmon resonance (SPR) experiments. The results revealed that compound IIB-2 can impede virus proliferation by blocking the export of influenza virus nucleoprotein. Thusly, our findings further emphasize influenza nuclear export as a viable target for designing novel chalcone-like derivatives with potential inhibitory properties that could be explored in future lead optimization studies.

Activation-induced cytidine deaminase an antibody diversification enzyme interacts with chromatin modifier UBN1 in B-cells

Activation-induced cytidine deaminase (AID) is the key mediator of antibody diversification in activated B-cells by the process of somatic hypermutation (SHM) and class switch recombination (CSR). Targeting AID to the Ig genes requires transcription (initiation and elongation), enhancers, and its interaction with numerous factors. Furthermore, the HIRA chaperon complex, a regulator of chromatin architecture, is indispensable for SHM. The HIRA chaperon complex consists of UBN1, ASF1a, HIRA, and CABIN1 that deposit H3.3 onto the DNA, the SHM hallmark. We explored whether UBN1 interacts with AID using computational and in-vitro experiments. Interestingly, our in-silico studies, such as molecular docking and molecular dynamics simulation results, predict that AID interacts with UBN1. Subsequently, co-immunoprecipitation and pull-down experiments established interactions between UBN1 and AID inside B-cells. Additionally, a double immunofluorescence assay confirmed that AID and UBN1 were co-localized in the human and chicken B-cell lines. Moreover, proximity ligation assay studies validated that AID interacts with UBN1. Ours is the first report on the interaction of genome mutator enzyme AID with UBN1. Nevertheless, the fate of interaction between UBN1 and AID is yet to be explored in the context of SHM or CSR.

Molecular level investigation for identifying potential inhibitors against thymidylate kinase of monkeypox through in silico approaches

The need for more advanced and effective monkeypox (Mpox) treatments has become evident with numerous Mpox virus (MPXV) outbreaks. Over the years, interest has increased in developing targeted medicines that are efficient, safe, and precise while avoiding adverse effects. Here, we screened 32409 compounds against thymidylate kinase (TMPK), an emerging target for Mpox treatment. We studied their pharmacological characteristics and analyzed those through all-atom molecular dynamics simulations followed by molecular mechanics Poisson Boltzmann surface area (MM-PBSA) based free energy calculations. According to our findings, the leads CID40777874 and CID28960001 had the highest binding affinities towards TMPK with ΔGbind of -8.04 and -5.58 kcal/mol, respectively, which outperformed our control drug cidofovir (ΔGbind = -2.92 kcal/mol) in terms of binding favourability. Additionally, we observed crucial TMPK dynamics brought on by ligand-binding and identified key residues such as Phe68 and Tyr101 as the critical points of the protein-ligand interaction. The DCCM analysis revealed the role of ligand binding in stabilizing TMPK's binding region, as indicated by residual correlation motions. Moreover, the PSN analysis revealed that the interaction with ligand induces changes in residual network properties, enhancing the stability of complexes. We successfully identified novel compounds that may serve as potential building blocks for constructing contemporary antivirals against MPXV and highlighted the molecular mechanisms underlying their binding with TMPK. Overall, our findings will play a significant role in advancing the development of new therapies against Mpox and facilitating a comprehensive understanding of their interaction patterns.Communicated by Ramaswamy H. Sarma.

Plant derived active compounds of ayurvedic neurological formulation, Saraswatharishta as a potential dual leucine zipper kinase inhibitor: an in-silico study

Recent findings have highlighted the essential role of dual leucine zipper kinase (DLK) in neuronal degeneration. Saraswatharishta (SWRT), an ayurvedic formulation utilized in traditional Indian medicine, has demonstrated effectiveness in addressing neurodegenerative diseases. Herein, we aim to delve into the atomistic details of the mode of action of phytochemicals present in SWRT against DLK. Our screening process encompassed over 500 distinct phytochemicals derived from the main ingredients of the SWRT formulation. Through a comparative analysis of docking scores and relative poses, we successfully identified four novel compounds, which underwent further investigation via 2 × 500 ns long molecular dynamics (MD) simulations. Among the top four compounds, CID16066851 sourced from the Acorus calamus displayed the most stable complex with DLK. The molecular mechanics Poisson - Boltzmann surface area (MM-PBSA) calculations highlighted the significance of electrostatic and van der Waals interactions in the binding recognition process. Additionally, we identified key residues, namely Phe192, Leu243, Val139, and Leu141, as hotspots that predominantly govern the DLK-inhibitor interaction. Notably, the leading compounds are sourced from the Acorus calamus, Syzygium aromaticum, Zingiber officinale, and Anethum sowa plants present in the SWRT formulation. Overall, the findings of our study hold promise for future drug development endeavors combating neurodegenerative conditions.Communicated by Ramaswamy H. Sarma.

Flexible Gaussian Accelerated Molecular Dynamics to Enhance Biological Sampling

Molecular dynamics simulations often struggle to obtain sufficient sampling to study complex molecular events due to high energy barriers separating the minima of interest. Multiple enhanced sampling techniques have been developed and improved over the years to tackle this issue. Gaussian accelerated molecular dynamics (GaMD) is a recently developed enhanced sampling technique that works by adding a biasing potential, lifting the energy landscape up, and decreasing the height of its barriers. GaMD allows one to increase the sampling of events of interest without the need of a priori knowledge of the system or the relevant coordinates. All required acceleration parameters can be obtained from a previous search run. Upon its development, several improvements for the methodology have been proposed, among them selective GaMD in which the boosting potential is selectively applied to the region of interest. There are currently four selective GaMD methods that have shown promising results. However, all of these methods are constrained on the number, location, and scenarios in which this selective boosting potential can be applied to ligands, peptides, or protein-protein interactions. In this work, we showcase a GROMOS implementation of the GaMD methodology with a fully flexible selective GaMD approach that allows the user to define, in a straightforward way, multiple boosting potentials for as many regions as desired. We show and analyze the advantages of this flexible selective approach on two previously used test systems, the alanine dipeptide and the chignolin peptide, and extend these examples to study its applicability and potential to study conformational changes of glycans and glycosylated proteins.

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.

Predicting the effects of rare genetic variants on oncogenic signaling pathways: A computational analysis of HRAS protein function

The HRAS gene plays a crucial role in regulating essential cellular processes for life, and this gene's misregulation is linked to the development of various types of cancers. Nonsynonymous single nucleotide polymorphisms (nsSNPs) within the coding region of HRAS can cause detrimental mutations that disrupt wild-type protein function. In the current investigation, we have employed in-silico methodologies to anticipate the consequences of infrequent genetic variations on the functional properties of the HRAS protein. We have discovered a total of 50 nsSNPs, of which 23 were located in the exon region of the HRAS gene and denoting that they were expected to cause harm or be deleterious. Out of these 23, 10 nsSNPs ([G60V], [G60D], [R123P], [D38H], [I46T], [G115R], [R123G], [P11OL], [A59L], and [G13R]) were identified as having the most delterious effect based on results of SIFT analysis and PolyPhen2 scores ranging from 0.53 to 69. The DDG values -3.21 kcal/mol to 0.87 kcal/mol represent the free energy change associated with protein stability upon mutation. Interestingly, we identified that the three mutations (Y4C, T58I, and Y12E) were found to improve the structural stability of the protein. We performed molecular dynamics (MD) simulations to investigate the structural and dynamic effects of HRAS mutations. Our results showed that the stable model of HRAS had a significantly lower energy value of -18756 kj/mol compared to the initial model of -108915 kj/mol. The RMSD value for the wild-type complex was 4.40 Å, and the binding energies for the G60V, G60D, and D38H mutants were -107.09 kcal/mol, -109.42 kcal/mol, and -107.18 kcal/mol, respectively as compared to wild-type HRAS protein had -105.85 kcal/mol. The result of our investigation presents convincing corroboration for the potential functional significance of nsSNPs in augmenting HRAS expression and adding to the activation of malignant oncogenic signalling pathways.

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.

Structure-based design and synthesis of a novel long-chain 4''-alkyl ether derivative of EGCG as potent EGFR inhibitor: in vitro and in silico studies

Herein, we report the discovery of a novel long-chain ether derivative of (-)-epigallocatechin-3-gallate (EGCG), a major green tea polyphenol as a potent EGFR inhibitor. A series of 4''-alkyl EGCG derivatives have been synthesized via regio-selectively alkylating the 4'' hydroxyl group in the D-ring of EGCG and tested for their antiproliferative activities against high (A431), moderate (HeLa), and low (MCF-7) EGFR-expressing cancer cell lines. The most potent compound, 4''-C14 EGCG showed the lowest IC50 values across all the tested cell lines. 4''-C14 EGCG was also found to be significantly more stable than EGCG under physiological conditions (PBS at pH 7.4). Further western blot analysis and imaging data revealed that 4''-C14 EGCG induced cell death in A431 cells with shrunken nuclei, nuclear fragmentation, membrane blebbing, and increased population of apoptotic cells where BAX upregulation and BCLXL downregulation were observed. In addition, autophosphorylation of EGFR and its downstream signalling proteins Akt and ERK were markedly inhibited by 4''-C14 EGCG. MD simulation and the MM/PBSA analysis disclosed the binding mode of 4''-C14 EGCG in the ATP-binding site of EGFR kinase domain. Taken together, our findings demonstrate that 4''-C14 EGCG can act as a promising potent EGFR inhibitor with enhanced stability.

Effect of Sulfation on the Conformational Dynamics of Dermatan Sulfate Glycosaminoglycan: A Gaussian Accelerated Molecular Dynamics Study

Glycosaminoglycans (GAGs) are anionic biopolymers present on cell surfaces as a part of proteoglycans. The biological activities of GAGs depend on the sulfation pattern. In our study, we have considered three octadecasaccharide dermatan sulfate (DS) chains with increasing order of sulfation (dp6s, dp7s, and dp12s) to illuminate the role of sulfation on the GAG units and its chain conformation through 10 μs-long Gaussian accelerated molecular dynamics simulations. DS is composed of repeating disaccharide units of iduronic acid (IdoA) and N-acetylgalactosamine (N-GalNAc). Here, N-GalNAc is linked to IdoA via β(1-4), while IdoA is linked to N-GalNAc through α(1-3). With the increase in sulfation, the DS structure becomes more rigid and linear, as is evident from the distribution of root-mean-square deviations (RMSDs) and end-to-end distances. The tetrasaccharide linker region of the main chain shows a rigid conformation in terms of the glycosidic linkage. We have observed that upon sulfation (i.e., dp12s), the ring flip between two chair forms vanished for IdoA. The dynamic cross-correlation analysis reveals that the anticorrelation motions in dp12s are reduced significantly compared to dp6s or dp7s. An increase in sulfation generates relatively more stable hydrogen-bond networks, including water bridging with the neighboring monosaccharides. Despite the favorable linear structures of the GAG chains, our study also predicts few significant bendings related to the different puckering states, which may play a notable role in the function of the DS. The relation between the global conformation with the micro-level parameters such as puckering and water-mediated hydrogen bonds shapes the overall conformational space of GAGs. Overall, atomistic details of the DS chain provided in this study will help understand their functional and mechanical roles, besides developing new biomaterials.

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.

Comparative Structural Dynamics of Isoforms of Helicobacter pylori Adhesin BabA Bound to Lewis b Hexasaccharide via Multiple Replica Molecular Dynamics Simulations

BabA of Helicobacter pylori is the ABO blood group antigen-binding adhesin. Despite considerable diversity in the BabA sequence, it shows an extraordinary adaptation in attachment to mucosal layers. In the current study, multiple replica molecular dynamics simulations were conducted in a neutral aqueous solution to elucidate the conformational landscape of isoforms of BabA bound to Lewis b (Le^b) hexasaccharide. In addition, we also investigated the underlying molecular mechanism of the BabA-glycan complexation using the MM/GBSA scheme. The conformational dynamics of Le^b in the free and protein-bound states were also studied. The carbohydrate-binding site across the four isoforms was examined, and the conformational variability of several vital loops was observed. The cysteine–cysteine loops and the two diversity loops (DL1 and DL2) were identified to play an essential role in recognizing the glycan molecule. The flexible crown region of BabA was stabilized after association with Le^b. The outward movement of the DL2 loop vanished upon ligand binding for the Spanish specialist strain (S381). Our study revealed that the S831 strain shows a stronger affinity to Le^b than other strains due to an increased favorable intermolecular electrostatic contribution. Furthermore, we showed that the α1-2-linked fucose contributed most to the binding by forming several hydrogen bonds with key amino acids. Finally, we studied the effect of the acidic environment on the BabA-glycan complexation via constant pH MD simulations, which showed a reduction in the binding free energy in the acidic environment. Overall, our study provides a detailed understanding of the molecular mechanism of Le^b recognition by four isoforms of H. pylori that may help the development of therapeutics targeted at inhibiting H. pylori adherence to the gastric mucosa.

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.

Gaussian accelerated molecular dynamics (GaMD): principles and applications

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.