Mechanistic and structural insight into R2R3-MYB transcription factor in plants: molecular dynamics based binding free energy analysis.

The MYB transcription factor (TF) family is essential for various plant growth and development processes, including responses to biotic and abiotic stresses. This study investigated the R2R3-MYB protein structure from five plants, including cereal crops. The R2R3-MYB protein structure was docked with the DNA structure, and the best complexes were selected for two runs of molecular dynamics (MD) simulations to investigate the key interacting residues and conformational changes in the R2R3-MYB proteins caused by DNA binding. The MM/PBSA method calculated the binding free energy for each R2R3-MYB protein-DNA complex, showing strong interaction. Hydrophobic and hydrogen bonds significantly stabilized the R2R3-MYB protein-DNA complexes. The principal component analysis showed high restrictions on the movement of protein atoms in the phase space. A similar MD simulation analysis was performed using the crystal structure of the R2R3-MYB protein-DNA complex from Arabidopsis thaliana, and the generated complexes resembled the X-ray crystal structure. This is the first detailed study on the R2R3-MYB protein-DNA complex in cereal crops, providing a cost-effective solution to identify the key interacting residues and analyze the conformational changes in the MYB domain before and after DNA binding.Communicated by Ramaswamy H. Sarma.

Genome-wide identification and expression analysis of the NHX gene family under salt stress in wheat (Triticum aestivum L).

Salt stress affects plant growth and development, resulting in the loss of crop yield across the world, and sodium-proton antiporters (NHXs) are one of the genes known to promote salt tolerance in transgenic plants. In this study, we conducted a comprehensive genome-wide analysis and expression profile of NHX genes in wheat under salinity stress. We identified 30 TaNHX genes in wheat based on the Na+/H+ exchanger domain, with all genes containing an amiloride motif except one, a known for inhibiting Na+ ions in plants. Phylogenetic analysis classified these genes into three classes with subfamilies: 12 were localized in vacuoles, while 18 were in the endoplasmic reticulum and plasma membrane. Promoter analysis revealed stress-related cis-acting elements, indicating their potential role in abiotic stress tolerance. The non-synonymous (Ka)/synonymous (Ks) ratios highlighted that the majority of TaNHX genes experienced robust purifying selection throughout their evolutionary history. Transcriptomis data analysis and qRT-PCR demonstrated distinct expression patterns for TaNHX genes across various tissues when subjected to salt stress. Additionally, we predicted 20 different miRNA candidates targeting the identified TaNHX genes. Protein-protein interaction prediction revealed NHX6's involvement in the SOS1 pathway, while NHX1 gene exhibit proton antiporter activity. Molecular dynamics (MD) simulations were also conducted to examine the interactions of TaNHX1, TaNHX2, and TaNHX3. These results represent a significant advancement in our understanding of the molecular mechanisms governing Na+ transporters. This may also offer promising avenues for future studies aimed at unraveling the intricate details of their biological roles and applications.

Unveiling the Wheat Microbiome under Varied Agricultural Field Conditions.

Wheat being the important staple food crop plays a significant role in nutritional security. A wide variety of microbial communities beneficial to plants and contributing to plant health and production are found in the rhizosphere. The wheat microbiome encompasses an extensive variety of microbial species playing a key role in sustaining the physiology of the crop, nutrient uptake, and biotic/abiotic stress resilience. This report presents wheat microbiome analysis under six different farm practices, namely, organic (Org), timely sown (TS), wheat after pulse crop (WAPC), temperature-controlled phenotyping facility (TCPF), maize-wheat cropping system (MW), and residue burnt field (Bur), using 16S rRNA sequencing methodology. The soil samples collected from either side of the wheat row were mixed to get a final sample set for DNA extraction under each condition. After the data preprocessing, microbial community analysis was performed, followed by functional analysis and annotation. An abundance of the phylum Proteobacteria was observed, followed by Acidobacteria, Actinobacteria, and Gemmatimonadetes in the majority of the samples, while relative abundance was found to vary at the genus level. Analysis against the Carbohydrate-Active Enzymes (CAZy) database showed a high number of glycoside hydrolase genes in the TS, TCPF, and WAPC samples, while the Org, MW, and Bur samples predominantly had glycosyltransferase genes and carbohydrate esterase genes were in the lowest numbers. Also, the Org and TCPF samples showed lower diversity, while rare and abundant species ranged from 12 to 25% and 20 to 32% of the total bacterial species in all the sets, respectively. These variations indicate that the different cropping sequence had a significant impact on soil microbial diversity and community composition, which characterizes its economic and environmental value as a sustainable agricultural approach to maintaining food security and ecosystem health. IMPORTANCE This investigation examined the wheat microbiome under six different agricultural field conditions to understand the role of cropping pattern on soil microbial diversity. This study also elaborated the community composition, which has importance in economic (role of beneficial community leading to higher production) and environmental (role of microbial diversity/community in safeguarding the soil health, etc.) arenas. This could lead to a sustainable farming approach for food security and improved ecosystem health. Also, the majority of the microbes are unculturable; hence, technology-based microcultivation will be a potential approach for harnessing other cultured microorganisms, leading to unique species for commercial production. The outcome of this research-accelerated work can provide an idea to the scientists/breeders/agronomists/pathologists under the mentioned field conditions regarding their influence over their crops.

Structure-function analysis reveals Trichoderma virens Tsp1 to be a novel fungal effector protein modulating plant defence.

Trichoderma virens colonizes roots and develops a symbiotic relationship with plants where the fungal partner derives nutrients from plants and offers defence, in return. Tsp1, a small secreted cysteine-rich protein, was earlier found to be upregulated in co-cultivation of T. virens with maize roots. Tsp1 is well conserved in Ascomycota division of fungi, but none of its homologs have been studied yet. We have expressed and purified recombinant Tsp1, and resolved its structure to 1.25 Å resolutions, from two crystal forms, using Se-SAD methods. The Tsp1 adopts a ß barrel fold and forms dimer in structure as well as in solution form. DALI based structure analysis revealed the structure similarity with two known fungal effector proteins: Alt a1 and PevD1. Structure and evolutionary analysis suggested that Tsp1 belongs to a novel effector protein family. Tsp1 acted as an inducer of salicylic acid mediated susceptibility in plants, rendering maize plants more susceptible to a necrotrophic pathogen Cochliobolus heterostrophus, as observed using plant defence assay and RT-qPCR analysis.

Structural insights into SARS-CoV-2 proteins.

The unprecedented scale of the ongoing COVID-19 pandemic has catalyzed an intense effort of the global scientific community to unravel different aspects of the disease in a short time. One of the crucial aspects of these developments is the determination of more than three hundred experimental structures of SARS-CoV-2 proteins in the last few months. These include structures of viral non-structural, structural, and accessory proteins and their complexes determined by either X-ray diffraction or cryo-electron microscopy. These structures elucidate the intricate working of different components of the viral machinery at the atomic level during different steps of the viral life cycle, including attachment to the host cell, viral genome replication and transcription, and genome packaging and assembly of the virion. Some of these proteins are also potential targets for drug development against the disease. In this review, we discuss important structural features of different SARS-CoV-2 proteins with their function, and their potential as a target for therapeutic interventions.

Computational analysis identifies druggable mutations in human rBAT mediated Cystinuria.

Culex quinquefasciatus Cqm1 protein acts as the receptor for Lysinibacillus sphaericus mosquito-larvicidal binary (BinAB) toxin that is used worldwide for mosquito control. We found amino acid transporter protein, rBAT, as phylogenetically closest Cqm1 homolog in humans. The present study reveals large evolutionary distance between Cqm1 and rBAT, and rBAT ectodomain lacks the sequence motif which serves as binding-site for the BinAB toxin. Thus, BinAB toxin can be expected to remain safe for humans. rBAT (heavy subunit; SLC3A1) and catalytic b0,+AT (light subunit; SLC7A9), linked by single disulfide bond, mediate renal reabsorption of cystine and dibasic amino acids in Na+ independent manner. Mutations in rBAT cause type I Cystinuria disease which shows global prevalence, and rBAT can be thought as an important pharmacological target. However, 3D structures of rBAT and b0,+AT, the two components of b0,+ heteromeric amino acid transporter systems, are not available. We constructed a reliable homology model of rBAT using Cqm1 coordinates and that of transmembrane b0,+AT subunit using LAT1 coordinates. Mapping of pathogenic mutations onto rBAT ectodomain revealed their scattered distribution throughout the rBAT protein. Further, our computational simulations-based scoring of several known deleterious mutations of rBAT revealed that mutations those do not compromise the protein fold and stability, are localized on the same face of the molecule. These residues are expected to interact with the b0,+AT transporter. The present study thus identifies druggable sites on rBAT that could be targeted for the treatment of type I Cystinuria.Communicated by Ramaswamy H. Sarma.

Analysis of mutations of defensin protein using accelerated molecular dynamics simulations.

Plant defensins possess diverse biological functions that include antifungal and antibacterial activities and α-amylase and trypsin inhibitory properties. Two mutations, G9R and V39R, were confirmed to increase the antifungal activity of Raphanus sativus antifungal protein 2 (RsAFP2). Accelerated Molecular Dynamics (aMD) were carried out to examine the conformational changes present in these RsAFP2 mutants, and its two closest homologs compared to the wild-type protein. Specifically, the root mean square fluctuation values for the eight cysteine amino acids involved in the four disulfide bonds were low in the V39R mutant compared to the wild-type. Additionally, analysis of the free energy change revealed that G9R and V39R mutations exert a neutral and stabilizing effect on RsAFP2 conformation, and this is supported by the observed lower total energy of mutants compared to the wild-type, suggesting that enhanced stability of the mutants. However, MD simulations to a longer time scale would aid in capturing more conformational state of the wild-type and mutants defensin protein. Furthermore, the aMD simulations on fungal mimic membranes with RsAFP2 and its mutants and homologs showed that the mutant proteins caused higher deformation and water diffusion than the native RsAFP2, especially the V39R mutant. The mutant variants seem to interact by specifically targeting the POPC and POPI lipids amongst others. This work highlights the stabilizing effect of mutations at the 9th and 39th positions of RsAFP2 and their increased membrane deformation activity.

Analysis of mutations leading to para-aminosalicylic acid resistance in Mycobacterium tuberculosis.

Thymidylate synthase A (ThyA) is the key enzyme involved in the folate pathway in Mycobacterium tuberculosis. Mutation of key residues of ThyA enzyme which are involved in interaction with substrate 2'-deoxyuridine-5'-monophosphate (dUMP), cofactor 5,10-methylenetetrahydrofolate (MTHF), and catalytic site have caused para-aminosalicylic acid (PAS) resistance in TB patients. Focusing on R127L, L143P, C146R, L172P, A182P, and V261G mutations, including wild-type, we performed long molecular dynamics (MD) simulations in explicit solvent to investigate the molecular principles underlying PAS resistance due to missense mutations. We found that these mutations lead to (i) extensive changes in the dUMP and MTHF binding sites, (ii) weak interaction of ThyA enzyme with dUMP and MTHF by inducing conformational changes in the structure, (iii) loss of the hydrogen bond and other atomic interactions and (iv) enhanced movement of protein atoms indicated by principal component analysis (PCA). In this study, MD simulations framework has provided considerable insight into mutation induced conformational changes in the ThyA enzyme of Mycobacterium.

Identification of GCC-box and TCC-box motifs in the promoters of differentially expressed genes in rice (Oryza sativa L.): Experimental and computational approaches.

The transcription factor selectively binds with the cis-regulatory elements of the promoter and regulates the differential expression of genes. In this study, we aimed to identify and validate the presence of GCC-box and TCC-box motifs in the promoters of upregulated differentially expressed genes (UR-DEGs) and downregulated differentially expressed genes (DR-DEGs) under anoxia using molecular beacon probe (MBP) based real-time PCR. The GCC-box motif was detected in UR-DEGs (DnaJ and 60S ribosomal protein L7 genes), whereas, the TCC-box was detected in DR-DEGs (DnaK and CPuORF11 genes). In addition, the mechanism of interaction of AP2/EREBP family transcription factor (LOC_Os03g22170) with GCC-box promoter motif present in DnaJ gene (LOC_Os06g09560) and 60S ribosomal protein L7 gene (LOC_Os08g42920); and TCC-box promoter motif of DnaK gene (LOC_Os02g48110) and CPuORF11 gene (LOC_Os02g01240) were explored using molecular dynamics (MD) simulations analysis including binding free energy calculations, principal component analyses, and free energy landscapes. The binding free energy analysis revealed that AP2/EREBP model residues such as Arg68, Arg72, Arg83, Lys87, and Arg90 were commonly involved in the formation of hydrogen bonds with GCC and TCC-box promoter motifs, suggesting that these residues are critical for strong interaction. The movement of the entire protein bound to DNA was restricted, confirming the stability of the complex. This study provides comprehensive binding information and a more detailed view of the dynamic interaction between proteins and DNA.

Structural basis for isoniazid resistance in KatG double mutants of Mycobacterium tuberculosis.

The failure of drugs for effective treatment against infectious diseases can be attributed to resistant forms of causative agents. The evasive nature of Mycobacterium tuberculosis is partly associated to its physical features, such as having a thick cell wall and incorporation of beneficial mutations leading to drug resistance. The pro drug Isoniazid (INH) interacts with an enzyme catalase peroxidase to get converted into its active form and upon activation stops the cell wall synthesis thus killing the Mycobacterium. The most common mutation i.e. S315T leads to high degree of drug resistance by virtue of its position in the active site. Here, we have characterized the prominent attributes of two double mutant isolates S315 T/D194G and S315T/M624V which are multi drug resistant and extremely drug resistant, respectively. Protein models were generated using the crystal structure which were then subjected to energy minimization and long term molecular dynamics simulations. Further, computational analysis showed decreasing ability of INH binding to the mutants in order of: Native > S315T/D194G > S315T/M624V. Also, a trend was observed that as the docking score and binding area decreased, there was a significant increase in the distortion of the 3D geometry of the mutants as observed by PCA analysis.

Novel missense mutations in gidB gene associated with streptomycin resistance in Mycobacterium tuberculosis: insights from molecular dynamics.

Streptomycin was the first antibiotic used for the treatment of tuberculosis by inhibiting translational proof reading. Point mutation in gidB gene encoding S-adenosyl methionine (SAM)-dependent 7-methylguanosine (m7G) methyltransferase required for methylation of 16S rRNA confers streptomycin resistance. As there was no structural substantiation experimentally, gidB protein model was built by threading algorithm. In this work, molecular dynamics (MD) simulations coupled with binding free energy calculations were performed to outline the mechanism underlying high-level streptomycin resistance associated with three novel missense mutants including S70R, T146M, and R187M. Results from dynamics analyses suggested that the structure distortion in the binding pocket of gidB mutants modulate SAM binding affinity. At the structural level, these conformational changes bring substantial decrease in the number of residues involved in hydrogen bonding and dramatically reduce thermodynamic stability of mutant gidB-SAM complexes. The outcome of comparative analysis of the MD simulation trajectories revealed lower conformational stability associated with higher flexibility in mutants relative to the wild-type, turns to be major factor driving the emergence of drug resistance toward antibiotic. This study will pave way toward design and development of resistant defiant gidB inhibitors as potent anti-TB agents.

Drug targeted virtual screening and molecular dynamics of LipU protein of Mycobacterium tuberculosis and Mycobacterium leprae.

The lipolytic protein LipU was conserved in mycobacterium sp. including M. tuberculosis (MTB LipU) and M. leprae (MLP LipU). The MTB LipU was identified in extracellular fraction and was reported to be essential for the survival of mycobacterium. Therefore to address the problem of drug resistance in pathogen, LipU was selected as a drug target and the viability of finding out some FDA approved drugs as LipU inhibitors in both the cases was explored. Three-dimensional (3D) model structures of MTB LipU and MLP LipU were generated and stabilized through molecular dynamics (MD). FDA approved drugs were screened against these proteins. The result showed that the top-scoring compounds for MTB LipU were Diosmin, Acarbose and Ouabain with the Glide XP score of -12.8, -11.9 and -11.7 kcal/mol, respectively, whereas for MLP LipU protein, Digoxin (-9.2 kcal/mol), Indinavir (-8.2 kcal/mol) and Travoprost (-8.2 kcal/mol) showed highest affinity. These drugs remained bound in the active site pocket of MTB LipU and MLP LipU structure and interaction grew stronger after dynamics. RMSD, RMSF and Rg were found to be persistent throughout the simulation period. Hydrogen bonds along with large number of hydrophobic interactions stabilized the complex structures. Binding free energies obtained through Prime/MM-GBSA were found in the significant range from -63.85 kcal/mol to -34.57 kcal/mol for MTB LipU and -71.33 kcal/mol to -23.91 kcal/mol for MLP LipU. The report suggested high probability of these drugs to demolish the LipU activity and could be probable drug candidates to combat TB and leprosy disease.

Dynamics of Dof domain-DNA interaction in wheat: Insights from atomistic simulations and free energy landscape.

DNA-binding one zinc finger protein (Dof) is a plant-specific transcription factor involved in numerous biological processes. In the current study, the plausible mechanism underlying Dof domain-DNA interaction in wheat was investigated using extensive molecular dynamics (MD) simulations analysis. We depicted that one key residue Lys29, possessing the ability to disturb the interaction between Dof domain-DNA upon substitution to Arg29. Frequent conformational changes were observed in Lys29Arg (K29R)-DNA complex during the entire MD simulation period, which significantly altered the interactions, thereby indicating the importance of Lys29 in complex stability. Principal component analysis and free energy landscape results also suggested strong affinity between wild-type Dof domain and DNA due to restricted atomic movement. Our study not only substantiates the structural and mechanistic insights of Dof transcription factor but also provides new avenues toward employment of these key amino acid residues in genetic engineering for development of abiotic stress tolerance crop plant.

Molecular dynamics simulations revealed structural differences among WRKY domain-DNA interaction in barley (Hordeum vulgare).

BACKGROUND: The WRKY transcription factors are a class of DNA-binding proteins involved in diverse plant processes play critical roles in response to abiotic and biotic stresses. Genome-wide divergence analysis of WRKY gene family in Hordeum vulgare provided a framework for molecular evolution and functional roles. So far, the crystal structure of WRKY from barley has not been resolved; moreover, knowledge of the three-dimensional structure of WRKY domain is pre-requisites for exploring the protein-DNA recognition mechanisms. Homology modelling based approach was used to generate structures for WRKY DNA binding domain (DBD) and its variants using AtWRKY1 as a template. Finally, the stability and conformational changes of the generated model in unbound and bound form was examined through atomistic molecular dynamics (MD) simulations for 100 ns time period. RESULTS: In this study, we investigated the comparative binding pattern of WRKY domain and its variants with W-box cis-regulatory element using molecular docking and dynamics (MD) simulations assays. The atomic insight into WRKY domain exhibited significant variation in the intermolecular hydrogen bonding pattern, leading to the structural anomalies in the variant type and differences in the DNA-binding specificities. Based on the MD analysis, residual contribution and interaction contour, wild-type WRKY (HvWRKY46) were found to interact with DNA through highly conserved heptapeptide in the pre- and post-MD simulated complexes, whereas heptapeptide interaction with DNA was missing in variants (I and II) in post-MD complexes. Consequently, through principal component analysis, wild-type WRKY was also found to be more stable by obscuring a reduced conformational space than the variant I (HvWRKY34). Lastly, high binding free energy for wild-type and variant II allowed us to conclude that wild-type WRKY-DNA complex was more stable relative to variants I. CONCLUSIONS: The results of our study revealed complete dynamic and structural information about WRKY domain-DNA interactions. However, no structure base information reported to date for WRKY variants and their mechanism of interaction with DNA. Our findings highlighted the importance of selecting a sequence to generate newer transgenic plants that would be increasingly tolerance to stress conditions.

Alanine mutation of the catalytic sites of Pantothenate Synthetase causes distinct conformational changes in the ATP binding region.

The enzyme Pantothenate synthetase (PS) represents a potential drug target in Mycobacterium tuberculosis. Its X-ray crystallographic structure has demonstrated the significance and importance of conserved active site residues including His44, His47, Asn69, Gln72, Lys160 and Gln164 in substrate binding and formation of pantoyl adenylate intermediate. In the current study, molecular mechanism of decreased affinity of the enzyme for ATP caused by alanine mutations was investigated using molecular dynamics (MD) simulations and free energy calculations. A total of seven systems including wild-type + ATP, H44A + ATP, H47A + ATP, N69A + ATP, Q72A + ATP, K160A + ATP and Q164A + ATP were subjected to 50 ns MD simulations. Docking score, MM-GBSA and interaction profile analysis showed weak interactions between ATP (substrate) and PS (enzyme) in H47A and H160A mutants as compared to wild-type, leading to reduced protein catalytic activity. However, principal component analysis (PCA) and free energy landscape (FEL) analysis revealed that ATP was strongly bound to the catalytic core of the wild-type, limiting its movement to form a stable complex as compared to mutants. The study will give insight about ATP binding to the PS at the atomic level and will facilitate in designing of non-reactive analogue of pantoyl adenylate which will act as a specific inhibitor for PS.

Insight into the inhibitor discrimination by FLT3 F691L.

Fms-like tyrosine kinase 3 (FLT3) belongs to the receptor tyrosine kinase family and expressed in hematopoietic progenitor cells. FLT3 gene mutations are reported in ~30% of acute myeloid leukemia cases. FLT3 kinase domain mutation F691L is one of the common causes of acquired resistance to the FLT3 inhibitors including quizartinib. MZH29 and crenolanib were previously reported to inhibit FLT3 F691L. However, crenolanib was reported for the moderate inhibition. We found that Glu661and Asp829 were the most significant residues to target the FLT3 F691L which contribute most significantly to the binding energy with MZH29 and crenolanib. These interactions were found absent with quizartinib. Further free energy landscape analysis revealed that FLT3 F691L bound to MZH29 and crenolanib was more stable as compared to quizartinib.

Dynamics of fluoroquinolones induced resistance in DNA gyrase of Mycobacterium tuberculosis.

DNA gyrase is a validated target of fluoroquinolones which are key components of multidrug resistance tuberculosis (TB) treatment. Most frequent occurring mutations associated with high level of resistance to fluoroquinolone in clinical isolates of TB patients are A90V, D94G, and A90V-D94G (double mutant [DM]), present in the larger subunit of DNA Gyrase. In order to explicate the molecular mechanism of drug resistance corresponding to these mutations, molecular dynamics (MD) and mechanics approach was applied. Structure-based molecular docking of complex comprised of DNA bound with Gyrase A (large subunit) and Gyrase C (small subunit) with moxifloxacin (MFX) revealed high binding affinity to wild type with considerably high Glide XP docking score of -7.88 kcal/mol. MFX affinity decreases toward single mutants and was minimum toward the DM with a docking score of -3.82 kcal/mol. Docking studies were also performed against 8-Methyl-moxifloxacin which exhibited higher binding affinity against wild and mutants DNA gyrase when compared to MFX. Molecular Mechanics/Generalized Born Surface Area method predicted the binding free energy of the wild, A90V, D94G, and DM complexes to be -55.81, -25.87, -20.45, and -12.29 kcal/mol, respectively. These complexes were further subjected to 30 ns long MD simulations to examine significant interactions and conformational flexibilities in terms of root mean square deviation, root mean square fluctuation, and strength of hydrogen bond formed. This comparative drug interaction analysis provides systematic insights into the mechanism behind drug resistance and also paves way toward identifying potent lead compounds that could combat drug resistance of DNA gyrase due to mutations.

Wild-type catalase peroxidase vs G279D mutant type: Molecular basis of Isoniazid drug resistance in Mycobacterium tuberculosis.

Mycobacterium tuberculosis katG gene is responsible for production of an enzyme catalase peroxidase that peroxidises and activates the prodrug Isoniazid (INH), a first-line antitubercular agent. INH interacts with catalase peroxidase enzyme within its heme pocket and gets converted to an active form. Mutations occurring in katG gene are often linked to reduced conversion rates for INH. This study is focussed on one such mutation occurring at residue 279, where glycine often mutates to aspartic acid (G279D). In the present study, several structural analyses were performed to study the effect of this mutation on functionality of KatG protein. On comparison, mutant protein exhibited a lower docking score, smaller binding cavity and reduced affinity towards INH. Molecular dynamics analysis revealed the mutant to be more rigid and less compact than the native protein. Essential dynamics analysis determined correlated motions of residues within the protein structure. G279D mutant was found to have many residues that showed related motions and an undesirable effect on the functionality of protein.

Structural insights into impact of Y134F mutation and discovery of novel fungicidal compounds against CYP51 in Puccinia triticina.

Sterol 14α-Demethylase Cytochrome P450 (CYP51) protein involved in ergosterol biosynthesis pathways turn out to be a crucial target for the fungicidal compound. However, the recognition mechanism and dynamic behavior of CYP51 in wheat leaf rust pathogen, Puccinia triticina, is still obscure. Previously, a mutation at position 134 (Y134F) was reported in five European isolates of P. triticina, conversely, structural basis of this mutation remains unclear. To address this problem, three-dimensional structure of CYP51 protein from P. triticina was successfully built using homology modeling approach. To assess the protein structure stability, wild and mutant-type CYP51 proteins bound with azole fungicide was subjected to 50 ns molecular dynamics (MD) simulations run. Observably, the comparative protein-ligand interaction analysis and binding free energy results revealed that impact of the mutation on the thermodynamics and conformational stability of the CYP51 protein was negligible. In addition, we carried out structure-based virtual screening and identified potent novel fungicidal compounds from four different databases and libraries. Consequently, through MD simulation and thermodynamic integration, four novel compounds such as CoCoCo54211 (CoCoCo database), ZINC04089470 (ZINC database), Allyl pyrocatechol 3,4 diacetate (Natural compound library), and 9-octadecenoic acid (Traditional Chinese Medicine database) has been predicted as potent fungicidal compound against CYP51 with XPGlide docking score of -11.41, -13.64, -7.40, and -6.55 kcal/mol, respectively. These compounds were found to form hydrogen bonds with heme group of CYP51, subsequently disturbing the stability and survival of fungus and can be used to control leaf rust in wheat.

Role of pncA gene mutations W68R and W68G in pyrazinamide resistance.

Mycobacterium tuberculosis (Mtb) resistance toward anti-tuberculosis drugs is a widespread problem. Pyrazinamide (PZA) is a first line antitubercular drug that kills semi-dormant bacilli when converted into its activated form, that is, pyrazinoic acid (POA) by Pyrazinamidase (PZase) enzyme coded by pncA gene. In this study, we conducted several analyses on native and mutant structures (W68R, W68G) of PZase before and after docking with the PZA drug to explore the molecular mechanism behind PZA resistance caused due to pncA mutations. Structural changes caused by mutations were studied with respect to their effects on functionality of protein. Docking was performed to analyze the protein-drug binding and comparative analysis was done to observe how the mutations affect drug binding affinity and binding site on protein. Native PZase protein was observed to have the maximum binding affinity in terms of docking score as well as shape complementarity in comparison to the mutant forms. Molecular dynamics simulation analyses showed that mutation in the 68th residue of protein results in a structural change at its active site which further affects the biological function of protein, that is, conversion of PZA to POA. Mutations in the protein thereby led to PZA resistance in the bacterium due to the inefficient binding.