Pyrazinamide resistance of novel mutations in pncA and their dynamic behavior

Mechanistic insights into the conformational changes and alterations in residual communications due to the mutations in the pncA Gene of Mycobacterium tuberculosis: A computational perspective for effective therapeutic solutions

Due to its emerging resistance to first-line anti-TB medications, tuberculosis (TB) is one of the most contagious illness in the world. According to reports, the effectiveness of treating TB is severely impacted by drug resistance, notably resistance caused by mutations in the pncA gene-encoded pyrazinamidase (PZase) to the front-line drug pyrazinamide (PZA). The present study focused on investigating the resistance mechanism caused by the mutations D12N, T47A, and H137R to better understand the structural and molecular events responsible for the resistance acquired by the pncA gene of Mycobacterium tuberculosis (MTB) at the structural level. Bioinformatics analysis predicted that all three mutations were deleterious and located near the active centre of the pncA, affecting its functional activity. Furthermore, molecular dynamics simulation (MDS) results established that mutations significantly reduced the structural stability and caused the rearrangement of FE2+ in the active centre of pncA. Moreover, essential dynamics analysis, including principal component analysis (PCA) and free energy landscape (FEL), concluded variations in the protein motion and decreased conformational space in the mutants. Additionally, the mutations potentially impacted the network topologies and altered the residual communications in the network. The complex simulation study results established the significant movement of the flap region from the active centre of mutant complexes, further supporting the flap region's significance in developing resistance to the PZA drug. This study advances our knowledge of the primary cause of the mechanism of PZA resistance and the structural dynamics of pncA mutants, which will help us to design new and potent chemical scaffolds to treat drug-resistant TB (DR-TB).

Molecular modeling of novel 2-aminopyridine derivatives as potential JAK2 inhibitors: a rational strategy for promising anticancer agents

Janus kinase 2(JAK2) is a potential target for anticancer drugs in the treatment of numerous myeloproliferative diseases due to its central role in the JAK/STAT signaling cascade. In this study, the binding behavior of 2 amino-pyridine derivatives as JAK2 inhibitors was investigated by using multifaceted strategies including 3D-QSAR, molecular docking, Fingerprint analysis, MD simulations, and MM-PBSA calculations. A credible COMFA (q2 = 0.606 and r2 = 0.919) and COMSIA (q2 = 0.641 and r2 = 0.992) model was developed, where the internal and external validation revealed that the obtained 3D-QSAR models could be capable of predicting bioactivities of JAK2 inhibitors. The structural criteria provided by the contour maps of model were used to computationally develop more potent 100 new JAK2 inhibitors. Docking studies were conducted on the model data set and newly developed compounds (in-house library) to demonstrate their binding mechanism and highlight the key interacting residues within JAK2 active site. The selected docked complexes underwent MD simulation (100 ns), which contributed in the further study of the binding interactions. Binding free energy analyses (MMGB/PBSA) revealed that key residues such as Glu930, Leu932 (hinge region), Asp939 (solvent accessible region), Arg980, Asn981and Asp994 (catalytic site) have a significantly facilitate ligand-protein interactions through H-bonding and van der Waals interactions. The preliminary in-silico ADMET evaluation revealed encouraging results for all the modeled and in-house library compounds. The findings of this research have the potential to offer valuable recommendations for the advancement of novel, potent, and efficacious JAK2 inhibitors. Overall, this work has successfully employed a wide range of computer-based methodologies to understand the interaction dynamics between 2-amino-pyridine derivatives and the JAK2 enzyme, which is a crucial target in myeloproliferative disorders.Communicated by Ramaswamy H. Sarma.

Free energy landscape and thermodynamics properties of novel mutations in PncA of pyrazinamide resistance isolates of Mycobacterium tuberculosis

Pyrazinamide (PZA) is one of the first-line antituberculosis therapy, active against non-replicating Mycobacterium tuberculosis (Mtb). The conversion of PZA into pyrazinoic acid (POA), the active form, required the activity of pncA gene product pyrazinamidase (PZase) activity. Mutations occurred in pncA are the primary cause behind the PZA resistance. However, the resistance mechanism is important to explore using high throughput computational approaches. Here we aimed to explore the mechanism of PZA resistance behind novel P62T, L120R, and V130M mutations in PZase using 200 ns molecular dynamics (MD) simulations. MD simulations were performed to observe the structural changes for these three mutants (MTs) compared to the wild types (WT). Root means square fluctuation, the radius of gyration, free energy landscape, root means square deviation, dynamic cross-correlation motion, and pocket volume were found in variation between WT and MTs, revealing the effects of P62T, L120R, and V130M. The free energy conformational landscape of MTs differs significantly from the WT system, lowering the binding of PZA. The geometric shape complementarity of the drug (PZA) and target protein (PZase) further confirmed that P62T, L120R, and V130M affect the protein structure. These effects on PZase may cause vulnerability to convert PZA into POA.Communicated by Ramaswamy H. Sarma.

Deciphering the mechanism of resistance by novel double mutations in pncA in Mycobacterium tuberculosis using protein structural graphs (PSG) and structural bioinformatic approaches

The evolution of MDR and XDR-TB is a growing concern and public health safety threat around the world. Gene mutations are the prime cause of drug resistance in tuberculosis, however the reports of double mutations further aggravated the situation. Despite the large-scale genomic sequencing and identification of novel mutations, structure investigation of the protein is still required to structurally and functionally characterize these novel mutations to design novel drugs for improved clinical outcome. Hence, we used structural bioinformatics approaches i.e. molecular modeling, residues communication and molecular simulation to understand the impact of novel double S59Y-L85P, D86G-V180F and S104G-V130 M mutation on the structure, function of pncA encoded Pyrazinamidase (PZase) and resistance of Pyrazinamide (PZA). Our results revealed that these mutations alter the binding paradigm and destabilize the protein to release the drug. Protein commination network (PCN) revealed variations in the hub residues and sub-networks which consequently alter the internal communication and signaling. The region 1-75 demonstrated higher flexibility in the mutant structures and minimal by the wild type which destabilize of the internally arranged beta-sheets which consequently reduce the binding of PZA and potentially Fe ion in the mutants. Hydrogen bonding analysis further validated the findings. The total binding free energy (ΔG) for each complex i.e. wild type -7.46 kcal/mol, S59Y-L85P -5.21 kcal/mol, S104G-V130 M -5.33 kcal/mol while for the D86G-V180F mutant the TBE was calculated to be -6.26 kcal/mol. This further confirms that these mutations reduce the binding energy of PZA for PZase and causes resistance in the effective therapy for TB. The trajectories motion was also observed to be affected by these mutations. In conclusion, these mutations use destabilizing approach to reduce the binding of PZA and causes resistance. These features can be used to design novel structure-based drugs against Tuberculosis.

Attenuation of Lipopolysaccharide-Induced Inflammatory Responses through Inhibition of the NF-κB Pathway and the Increased NRF2 Level by a Flavonol-Enriched n-Butanol Fraction from Uvaria alba

Pathologic hyperreactive inflammatory responses occur when there is excessive activation of a proinflammatory NF-κB pathway and a reduced cytoprotective NRF2 cascade. The noncytotoxic, highly selective COX-2 inhibitory flavonol-enriched butanol fraction (UaB) from Uvaria alba (U. alba) was investigated for its inflammatory modulating potential by targeting NF-κB activation and NRF2 activity. Enzyme-linked immunosorbent assay was initially performed to measure levels of proinflammatory mediators [nitric oxide (NO), prostaglandin E₂, and reactive oxygen species (ROS)] and cytokines [tumor necrosis factor-alpha (TNF-α), IL-1β, and IL-6], followed by reverse transcription-polymerase chain reaction and western blotting to determine mRNA and protein expression, respectively. Using immunofluorescence staining combined with western blot analysis, the activation of NF-κB was further investigated. NRF2 activity was also measured using a luciferase reporter assay. UaB abrogated protein and mRNA expressions of inducible nitric oxide synthase (iNOS), COX-2, TNF-α, IL-1β, and IL-6 in RAW 264.7 macrophages, thereby suppressing the production of proinflammatory mediators and cytokines. This was further validated when a concentration-dependent decrease in NO and ROS production was observed in zebrafish (Danio rerio) larvae. UaB also increased NRF2 activity in HaCaT/ARE cell line and attenuated NF-κB activation by inhibiting the nuclear translocation of transcription factor p65 in RAW 264.7 macrophages. Nontargeted LC-MS analysis of UaB revealed the presence of the flavonols quercitrin (1), quercetin (2), rutin (3), kaempferol (4), and kaempferol 3-O-rutinoside (5). Molecular docking indicates that major flavonol aglycones have high affinity toward COX-2 NSAID-binding sites, TNF-α, and TNF-α converting enzyme, while the glycosylated flavonoids showed strong binding toward iNOS and IKK-all possessing dynamic stability when performing molecular dynamics simulations at 140 ns. This is the first report to have elucidated the mechanistic anti-inflammatory potential of the Philippine endemic plant U. alba.

Investigation of the binding and dynamic features of A.30 variant revealed higher binding of RBD for hACE2 and escapes the neutralizing antibody: A molecular simulation approach

With the emergence of Delta and Omicron variants, many other important variants of SARS-CoV-2, which cause Coronavirus disease-2019, including A.30, are reported to increase the concern created by the global pandemic. The A.30 variant, reported in Tanzania and other countries, harbors spike gene mutations that help this strain to bind more robustly and to escape neutralizing antibodies. The present study uses molecular modelling and simulation-based approaches to investigate the key features of this strain that result in greater infectivity. The protein-protein docking results for the spike protein demonstrated that additional interactions, particularly two salt-bridges formed by the mutated residue Lys484, increase binding affinity, while the loss of key residues at the N terminal domain (NTD) result in a change to binding conformation with monoclonal antibodies, thus escaping their neutralizing effects. Moreover, we deeply studied the atomic features of these binding complexes through molecular simulation, which revealed differential dynamics when compared to wild type. Analysis of the binding free energy using MM/GBSA revealed that the total binding free energy (TBE) for the wild type receptor-binding domain (RBD) complex was -58.25 kcal/mol in contrast to the A.30 RBD complex, which reported -65.59 kcal/mol. The higher TBE for the A.30 RBD complex signifies a more robust interaction between A.30 variant RBD with ACE2 than the wild type, allowing the variant to bind and spread more promptly. The BFE for the wild type NTD complex was calculated to be -65.76 kcal/mol, while the A.30 NTD complex was estimated to be -49.35 kcal/mol. This shows the impact of the reported substitutions and deletions in the NTD of A.30 variant, which consequently reduce the binding of mAb, allowing it to evade the immune response of the host. The reported results will aid the development of cross-protective drugs against SARS-CoV-2 and its variants.

The Omicron (B.1.1.529) variant of SARS-CoV-2 binds to the hACE2 receptor more strongly and escapes the antibody response: Insights from structural and simulation data

As SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) continues to inflict chaos globally, a new variant officially known as B.1.1.529 was reported in South Africa and was found to harbor 30 mutations in the spike protein. It is too early to speculate on transmission and hospitalizations. Hence, more analyses are required, particularly to connect the genomic patterns to the phenotypic attributes to reveal the binding differences and antibody response for this variant, which can then be used for therapeutic interventions. Given the urgency of the required analysis and data on the B.1.1.529 variant, we have performed a detailed investigation to provide an understanding of the impact of these novel mutations on the structure, function, and binding of RBD to hACE2 and mAb to the NTD of the spike protein. The differences in the binding pattern between the wild type and B.1.1.529 variant complexes revealed that the key substitutions Asn417, Ser446, Arg493, and Arg498 in the B.1.1.529 RBD caused additional interactions with hACE2 and the loss of key residues in the B.1.1.529 NTD resulted in decreased interactions with three CDR regions (1-3) in the mAb. Further investigation revealed that B.1.1.529 displayed a stable dynamic that follows a global stability trend. In addition, the dissociation constant (K_D), hydrogen bonding analysis, and binding free energy calculations further validated the findings. Hydrogen bonding analysis demonstrated that significant hydrogen bonding reprogramming took place, which revealed key differences in the binding. The total binding free energy using MM/GBSA and MM/PBSA further validated the docking results and demonstrated significant variations in the binding. This study is the first to provide a basis for the higher infectivity of the new SARS-CoV-2 variants and provides a strong impetus for the development of novel drugs against them.

Structural and Dynamic Insights into the W68L, L85P, and T87A Mutations of Mycobacterium tuberculosis Inducing Resistance to Pyrazinamide

Tuberculosis (TB), the most frequent bacterium-mediated infectious disease caused by Mycobacterium tuberculosis, has been known to infect humans since ancient times. Although TB is common worldwide, the most recent report by the WHO (World Health Organization) listed the three countries of India, China, and Russia with 27%, 14%, and 8% of the global burden of TB, respectively. It has been reported that resistance to TB drugs, particularly by the pncA gene to the pyrazinamide drug due to mutations, significantly affects the effective treatment of TB. Understanding the mechanism of drug resistance using computational methods is of great interest to design effective TB treatment, exploring the structural features with these tools. Thus, keeping in view the importance of these methods, we employed state-of-the-art computational methods to study the mechanism of resistance caused by the W68L, L85P, and T87A mutations recently reported in 2021. We employed a molecular docking approach to predict the binding conformation and studied the dynamic properties of each complex using molecular dynamics simulation approaches. Our analysis revealed that compared to the wildtype, these three mutations altered the binding pattern and reduced the binding affinity. Moreover, the structural dynamic features also showed that these mutations significantly reduced the structural stability and packing, particularly by the W68L and L85P mutations. Moreover, principal component analysis, free energy landscape, and the binding free energy results revealed variation in the protein’s motion and the binding energy. The total binding free energy was for the wildtype −9.61 kcal/mol, W68L −7.57 kcal/mol, L85P −6.99 kcal/mol, and T87A −7.77 kcal/mol. Our findings can help to design a structure-based drug against the MDR (multiple drug-resistant) TB.

Structural and Genomic Insights Into Pyrazinamide Resistance in Mycobacterium tuberculosis Underlie Differences Between Ancient and Modern Lineages

Resistance to drugs used to treat tuberculosis disease (TB) continues to remain a public health burden, with missense point mutations in the underlying Mycobacterium tuberculosis bacteria described for nearly all anti-TB drugs. The post-genomics era along with advances in computational and structural biology provide opportunities to understand the interrelationships between the genetic basis and the structural consequences of M. tuberculosis mutations linked to drug resistance. Pyrazinamide (PZA) is a crucial first line antibiotic currently used in TB treatment regimens. The mutational promiscuity exhibited by the pncA gene (target for PZA) necessitates computational approaches to investigate the genetic and structural basis for PZA resistance development. We analysed 424 missense point mutations linked to PZA resistance derived from ∼35K M. tuberculosis clinical isolates sourced globally, which comprised the four main M. tuberculosis lineages (Lineage 1-4). Mutations were annotated to reflect their association with PZA resistance. Genomic measures (minor allele frequency and odds ratio), structural features (surface area, residue depth and hydrophobicity) and biophysical effects (change in stability and ligand affinity) of point mutations on pncA protein stability and ligand affinity were assessed. Missense point mutations within pncA were distributed throughout the gene, with the majority (>80%) of mutations with a destabilising effect on protomer stability and on ligand affinity. Active site residues involved in PZA binding were associated with multiple point mutations highlighting mutational diversity due to selection pressures at these functionally important sites. There were weak associations between genomic measures and biophysical effect of mutations. However, mutations associated with PZA resistance showed statistically significant differences between structural features (surface area and residue depth), but not hydrophobicity score for mutational sites. Most interestingly M. tuberculosis lineage 1 (ancient lineage) exhibited a distinct protein stability profile for mutations associated with PZA resistance, compared to modern lineages.