Search

The last decade has witnessed the rise of an extremely threatening healthcare-associated multidrug-resistant non-albicans Candida (NAC) species, Candida auris. Since besides target alterations, efflux mechanisms contribute maximally to antifungal resistance, it is imperative to investigate their contributions in this pathogen. Of note, within the major facilitator superfamily (MFS) of efflux pumps, drug/H

The peptide transport (PTR) or proton-dependent oligopeptide transporter (POT) family exploits the inwardly directed proton motive force to facilitate the cellular uptake of di/tripeptides. Interestingly, some representatives are also shown to import peptide-based antifungals in certain Candida species. Thus, the identification and characterization of PTR transporters serve as an essential first step for their potential usage as antifungal peptide uptake systems. Herein, we present a genome-wide inventory of the PTR transporters in five prominent Candida species. Our study identifies 2 PTR transporters each in C. albicans and C. dubliniensis, 1 in C. glabrata, 4 in C. parapsilosis, and 3 in C. auris. Notably, despite all representatives retaining the conserved features seen in the PTR family, there exist two distinct classes of PTR transporters that differ in terms of their sequence identities and lengths of certain extracellular and intracellular segments. Further, we also evaluated the contribution of each PTR protein of the newly emerged multi-drug-resistant C. auris in di/tripeptide uptake. Notably, deletion of two PTR genes BNJ08_003830 and BNJ08_005124 led to a marked reduction in the transport capabilities of several tested di/tripeptides. However, all three genes could complement the role of native PTR2 gene of Saccharomyces cerevisiae, albeit to varied levels. Besides, BNJ08_005124 deletion also resulted in increased resistance toward the peptide-nucleoside drug Nikkomycin Z as well as the glucosamine-6-phosphate synthase inhibitor, L-norvalyl-N3-(4-methoxyfumaroyl)-L-2,3-diaminopropionoic acid (Nva-FMDP), pointing toward its predominant role in their uptake mechanism. Altogether, the study provides an important template for future structure-function investigations of PTR transporters in Candida species. KEY POINTS: • Candida genome encodes for two distinct classes of PTR transporters. • Candida auris encodes for 3 PTR transporters with different specificities. • BNJ08_005124 in C. auris is involved in the uptake of Nikkomycin Z and Nva-FMDP.

Antifungal resistance mediated by overexpression of ABC transporters is one of the primary roadblocks to effective therapy against Candida infections. Thus, identification and characterization of the ABC transporter repertoire in Candida species are of high relevance. The method described in the chapter is based on our previously developed bioinformatic pipeline for identification of ABC proteins in Candida species. The methodology essentially involves the utilization of a hidden Markov model (HMM) profile of the nucleotide-binding domain (NBD) of ABC proteins to mine these proteins from the proteome of Candida species. Further, a widely used tool to predict membrane protein topology is exploited to identify the true transporter candidates out of the ABC proteins. Even though the chapter specifically focuses on a method to identify ABC transporters in Candida auris , the same can also be applied to any other Candida species.

Known crystal structures of 10 cytochromes (cyts) c from different sources led to the conclusion that natural mutations in these proteins does not affect their 3D structure, hence evolution preserved structure for function. A sequence alignment of horse cyt c with all other 284 cyts c led to two important conclusions: (i) Leu at position 94 is conserved in all 30 mammalian known sequences, and (ii) there are 14 other species which have either Val or Ile at 94th position. We asked a question: Is the avoidance of substitution by Val or Ile at position 94 in the mammalian cyts c by design or by chance? To answer this question, we introduced natural substitutes of Leu94 by Val and Ile in horse cyt c using site-directed mutagenesis. Here, from our in vitro and molecular dynamic simulation studies on L94V and L94I mutants, we concluded that (i) although the natural mutations destabilize the wild type cyt c, it does not significantly affect the mechanism of folding of the protein, (ii) urea-induced denaturation of WT cyt c and its mutants is a two-state process, and (iii) denaturation of WT cyt c and its mutants by guanidinium chloride is not a two-state process.

Aberrant misfolding and aggregation of TAR DNA-binding protein 43 (TDP-43) and its fragments have been implicated in amyotrophic lateral sclerosis (ALS) and other neurodegenerative diseases. Within the protein, the second RNA recognition motif (RRM2) has recently been demonstrated to be a major contributor towards aggregation and the resultant toxicity. However, the physicochemical mechanism of its misfolding from the functional folded state is poorly understood. In the present work, we have used a cumulative ∼2µs of molecular dynamics (MD) simulation to study the structural and thermodynamic characteristics of different unfolded intermediates of RRM2 domain of TDP-43. In 6 M GdmCl at 400 K, at RMSD around 1.5 nm, part of the secondary structure i.e. helix still does not melt without significant change in solvent accessibility and intra-protein hydrogen bonds. However, hydrophobic contacts disrupt significantly suggesting that unfolding proceeds through disruption of hydrophobic core of the protein.The temperature dependent free-energy landscapes (FELs) reveal the presence of multiple metastable intermediate states stabilized by hydrophobic (ILV) contacts and hydrogen bonds. These conformational states have all the native helices intact with significant loss of β-sheets. These partially unfolded states are quite compact and characterized by the exposure of aggregation-prone β-sheets, suggesting the increased aggregation propensity of the partially unfolded states. Our results will thus serve to uncover the structural properties of partially unfolded intermediate states that drive TDP-43 misfolding and aggregation. Elucidating the structural characterization of the misfolding and aggregation prone intermediate states of TDP-43 are important to understand its role in ALS and other neurodegenerative diseases. Communicated by Ramaswamy H. Sarma

A novel series of benzothiazole–piperazine hybrids were rationally designed, synthesized, and evaluated as multifunctional ligands against Alzheimer's disease (AD). The synthesized hybrid molecules illustrated modest to strong inhibition of acetylcholinesterase (AChE) and Aβ1-42 aggregation. Compound 12 emerged as the most potent hybrid molecule exhibiting balanced functions with effective, uncompetitive and selective inhibition against AChE (IC50 = 2.31 μM), good copper chelation, Aβ1-42 aggregation inhibition (53.30%) and disaggregation activities. Confocal laser scanning microscopy and TEM analysis also validate the Aβ fibril inhibition ability of this compound. Furthermore, this compound has also shown low toxicity and is capable of impeding loss of cell viability elicited by H2O2 neurotoxicity in SHSY-5Y cells. Notably, compound 12 significantly improved cognition and spatial memory against scopolamine-induced memory deficit in a mouse model. Hence, our results corroborate the multifunctional nature of novel hybrid molecule 12 against AD and it may be a suitable lead for further development as an effective therapeutic agent for therapy in the future.

Thromboembolic diseases are a major cause of mortality in human and the currently available anticoagulants are associated with various drawbacks, therefore the search for anticoagulants that have better safety profile is highly desirable. Compounds that are part of the dietary routine can be modified to possibly increase their anticoagulant potential. We show mannose 2,3,4,5,6

The PTR or proton-dependent oligopeptide transporter (POT) family exploits the inwardly directed proton motive force to facilitate the cellular uptake of di- and tripeptides. Interestingly, representatives from the family have shown efficacy in delivering peptide-based antifungal derivatives in certain Candida species. Given the increased incidences of fungal infections by Candida species and the associated escalating orders of resistance against frontline antifungals, peptide derivatives are attractive therapeutic options. In that direction, the identification and characterization of PTR transporters serve as an essential first step in the translation of peptide-based antifungals as next-gen therapeutics. Herein, we present a genome-wide inventory of the PTR transporters in five prominent Candida species. Our study identifies 2 PTR transporters each in C. albicans and C. dubliniensis, 1 in C. glabrata, 4 in C. parapsilosis, and 3 in C. auris. Notably, despite all representatives retaining the conserved features seen in the PTR family, there exist two distinct classes of PTR transporters that differ in terms of their sequence identities and membrane topology. Further, we also evaluated the contribution of each PTR protein of the newly emerged multidrug-resistant C. auris in di-/tripeptide uptake. Notably, deletion of the PTR transporters encoded by BNJ08_003830 and BNJ08_005124 led to a marked reduction in the transport capabilities of several tested di-/tripeptides. Besides, BNJ08_005124 deletion also resulted in increased resistance towards the peptide-nucleoside drug Nikkomycin Z, pointing towards its predominant role in the uptake mechanism. Altogether, the study provides an important template for future structure-function investigations of PTR transporters in Candida species.

The last decade has witnessed the rise of extremely threatening healthcare-associated multidrug resistant non-albicans Candida (NAC) species, Candida auris. Thus, understanding the molecular basis of antifungal resistance has emerged as the single most important goal amongst the research community. Since besides target alterations, efflux mechanisms contribute maximally to antifungal resistance, it is imperative to investigate their contributions in this pathogen. Of note, within the Major Facilitator Superfamily (MFS) of efflux pumps, Drug/H+ antiporter family 1 (DHA1) has been established as a predominant contributor towards xenobiotic efflux. Our study provides a complete landscape of DHA1 transporters encoded in the genome of C. auris. This study identifies 14 DHA1 transporters encoded in the genome of the pathogen. We also construct deletion and heterologous overexpression strains for the most important DHA1 drug transporter, viz., CauMdr1 to map the spectrum of its substrates. While the knockout strain did not show any significant changes in the resistance patterns against the majority of the tested substrates, the ortholog when overexpressed in a minimal background S. cerevisiae strain, AD1-8u-, showed significant enhancement in the Minimum inhibitory concentrations (MICs) against a large panel of antifungal molecules. Altogether, the present study provides a comprehensive template for investigating the role of DHA1 members of C. auris in antifungal resistance mechanisms.

Thromboembolic diseases are a major cause of mortality in human and the currently available anticoagulants are associated with various drawbacks, therefore the search for anticoagulants that have better safety profile is highly desirable. Compounds that are part of the dietary routine can be modified to possibly increase their anticoagulant potential. We show mannose 2,3,4,5,6-O-pentasulfate (MPS) as a synthetically modified form of mannose that has appreciable anticoagulation properties. An in silico study identified that mannose in sulfated form can bind effectively to the heparin-binding site of antithrombin (ATIII) and heparin cofactor II (HCII). Mannose was sulfated using a simple sulfation strategy-involving triethylamine-sulfur trioxide adduct. HCII and ATIII were purified from human plasma and the binding analysis using fluorometer and isothermal calorimetry showed that MPS binds at a unique site. A thrombin inhibition analysis using the chromogenic substrate showed that MPS partially enhances the activity of HCII. Further an assessment of in vitro blood coagulation assays using human plasma showed that the activated partial thromboplastin time (APTT) and prothrombin time (PT) were prolonged in the presence of MPS. A molecular dynamics simulation analysis of the HCII-MPS complex showed fluctuations in a N-terminal loop and the cofactor binding site of HCII. The results indicate that MPS is a promising lead due to its effect on the in vitro coagulation rate. Communicated by Ramaswamy H. Sarma

TDP-43 proteinopathies is a disease hallmark that characterizes amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD). The N-terminal domain of TDP-43 (NTD) is important to both TDP-43 physiology and TDP-43 proteinopathy. However, its folding and dimerization process is still poorly characterized. In the present study, we have investigated the folding/unfolding of NTD employing all-atom molecular dynamics (MD) simulations in 8 M dimethylsulfoxide (DMSO) at high temperatures. The MD results showed that the unfolding of the NTD at high temperature evolves through the formation of a number of conformational states differing in their stability and free energy. The presence of structurally heterogeneous population of intermediate ensembles was further characterized by the different extents of solvent exposure of Trp80 during unfolding. We suggest that these non-natives unfolded intermediate ensembles may facilitate NTD oligomerization and subsequently TDP-43 oligomerization, which might lead to the formation of irreversible pathological aggregates, characteristics of disease pathogenesis.

Mycolic acids are the key constituents of mycobacterial cell wall, which protect the bacteria from antibiotic susceptibility, helping to subvert and escape from the host immune system. Thus, the enzymes involved in regulating and biosynthesis of mycolic acids can be explored as potential drug targets to kill

ABC transporter proteins are involved in active transport, both in prokaryotes and eukaryotes. Sequence analysis of nucleotide binding domains (NBDs) of ABC proteins from all taxa revealed a well-conserved new motif having the signature: xT/ShxE/DNhxF, located between Q-loop and ABC signature sequence. A recent structure of an ABC transporter, ABCG5/G8 highlighted the motif as an essential structural determinant of inter-domain crosstalk and termed it as E-helix. We carried out an extensive computational analysis to unravel important structural entities alongside E-helix which plausibly play role in the interlocking mechanism of NBD with TMD. We identified E-helix to be a central structural moiety which interacts with three helices and an intracellular loop that leads to the transmembrane domain. Considering its wide occurrence, we examined the importance of this motif in one representative multidrug ABC transporter of Candida albicans, Cdr1p. The motif residues were replaced by alanines both individually as well as in combinations. The GFP-tagged versions of mutant proteins were overexpressed in Saccharomyces cerevisiae. Overall, our mutational data suggested that this motif plays a role in the maintenance of proper structural fold and/or inter-domain contacts in Cdr1p. We, thus, unveil an essential structural motif in ABC superfamily transporters.

Background Growing resistance to antimicrobials has become an important health issue of the 21st century. Many international, national and local approaches are being employed for the control and prevention of antimicrobial resistance (AMR). Among them, surveillance is reported to be the best method to reduce the spread of infection and thereby AMR. An integral component of AMR surveillance is the informatics suite for collection, storage and analysis of surveillance data. Methods Considering the traits of an optimal surveillance tool and constraints with existing tools, Indian Council of Medical Research (ICMR) initiated the design and development of ICMR’s Antimicrobial Resistance Surveillance system (i-AMRSS). i-AMRSS is a web-based tool built using modular architecture. It is capable of collecting standardized data from small laboratories to generate local and nationwide reports. Results i-AMRSS is a robust, comprehensive, modular, extensible and intelligent open-source tool piloted in ICMR’s AMR Network (31 hospitals and laboratories across India) since 2016. The developed tool has collected more than 280 000 patient records to date. Conclusions The standardized data collected through i-AMRSS would be valuable for various collaborators to monitor outbreaks and infection control practices, evaluate transmission dynamics and formulate antibiotic use and selling policies. The tool is presently being used to capture human testing and consumption data, however, it can be extended for AMR surveillance using a ‘One Health’ approach.

The success of a pathogen within the host depends on various extrinsic factors that work in a synergistic mechanism to promote pathogenesis. One such factor is driven by the changes observed within the host genome, providing survival and establishment of pathogens inside the host. Pathogens are also known for establishing their intracellular niche within the host by mimicking the host enzymes and immune system for survival. Understanding the strategies used by pathogens to intervene in host genetic machinery for pathogenesis is important for creating successful targets and personalized drugs to counterbalance their effects. Accumulation of omics data and simultaneous development of bioinformatics analysis tools have allowed researchers to understand the interplay between prokaryotic and eukaryotic cells through the multi-omics approach. This permits a better understanding of diseases associated with host-parasite interactions and subsequent development of personalized medicines as therapeutics.

The chapter describes the development of a virtual screening workflow involving open-source docking algorithms, consensus scoring, and rescoring using binding energy estimation to identify hits in the drug discovery process. Plasmodium falciparum dihydrofolate reductase (PfDHFR-TS, PDB 1J3I) is used as an example to describe the Virtual Screening (VS) workflow.

Skin cancer is one of the major public health problems worldwide. Its early diagnosis and timely therapy are immensely important in improving patient health. Thus, the improved and accessible diagnostic systems for skin cancer are the most potent determinant of getting the right treatment at the right time. The last two decades have seen unprecedented growth in the application of artificial intelligence (AI) for skin cancer research. Recent advancement in the computational power, digitization of medical imaging, rise of -omics data have accumulated a new opportunity. The ability of AI methods to detect hidden or unknown patterns from such complex datasets reveals their importance. The AI approach in skin cancer has helped to improve the diagnostic and therapeutic strategies, from risk assessment using genomic sequences, accessible smartphone-based “apps” for diagnosis to predict the likelihood of therapy response amidst others. This technology holds a promising potential to automate and assist primary clinicians in improving patient health outcomes through effective diagnostic and therapy strategies using the complex healthcare data paving in the way for precision medicine.

1.ABSTRACTUnderstanding pH-dependent protein stability is important in biological - transport, storage, and delivery, in vivo conditions such as aggregation potential in neurodegenerative disease, and in studying the folding/unfolding of proteins. Using computer simulations, we can replace complex experimental determination and provide an atomistic-level interpretation of the cause and effect of pH on protein stability. Here, we standardize a method that provides a framework through which we examined pH-dependent transient conformations during unfolding simulations of proteins. Constant pH simulations utilized in the prediction of pKa values of charged groups of the peptide. The calculated pKa values employed to fix the appropriate protonation state of the amino acid to simulate the effect of pH on the system. Trajectories from multiple high-temperature MD simulations of the protein sample the conformational space during unfolding for a given pH state. The ensemble of conformations is analyzed from its free energy landscape to identify transient and stable conformations both at a given pH and between different pH. As a test system RN80, a protein fragment analog of the C-peptide from bovine pancreatic ribonuclease-A used to measure the accuracy of the predictions from simulations. Experimental measures of the helix content determined as a function of pH display a bell-shaped curve, i.e. RN80 alpha-helix formation is maximum at pH5 with a subsequent loss in helicity at higher and lower pH. The main forces stabilizing the alpha-helix are a salt-bridge formed between Glu-2 and Arg-10 and cation-pi-interaction between Tyr-8 and His-12. Our protocol includes constant pH calculations, optimal high-temperature simulations, and Free Energy landscape analysis exhibited the agreement with the experimental observations.

The global pandemic crisis, coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has claimed the lives of millions of people across the world. Development and testing of anti-SARS-CoV-2 drugs or vaccines have not turned to be realistic within the timeframe needed to combat this pandemic. Here, we report a comprehensive computational approach to identify the multi-targeted drug molecules against the SARS-CoV-2 proteins, whichare crucially involved in the viral–host interaction, replication of the virus inside the host, disease progression and transmission of coronavirus infection. Virtual screening of 75 FDA-approved potential antiviral drugs against the target proteins, spike (S) glycoprotein, human angiotensin-converting enzyme 2 (hACE2), 3-chymotrypsin-like cysteine protease (3CLpro), cathepsin L (CTSL), nucleocapsid protein, RNA-dependent RNA polymerase (RdRp) and non-structural protein 6 (NSP6), resulted in the selection of seven drugs which preferentially bind to the target proteins. Further, the molecular interactions determined by molecular dynamics simulation revealed that among the 75 drug molecules, catechin can effectively bind to 3CLpro, CTSL, RBD of S protein, NSP6 and nucleocapsid protein. It is more conveniently involved in key molecular interactions, showing binding free energy (ΔGbind) in the range of −5.09 kcal/mol (CTSL) to −26.09 kcal/mol (NSP6). At the binding pocket, catechin is majorly stabilized by the hydrophobic interactions, displays ΔEvdW values: −7.59 to −37.39 kcal/mol. Thus, the structural insights of better binding affinity and favorable molecular interaction of catechin toward multiple target proteins signify that catechin can be potentially explored as a multi-targeted agent against COVID-19., Graphical Abstract Graphical Abstract

In view of the multiple clinical and physiological implications of ABC transporter proteins, there is a considerable interest among researchers to characterize them functionally. However, such characterizations are based on the premise that ABC proteins are accurately identified in the proteome of an organism, and their topology is correctly predicted. With this objective, we have developed ABC-finder, i.e., a Docker-based package for the identification of ABC proteins in all organisms, and visualization of the topology of ABC proteins using a web browser. ABC-finder is built and deployed in a Linux container, making it scalable for many concurrent users on our servers and enabling users to download and run it locally. Overall, ABC-finder is a convenient, portable, and platform-independent tool for the identification and topology prediction of ABC proteins.ABC-finder is accessible at http://abc-finder.osdd.jnu.ac.in

Background A structural phenomenon seen in certain lineages of angiosperms that has captivated many scholars including Charles Darwin is the evolution of plant carnivory. Evidently, these structural features collectively termed carnivorous syndrome, evolved to aid nutritional acquisition from attracted, captured and digested prey. We now understand why plant carnivory evolved but how carnivorous plants acquired these attributes remains a mystery. In an attempt to understand the evolution of Nepenthes pitcher and to shed more light on its role in prey digestion, we analyzed the transcriptome data of the highly specialized Nepenthes khasiana leaf comprising the leaf base lamina, tendril and the different parts/zones of the pitcher tube viz. digestive zone, waxy zone and lid. Results In total, we generated around 262 million high-quality Illumina reads. Reads were pooled, normalized and de novo assembled to generate a reference transcriptome of about 412,224 transcripts. We then estimated transcript abundance along the N. khasiana leaf by mapping individual reads from each part/zone to the reference transcriptome. Correlation-based hierarchical clustering analysis of 27,208 commonly expressed genes indicated functional relationship and similar cellular processes underlying the development of the leaf base and the pitcher, thereby implying that the Nepenthes pitcher is indeed a modified leaf. From a list of 2386 differentially expressed genes (DEGs), we identified transcripts encoding key enzymes involved in prey digestion and protection against pathogen attack, some of which are expressed at high levels in the digestive zone. Interestingly, many of these enzyme-encoding genes are also expressed in the unopened N. khasiana pitcher. Transcripts showing homology to both bacteria and fungi were also detected; and in the digestive zone, fungi are more predominant as compared to bacteria. Taking cues from histology and scanning electron microscopy (SEM) photomicrographs, we found altered expressions of key regulatory genes involved in leaf development. Of particular interest, the expression of class III HOMEODOMAIN-LEUCINE ZIPPER (HD-ZIPIII) and ARGONAUTE (AGO) genes were upregulated in the tendril. Conclusions Our findings suggest that N. khasiana pitchers employ a wide range of enzymes for prey digestion and plant defense, harbor microbes and probably evolved through altered expression of leaf polarity genes.

Background Growing resistance to antimicrobials has become an important health issue of the 21st century. Many international, national and local approaches are being employed for the control and prevention of antimicrobial resistance (AMR). Among them, surveillance is reported to be the best method to reduce the spread of infection and thereby AMR. An integral component of AMR surveillance is the informatics suite for collection, storage and analysis of surveillance data. Methods Considering the traits of an optimal surveillance tool and constraints with existing tools, Indian Council of Medical Research (ICMR) initiated the design and development of ICMR’s Antimicrobial Resistance Surveillance system (i-AMRSS). i-AMRSS is a web-based tool built using modular architecture. It is capable of collecting standardized data from small laboratories to generate local and nationwide reports. Results i-AMRSS is a robust, comprehensive, modular, extensible and intelligent open-source tool piloted in ICMR’s AMR Network (31 hospitals and laboratories across India) since 2016. The developed tool has collected more than 280 000 patient records to date. Conclusions The standardized data collected through i-AMRSS would be valuable for various collaborators to monitor outbreaks and infection control practices, evaluate transmission dynamics and formulate antibiotic use and selling policies. The tool is presently being used to capture human testing and consumption data, however, it can be extended for AMR surveillance using a ‘One Health’ approach.

Background Homology based methods are one of the most important and widely used approaches for functional annotation of high-throughput microbial genome data. A major limitation of these methods is the absence of well-characterized sequences for certain functions. The non-homology methods based on the context and the interactions of a protein are very useful for identifying missing metabolic activities and functional annotation in the absence of significant sequence similarity. In the current work, we employ both homology and context-based methods, incrementally, to identify local holes and chokepoints, whose presence in the Mycobacterium tuberculosis genome is indicated based on its interaction with known proteins in a metabolic network context, but have not been annotated. We have developed two computational procedures using network theory to identify orphan enzymes (‘Hole finding protocol’) coupled with the identification of candidate proteins for the predicted orphan enzyme (‘Hole filling protocol’). We propose an integrated interaction score based on scores from the STRING database to identify candidate protein sequences for the orphan enzymes from M. tuberculosis, as a case study, which are most likely to perform the missing function. Results The application of an automated homology-based enzyme identification protocol, ModEnzA, on M. tuberculosis genome yielded 56 novel enzyme predictions. We further predicted 74 putative local holes, 6 choke points, and 3 high confidence local holes in the genome using ‘Hole finding protocol’. The ‘Hole-filling protocol’ was validated on the E. coli genome using artificial in-silico enzyme knockouts where our method showed 25% increased accuracy, compared to other methods, in assigning the correct sequence for the knocked-out enzyme amongst the top 10 ranks. The method was further validated on 8 additional genomes. Conclusions We have developed methods that can be generalized to augment homology-based annotation to identify missing enzyme coding genes and to predict a candidate protein for them. For pathogens such as M. tuberculosis, this work holds significance in terms of increasing the protein repertoire and thereby, the potential for identifying novel drug targets.

Schizosaccharomyces pombe RNA polymerase II comprises twelve different subunits. Its Rpb9 subunit comprises 113 amino acids, and is the only non-essential subunit of S. pombe RNA polymerase II. However, its functions have not been studied in S. pombe. The results presented in this study demonstrate that Rpb9 is involved in regulating growth under optimum and certain stress conditions in S. pombe. To further address the role (s) of various domains of this subunit in regulating these phenotypes, deletion mutant analysis was done. We observed that the region spanning 1-74 amino acids, encompassing the amino-terminal zinc finger domain and the linker region of Rpb9 was able to rescue the phenotypes associated with rpb9+deletion. We also demonstrate that the functions of this subunit are only partially conserved among yeast and humans. Our computational biology approaches provide a structural basis for the differential role of various Rpb9 domains in S. pombe. Furthermore, using these tools we show that there has been a co-evolution of the interaction residues between the Rpb9 subunit and the two largest subunits of RNA polymerase II, allowing for a more stringent organism-specific packing. Taken together, our results have provided functional and structural insights into the Rpb9 subunit of S. pombe.

In Mycobacterium tuberculosis (MTB), the cell wall synthesis flavoenzyme decaprenylphosphoryl-β- d -ribose 2′-epimerase (DprE1) plays a crucial role in host pathogenesis, virulence, lethality and survival under stress. The emergence of different variants of drug resistant MTB are a major threat worldwide which essentially requires more effective new drug molecules with no major side effects. Here, we used structure based virtual screening of bioactive molecules from the ChEMBL database targeting DprE1, having bioactive 78,713 molecules known for anti-tuberculosis activity. An extensive molecular docking, binding affinity and pharmacokinetics profile filtering results in the selection four compounds, C5 (ChEMBL2441313), C6 (ChEMBL2338605), C8 (ChEMBL441373) and C10 (ChEMBL1607606) which may explore as potential drug candidates. The obtained results were validated with thirteen known DprE1 inhibitors. We further estimated the free-binding energy, solvation and entropy terms underlying the binding properties of DprE1-ligand interactions with the implication of MD simulation, MM/GBSA, MM/PBSA and MM/3D-RISM. Interestingly, we find that C6 shows the highest ΔG scores (−41.28 ± 3.51, −22.36 ± 3.17, −10.33 ± 5.70 kcal mol−1) in MM/GBSA, MM/PBSA and MM/3D-RISM assay, respectively. Whereas, the lowest ΔG scores (−35.31 ± 3.44, −13.67 ± 2.65, −3.40 ± 4.06 kcal mol−1) observed for CT319, the inhibitor co-crystallized with DprE1. Collectively, the results demonstrated that hit-molecules: C5, C6, C8 and C10 having better binding free energy and molecular affinity as compared to CT319. Thus, we proposed that selected compounds may be explored as lead molecules in MTB therapy.

The global pandemic crisis, COVID-19 caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has claimed the lives of millions of people across the world. Development and testing of anti-SARS-CoV-2 drugs or vaccines, are not turned to be realistic in the timeframe needed to combat this pandemic. Thus, rigorous efforts are still ongoing for the drug repurposing as a clinical treatment strategy to control COVID-19. Here we report a comprehensive computational approach to identify the multi-targeted drug molecules against the SARS-CoV-2 proteins, which are crucially involved in the viral-host interaction, replication of the virus inside the host, disease progression and transmission of coronavirus infection. Virtual screening of 72 FDA approved potential antiviral drugs against the target proteins: Spike (S) glycoprotein, human angiotensin-converting enzyme 2 (hACE2), 3-chymotrypsin-like cysteine protease (3CLpro), Cathepsin L, Nucleocapsid protein, RNA-dependent RNA polymerase (RdRp) and nonstructural protein 6 (NSP6) resulted in the selection of seven drugs which preferentially binds to the target proteins. Further, the molecular interactions determined by MD simulation, free energy landscape and the binding free energy estimation, using MM-PBSA revealed that among 72 drug molecules, catechin (flavan-3-ol) can effectively bind to 3CLpro, Cathepsin L, RBD of S protein, NSP-6, and Nucleocapsid protein. It is more conveniently involved in key molecular interactions, showing binding free energy (ΔGbind) in the range of -5.09 kcal/mol (Cathepsin L) to -26.09 kcal/mol (NSP6). At the binding pocket, catechin is majorly stabilized by the hydrophobic interactions, displays ΔEvdW values -7.59 to -37.39 kcal/mol. Thus, the structural insights of better binding affinity and favourable molecular interaction of catechin towards multiple target proteins, signifies that catechin can be potentially explored as a multitargeted agent in the rational design of effective therapies against COVID-19.

A novel series of benzothiazole-piperazine hybrids were rationally designed, synthesized, and evaluated as multifunctional ligands against Alzheimer's disease (AD). The synthesized hybrid molecules illustrated modest to strong inhibition of acetylcholinesterase (AChE) and Aβ

ATP-binding cassette (ABC) is one of the largest superfamily of proteins, which are ubiquitously present, performing variety of cellular functions. These proteins as drug transporters have been enticing substantial consideration because of their clinical importance. The present study focuses on genome wide identification of ABC proteins of an important halotolerant yeast Debaryomyces hansenii and explores their role in salt and drug tolerance. Our bioinformatics analysis identified a total of 30 putative ABC protein-coding genes whose expression at transcript level was confirmed by qRT-PCR. Our comparative phylogenetic analysis of nucleotide binding domains of D. hansenii and topology prediction categorized these proteins into six subfamilies; ABCB/MDR, ABCC/MRP, ABCD/ALDP, ABCF/YEF3, ABCE/RLI, and ABCG/PDR based on the nomenclature adopted by the Human Genome Organization (HUGO). Further, our transmembrane domain (TMD) predictions suggest that out of 30 ABC proteins, only 22 proteins possess either two or one TMD and hence are considered as membrane localized ABC proteins. Notably, our transcriptional dynamics of ABC proteins encoding genes following D. hansenii cells treatment with different salts and drugs concentrations illustrated variable transcriptional response of some of the genes, pointing to their role in salt and drug tolerance. This study first time provides a comprehensive inventory of the ABC proteins of a haploid D. hansenii which will be helpful for exploring their functional relevance.

Proteins of the SNARE (Soluble N-ethylmaleimide-sensitive factor attachment protein receptors) family play a significant role in all vesicular fusion events involved in endocytic and exocytic pathways. These proteins act as molecular machines that assemble into tight four-helix bundle complex, bridging the opposing membranes into close proximity forming membrane fusion. Almost all SNARE proteins share a 53 amino acid coiled-coil domain, which is mostly linked to the transmembrane domain at the C-terminal end. Despite significant variations between SNARE sequences across species, the SNARE mediated membrane fusion is evolutionary conserved in all eukaryotes. It is of interest to compare the functional divergence of SNARE proteins across various eukaryotic groups during evolution. Here, we report an exhaustive phylogeny of the SNARE proteins retrieved from SNARE database including plants, animals, fungi and protists. The Initial phylogeny segregated SNARE protein sequences into five well-supported clades Qa, Qb, Qc, Qbc and R reflective of their positions in the four-helix SNARE complex. Further to improve resolution the Qa, Qb, Qc and R family specific trees were reconstructed, each of these were further segregated into organelle specific clades at first and later diverged into lineage specific subgroups. This revealed that most of the SNARE orthologs are conserved at subcellular locations or at trafficking pathways across various species during eukaryotic evolution. The paralogous expansion in SNARE repertoire was observed at metazoans (animals) and plants independently during eukaryotic evolution. However, results also show that the multi-cellular and saprophytic fungi have limited SNAREs.

PURPOSE: Transient receptor potential cation channel subfamily M member 4 (TRPM4) has been shown to play a key role in detrusor contractility under physiological conditions. In this study we investigated the potential role of TRPM4 in detrusor overactivity following spinal cord transection (SCT) in mice. METHODS: TRPM4 expression and function were evaluated in bladder tissue with or without the mucosa from spinal intact (SI) and SCT female mice (T8-T9 vertebra; 1–28 days post SCT) using PCR, western blots, immunohistochemistry and muscle strip contractility techniques. RESULTS: TRPM4 was expressed in the urothelium (UT) and detrusor smooth muscle (DSM) and was upregulated after SCT. Expression levels peaked 3–7 days post SCT in both the UT and DSM. Pharmacological block of TRPM4 with the antagonist, 9-Phenanthrol (30 µM) greatly reduced spontaneous phasic activity that developed after SCT, regardless of the presence or absence of the mucosa. CONCLUSIONS: Detrusor overactivity following spinal cord injury leads to incontinence and/or renal impairment, and represents a major health problem for which current treatments are not satisfactory. Augmented TRPM4 expression in the bladder after chronic SCT supports the hypothesis that TRPM4 channels play a role in DSM overactivity following SCT. Inhibition of TRPM4 may be beneficial for improving detrusor overactivity in SCI.

Surface localized microbial enolases' binding with human plasminogen has been increasingly proven to have an important role in initial infection cycle of several human pathogens.Likewise, surface localized Mycobacterium tuberculosis (Mtb) enolase also binds to human plasminogen, and this interaction may entail crucial consequences for granuloma stability. The current study is the first attempt to explore the plasminogen interacting residues of enolase from Mtb. Beginning with the structural modeling of Mtbenolase, the binding pose of Mtbenolase and human plasminogen was predicted using protein-protein docking simulations. The binding pose revealed the interface region with interacting residues and molecular interactions. Next, the interacting residues were refined and ranked by using various criteria. Finally, the selected interacting residues were tested experimentally for their involvement in plasminogen binding. The two consecutive lysine residues, Lys-193 and Lys-194, turned out to be active residues for plasminogen binding. These residues when substituted for alanine along with the most active residue Lys-429 i.e., the triple mutant (K193A + K194A + K429A) Mtb enolase, exhibited40% reduction in plasminogen binding. It is worth noting that Mtb enolase lost nearly half of the plasminogen binding activity with only three simultaneous substitutions, without any significant secondary structure perturbation. Further, the sequence comparison between Mtb and human enolase isoforms suggests the possibility of selective targeting of Mtb enolase to obstruct binding of human plasminogen. This article is protected by copyright. All rights reserved

Human carbonic anhydrase IX (CAIX) has evolved as a promising biomarker for cancer prognosis, due to its overexpression in various cancers and restricted expression in normal tissue. However, limited information is available on its biophysical behavior. The unfolding of CAIX in aqueous urea solution was studied using all-atom molecular dynamics simulation approach. The results of this study revealed a stable intermediate state along the unfolding pathway of CAIX. At intermediate concentrations of urea (2.0–4.0 M), the protein displays a native-like structure with a large population of its secondary structure and hydrophobic contacts remaining intact in addition to small confined overall motions. Beyond 4.0 M urea, the unfolding is more gradual and at 8.0 M urea the structure is largely collapsed due to the solvent effect. The hydrophobic contact analysis suggests that the contact in terminal α-helices is separated initially which propagates in the loss of contacts from centrally located β-sheets. The reduction of 60–65% tertiary contacts in 7.0–8.0 M urea suggested the presence of residual structure in unfolded state and is confirmed with structural snap shot. Free energy landscape analysis suggested that unfolding of CAIX exists through the different intermediate states.

α-Isopropylmalate Synthase (α-IPMS) encoded by leuA in Mycobacterium tuberculosis (M.tb) is involved in the leucine biosynthesis pathway and is extremely critical for the synthesis of branched-chain amino acids (leucine, isoleucine and valine). α-IPMS activity is required not only for the proliferation of M.tb but is also indispensable for its survival during the latent phase of infection. It is absent in humans and is widely regarded as one of the validated drug targets against Tuberculosis (TB). Despite its essentiality, any study on designing of potential chemical inhibitors against α-IPMS has not been reported so far. In the present study, in silico identification of putative inhibitors against α-IPMS exploring three chemical databases i.e. NCI, DrugBank and ChEMBL is reported through structurebased drug design and filtering of ligands based on the pharmacophore feature of the actives. In the absence of experimental results of any inhibitor against α-IPMS, a stringent validation of docking results is done by comparing with molecular mechanics/Poisson- Boltzmann surface area (MM/PBSA) calculations by investigating two more proteins for which experimental results are known.