Structural Prediction and Mutational Analysis of Rv3906c Gene of Mycobacterium tuberculosis H37Rv to Determine Its Essentiality in Survival

Elucidation of furanone as ergosterol pathway inhibitor in Cryptococcus neoformans

In the era of antiretroviral therapy, the prevalence of Cryptococcal infection among HIV patients in developed countries has decreased considerably. However, C. neoformans ranks top among the critical priority pathogen that affects a wide range of immunocompromised individuals. The threat of C. neoformans is because of its incredibly multifaceted intracellular survival capabilities. Cell membrane sterols especially ergosterol and enzymes of its biosynthetic pathway are considered fascinating drug targets because of their structural stability. In this study, the ergosterol biosynthetic enzymes were modeled and docked with furanone derivatives. Among the tested ligands Compound 6 has shown a potential interaction with Lanosterol 14 α-demethylase. This best-docked protein-ligand complex was taken further to molecular dynamics simulation. In addition, Compound 6 was synthesized and an in vitro study was conducted to quantify the ergosterol in Compound 6 treated cells. Altogether the computational and in vitro study demonstrates that Compound 6 has anticryptococcal activity by targeting the biosynthetic pathway of ergosterol.Communicated by Ramaswamy H. Sarma.

Structural and functional characterization of mycobacterial PhoH2 and identification of potential inhibitor of its enzymatic activity

Mycobacterium tuberculosis is composed of a cumbersome signaling and protein network which partakes in bacterial survival and augments its pathogenesis. Mycobacterial PhoH2 (Mt-PhoH2) is a signaling element and a predictive phosphate starvation protein that works in an ATP-dependent manner. Here, we elaborated the characterization of Mt-PhoH2 through biophysical, biochemical, and computational methods. In addition to its intrinsic ATPase activity, the biochemical experiments revealed its GTPase activity and both activities are metal ion dependent. Magnesium, manganese, copper, iron, nickel, zinc, cesium, calcium, and lithium were examined for their effect on activity, and the optimum activity was found with 10 mM of Mg2+ ions. The kinetic parameters of 3 µM Mt-PhoH2 were observed as Km 4.873 ± 0.44 µM, Vmax 12.3817 ± 0.084 µM/min/mg, Kcat 0.0075 ± 0.00005 s-1, and Kcat/Km 0.0015 ± 0.000001 µM-1 s-1 with GTP. In the case of GTP as a substrate, a 20% decrease in enzymatic activity and a 50% increase in binding affinity of Mt-PhoH2 were observed. The substrates ADP and GDP inhibit the ATPase and GTPase activity of Mt-PhoH2. CD spectroscopy showed the dominance of alpha helix in the secondary structure of Mt-PhoH2, and this structural pattern was altered upon addition of ATP and GTP. In silico inhibitor screening revealed ML141 and NAV_2729 as two potential inhibitors of the catalytic activity of Mt-PhoH2. Mt-PhoH2 is essential for mycobacterial growth as its knockdown strain showed a decreased growth effect. Overall, the present article emphasizes the factors essential for the proper functioning of Mt-PhoH2 which is a participant in the toxin-antitoxin machinery and may also play an important role in phosphate starvation.

CDH1 gene as biomarker towards breast cancer prediction

Breast cancer is considered to be happened due to genetic aberration. Out of several genes expressed, it is found that cadherin 1, type 1 (CDH1) is responsible in several ways to control the metabolic order in human. Deregulation of the function of protein E-cadherin, expressed from CDH1 plays an important role in lobular breast cancer. In order to understand the root cause of this recent claim, we focus on CDH1 gene: whether the genetic information translated due to any deviation/alteration/modification in its sequence is related to the occurrence of the different types breast cancer. Towards this end, quantitative analysis of different biophysical and bio-chemical properties of CDH1 gene in genomic and proteomic levels from the available genomic (cDNA) sequences of CDH1 gene (obtained from the COSMIC Database for 78 patients, suffering from various types of breast cancer) clearly emphasizes that alternation/modification in the sequence of the CDH1 gene can be detrimental. Furthermore, Random forest, K-nearest neighbour and stochastic gradient descent (SGD) algorithms are applied on the derived dataset to classify the types of breast cancer, and to validate our hypothesis regarding the acute role of CDH1 as potential bio marker for breast cancer. Analysis of the mutated CDH1 gene sequences, and their related parameters using aforesaid machine learning techniques clearly establish that CDH1 gene can take the deterministic role in predicting the chances of occurrences of different types of breast cancer with an accuracy of > 90 % . Such an observation opens a new paradigm in diagnostic approach of breast cancer.Communicated by Ramaswamy H. Sarma.

Structure prediction, docking studies and molecular cloning of novel Pichia kudriavzevii YK46 metalloprotease (MetPr) for improvement of feather waste biodegradation

This study addresses the environmental risks associated with the accumulation of keratin waste from poultry, which is resistant to conventional protein degradation methods. To tackle this issue, microbial keratinases have emerged as promising tools for transforming resilient keratin materials into valuable products. We focus on the Metalloprotease (MetPr) gene isolated from novel Pichia kudriavzevii YK46, sequenced, and deposited in the NCBI GenBank database with the accession number OQ511281. The MetPr gene encodes a protein consisting of 557 amino acids and demonstrates a keratinase activity of 164.04 U/ml. The 3D structure of the protein was validated using Ramachandran's plot, revealing that 93% and 97.26% of the 557 residues were situated within the most favoured region for the MetPr proteins of template Pichia kudriavzevii strain 129 and Pichia kudriavzevii YK46, respectively. Computational analyses were employed to determine the binding affinities between the deduced protein and beta keratin. Molecular docking studies elucidated the optimal binding affinities between the metalloprotease (MetPr) and beta-keratin, yielding values of - 260.75 kcal/mol and - 257.02 kcal/mol for the template strains Pichia kudriavzevii strain 129 and Pichia kudriavzevii YK46, respectively. Subsequent molecular cloning and expression of the MetPr gene in E. coli DH5α led to a significantly higher keratinase activity of 281 ± 12.34 U/ml. These findings provide valuable insights into the potential of the MetPr gene and its encoded protein for keratin waste biotransformation, with implications for addressing environmental concerns related to keratinous waste accumulation.

Structural-functional analysis and molecular characterization of arsenate reductase from Enterobacter cloacae RSC3 for arsenic biotransformation

Enterobacter cloacae RSC3 isolated from an industrial pesticide site transformed arsenate into arsenite. The arsenate is transported by membrane-bound phosphate transporter and transformed to arsenite by arsenate reductase (arsC). E. cloacae RSC3 produced an arsenate reductase enzyme with a maximum activity of 354 U after 72 h of incubation. Arsenate reductase was found to be active and stable at a wide range of temperatures (20 and 45 °C) and pH (5-10), with maximum activity at 35 °C and pH 7.0. The arsenate reductase protein was further characterised molecularly using different bioinformatics tools. The 3D structure of ArsC protein was predicted by homology modelling and validated by the Ramachandran plot with 91.9% residues in the most favoured region. ArsC protein of E. cloacae RSC3 revealed structural homology with ArsC from PDB ID: 1S3C. The gene ontology results also showed that the ArsC protein had a molecular functionality of the arsenate reductase (glutaredoxin) activity and the biological function of cellular response to DNA damage stimulus. Molecular docking analysis of 3D structures using AutoDock vina-1.5.7 server predicted four ligand binding active site residues at Gln70, Asp68, Leu68, and Leu63. Strong ArsC-arsenate ion interaction was observed with binding energy -1.03 kcal/mol, indicating significant arsenate reductase activity and specificity of ArsC protein. On the basis of molecular dynamics simulation analysis, the RMSD and RMSF values revealed the stability of ArsC protein from E. cloacae RSC3.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03730-9.

In silico exploration and molecular dynamics of deleterious SNPs on the human TERF1 protein triggering male infertility

By limiting chromosome erosion and end-to-end fusions, telomere integrity is critical for chromosome stability and cell survival. During mitotic cycles or due to environmental stresses, telomeres become progressively shorter and dysfunctional, thus triggering cellular senescence, genomic instability and cell death. To avoid such consequences, the telomerase action, as well as the Shelterin and CST complexes, assure the telomere's protection. Telomeric repeat binding factor 1 (TERF1), which is one of the primary components of the Shelterin complex, binds directly to the telomere and controls its length and function by regulating the telomerase activity. Several reports about TERF1 gene variations have been associated with different diseases, and some of them have linked these variations to male infertility. Hence, this paper can be advantageous to investigate the association between the missense variants of the TERF1 gene and the susceptibility to male infertility. The stepwise prediction of SNPs pathogenicity followed in this study was based on stability and conservation analysis, post-translational modification, secondary structure, functional interaction prediction, binding energy evaluation and finally molecular dynamic simulation. Prediction matching among the tools revealed that out of 18 SNPs, only four (rs1486407144, rs1259659354, rs1257022048 and rs1320180267) were predicted as the most damaging and highly deleterious SNPs affecting the TERF1 protein and its molecular dynamics when interacting with the TERB1 protein by influencing the function, structural stability, flexibility and compaction of the overall complex. Interestingly, these polymorphisms should be considered during genetic screening so they can be used effectively as genetic biomarkers for male infertility diagnosis.Communicated by Ramaswamy H. Sarma.

Fludarabine, a Potential DNA-Dependent RNA Polymerase Inhibitor, as a Prospective Drug against Monkeypox Virus: A Computational Approach

Monkeypox is a zoonotic contagious disease that has recently re-emerged in different countries worldwide. Due to the lack of an effective treatment that eliminates the virus, there is an urgent need to find effective drugs to stop the spread of the multi-country outbreak. The current study aimed to use computational methods to quickly identify potentially effective drugs against the Monkeypox virus (MPXV). Three MPXV proteins were targeted in this study due to their essential role in viral replication (a DNA-Dependent RNA Polymerase subunit (A6R)), a protein involved in cell entry (D8L), and a protein catalyzing the envelopment of intracellular mature virus particles (F13L). We virtually screened a library of 1615 FDA-approved compounds, utilizing different in-silico approaches including computational modeling, molecular docking, molecular dynamic (MD) simulation, and MM-GBSA. The compound Fludarabine was found to have the best docking score (−7.53 kcal/mol) in relation to the MPXV A6R protein. Additionally, Fludarabine showed in-silico activity on the D8L and F13L proteins. During the whole period of the 100 ns MD simulation, the complex of A6R and Fludarabine exhibited the best stability. This stability was reflected in a good score of MM-GBSA, with an average value of −44.62 kcal/mole in a range between −53.26 and −35.49 and a low value of standard deviation (3.76). Furthermore, Fludarabine blocked efficiently the Asn175 residue which has an important role in the attachment of the virus to a host cell. The results of this study recommend more in vitro studies on this compound, as a starting point to develop a novel treatment against MPXV.

Mechanistic Insight into the Enzymatic Inhibition of β-Amyrin against Mycobacterial Rv1636: In Silico and In Vitro Approaches

Simple Summary

Rv1636 is a mycobacterial universal stress protein whose expression level increases in different type of stress conditions. This protein promotes the growth of Mycobacterium tuberculosis in the host derived stress conditions generated during infection. Therefore in this manuscipt, we are trying to target Rv1636 using natural inhibitor. Targeting essential Mycobacterial protein using natural prodect was hypothesized to generate a molecule with low toxic effects and high inhibitory activity. It was found that Rv1636 contains ATPase activity and its ATPase activity gets disturbed by addition of β-Amyrin in the reaction. β-Amyrin was forund to interfere with the ATP binding site of Rv1636 which was confirmed by molecular docking anad dynamic studies. In addition to the ATPase activity, Rv1636 was also contain the cAMP binding capacity and also involved in balancing the cAMP levels inside cells. So, targeting Rv1636 using β-Amyrin disrupts its ATPase activity and cAMP regulatory activity and these conditions might make Mycobacterium tuberculosis more susceptible to the host derived stress conditions.

Abstract

Mycobacterium tuberculosis has seen tremendous success as it has developed defenses to reside in host alveoli despite various host-related stress circumstances. Rv1636 is a universal stress protein contributing to mycobacterial survival in different host-derived stress conditions. Both ATP and cAMP can be bound with the Rv1636, and their binding actions are independent of one another. β-Amyrin, a triterpenoid compound, is abundant in medicinal plants and has many pharmacological properties and broad therapeutic potential. The current study uses biochemical, biophysical, and computational methods to define the binding of Rv1636 with β-Amyrin. A substantial interaction between β-Amyrin and Rv1636 was discovered by molecular docking studies, which helped decipher the critical residues involved in the binding process. VAL60 is a crucial residue found in the complexes of both Rv1636_β-Amyrin and Rv1636-ATP. Additionally, the Rv1636_β-Amyrin complex was shown to be stable by molecular dynamics simulation studies (MD), with minimal changes observed during the simulation. In silico observations were further complemented by in vitro assays. Successful cloning, expression, and purification of Rv1636 were accomplished using Ni-NTA affinity chromatography. The results of the ATPase activity assay indicated that Rv1636’s ATPase activity was inhibited in the presence of various β-Amyrin concentrations. Additionally, circular dichroism spectroscopy (CD) was used to examine modifications to Rv1636 secondary structure upon binding of β-Amyrin. Finally, isothermal titration calorimetry (ITC) advocated spontaneous binding of β-Amyrin with Rv1636 elucidating the thermodynamics of the Rv1636_β-Amyrin complex. Thus, the study establishes that β-Amyrin binds to Rv1636 with a significant affinity forming a stable complex and inhibiting its ATPase activity. The present study suggests that β-Amyrin might affect the functioning of Rv1636, which makes the bacterium vulnerable to different stress conditions.

Computational Analysis of Missense Variants in the Human Transmembrane Protease Serine 2 (TMPRSS2) and SARS-CoV-2

The human transmembrane protease serine 2 (TMPRSS2) protein plays an important role in prostate cancer progression. It also facilitates viral entry into target cells by proteolytically cleaving and activating the S protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In the current study, we used different available tools like SIFT, PolyPhen2.0, PROVEAN, SNAP2, PMut, MutPred2, I-Mutant Suite, MUpro, iStable, ConSurf, ModPred, SwissModel, PROCHECK, Verify3D, and TM-align to identify the most deleterious variants and to explore possible effects on the TMPRSS2 stability, structure, and function. The six missense variants tested were evaluated to have deleterious effects on the protein by SIFT, PolyPhen2.0, PROVEAN, SNAP2, and PMut. Additionally, V160M, G181R, R240C, P335L, G432A, and D435Y variants showed a decrease in stability by at least 2 servers; G181R, G432A, and D435Y are highly conserved and identified posttranslational modifications sites (PTMs) for proteolytic cleavage and ADP-ribosylation using ConSurf and ModPred servers. The 3D structure of TMPRSS2 native and mutants was generated using 7 meq as a template from the SwissModeller group, refined by ModRefiner, and validated using the Ramachandran plot. Hence, this paper can be advantageous to understand the association between these missense variants rs12329760, rs781089181, rs762108701, rs1185182900, rs570454392, and rs867186402 and susceptibility to SARS-CoV-2.

Characterization of Mercury-Resistant Rhizobacteria for Plant Growth Promotion: An In Vitro and In Silico Approach

In this study, a total 30 rhizobacterial isolates were screened out based on resistance against different concentrations of mercuric chloride (HgCl₂), growth on nitrogen-free mannitol (NFM) and production of indole-3-acetic acid (IAA). The biochemical and plant growth promoting characterization of selected isolates was performed by different biochemical tests. Out of 30, six isolates, UM-3, AZ-5, UM-7, UM-11, UM-26, and UM-28 showed resistance at 30 µg/ml HgCl₂, pronounced growth on NFM and high production of IAA as 18.6, 16.7, 16, 18.7, 14, and 16 µg/ml, respectively (P < 0.05). The 16S rDNA ribotyping and phylogenetic analysis of selected bacterial isolates were performed and characterized as Exiguobacterium sp. UM-3 (KJ736011), Bacillus thuringiensis AZ-5 (KJ675627), Bacillus subtilis UM-7 (KJ736013), Enterobacter cloacae UM-11 (KJ736014), Pseudomonas aeruginosa UM-26 (KJ736016), P. aeruginosa UM-28 (KJ736017) and Bacillus pumilus UM-16 (KJ736015) used as negative control. B. thuringiensis AZ-5 showed high resistance against 30 µg/ml of HgCl₂ due to the presence of merB gene. The structural determination of MerB protein was carried out using bioinformatics tools, i.e., Protparam, Pfam, InterProScan, STRING, Jpred4, PSIPRED, I-TASSER, COACH server and ERRAT. These tools predicted the structural based functional homology of MerB protein (organomercuric lyase) in association with MerA (mercuric reductase) in bacterial Hg-detoxification system.

Revealing new therapeutic opportunities in hypertension through network-driven integrative genetic analysis and drug target prediction approach

Epidemiological studies have established that untreated hypertension (HTN) is a major independent risk factor for developing cardiovascular diseases (CVD), stroke, renal failure, and other conditions. Several important studies have been published to prevent and manage HTN; however, antihypertensive agents' optimal choice remains controversial. Therefore, the present study is undertaken to update our knowledge in the primary treatment of HTN, specifically in the setting of other three important diseases. MicroRNAs (miRNAs) are remarkably stable short endogenous conserved non-coding RNAs that bind to the mRNA at its (3' UTR) to regulate its gene expression by causing translational repression or mRNA degradation. Through their coordinated activities on different pathways and networks, individual miRNAs control normal and pathological cellular processes. Therefore, to identify the critical miRNA-mRNA-TF interactions, we performed systematic bioinformatics analysis. We have also employed the molecular modelling and docking approach to identify the therapeutic target that delivers novel empathies into Food and Drug Administration approved and herbal drug response physiology. Gene Expression Omnibus (GEO) was employed to identify the differentially expressed genes (DEGs) and hub genes- KNG1, HLA-DPB1, CXCL8, IL1B, and BCL2. The HTN associated feed-forward loop (FFL) network included miR-9-5p, KNG1 and AR. We employed high throughput screening to get the best interacting compounds, telmisartan and limonin, that provided a significant docking score (-13.3 and -12.0 kcal/mol) and a potential protective effect that may help to combat the impact of HTN. The present study provides novel insight into HTN etiology through the identification of mRNAs and miRNAs and associated pathways.

Computational-approach understanding the structure-function prophecy of Fibrinolytic Protease RFEA1 from Bacillus cereus RSA1

Microbial fibrinolytic proteases are therapeutic enzymes responsible to ameliorate thrombosis, a fatal cardiac-disorder which effectuates due to excessive fibrin accumulation in blood vessels. Inadequacies such as low fibrin specificity, lethal after-effects and short life-span of available fibrinolytic enzymes stimulates an intensive hunt for novel, efficient and safe substitutes. Therefore, we herewith suggest a novel and potent fibrinolytic enzyme RFEA1 from Bacillus cereus RSA1 (MK288105). Although, attributes such as in-vitro purification, characterization and thrombolytic potential of RFEA1 were successfully accomplished in our previous study. However, it is known that structure-function traits and mode of action significantly aid to commercialization of an enzyme. Also, predicting structural model of a protein from its amino acid sequence is challenging in computational biology owing to intricacy of energy functions and inspection of vast conformational space. Our present study thus reports In-silico structural-functional analysis of RFEA1. Sequence based modelling approaches such as-Iterative threading ASSEmbly Refinement (I-TASSER), SWISS-MODEL, RaptorX and Protein Homology/analogY Recognition Engine V 2.0 (Phyre2) were employed to model three-dimensional structure of RFEA1 and the modelled RFEA1 was validated by structural analysis and verification server (SAVES v6.0). The modelled crystal structure revealed the presence of high affinity Ca1 binding site, associated with hydrogen bonds at Asp147, Leu181, Ile185 and Val187residues. RFEA1 is structurally analogous to Subtilisin E from Bacillus subtilis 168. Molecular docking analysis using PATCH DOCK and FIRE DOCK servers was performed to understand the interaction of RFEA1 with substrate fibrin. Strong RFEA1-fibrin interaction was observed with high binding affinity (-21.36 kcal/mol), indicating significant fibrinolytic activity and specificity of enzyme RFEA1. Overall, the computational research suggests that RFEA1 is a subtilisin-like serine endopeptidase with proteolytic potential, involved in thrombus hydrolysis.

Computational modeling and bioinformatic analyses of functional mutations in drug target genes in Mycobacterium tuberculosis

•

MycoTRAP-DB, a database of mutations and their impact on normal functionality of protein in M.tb genes.

•

Several secondary mutations were identified with significant impact on protein structure and function.

•

Comprehensive information gives insight for screening of suspected hotspots in advance to combat drug resistant TB.

Tuberculosis (TB) continues to be the leading cause of deaths due to its persistent drug resistance and the consequent ineffectiveness of anti-TB treatment. Recent years witnessed huge amount of sequencing data, revealing mutations responsible for drug resistance. However, the lack of an up-to-date repository remains a barrier towards utilization of these data and identifying major mutations-associated with resistance. Amongst all mutations, non-synonymous mutations alter the amino acid sequence of a protein and have a much greater effect on pathogenicity. Hence, this type of gene mutation is of prime interest of the present study. The purpose of this study is to develop an updated database comprising almost all reported substitutions within the Mycobacterium tuberculosis (M.tb) drug target genes rpoB, inhA, katG, pncA, gyrA and gyrB. Various bioinformatics prediction tools were used to assess the structural and biophysical impacts of the resistance causing non-synonymous single nucleotide polymorphisms (nsSNPs) at the molecular level. This was followed by evaluating the impact of these mutations on binding affinity of the drugs to target proteins. We have developed a comprehensive online resource named MycoTRAP-DB (Mycobacterium tuberculosis Resistance Associated Polymorphisms Database) that connects mutations in genes with their structural, functional and pathogenic implications on protein. This database is accessible at http://139.59.12.92. This integrated platform would enable comprehensive analysis and prioritization of SNPs for the development of improved diagnostics and antimycobacterial medications. Moreover, our study puts forward secondary mutations that can be important for prognostic assessments of drug-resistance mechanism and actionable anti-TB drugs.

In Silico Selection and Evaluation of Pugnins with Antibacterial and Anticancer Activity Using Skin Transcriptome of Treefrog (Boana pugnax)

In order to combat bacterial and cancer resistance, we identified peptides (pugnins) with dual antibacterial l-anticancer activity from the Boana pugnax (B. pugnax) skin transcriptome through in silico analysis. Pugnins A and B were selected owing to their high similarity to the DS4.3 peptide, which served as a template for their alignment to the B. pugnax transcriptome, as well as their function as part of a voltage-dependent potassium channel protein. The secondary peptide structure stability in aqueous medium was evaluated as well, and after interaction with the Escherichia coli (E. coli) membrane model using molecular dynamics. These pugnins were synthesized via solid-phase synthesis strategy and verified by Reverse phase high-performance liquid chromatography (RP-HPLC) and mass spectrometry. Subsequently, their alpha-helix structure was determined by circular dichroism, after which antibacterial tests were then performed to evaluate their antimicrobial activity. Cytotoxicity tests against cancer cells also showed selectivity of pugnin A toward breast cancer (MFC7) cells, and pugnin B toward prostate cancer (PC3) cells. Alternatively, flow cytometry revealed necrotic cell damage with a major cytotoxic effect on human keratinocytes (HaCaT) control cells. Therefore, the pugnins found in the transcriptome of B. pugnax present dual antibacterial-anticancer activity with reduced selectivity to normal eukaryotic cells.

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.

Computational evaluation of anticipated PE_PGRS39 protein involvement in host-pathogen interplay and its integration into vaccine development

Mycobacterium tuberculosis causes more than 1 million deaths every year, which is higher than any other bacterial pathogen. Its success depends on its interaction with the host and its ability to regulate the host's immune system for its own survival. Mycobacterium tuberculosis H37Rv (Mtb) proteome consists of unique PE_PGRS family proteins, which present a significant role in bacterial pathogenesis over the past years. Earlier evidence suggests that some PE_PGRS proteins display fibronectin-binding activity. In this manuscript, computational characterization of the PE_PGRS39 protein has indicated something peculiar about this protein. Investigation showed that PE_PGRS39 is an extracellular protein that, instead of acting as fibronectin-binding protein, might mimic fibronectin which binds to alpha-5 beta-1 (α5β1) integrin. PE_PGRS39 protein additionally turned into proven pieces of evidence to have motifs such as DXXG and GGXGXD and PXXP that bind with guanosine triphosphate (GTP), calcium, and host Src homology 3 (SH3) domains, respectively, in conjunction with RGD-integrin binding. These interactions designate the direct role of PE_PGRS39 in bacterial pathogenesis via cell adhesion and signaling. Additionally, the analysis showed that PE_PGRS39 is an antigenic protein and epitope prediction provided functional regions of the protein that trigger a cellular immune response facilitated by T or B cells. Further, an experimental analysis could also open up new avenues for developing novel drugs by targeting signaling motifs or novel vaccines using functional epitopes that could evoke an immune response in the host.

In silico studies of selected multi-drug targeting against 3CLpro and nsp12 RNA-dependent RNA-polymerase proteins of SARS-CoV-2 and SARS-CoV

An outbreak of a cluster of viral pneumonia cases, subsequently identified as coronavirus disease 2019 (COVID-19), due to a novel SARS-CoV-2 necessitates an urgent need for a vaccine to prevent infection or an approved medication for a cure. In our in silico molecular docking study, a total of 173 compounds, including FDA-approved antiviral drugs, with good ADME descriptors, and some other nucleotide analogues were screened. The results show that these compounds demonstrate strong binding affinity for the residues at the active sites of RNA-dependent RNA-polymerase (RdRp) modelled structures and Chymotrypsin-like cysteine protease (3CLpro) of the HCoV proteins. Free energies (ΔG's) of binding for SARS-CoV-2 and SARS-CoV RdRp range from - 5.4 to - 8.8 kcal/mol and - 4.9 to - 8.7 kcal/mol, respectively. Also, SARS-CoV-2 and SARS-CoV 3CLpro gave ΔG values ranging from - 5.1 to - 8.4 kcal/mol and - 5.5 to - 8.6 kcal/mol, respectively. Interesting results are obtained for ivermectin, an antiparasitic agent with broad spectrum activity, which gave the highest binding energy value (- 8.8 kcal/mol) against the 3CLpro of SARS-CoV-2 and RdRps of both SARS-CoV and SARS-CoV-2. The reason for such high binding energy values is probably due to the presence of hydroxy, methoxy and sugar moieties in its structure. The stability of the protein-ligand complexes of polymerase inhibitors considered in this investigation, such as Sofosbuvir, Remdesivir, Tenofovir, Ribavirin, Galidesivir, 5c3, 5h1 and 7a1, show strong to moderate hydrogen bonding and hydrophobic interactions (π-π stacked, π-π T-shaped, π-sigma and π-alkyl). The stability provided from such interactions translate into greater antiviral activity or inhibitory effect of the ligands. Assessment of the average free energies of binding of the FDA approved drugs are highly comparable for conformers of a particular inhibitor, indicating similar modes of binding within the pockets.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13721-021-00299-2.

Domain-wise differentiation of Mycobacterium tuberculosis H37 Rv hypothetical proteins: A roadmap to discover bacterial survival potentials

Proteomic information revealed approximately 3,923 proteins in Mycobacterium tuberculosis H₃₇ Rv genome of which around ∼25% of proteins are hypothetical proteins (HPs). The present work comprises computational approaches to identify and characterize the HPs of M. tuberculosis that symbolize the putative target for rationale development of a drug or antituberculosis strategy. Proteins were primarily classified based on motif and domain information, which were further analyzed for the presence of virulence factors (VFs), determination of localization, and signal peptide/enzymatic cleavage sites. 863 HPs were found, and 599 HPs were finalized based on motifs, that is, GTP (525), Trx (47), SAM (14), PE-PGRS (5), and CBD (8). 80 HPs contain virulence factor (VF), 24 HPs localized in membrane region, and 4 HPs contain signal peptide/enzymatic cleavage sites. The overall parametric study finalizes four HPs Rv0679c, Rv0906, Rv3627c, and Rv3811 that also comprise GTPase domain. Structure prediction, structure-based function prediction, molecular docking and mutation analysis of selected proteins were done. Docking studies revealed that GTP and GTPase inhibitor (mac0182344) were docked with all four proteins with high affinities. In silico point mutation studies showed that substitution of aspartate with glycine within a GTPase motif showed the largest decrease in stability and pH differentiation also affects protein's stability. This analysis thus fixes a roadmap in the direction of finding potential target of this bacterium for drug development and enlightens the efficacy of GTP as a major regulator of Mycobacterial cellular pathways.

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.

Essential biochemical, biophysical and computational inputs on efficient functioning of Mycobacterium tuberculosis H37Rv FtsY

Mycobacterium tuberculosis (M. tuberculosis H₃₇Rv) utilizes the signal recognition particle pathway (SRP pathway) system for secretion of various proteins from ribosomes to the extracellular surface which plays an important role in the machinery running inside the bacterium. This system comprises of three major components FtsY, FfH and 4.5S rRNA. This manuscript highlights essential factors responsible for the optimized enzymatic activity of FtsY. Kinetic parameters include V_max and K_m for the hydrolysis of GTP by ftsY which were 20.25±5.16 μM/min/mg and 39.95±7.7 μM respectively. k_cat and catalytic efficiency of the reaction were 0.012±0.003 s^-1 and 0.00047±0.0001 μM/s^-1 respectively. These values were affected upon changing the standard conditions. Cations (Mg²⁺ and Mn²⁺) play important role in FtsY enzymatic activity as increasing Mg²⁺ decrease the activity. Mn²⁺on the other hand is required at higher concentration around 60 mM for carrying optimum GTPase activity. FtsY is hydrolyzing ATP and GDP as well and GDP acts as an inhibitor of the reaction. MD simulation shows effective binding and stabilization of the FtsY complexed structure with GTP, GDP and ATP. Mutational analysis was done at two important residues of GTP binding motif of FtsY, namely, GXXXXGK (K236) and DXXG (D367) and showed that these mutations significantly decrease FtsY GTPase activity. FtsY is comprised of α helices, but this structural pattern was shown to change with increasing concentrations of GTP and ATP which symbolize that these ligands cause significant conformational change by variating the secondary structure to transduce signals required by downstream effectors. This binding favors the functional stabilization of FtsY by destabilization of α-helix integrity. Revealing the hidden aspects of the functioning of FtsY might be an essential part for the understanding of the SRP pathway which is one of the important contributors of M. tuberculosis virulence.

Structure-Based Virtual Screening to Discover Potential Lead Molecules for the SARS-CoV-2 Main Protease

The COVID-19 disease is caused by a new strain of the coronavirus family (SARS-CoV-2), and it has affected at present millions of people all over the world. The indispensable role of the main protease (Mpro) in viral replication and gene expression makes this enzyme an attractive drug target. Therefore, inhibition of SARS-CoV-2 Mpro as a proposition to halt virus ingression is being pursued by scientists globally. Here we carried out a study with two objectives: the first being to perform comparative protein sequence and 3D structural analysis to understand the effect of 12 point mutations on the active site. Among these, two mutations, viz., Ser46 and Phe134, were found to cause a significant change at the active sites of SARS-CoV-2. The Ser46 mutation present at the entrance of the S5 subpocket of SARS-CoV-2 increases the contribution of other two hydrophilic residues, while the Phe134 mutation, present in the catalytic cysteine loop, can cause an increase in catalytic efficiency of Mpro by facilitating fast proton transfer from the Cys145 to His41 residue. It was observed that active site remained conserved among Mpro of both SARS-CoVs, except at the entrance of the S5 subpocket, suggesting sustenance of substrate specificity. The second objective was to screen the inhibitory effects of three different data sets (natural products, coronaviruses main protease inhibitors, and FDA-approved drugs) using a structure-based virtual screening approach. A total of 73 hits had a combo score >2.0. Eight different structural scaffold classes were identified, such as one/two tetrahydropyran ring(s), dipeptide/tripeptide/oligopeptide, large (approximately 20 atoms) cyclic peptide, and miscellaneous. The screened hits showed key interactions with subpockets of the active site. Further, molecular dynamics studies of selected screened compounds confirmed their perfect fitting into the subpockets of the active site. This study suggests promising structures that can fit into the SARS-CoV-2 Mpro active site and also offers direction for further lead optimization and rational drug design.

Prediction and Characterization of Cationic Arginine-Rich Plant Antimicrobial Peptide SM-985 From Teosinte (Zea mays ssp. mexicana)

Antimicrobial peptides (AMPs) are effective against different plant pathogens and newly considered as part of plant defense systems. From prokaryotes to eukaryotes, AMPs can exist in all forms of life. SM-985 is a cationic AMP (CAMP) isolated from the cDNA library of Mexican teosinte (Zea mays ssp. mexicana). A computational prediction server running with different algorithms was used to screen the teosinte cDNA library for AMPs, and the SM-985 peptide was predicted as an AMP with high probability prediction values. SM-985 is an arginine-rich peptide and composed of 21 amino acids (MW: 2671.06 Da). The physicochemical properties of SM-985 are very promising as an AMP, including the net charge (+8), hydrophobicity ratio of 23%, Boman index of 5.19 kcal/mol, and isoelectric point of 12.95. The SM-985 peptide has amphipathic α-helix conformations. The antimicrobial activity of SM-985 was confirmed against six bacterial plant pathogens, and the MIC of SM-985 against Gram-positive indicators was 8 μM, while the MIC of SM-985 against Gram-negative indicators was 4 μM. The SM-985 interacting with the bacterial membrane and this interaction were examined by treatment of the bacterial indicators with FITC-SM-985 peptide, which showed a high binding affinity of SM-985 to the bacterial membrane (whether Gram-positive or Gram-negative). Scanning electron microscopy (SEM) and transmission electron microscopy (TEM) images of the treated bacteria with SM-985 demonstrated cell membrane damage and cell lysis. In vivo antimicrobial activity was examined, and SM-985 prevented leaf spot disease infection caused by Pst DC3000 on Solanum lycopersicum. Moreover, SM-985 showed sensitivity to calcium chloride salt, which is a common feature of CAMPs.

Increases in Hydrophilicity and Charge on the Polar Face of Alyteserin 1c Helix Change its Selectivity towards Gram-Positive Bacteria

Recently, resistance of pathogens towards conventional antibiotics has increased, representing a threat to public health globally. As part of the fight against this, studies on alternative antibiotics such as antimicrobial peptides have been performed, and it has been shown that their sequence and structure are closely related to their antimicrobial activity. Against this background, we here evaluated the antibacterial activity of two peptides developed by solid-phase synthesis, Alyteserin 1c (WT) and its mutant derivative (ΔM), which shows increased net charge and reduced hydrophobicity. These structural characteristics were modified as a result of amino acid substitutions on the polar face of the WT helix. The minimum inhibitory concentration (MIC) of both peptides was obtained in Gram-positive and Gram-negative bacteria. The results showed that the rational substitutions of the amino acids increased the activity in Gram-positive bacteria, especially against Staphylococcus aureus, for which the MIC was one-third of that for the WT analog. In contrast to the case for Gram-positive bacteria, these substitutions decreased activity against Gram-negative bacteria, especially in Escherichia coli, for which the MIC was eight-fold higher than that exhibited by the WT peptide. To understand this, models of the peptide behavior upon interacting with membranes of E. coli and S. aureus created using molecular dynamics were studied and it was determined that the helical stability of the peptide is indispensable for antimicrobial activity. The hydrogen bonds between the His20 of the peptides and the phospholipids of the membranes should modulate the selectivity associated with structural stability at the carboxy-terminal region of the peptides.

Molecular docking studies of chloroquine and its derivatives against P23pro-zbd domain of chikungunya virus: Implication in designing of novel therapeutic strategies

The arthropod-transmitted chikungunya virus has emerged as an epidemic menace that causes debilitating polyarthritis. With this life-threatening impact on humans, the possible treatment requires to cure the viral infectivity. But, devoid of any vaccine against the chikungunya virus (CHIKV), there is a need to develop a novel chemotherapeutic strategy to treat this noxious infection. CHIKV carries highly compact P23^pro-zbd structure that possesses potential RNA-binding surface domains which extremely influences the use of RNA template during genome replication at the time of infection and pathogenesis. Therefore, computational approaches were used to explore the novel small molecule inhibitors targeting P23^pro-zbd domain. The tertiary structure was modeled and optimized using in silico approaches. The results obtained from PROCHECK (93.1% residues in favored regions), ERRAT (87.480 overall model quality) and ProSA (Z-score: -11.72) revealed the reliability of the proposed model. Interestingly, a previously reported inhibitor, chloroquine possesses good binding affinities with the target domain. In-depth analysis revealed that chloroquine derivatives such as didesethyl chloroquine hydroxyacetamide, cletoquine, hydroxychloroquine exhibited a better binding affinity. Notably, MD simulation analysis exhibited that Thr1312, Ala1355, Ala1356, Asn1357, Asp1364, Val1366, Cys1367, Ala1401, Gly1403, Ser1443, Tyr1444, Gly1445, Asn1459, and Thr1463 residues are the key amino acid responsible for stable ligand-protein interaction. The results obtained from this study provide new insights and advances the understanding to develop a new approach to consider effective and novel drug against chikungunya. However, a detailed in vivo study is required to explore its drug likeliness against this life-threatening disease.