Insights into the structure-function relationship of disease resistance protein HCTR in maize (Zea mays L.): a computational structural biology approach

Spectroscopic insights with molecular docking and molecular dynamic simulation studies of anticancer drug 5-Fluorouracil targeting human pyruvate kinase m2

This study was conducted to test the efficacy of 5-fluorouracil (5-FU) as an anticancer drug against the human pyruvate kinase isozyme M2 (PKM2) using spectroscopic, molecular docking and molecular dynamic simulation studies. PKM2 fluorescence quenching studies in the presence of 5-FU performed at three different temperatures indicates dynamic quenching processes with single-set of binding (n ≈ 1) profile. The biomolecular quenching constants (kq) and the effective binding constants (Kb) obtained are shown to increase with temperature. The calculated enthalpy (ΔH) and entropy changes (ΔS) are estimated to be -118.06 kJ/mol and 146.14 kJ/mol/K respectively, which suggest the possible mode of interaction as electrostatic and hydrogen bonding. Further, these values were used to estimate the free energy changes (ΔG) and that increases with temperature. The negative ΔG values clearly indicates spontaneous binding process that stabilizes the complex formed between 5-FU and PKM2. Far-UV CD spectra of PKM2 in the presence of 5-FU shows decrease in α-helix contents which point towards the destabilization of secondary structure that weakens the biological activity of PKM2. The intrinsic fluorescence study and circular dichroism (CD) spectra showed minor conformational changes of PKM2 in the presence of 5-FU. Additionally, the results obtained from molecular docking and all-atom molecular dynamic simulation study supports the insight of the spectroscopic binding studies, and strengthens the dynamic stability of the complex between 5-FU and PKM2 through H-bonding. This study establishes a paradigm of 5-FU-PKM2 complexation and the efficacy of 5-FU that compromises the biological activity of the targeted PKM2.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulation, synthesis and topoisomerase inhibitory actions of vanillin derivatives: a systematic computational structural integument

A series of 4-hydroxy-3-methoxy benzaldehyde (vanillin) derivatives (3a-3r) was designed for the principle of Schiff base condensation with several individual sulfanilamide analogues. The inhibitory potencies of the designed compounds were evaluated through molecular docking simulation studies against the targets, breast cancer-topo isomerase-IIα and estrogen receptor-α; and the top scoring poses with higher binding energy were selected to assess the mode of binding and stability of each complex through molecular dynamics simulations. Compounds that remained stable in the active sites of the both target receptors through a number of strong H-bonds and hydrophobic contacts were selected. Based on the computational results, these selected compounds, 3b, 3e and 3f were synthesized and were followed up for structural elucidation attempts, by FT/ATR, ¹H NMR and ¹³C NMR. From the experimental in vitro studies on 3b, 3e and 3f, the following remarkable activities against breast cancer cell line were done; IC₅₀ values of 3b, 3e and 3f were noted, 6.7, 4.3 and 11 ng/mL, respectively. These newly synthesized compounds may be used as novel inhibitors of nuclear receptors with potential therapeutic applications in control of cancer.Communicated by Ramaswamy H. Sarma.

Novel Insights into Understanding the Molecular Dialogues between Bipolaroxin and the Gα and Gβ Subunits of the Wheat Heterotrimeric G-Protein during Host–Pathogen Interaction

Spot blotch disease of wheat, caused by the fungus Bipolaris sorokiniana (Sacc.) Shoem., produces several toxins which interact with the plants and thereby increase the blightening of the wheat leaves, and Bipolaroxin is one of them. There is an urgent need to decipher the molecular interaction between wheat and the toxin Bipolaroxin for in-depth understanding of host–pathogen interactions. In the present study, we have developed the three-dimensional structure of G-protein alpha subunit from Triticum aestivum. Molecular docking studies were performed using the active site of the modeled G-protein alpha and cryo-EM structure of beta subunit from T. aestivum and ‘Bipolaroxin’. The study of protein–ligand interactions revealed that six H-bonds are mainly formed by Glu29, Ser30, Lys32, and Ala177 of G-alpha with Bipolaroxin. In the beta subunit, the residues of the core beta strand domain participate in the ligand interaction where Lys256, Phe306, and Leu352 formed seven H-bonds with the ligand Bipolaroxin. All-atoms molecular dynamics (MD) simulation studies were conducted for G-alpha and -beta subunit and Bipolaroxin complexes to explore the stability, conformational flexibility, and dynamic behavior of the complex system. In planta studies clearly indicated that application of Bipolaroxin significantly impacted the physio-biochemical pathways in wheat and led to the blightening of leaves in susceptible cultivars as compared to resistant ones. Further, it interacted with the Gα and Gβ subunits of G-protein, phenylpropanoid, and MAPK pathways, which is clearly supported by the qPCR results. This study gives deeper insights into understanding the molecular dialogues between Bipolaroxin and the Gα and Gβ subunits of the wheat heterotrimeric G-protein during host–pathogen interaction.

Transcriptome-Based Analysis of Tomato Genotypes Resistant to Bacterial Spot (Xanthomonas perforans) Race T4

Bacterial spot (BS) is one of the most devastating foliar bacterial diseases of tomato and is caused by multiple species of Xanthomonas. We performed the RNA sequencing (RNA-Seq) analysis of three tomato lines with different levels of resistance to Xanthomonas perforans race T4 to study the differentially expressed genes (DEGs) and transcript-based sequence variations. Analysis between inoculated and control samples revealed that resistant genotype Solanum pimpinellifolium accession PI 270443 had more DEGs (834), followed by susceptible genotype tomato (S. lycopersicum L) breeding line NC 714 (373), and intermediate genotype tomato breeding line NC 1CELBR (154). Gene ontology (GO) terms revealed that more GO terms (51) were enriched for upregulated DEGs in the resistant genotype PI 270443, and more downregulated DEGs (67) were enriched in the susceptible genotype NC 714. DEGs in the biotic stress pathway showed more upregulated biotic stress pathway DEGs (67) for PI 270443 compared to more downregulated DEGs (125) for the susceptible NC 714 genotype. Resistant genotype PI 270443 has three upregulated DEGs for pathogenesis-related (PR) proteins, and susceptible genotype NC 714 has one downregulated R gene. Sequence variations called from RNA-Seq reads against the reference genome of susceptible Heinz 1706 showed that chr11, which has multiple reported resistance quantitative trait loci (QTLs) to BS race T4, is identical between two resistant lines, PI 270443 and NC 1CELBR, suggesting that these two lines share the same resistance QTLs on this chromosome. Several loci for PR resistance proteins with sequence variation between the resistant and susceptible tomato lines were near the known Rx4 resistance gene on chr11, and additional biotic stress associated DEGs near to the known Rx4 resistance gene were also identified from the susceptible NC 714 line.

A combination of 2D similarity search, pharmacophore, and molecular docking techniques for the identification of vascular endothelial growth factor receptor-2 inhibitors

The human vascular endothelial growth factor receptor-2 (VEGFR-2) has been an attractive target for the inhibition of angiogenesis. In the current study, we used a hybrid protocol of virtual screening methods to retrieve new VEGFR-2 inhibitors from the Zinc-Specs Database (441 574 compounds). The hybrid protocol included the initial screening of candidates by comparing the 2D similarity to five reported top active inhibitors of 13 VEGFR-2 X-ray crystallography structures, followed by the pharmacophore modeling of virtual screening on the basis of receptor-ligand interactions and further narrowing by LibDOCK to obtain the final hits. Two compounds (AN-919/41439526 and AK-968/40939851) with a high libscore were selected as the final hits for a subsequent cell cytotoxicity study. The two compounds screened exerted significant inhibitory effects on the proliferation of cancer cells (U87 and MCF-7). The results indicated that the hybrid procedure is an effective approach for screening specific receptor inhibitors.

The blast resistance gene Pi54of cloned from Oryza officinalis interacts with Avr-Pi54 through its novel non-LRR domains

The dominant rice blast resistance gene Pi54 cloned by map-based cloning approach from indica rice cultivar Tetep confers broad spectrum resistance to Magnaporthe oryzae. In this investigation, an orthologue of Pi54 designated as Pi54of was cloned from Oryza officinalis conferring high degree of resistance to M. oryzae and is functionally validated. We have also characterized the Pi54of protein and demonstrate its interaction with AVR-Pi54 protein. The Pi54of encoded ∼43 kDa small and unique cytoplasmic LRR family of disease resistance protein having unique Zinc finger domain overlapped with the leucine rich repeat regions. Pi54of showed Magnaporthe-induced expression. The phylogenetic and western blot analysis confirmed orthologous nature of Pi54 and Pi54of genes, whereas the identity of protein was confirmed through MALDI-TOF analysis. The in silico analysis showed that Pi54of is structurally more stable than other cloned Pi54 proteins. The molecular docking revealed that Pi54of protein interacts with AVR-Pi54 through novel non-LRR domains such as STI1 and RhoGEF. The STI1 and GEF domains which interact with AVR-Pi54 are also components of rice defensome complex. The Pi54of protein showed differential domain specificity while interacting with the AVR protein. Functional complementation revealed that Pi54of transferred in two rice lines belonging to indica and japonica background imparts enhanced resistance against three highly virulent strains of M. oryzae. In this study, for the first time, we demonstrated that a rice blast resistance gene Pi54of cloned from wild species of rice provides high degree of resistance to M. oryzae and might display different molecular mechanism involved in AVRPi54-Pi54of interaction.

Structure-based computational study of two disease resistance gene homologues (Hm1 and Hm2) in maize (Zea mays L.) with implications in plant-pathogen interactions

The NADPH-dependent HC-toxin reductases (HCTR1 and 2) encoded by enzymatic class of disease resistance homologous genes (Hm1 and Hm2) protect maize by detoxifying a cyclic tetrapeptide, HC-toxin, secreted by the fungus Cochliobolus carbonum race 1(CCR1). Unlike the other classes' resistance (R) genes, HCTR-mediated disease resistance is an inimitable mechanism where the avirulence (Avr) component from CCR1 is not involved in toxin degradation. In this study, we attempted to decipher cofactor (NADPH) recognition and mode of HC-toxin binding to HCTRs through molecular docking, molecular dynamics (MD) simulations and binding free energy calculation methods. The rationality and the stability of docked complexes were validated by 30-ns MD simulation. The binding free energy decomposition of enzyme-cofactor complex was calculated to find the driving force behind cofactor recognition. The overall binding free energies of HCTR1-NADPH and HCTR2-NADPH were found to be −616.989 and −16.9749 kJ mol⁻¹ respectively. The binding free energy decomposition revealed that the binding of NADPH to the HCTR1 is mainly governed by van der Waals and nonpolar interactions, whereas electrostatic terms play dominant role in stabilizing the binding mode between HCTR2 and NADPH. Further, docking analysis of HC-toxin with HCTR-NADPH complexes showed a distinct mode of binding and the complexes were stabilized by a strong network of hydrogen bond and hydrophobic interactions. This study is the first in silico attempt to unravel the biophysical and biochemical basis of cofactor recognition in enzymatic class of R genes in cereal crop maize.