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

The conformational dynamics of Hepatitis C Virus E2 glycoprotein with the increasing number of N-glycosylation unraveled by molecular dynamics simulations

The Hepatitis C Virus (HCV), responsible for causing hepatitis and a significant contributor to liver disorders, presents a challenge for treatment due to its high genetic variability. Despite efforts, there is still no effective medication available for this virus. One of the promising targets for drug development involves targeting glycoprotein E2. However, our understanding of the dynamic behavior of E2 and its associated glycans remains limited. In this study, we investigated the dynamic characteristics of E2 with varying degrees of glycosylation using all-atom molecular dynamics simulations. We also explored glycan's interactions with the protein and among themselves. An overall increase in correlation between the vital protein regions was observed with an increase in glycan number. The protein dynamics is followed by the analysis of glycan dynamics, where the flexibility of the individual glycans was analyzed in their free and bound state, which revealed a decrease in their fluctuation in some cases. Furthermore, we generated the free energy landscape of individual N-glycan linkages in both free and bound states and observed both increases and decreases in flexibility, which can be attributed to the formation and breakage of hydrogen bonds with amino acids. Finally, we found that for a high glycosylation system, glycans interact with glycoprotein and form hydrogen bonds among themselves. Moreover, the hydrogen bond profiles of a given glycan can vary when influenced by other glycans. In summary, our study provides valuable insights into the dynamics of the core region of HCV E2 glycoprotein and its associated glycans.Communicated by Ramaswamy H. Sarma.

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.

Elucidating the conformational dynamics of histo-blood group antigens and their interactions with the rotavirus spike protein through computational lens

In the present study, we investigated the conformational dynamics of histo-blood group antigens (HBGAs) and their interactions with the VP8* domain of four rotavirus genotypes (P[4], P[6], P[19], and P[11]) utilizing all-atom molecular dynamics simulations in explicit water. Our study revealed distinct changes in the dynamic behavior of the same glycan due to linkage variations. We observed that LNFPI HBGA having a terminal β linkage shows two dominant conformations after complexation, whereas only one was obtained for LNFPI with a terminal α linkage. Interestingly, both variants displayed a single dominant structure in the free state. Similarly, LNT and LNnT show a shift in their dihedral linkage profile between their two terminal monosaccharides because of a change in the linkage from β(1-3) to β(1-4). The molecular mechanics generalized Born surface area (MM/GBSA) calculations yielded the highest binding affinity for LNFPI(β)/P[6] (-13.93 kcal/mol) due to the formation of numerous hydrogen bonds between VP8* and HBGAs. LNnT binds more strongly to P[11] (-12.88 kcal/mol) than LNT (-4.41 kcal/mol), suggesting a single change in the glycan linkage might impact its binding profile significantly. We have also identified critical amino acids and monosaccharides (Gal and GlcNAc) that contributed significantly to the protein-ligand binding through the per-residue decomposition of binding free energy. Moreover, we found that the interaction between the same glycan and different protein receptors within the same rotavirus genogroup influenced the micro-level dynamics of the glycan. Overall, our study helps a deeper understanding of the H-type HBGA and rotavirus spike protein interaction.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.

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.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

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′′-C₁₄ EGCG showed the lowest IC₅₀ values across all the tested cell lines. 4′′-C₁₄ 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′′-C₁₄ EGCG induced cell death in A431 cells with shrunken nuclei, nuclear fragmentation, membrane blebbing, and increased population of apoptotic cells where BAX upregulation and BCL_XL downregulation were observed. In addition, autophosphorylation of EGFR and its downstream signalling proteins Akt and ERK were markedly inhibited by 4′′-C₁₄ EGCG. MD simulation and the MM/PBSA analysis disclosed the binding mode of 4′′-C₁₄ EGCG in the ATP-binding site of EGFR kinase domain. Taken together, our findings demonstrate that 4′′-C₁₄ EGCG can act as a promising potent EGFR inhibitor with enhanced stability.

Among the synthesized 4′′-alkyl EGCG derivatives, 4′′-C₁₄ EGCG inhibited EGF stimulated phosphorylation of EGFR and its downstream signaling pathways, ERK and Akt. 4′′-C₁₄ EGCG showed significantly improved stability than EGCG and induced apoptosis.

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.

Mechanistic insights of ABC importer HutCD involved in heme internalization by Vibrio cholerae

Heme internalization by pathogenic bacteria inside a human host to accomplish the requirement of iron for important cellular processes is of paramount importance. Despite this, the mechanism of heme import by the ATP-binding-cassette (ABC) transporter HutCD in Vibrio cholerae remains unexplored. We have performed biochemical studies on ATPase HutD and its mutants, along with molecular modelling, docking and unbiased all-atom MD simulations on lipid-solvated models of permease-ATPase complex HutCD. The results demonstrated mechanisms of ATP binding/hydrolysis and trapped transient and global conformational changes in HutCD, necessary for heme internalization. ATPase HutD forms a dimer, independent of the permease HutC. Each HutD monomer canonically binds ATP in a 1:1 stoichiometry. MD simulations demonstrated that a rotational motion of HutC dimer occurs synchronously with the inter-dimeric D-loop interactions of HutDs. F151 of TM4–TM5 loop of HutC, packs with ATP and Y15 of HutD, initiating ‘cytoplasmic gate opening’ which mimics an ‘outward-facing’ to ‘inward-facing’ conformational switching upon ATP hydrolysis. The simulation on ‘inward-facing’ HutCD culminates to an ‘occluded’ state. The simulation on heme-docked HutCD indicated that the event of heme release occurs in ATP-free ‘inward-facing’ state. Gradual conformational changes of the TM5 helices of HutC towards the ‘occluded’ state facilitate ejection of heme.

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.

Specific interactions of peripheral membrane proteins with lipids: what can molecular simulations show us?

Peripheral membrane proteins (PMPs) can reversibly and specifically bind to biological membranes to carry out functions such as cell signalling, enzymatic activity, or membrane remodelling. Structures of these proteins and of their lipid-binding domains are typically solved in a soluble form, sometimes with a lipid or lipid headgroup at the binding site. To provide a detailed molecular view of PMP interactions with the membrane, computational methods such as molecular dynamics (MD) simulations can be applied. Here, we outline recent attempts to characterise these binding interactions, focusing on both intracellular proteins, such as phosphatidylinositol phosphate (PIP)-binding domains, and extracellular proteins such as glycolipid-binding bacterial exotoxins. We compare methods used to identify and analyse lipid-binding sites from simulation data and highlight recent work characterising the energetics of these interactions using free energy calculations. We describe how improvements in methodologies and computing power will help MD simulations to continue to contribute to this field in the future.

Conformational preferences of triantennary and tetraantennary hybrid N-glycans in aqueous solution: Insights from 20 μs long atomistic molecular dynamic simulations

In the current study, we have investigated the conformational dynamics of a triantennary (N-glycan1) and tetraantennary (N-glycan2) hybrid N-glycans found on the surface of the HIV glycoprotein using 20 μs long all-atom molecular dynamics (MD) simulations. The main objective of the present study is to elucidate the influence of adding a complex branch on the overall glycan structural dynamics. Our investigation suggests that the average RMSD value increases when a complex branch is added to N-glycan1. However, the RMSD distribution is relatively wider in the case of N-glycan1 compared to N-glycan2, which indicates that multiple complex branches restrict the conformational variability of glycans. A similar observation is obtained from the principal component analysis of both glycans. All the puckering states (⁴C₁ to ¹C₄) of each monosaccharide except mannose are sampled in our simulations, although the ⁴C₁ chair form is energetically more favorable than ¹C₄. In N-glycan1, the 1-6 linkage in the mannose branch [Man(9)-α(1-6)-Man(5)] stays in the gauche-gauche cluster, whereas it moves towards trans-gauche in N-glycan2. For both glycans, mannose branches are more flexible than the complex branches, and adding a complex branch does not influence the dynamics of the mannose branches. We have noticed that the end-to-end distance of the complex branch shortens by ∼ 10 Å in the presence of another complex branch. This suggests that in the presence of an additional complex branch, the other complex branch adopts a close folded structure. All these conformational changes involve the selective formation of inter-residue and water-mediated hydrogen-bond networks.

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.

Identification of Food Compounds as inhibitors of SARS-CoV-2 main protease using molecular docking and molecular dynamics simulations

SARS-CoV-2 has rapidly emerged as a global pandemic with high infection rate. At present, there is no drug available for this deadly disease. Recently, M^pro (Main Protease) enzyme has been identified as essential proteins for the survival of this virus. In the present work, Lipinski's rules and molecular docking have been performed to identify plausible inhibitors of M^pro using food compounds. For virtual screening, a database of food compounds was downloaded and then filtered using Lipinski's rule of five. Then, molecular docking was accomplished to identify hits using M^pro protein as the target enzyme. This led to identification of a Spermidine derivative as a hit. In the next step, Spermidine derivatives were collected from PubMed and screened for their binding with M^pro protein. In addition, molecular dynamic simulations (200 ns) were executed to get additional information. Some of the compounds are found to have strong affinity for M^pro, therefore these hits could be used to develop a therapeutic agent for SARS-CoV-2.

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.

Elucidating biophysical basis of binding of inhibitors to SARS-CoV-2 main protease by using molecular dynamics simulations and free energy calculations

The recent outbreak of novel "coronavirus disease 2019" (COVID-19) has spread rapidly worldwide, causing a global pandemic. In the present work, we have elucidated the mechanism of binding of two inhibitors, namely α-ketoamide and Z31792168, to SARS-CoV-2 main protease (Mpro or 3CLpro) by using all-atom molecular dynamics simulations and free energy calculations. We calculated the total binding free energy (ΔGbind) of both inhibitors and further decomposed ΔGbind into various forces governing the complex formation using the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) method. Our calculations reveal that α-ketoamide is more potent (ΔGbind= - 9.05 kcal/mol) compared to Z31792168 (ΔGbind= - 3.25 kcal/mol) against COVID-19 3CLpro. The increase in ΔGbind for α-ketoamide relative to Z31792168 arises due to an increase in the favorable electrostatic and van der Waals interactions between the inhibitor and 3CLpro. Further, we have identified important residues controlling the 3CLpro-ligand binding from per-residue based decomposition of the binding free energy. Finally, we have compared ΔGbind of these two inhibitors with the anti-HIV retroviral drugs, such as lopinavir and darunavir. It is observed that α-ketoamide is more potent compared to lopinavir and darunavir. In the case of lopinavir, a decrease in van der Waals interactions is responsible for the lower binding affinity compared to α-ketoamide. On the other hand, in the case of darunavir, a decrease in the favorable intermolecular electrostatic and van der Waals interactions contributes to lower affinity compared to α-ketoamide. Our study might help in designing rational anti-coronaviral drugs targeting the SARS-CoV-2 main protease.Communicated by Ramaswamy H. Sarma.

Computational modeling - an approach to the development of blood grouping reagents

Introduction: Blood group antigens are defined by an immune response that generates antibodies against a red blood cell molecule. Antibodies against these antigens can be associated with hemolytic transfusion reactions. However, difficulties can arise when developing antibodies against antigens through the use of peptide sequences alone. Three-dimensional representations (models) of the molecular structure of antigen-bearing proteins can provide valuable insights into tertiary structures and their consequent antigenicity. This can be achieved through predictive computational modeling to produce both structural and molecular dynamics models of blood group proteins.Areas covered: Authors discuss the use of molecular dynamic simulations on existing structures, as well as the use of computational modeling techniques in the development of protein models lacking preexisting data. Finally, the authors discuss specific examples of the use of computationally derived models of the MNS blood group system and its use in attempts to produce antibodies against MNS proteins.Expert opinion: Although in silico techniques have limitations, computer-based predictive models can inform the direction of research into blood group proteins. It is to be expected that as computer-based techniques grow more powerful these contributions will be even more significant.

Characterizing an allosteric inhibitor-induced inactive state in with-no-lysine kinase 1 using Gaussian accelerated molecular dynamics simulations

The binding of an allosteric inhibitor in WNK1 leads to the inactive state.

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.

An in-silico study on selected organosulfur compounds as potential drugs for SARS-CoV-2 infection via binding multiple drug targets

•

Lurasidone and Lurasidone exo inhibit SARS-CoV-2 via ‘one drug multiple targets’ strategy.

•

MD simulation revealed their strong binding affinity against multiple targets of SARS-CoV-2.

•

Both the drugs exhibit favourable pharmacokinetic properties as predicted by ADMET parameters.

•

Their unique multitargeting feature warrants further in-vitro and in-vivo experiments.

The emerging paradigm shift from ‘one molecule, one target, for one disease’ towards ‘multi-targeted small molecules’ has paved an ingenious pathway in drug discovery in recent years. We extracted this idea for the investigation of drugs for COVID-19. Perceiving the importance of organosulfur compounds, seventy-six known organosulfur compounds were screened and studied for the interaction with multiple SARS-CoV-2 target proteins by molecular dynamics simulation. Lurasidone and its derivatives displayed substantial binding affinity against five proteins (Mpro, PLpro, Spro, helicase and RdRp). The pharmacokinetics, ADMET properties and target prediction studies performed in this work further potentiates the effectiveness against SARS-CoV-2.

Plant-derived natural polyphenols as potential antiviral drugs against SARS-CoV-2 via RNA-dependent RNA polymerase (RdRp) inhibition: an in-silico analysis

The sudden outburst of Coronavirus disease (COVID-19) caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) poses a massive threat to global public health. Currently, no therapeutic drug or vaccine exists to treat COVID-19. Due to the time taking process of new drug development, drug repurposing might be the only viable solution to tackle COVID-19. RNA-dependent RNA polymerase (RdRp) catalyzes SARS-CoV-2 RNA replication and hence, is an obvious target for antiviral drug design. Interestingly, several plant-derived polyphenols effectively inhibit the RdRp of other RNA viruses. More importantly, polyphenols have been used as dietary supplementations for a long time and played beneficial roles in immune homeostasis. We were curious to study the binding of polyphenols with SARS-CoV-2 RdRp and assess their potential to treat COVID-19. Herein, we made a library of polyphenols that have shown substantial therapeutic effects against various diseases. They were successfully docked in the catalytic pocket of RdRp. The investigation reveals that EGCG, theaflavin (TF1), theaflavin-3'-O-gallate (TF2a), theaflavin-3'-gallate (TF2b), theaflavin 3,3'-digallate (TF3), hesperidin, quercetagetin, and myricetin strongly bind to the active site of RdRp. Further, a 150-ns molecular dynamic simulation revealed that EGCG, TF2a, TF2b, TF3 result in highly stable bound conformations with RdRp. The binding free energy components calculated by the MM-PBSA also confirm the stability of the complexes. We also performed a detailed analysis of ADME prediction, toxicity prediction, and target analysis for their druggability. Overall, our results suggest that EGCG, TF2a, TF2b, TF3 can inhibit RdRp and represent an effective therapy for COVID-19.Communicated by Ramaswamy H. Sarma.