Structural dynamics of Casein Kinase I (CKI) from malarial parasite Plasmodium falciparum (Isolate 3D7): Insights from theoretical modelling and molecular simulations

Engineered β-glycosidase from Hyperthermophilic Sulfolobus solfataricus with Improved Rd-hydrolyzing Activity for Ginsenoside Compound K Production

Hyperthermophilic Sulfolobus solfataricus β-glycosidase (SS-βGly), with higher stability and activity than mesophilic enzymes, has potential for industrial ginsenosides biotransformation. However, its relatively low ginsenoside Rd-hydrolyzing activity limits the production of pharmaceutically active minor ginsenoside compound K (CK). In this study, first, we used molecular docking to predict the key enzyme residues that may hypothetically interact with ginsenoside Rd. Then, based on sequence alignment and alanine scanning mutagenesis approach, key variant sites were identified that might improve the enzyme catalytic efficiency. The enzyme catalytic efficiency (kcat/Km) and substrate affinity (Km) of the N264D variant enzyme for ginsenoside Rd increased by 60% and decreased by 17.9% compared with WT enzyme, respectively, which may be due to a decrease in the binding free energy (∆G) between the variant enzyme and substrate Rd. In addition, Markov state models (MSM) analysis during the whole 1000-ns MD simulations indicated that altering N264 to D made the variant enzyme achieve a more stable SS-βGly conformational state than the wild-type (WT) enzyme and corresponding Rd complex. Under identical conditions, the relative activities and the CK conversion rates of the N264D enzyme were 1.7 and 1.9 folds higher than those of the WT enzyme. This study identified an excellent hyperthermophilic β-glycosidase candidate for industrial biotransformation of ginsenosides.

Post-Translational Modifications of Proteins of Malaria Parasites during the Life Cycle

Post-translational modifications (PTMs) are essential for regulating protein functions, influencing various fundamental processes in eukaryotes. These include, but are not limited to, cell signaling, protein trafficking, the epigenetic control of gene expression, and control of the cell cycle, as well as cell proliferation, differentiation, and interactions between cells. In this review, we discuss protein PTMs that play a key role in the malaria parasite biology and its pathogenesis. Phosphorylation, acetylation, methylation, lipidation and lipoxidation, glycosylation, ubiquitination and sumoylation, nitrosylation and glutathionylation, all of which occur in malarial parasites, are reviewed. We provide information regarding the biological significance of these modifications along all phases of the complex life cycle of Plasmodium spp. Importantly, not only the parasite, but also the host and vector protein PTMs are often crucial for parasite growth and development. In addition to metabolic regulations, protein PTMs can result in epitopes that are able to elicit both innate and adaptive immune responses of the host or vector. We discuss some existing and prospective results from antimalarial drug discovery trials that target various PTM-related processes in the parasite or host.

Classical molecular dynamics simulation identifies catechingallate as a promising antiviral polyphenol against MPOX palmitoylated surface protein

Cumulative global prevalence of the emergent monkeypox (MPX) infection in the non-endemic countries has been professed as a global public health predicament. Lack of effective MPX-specific treatments sets the baseline for designing the current study. This research work uncovers the effective use of known antiviral polyphenols against MPX viral infection, and recognises their mode of interaction with the target F13 protein, that plays crucial role in formation of enveloped virions. Herein, we have employed state-of-the-art machine learning based AlphaFold2 to predict the three-dimensional structure of F13 followed by molecular docking and all-atoms molecular dynamics (MD) simulations to investigate the differential mode of F13-polyphenol interactions. Our extensive computational approach identifies six potent polyphenols Rutin, Epicatechingallate, Catechingallate, Quercitrin, Isoquecitrin and Hyperoside exhibiting higher binding affinity towards F13, buried inside a positively charged binding groove. Intermolecular contact analysis of the docked and MD simulated complexes divulges three important residues Asp134, Ser137 and Ser321 that are observed to be involved in ligand binding through hydrogen bonds. Our findings suggest that ligand binding induces minor conformational changes in F13 to affect the conformation of the binding site. Concomitantly, essential dynamics of the six-MD simulated complexes reveals Catechin gallate, a known antiviral agent as a promising polyphenol targeting F13 protein, dominated with a dense network of hydrophobic contacts. However, assessment of biological activities of these polyphenols need to be confirmed through in vitro and in vivo assays, which may pave the way for development of new novel antiviral drugs.

Structural insights of Labeo catla (catla) myxovirus resistance protein,GTP binding recognition and constitutive expression induced with Poly I:C

The Myxovirus resistance (Mx) proteins are critical effectors belonging to the super-family of guanidine triphosphatase, often stimulated by type I interferon (IFN) and mediates antiviral responses to restrict the replication of numerous viral genes in fishes. In teleosts, Mx proteins display diverse and complicated antiviral activity in different species. The present investigation seeks to characterize the Mx gene from Labeo catla upon induction by double-stranded (ds) RNA, polyinosinic-polycytidylic acid, (poly I: C). Molecular modeling and all-atoms molecular dynamics (MD) simulations were employed to understand the architecture of the GTPase domain and its plausible mode of GTP recognition in Mx protein. The full-length L. catla Mx (LcMx) gene sequence (1821 bp nucleotides) encodes an open reading frame of 606 amino acids. Domain search indicated conserved tripartite domain architecture of LcMx and forms a major cluster with the Mx from other teleosts. The positively charged Arginine and polar Glutamine residues from helix 3 and 4 of stalk region LcMx aid in homo-oligomerization. MD simulation portrayed the role of conserved critical residues aid in GTP recognition by the GTPase domain which perfectly corroborates with experimental findings and prior MD studies. After injection of poly I:C, the temporal mRNA profile showed that LcMx expression was significantly elevated in the spleen, brain, kidney, liver, muscle, heart, intestine, and gill tissues. Collectively, these results suggest that the elevated expression of the major innate immune defense gene Mx was able to inhibit the poly I: C mediated virulence in fish.Communicated by Ramaswamy H. Sarma.

All-atoms molecular dynamics study to screen potent efflux pump inhibitors against KpnE protein of Klebsiella pneumoniae

The Small Multidrug Resistance efflux pump protein KpnE, plays a pivotal role in multi-drug resistance in Klebsiella pneumoniae. Despite well-documented study of its close homolog, EmrE, from Escherichia coli, the mechanism of drug binding to KpnE remains obscure due to the absence of a high-resolution experimental structure. Herein, we exclusively elucidate its structure-function mechanism and report some of the potent inhibitors through drug repurposing. We used molecular dynamics simulation to develop a dimeric structure of KpnE and explore its dynamics in lipid-mimetic bilayers. Our study identified both semi-open and open conformations of KpnE, highlighting its importance in transport process. Electrostatic surface potential map suggests a considerable degree of similarity between KpnE and EmrE at the binding cleft, mostly occupied by negatively charged residues. We identify key amino acids Glu14, Trp63 and Tyr44, indispensable for ligand recognition. Molecular docking and binding free energy calculations recognizes potential inhibitors like acarbose, rutin and labetalol. Further validations are needed to confirm the therapeutic role of these compounds. Altogether, our membrane dynamics study uncovers the crucial charged patches, lipid-binding sites and flexible loop that could potentiate substrate recognition, transport mechanism and pave the way for development of novel inhibitors against K. pneumoniae.Communicated by Ramaswamy H. Sarma.

Computational study on Schiff base derived salicylaldehyde and furfuraldehyde derivatives as potent anti-tubercular agents: prospect to dihydropteroate synthase inhibitors

Nowadays, bacterial multidrug resistance has become a commonplace problem in clinics due to several intrinsic factors mediated through resistance to antibacterials obtained via bacterial consortia and extrinsic factors, such as non-uniform antibacterial policy and migration of resistant bacteria through human and other routes. The development of newer, effective anti-mycobacterial candidate(s) is coveted by clinics. Hybrid molecules would be comparatively more emulating against invasive bacterial strains; nevertheless, newer antibiotics are continually added. Herein, designing and developments of two series of Schiff-based salicylaldehyde S1-S7 and furfuraldehyde F1-F7 molecules individually bearing sulfonamide group are described; and those were synthesized and their structures by spectral characterization were confirmed. Concomitantly, molecule dynamic simulations of all atoms had been performed to fathom the mechanism of the action with these leading complexes. These data imply that the synthesized Schiff-based salicylaldehyde hybrids would be promising anti-tubercular compounds, which further need potent pharmacological evaluations.Communicated by Ramaswamy H. Sarma.

Spectroscopic insights with molecular docking and molecular dynamic simulation studies of anticancer drug 5-Fluorouracil targeting human pyruvate kinase m2

This study was conducted to test the efficacy of 5-fluorouracil (5-FU) as an anticancer drug against the human pyruvate kinase isozyme M2 (PKM2) using spectroscopic, molecular docking and molecular dynamic simulation studies. PKM2 fluorescence quenching studies in the presence of 5-FU performed at three different temperatures indicates dynamic quenching processes with single-set of binding (n ≈ 1) profile. The biomolecular quenching constants (kq) and the effective binding constants (Kb) obtained are shown to increase with temperature. The calculated enthalpy (ΔH) and entropy changes (ΔS) are estimated to be -118.06 kJ/mol and 146.14 kJ/mol/K respectively, which suggest the possible mode of interaction as electrostatic and hydrogen bonding. Further, these values were used to estimate the free energy changes (ΔG) and that increases with temperature. The negative ΔG values clearly indicates spontaneous binding process that stabilizes the complex formed between 5-FU and PKM2. Far-UV CD spectra of PKM2 in the presence of 5-FU shows decrease in α-helix contents which point towards the destabilization of secondary structure that weakens the biological activity of PKM2. The intrinsic fluorescence study and circular dichroism (CD) spectra showed minor conformational changes of PKM2 in the presence of 5-FU. Additionally, the results obtained from molecular docking and all-atom molecular dynamic simulation study supports the insight of the spectroscopic binding studies, and strengthens the dynamic stability of the complex between 5-FU and PKM2 through H-bonding. This study establishes a paradigm of 5-FU-PKM2 complexation and the efficacy of 5-FU that compromises the biological activity of the targeted PKM2.Communicated by Ramaswamy H. Sarma.

Lipid and cholesterols modulate the dynamics of SARS-CoV-2 viral ion channel ORF3a and its pathogenic variants

SARS-CoV-2 accessory protein, ORF3a is a putative ion channel which immensely contributes to viral pathogenicity by modulating host immune responses and virus-host interactions. Relatively high expression of ORF3a in diseased individuals and implication with inflammasome activation, apoptosis and autophagy inhibition, ratifies as an effective target for developing vaccines and therapeutics. Herein, we present the elusive dynamics of ORF3a-dimeric state using all-atoms molecular dynamics (MD) simulations at μ-seconds scale in a heterogeneous lipid-mimetic system in multiple replicates. Additionally, we also explore the effect of non-synonymous pathogenic mutations on ORF3a ion channel activity and viral pathogenicity in different SARS-CoV-2 variants using various structure-based protein stability (ΔΔG) tools and computational saturation mutagenesis. Our study ascertains the role of phosphatidylcholines and cholesterol in modulating the structure of ORF3a, which perturbs the size and flexibility of the polar cavity that allows permeation of large cations. Discrete trend in ion channel pore radius and area per lipid arises the premise that presence of lipids might also affect the overall conformation of ORF3a. MD structural-ensembles, in some replicates rationalize the crucial role of TM2 in maintaining the native structure of ORF3a. We also infer that loss of structural stability primarily grounds for pathogenicity in more than half of the pathogenic variants of ORF3a. Overall, the effect of mutation on alteration of ion permeability of ORF3a, proposed in this study brings mechanistic insights into variant consequences on viral membrane proteins of SARS-CoV-2, which can be utilized for the development of novel therapeutics to treat COVID-19 and other coronavirus diseases.

Effect of Fragment 1 on the Binding of Epigallocatechin Gallate to the PD-L1 Dimer Explored by Molecular Dynamics

Blocking the interaction between programmed cell death-1 (PD-1) and programmed cell death-ligand 1 (PD-L1) by directly targeting the PD-L1 dimer has emerged as a hot topic in the field of cancer immunotherapy. Epigallocatechin gallate (EGCG), a natural product, has been demonstrated binding to the PD-L1 dimer in our previous study, but has a weaker binding capacity, moreover, EGCG is located at the end of the binding pocket of the PD-L1 dimer. The inhibitor fragment 1 (FRA) lies at the other end. So, we proposed that the introduction of FRA might be able to improve the binding ability. To illuminate this issue, molecular dynamics (MD) simulation was performed in the present study. Binding free energy calculations show that the binding affinity is significantly increased by 17 kcal/mol upon the introduction of FRA. It may be due to the energy contributions of emerging key residues ATyr56, AMet115, BTyr123, AIle54 and the enhanced contributions of initial key residues ATyr123 and BVal68. Binding mode and non-bonded interaction results indicate that FRA_EGCG (EGCG in combination with FRA) binds to the C-, F- and G-sheet of the PD-L1 dimer. Importantly, the introduction of FRA mainly strengthened the nonpolar interactions. The free energy landscape and secondary structure results further show that FRA_EGCG can interact with the PD-L1 dimer more stably. These data demonstrated here provide the theoretical basis for screening two or more natural products with additive inhibitory effect on this pathway and therefore exerting more effective anticancer immunity.

Identification of PPA1 inhibitor candidates for potential repurposing in cancer medicine

Inorganic pyrophosphatase 1 (PPA1) is pivotal to cellular metabolism as it facilitates the hydrolysis of PPi-a by-product of various metabolic processes that influence cell growth and differentiation. Overexpression of PPA1 enzyme has been linked to diminished patient survival and was shown to influence tumor cell dynamics, thereby positioning it as a potential therapy target for a variety of cancers including colorectal cancer, diffuse large B-cell lymphoma, and lung adenocarcinoma. Despite this therapeutic promise, there are no known inhibitors of PPA1 as of today. In this study, we searched for potential PPA1 inhibitors using a molecular docking screen of 30 470 compounds with a history of clinical trials and/or US Food and Drug Administration approval. We specifically targeted the active pocket that coincides with the established catalytic domain. Our screen identified promising hits, which we further subjected to ADMET (absorption, distribution, metabolism, excretion, and toxicity) filtering. Subsequent molecular dynamics (MD) analyses were conducted on devazepide, quinotolast, and tarazepide-the three substances that successfully navigated all filters. MD analyses reinforced the stability of the protein-ligand complexes and confirmed ligand binding, as substantiated by our root mean square deviation, radius of gyration and secondary structures of proteins analyses. Furthermore, Molecular Mechanics Poisson-Boltzmann Surface Area calculations post-MD identified devazepide and quinotolast as showing higher binding affinities; being supported by principal component analysis, free energy landscape, and dynamic cross-correlation matrix results. Overall, our study reveals devazepide and quinotolast as potential candidates for PPA1 inhibition which could be considered for repurposing studies that need further experimental validation. These results not only reveal a potential for clinical repurposing for PPA1 inhibition but they also offer valuable insights into the development of future compounds for targeting the crucial PPA1 enzyme.

Exploiting the potential of natural polyphenols as antivirals against monkeypox envelope protein F13 using machine learning and all-atoms MD simulations

The re-emergence of monkeypox (MPX), in the era of COVID-19 pandemic is a new global menace. Regardless of its leniency, there are chances of MPX expediting severe health deterioration. The role of envelope protein, F13 as a critical component for production of extracellular viral particles makes it a crucial drug target. Polyphenols, exhibiting antiviral properties have been acclaimed as an effective alternative to the traditional treatment methods for management of viral diseases. To facilitate the development of potent MPX specific therapeutics, herein, we have employed state-of-the-art machine learning techniques to predict a highly accurate 3-dimensional structure of F13 as well as identify binding hotspots on the protein surface. Additionally, we have effectuated high-throughput virtual screening methodology on 57 potent natural polyphenols having antiviral activities followed by all-atoms molecular dynamics (MD) simulations, to substantiate the mode of interaction of F13 protein and polyphenol complexes. The structure-based virtual screening based on Glide SP, XP and MM/GBSA scores enables the selection of six potent polyphenols having higher binding affinity towards F13. Non-bonded contact analysis, of pre- and post- MD complexes propound the critical role of Glu143, Asp134, Asn345, Ser321 and Tyr320 residues in polyphenol recognition, which is well supported by per-residue decomposition analysis. Close-observation of the structural ensembles from MD suggests that the binding groove of F13 is mostly hydrophobic in nature. Taken together, this structure-based analysis from our study provides a lead on Myricetin, and Demethoxycurcumin, which may act as potent inhibitors of F13. In conclusion, our study provides new insights into the molecular recognition and dynamics of F13-polyphenol bound states, offering new promises for development of antivirals to combat monkeypox. However, further in vitro and in vivo experiments are necessary to validate these results.

Structural insights into SARS-CoV-2 main protease conformational plasticity

The spread of different severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants underscores the need for insights into the structural properties of its structural and non-structural proteins. The highly conserved homo-dimeric chymotrypsin-like protease (3CL M^PRO ), belonging to the class of cysteine hydrolases, plays an indispensable role in processing viral polyproteins that are involved in viral replication and transcription. Studies have successfully demonstrated the role of M^PRO as an attractive drug target for designing antiviral treatments because of its importance in the viral life cycle. Herein, we report the structural dynamics of six experimentally solved structures of M^PRO (i.e., 6LU7, 6M03, 6WQF, 6Y2E, 6Y84, and 7BUY including both ligand-free and ligand-bound states) at different resolutions. We have employed a structure-based balanced forcefield, CHARMM36m through state-of-the-art all-atoms molecular dynamics simulations at µ-seconds scale at room temperature (303K) and pH 7.0 to explore their structure-function relationship. The helical domain-III responsible for dimerization mostly contributes to the altered conformational states and destabilization of M^PRO . A keen observation of the high degree of flexibility in the P5 binding pocket adjoining domain II-III highlights the reason for observation of conformational heterogeneity among the structural ensembles of M^PRO . We also observe a differential dynamics of the catalytic pocket residues His41, Cys145, and Asp187, which may lead to catalytic impairment of the monomeric proteases. Among the highly populated conformational states of the six systems, 6LU7 and 7M03 forms the most stable and compact M^PRO conformation with intact catalytic site and structural integrity. Altogether, our findings from this extensive study provides a benchmark to identify physiologically relevant structures of such promising drug targets for structure-based drug design and discovery of potent drug-like compounds having clinical potential.

Identification of core therapeutic targets for Monkeypox virus and repurposing potential of drugs against them: An in silico approach

Monkeypox virus (mpox virus) outbreak has rapidly spread to 82 non-endemic countries. Although it primarily causes skin lesions, secondary complications and high mortality (1-10%) in vulnerable populations have made it an emerging threat. Since there is no specific vaccine/antiviral, it is desirable to repurpose existing drugs against mpox virus. With little knowledge about the lifecycle of mpox virus, identifying potential inhibitors is a challenge. Nevertheless, the available genomes of mpox virus in public databases represent a goldmine of untapped possibilities to identify druggable targets for the structure-based identification of inhibitors. Leveraging this resource, we combined genomics and subtractive proteomics to identify highly druggable core proteins of mpox virus. This was followed by virtual screening to identify inhibitors with affinities for multiple targets. 125 publicly available genomes of mpox virus were mined to identify 69 highly conserved proteins. These proteins were then curated manually. These curated proteins were funnelled through a subtractive proteomics pipeline to identify 4 highly druggable, non-host homologous targets namely; A20R, I7L, Top1B and VETFS. High-throughput virtual screening of 5893 highly curated approved/investigational drugs led to the identification of common as well as unique potential inhibitors with high binding affinities. The common inhibitors, i.e., batefenterol, burixafor and eluxadoline were further validated by molecular dynamics simulation to identify their best potential binding modes. The affinity of these inhibitors suggests their repurposing potential. This work can encourage further experimental validation for possible therapeutic management of mpox.

Molecular dynamics simulation, synthesis and topoisomerase inhibitory actions of vanillin derivatives: a systematic computational structural integument

A series of 4-hydroxy-3-methoxy benzaldehyde (vanillin) derivatives (3a-3r) was designed for the principle of Schiff base condensation with several individual sulfanilamide analogues. The inhibitory potencies of the designed compounds were evaluated through molecular docking simulation studies against the targets, breast cancer-topo isomerase-IIα and estrogen receptor-α; and the top scoring poses with higher binding energy were selected to assess the mode of binding and stability of each complex through molecular dynamics simulations. Compounds that remained stable in the active sites of the both target receptors through a number of strong H-bonds and hydrophobic contacts were selected. Based on the computational results, these selected compounds, 3b, 3e and 3f were synthesized and were followed up for structural elucidation attempts, by FT/ATR, ¹H NMR and ¹³C NMR. From the experimental in vitro studies on 3b, 3e and 3f, the following remarkable activities against breast cancer cell line were done; IC₅₀ values of 3b, 3e and 3f were noted, 6.7, 4.3 and 11 ng/mL, respectively. These newly synthesized compounds may be used as novel inhibitors of nuclear receptors with potential therapeutic applications in control of cancer.Communicated by Ramaswamy H. Sarma.

Novel Insights into Understanding the Molecular Dialogues between Bipolaroxin and the Gα and Gβ Subunits of the Wheat Heterotrimeric G-Protein during Host–Pathogen Interaction

Spot blotch disease of wheat, caused by the fungus Bipolaris sorokiniana (Sacc.) Shoem., produces several toxins which interact with the plants and thereby increase the blightening of the wheat leaves, and Bipolaroxin is one of them. There is an urgent need to decipher the molecular interaction between wheat and the toxin Bipolaroxin for in-depth understanding of host–pathogen interactions. In the present study, we have developed the three-dimensional structure of G-protein alpha subunit from Triticum aestivum. Molecular docking studies were performed using the active site of the modeled G-protein alpha and cryo-EM structure of beta subunit from T. aestivum and ‘Bipolaroxin’. The study of protein–ligand interactions revealed that six H-bonds are mainly formed by Glu29, Ser30, Lys32, and Ala177 of G-alpha with Bipolaroxin. In the beta subunit, the residues of the core beta strand domain participate in the ligand interaction where Lys256, Phe306, and Leu352 formed seven H-bonds with the ligand Bipolaroxin. All-atoms molecular dynamics (MD) simulation studies were conducted for G-alpha and -beta subunit and Bipolaroxin complexes to explore the stability, conformational flexibility, and dynamic behavior of the complex system. In planta studies clearly indicated that application of Bipolaroxin significantly impacted the physio-biochemical pathways in wheat and led to the blightening of leaves in susceptible cultivars as compared to resistant ones. Further, it interacted with the Gα and Gβ subunits of G-protein, phenylpropanoid, and MAPK pathways, which is clearly supported by the qPCR results. This study gives deeper insights into understanding the molecular dialogues between Bipolaroxin and the Gα and Gβ subunits of the wheat heterotrimeric G-protein during host–pathogen interaction.

Structural insights into the RNA interaction with Yam bean Mosaic virus (coat protein) from Pachyrhizus erosus using bioinformatics approach

Plants are constantly threatened by a virus infection, i.e., Potyviruses, the second largest genus of plant viruses which results in several million-dollar losses in various essential crops globally. Yam bean (Pachyrhizus erosus) is considered to be one of the essential tuberous legume crops holding a great potential source of starch. Yam Bean Mosaic Virus (YBMV) of Potyvirus group belonging to the family potyviridae affects Yam bean and several angiosperms both in the tropical and sub-tropical regions causing large economical losses in crops. In this study, we attempted to understand the sequence-structure relationship and mode of RNA binding mechanism in YBMV CP using in silico integrative modeling and all-atoms molecular dynamics (MD) simulations. The assembly of coat protein (CP) subunits from YBMV and the plausible mode of RNA binding were compared with the experimental structure of CP from Watermelon mosaic virus potyvirus (5ODV). The transmembrane helix region is present in the YBMV CP sequence ranging from 76 to 91 amino acids. Like the close structural-homolog, 24 CPs monomeric sub-units formed YBMV a conserved fold. Our computational study showed that ARG124, ARG155, and TYR151 orient towards the inner side of the virion, while, THR122, GLN125, SER92, ASP94 reside towards the outer side of the virion. Despite sharing very low sequence similarity with CPs from other plant viruses, the strongly conserved residues Ser, Arg, and Asp within the RNA binding pocket of YBMV CP indicate the presence of a highly conserved RNA binding site in CPs from different families. Using several bioinformatics tools and comprehensive analysis from MD simulation, our study has provided novel insights into the RNA binding mechanism in YBMV CP. Thus, we anticipate that our findings from this study will be useful for the development of new therapeutic agents against the pathogen, paving the way for researchers to better control this destructive plant virus.

Optimization of the Urea Removal in a Wearable Dialysis Device Using Nitrogen-Doped and Phosphorus-Doped Graphene

Dialysis has been recognized as an essential treatment for end-stage renal disease (ESRD). This therapy, however, suffers from several limitations leading to numerous complications in the patients. As dialysis cannot completely substitute healthy kidney functions, the health condition of an ESRD patient is ultimately affected. Wearable artificial kidney (WAK) can resolve the restrictions of blood purification by the dialysis method. However, absorbing large amounts of urea produced in the body is one of the main challenges of these WAK and overcoming this is necessary to improve both functionality and footprint of the device. This study investigates the adsorption capabilities of N- and P-doped graphene nanosorbents for the first time by using molecular dynamic simulation. Urea removal on carbon nanosheets was simulated with different percentages of phosphorus and nitrogen dopants along with the pristine graphene. Specifically, the effects of interaction energy, adsorption percentage, gyration radius, hydrogen bonding, and other molecular dynamic analyses on urea removal were also investigated. The results from this study match well with the existing research, demonstrating the accuracy of the model. The results further suggest that graphene nanosheets doped by 10% nitrogen are likely the most effective in removing urea given that it is associated with the maximum radial distribution function (RDF), the maximum reduction in gyration radius, a high number of hydrogen bonds, and the most negative adsorption energy. This molecular study offers attractive suggestions for the novel adsorbents of artificial kidney devices and paves the way for the development of novel and enhanced urea adsorbents.

Naltrexone a potential therapeutic candidate for COVID-19

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is the cause of Coronavirus Disease (COVID-19) that has resulted in a global pandemic. At the time of writing, approximately 16.06 million cases have been reported worldwide. Like other coronaviruses, SARS-CoV-2 relies on the surface Spike glycoprotein to access the host cells, mainly through the interaction of its Receptor Binding Domain (RBD) with the host receptor Angiotensin-Converting Enzyme2 (ACE2). SARS-CoV-2 infection induces a profound downstream pro-inflammatory cytokine storm. This release of the pro-inflammatory cytokines is underpinning lung tissue damage, respiratory failure, and eventually multiple organ failure in COVID-19 patients. The phosphorylation status of ERK1/2 is positively correlated with virus load and ERK1/2 inhibition suppressed viral replication and viral infectivity. Therefore, molecular entities able to interfere with binding of the SARS-CoV-2 Spike protein to ACE2, or damping hyperinflammatory cytokines storm, blocking ERK1/2 phosphorylation have a great potential to inhibit viral entry along with viral infectivity. Herein, we report that the FDA-approved non-peptide opioid antagonist drug, naltrexone suppresses high fat/LPS induced pro-inflammatory cytokine release both from macrophage cells and Adipose Tissue Macrophage. Moreover, Low Dose Naltrexone (LDN) also showed its activity as an ERK1/2 inhibitor. Notably, virtual docking and simulation data also suggest LDN may disrupt the interaction of ACE2 with RBD. LDN may be considered as a target as the treatment and (or) adjuvant therapy for coronavirus infection. Clinical toxicity measurements may not be required for LDN since naltrexone was previously tested and is an approved drug by the FDA.Communicated by Ramaswamy H. Sarma.

Atomistic insight into 2D COFs as antiviral agents against SARS-CoV-2

The recent pandemic of COVID-19 has raised global health concerns. Preventing severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) activity in the body is a very promising method to overcome the COVID-19 pandemic. One of the prevention methods is constraining the binding process among the human cell receptor-ACE2 and coronavirus spike protein. In the research done, the effect of deformation of the spike protein structure, due to the covalent organic frameworks (COFs), in reducing the interactions of ACE2 and the spike protein by the computational method was investigated. In this regard, atomic analysis of the interactions of ACE2 and the spike protein is provided using a molecular dynamics simulation. First, we investigated the interactions of the three different COFs, including COF-78, DAAQ-TFP, and COF-OEt, with the spike protein by analyzing the bond energies, as well as structural changes of the spike protein. Then, intermolecular interactions of the deformed spike protein along with ACE2 were assessed to clarify the protein's fusion after the deformation. As indicated by the results, although all introduced COFs deformed the spike protein in an effective way, COF-78 showed the best performance in the prevention of spike protein-ACE2 interactions by changing the molecular structure of the protein. Indeed, the interaction analysis of the deformed spike protein by COF-78 with the ACE2 showed that their interactions had the lowest absolute value of energy, along with the least amount of hydrogen bonds, in which the compaction of the protein was lower compared to the other deformed proteins. Moreover, having a high contact area with an aqueous media as well as severe fluctuations during the simulation time confirmed the positive performance of COF-78. In the current study, we aimed to introduce novel materials and COVID-19 prevention methodology that can be used in face masks and for surface disinfection.

Molecular Mechanism of Food-Derived Polyphenols on PD-L1 Dimerization: A Molecular Dynamics Simulation Study

In cancer immunotherapy, an emerging approach is to block the interactions of programmed cell death-1 (PD-1) and programmed cell death-ligand 1 (PD-L1) using small-molecule inhibitors. The food-derived polyphenols curcumin (CC), resveratrol (RSV) and epigallocatechin gallate (EGCG) have anticancer immunologic functions, which, recently, have been proposed to act via the downregulation of PD-L1 expression. However, it remains unclear whether they can directly target PD-L1 dimerization and, thus, interrupt the PD-1/PD-L1 pathway. To elucidate the molecular mechanism of such compounds on PD-L1 dimerization, molecular docking and nanosecond molecular dynamics simulations were performed. Binding free energy calculations show that the affinities of CC, RSV and EGCG to the PD-L1 dimer follow a trend of CC > RSV > EGCG. Hence, CC is the most effective inhibitor of the PD-1/PD-L1 pathway. Analysis on contact numbers, nonbonded interactions and residue energy decomposition indicate that such compounds mainly interact with the C-, F- and G-sheet fragments of the PD-L1 dimer, which are involved in interactions with PD-1. More importantly, nonpolar interactions between these compounds and the key residues Ile54, Tyr56, Met115, Ala121 and Tyr123 play a dominant role in binding. Free energy landscape and secondary structure analyses further demonstrate that such compounds can stably interact with the binding domain of the PD-L1 dimer. The results provide evidence that CC, RSV and EGCG can inhibit PD-1/PD-L1 interactions by directly targeting PD-L1 dimerization. This provides a novel approach to discovering food-derived small-molecule inhibitors of the PD-1/PD-L1 pathway with potential applications in cancer immunotherapy.

Molecular insight into optimizing the N- and P-doped fullerenes for urea removal in wearable artificial kidneys

Urea is the result of the breakdown of proteins in the liver, the excess of which circulates in the blood and is adsorbed by the kidneys. However, in the case of kidney diseases, some products, specifically urea, cannot be removed from the blood by the kidneys and causes serious health problems. The end-stage renal disease (ESRD) patients are not able to purify their blood, which endangers their life. ESRD patients require dialysis, a costly and difficult method of urea removal from the blood. Wearable artificial kidneys (WAKs) are consequently designed to remove the waste from blood. Regarding the great amount of daily urea production in the body, WAKs should contain strong and selective urea adsorbents. Fullerenes-which possess fascinating chemical properties-have been considered herein to develop novel urea removal adsorbents. Molecular dynamics (MD) has enabled researchers to study the interaction of different materials and can pave the way toward facilitating the development of wearable devices. In this study, urea adsorption by N-doped fullerenes and P-doped fullerenes were assessed through MD simulations. The urea adsorption was simulated by five samples of fullerenes, with phosphorous and different nitrogen dopant contents. For comparing the urea adsorption capacity in the performed simulations, detailed characteristics-including the energy analysis, radius of gyration, radial distribution function (RDF), root-mean-square fluctuation (RMSD), and H-bond analyses were investigated. It had been determined that the fullerene containing 8% nitrogen-with the highest reduction in the radius of gyration, the maximum RDF, a high adsorption energy, and a high number of hydrogen bonds-adsorbs urea more efficiently.

Kanamycin-Mediated Conformational Dynamics of Escherichia coli Outer Membrane Protein TolC

TolC is a member of the outer membrane efflux proteins (OEPs) family and acts as an exit duct to export proteins, antibiotics, and substrate molecules across the Escherichia coli cell membrane. Export of these molecules is evidenced to be brought about through the reversible interactions and binding of substrate-specific drug molecules or antibiotics with TolC and by being open for transport, which afterward leads to cross-resistance. Hence, the binding of kanamycin with TolC was monitored through molecular docking (MD), the structural fluctuations and conformational changes to the atomic level. The results were further supported from the steady-state fluorescence binding and isothermal titration calorimetry (ITC) studies. Binding of kanamycin with TolC resulted in a concentration dependent fluorescence intensity quenching with 7 nm blue shift. ITC binding data maintains a single binding site endothermic energetic curve with binding parameters indicating an entropy driven binding process. The confirmational changes resulting from this binding were monitored by a circular dichroism (CD) study, and the results showed insignificant changes in the α-helix and β-sheets secondary structure contents, but the tertiary structure shows inclusive changes in the presence of kanamycin. The experimental data substaintially correlates the RMSD, R _g, and RMSF results. The resulting conformational changes of the TolC-kanamycin complexation was stabilized through H-bonding and other interactions.

iTRAQ-based comparative proteomic analysis in different developmental stages of Echinococcus granulosus

Cystic echinococcosis, caused by infection with the larval stage of the cestode Echinococcus granulosus, is a chronic zoonosis. The lifecycle of the E. granulosus parasite includes three consecutive stages that require specific gene regulation or protein expression to survive environmental shifts between definitive hosts and intermediate hosts. The aim of the present study is to screen and analyze the stage differential antigens to be considered for vaccine development against E. granulosus. By using the iTRAQ (isobaric tags for relative and absolute quantification) method, the differentially expressed proteins were selected from the three consecutive developmental stages of E. granulosus: oncosphere, adult tapeworms, and protoscolex. Through a bioinformatics analysis including Clusters of Orthologous Groups (COG), Gene Ontology (GO), and pathway metabolic annotation, we identified some proteins of interest from each stage. The results showed that a large number of differentially expressed proteins (375: oncosphere vs. adult, 346: oncosphere vs. protoscolex, and 391: adult vs. protoscolex) were identified from the three main lifecycle stages. Analysis of the differential protein pathways showed that these differential proteins are mainly enriched in metabolic pathways, Huntington’s diseases, Alzheimer’s diseases, and ribosome metabolic pathways. Interestingly, among these differential proteins, expression levels of paramyosin, HSP60, HSP70, HSP90, cathepsin L1, cathepsin D, casein kinase, and calmodulin were significantly higher in the oncosphere than in the adult or protoscolex (p < 0.05). We hope our findings will help to identify potential targets for diagnosis or for therapeutic and prophylactic intervention.

Sequential Dihedral Angles (SDAs): A Method for Evaluating the 3D Structure of Proteins

One of the most important steps in modeling three-dimensional (3D) structures of proteins is the evaluation of the constructed models. The present study suggests that the correctness of a structure may be tested by using the characteristics of sequential dihedral angles (SDAs) between adjacent alpha-carbons (C_α) in the main chains of proteins. From our studies on protein structures in the protein data bank (PDB), the SDAs between the C_α in the main chains are limited in their values. In addition, the sum of the absolute values of the three sequential dihedral angles (SDAs) can never be 0 degree. Moreover, 48 degrees is the lowest value existing for the sum of the absolute values of three sequential dihedral angles (SDAs). Thus, the SDAs between the alpha-carbons along the main chains of proteins may be a useful parameter for evaluating anomalies in protein structures.

Insights into structure and dynamics of extracellular domain of Toll-like receptor 5 in Cirrhinus mrigala (mrigala): A molecular dynamics simulation approach

The toll-like receptor 5 (TLR5) is the most conserved important pattern recognition receptors (PRRs) often stimulated by bacterial flagellins and plays a major role in the first-line defense against invading pathogenic bacteria and in immune homeostasis. Experimental crystallographic studies have shown that the extracellular domain (ECD) of TLR5 recognizes flagellin of bacteria and functions as a homodimer in model organism zebrafish. However, no structural information is available on TLR5 functionality in the major carp Cirrhinus mrigala (mrigala) and its interaction with bacterial flagellins. Therefore, the present study was undertaken to unravel the structural basis of TLR5-flagellin recognition in mrigala using structural homodimeric TLR5-flagellin complex of zebrafish as reference. Integrative structural modeling and molecular dynamics simulations were employed to explore the structural and mechanistic details of TLR5 recognition. Results from structural snapshots of MD simulation revealed that TLR5 consistently formed close interactions with the three helices of the D1 domain in flagellin on its lateral side mediated by several conserved amino acids. Results from the intermolecular contact analysis perfectly substantiate with the findings of per residue-free energy decomposition analysis. The differential recognition mediated by flagellin to TLR5 in mrigala involves charged residues at the interface of binding as compared to the zebrafish complex. Overall our results shows TLR5 of mrigala involved in innate immunity specifically recognized a conserved site on flagellin which advocates the scientific community to explore host-specific differences in receptor activation.

Graphene-Based Nanoparticles as Potential Treatment Options for Parkinson's Disease: A Molecular Dynamics Study

INTRODUCTION

The study of abnormal aggregation of proteins in different tissues of the body has recently earned great attention from researchers in various fields of science. Concerning neurological diseases, for instance, the accumulation of amyloid fibrils can contribute to Parkinson's disease, a progressively severe neurodegenerative disorder. The most prominent features of this disease are the degeneration of neurons in the substantia nigra and accumulation of α-synuclein aggregates, especially in the brainstem, spinal cord, and cortical areas. Dopamine replacement therapies and other medications have reduced motor impairment and had positive consequences on patients' quality of life. However, if these medications are stopped, symptoms of the disease will recur even more severely. Therefore, the improvement of therapies targeting more basic mechanisms like prevention of amyloid formation seems to be critical. It has been shown that the interactions between monolayers like graphene and amyloids could prevent their fibrillation.

METHODS

For the first time, the impact of four types of last-generation graphene-based nanostructures on the prevention of α-synuclein amyloid fibrillation was investigated in this study by using molecular dynamics simulation tools.

RESULTS

Although all monolayers were shown to prevent amyloid fibrillation, nitrogen-doped graphene (N-Graphene) caused the most instability in the secondary structure of α-synuclein amyloids. Moreover, among the four monolayers, N-Graphene was shown to present the highest absolute value of interaction energy, the lowest contact level of amyloid particles, the highest number of hydrogen bonds between water and amyloid molecules, the highest instability caused in α-synuclein particles, and the most significant decrease in the compactness of α-synuclein protein.

DISCUSSION

Ultimately, it was concluded that N-Graphene could be the most effective monolayer to disrupt amyloid fibrillation, and consequently, prevent the progression of Parkinson's disease.

Effect of mutation on structure, function and dynamics of receptor binding domain of human SARS-CoV-2 with host cell receptor ACE2: a molecular dynamics simulations study

Recent studies have pointed the role of angiotensin-converting enzyme-II (ACE2) in mediating the entry of SARS-CoV-2 to the host cell by binding to the receptor-binding domain (RBD) of viral spike protein, and successive priming by cellular proteases initiates the infection. SARS-CoV replication rate and disease severity is controlled by the binding affinity of RBD with ACE2. To understand, how mutations in the conserved residues of RBD affect the molecular interaction with ACE2, we generated five alanine mutants i.e. Y449A, N487A, Y489A, N501A and Y505A in the receptor binding motif (RBM) of the ACE2-RBD SARS-CoV-2 complex (PDB: 6M0J). Computational site directed mutagenesis induced dynamics in wild-type and mutant complexes were extensively studied through all-atoms molecular dynamics (MD) simulations of 150 ns. In silico mutational analysis revealed loss of important intermolecular hydrogen bonds and other non-bonded contacts, critical for molecular recognition of SARS-CoV-2 RBD to ACE2, which is well supported by saturation mutagenesis study of binding interface residues. MD simulations results showed that RBM motif is flexible, where mutant residues are relatively more mobile than corresponding wild-type residues. Global motion analysis through principal component studies revealed that RBD exhibits protuberant in-ward motion towards the human ACE2 binding interface which may be crucial for molecular interaction. Conclusively, the present finding are in congruence with previous experimental reports and provides detailed information on the structural basis of receptor binding by human SARS-CoV-2, which will crucial for the development of novel inhibitors or drugs to combat against SARS-CoV-2.

Communicated by Ramaswamy H. Sarma

Analysis of polyphenolic metabolites from in vitro gastrointestinal digested soft fruit extracts identify malvidin-3-glucoside as an inhibitor of PTP1B

Protein-tyrosine phosphatase 1B (PTP1B, EC 3.1.3.48) is an important regulator of insulin signalling. Herein, we employed experimental and computational biology techniques to investigate the inhibitory properties of phenolics, identified from four in vitro gastrointestinal digested (IVGD) soft fruits, on PTP1B. Analysis by LC-MS/MS identified specific phenolics that inhibited PTP1B in vitro. Enzyme kinetics identified the mode of inhibition, while dynamics, stability and binding mechanisms of PTP1B-ligand complex were investigated through molecular modelling, docking, molecular dynamics (MD) simulations, and MM/PBSA binding free energy estimation. IVGD extracts and specific phenolics identified from the four soft fruits inhibited PTP1B (P < 0.0001) activity. Among the phenolics tested, the greatest inhibition was shown by malvidin-3-glucoside (P < 0.0001) and gallic acid (P < 0.0001). Malvidin-3-glucoside (Ki = 3.8 µg/mL) was a competitive inhibitor and gallic acid (Ki = 33.3 µg/mL) a non-competitive inhibitor of PTP1B. Malvidin-3-glucoside exhibited better binding energy than gallic acid and the synthetic inhibitor Dephostatin (-7.38 > -6.37 > -5.62 kcal/mol) respectively. Principal component analysis demonstrated malvidin-3-glucoside PTP1B-complex occupies more conformational space where critical WPD-loop displayed a higher degree of motion. MM/PBSA binding free energy for malvidin-3-glucoside to PTP1B was found to be higher than other complexes mediated by Van der Waals energy rather than electrostatic interaction for the other two inhibitors (-80.32 ± 1.25 > -40.64 ± 1.43 > -21.63 ± 1.73 kcal/mol) respectively. Altogether, we have established novel insights into the specific binding of dietary phenolics and have identified malvidin-3-glucoside as an PTP1B inhibitor, which may be further industrially developed for the treatment of type-2 diabetes.

Computer-aided synthesis of dapsone-phytochemical conjugates against dapsone-resistant Mycobacterium leprae

Leprosy continues to be the belligerent public health hazard for the causation of high disability and eventual morbidity cases with stable prevalence rates, even with treatment by the on-going multidrug therapy (MDT). Today, dapsone (DDS) resistance has led to fear of leprosy in more unfortunate people of certain developing countries. Herein, DDS was chemically conjugated with five phytochemicals independently as dapsone-phytochemical conjugates (DPCs) based on azo-coupling reaction. Possible biological activities were verified with computational chemistry and quantum mechanics by molecular dynamics simulation program before chemical synthesis and spectral characterizations viz., proton-HNMR, FTIR, UV and LC-MS. The in vivo antileprosy activity was monitored using the 'mouse-foot-pad propagation method', with WHO recommended concentration 0.01% mg/kg each DPC for 12 weeks, and the host-toxicity testing of the active DPC4 was seen in cultured-human-lymphocytes in vitro. One-log bacilli cells in DDS-resistant infected mice footpads decreased by the DPC4, and no bacilli were found in the DDS-sensitive mice hind pads. Additionally, the in vitro host toxicity study also confirmed that the DCP4 up to 5,000 mg/L level was safety for oral administration, since a minor number of dead cells were found in red color under a fluorescent microscope. Several advanced bioinformatics tools could help locate the potential chemical entity, thereby reducing the time and resources required for in vitro and in vitro tests. DPC4 could be used in place of DDS in MDT, evidenced from in vivo antileprosy activity and in vitro host toxicity study.

Structural Characterization of Open Reading Frame-Encoded Functional Genes from Tilapia Lake Virus (TiLV)

In recent years, large-scale mortalities are observed in tilapia due to infection with a novel orthomyxo-like virus named, tilapia lake virus (TiLV) which is marked to be a severe threat to universal tilapia industry. Currently, there are knowledge gaps relating to the antiviral peptide as well as there are no affordable vaccines or drugs available against TiLV yet. To understand the spreading of infection of TiLV in different organs of Oreochromis niloticus, RT-PCR analysis has been carried out. The gene segments of TiLV were retrieved from the NCBI database for computational biology analysis. The 14 functional genes were predicted from the 10 gene segments of TiLV. Phylogenetic analysis was employed to find out a better understanding for the evolution of tilapia lake virus genes. Out of 14 proteins, only six proteins show transmembrane helix region. Moreover, molecular modeling and molecular dynamics simulations of the predicted proteins revealed structural stability of the protein stabilized after 10-ns simulation. Overall, our study provided a basic bioinformatics on functional proteome of TiLV. Further, this study could be useful for development of novel peptide-based therapeutics to control TiLV infection.

Molecular dynamics and computational study of Mannich-based coumarin derivatives: potent tyrosine kinase inhibitor

The manifestation of bacterial UTI (Urinary Tract Infection) has been predominantly endemic, globally; eventually, the development of new UTI antibacterial agent(s) remains the call of the day. Herein, two series of Mannich-based 4-hydroxy coumarin derivatives, 7a-m and 8a-m were designed by suitable heterocyclic amines condensed with aldehydes. The synthesised molecules were interpreted by ¹H-NMR and ¹³C-NMR spectral analyses with in vitro antibacterial studies. The compound, 4-hydroxy-3-((4-hydroxy-3-methoxyphenyl)(morpholino)methyl)-2H-chromen-2-one 8l was the significant derivative against pathogenic bacteria Staphylococcus aureus and Escherichia coli with MIC values 12.50 and 25 µM, respectively. Computational assessments with the Lipinski's rule of five, ADMET properties and molecular docking studies revealed that analogues, 7f, 7l, 8d, 8j and 8k could be potent druggable molecules with significant binding affinity towards bacterial tyrosine kinase, as target. To understand the mode of binding and intrinsic stabilities of potent receptor-ligand complexes, each system was subjected molecular dynamics simulations for 100 ns. Inter-molecular contact analysis and intrinsic hydrogen-bond stability portrayed the analogues 8l form a number of non-bonded contacts with the receptor tyrosine kinase being mostly dominated by electrostatic and hydrophobic contacts. The results from the present structure-based designing approach might be a valuable tool towards identification of a new antibacterial drug candidate(s) against UTI in near future.Abbreviations4-HC4-hydroxy coumarinADMETabsorption, distribution, metabolism, excretion, toxicityDCE1,2-dichloroethaneDCMdichloromethaneDMSOdimethyl sulfoxideMDmolecular dynamicsMICminimum inhibitory concentrationNMRnuclear magnetic resonanceNPTisothermal isobaric ensemblesNVTcanonical ensemblesPCAprincipal component analysisPDBprotein data bankRgradius gyrationRMSFroot-mean-square fluctuationsRSMDroot-mean-square deviationUTIurinary Tract InfectionCommunicated by Ramaswamy H. Sarma.

Virtual screening to identify novel potential inhibitors for Glutamine synthetase of Mycobacterium tuberculosis

Glutamine synthetase (GS) of Mycobacterium tuberculosis (Mtb) is an essential enzyme which is involved in nitrogen metabolism and cell wall synthesis. It is involved in the inhibition of phagosome-lysosome fusion by preventing acidification. Targeting GS can be helpful to control the infection of Mtb. In order to identify potential inhibitors, we screened chemical libraries (56,400 compounds of ChEMBL anti-mycobacterial, 1596 FDA approved drugs, 419 Natural product and 916 phytochemical) against this target using the virtual screening approach. Screening by molecular docking identified ten top-ranked compounds as GSMtb inhibitors and they were compared with known inhibitors (as control). Since GS enzyme (GSHs) is also present in human. We have compared the protein sequence of GS from Mtb and human using the P-BLAST in NCBI. We found ∼27% identity in between these two sequences, so we also compared the binding affinity of inhibitor between Mtb and human. Finally, we identified top two compounds namely CHEMBL387509, CHEMBL226198 from ChEMBL anti-mycobacterial dataset, and Eriocitrin and Malvidin from phytochemical dataset which showed lees binding affinity towards GSHs whereas Pamidronate, and Phentermine from FDA approved drugs and (-)-Quinic Acid, Hexopyranuronic acid, Quebrachit, and Castanospermine from natural product showed protein-ligand interaction with Mtb protein while no interaction with GSHs. The top two docked complexes were subjected to molecular dynamic simulation to understand the stability of the molecule. Further, we calculated the binding free energy of the docked complex and analyzed hydrogen bond, salt bridge, pie stacking, and hydrophobic interaction in the docking region. These ligands exhibited very good binding affinity GSMtb enzymes. Therefore, these ligands are novel and drug-likeness compounds, and they may be potential inhibitors of M tuberculosis.Communicated by Ramaswamy H. Sarma.

Molecular characterization, constitutive expression and GTP binding mechanism of Cirrhinus mrigala (Hamilton, 1822) Myxovirus resistance (Mx) protein

Myxovirus resistance (Mx) proteins represents the subclass of the dynamin superfamily of large Guanosine triphosphates (GTPases), play esential role in intracellular vesicle trafficking, endocytosis, organelle homeostasis and mitochondria distribution. These proteins are key players of the vertebrate immune system, induced by type-I and type-III interferons (IFN) of infected host and inhibit viral replication by sequestering its nucleoprotein. In the present study, we report the sequencing and characterization of Cirrhinus mrigala Mx protein (CmMx) for the first time and observed its constitutive expression in different tissues for a period of fourteen days. The synthetic peptide, LSGVALPRGTGI, was dissolved in PBS and injected into a rabbit and the antibody raised against CmMx was used to study the level of its expression. The full length of the CmMx cDNA is 2244 bp with a molecular mass of 70.9 kDa and a predicted isoelectric point of 8.25. The 627 amino acids polypeptide formed of three main functional domains: N-terminal GTPase domain (GD), a middle domain (MD) and GTPase effector domain (GED) with carboxy terminal leucine zipper motif. The 3D models of CmMx protein was modeled based on available close structural homologs and further validated through molecular dynamics (MD) simulations. MD study revealed the importance of G-domain responsible for recognition of GTP, which perfectly corroborate with earlier studies. MM/PBSA binding free energy analysis displayed that van der Waals and electrostatic energy were the key driving force behind molecular recognition of GTP by CmMx protein. The results from this study will illuminate more lights into the ongoing research on myxovirus resistance protein and its role in inhibition of viral replication in other eukaryotic system as well.

Insights into the mode of flavin mononucleotide binding and catalytic mechanism of bacterial chromate reductases: A molecular dynamics simulation study

Enzymes from natural sources protect the environment via complex biological mechanisms, which aid in reductive immobilization of toxic metals including chromium. Nevertheless, progress was being made in elucidating high-resolution crystal structures of reductases and their binding with flavin mononucleotide (FMN) to understand the underlying mechanism of chromate reduction. Therefore, herein, we employed molecular dynamics (MD) simulations, principal component analysis (PCA), and binding free energy calculations to understand the dynamics behavior of these enzymes with FMN. Six representative chromate reductases in monomeric and dimeric forms were selected to study the mode, dynamics, and energetic component that drive the FMN binding process. As evidenced by MD simulation, FMN prefers to bind the cervix formed between the catalytic domain surrounded by strong conserved hydrogen bonding, electrostatic, and hydrophobic contacts. The slight movement and reorientation of FMN resulted in breakage of some crucial H-bonds and other nonbonded contacts, which were well compensated with newly formed H-bonds, electrostatic, and hydrophobic interactions. The critical residues aiding in tight anchoring of FMN within dimer were found to be strongly conserved in the bacterial system. The molecular mechanics combined with the Poisson-Boltzmann surface area binding free energy of the monomer portrayed that the van der Waals and electrostatic energy contribute significantly to the total free energy, where, the polar solvation energy opposes the binding of FMN. The proposed proximity relationships between enzyme and FMN binding site presented in this study will open up better avenues to engineer enzymes with optimized chromate reductase activity for sustainable bioremediation of heavy metals.

Molecular cloning, GTP recognition mechanism and tissue-specific expression profiling of myxovirus resistance (Mx) protein in Labeo rohita (Hamilton) after Poly I:C induction

The myxovirus resistance (Mx) proteins belong to interferon-induced dynamin GTPase and play pivotal role in the inhibition of replication of numerous viruses. These antiviral proteins are released in usual or diseased condition to prevent the viral attack and to carry regular cellular activities like endocytosis and trafficking of nucleoproteins into the nucleus. The invasion of virus up-regulates the expression of Mx transcripts and double-stranded RNA mimic like polyinosinic polycytidyilic acid (Poly I:C). To understand the tissue-specific expression profiling and mechanism of GTP recognition of Mx protein from Labeo rohita (rohu), the full-length gene was cloned, sequenced and characterized through various Bioinformatics tools for the first time. The Mx cDNA was comprised of 2297 bp, and the open reading frame of 1938 bp encodes polypeptide of 631 amino acids. The coding sequence of Mx protein possess the signature motif of dynamin superfamily, LPRG(S/K)GIVTR, the tripartite guanosine-5/triphosphate (GTP)-binding motif (GXXXSGKS/T, DXXG and T/NKXD) and the leucine zipper motifs at the C-terminal end, well conserved in all interferon-induced Mx protein in vertebrates. Western blotting confirmed the molecular weight of Mx protein to be 72 kDa. After the intraperitoneal challenge of L. rohita with a Poly I:C, up-regulation of Mx protein was observed in brain, spleen, liver, kidney, intestine, heart, muscle, and gill. Ontogeny study displayed pronounced expression of Mx protein in all stages of the developmental of Rohu after Poly I:C induction. However a persistent expression of Mx transcript was also observed in Rohu egg as well as milt without induction with Poly I:C. Higher expression of Mx gene was observed on 96 h where it was 6.4 folds higher than the control. The computational modelling of Mx protein portrayed the tripartite N-terminal G-domain that binds to GTP, the bundle-signaling element (BSE) which interconnects the G-domain to the elongated stalk domain and C-terminal helical stalk domain. In agreement with the experimental studies, a series of conserved residues viz., Gln52, Ser53, Ser54, Leu68, Pro69, Gly71, Gly73, Thr76, Asp151, Gly154, Thr220, Lys221, Val251, Cys253, Arg254, and Gly255 were computed to be indispensable for tight anchoring of GTP within binding cavity of G-domain. The binding free energy calculation study depicted that the van der Waals and electrostatic terms contributs significantly to molecular recognition of GTP. Collectively, our study provides mechanistic insights into the tissue-specific expression profiling and GTP binding mechanism of Mx protein from Labeo rohita, which is expected to drive further research on several cellular events including viral resistance and endocytosis in the near future.

A theoretical and experimental investigation of the effect of sodium dodecyl sulfate on the structural and conformational properties of bovine β-casein

A predicted three-dimensional structure of bovine β-casein was constructed using homology modeling with the aid of MODELLER and I-TASSER programs, with the validity and reliability of the models evaluated according to stereochemical qualities and small angle X-ray scattering. By comparing the results obtained from the two models using the CRYSOL program, an optimal model of the β-casein structure derived from I-TASSER was selected and used in subsequent molecular dynamics (MD) analysis. 300 ns MD simulations of β-casein in water and in the presence of different SDS concentrations at 300 K were performed. The results of the MD simulations indicated that SDS molecules played a dual role in modifying the conformation of β-casein at 300 K. Concentrations of SDS below its CMC (1 mM), at which only the monomer form of SDS was present, induced β-casein to lose its secondary structure by converting helices into random coils; however the conformation of the complex was still comparable with that of native β-casein. In the presence of 10 mM SDS (above its CMC), the helical content of β-casein was increased along with reduced random coils, and the structural rearrangement led to a more compact conformation. The latter change is likely related to the hydrophobic interactions that dominate the binding of the C-terminal region, along with the anchoring of sulfate groups of SDS on the positively charged N-terminal portion via electrostatic attraction. Hydrogen bonding supplemented the SDS-induced stabilization of β-casein. A correlated "necklace and bead" model, in which the micelles nucleate on the protein hydrophobic sites, was proposed for the structure of β-casein-SDS complexes.

Structural bioinformatics insights into ATP binding mechanism in zebrafish (Danio rerio) cyclin-dependent kinase-like 5 (zCDKL5) protein

In mammalian systems, the conserved cyclin-dependent protein kinases (CDKs) control the process of cell division and curb the transcription mechanism in response to diverse signaling events that are essential for the catalytic activity. In zebrafish, zCDKL5 portrays differential expression profiling in several tissues and presumed to play a vital role in the neuronal development. In this present study, the sequence-structure relationship and mode of ATP binding in zCDKL5 was unveiled through theoretical modeling, molecular docking, and MD simulations. Like human CDKs, the modeled zCDKL5 was found to be bipartite in nature, where, ATP binds to the central cavity of the catalytic domain through a strong network of H-bonding, electrostatic, and hydrophobic interactions. MD simulation portrayed that conserved residues, viz, Ile10, Gly11, Glu12, Val18, Val64, Glu81, Cys143, and Asp144 were indispensable for tight anchoring of ATP and contribute to the stability of the zCDKL5-ATP complex. MM/PBSA binding free energy analysis displayed that van der Waal energy (ΔG _vwd ) and Electrostatic energy (ΔG _ele ) were the major contributors towards the overall binding free energy. Thus, the comparative structural bioinformatics approach has shed new insights into the dynamics and ATP binding mechanism of zCDKL5. The results from the study will help to undertake further research on the role of phosphorylated CDKL5 in the onset of neurodevelopmental disorders caused by mutations in higher eukaryotic systems.

The potential HLA Class I-restricted epitopes derived from LeIF and TSA of Leishmania donovani evoke anti-leishmania CD8+ T lymphocyte response

To explore new protective measure against visceral leishmaniasis, reverse vaccinology approach was employed to identify key immunogenic regions which can mediate long-term immunity. In-depth computational analysis revealed nine promiscuous epitopes which can possibly be presented by 46 human leukocyte antigen, thereby broadening the worldwide population up to 94.16%. This is of reasonable significance that most of the epitopes shared 100% sequence homology with other Leishmania species and could evoke a common pattern of protective immune response. Transporter associated with antigen processing binding affinity, molecular docking approach followed by dynamics simulation and human leukocyte antigen stabilization assay suggested that the best five optimal set of epitopes bind in between α1 and α2 binding groove with sufficient affinity and stability which allows the translocation of intact epitope to the cell surface. Fascinatingly, the human leukocyte antigen stabilization assay exhibited a modest correlation with the positive immunogenicity score predicted by class I pMHC immunogenicity predictor. A support for this notion came from ELISA and FACS analysis where the epitopes as a cocktail induced CD8+ IFN-γ and Granzyme B levels significantly in treated visceral leishmaniasis subject which suggests the immunogenic ability of the selected epitopes.

Deep insights into the mode of ATP-binding mechanism in Zebrafish cyclin-dependent protein kinase-like 1 (zCDKL1): A molecular dynamics approach

In eukaryotes, the serine/threonine kinases (STKs) belonging to cyclin-dependent protein kinases (CDKs) play significant role in control of cell division and curb transcription in response to several extra and intra-cellular signals indispensable for enzymatic activity. The zebrafish cyclin-dependent protein kinase-like 1 protein (zCDKL1) shares a high degree of sequence and structural similarity with mammalian orthologs and express in brain, ovary, testis, and low levels in other tissues. Regardless of its importance in the developmental process, the structure, function and mode of ATP recognition have not been investigated yet due to lack of experimental data. Henceforth, to gain atomistic insights in to the structural dynamics and mode of ATP binding, a series of computational techniques involving theoretical modeling, docking, molecular dynamics (MD) simulations and MM/PBSA binding free energies were employed. The modeled bi-lobed zCDKL1 shares a high degree of secondary structure topology with human orthologs where ATP prefers to lie in the central cavity of the bi-lobed catalytic domain enclosed by strong hydrogen bonding, electrostatic and hydrophobic contacts. Long range MD simulation portrayed that catalytic domain of zCDKL1 to be highly rigid in nature as compared to the complex (zCDKL1-ATP) form. Comparative analysis with its orthologs revealed that conserved amino acids i.e., Ile10, Gly11, Glu12, Val18, Arg31, Phe80, Glu 130, Cys143 and Asp144 were crucial for ATP binding mechanism, which needs further investigation for legitimacy. MM/PBSA method revealed that van der Waals, electrostatic and polar solvation energy mostly contributes towards negative free energy. The implications of ATP binding mechanism inferred through these structural bioinformatics approaches will help in understanding the catalytic mechanisms of important STKs in eukaryotic system.

Identification of Potential MHC Class-II-Restricted Epitopes Derived from Leishmania donovani Antigens by Reverse Vaccinology and Evaluation of Their CD4+ T-Cell Responsiveness against Visceral Leishmaniasis

Visceral leishmaniasis (VL) is one of the most neglected tropical diseases for which no vaccine exists. In spite of extensive efforts, no successful vaccine is available against this dreadful infectious disease. To support vaccine development, an immunoinformatics approach was applied to screen potential MHC class-II-restricted epitopes that can activate the immune cells. Initially, 37 epitopes derived from six stage-dependent, overexpressed antigens were predicted, which were presented by at least 26 diverse MHC class-II allele. Based on a population coverage analysis and human leukocyte antigen cross-presentation ability, six of the 37 epitopes were selected for further analysis. Stimulation with synthetic peptide alone or as a cocktail triggered intracellular IFN-γ production. Moreover, specific IgG antibodies were detected in the serum of active VL cases against P1, P4, P5, and P6 in order to evaluate the peptide effect on the humoral immune response. Additionally, most of the peptides, except P2, were found to be non-inducers of CD4+ IL-10 against both active VL as well as treated VL subjects. This finding suggests there is no role of these peptides in the pathogenesis of Leishmania. Peptide immunogenicity was validated in BALB/c mice immunized with a cocktail of synthetic peptide emulsified in complete Freund’s adjuvant/incomplete Freund’s adjuvant. The immunized splenocytes induced strong spleen cell proliferation upon parasite re-stimulation. Furthermore, increased IFN-γ, interleukin-12, IL-17, and IL-22 production augmented with elevated nitric oxide (NO) synthesis is thought to play a crucial role in macrophage activation. In this investigation, we identified six MHC class-II-restricted epitope hotspots of Leishmania antigens that induce CD4+ Th1 and Th17 responses, which could be used to potentiate a human universal T-epitope vaccine against VL.