Investigating Phosphorylation-Induced Conformational Changes in WNK1 Kinase by Molecular Dynamics Simulations

Macrophage WNK1 senses intracellular hypo-chlorine to regulate vulnerability to sepsis attack during hypochloremia

Sepsis is one of the leading causes of death in critical patients worldwide and its occurrence is related to the excessive activation of macrophages. Chloride loss worsens the prognosis of patients with sepsis but the underlying mechanism is currently unclear. In this study, we founded that macrophages deficient in intracellular Cl- secrete more inflammatory cytokines such as IL-1β, IL-6 and TNF-α compared with control group. The intracellular chloride level decreased in WNK1 deficiency or activity inhibited macrophages with more severe inflammatory response after LPS treatment. Remimazolam, as classic GABAa receptor agonist, alleviates excessive inflammation cascade by promoting macrophage chloride influx during sepsis progression. Collectively, this study proves that macrophage WNK1 acts as a negative regulator of inflammatory response by sensing chloride to maintain intracellular chloride balance during sepsis coupled with hypochloremia.

Protein Oxidative Modifications in Neurodegenerative Diseases: From Advances in Detection and Modelling to Their Use as Disease Biomarkers

Oxidation-reduction post-translational modifications (redox-PTMs) are chemical alterations to amino acids of proteins. Redox-PTMs participate in the regulation of protein conformation, localization and function, acting as signalling effectors that impact many essential biochemical processes in the cells. Crucially, the dysregulation of redox-PTMs of proteins has been implicated in the pathophysiology of numerous human diseases, including neurodegenerative diseases such as Alzheimer's disease and Parkinson's disease. This review aims to highlight the current gaps in knowledge in the field of redox-PTMs biology and to explore new methodological advances in proteomics and computational modelling that will pave the way for a better understanding of the role and therapeutic potential of redox-PTMs of proteins in neurodegenerative diseases. Here, we summarize the main types of redox-PTMs of proteins while providing examples of their occurrence in neurodegenerative diseases and an overview of the state-of-the-art methods used for their detection. We explore the potential of novel computational modelling approaches as essential tools to obtain insights into the precise role of redox-PTMs in regulating protein structure and function. We also discuss the complex crosstalk between various PTMs that occur in living cells. Finally, we argue that redox-PTMs of proteins could be used in the future as diagnosis and prognosis biomarkers for neurodegenerative diseases.

Design of novel DABO derivatives as HIV-1 RT inhibitors using molecular docking, molecular dynamics simulations and ADMET properties

HIV-1 reverse transcriptase is an important target for developing effective anti-HIV-1 inhibitors. Different types of small molecules have been designed based on this target, showing different levels of inhibitory activity against various types of HIV-1 strains. The relationship between structure and activity of DABO derivatives was investigated by means of 3D-QSAR molecular model, molecular docking, molecular dynamics and ADMET properties. The statistical results of molecular models show that the CoMFA and CoMSIA models have good internal stability (CoMFA: q2 = 0.623, r2 = 0.946; CoMSIA: q2 = 0.668, r2 = 0.983) and external prediction ability (CoMFA: rpred2 = 0.961; CoMSIA: rpred2 = 0.961). In addition, molecular docking has explored the mechanism of action between small molecules and receptor proteins, and the results show that hydrogen bonding between amino acid Lys101 and small molecules can improve the affinity of ligands to receptor binding. A total of 12 novel molecules were designed and their activities were predicted based on the 3D-QSAR model and molecular docking results. The results showed that the designed molecules had higher predictive activity. Subsequently, 100 ns MD simulation and binding free energy verified the stability of molecular docking results. Finally, the pharmacokinetic properties of the novel designed molecule were verified by using ADMET to predict its properties. These results can provide reference for the design and development of novel and effective HIV-1 RT inhibitors, and provide new ideas for the design of subsequent drugs.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulations reveal phosphorylation-induced conformational dynamics of the fibroblast growth factor receptor 1 kinase

The Fibroblast Growth Factor Receptor1 (FGFR1) kinase wields exquisite control on cell fate, proliferation, differentiation, and homeostasis. An imbalance of FGFR1 signaling leads to several pathogeneses of diseases ranging from multiple cancers to allergic and neurodegenerative disorders. In this study, we investigated the phosphorylation-induced conformational dynamics of FGFR1 in apo and ATP-bound states via all-atom molecular dynamics simulations. All simulations were performed for 2 × 2 µs. We have also investigated the energetics of the binding of ATP to FGFR1 using the molecular mechanics Poisson-Boltzmann scheme. Our study reveals that the FGFR1 kinase can reach a fully active configuration through phosphorylation and ATP binding. A 3-10 helix formation in the activation loop signifies its rearrangement leading to stability upon ATP binding. The interaction of phosphorylated tyrosine (pTyr654) with positively charged residues forms strong salt-bridge interactions, driving the compactness of the structure. The dynamic cross-correlation map reveals phosphorylation enhances correlated motions and reduces anti-correlated motions between different domains. We believe that the mechanistic understanding of large-conformational changes upon the activation of the FGFR1 kinase will aid the development of novel targeted therapeutics.Communicated by Ramaswamy H. Sarma.

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.

Computer-aided Affinity Enhancement of a Cross-reactive Antibody against Dengue Virus Envelope Domain III

The dengue virus (DENV), composed of four distinct but serologically related Flaviviruses, causes the most important emerging viral disease, with nearly 400 million infections yearly. Currently, there are no approved therapies. Although DENV infection induces lifelong immunity against the same serotype, the antibodies raised contribute to severe disease in heterotypic infections. Therefore, understanding the mechanism of DENV neutralization by antibodies is crucial in the design of vaccines against all serotypes. This study reports a comparative structural and energetic analysis of the monoclonal antibody (mAb) 4E11 in complex with its target domain III of the envelope protein for all four DENV serotypes. We use extensive replica molecular dynamics simulations in conjunction with the binding free energy calculations. Further single point and double mutations were designed through computational site-directed mutagenesis and observed that the re-engineered antibody exhibits high affinity to binding and broadly neutralizing activity against serotypes. Our results showed improved binding affinity by the gain of enthalpy, which could be attributed to the stabilization of salt-bridge and hydrogen bond interactions at the antigen-antibody interface. The findings provide valuable results in understanding the structural dynamics and energetic contributions that will be helpful to the design of high-affinity antibodies against dengue infections.

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.

The Integral Role of Chloride & With-No-Lysine Kinases in Cell Volume Regulation & Hypertension

Chloride anions are the most abundant in humans. For many years, it has been believed that chloride is simply a counterion of all other cations, ensuring the electroneutrality of the extracellular space. Recent data suggests that chloride anions possess a broad spectrum of important activities that regulate vital cellular functions. It is now evident that, apart from its contribution to the electroneutrality of the extracellular space, it acts as an osmole and contributes to extracellular and intracellular volume regulation. Its anionic charge also contributes to the generation of cell membrane potential. The most interesting action of chloride anions is their ability to regulate the activity of with-no-lysine kinases, which in turn regulate the activity of sodium chloride and potassium chloride cotransporters and govern the reabsorption of salt and excretion of potassium by nephron epithelia. Chloride anions seem to play a crucial role in cell functions, such as cell volume regulation, sodium reabsorption in the distal nephron, potassium balance, and sodium sensitivity, which lead to hypertension. All of these functions are accomplished on a molecular level via complicated metabolic pathways, many of which remain poorly defined. We attempted to elucidate some of these pathways in light of recent advances in our knowledge, obtained mainly from experimental studies.

Structural Insights into Phosphorylation-Mediated Polymerase Function Loss for DNA Polymerase β Bound to Gapped DNA

DNA polymerase β is a member of the X-family of DNA polymerases, playing a critical role in the base excision repair (BER) pathway in mammalian cells by implementing the nucleotide gap-filling step. In vitro phosphorylation of DNA polymerase β with PKC on S44 causes loss in the enzyme's DNA polymerase activity but not single-strand DNA binding. Although these studies have shown that single-stranded DNA binding is not affected by phosphorylation, the structural basis behind the mechanism underlying phosphorylation-induced activity loss remains poorly understood. Previous modeling studies suggested phosphorylation of S44 was sufficient to induce structural changes that impact the enzyme's polymerase function. However, the S44 phosphorylated-enzyme/DNA complex has not been modeled so far. To address this knowledge gap, we conducted atomistic molecular dynamics simulations of pol β complexed with gapped DNA. Our simulations, which used explicit solvent and lasted for microseconds, revealed that phosphorylation at the S44 site, in the presence of Mg ions, induced significant conformational changes in the enzyme. Specifically, these changes led to the transformation of the enzyme from a closed to an open structure. Additionally, our simulations identified phosphorylation-induced allosteric coupling between the inter-domain region, suggesting the existence of a putative allosteric site. Taken together, our results provide a mechanistic understanding of the conformational transition observed due to phosphorylation in DNA polymerase β interactions with gapped DNA. Our simulations shed light on the mechanisms of phosphorylation-induced activity loss in DNA polymerase β and reveal potential targets for the development of novel therapeutics aimed at mitigating the effects of this post-translational modification.

Identification of novel nematode succinate dehydrogenase inhibitors: Virtual screening based on ligand-pocket interactions

To discover new nematicidal succinate dehydrogenase (SDH) inhibitors with novel structures, we conducted a virtual screening of the ChemBridge library with 1.7 million compounds based on ligand-pocket interactions. The homology model of Caenorhabditis elegans SDH was established, along with a pharmacophore model based on ligand-pocket interactions. After the pharmacophore-based and docking-based screening, 19 compounds were selected for the subsequent enzymatic assays. The results showed that compound 1 (ID: 7607321) exhibited inhibitory activity against SDH with a determined IC₅₀ value of 19.6 μM. Structural modifications and nematicidal activity studies were then carried out, which provided further evidence that compound 1 exhibited excellent nematicidal activity. Molecular dynamics simulations were then conducted to investigate the underlying molecular basis for the potency of these inhibitors against SDH. This work provides a reliable strategy and useful information for the future design of nematode SDH inhibitors.

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.

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.

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.

Functional Loop Dynamics and Characterization of the Inactive State of the NS2B-NS3 Dengue Protease due to Allosteric Inhibitor Binding

Dengue virus, a flavivirus that causes dengue shock syndrome and dengue hemorrhagic fever, is currently prevalent worldwide. A two-component protease (NS2B-NS3) is essential for maturation, representing an important target for designing anti-flavivirus drugs. Previously, consideration has been centered on developing active-site inhibitors of NS2B-NS3pro. However, the flat and charged nature of its active site renders difficulties in developing inhibitors, suggesting an alternative strategy for identifying allosteric inhibitors. The allosterically sensitive site of the dengue protease is located near Ala125, between the 120s loop and 150s loop. Using atomistic molecular dynamics simulations, we have explored the protease's conformational dynamics upon binding of an allosteric inhibitor. Furthermore, characterization of the inherent flexible loops (71-75s loop, 120s loop, and 150s loop) is carried out for allosteric-inhibitor-bound wild-type and mutant A125C variants and a comparison is performed with its unbound state to extract the structural changes describing the inactive state of the protease. Our study reveals that compared to the unliganded system, the inhibitor-bound system shows large structural changes in the 120s loop and 150s loop in contrast to the rigid 71-75s loop. The unliganded system shows a closed-state pocket in contrast to the open state for the wild-type complex that locks the protease into the open and inactive-state conformations. However, the mutant complex fluctuates between open and closed states. Also, we tried to see how mutation and binding of an allosteric inhibitor perturb the connectivity in a protein structure network (PSN) at contact levels. Altogether, our study reveals the mechanism of conformational rearrangements of loops at the molecular level, locking the protein in an inactive conformation, which may be useful for developing allosteric inhibitors.

Phosphorylation Induced Conformational Transitions in DNA Polymerase β

DNA polymerase β (pol β) is a member of the X- family of DNA polymerases that catalyze the distributive addition of nucleoside triphosphates during base excision DNA repair. Previous studies showed that the enzyme was phosphorylated in vitro with PKC at two serines (44 and 55), causing loss of DNA polymerase activity but not DNA binding. In this work, we have investigated the phosphorylation-induced conformational changes in DNA polymerase β in the presence of Mg ions. We report a comprehensive atomic resolution study of wild type and phosphorylated DNA polymerase using molecular dynamics (MD) simulations. The results are examined via novel methods of internal dynamics and energetics analysis to reveal the underlying mechanism of conformational transitions observed in DNA pol β. The results show drastic conformational changes in the structure of DNA polymerase β due to S44 phosphorylation. Phosphorylation-induced conformational changes transform the enzyme from a closed to an open structure. The dynamic cross-correlation shows that phosphorylation enhances the correlated motions between the different domains. Centrality network analysis reveals that the S44 phosphorylation causes structural rearrangements and modulates the information pathway between the Lyase domain and base pair binding domain. Further analysis of our simulations reveals that a critical hydrogen bond (between S44 and E335) disruption and the formation of three additional salt bridges are potential drivers of these conformational changes. In addition, we found that two of these additional salt bridges form in the presence of Mg ions on the active sites of the enzyme. These results agree with our previous study of DNA pol β S44 phosphorylation without Mg ions which predicted the deactivation of DNA pol β. However, the phase space of structural transitions induced by S44 phosphorylation is much richer in the presence of Mg ions.

Phosphorylation-Induced Conformational Dynamics and Inhibition of Janus Kinase 1 by Suppressors of Cytokine Signaling 1

The dysfunction of the JAK/STAT (Janus kinase/signal transducers and activators of transcription) pathway results in several pathophysiological conditions, including autoimmune disorders. The negative feedback regulators of the JAK/STAT signaling pathway, suppressors of cytokine signaling (SOCS), act as a natural inhibitor of JAK and inhibit aberrant activity. SOCS1 is the most potent member of the SOCS family, whose kinase inhibitory region targets the substrate-binding groove of JAK with high affinity and blocks the phosphorylation of JAK kinases. Overall, we performed an aggregate of 13 μs molecular dynamics simulations on the activation loop's three different phosphorylation (double and single) states. Results from our simulations show that the single Tyr¹⁰³⁴ phosphorylation could stabilize the JAK1/SOCS1 complex as well as the flexible activation segment. The phosphate-binding loop (P-loop) shows conformational variability at dual and single phosphorylated states. Principal component analysis and protein structure network (PSN) analysis reveal that the different phosphorylation states and SOCS1 binding induce intermediate inactive conformations of JAK1, which could be a better target for future JAK1 selective drug design. PSN analysis suggests that the com-pY1034 system is stabilized due to higher values of network hubs than the other two complex systems. Moreover, the binding free energy calculations suggest that pTyr¹⁰³⁴ states show a higher affinity toward SOCS1 than the dual and pTyr¹⁰³⁵ states. We believe that the mechanistic understanding of JAK1/SOCS1 complexation will aid future studies related to peptide inhibitors based on SOCS1.

Structural Rearrangement of the Serotonin Transporter Intracellular Gate Induced by Thr276 Phosphorylation

The reuptake of the neurotransmitter serotonin from the synaptic cleft by the serotonin transporter, SERT, is essential for proper neurological signaling. Biochemical studies have shown that Thr276 of transmembrane helix 5 is a site of PKG-mediated SERT phosphorylation, which has been proposed to shift the SERT conformational equilibria to promote inward-facing states, thus enhancing 5-HT transport. Recent structural and simulation studies have provided insights into the conformation transitions during substrate transport but have not shed light on SERT regulation via post-translational modifications. Using molecular dynamics simulations and Markov state models, we investigate how Thr276 phosphorylation impacts the SERT mechanism and its role in enhancing transporter stability and function. Our simulations show that Thr276 phosphorylation alters the hydrogen-bonding network involving residues on transmembrane helix 5. This in turn decreases the free energy barriers for SERT to transition to the inward-facing state, thus facilitating 5-HT import. The results provide atomistic insights into in vivo SERT regulation and can be extended to other pharmacologically important transporters in the solute carrier family.

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.

Discovery of Novel Tubulin Polymerization Inhibitors by Utilizing 3D-QSAR, Molecular Docking and Molecular Dynamics Simulation

Tubulin is a potential therapeutic target for cancer. Compounds inhibit the polymerization of tubulin or promote the polymerization of tubulin to interfere with the mitotic process of cells, resulting in...

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.

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.

Identification of Potential Inhibitors against Epstein-Barr Virus Nuclear Antigen 1 (EBNA1): An Insight from Docking and Molecular Dynamic Simulations

Epstein-Barr virus (EBV), a known tumorigenic virus, is associated with various neuropathies, including multiple sclerosis (MS). However, there is no anti-EBV FDA-approved drug available in the market. Our study targeted EBV protein EBV nuclear antigen 1 (EBNA1), crucial in virus replication and expressed in all the stages of viral latencies. This dimeric protein binds to an 18 bp palindromic DNA sequence and initiates the process of viral replication. We chose phytochemicals and FDA-approved MS drugs based on literature survey followed by their evaluation efficacies as anti-EBNA1 molecules. Molecular docking revealed FDA drugs ozanimod, siponimod, teriflunomide, and phytochemicals; emodin; protoapigenone; and EGCG bound to EBNA1 with high affinities. ADMET and Lipinski's rule analysis of the phytochemicals predicted favorable druggability. We supported our assessments of pocket druggability with molecular dynamics simulations and binding affinity predictions by the molecular mechanics generalized Born surface area (MM/GBSA) method. Our results establish a stable binding for siponimod and ozanimod with EBNA1 mainly via van der Waals interactions. We identified hot spot residues like I481', K477', L582', and K586' in the binding of ligands. In particular, K477' at the amino terminal of EBNA1 is known to establish interaction with two bases at the major groove of the DNA. Siponimod bound to EBNA1 engaging K477', thus plausibly making it unavailable for DNA interaction. Computational alanine scanning further supported the significant roles of K477', I481', and K586' in the binding of ligands with EBNA1. Conclusively, the compounds showed promising results to be used against EBNA1.

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.

Phosphorylation-dependent subfunctionalization of the calcium-dependent protein kinase CPK28

Calcium (Ca2+)-dependent protein kinases (CDPKs or CPKs) are a unique family of Ca2+ sensor/kinase-effector proteins with diverse functions in plants. In Arabidopsis thaliana, CPK28 contributes to immune homeostasis by promoting degradation of the key immune signaling receptor-like cytoplasmic kinase BOTRYTIS-INDUCED KINASE 1 (BIK1) and additionally functions in vegetative-to-reproductive stage transition. How CPK28 controls these seemingly disparate pathways is unknown. Here, we identify a single phosphorylation site in the kinase domain of CPK28 (Ser318) that is differentially required for its function in immune homeostasis and stem elongation. We show that CPK28 undergoes intermolecular autophosphorylation on Ser318 and can additionally be transphosphorylated on this residue by BIK1. Analysis of several other phosphorylation sites demonstrates that Ser318 phosphorylation is uniquely required to prime CPK28 for Ca2+ activation at physiological concentrations of Ca2+, possibly through stabilization of the Ca2+-bound active state as indicated by intrinsic fluorescence experiments. Together, our data indicate that phosphorylation of Ser318 is required for the activation of CPK28 at low intracellular [Ca2+] to prevent initiation of an immune response in the absence of infection. By comparison, phosphorylation of Ser318 is not required for stem elongation, indicating pathway-specific requirements for phosphorylation-based Ca2+-sensitivity priming. We additionally provide evidence for a conserved function for Ser318 phosphorylation in related group IV CDPKs, which holds promise for biotechnological applications by generating CDPK alleles that enhance resistance to microbial pathogens without consequences to yield.

Mining of Ebola virus genome for the construction of multi-epitope vaccine to combat its infection

Ebola virus is the primary causative agent of viral hemorrhagic fever that is an epidemic disease and responsible for the massive premature deaths in humans. Despite knowing the molecular mechanism of its pathogenesis, to date, no commercial or FDA approved multiepitope vaccine is available against Ebola infection. The current study focuses on designing a multi-epitope subunit vaccine for Ebola using a novel immunoinformatic approach. The best predicted antigenic epitopes of Cytotoxic-T cell (CTL), Helper-T cells (HTL), and B-cell epitopes (BCL) joined by various linkers were selected for the multi-epitope vaccine designing. For the enhanced immune response, two adjuvants were also added to the construct. Further analysis showed the vaccine to be immunogenic and non-allergenic, forming a stable and energetically favorable structure. The stability of the unbound vaccine construct and vaccine/TLR4 was elucidated via atomistic molecular dynamics simulations. The binding free energy analysis (ΔG_Bind = -194.2 ± 0.5 kcal/mol) via the molecular mechanics Poisson-Boltzmann docking scheme revealed a strong association and thus can initiate the maximal immune response. Next, for the optimal expression of the vaccine construct, its gene construct was cloned in the pET28a + vector system. In summary, the Ebola viral proteome was screened to identify the most potential HTLs, CTLs, and BCL epitopes. Along with various linkers and adjuvants, a multi-epitope vaccine is constructed that showed a high binding affinity with the immune receptor, TLR4. Thus, the current study provides a highly immunogenic multi-epitope subunit vaccine construct that may induce humoral and cellular immune responses against the Ebola infection.Communicated by Ramaswamy H. Sarma.

Dimeric phosphorylation of glyoxalase I alters its symmetry and substrate binding mechanism: simulation studies

Glyoxalase I (GLO1) is a dimeric esterase of the glyoxalase system. Phosphorylation of the residue T106 has been found to inhibit GLO1 activity, and contribute to the onset of oxidative stress and cellular damage. This research uses multiple molecular dynamics simulations and automated docking of both GLO1 and dimerically phosphorylated GLO1 (p2-GLO1) to predict the initial structural differences induced by phosphorylation, and their interaction with the intermediate substrate Hemimercaptal. This research indicates that immediately following phosphorylation, GLO1 exhibits reduced sphericity, partly caused by outward splaying of the loop region surrounding T106. Phosphorylation induces enhanced concerted motions in the loop composed of residues immediately surrounding T106, which are correlated with motions at the active site pocket at the distant, opposite end of the dimer. These T106 region loop motions result in the distortion of the shape of the active site, and potentially alter its accessibility. Phosphorylation alters the manner in which GLO1 interacts with Hemimercaptal. For GLO1, Hemimercaptal is predicted to bind to T106, which we propose constitutes a novel, highly accessible 'capture site' responsible for initial contact with the substrate. In contrast, for p2-GLO1, Hemimercaptal is unable to bind favourably to (phosphorylated) position T106, suggesting that this proposed transient 'capture site' is abolished upon phosphorylation of GLO1. Hence, a novel physiological role is here proposed for the known essential GLO1 residue T106. These results may further contribute to understanding the inhibition mechanism of GLO1 upon phosphorylation.Communicated by Ramaswamy H. Sarma.

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.

Activation Loop Dynamics Are Coupled to Core Motions in Extracellular Signal-Regulated Kinase-2

The activation loop segment in protein kinases is a common site for regulatory phosphorylation. In extracellular signal-regulated kinase 2 (ERK2), dual phosphorylation and conformational rearrangement of the activation loop accompany enzyme activation. X-ray structures show the active conformation to be stabilized by multiple ion pair interactions between phosphorylated threonine and tyrosine residues in the loop and six arginine residues in the kinase core. Despite the extensive salt bridge network, nuclear magnetic resonance Carr-Purcell-Meiboom-Gill relaxation dispersion experiments show that the phosphorylated activation loop is conformationally mobile on a microsecond to millisecond time scale. The dynamics of the loop match those of previously reported global exchange within the kinase core region and surrounding the catalytic site that have been found to facilitate productive nucleotide binding. Mutations in the core region that alter these global motions also alter the dynamics of the activation loop. Conversely, mutations in the activation loop perturb the global exchange within the kinase core. Together, these findings provide evidence for coupling between motions in the activation loop and those surrounding the catalytic site in the active state of the kinase. Thus, the activation loop segment in dual-phosphorylated ERK2 is not held statically in the active X-ray conformation but instead undergoes exchange between conformers separated by a small energetic barrier, serving as part of a dynamic allosteric network controlling nucleotide binding and catalytic function.

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.

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

Acquired immune deficiency syndrome (AIDS) is caused by the human immunodeficiency virus (HIV), type 1 and 2. Further, the diversity in HIV-1 has given rise to many serotypes and recombinant strains. The currently used protease inhibitors have been developed for subtype B, although non-B subtype strains account for ∼ 90% of the global HIV infections. Subtype D is spreading rapidly and infecting a large population in North Africa and the Middle East. In the current study, molecular dynamics simulations in conjunction with the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) scheme was used to investigate the potency of four drugs, namely atazanavir (ATV), darunavir (DRV), lopinavir (LPV) and tipranavir (TPV) against the subtype D variant. Our calculations predicted that the potency of the inhibitors decreased in the order TPV > ATV > DRV > LPV. TPV was found to be the most potent against subtype D due to an increase in van der Waals and electrostatic interactions and reduction in the desolvation energy compared to other inhibitors. This result is further supported by the hydrogen bond interactions between inhibitors and protease. Furthermore, our analyses suggested that the binding of TPV induced a more closed conformation of the flap compared to apo or other complexes. It was observed that TPV/PR^D has a lower cavity volume relative to the other three complexes leading to a tighter binding. The open conformation of the flap was observed for LPV/PR^D. We expect that this study might be useful for designing more potent inhibitors against HIV-1 subtype D. Communicated by Ramaswamy H. Sarma.