Understanding the mechanism of the SARS CoV-2 coinfection with other respiratory viruses

Coinfections have a potential role in increased morbidity and mortality rates during pandemics. Our investigation is aimed at evaluating the viral coinfection prevalence in COVID-19 patients. Rapid diagnostic tests are tools with a paramount impact both on improving patient care. Particularly in the case of respiratory infections, it is of great importance to quickly confirm/exclude the involvement of pathogens. The COVID-19 pandemic has been associated with changes in respiratory virus infections worldwide, which have differed between virus types. In this paper, we systematically searched the percentage of coinfection of various respiratory viruses in COVID-19-positive samples. We included patients of all ages, in all settings. The main outcome was the proportion of patients with viral coinfection. By describing the differences in changes between viral species across different geographies over the course of the COVID-19 pandemic, we may better understand the complex factors involved in the community cocirculation of respiratory viruses.

Genomic characterization and Epidemiology of an emerging SARS-CoV-2 variant in Delhi, India

Delhi, the national capital of India, has experienced multiple SARS-CoV-2 outbreaks in 2020 and reached a population seropositivity of over 50% by 2021. During April 2021, the city became overwhelmed by COVID-19 cases and fatalities, as a new variant B.1.617.2 (Delta) replaced B.1.1.7 (Alpha). A Bayesian model explains the growth advantage of Delta through a combination of increased transmissibility and partial reduction of immunity elicited by prior infection (median estimates; x1.5-fold, 20% reduction). Seropositivity of an employee and family cohort increased from 42% to 86% between March and July 2021, with 27% reinfections, as judged by increased antibody concentration after previous decline. The likely high transmissibility and partial evasion of immunity by the Delta variant contributed to an overwhelming surge in Delhi. One-Sentence SummaryDelhi experienced an overwhelming surge of COVID-19 cases and fatalities peaking in May 2021 as the highly transmissible and immune evasive Delta variant replaced the Alpha variant.

SARS-CoV-2 B.1.617 emergence and sensitivity to vaccine-elicited antibodies

The SARS-CoV-2 B.1.617.2 (Delta) variant was first identified in the state of Maharashtra in late 2020 and spread throughout India, outcompeting pre-existing lineages including B.1.617.1 (Kappa) and B.1.1.7 (Alpha). In vitro, B.1.617.2 is 6-fold less sensitive to serum neutralising antibodies from recovered individuals, and 8-fold less sensitive to vaccine-elicited antibodies as compared to wild type Wuhan-1 bearing D614G. Serum neutralising titres against B.1.617.2 were lower in ChAdOx-1 versus BNT162b2 vaccinees. B.1.617.2 spike pseudotyped viruses exhibited compromised sensitivity to monoclonal antibodies against the receptor binding domain (RBD) and N-terminal domain (NTD), in particular to the clinically approved bamlavinimab and imdevimab monoclonal antibodies. B.1.617.2 demonstrated higher replication efficiency in both airway organoid and human airway epithelial systems as compared to B.1.1.7, associated with B.1.617.2 spike being in a predominantly cleaved state compared to B.1.1.7. Additionally we observed that B.1.617.2 had higher replication and spike mediated entry as compared to B.1.617.1, potentially explaining B.1.617.2 dominance. In an analysis of over 130 SARS-CoV-2 infected healthcare workers across three centres in India during a period of mixed lineage circulation, we observed substantially reduced ChAdOx-1 vaccine efficacy against B.1.617.2 relative to non-B.1.617.2. Compromised vaccine efficacy against the highly fit and immune evasive B.1.617.2 Delta variant warrants continued infection control measures in the post-vaccination era.

Intra-host variability in global SARS-CoV-2 genomes as signatures of RNA editing: implications in viral and host response outcomes

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying Variants Of Concern (VOC). Besides, viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host Single Nucleotide Variations (iSNVs). Analysing 1, 347 samples collected till June 2020, we recorded 18, 146 iSNV sites throughout the SARS-CoV-2 genome. Both, mutations in RdRp as well as APOBEC and ADAR mediated RNA editing seem to contribute to the differential prevalence of iSNVs in hosts. Noteworthy, 41% of all unique iSNVs were reported as SNVs by 30th September 2020 in samples submitted to GISAID, which increased to [~]80% by 30th June 2021. Following this, analysis of another set of 1, 798 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) variations appeared as iSNVs before getting fixed in the population. We also observe hyper-editing events at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=177 SRC="FIGDIR/small/417519v3_ufig1.gif" ALT="Figure 1"> View larger version (41K): org.highwire.dtl.DTLVardef@12b6ac2org.highwire.dtl.DTLVardef@16df897org.highwire.dtl.DTLVardef@dbbec2org.highwire.dtl.DTLVardef@c8de14_HPS_FORMAT_FIGEXP M_FIG C_FIG

Multi-Omics and Integrated Network Approach to Unveil Evolutionary Patterns, Mutational Hotspots, Functional Crosstalk and Regulatory Interactions in SARS-CoV-2

SARS-CoV-2 pandemic resulted in 92 million cases in a span of one year. The study focuses on understanding population specific variations attributing its high rate of infections in specific geographical regions particularly in USA. Rigorous phylogenomic network analysis of complete SARS-CoV-2 genomes (245) inferred five central clades named a (ancestral), b, c, d and e (subtype e1 & e2). The clade d & e2 were found exclusively comprising of USA. Clades were distinguished by 10 co-mutational combinations in Nsp3, ORF8, Nsp13, S, Nsp12, Nsp2 and Nsp6. Our analysis revealed that only 67.46% of SNP mutations were at amino acid level. T1103P mutation in Nsp3 was predicted to increase protein stability in 238 strains except 6 strains which were marked as ancestral type; whereas co-mutation (P409L & Y446C) in Nsp13 were found in 64 genomes from USA highlighting its 100% co-occurrence. Docking highlighted mutation (D614G) caused reduction in binding of Spike proteins with ACE2, but it also showed better interaction with TMPRSS2 receptor contributing to high transmissibility among USA strains. We also found host proteins, MYO5A, MYO5B, MYO5C had maximum interaction with viral proteins (N, S, M). Thus, blocking the internalization pathway by inhibiting MYO5 proteins which could be an effective target for COVID-19 treatment. The functional annotations of the HPI network were found to be closely associated with hypoxia and thrombotic conditions confirming the vulnerability and severity of infection. We also screened CpG islands in Nsp1 & N conferring ability of SARS-CoV-2 to enter and trigger ZAP activity inside host cell. ImportanceIn the current study we presented a global view of mutational pattern observed in SARS-CoV-2 virus transmission. This provided a who-infect-whom geographical model since the early pandemic. This is hitherto the most comprehensive comparative genomics analysis of full-length genomes for co-mutations at different geographical regions specially in USA strains. Compositional structural biology results suggested that mutations have balance of contrary forces effect on pathogenicity suggesting only few mutations to effective at translation level but not all. Novel HPI analysis and CpG predictions elucidates the proof of concept of hypoxia and thrombotic conditions in several patients. Thus, the current study focuses the understanding of population specific variations attributing high rate of SARS-CoV-2 infections in specific geographical regions which may eventually be vital for the most severely affected countries and regions for sharp development of custom-made vindication strategies.

ACE2 Homo-dimerization, Human Genomic variants and Interaction of Host Proteins Explain High Population Specific Differences in Outcomes of COVID19

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a positive single-stranded RNA virus that causes a highly contagious Corona Virus Disease (COVID19). Entry of SARS-CoV-2 in human cells depends on binding of the viral spike (S) proteins to cellular receptor Angiotensin-converting enzyme 2 (ACE2) and on S-protein priming by host cell serine protease TMPRSS2. Recently, COVID19 has been declared pandemic by World Health Organization (WHO) yet high differences in disease outcomes across countries have been seen. We provide evidences to explain these population-level differences. One of the key factors of entry of the virus in host cells presumably is because of differential interaction of viral proteins with host cell proteins due to different genetic backgrounds. Based on our findings, we conclude that a higher expression of ACE2 is facilitated by natural variations, acting as Expression quantitative trait loci (eQTLs), with different frequencies in different populations. We suggest that high expression of ACE2 results in homo-dimerization, proving disadvantageous for TMPRSS2 mediated cleavage of ACE2; whereas, the monomeric ACE2 has higher preferential binding with SARS-CoV-2 S-Protein vis-a-vis its dimerized counterpart. Further, eQTLs in TMPRSS2 and natural structural variations in the gene may also result in differential outcomes towards priming of viral S-protein, a critical step for entry of the Virus in host cells. In addition, we suggest that several key host genes, like SLC6A19, ADAM17, RPS6, HNRNPA1, SUMO1, NACA, BTF3 and some other proteases as Cathepsins, might have a critical role. To conclude, understanding population specific differences in these genes may help in developing appropriate management strategies for COVID19 with better therapeutic interventions.

Assessment of epigenetic alterations and in silico analysis of mutation affecting PTEN expression among Indian cervical cancer patients.

Genetic and epigenetic anomalies accountable for genetic dysregulation are the most common aberrations that determine the underlying heterogeneity of the tumor cells. Currently, phosphatase and tensin homolog (PTEN) incongruity has emerged as potent and persuasive malfunctioning in varied human malignancies. In this study, we have analysed the promoter hypermethylation and expression status of PTEN. We identified different mutations in the exonic region of PTEN. Functional consequences of these mutations were explored using in silico techniques. Promoter hypermethylation of PTEN was detected using methylation-specific polymerase chain reaction (MS-PCR), expression analysis was performed with immunohistochemistry (IHC) and mutation by direct sequencing in a total of 168 uterine cervix tumor cases. The findings were statistically correlated with the clinical parameters. In addition, the effect of nonsynonymous mutations was studied with molecular dynamics simulations. PTEN promoter hypermethylation (45.8%) was found to be significantly associated with the of PTEN loss (57.14%, P < 0.0001). Tumor stages, tumor size, lymph node (LN) were found to be significantly correlated with both PTEN promoter hypermethylation and PTEN loss. Histological grade, however, showed a significant association with only PTEN loss. In total, 11.76% of tumors exhibited mutations in exon 5 and 7, out of which E150K of exon 5 showed the highest deviations in the crystal structure of PTEN by in silico analysis. This study provides valuable insights into oncology and paves the path in the development of efficient biomarker and/or imperative therapeutic tool for cervical cancer treatment.

MicroRNA (hsa-miR-19b-2-5p) targets key mitochondrial biogenesis genes-a bioinformatics analysis.

The present study on the basis of a detailed bioinformatics analysis proposed a potential role of a miRNA, hsa-miR-19b-2-5p, in regulating the mitochondrial biogenesis. The miRNA has shown to be involved in important biological processes of cellular metabolic, cellular macromolecule biosynthetic processes and gene expression pathways. The miRNA, hsa-miR-19b-2-5p, was predicted to regulate the molecular function of nucleic acid, organic/heterocyclic compound, nucleic acid binding transcription factor activity. The pathway enrichment analysis suggested that this miRNA participated in several metabolic pathways which could be a key to the regulation of the mitochondrial gene expression and biogenesis. In addition, this miRNA targets a total of 112 mitochondria-related genes, establishing further the crucial role of the candidate miRNA in mitochondrial biology.

Identification of novel inhibitors against UDP-galactopyranose mutase to combat leishmaniasis.

Leishmania, a protozoan parasite that causes leishmaniasis, affects 1-2 million people every year worldwide. Leishmaniasis is a vector born disease and characterized by a diverse group of clinical syndromes. Current treatment is limited because of drug resistance, high cost, poor safety, and low efficacy. The urgent need for potent agents against Leishmania has led to significant advances in the development of novel antileishmanial drugs. ß-galactofuranose (ß-Galf) is an important component of Leishmanial cell surface matrix and plays a critical role in the pathogenesis of parasite. UDP-galactopyranose mutase (UGM) converts UDP-galactopyranose (UDP-Galp) to UDP-galactofuranose (UDP-Galf) which acts as the precursor for ß-Galf synthesis. Due to its absence in human, this enzyme is selected as the potential target in search of new antileishmanial drugs. Three dimensional protein structure model of Leishmania major UGM (LmUGM) has been homology modeled using Trypanosoma cruzi UGM (TcUGM) as a template. The stereochemistry was validated further. We selected already reported active compounds from PubChem database to target the LmUGM. Three compounds (6064500, 44570814, and 6158954) among the top hit occupied the UDP binding site of UGM suggested to work as a possible inhibitor for it. In vitro antileishmanial activity assay was performed with the top ranked inhibitor, 6064500. The 6064500 molecule has inhibited the growth of Leishmania donovani promastigotes significantly. Further, at similar concentrations it has exhibited significantly lesser toxicity than standard drug miltefosine hydrate in mammalian cells.

Genetic and epigenetic alterations affecting PARK-2 expression in cervical neoplasm among North Indian patients.

The recent investigation on PARK-2, a putative tumor suppressor gene, has found that it has been altered in multiple human malignancies. However, the clinical impact of PARK-2 alteration in uterine cervix carcinoma has not yet been studied. Therefore, we aimed to examine mutations, promoter hypermethylation, and protein expression of PARK-2 among the North Indian patients and their association with clinical parameters to evaluate the implication of PARK-2 in the genesis of cervical cancer. A total of 168 patient samples were processed for mutational analysis by single-strand conformation polymorphism, sequencing, and further in silico analysis of the identified mutations. Promoter hypermethylation by methylation-specific polymerase chain reaction and expression of PARK-2 were performed using immunohistochemistry. Statistical correlation between molecular findings and the clinicopathological parameters was taken to figure out the meaningful outcome. As per our findings, 3.5% (6/168) tumors showed novel missense mutations in exon 11 of PARK-2. In silico analysis showed high structural deviations manifested by mutations, A398D and Y391N, in both mutant proteins as compared to wild type. Promoter hypermethylation was observed in total of 29% of (48/168) tumor samples. Furthermore, 46.43% tumors (78/168) exhibited loss of PARK-2 expression in cervical carcinoma. The loss of expression of PARK-2 when correlated with clinical parameters resulted in significant association with tumor stage (p = 0.002) and with histological grade (p = 0.025). However, only clinical stage remained significant after Bonferroni correction (p < 0.007). A trend was observed between PARK-2 promoter hypermethylation and its protein expression. Our study provided sufficient information and insight for investigation of PARK-2 and highlighted its role as a tumor suppressor gene in cervical cancer in North Indian population.

Mitochondrial ND5 mutation mediated elevated ROS regulates apoptotic pathway epigenetically in a P53 dependent manner for generating pro-cancerous phenotypes.

We have previously observed concomitant events of mutations in mitochondrial and nuclear genes, along with elevated reactive oxygen species (ROS) and differential methylation within the promoters of nuclear genes in tumors and in vitro experiments of tumorigenesis. These observations have made it pertinent to replicate and understand the role of acquired mitochondrial condition in tuning a cell to accomplish a pro-cancerous state. Using a codon optimized vector system for exogenous over-expression and mitochondrial localization; we have characterized here the role of over-expressed wild type mtND5 and one of its non-synonymous somatic mutation, ND5:P265H. The ectopically over-expressed ND5:P265H in mitochondria resulted in a reduced Complex I activity, generation of higher ADP/ATP ratio, reactive oxygen species (ROS) and carbonylation of proteins as compared to mock-transfected cells. Cells over-expressing mtND5 variant produced both peroxide as well as super-oxide ROS; the generation of which was dependent on the functional status of P53; modulating epigenetically the expression of key apoptosis pathway genes. The pro-cancerous phenotypes, of anchorage dependent and independent growth; increased glucose uptake and lactate production, were selectively observed only in P53 non-functional cells over-expressing mutant ND5:P265H. We propose that somatic mutation in mtND5 resulting in down-regulated complex I enzyme activity, elevated ROS and up-regulation of a set of nuclear anti-apoptotic genes epigenetically in the P53 dysfunctional cellular background, has provided a unique understanding of the molecular mechanism of mitochondrial mutation; and the concomitant existence of somatically acquired mitochondrial and nuclear p53 mutations, in cancer progression and promotion.

Understanding Genetic Heterogeneity in Type 2 Diabetes by Delineating Physiological Phenotypes: SIRT1 and its Gene Network in Impaired Insulin Secretion.

Type 2 diabetes (T2D) is a chronic metabolic disease which shows an exponential increase in all parts of the world. However, the disease is controllable by early detection and modified lifestyle. A series of factors have been associated with the pathogenesis of diabetes, and genes are considered to play a critical role. The individual risk of developing T2D is determined by an altered genetic background of the en-zymes involved in several metabolism-related biological mechanisms, including glucose homeostasis, insulin metab-olism, the glucose and ion transporters involved in glucose uptake, transcription factors, signaling intermediates of insulin signaling pathways, insulin production and secretion, pancreatic tissue development, and apoptosis. However, many candidate genes have shown heterogeneity of associations with the disease in different populations. A possible approach to resolving this complexity and under-standing genetic heterogeneity is to delineate the physiological phenotypes one by one as studying them in combination may cause discrepancies in association studies. A systems biology approach involving regulatory proteins, transcription factors, and microRNAs is one way to understand and identify key factors in complex diseases such as T2D. Our earlier studies have screened more than 100 single nucleotide polymorphisms (SNPs) belonging to more than 60 globally known T2D candidate genes in the Indian population. We observed that genes invariably involved in the activity of pancreatic ß-cells provide susceptibility to type 2 diabetes (T2D). Encouraged by these results, we attempted to delineate in this review one of the commonest physiological phenotypes in T2D, namely impaired insulin secretion, as the cause of hyperglycemia. This review is also intended to explain the genetic basis of the pathophysiology of insulin secretion in the context of variations in the SIRT1 gene, a major switch that modulates insulin secretion, and a set of other genes such as HHEX, PGC-α, TCF7L2, UCP2, and ND3 which were found to be in association with T2D. The review aims to look at the genotypic and transcriptional regulatory relationships with the disease phenotype.

Biophysical insight into the anti-amyloidogenic behavior of taurine.

In this work, we investigated the inhibitory ability of taurine on the aggregation of Human serum albumin (HSA) and also examined how it controls the kinetic parameters of the aggregation process. We demonstrated the structural alterations in the HSA after binding to the taurine at 65 °C by exploiting various biophysical techniques. UV-vis spectroscopy was used to check the turbidometric changes in the protein. Thioflavin T fluorescence kinetics was subjected to explore kinetic parameters comparing the amyloid formation in the presence of varying concentration of taurine. Further, Congo red binding and ANS binding assays were performed to determine the inhibitory effect of taurine on HSA fibrillation process and surface hydrophobicity modifications occurring before and after the addition of taurine with protein, respectively. Far UV CD and Dynamic Light Scattering (DLS) confirmed that taurine stabilized the protein α-helical structure and formed complex with HSA which is further supported by differential scanning calorimetry (DSC). Moreover, microscopic imaging techniques were also done to analyze the morphology of aggregation formed. Taurine is also capable of altering the cytotoxicity of the proteinaceous aggregates. Molecular docking study also deciphered the possible residues involved in protein and drug interaction.

In silico screening, genotyping, molecular dynamics simulation and activity studies of SNPs in pyruvate kinase M2.

Role of, 29-non-synonymous, 15-intronic, 3-close to UTR, single nucleotide polymorphisms (SNPs) and 2 mutations of Human Pyruvate Kinase (PK) M2 were investigated by in-silico and in-vitro functional studies. Prediction of deleterious substitutions based on sequence homology and structure based servers, SIFT, PANTHER, SNPs&GO, PhD-SNP, SNAP and PolyPhen, depicted that 19% emerged common between all the mentioned programs. SNPeffect and HOPE showed three substitutions (C31F, Q310P and S437Y) in-silico as deleterious and functionally important. In-vitro activity assays showed C31F and S437Y variants of PKM2 with reduced activity, while Q310P variant was catalytically inactive. The allosteric activation due to binding of fructose 1-6 bisphosphate (FBP) was compromised in case of S437Y nsSNP variant protein. This was corroborated through molecular dynamics (MD) simulation study, which was also carried out in other two variant proteins. The 5 intronic SNPs of PKM2, associated with sporadic breast cancer in a case-control study, when subjected to different computational analyses, indicated that 3 SNPs (rs2856929, rs8192381 and rs8192431) could generate an alternative transcript by influencing splicing factor binding to PKM2. We propose that these, potentially functional and important variations, both within exons and introns, could have a bearing on cancer metabolism, since PKM2 has been implicated in cancer in the recent past.

Apoptosis regulatory protein-protein interaction demonstrates hierarchical scale-free fractal network.

Dysregulation or inhibition of apoptosis favors cancer and many other diseases. Understanding of the network interaction of the genes involved in apoptotic pathway, therefore, is essential, to look for targets of therapeutic intervention. Here we used the network theory methods, using experimentally validated 25 apoptosis regulatory proteins and identified important genes for apoptosis regulation, which demonstrated a hierarchical scale-free fractal protein-protein interaction network. TP53, BRCA1, UBIQ and CASP3 were recognized as a four key regulators. BRCA1 and UBIQ were also individually found to control highly clustered modules and play an important role in the stability of the overall network. The connection among the BRCA1, UBIQ and TP53 proteins was found to be important for regulation, which controlled their own respective communities and the overall network topology. The feedback loop regulation motif was identified among NPM1, BRCA1 and TP53, and these crucial motif topologies were also reflected in high frequency. The propagation of the perturbed signal from hubs was found to be active upto some distance, after which propagation started decreasing and TP53 was the most efficient signal propagator. From the functional enrichment analysis, most of the apoptosis regulatory genes associated with cardiovascular diseases and highly expressed in brain tissues were identified. Apart from TP53, BRCA1 was observed to regulate apoptosis by influencing motif, propagation of signals and module regulation, reflecting their biological significance. In future, biochemical investigation of the observed hub-interacting partners could provide further understanding about their role in the pathophysiology of cancer.

Insight into the effect of inhibitor resistant S130G mutant on physico-chemical properties of SHV type beta-lactamase: a molecular dynamics study.

Bacterial resistance is a serious threat to human health. The production of ß-lactamase, which inactivates ß-lactams is most common cause of resistance to the ß-lactam antibiotics. The Class A enzymes are most frequently encountered among the four ß-lactamases in the clinic isolates. Mutations in class A ß-lactamases play a crucial role in substrate and inhibitor specificity. SHV and TEM type are known to be most common class A ß-lactamases. In the present study, we have analyzed the effect of inhibitor resistant S130G point mutation of SHV type Class-A ß-lactamase using molecular dynamics and other in silico approaches. Our study involved the use of different in silico methods to investigate the affect of S130G point mutation on the major physico-chemical properties of SHV type class A ß-lactamase. We have used molecular dynamics approach to compare the dynamic behaviour of native and S130G mutant form of SHV ß-lactamase by analyzing different properties like root mean square deviation (RMSD), H-bond, Radius of gyration (Rg) and RMS fluctuation of mutation. The results clearly suggest notable loss in the stability of S130G mutant that may further lead to decrease in substrate specificity of SHV. Molecular docking further indicates that S130G mutation decreases the binding affinity of all the three inhibitors in clinical practice.

mtDNA germ line variation mediated ROS generates retrograde signaling and induces pro-cancerous metabolic features.

mtDNA non-synonymous germ line variation (G10398A; p.A114T) has remained equivocal with least mechanistic understanding in showing an association with cancer. This has necessitated showing in-vitro how an over-expression within mitochondria of either of the variants produces higher intracellular ROS, resulting in differential anchorage dependent and independent growth. Both these features were observed to be relatively higher in ND3:114T variant. An elevated amount of intracellular carbonylated proteins and a reduced activity of a key glycolytic enzyme, Pyruvate kinase M2, along with high glucose uptake and lactate production were other pro-cancerous features observed. The retrograde signaling through surplus ROS was generated by post-ND3 over-expression regulated nuclear gene expression epigenetically, involving selectively the apoptotic-DDR-pathways. The feature of ND3 over-expression, inducing ROS mediated pro-cancerous features in the cells in in vitro, was replicated in a pilot study in a limited number of sporadic breast tumors, suggesting the importance of mitochondrial germ-line variant(s) in enabling the cells to acquire pro-cancerous features.

Molecular simulation of Tyr105 phosphorylated pyruvate kinase M2 to understand its structure and dynamics.

Tyrosine phosphorylation (p-Y105) of pyruvate kinase (PK) M2, in recent years, has been suggested to facilitate Warburg effect and tumor cell growth. However, a comparison of the structural dynamics of the un-phosphorylated, the active, and the phosphorylated-at-Y105, the inactive-states, is not clear. We studied molecular dynamics of the two states to unravel these features, where phosphorylated PKM2 showed a rapid global conformation change in the initial stages of the simulation. The overall simulation identified that the phosphorylation event results in more buried and less flexible PKM2 conformation, as compared to the un-phosphorylated form, resulting in an open and closed conformation of the active site in un-phosphorylated and phosphorylated forms, respectively, due to the movement of B domain. This conformational shift in Y105-phosphorylated-PKM2 (p-Y105-PKM2) with closed active site, responsible for inhibition of PKM2 activity, was an outcome of the bending residues (117-118, 218-219, 296-297, and 301-308) within the loop connecting A and B domains and the presence of helix-loop-helix motif in A domain. The un-phosphorylated PKM2 formed a helix bend (H4) due to less fluctuation of the residue S-100; where the other end of the helix (H4) was connected to the substrate binding pocket. Further, simulation analysis showed that phosphorylation did not affect the FBP binding predominantly. We propose that p-Y105 inhibits the activity of PKM2 without influencing FBP binding directly and not allowing the open binding conformation by influencing G128, S100, G506 and gamma turn, G126 and S127 residues. Phosphorylated PKM2 was also identified to gain the transcriptional factor function which was not the case with un-phosphorylated form. These structurally important residues in PKM2 could have a bearing on cancer metabolism, since PKM2 has been implicated in the promotion of cancer in the recent past.

PARK2 and proinflammatory/anti-inflammatory cytokine gene interactions contribute to the susceptibility to leprosy: a case-control study of North Indian population.

OBJECTIVES: Cytokines and related molecules in immune-response pathways seem important in deciding the outcome of the host-pathogen interactions towards different polar forms in leprosy. We studied the role of significant and functionally important single-nucleotide polymorphisms (SNPs) in these genes, published independently from our research group, through combined interaction with an additional analysis of the in silico network outcome, to understand how these impact the susceptibility towards the disease, leprosy. DESIGN: The study was designed to assess an overall combined contribution of significantly associated individual SNPs to reflect on epistatic interactions and their outcome in the form of the disease, leprosy. Furthermore, in silico approach was adopted to carry out protein-protein interaction study between PARK2 and proinflammatory/anti-inflammatory cytokines. SETTING: Population-based case-control study involved the data of North India. Protein-protein interaction networks were constructed using cytoscape. PARTICIPANTS: Study included the data available from 2305 Northern Indians samples (829 patients with leprosy; 1476 healthy controls), generated by our research group. PRIMARY AND SECONDARY OUTCOME MEASURES: For genotype interaction analysis, all possible genotype combinations between selected SNPs were used as an independent variable, using binary logistic regression with the forward likelihood ratio method, keeping the gender as a covariate. RESULTS: Interaction analysis between PARK2 and significant SNPs of anti-inflammatory/proinflammatory cytokine genes, including BAT1 to BTNL2-DR spanning the HLA (6p21.3) region in a case-control comparison, showed that the combined analysis of: (1) PARK2, tumour necrosis factor (TNF), BTNL2-DR, interleukin (IL)-10, IL-6 and TGFBR2 increased the risk towards leprosy (OR=2.54); (2) PARK2, BAT1, NFKBIL1, LTA, TNF-LTB, IL12B and IL10RB provided increased protection (OR=0.26) in comparison with their individual contribution. CONCLUSIONS: Epistatic SNP-SNP interactions involving PARK2 and cytokine genes provide an additive risk towards leprosy susceptibility. Furthermore, in silico protein-protein interaction of PARK2 and important proinflammatory/anti-inflammatory molecules indicate that PARK2 is central to immune regulation, regulating the production of different cytokines on infection.

Structure-based screening of inhibitors against KPC-2: designing potential drug candidates against multidrug-resistant bacteria.

KPC-2 ß-lactamase demonstrates a wide substrate spectrum that includes carbapenamases, oxyimino-cephalosporins, and cephamycins. In addition, strains harboring KPC-type ß-lactamases are often identified as resistant to standard ß-lactamase inhibitors. Thus, KPC-2 carbapenems present a significant clinical challenge, as the mechanistic bases for KPC-2-associated phenotypes remain mysterious. Inhibiting the function of these resistance enzymes could control the hydrolysis of antibiotics. In the present study, we have reported two novel (non-ß-lactatam) compounds that inhibit the activity of the KPC-2 enzyme. These compounds were identified by structure-based virtual screening using computational docking programs and molecular dynamics simulations with the solved crystal structure. Two compounds (ZINC01807204 and ZINC02318494) were selected on the basis of fitness scores from docking program and 5 ns molecular dynamics simulations. These commercially available compounds have been procured and their biological activity was experimentally evaluated on the E. coli strain carrying recombinant KPC-2. These new compounds in combination with ceftazidime and cefoxitin exhibited the Minimum Inhibitory Concentration (MIC) values of 2 and 8 µg/ml respectively, which were found to be lower as compared to known ß-lactamase inhibitors. Moreover, these compounds were also found to have comparable MICs values being 64 µg/ml in combination with ceftriaxone. This study explored novel inhibitors against KPC-2, a class A ß-lactamase, which may be putative drug candidates against KPC-2 producing bacterial infection.