Insight into the drug resistance whole genome of Mycobacterium tuberculosis isolates from Khyber Pakhtunkhwa, Pakistan.

Whole genome sequencing (WGS) is one of the most reliable methods for detection of drug resistance, genetic diversity in other virulence factor and also evolutionary dynamics of Mycobacterium tuberculosis complex (MTBC). First-line anti-tuberculosis drugs are the major weapons against Mycobacterium tuberculosis (MTB). However, the emergence of drug resistance remained a major obstacle towards global tuberculosis (TB) control program 2030, especially in high burden countries including Pakistan. To overcome the resistance and design potent drugs, genomic variations in drugs targets as well as in the virulence and evolutionary factors might be useful for better understanding and designing potential inhibitors. Here we aimed to find genomic variations in the first-line drugs targets, along with other virulence and evolutionary factors among the circulating isolates in Khyber Pakhtunkhwa, Pakistan. Samples were collected and drug susceptibility testing (DST) was performed as per WHO standard. The resistance samples were subjected to WGS. Among the five whole genome sequences, three samples (NCBI BioProject Accession: PRJNA629298, PRJNA629388) harbored 1997, 1162, and 2053 mutations. Some novel mutations have been detected in drugs targets. Similarly, numerous novel variants have also been detected in virulency and evolutionary factors, PE, PPE, and secretory system of MTB isolates. Exploring the genomic variations among the circulating isolates in geographical specific locations might be useful for future drug designing. To the best of our knowledge, this is the first study that provides useful data regarding the insight genomic variations in virulency, evolutionary factors including ESX and PE/PPE as well as drug targets, for better understanding and management of TB in a WHO declared high burden country.

Structures of SARS-CoV-2 RNA-Binding Proteins and Therapeutic Targets.

BACKGROUND: The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) epidemic has resulted in thousands of infections and deaths worldwide. Several therapies are currently undergoing clinical trials for the treatment of SARS-CoV-2 infection. However, the development of new drugs and the repositioning of existing drugs can only be achieved after the identification of potential therapeutic targets within structures, as this strategy provides the most precise solution for developing treatments for sudden epidemic infectious diseases. SUMMARY: In the current investigation, crystal and cryo-electron microscopy structures encoded by the SARS-CoV-2 genome were systematically examined for the identification of potential drug targets. These structures include nonstructural proteins (Nsp-9; Nsp-12; and Nsp-15), nucleocapsid (N) proteins, and the main protease (Mpro). Key Message: The structural information reveals the presence of many potential alternative therapeutic targets, primarily involved in interaction between N protein and Nsp3, forming replication-transcription complexes (RTCs) which might be a potential drug target for effective control of current SARS-CoV-2 pandemic. RTCs consist of 16 nonstructural proteins (Nsp1-16) that play the most essential role in the synthesis of viral RNA. Targeting the physical linkage between the envelope and single-stranded positive RNA, a process facilitated by matrix proteins may provide a good alternative strategy. Our current study provides useful information for the development of new lead compounds against SARS-CoV-2 infections.

Marine natural compounds as potents inhibitors against the main protease of SARS-CoV-2-a molecular dynamic study.

Sever acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a single-stranded RNA (ssRNA) virus, responsible for severe acute respiratory disease (COVID-19). A large number of natural compounds are under trial for screening compounds, possessing potential inhibitory effect against the viral infection. Keeping in view the importance of marine compounds in antiviral activity, we investigated the potency of some marine natural products to target SARS-CoV-2 main protease (Mpro) (PDB ID 6MO3). The crystallographic structure of Mpro in an apo form was retrieved from Protein Data Bank and marine compounds from PubChem. These structures were prepared for docking and the complex with good docking score was subjected to molecular dynamic (MD) simulations for a period of 100 ns. To measure the stability, flexibility, and average distance between the target and compounds, root mean square deviations (RMSD), root mean square fluctuation (RMSF), and the distance matrix were calculated. Among five marine compounds, C-1 (PubChem CID 11170714) exhibited good activity, interacting with the active site and surrounding residues, forming many hydrogen and hydrophobic interactions. The C-1 also attained a stable dynamic behavior, and the average distance between compound and target remains constant. In conclusion, marine natural compounds may be used as a potential inhibitor against SARS-CoV-2 for better management of COVID-19.

Molecular dynamics simulation and binding free energy calculations of microcin J25 binding to the FhuA receptor.

Microcins are low-molecular weight, highly stable, ribosomally produced, bacterial inhibitory molecules involved in competitive research. Microcin J25 consists of 21 amino acids and has a lasso-like structure. The first step is bacteria binding to the ion receptor of FhuA on the bacterial surface. In this study, molecular dynamics simulation was implemented to study the binding mechanism of MJ25 and three mutants, to find the changing of individual amino acids in the ß-hair region of MJ25. The binding-free energy calculation was subjected to MJ25 and FhuA. In addition, computational mutation analysis was revealed for association between MJ25 and FhuA. However, we found that the mutating 12th amino acid of the ß-hair region into histidine is extremely important for the binding of MJ25 to FhuA. In addition, the number of hydrogen bonds is essential for binding of MJ25 to FhuA. The overall results show that the key guiding significance is improving the sterilizing efficiency of MJ25 and the drug design of MJ25.Communicated by Ramaswamy H. Sarma.

Mechanistic analysis of A46V, H57Y, and D129N in pyrazinamidase associated with pyrazinamide resistance.

Pyrazinamide (PZA) is a component of first-line drugs, active against latent Mycobacterium tuberculosis (MTB) isolates. The prodrug is activated into the active form, pyrazinoic acid (POA) via pncA gene-encoded pyrazinamidase (PZase). Mutations in pncA have been reported, most commonly responsible for PZA-resistance in more than 70% of the resistant cases. In our previous study, we detected many mutations in PZase among PZA-resistance MTB isolates including A46V, H71Y, and D129N. The current study was aimed to investigate the molecular mechanism of PZA-resistance behind mutants (MTs) A46V, H71Y, and D129N in comparison with the wild type (WT) through molecular dynamic (MD) simulation. MTB positive samples were subjected to PZA drug susceptibility testing (DST) against critical concentration (100ug/ml). The resistant samples were subjected to pncA sequencing. Thirty-six various mutations have been observed in the coding region of pncA of PZA-resistant isolates (GenBank accession No. MH461111) including A46V, H71Y, and D129N. The post-simulation analysis revealed a significant variation in MTs structural dynamics as compared to the WT. Root means square deviations (RMSD) and Root means square fluctuation (RMSF) has been found in variation between WT and MTs. Folding effect and pocket volume were altered in MTs when compared with WT. Geometric matching supports the effect of mutation A46V, H71Y, and D129N on PZase structure that may have an insight effect on PZase dynamics, making them vulnerable to convert pro-PZA into active form, POA. In conclusion, the current analyses will provide useful information behind PZA-resistance for better management of drug-resistant TB.

Pyrazinamide resistance of novel mutations in pncA and their dynamic behavior.

Pyrazinamide (PZA) is one of the essential anti-mycobacterium drugs, active against non-replicating Mycobacterium tuberculosis (MTB) isolates. PZA is converted into its active state, called pyrazinoic acid (POA), by action of pncA encoding pyrazinamidase (PZase). In the majority of PZA-resistance isolates, pncA harbored mutations in the coding region. In our recent report, we detected a number of novel variants in PZA-resistance (PZAR) MTB isolates, whose resistance mechanisms were yet to be determined. Here we performed several analyses to unveil the PZAR mechanism of R123P, T76P, G150A, and H71R mutants (MTs) through molecular dynamics (MD) simulations. In brief, culture positive MTB isolates were subjected to PZA susceptibility tests using the WHO recommended concentration of PZA (100 µg ml-1). The PZAR samples were screened for mutations in pncA along sensitive isolates through polymerase chain reactions and sequencing. A large number of variants (GeneBank accession no. MH461111), including R123P, T76P, G150A, and H71R, have been spotted in more than 70% of isolates. However, the mechanism of PZAR for mutants (MTs) R123P, T76P, G150A, and H71R was unknown. For the MTs and native PZase structures (WT), thermodynamic properties were compared using molecular dynamics simulations for 100 ns. The MTs structural activity was compared to the WT. Folding effect and pocket volume variations have been detected when comparing between WT and MTs. Geometric matching further confirmed the effect of R123P, T76P, G150A, and H71R mutations on PZase dynamics, making them vulnerable for activating the pro-drug into POA. This study offers a better understanding for management of PZAR TB. The results may be used as alternative diagnostic tools to infer PZA resistance at a structural dynamics level.

Combining in silico and in vitro approaches to identification of potent inhibitor against phospholipase A2 (PLA2).

Phospholipase A2 (PLA2) is the main constituent of snake venom. PLA2 enzymes catalyze the Ca2+ dependent hydrolysis of 2-acyl ester bonds of 3-sn-phospholipids, releasing fatty acids and lysophospholipids. Inside the body of the victim, PLA2 from snake venom induces either direct or indirect pathophysiological effects, including anticoagulant, inflammatory, neurotoxic, cardiotoxic, edematogenic, and myotoxic activities. Therefore, there is a need to find the potential inhibitors against PLA2 responsible for snakebite. In this study, we employed in silico and in vitro methods to identify the potential inhibitor against PLA2. Virtual screening and molecular docking studies were performed to find potent inhibitor against PLA2 using Traditional Chinese Medicine Database (TCM). Based on these studies, Scutellarin (TCM3290) was selected and calculated by density functional theory calculation at B3LYP/6-31G**++ level to explore the stereo-electronic features of the molecule. Further, molecular docking and DFT of Minocycline was carried out. Quantum polarized ligand docking was performed to optimize the geometry of the protein-ligand complexes. The protein-ligand complexes were subjected to molecular dynamics simulation and binding free energy calculations. The residence time of a protein-ligand complex is a critical parameter affecting natural influences in vitro. It is nonetheless a challenging errand to expect, regardless of the accessibility of incredible PC assets and a large variety of computing procedures. In this metadynamics situation, we used the conformational flooding technique to deal with rank inhibitors constructions. The systematic free energy perturbation (FEP) protocol and calculate the energy of both complexes. Finally, the selected compound of TCM3290 was studied in vitro analysis such as inhibition of PLA2 activity, hyaluronidase activity and fibrinogenolytic activity. The TCM3290 had a more binding affinity compare to Minocycline, and interacted with the key residues of TYR63 and GLY31. DFT represented the highest HOMO and LUMO energy of 0.15146 eV. MD simulation with 100 ns proved that an inhibitor binding mode is more stable inside the binding site of PLA2. In vitro analysis shows that TCM3290 significantly neutralized by PLA2. The above observations confirmed that Scutellarin (TCM3290) had a potent snake venom neutralizing capacity and could hypothetically be used for therapeutic drives of snakebite envenomation.

Exploring the Papillomaviral Proteome to Identify Potential Candidates for a Chimeric Vaccine against Cervix Papilloma Using Immunomics and Computational Structural Vaccinology.

The human papillomavirus (HPV) 58 is considered to be the second most predominant genotype in cervical cancer incidents in China. HPV type-restriction, non-targeted delivery, and the highcost of existing vaccines necessitate continuing research on the HPV vaccine. We aimed to explore the papillomaviral proteome in order to identify potential candidates for a chimeric vaccine against cervix papilloma using computational immunology and structural vaccinology approaches. Two overlapped epitope segments (23⁻36) and (29⁻42) from the N-terminal region of the HPV58 minor capsid protein L2 are selected as capable of inducing both cellular and humoral immunity. In total, 318 amino acid lengths of the vaccine construct SGD58 contain adjuvants (Flagellin and RS09), two Th epitopes, and linkers. SGD58 is a stable protein that is soluble, antigenic, and non-allergenic. Homology modeling and the structural refinement of the best models of SGD58 and TLR5 found 96.8% and 93.9% favored regions in Rampage, respectively. The docking results demonstrated a HADDOCK score of -62.5 ± 7.6, the binding energy (-30 kcal/mol) and 44 interacting amino acid residues between SGD58-TLR5 complex. The docked complex are stable in 100 ns of simulation. The coding sequences of SGD58 also show elevated gene expression in Escherichia coli with 1.0 codon adaptation index and 59.92% glycine-cysteine content. We conclude that SGD58 may prompt the creation a vaccine against cervix papilloma.

Computational and in vitro insights on snake venom phospholipase A2 inhibitor of phytocompound ikshusterol3-O-glucoside of Clematis gouriana Roxb. ex DC.

Ikshusterol3-O-glucoside was isolated from Clematis gouriana Roxb. ex DC. root. A structure of the isolated compound was determined on the basis of various spectroscopic interpretations (UV, NMR, FTIR, and GC-MS-EI). This structure was submitted in the PubChem compound database (SID 249494133). SID 249494133 was carried out by density functional theory calculation to observe the chemical stability and electrostatic potential of this compound. The absorption, distribution, metabolism, and excretion property of this compound was predicted to evaluate the drug likeness and toxicity. In addition, molecular docking, quantum polarized ligand docking, prime MMGBSA calculation, and induced fit docking were performed to predict the binding status of SID 249494133 with the active site of phospholipase A2 (PLA2) (PDB ID: 1A3D). The stability of the compound in the active site of PLA2 was carried out using molecular dynamics simulation. Further, the anti-venom activity of the compound was assessed using the PLA2 assay against Naja naja (Indian cobra) crude venom. The results strongly show that Ikshusterol3-O-glucoside has a potent snake-venom neutralizing capacity and it might be a potential molecule for the therapeutic treatment for snakebites.

Isolation and characterization of bioactive compounds of Clematis gouriana Roxb. ex DC against snake venom phospholipase A2 (PLA2) computational and in vitro insights.

Bioactive compounds were isolated from Clematis gouriana Roxb. ex DC. The compounds were separated, characterized, the structures elucidated and submitted to the PubChem Database. The PubChem Ids SID 249494134 and SID 249494135 were tested against phospholipases A2 (PLA2) of Naja naja (Indian cobra) venom for PLA2 activity. Both the compounds showed promising inhibitory activity; computational data also substantiated the results. The two compounds underwent density functional theory calculation to observe the chemical stability and electrostatic potential profile. Molecular interactions between the compounds and PLA2 were observed at the binding pocket of the PLA2 protein. Further, this protein-ligand complexes were simulated for a timescale of 100 ns of molecular dynamics simulation. Experimental and computational results showed significant PLA2 inhibition activity.

High-affinity selective inhibitor against phospholipase A2 (PLA2): a computational study.

Phospholipase A2 (PLA2) is the most abundant protein found in snake venom. PLA2 induces a variety of pharmacological effects such as neurotoxicity, myotoxicity and cardiotoxicity as well as anticoagulant, hemolytic, anti-platelet, hypertensive, hemorrhagic and edema inducing effects. In this study, the three dimensional structure of PLA2 of Naja sputatrix (Malayan spitting cobra) was modeled by I-TASSER, SWISS-MODEL, PRIME and MODELLER programs. The best model was selected based on overall stereo-chemical quality. Further, molecular dynamics simulation was performed to know the stability of the modeled protein using Gromacs software. Average structure was generated during the simulation period of 10 ns. High throughput virtual screening was employed through different databases (Asinex, Hit finder, Maybridge, TOSLab and ZINC databases) against PLA2. The top seven compounds were selected based on the docking score and free energy binding calculations. These compounds were analyzed by quantum polarized ligand docking, induced fit docking and density functional theory calculation. Furthermore, the stability of lead molecules in the active site of PLA2 was employed by MD simulation. The results show that selected lead molecules were highly stable in the active site of PLA2.

Identification of potent inhibitors against snake venom metalloproteinase (SVMP) using molecular docking and molecular dynamics studies.

Snake venom metalloproteinase (SVMP) (Echis coloratus (Carpet viper) is a multifunctional enzyme that is involved in producing several symptoms that follow a snakebite, such as severe local hemorrhage, nervous system effects and tissue necrosis. Because the three-dimensional (3D) structure of SVMP is not known, models were constructed, and the best model was selected based on its stereo-chemical quality. The stability of the modeled protein was analyzed through molecular dynamics (MD) simulation studies. Structure-based virtual screening was performed, and 15 potential molecules with the highest binding energies were selected. Further analysis was carried out with induced fit docking, Prime/MM-GBSA (ΔGBind calculations), quantum-polarized ligand docking, and density functional theory calculations. Further, the stability of the lead molecules in the SVMP-active site was examined using MD simulation. The results showed that the selected lead molecules were highly stable in the active site of SVMP. Hence, these molecules could potentially be selective inhibitors of SVMP. These lead molecules can be experimentally validated, and their backbone structural scaffold could serve as building blocks in designing drug-like molecules for snake antivenom.

High throughput virtual screening and E-pharmacophore filtering in the discovery of new BACE-1 inhibitors.

Alzheimer's disease is a progressive neurodegenerative disorder, which is characterized by amyloid ß peptide deposition in the brain. Aß peptide, the major component of amyloid plaques is generated by the sequential processing of a larger protein called amyloid Precursor Protein by ß-amyloid cleaving enzyme (BACE-1). In this study, we appllied computer assisted methodology unifying molecular docking and pharmacophore filtering to identify potent inhibitors against BACE-1. In order to inspect the pharmacophore region and binding mode of BACE-1 135 reported co-crystallized ligands of BACE-1 were docked into the active site using Glide XP. The present molecular docking studies provided critical information on protein ligand interactions that revealed imminent information on chemical features essential to inhibiting BACE-1. Based on the docking results we proposed structure based pharmacophore features that hold well as potent BACE-1 inhibitors. A huge set of compounds was docked into the active site of BACE-1 and the hits from the docking were filtered to match the chemical features of the pharmacophore model. The compounds resulting from the pharmacophore filtering were again re-docked into the active site of BACE-1 and the three hits bound well into the active sites and matched the pharmacophore models which were identified as possible potential inhibitors of BACE-1. Molecular dynamics simulation reveals that lead 3 shows constant RMSD and the number of hydrogen bonding with the protein among the identified three lead molecules.