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

MP-13, a novel chimeric peptide of morphiceptin and pepcan-9, produces potent antinociception with limited side effects

Pharmacological investigations have substantiated the potential of bifunctional opioid/cannabinoid agonists in delivering potent analgesia while minimizing adverse reactions. Peptide modulators of cannabinoid receptors, known as pepcans, have been investigated before. In this study, we designed a series of chimeric peptides based on pepcans and morphiceptin (YPFP-NH2). Here, we combined injections of pepcans and morphiceptin to investigate the combination treatment of opioids and cannabis and compared the analgesic effect with chimeric compounds. Subsequently, we employed computational docking to screen the compounds against opioid and cannabinoid receptors, along with an acute pain model, to identify the most promising peptide. Among these peptides, MP-13, a morphiceptin and pepcan-9 (PVNFKLLSH) construct, exhibited superior supraspinal analgesic efficacy in the tail-flick test, with an ED50 value at 1.43 nmol/mouse, outperforming its parent peptides and other chimeric analogs. Additionally, MP-13 displayed potent analgesic activity mediated by mu-opioid receptor (MOR), delta-opioid receptor (DOR), and cannabinoid type 1 (CB1) receptor pathways. Furthermore, MP-13 did not induce psychological dependence and gastrointestinal motility inhibition at the effective analgesic doses, and it maintained non-tolerance-forming antinociception throughout a 7-day treatment regimen, with an unaltered count of microglial cells in the periaqueductal gray region, supporting this observation. Moreover, intracerebroventricular administration of MP-13 demonstrated dose-dependent antinociception in murine models of neuropathic, inflammatory, and visceral pain. Our findings provide promising insights for the development of opioid/cannabinoid peptide agonists, addressing a crucial gap in the field and holding significant potential for future research and development. PERSPECTIVE: This article offers insights into the combination treatment of pepcans with morphiceptin. Among the chimeric peptides, MP-13 exhibited potent analgesic effects in a series of preclinical pain models with a favorable side-effect profile.

Structure and dynamics of whole-sequence homology model of ORF3a protein of SARS-CoV-2: An insight from microsecond molecular dynamics simulations

The ORF3a is a large accessory protein in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which plays an important role in virulence and viral replication; especially in inflammasome activation and apoptosis. However,, the existing cryo-EM structure of SARS-CoV-2 ORF3a is incomplete, . making it challenging to understand its structural and functional features. The aim of this study is to investigate the dynamic behaviors of the full-sequence homology model of ORF3a and compare it with the cryo-EM structure using microsecond molecular dynamics simulations. The previous studies indicated that the unresolved residues of the cryo-EM structure are not only involved in the pathogenesis of the SARS-CoV-2 but also exhibit a significant antigenicity. The dynamics scenario of homology model revealed higher RMSD, Rg, and SASA values with stable pattern when compared to the cryo-EM structure. Moreover, the RMSF analysis demonstrated higher fluctuations at specific positions (1-43, 97-110, 172-180, 219-243) in the model structure, whereas the cryo-EM structure displayed lower overall drift (except 1-43) in comparison to the model structure.Secondary structural features indicated that a significant unfolding in the transmembrane domains and β-strand at positions 166 to 172, affecting the stability and compactness of the cryo-EM structure , whereas the model exhibited noticeable unfolding in transmembrane domains and small-coiled regions in the N-terminal. , The results from molecular docking and steered molecular dynamics investigations showed the model structure had a greater number of non-bonding interactions, leading to enhanced stability when compared to the cryo-EM structure. Consequently, higher forces were necessary for unbinding of the baricitinib and ruxolitinib inhibitors from the model structure.. Our findings can help better understanding of the significance of unresolved residues at the molecular level. Additionally, this information can guide researchers for experimental endeavors aimed at completing the full-sequence structure of the ORF3a.Communicated by Ramaswamy H. Sarma.

Genome-Wide Identification of the WRKY Gene Family in Four Cotton Varieties and the Positive Role of GhWRKY31 in Response to Salt and Drought Stress

The WRKY gene family is ubiquitously distributed in plants, serving crucial functions in stress responses. Nevertheless, the structural organization and evolutionary dynamics of WRKY genes in cotton have not been fully elucidated. In this study, a total of 112, 119, 217, and 222 WRKY genes were identified in Gossypium arboreum, Gossypium raimondii, Gossypium hirsutum, and Gossypium barbadense, respectively. These 670 WRKY genes were categorized into seven distinct subgroups and unequally distributed across chromosomes. Examination of conserved motifs, domains, cis-acting elements, and gene architecture collectively highlighted the evolutionary conservation and divergence within the WRKY gene family in cotton. Analysis of synteny and collinearity further confirmed instances of expansion, duplication, and loss events among WRKY genes during cotton evolution. Furthermore, GhWRKY31 transgenic Arabidopsis exhibited heightened germination rates and longer root lengths under drought and salt stress. Silencing GhWRKY31 in cotton led to reduced levels of ABA, proline, POD, and SOD, along with downregulated expression of stress-responsive genes. Yeast one-hybrid and molecular docking assays confirmed the binding capacity of GhWRKY31 to the W box of GhABF1, GhDREB2, and GhRD29. The findings collectively offer a systematic and comprehensive insight into the evolutionary patterns of cotton WRKYs, proposing a suitable regulatory framework for developing cotton cultivars with enhanced resilience to drought and salinity stress.

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.

Studying Conformational Properties of Transmembrane Domain of KCNE3 in a Lipid Bilayer Membrane Using Molecular Dynamics Simulations

KCNE3 is a single-pass integral membrane protein that regulates numerous voltage-gated potassium channel functions such as KCNQ1. Previous solution NMR studies suggested a moderate degree of curved α-helical structure in the transmembrane domain (TMD) of KCNE3 in lyso-myristoylphosphatidylcholine (LMPC) micelles and isotropic bicelles with the residues T71, S74 and G78 situated along the concave face of the curved helix. During the interaction of KCNE3 and KCNQ1, KCNE3 pushes its transmembrane domain against KCNQ1 to lock the voltage sensor in its depolarized conformation. A cryo-EM study of KCNE3 complexed with KCNQ1 in nanodiscs suggested a deviation of the KCNE3 structure from its independent structure in isotropic bicelles. Despite the biological significance of KCNE3 TMD, the conformational properties of KCNE3 are poorly understood. Here, all atom molecular dynamics (MD) simulations were utilized to investigate the conformational dynamics of the transmembrane domain of KCNE3 in a lipid bilayer containing a mixture of POPC and POPG lipids (3:1). Further, the effect of the interaction impairing mutations (V72A, I76A and F68A) on the conformational properties of the KCNE3 TMD in lipid bilayers was investigated. Our MD simulation results suggest that the KCNE3 TMD adopts a nearly linear α helical structural conformation in POPC-POPG lipid bilayers. Additionally, the results showed no significant change in the nearly linear α-helical conformation of KCNE3 TMD in the presence of interaction impairing mutations within the sampled time frame. The KCNE3 TMD is more stable with lower flexibility in comparison to the N-terminal and C-terminal of KCNE3 in lipid bilayers. The overall conformational flexibility of KCNE3 also varies in the presence of the interaction-impairing mutations. The MD simulation data further suggest that the membrane bilayer width is similar for wild-type KCNE3 and KCNE3 containing mutations. The Z-distance measurement data revealed that the TMD residue site A69 is close to the lipid bilayer center, and residue sites S57 and S82 are close to the surfaces of the lipid bilayer membrane for wild-type KCNE3 and KCNE3 containing interaction-impairing mutations. These results agree with earlier KCNE3 biophysical studies. The results of these MD simulations will provide complementary data to the experimental outcomes of KCNE3 to help understand its conformational dynamic properties in a more native lipid bilayer environment.

Review: Structural-functional relationship of WRKY transcription factors: Unfolding the role of WRKY in plants

WRKY as the name suggests, are the transcription factors (TFs) that contain the signature WRKY domains, hence named after it. Since their discovery in 1994, they have been well studied in plants with exploration of approximately 74 WRKY genes in the model plant, Arabidopsis alone. However, the study of these transcription factors (TFs) is not just limited to model plant now. They have been studied widely in crop plants as well, because of their tremendous contribution in stress as well as in growth and development. Here, in this review, we describe the story of WRKY TFs from their identification to their origin, the binding mechanisms, structure and their contribution in regulating plant development and stress physiology. High throughput transcriptomics-based data also opened a doorway to understand the comprehensive and detailed functioning of WRKY TFs in plants. Indeed, the detailed functional role of each and every WRKY member in regulating the gene expression is required to pave the path to develop holistic understanding of their role in stress physiology and developmental processes in plants.

Evaluation of mutations on O6-methylguanine methyl transferase structure and its interactions: molecular dynamics simulation study

O6-methylguanine DNA methyl transferase (MGMT) is a significant vehicle for the cellular clearance of alkyl lesions, particularly the methyl group of the O-6 and O-4 positions of guanine and thymine, respectively. Many publications have studied the correlation between polymorphisms in MGMT and susceptibility to various cancers. In the present study, we investigated the consequence of L84F, common single-nucleotide polymorphism, K125E, site-specific mutagenesis, and L84F/K125E on conformation, stability, and behavior of MGMT in the free form and interaction with proliferating cell nuclear antigen (PCNA) and DNA as partners in the biochemical network by using molecular dynamics simulation method. Our results showed that all free variants of MGMT differed from the native form. However, among all free variants of MGMT, the L84F/K125E variant exhibited similar properties compared with the wild-type. In contrast, in complex modes, only amino acid residues of the L84F variant are involved in the interactions with PCNA and DNA somewhat differently relative to the wild-type. Furthermore, L84F SNP showed the highest binding free energy compared to other variants and native forms. These alterations in the amino acids and binding free energy of L84F relative to the native are the reasons for changing its region connection compared to the native form. Therefore, we propose conducting further investigations into the impact of inhibitors or chemotherapeutic agents to assess their effectiveness on MGMT variants compared to the wild-type, aiming to reduce the cost of cancer treatment that will depend on inhibiting native MGMT protein.Communicated by Ramaswamy H. Sarma.

Multifaceted roles of WRKY transcription factors in abiotic stress and flavonoid biosynthesis

Increasing biotic and abiotic stresses are seriously impeding the growth and yield of staple crops and threatening global food security. As one of the largest classes of regulators in vascular plants, WRKY transcription factors play critical roles governing flavonoid biosynthesis during stress responses. By binding major W-box cis-elements (TGACCA/T) in target promoters, WRKYs modulate diverse signaling pathways. In this review, we optimized existing WRKY phylogenetic trees by incorporating additional plant species with WRKY proteins implicated in stress tolerance and flavonoid regulation. Based on the improved frameworks and documented results, we aim to deduce unifying themes of distinct WRKY subfamilies governing specific stress responses and flavonoid metabolism. These analyses will generate experimentally testable hypotheses regarding the putative functions of uncharacterized WRKY homologs in tuning flavonoid accumulation to enhance stress resilience.

CSC01 shows promise as a potential inhibitor of the oncogenic G13D mutant of KRAS: an in silico approach

About 30% of malignant tumors include KRAS mutations, which are frequently required for the development and maintenance of malignancies. KRAS is now a top-priority cancer target as a result. After years of research, it is now understood that the oncogenic KRAS-G12C can be targeted. However, many other forms, such as the G13D mutant, are yet to be addressed. Here, we used a receptor-based pharmacophore modeling technique to generate potential inhibitors of the KRAS-G13D oncogenic mutant. Using a comprehensive virtual screening workflow model, top hits were selected, out of which CSC01 was identified as a promising inhibitor of the oncogenic KRAS mutant (G13D). The stability of CSC01 upon binding the switch II pocket was evaluated through an exhaustive molecular dynamics simulation study. The several post-simulation analyses conducted suggest that CSC01 formed a stable complex with KRAS-G13D. CSC01, through a dynamic protein-ligand interaction profiling analysis, was also shown to maintain strong interactions with the mutated aspartic acid residue throughout the simulation. Although binding free energy analysis through the umbrella sampling approach suggested that the affinity of CSC01 with the switch II pocket of KRAS-G13D is moderate, our DFT analysis showed that the stable interaction of the compound might be facilitated by the existence of favorable molecular electrostatic potentials. Furthermore, based on ADMET predictions, CSC01 demonstrated a satisfactory drug likeness and toxicity profile, making it an exemplary candidate for consideration as a potential KRAS-G13D inhibitor.

Deciphering mitogen activated protein kinase pathway activated during insect attack in Nicotiana attenuata

Plant yields are compromised due to abiotic and biotic stresses. A crucial biotic stress instigated by insect attack, is a major concern that limits crop production. To overcome the deleterious effect of herbivory, pesticides are used but long-term usage of pesticides can be harmful to the environment and human health. Understanding the plants' inherent defense mechanism by interpreting the interaction pattern of defense-related proteins and signalling components and manipulating them to strengthen defense status, is one of the alternative approaches of green biotechnology. During insect attack, host plants initiate innumerable signalling pathways to activate defense response; Mitogen Activated Protein Kinase (MAPK) Pathway is a crucial component of signalling pathway that regulate the expression of downstream defense-related genes. MAPK pathway has three components: MAPKKK, MAPKK and MAPK. Earlier studies have shown participation of SIPK and WIPK (MAPKs) as well as MEK2 (MAPKK) during insect infestation and its association with plant defense. However, information on the third component and elucidation of the complete MAPK pathway are still elusive. Therefore, this study aims to identify the unknown component and decipher MAPK pathway in Nicotiana attenuata involved in plant defense against herbivory by identifying herbivory-inducible MAPKKKs and and their interaction with known partners of the MAPK pathway by docking and MD simulation. The possible pathway was predicted to be MAPKKK Na12134/Na04522-MEK2-SIPK/WIPK. Further, validation of the above interaction by in vitro and in vivo methods is highly recommended.Communicated by Ramaswamy H. Sarma.

Binding of stress-responsive OsWRKY proteins through WRKYGQK heptapeptide residue with the promoter region of two rice blast disease resistance genes Pi2 and Pi54 is important for development of blast resistance

Molecular docking was done to investigate the interactions between five differentially expressed rice WRKY proteins when challenged with the rice blast disease caused by Magnaporthe oryzae and drought stresses applied either individually or overlapped, with the promoter region of two blast resistance genes (Pi2 and Pi54). Molecular docking was performed using the HDOCK server. Initially, the homology models for each of the five rice WRKY proteins were prepared using I-TASSER server, and then the secondary structure as well as the DNA-binding pockets were predicted using PSIPRED and BindUP servers, respectively. The molecular docking study revealed a differential binding pattern of the rice WRKYs with the two blast resistance genes. The WRKY proteins (OsWRKY88 and OsWRKY102), whose transcript levels decrease when drought and blast stresses are overlapped, interact with the two resistance genes mostly involving the residues of the zinc finger structure. On the other hand, the WRKY proteins (OsWRKY53-1 and OsWRKY113), whose transcript levels did not reduce significantly when challenged by drought and blast overlapped condition compared to individual treatment of blast, interact mostly involving the residues of the conserved WRKYGQK heptapeptide sequence. Interestingly, the protein OsWRKY74 whose transcript levels are unaffected in both individual and overlapped stresses, interacts with both the blast resistance genes involving few residues of both WRKYGQK heptapeptide and the zinc finger structure. The findings thus indicate that the interaction of OsWRKY proteins involving the conserved WRKYGQK heptapeptide sequence with the blast resistance genes Pi2 and Pi54 is important to mitigate the blast challenge in rice even during overlapping challenges of drought.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03711-y.

Rice WRKY13 TF protein binds to motifs in the promoter region to regulate downstream disease resistance-related genes

In plants, pathogen resistance is brought about by the binding of certain transcription factor (TF) proteins to the cis-elements of certain target genes. These cis-elements are present upstream in the motif of the promoters of each gene. This ensures the binding of a specific TF to a specific promoter, therefore regulating the expression of that gene. Therefore, the study of each promoter sequence of all the rice genes would help identify the target genes of a specific TF. Rice 1 kb upstream promoter sequences of 55,986 annotated genes were analyzed using the Perl program algorithm to detect WRKY13 binding motifs (bm). The resulting genes were grouped using Gene Ontology and gene set enrichment analysis. A gene with more than 4 TF bm in their promoter was selected. Ten genes reported to have a role in rice disease resistance were selected for further analysis. Cis-acting regulatory element analysis was carried out to find the cis-elements and confirm the presence of the corresponding motifs in the promoter sequences of these genes. The 3D structure of WRKY13 TF and the corresponding ten genes were built, and the interacting residues were determined. The binding capacity of WRKY13 to the promoter of these selected genes was analyzed using docking studies. WRKY13 was considered for docking analysis based on the prior reports of autoregulation. Molecular dynamic simulations provided more details regarding the interactions. Expression data revealed the expression of the genes that helped provide the mechanism of interaction. Further co-expression network helped to characterize the interaction of these selected disease resistance-related genes with the WRKY13 TF protein. This study suggests downstream target genes that are regulated by the WRKY13 TF. The molecular mechanism involving the gene network regulated by WRKY13 TF in disease resistance against rice fungal pathogens is explored.

Molecular dynamics modelling of the interaction of a synthetic zinc-finger miniprotein with DNA

We report the modelling of the DNA complex of an artificial miniprotein composed of two zinc finger modules and an AT-hook linking peptide. The computational study provides for the first time a structural view of these types of complexes, dissecting interactions that are key to modulate their stability. The relevance of these interactions was validated experimentally. These results confirm the potential of this type of computational approach for studying peptide-DNA complexes and suggest that they could be very useful for the rational design of non-natural, DNA binding miniproteins.

Application of molecular dynamics simulation for exploring the roles of plant biomolecules in promoting environmental health

Understanding the dynamic changes of plant biomolecules is vital for exploring their mechanisms in the environment. Molecular dynamics (MD) simulation has been widely used to study structural evolution and corresponding properties of plant biomolecules at the microscopic scale. Here, this review (i) outlines structural properties of plant biomolecules, and the crucial role of MD simulation in advancing studies of the biomolecules; (ii) describes the development of MD simulation in plant biomolecules, determinants of simulation, and analysis parameters; (iii) introduces the applications of MD simulation in plant biomolecules, including the response of the biomolecules to multiple stresses, their roles in corrosive environments, and their contributions in improving environmental health; (iv) reviews techniques integrated with MD simulation, such as molecular biology, quantum mechanics, molecular docking, and machine learning modeling, which bridge gaps in MD simulation. Finally, we make suggestions on determination of force field types, investigation of plant biomolecule mechanisms, and use of MD simulation in combination with other techniques. This review provides comprehensive summaries of the mechanisms of plant biomolecules in the environment revealed by MD simulation and validates it as an applicable tool for bridging gaps between macroscopic and microscopic behavior, providing insights into the wide application of MD simulation in plant biomolecules.

In vitro and in silico characterization of a novel glutamate carboxypeptidase from Cohnella sp. A01

Glutamate carboxypeptidase is a bacterial enzyme of metallopeptidase superfamily. This enzyme is an exo-peptidase that catalyzes the hydrolysis of glutamate residues at the C-terminus of folic acid. The rCP302 is a novel zinc ion-dependent recombinant glutamate carboxypeptidase derived from a thermophilic bacterium, Cohnella sp. A01 (PTCC No: 1921). By simulating the structure of rCP302, analyzing its activity in various environmental settings, and contrasting it with that of related enzymes, we wanted to evaluate the heterologous production, purification, and characterization of this enzyme. The bioinformatics study showed that rCP302 had maximum similarity to M20 family of metallopeptidases. The purified rCP302 molecular weight was about 41.6 kDa. The optimum temperature and pH for the catalytic activity of rCP302 were 50 °C and 7.2, respectively. Fluorescence spectroscopy data elucidated the secondary structure of rCP302 and determined conformational changes caused by alterations in ambient conditions. Using folate as a substrate, K_m and specific activity values were calculated as 0.108 μM and 687 μmol/min/mg, respectively. The enzyme activity was strongly inhibited when EDTA sequestered zinc ions. The half-life of this enzyme at 30 °C was 2012 min. Regarding the ability of rCP302 to degrade folic acid, and its long half-life at 37 °C, the normal temperature of many mammals, this enzyme can be introduced for further study for use in the pharmaceutical industry.

WRKY transcription factors: evolution, regulation, and functional diversity in plants

The recent advancements in sequencing technologies and informatic tools promoted a paradigm shift to decipher the hidden biological mysteries and transformed the biological issues into digital data to express both qualitative and quantitative forms. The transcriptomic approach, in particular, has added new dimensions to the versatile essence of plant genomics through the large and deep transcripts generated in the process. This has enabled the mining of super families from the sequenced plants, both model and non-model, understanding their ancestry, diversity, and evolution. The elucidation of the crystal structure of the WRKY proteins and recent advancement in computational prediction through homology modeling and molecular dynamic simulation has provided an insight into the DNA-protein complex formation, stability, and interaction, thereby giving a new dimension in understanding the WRKY regulation. The present review summarizes the functional aspects of the high volume of sequence data of WRKY transcription factors studied from different species, till date. The review focuses on the dynamics of structural classification and lineage in light of the recent information. Additionally, a comparative analysis approach was incorporated to understand the functions of the identified WRKY transcription factors subjected to abiotic (heat, cold, salinity, senescence, dark, wounding, UV, and carbon starvation) stresses as revealed through various sets of studies on different plant species. The review will be instrumental in understanding the events of evolution and the importance of WRKY TFs under the threat of climate change, considering the new scientific evidences to propose a fresh perspective.

Structural and functional characterisation of HepTH1-5 peptide as a potential hepcidin replacement

Hepcidin is a principal regulator of iron homeostasis and its dysregulation has been recognised as a causative factor in cancers and iron disorders. The strategy of manipulating the presence of hepcidin peptide has been used for cancer treatment. However, this has demonstrated poor efficiency and has been short-lived in patients. Many studies reported using minihepcidin therapy as an alternative way to treat hepcidin dysregulation, but this was only applied to non-cancer patients. Highly conserved fish hepcidin protein, HepTH1-5, was investigated to determine its potential use in developing a hepcidin replacement for human hepcidin (Hepc25) and as a therapeutic agent by targeting the tumour suppressor protein, p53, through structure-function analysis. The authors found that HepTH1-5 is stably bound to ferroportin, compared to Hepc25, by triggering the ferroportin internalisation via Lys42 and Lys270 ubiquitination, in a similar manner to the Hepc25 activity. Moreover, the residues Ile24 and Gly24, along with copper and zinc ligands, interacted with similar residues, Lys24 and Asp1 of Hepc25, respectively, showing that those molecules are crucial to the hepcidin replacement strategy. HepTH1-5 interacts with p53 and activates its function through phosphorylation. This finding shows that HepTH1-5 might be involved in the apoptosis signalling pathway upon a DNA damage response. This study will be very helpful for understanding the mechanism of the hepcidin replacement and providing insights into the HepTH1-5 peptide as a new target for hepcidin and cancer therapeutics.Communicated by Ramaswamy H. Sarma.

pH modulates the role of SP6 RNA polymerase in transcription process: an in silico study

SP6 RNA polymerase (SP6 RNAP) is an essential enzyme for the transcription process in SP6 bacteriophage. SP6 RNAP plays a vital role in mRNA vaccine designing technology and other translational biotechnology research due to the high specificity towards its promoter. The self-replicating performance also put this polymerase to study extensively. Despite of the reports emphasizing the function of this enzyme, a detailed structural and functional understanding of RNA polymerase is not reported so far. Here, we report the first-ever information about SP6RNAP structure and its effect on promoter binding at different pH environments using molecular docking and molecular dynamics simulation (MDS) study. We also report the changes in polymerase conformations in different pH conditions using in-silico approach. The docking study was also performed for SP6 RNAP with SP6 promoter at different pH environments using the in-silico docking tools and conducted the MDS study for complexes. MM/PBSA and per residue energy contribution has been performed at three different pH environments. The structural aspects confirmed that the pH 7.9 state favors the polymerase functional activity in the transcription process which was in the range reported using transcription assay. This polymerase's unique features may play its emerging role as an efficient transcription factor in translational biological research.Communicated by Ramaswamy H. Sarma.

Integrated computational approach towards identification of HSPG and ACE2 mimicking moieties for SARS-CoV-2 inhibition

A key step to inhibit severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is to prevent the entry of the virus into the host cells. The receptor-binding domains (RBDs) of spike proteins of SARS-CoV and other human coronaviruses utilize heparan sulfate proteoglycans (HSPGs) as the primary receptors for their accumulation on the cell surface and then scan for binding to the main entry receptor angiotensin-converting enzyme 2 (ACE2). SARS-CoV and SARS-CoV-2 share structurally similar RBDs and therefore, it is possible that SARS-COV-2 primarily binds to HSPGs followed by binding to the ACE2 receptors. A promising strategy to inhibit virus infection is to circulate exogenous bioactive moieties structurally mimicking cellular HSPG and ACE2 which act as decoy receptors binding to SARS-CoV-2 and competitively inhibit virus entry to the host cells mediated by cellular-bound HSPG and ACE2. Using a molecular docking tool, we identified carboxymethyl benzyl amide sulfonate (CMBS) and polyanetholesulfonic acid (PAS) as the suitable HSPG mimicking ligands, and Paenibacillus sp. B38-derived carboxypeptidase (B38-CAP) and Bacillus subtilis-derived carboxypeptidase (BS-CAP) as the potential ACE2-like enzymes having a strong binding affinity to the spike proteins as that of cellular HSPG and ACE2. Further, the binding stability and compactness of these moieties with SARS-CoV-2 were analyzed through molecular dynamics (MD) simulations, and the results indicated that these moieties form well-stable complexes with the RBD of spike proteins. The identified moieties could be conjugated to the surfaces of non-toxic nanoparticles to provide multiple interactions to efficiently shield SARS-CoV-2, and inhibit viral entry to the host cells.

Liposomal Formulation of Botanical Extracts may Enhance Yield Triggering PR Genes and Phenylpropanoid Pathway in Barley (Hordeum vulgare)

This work aimed to study the plant conditioning effect and mode of action of a plant-based biostimulant used in organic farming. This new generation plant biostimulant, named ELICE16INDURES®, is rich in plant bio-active ingredients containing eleven supercritical botanical extracts encapsulated in nano-scale liposomes. The dose-response (10 to 240 g ha-1) relationship was tested in a field population of autumn barley (Hordeum vulgare) test crop, and underlying molecular mechanisms were studied. Applying nanotechnology, cell-identical nanoparticles may help the better uptake and delivery of active ingredients increasing resilience, vitality, and crop yield. The amount of harvested crops showed a significant increase of 27.5% and 39.9% interconnected to higher normalized difference vegetation index (NDVI) of 20% and 25% after the treatment of low and high dosages (20 and 240 g ha-1), respectively. Illumina NextSeq 550 sequencing, gene expression profiling, and KEGG-pathway analysis of outstanding dosages indicated the upregulation of pathogenesis-related (PR) and other genes-associated with induced resistance-which showed dose dependency as well.

A drug repurposing endeavor to discover a multi-targeting ligand against RhlR and LasR proteins from opportunistic human pathogen Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic human pathogen. It synthesizes the poison called Hydrogen Cyanide (HCN). The synthesis of HCN is mediated by the enzyme HCN synthase which is obtained from the hcnABC operon and the transcription of the hcnABC operon is mediated by three proteins LasR, RhlR, and ANR. In our previous works, we analyzed the activation process of RhlR and LasR proteins by their cognate auto-inducer ligands (N-butanoyl-L-homoserine lactone and N-(3-oxododecanoyl)-homoserine lactone respectively). In this work, we attempted to identify some multi-targeting ligands which would be able to destroy the structural integrity of both the RhlR and LasR proteins using steered MD simulations. We used the virtual screening of ligand libraries, and for that purpose, we used the NCI drug database. We selected the top 4 ligands from our virtual screening experiments. We then tried to check their relative binding affinities with the LasR and RhlR proteins in comparison to their native auto-inducer ligands. Through this work, we were able to identify 4 such ligands which were capable of binding to both the RhlR and LasR proteins in a better way than their native auto-inducer ligands. The efficacies of these ligands to actually perturb the structural integrity of RhlR and LasR proteins could be tested in wet lab. The work is the first work in the field of structure-based drug design to come up with possible multi-targeting drug-like structures against the RhlR and LasR proteins from Pseudomonas aeruginosa.

Binding energy analysis and molecular dynamic simulation studies of the designed orally active, non-toxic GABARAP modulators

Epilepsy is a severe neurological disorder that occurs when the communication between the neurons is disturbed. Gamma-amino butyric acid-associated protein (GABARAP) plays a key role in balancing Gamma-aminobutyric acid-A (GABA(A)) receptor functions of inhibiting the neurotransmission and controlling the seizure. In this study, we introduce the derivatives of the selected anti-epileptic drugs, namely Felbamate and Clobazam, by substituting different hydrophilic and hydrophobic groups at the specified positions. Molecular docking studies between the derivatives and GABARAP were carried out using PyRx software. The interacting residues were identified from LigPlot⁺. Drug-likeness, drug-related properties, and toxic endpoints of each derivative were analyzed using the SwissADME, Osiris property explorer, and ProTox-II servers. After analyzing the binding energy, drug-properties, and toxicity, the best five derivatives of Felbamate and Clobazam were selected. Molecular Dynamic simulation studies involving the target-ligand interaction were carried out for 100 nanoseconds using GROMACS 2018. The root mean square deviation, root mean square fluctuation, radius of gyration, Solvent accessible area, Energy plots and trajectories of the ten GABARAP complexes of the derivatives, and two GABARAP complexes of parent drugs were compared and critically analyzed. Among the five Felbamate derivatives, F7 formed the most stable complex with GABARAP. Among the five Clobazam derivatives, C27, C33 and C32 showed stable GABARAP interaction. In light of the above systematic computational analysis, we propose F7, C27, C33, and C32 as the potential anti-epileptic drug candidates for developing novel therapeutics. The substitution of hydrophobic groups at para position on benzene ring has promoted strong binding to GABARAP. Communicated by Ramaswamy H. Sarma.

In-silico study for the screening and preparation of ionic liquid-AVDs conjugate to combat COVID-19 surge

The pandemic due to COVID-19 caused by SARS-CoV-2 has led to the recorded deaths worldwide and is still a matter of concern for scientists to find an effective counteragent. The combination therapy is always been a successful attempt in treating various threatful diseases. Recently, Ionic liquids (ILs) are known for their antiviral activity. Fascinating tunable properties of ILs make them a potential candidate for designing the therapeutic agent. The concern while using ILs in biomedical field remains is toxicity therefore, choline-based ILs were used in the study as they are considered to be greener as compared to other ILs. In the present study strategically, we performed the blind molecular docking of antiviral drug (Abacavir, Acyclovir, and Galidesivir)-choline based ILs conjugates with the target protein (Mpro protease). The molecules were screened on the basis of binding energy. The data suggested that the combination of AVDs-ILs have greater antiviral potential as compared to the drugs and ILs alone. Further, the ADME properties and toxicity analysis of the screened conjugates was done which revealed the non-toxicity of the conjugates. Additionally, the energetic profiling of the ILs drugs and their conjugates was done using DFT calculations which revealed the stability of the conjugates and have a better option to be developed as a therapeutic agent. Also, from molecular dynamic simulation was done and results showed the stability of the complex formed between target protein and the designed conjugates of AVDs and ILs.

Investigating Structural Dynamics of KCNE3 in Different Membrane Environments Using Molecular Dynamics Simulations

KCNE3 is a potassium channel accessory transmembrane protein that regulates the function of various voltage-gated potassium channels such as KCNQ1. KCNE3 plays an important role in the recycling of potassium ion by binding with KCNQ1. KCNE3 can be found in the small intestine, colon, and in the human heart. Despite its biological significance, there is little information on the structural dynamics of KCNE3 in native-like membrane environments. Molecular dynamics (MD) simulations are a widely used as a tool to study the conformational dynamics and interactions of proteins with lipid membranes. In this study, we have utilized all-atom molecular dynamics simulations to characterize the molecular motions and the interactions of KCNE3 in a bilayer composed of: a mixture of POPC and POPG lipids (3:1), POPC alone, and DMPC alone. Our MD simulation results suggested that the transmembrane domain (TMD) of KCNE3 is less flexible and more stable when compared to the N- and C-termini of KCNE3 in all three membrane environments. The conformational flexibility of N- and C-termini varies across these three lipid environments. The MD simulation results further suggested that the TMD of KCNE3 spans the membrane width, having residue A69 close to the center of the lipid bilayers and residues S57 and S82 close to the lipid bilayer membrane surfaces. These results are consistent with previous biophysical studies of KCNE3. The outcomes of these MD simulations will help design biophysical experiments and complement the experimental data obtained on KCNE3 to obtain a more detailed understanding of its structural dynamics in the native membrane environment.

A tau class GST, OsGSTU5, interacts with VirE2 and modulates the Agrobacterium-mediated transformation in rice

OsGSTU5 interacts and glutathionylates the VirE2 protein of Agrobacterium and its (OsGSTU5) overexpression and downregulation showed a low and high AMT efficiency in rice, respectively. During Agrobacterium-mediated transformation (AMT), T-DNA along with several virulence proteins such as VirD2, VirE2, VirE3, VirD5, and VirF enter the plant cytoplasm. VirE2 serves as a single-stranded DNA binding (SSB) protein that assists the cytoplasmic trafficking of T-DNA inside the host cell. Though the regulatory roles of VirE2 have been established, the cellular reaction of their host, especially in monocots, has not been characterized in detail. This study identified a cellular interactor of VirE2 from the cDNA library of rice. The identified plant protein encoded by the gene cloned from rice was designated OsGSTU5, it interacted specifically with VirE2 in the host cytoplasm. OsGSTU5 was upregulated during Agrobacterium infection and involved in the post-translational glutathionylation of VirE2 (gVirE2). Interestingly, the in silico analysis showed that the 'gVirE2 + ssDNA' complex was structurally less stable than the 'VirE2 + ssDNA' complex. The gel shift assay also confirmed the attenuated SSB property of gVirE2 over VirE2. Moreover, knock-down and overexpression of OsGSTU5 in rice showed increased and decreased T-DNA expression, respectively after Agrobacterium infection. The present finding establishes the role of OsGSTU5 as an important target for modulation of AMT efficiency in rice.

The Efficacy Mechanism of Epigallocatechin Gallate against Pre-Eclampsia based on Network Pharmacology and Molecular Docking

Pre-eclampsia (PE), a pregnancy complication, affects 3-5% of all pregnancies worldwide and is the main cause of maternal and perinatal morbidity. However, there is no drug which can clearly slow this disease progression. Epigallocatechin gallate (EGCG), a natural compound extracted from green tea, has been found to enhance the treatment efficacy of oral nifedipine against pregnancy-induced severe PE. This study aims to clarify the potential targets and pharmacological mechanisms of EGCG in treatment of PE. We used Traditional Chinese Medicine Systems Pharmacology database and Gene Cards database to obtain 179 putative target proteins of EGCG, 550 PE-related hub genes and 39 intersecting targets between EGCG and PE. By using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses, we got the gene entries and enrichment pathways closely related to the intersecting targets. The top 10 enrichment pathways were pathway in cancer, proteoglycans in cancer, HIF-1 signaling pathway, AGE-RAGE signaling pathway in diabetic complications, TNF signaling pathway, bladder cancer, hepatitis B, IL-17 signaling pathway, toxoplasmosis, PI3K-Akt signaling pathway. Furthermore, compound-target-pathway (CTP) and protein-protein interaction (PPI) network analysis were employed to explore the interaction of the top twelve targets for EGCG in treating PE. Molecular docking analysis showed combinations between these targets and EGCG, and the interaction between EGCG and the targets IL-6 and EGFR was confirmed by using molecular dynamic simulation. In conclusion, these findings hint the underlying mechanism of EGCG in the treatment of PE and point out directions in further studies on PE.

In silico molecular docking study of Andrographis paniculata phytochemicals against TNF-α as a potent anti-rheumatoid drug

Tumor necrosis factor-α (TNF-α) is a proinflammatory cytokine which plays a crucial role in controlling inflammatory responses. The pathway of Rheumatoid arthritis (RA) leading to TNF-alpha is activated by macrophages and quite often by natural killer cells and lymphocytes. In the inflammatory phase, it is believed to be the main mediator and to be anchored with the progression of different diseases such as ankylosing spondylitis, Crohn's disease, and Rheumatoid arthritis (RA). The major goal of this study is to use in silico docking studies to investigate the anti-inflammatory potential of a bioactive molecule from the medicinal plant Andrographis paniculata. The three-dimensional structures of different phytochemicals of A. paniculata were obtained from PubChem database, and the receptor protein was derived from PDB database. Docking analysis was executed using AutoDock vina, and the binding energies were compared. Bisandrographolide A and Andrographidine C revealed the highest score of -8.6 Kcal/mol, followed by, Neoandrographolide (-8.5 Kcal/mol). ADME and toxicity parameters were evaluated for these high scoring ligands and results showed that Andrographidine C could be a potent drug, whereas Neoandrographolide and Bisandrographolide A can be modified in in vitro and can lead to a promising drug. Further, the top scorer (Andrographidine C) and control drug (Leflunomide) were subjected to 100 ns MD Simulation. The protein complex with Andrographidine C had more stable confirmation with lower RMSD (0.28 nm) and higher binding energy (-133.927 +/- 13.866 kJ/mol). In conclusion, Andrographidine C may be a potent surrogate to the disease-modifying anti-rheumatic drugs (DMARD's) & Non-steroidal anti-inflammatory drugs (NSAID's) that has fewer or minor adverse effects and can aid in RA management.

Comparative molecular docking and molecular-dynamic simulation of wild-type- and mutant carboxylesterase with BTA-hydrolase for enhanced binding to plastic

According to the literature review, microbial degradation of polyethylene terephthalate by PETases has been detected effective and eco-friendly. However, the number of microorganisms capable of such feats is limited with some undesirable bioprospecting results. BTA-hydrolase has been already reported capable of degrading polyethylene terephthalate. Therefore, mutation by in silico site-directed mutagenesis means to introduce current isomer of PETase for polyethylene terephthalate degradative capability as a better approach to resolve this issue. This study aimed to use in silico site-directed mutagenesis to convert a carboxylesterase from Archaeoglobus fulgidus to BTA-hydrolase from Thermobifida fusca by replacing six amino acids in specific locations. This work was followed by molecular docking analysis with polyethylene terephthalate and polypropylene to compare their interactions. The best-docked enzyme-substrate complex was further subjected to molecular dynamics simulation to gauge the binding quality of the BTA-hydrolase, wild-type and mutant-carboxylesterase with only polyethylene terephthalate as a substrate. Results of molecular docking revealed lowest binding energy for the wild-type carboxylesterase-polypropylene complex (-7.5 kcal/mol). The root-mean-square deviation value was observed stable for BTA-hydrolase. Meanwhile, root-mean-square fluctuation was assessed with higher fluctuation for the mutated residue Lys178. Consequently, the Rg value for BTA-hydrolase-ligand complex (∼1.68 nm) was the lowest compared to the mutant and wild-type carboxylesterase. The collective data conveyed that mutations imparted a minimal change in the ability of the mutant carboxylesterase to bind to polyethylene terephthalate.

Predicting Diagnostic Potential of Cathepsin in Epithelial Ovarian Cancer: A Design Validated by Computational, Biophysical and Electrochemical Data

Background: Epithelial ovarian cancer remains one of the leading variants of gynecological cancer with a high mortality rate. Feasibility and technical competence for screening and detection of epithelial ovarian cancer remain a major obstacle and the development of point of care diagnostics (POCD) may offer a simple solution for monitoring its progression. Cathepsins have been implicated as biomarkers for cancer progression and metastasis; being a protease, it has an inherent tendency to interact with Cystatin C, a cysteine protease inhibitor. This interaction was assessed for designing a POCD module. Methods: A combinatorial approach encompassing computational, biophysical and electron-transfer kinetics has been used to assess this protease-inhibitor interaction. Results: Calculations predicted two cathepsin candidates, Cathepsin K and Cathepsin L based on their binding energies and structural alignment and both predictions were confirmed experimentally. Differential pulse voltammetry was used to verify the potency of Cathepsin K and Cathepsin L interaction with Cystatin C and assess the selectivity and sensitivity of their electrochemical interactions. Electrochemical measurements indicated selectivity for both the ligands, but with increasing concentrations, there was a marked difference in the sensitivity of the detection. Conclusions: This work validated the utility of dry-lab integration in the wet-lab technique to generate leads for the design of electrochemical diagnostics for epithelial ovarian cancer.

Hydration modulates oxygen channel residues for oxygenation of cysteine dioxygenase: Perspectives from molecular dynamics simulations

Cysteine dioxygenase (CDO) regulates the concentration of l-cysteine substrate by its oxidation in the body to prevent different diseases, including neurodegenerative and autoimmune diseases. CDO catalyzes the oxidation of thiol group of l-cysteine to l-cysteine sulfinic acid using molecular oxygen. In this study, molecular dynamics simulations were applied to ligand-free CDO, cysteine-bound CDO, and oxygen-bound CDO-cysteine complex which were primarily subjected to the evaluation of their structural and dynamical properties. The simulation data provided significant information not only on the conformational changes of the enzyme after its ligation but also on the co-ligation by sequential binding of l-cysteine and molecular oxygen. It was found that the ligation and co-ligation perturbed the active site region as well as the overall protein dynamics which were analyzed in terms of root mean square deviation, root mean square fluctuation and dynamic cross correlation matrices as well as principal component analysis. Furthermore, oxygen transport pathways were successfully explored by taking various tunnel clusters into account and one of those clusters was given preference based on the throughput value. The bottleneck formed by different amino acid residues was examined to figure out their role in the oxygenation process of the enzyme. The residues forming the tunnel's bottleneck and their dynamics mediated by water molecules were further investigated using radial distribution functions which gave insights into the hydration behavior of these residues. The findings based on the hydration behavior in turn served to explore the water-mediated dynamics of these residues in the modulation of the pathway, including tunnel gating for the oxygen entry and diffusion to the active site, which is essential for the CDO's catalytic function.

Thermodynamic stability, in-vitro permeability, and in-silico molecular modeling of the optimal Elaeis guineensis leaves extract water-in-oil nanoemulsion

Nanoemulsion is a delivery system used to enhance bioavailability of plant-based compounds across the stratum corneum. Elaeis guineensis leaves are rich source of polyphenolic antioxidants, viz. gallic acid and catechin. The optimal E. guineensis leaves extract water-in-oil nanoemulsion was stable against coalescence, but it was under significant influence of Ostwald ripening over 90 days at 25 °C. The in-vitro permeability revealed a controlled and sustained release of the total phenolic compounds (TPC) of EgLE with a cumulative amount of 1935.0 ± 45.7 µgcm-2 after 8 h. The steady-state flux and permeation coefficient values were 241.9 ± 5.7 µgcm-2 h-1 and 1.15 ± 0.03 cm.h-1, respectively. The kinetic release mechanism for TPC of EgLE was best described by the Korsmeyer-Peppas model due to the highest linearity of R2 = 0.9961, indicating super case II transport mechanism. The in-silico molecular modelling predicted that the aquaporin-3 protein in the stratum corneum bonded preferably to catechin over gallic acid through hydrogen bonds due to the lowest binding energies of - 57.514 kcal/mol and - 8.553 kcal/mol, respectively. Thus, the in-silico study further verified that catechin could improve skin hydration. Therefore, the optimal nanoemulsion could be used topically as moisturizer to enhance skin hydration based on the in-silico prediction.

Identification of Novel Mutations in Arabidopsis thaliana DOF 4.2 Coding Gene

Purpose: DOF (DNA-binding with One Finger) proteins are plant-specific transcription factors which mediate numerous biological processes. The purpose of the current study is to report new naturally occurring mutations in the gene encoding for one of the members of DOF proteins named DOF 4.2.

Methods: The expression of zinc finger domain of DOF 4.2 (DOF 4.2-ZF) was investigated by first synthesis of cDNA library using different parts of Arabidopsis thaliana plant. Then the coding sequence for zinc finger domain of DOF 4.2 protein was prepared using nested PCR experiment and cloned into pGEX-6P-1 expression vector. Finally, the prepared construct was used for protein expression. Furthermore, molecular dynamics (MD) simulation was carried out to predict DNA binding affinity of DOF 4.2-ZF using AMBER package.

Results: For the first time a new variant of DOF 4.2-ZF protein with three mutations was detected. One of the mutations is silent while the other two mutations lead to amino acid replacement (S18G) as well as introduction of a stop codon ultimately resulting in a truncated protein production. In order to investigate whether the truncated form is able to recognize DNA binding motif, MD simulations were carried out and the results showed that the chance of binding of DOF 4.2-ZF protein to cognate DNA in its truncated form is very low.

Conclusion: The findings demonstrated that the observed mutations adversely affect the DNA binding ability of the truncated form of DOF4.2 if it is expressed in the mutant variant of A. thaliana used in this study.

A novel variant fibrinogen, AαE11del, demonstrating the importance of AαE11 residue in thrombin binding

INTRODUCTION

We identified a novel heterozygous AαE11del variant in a patient with congenital dysfibrinogenemia. This mutation is located in fibrinopeptide A (FpA). We analyzed the effect of AαE11del on the catalyzation of thrombin and batroxobin and simulated the stability of the complex structure between the FpA fragment (AαG6-V20) peptide and thrombin.

MATERIALS AND METHODS

We performed fibrin polymerization and examined the kinetics of FpA release catalyzed by thrombin and batroxobin using purified plasma fibrinogen. To clarify the association between the AαE11 residue and thrombin, we calculated binding free energy using molecular dynamics simulation trajectories.

RESULTS

Increasing the thrombin concentration improved release of FpA from the patient's fibrinogen to approximately 90%, compared to the previous 50% of that of normal fibrinogen. Fibrin polymerization of variant fibrinogen also improved. In addition, greater impairment of variant FpA release from the patient's fibrinogen was observed with thrombin than with batroxobin. Moreover, the calculated binding free energy showed that the FpA fragment-thrombin complex became unstable due to the missing AαE11 residue.

CONCLUSIONS

Our findings indicate that the AαE11 residue is involved in FpA release in thrombin catalyzation more than in batroxobin catalyzation, and that the AαE11 residue stabilizes FpA fragment-thrombin complex formation.

Probe into the Molecular Mechanism of Ibuprofen Interaction with Warfarin Bound to Human Serum Albumin in Comparison to Ascorbic and Salicylic Acids: Allosteric Inhibition of Anticoagulant Release

The release of anticoagulant drugs such as warfarin from human serum albumin (HSA) has been important not only mechanistically but also clinically for patients who take multiple drugs simultaneously. In this study, the role of some commonly used drugs, including s-ibuprofen, ascorbic acid, and salicylic acid, was investigated in the release of warfarin bound to HSA in silico. The effects of the aforementioned drugs on the HSA-warfarin complex were investigated with molecular dynamics (MD) simulations using two approaches; in the first perspective, molecular docking was used to model the interaction of each drug with the HSA-warfarin complex, and in the second approach, drugs were positioned randomly and distant from the binary complex (HSA-warfarin) in a physiologically relevant concentration. The results obtained from both approaches indicated that s-ibuprofen and ascorbic acid both displayed allosteric effects on the release of warfarin from HSA. Although ascorbic acid aided in warfarin release, leading to destabilization of HSA, ibuprofen demonstrated a stabilizing effect on releasing the anticoagulant drug through several noncovalent interactions, including hydrophobic, electrostatic, and hydrogen-bonding interactions with the protein. The calculated binding free energy and energy contribution of involved residues using the molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) method, along with root mean square deviation (RMSD) values, protein gyration, and free energy surface (FES) mapping of the protein, provided valuable details on the nature of the interactions of each drug on the release of warfarin from HSA. These results can provide important information on the mechanisms of anticoagulant release that has not been revealed in molecular details previously.

Revealing atomic-scale molecular diffusion of a plant-transcription factor WRKY domain protein along DNA

In transcription factors’ search for target genes, one-dimensional diffusion of the protein along DNA is essential. Experimentally, it remains challenging to resolve the individual diffusional steps of protein on DNA. Here, we report mainly all-atom equilibrium simulations of a WRKY domain protein in association with and diffusion along DNA. We demonstrate a complete stepping cycle of the protein for one base pair on DNA within microseconds, along with stochastic motions. Processive protein diffusions on DNA have been further sampled in a coarse-grained model. We have also found preferential DNA-strand association of the domain protein, which becomes most prominent at specific DNA binding, and it can be common for small-domain proteins to balance movements on the DNA with the sequence recognition.

Transcription factor (TF) target search on genome is highly essential for gene expression and regulation. High-resolution determination of TF diffusion along DNA remains technically challenging. Here, we constructed a TF model system using the plant WRKY domain protein in complex with DNA from crystallography and demonstrated microsecond diffusion dynamics of WRKY on DNA by employing all-atom molecular-dynamics (MD) simulations. Notably, we found that WRKY preferentially binds to one strand of DNA with significant energetic bias compared with the other, or nonpreferred strand. The preferential DNA-strand binding becomes most prominent in the static process, from nonspecific to specific DNA binding, but less distinct during diffusive movements of the domain protein on the DNA. Remarkably, without employing acceleration forces or bias, we captured a complete one-base-pair stepping cycle of the protein tracking along major groove of DNA with a homogeneous poly-adenosine sequence, as individual hydrogen bonds break and reform at the protein–DNA binding interface. Further DNA-groove tracking motions of the protein forward or backward, with occasional sliding as well as strand crossing to minor groove of DNA, were also captured. The processive diffusion of WRKY along DNA has been further sampled via coarse-grained MD simulations. The study thus provides structural dynamics details on diffusion of a small TF domain protein, suggests how the protein approaches a specific recognition site on DNA, and supports further high-precision experimental detection. The stochastic movements revealed in the TF diffusion also provide general clues about how other protein walkers step and slide along DNA.

Molecular Mechanism of the Intercalation of the SOX-4 Protein into DNA Inducing Bends and Kinks

DNA-protein interactions regulate several biophysical functions, yet the mechanism of only a few is investigated in molecular detail. An important example is the intercalation of transcription factor proteins into DNA that produce bent and kinked DNA. Here, we have studied the molecular mechanism of the intercalation of a transcription factor SOX4 into DNA with a goal to understand the sequence of molecular events that precede the bending and kinking of the DNA. Our long well-tempered metadynamics and molecular dynamics (MD) simulations show that the protein primarily binds to the backbone of DNA and rotates around it to form an intercalative native state. We show that although there are multiple pathways for intercalation, the deintercalation pathway matches with the most probable intercalation pathway. In both cases, bending and kinking happen simultaneously, driven by the onset of the intercalation of the amino acid.

Evolution and identification of DREB transcription factors in the wheat genome: modeling, docking and simulation of DREB proteins associated with salt stress

Soil salinity and the resulting salt stress it imposes on crop plants is a major problem for modern agriculture. Understanding how salt tolerance mechanisms in plants are regulated is therefore important. One regulatory mechanism is the APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription factor family, including dehydration responsive element binding (DREB) transcription factors. By binding to DNA, specifically upstream of genes that play roles in salt tolerance pathways, DREB proteins upregulate expression of these genes. DREB in Triticum aestivum (wheat) cluster in sub-groups and in this study by scanning the recently extended predicted proteome of wheat for DREB, we increased the number of members of this sub-family. Using the wheat genome, we identified 576 genes coding for the AP2 domain of which 508 were identified to have one AP2 domain, a characteristic of the DREB/ERF subfamily. We confirmed the existing four sub-groups by sequence-based phylogenetic analyses but also identified 32 new DREB subfamily members, not belonging to any known sub-group. Transcription factor profile inference analysis identified two genes, TraesCS2B02G002700 and TraesCS2D02G015200, being homologous to DREB1A of Arabidopsis thaliana. Based on molecular simulation (25 ns) analysis, TraesCS2B02G002700 with a CCGAC motif was observed to interact very stably with DNA. In silico mutational analysis at the 19^th position in the DREB domain of TraesCS2B02G002700-DNA complex indicated this as a stable part for recognizing and forming interaction with DNA. Moreover, six target genes were predicted having an upstream CCGAC motif regulated by TraesCS2B02G002700. Our study provides an overall framework for exploring the transcription factors in plants and identifying e.g. potential salt stress target genes.Communicated by Ramaswamy H. Sarma.

Highlighting the potential utility of MBP crystallization chaperone for Arabidopsis BIL1/BZR1 transcription factor-DNA complex

The maltose-binding protein (MBP) fusion tag is one of the most commonly utilized crystallization chaperones for proteins of interest. Recently, this MBP-mediated crystallization technique was adapted to Arabidopsis thaliana (At) BRZ-INSENSITIVE-LONG (BIL1)/BRASSINAZOLE-RESISTANT (BZR1), a member of the plant-specific BZR TFs, and revealed the first structure of AtBIL1/BZR1 in complex with target DNA. However, it is unclear how the fused MBP affects the structural features of the AtBIL1/BZR1-DNA complex. In the present study, we highlight the potential utility of the MBP crystallization chaperone by comparing it with the crystallization of unfused AtBIL1/BZR1 in complex with DNA. Furthermore, we assessed the validity of the MBP-fused AtBIL1/BZR1-DNA structure by performing detailed dissection of crystal packings and molecular dynamics (MD) simulations with the removal of the MBP chaperone. Our MD simulations define the structural basis underlying the AtBIL1/BZR1-DNA assembly and DNA binding specificity by AtBIL1/BZR1. The methodology employed in this study, the combination of MBP-mediated crystallization and MD simulation, demonstrates promising capabilities in deciphering the protein-DNA recognition code.

Elucidation of the hetero-dimeric binding activity of LasR and RhlR proteins with the promoter DNA and the role of a specific Phe residue during the biosynthesis of HCN synthase from opportunistic pathogen Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic human pathogen. It causes secondary infections in patients suffering from cancer and other immunological disorders. The pathogenicity of the organism is dependent on the ability of the organism to code for hydrogen cyanide (HCN), the synthesis of which is mediated by HCN synthase enzyme. HCN synthase is encoded by hcnABC operon. The transcription of the operon is controlled by a complex interplay between the proteins LasR and RhlR. Till date, there is no report that deals with the binding interactions of the RhlR-LasR heterodimer with the promoter DNA region of the hcnABC operon. We, for the first time, tried to analyse the binding modes of the RhlR-LasR heterodimer with the promoter DNA regions. From our work, we could predict the importance of a specific amino acid residue Phe214 from RhlR which might be considered to have the desired specificity to bind to the promoter DNA. Therefore, the amino acid Phe214 may be targeted to develop suitable ligands to eradicate the spread of secondary infections by Pseudomonas aeruginosa.

In silico profiling and structural insights of zinc metal ion on O6-methylguanine methyl transferase and its interactions using molecular dynamics approach

O6-methylguanine DNA methyl transferase (MGMT) is a metalloenzyme participating in the repair of alkylated DNA. In this research, we performed a comparative study for evaluating the impact of zinc metal ion on the behavior and interactions of MGMT in the both enzymatic forms of apo MGMT and holo MGMT. DNA and proliferating cell nuclear antigen (PCNA), as partners of MGMT, were utilized to evaluate molecular interactions by virtual microscopy of molecular dynamics simulation. The stability and conformational alterations of each forms (apo and holo) MGMT-PCNA, and (apo and holo) MGMT-DNA complexes were calculated by MM/PBSA method. A total of seven systems including apo MGMT, holo MGMT, free PCNA, apo MGMT-PCNA, holo MGMT-PCNA, apo MGMT-DNA, and holo MGMT-DNA complexes were simulated. In this study, we found that holo MGMT was more stable and had better folding and functional properties than that of apo MGMT. Simulation analysis of (apo and holo) MGMT-PCNA complexes displayed that the sequences of the amino acids involved in the interactions were different in the two forms of MGMT. The important amino acids of holo MGMT involved in its interaction with PCNA included E92, K101, A119, G122, N123, P124, and K125, whereas the important amino acids of apo MGMT included R128, R135, S152, N157, Y158, and L162. Virtual microscopy of molecular dynamics simulation showed that the R128 and its surrounding residues were important amino acids involved in the interaction of holo MGMT with DNA that was exactly consistent with X-ray crystallography structure. In the apo form of the protein, the N157 and its surrounding residues were important amino acids involved in the interaction with DNA. The binding free energies of - 387.976, - 396.226, - 622.227, and - 617.333 kcal/mol were obtained for holo MGMT-PCNA, apo MGMT-PCNA, holo MGMT-DNA, and apo MGMT-DNA complexes, respectively. The principle result of this research was that the area of molecular interactions differed between the two states of MGMT. Therefore, in investigations of metalloproteins, the metal ion must be preserved in their structures. Finally, it is recommended to use the holo form of metalloproteins in in vitro and in silico researches.

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 9^th and 39^th positions of RsAFP2 and their increased membrane deformation activity.

An in silico approach to identify potential medicinal plants for treating Alzheimer disease: a case study with acetylcholinesterase

Alzheimer's disease (AD) is a progressive neurological disorder affecting an estimated 10 million people worldwide. There is no cure for AD, and only a handful of drugs are known to provide some relief of the symptoms. The prescription drug donepezil has been widely used to treat to slow the progression and onset of the disease; however, the unpleasant side effects have paved the way to find alternative medicines. Many herbs are known to improve brain function, but evidence of medicinal plants that can treat AD is limited due to the lack of concrete rational evidences. Moreover, the traditional method of randomly screening plant extract against AD targets takes time and resources. In this study, a receptor-based in silico method has been implemented which serves to accelerate the process of identification of medicinal plants useful for treatment of AD. A database of natural compounds was compiled to identify hits against acetylcholinesterase (AChE). Receptor-based pharmacophore screening was performed, and selected hits were subjected to docking and molecular dynamics simulations. Molecular Mechanics/Generalized Born surface area (MM/GBSA) calculations were carried out to identify the best scoring hits further. In vitro assays were done for the plant extracts containing the top-scoring hits against AChE. Three plant extracts showed favorable inhibitory activity.Communicated by Ramaswamy H. Sarma.

Kinetics and Mechanism of Solvent Influence on the Lipase-Catalyzed 1,3-Diolein Synthesis

1,3-Diacylglycerol preparation has roused increasing attention in recent years as the 1,3-diacylglycerol-rich oils can suppress the deposition of visceral fat and prevent the body weight increasing. Lipozyme TL IM-mediated esterification of oleic acid with monoolein was effective for 1,3-diacylglycerol production. During the esterification process, the solvent shows obvious influence on the diolein synthesis as well as the 1,3-diolein production. This work investigated the related kinetics and mechanism of the solvent effect on the esterification and Lipozyme TL IM performance. The results indicated that both the esterification rate constant and the acyl migration rate constant positively correlated with the logP of the solvent, while the site specificity of lipase has negative correlation with solvent logP. The acylation toward the 2-position of 1-monoolein was more sensitive to the solvent logP compared to the 1-position of glycerides. Molecular dynamics simulation revealed that solvents with different logP influenced the structure of Lipozyme TL IM including RMSD, hydrogen bond, and radial distribution function to a large extent, which subsequently led to the catalytic activity and selectivity variation of the lipase.

Molecular simulation strategies for the discovery of selective inhibitors of β-catenin

Tumor dissemination and relapse in lung cancer were found to be due to the existence of cancer stem cells. In particular, the [Formula: see text]-catenin pathway is found to be one of the crucial pathways in maintaining the stem-like properties of the cells. Thus, targeting the [Formula: see text]-catenin family of proteins is a significant therapeutic route in the treatment of lung cancer. Therefore, in the present study, a pharmacophore-based drug repurposing approach was accomplished to pinpoint potent [Formula: see text]-catenin inhibitors from the DrugBank database. Primarily, ligand-based pharmacophore hypothesis (AAHHR) was generated using existing [Formula: see text]-catenin inhibitors available in the literature and utilized for library screening. Subsequently, the inhibitory activity of the screened compounds was examined by the hierarchical docking process and the Prime MM-GBSA algorithm. Moreover, quantum chemical calculations and molecular dynamics simulations were executed to analyze the inhibitory effects of the screened hit molecule. The results indicate that hit molecule, DB08047 was found to possess better binding free energy, favorable ligand strain energy, satisfactory pharmacokinetic properties and superior free energy landscape profile. Eventually, the pIC[Formula: see text] values of the lead compounds were predicted by the AutoQSAR algorithm. It is noteworthy to mention that DB08047 was found to possess pyrazole scaffolds which could downregulate [Formula: see text]-catenin pathway and thus facilitate the controlled cell growth/inhibit tumor growth.

WRKY genes family study reveals tissue-specific and stress-responsive TFs in wild potato species

Wild potatoes, as dynamic resource adapted to various environmental conditions, represent a powerful and informative reservoir of genes useful for breeding efforts. WRKY transcription factors (TFs) are encoded by one of the largest families in plants and are involved in several biological processes such as growth and development, signal transduction, and plant defence against stress. In this study, 79 and 84 genes encoding putative WRKY TFs have been identified in two wild potato relatives, Solanum commersonii and S. chacoense. Phylogenetic analysis of WRKY proteins divided ScWRKYs and SchWRKYs into three Groups and seven subGroups. Structural and phylogenetic comparative analyses suggested an interspecific variability of WRKYs. Analysis of gene expression profiles in different tissues and under various stresses allowed to select ScWRKY045 as a good candidate in wounding-response, ScWRKY055 as a bacterial infection triggered WRKY and ScWRKY023 as a multiple stress-responsive WRKY gene. Those WRKYs were further studied through interactome analysis allowing the identification of potential co-expression relationships between ScWRKYs/SchWRKYs and genes of various pathways. Overall, this study enabled the discrimination of WRKY genes that could be considered as potential candidates in both breeding programs and functional studies.