Identifying key interactions stabilizing DOF zinc finger-DNA complexes using in silico approaches

Plant-TFClass: a structural classification for plant transcription factors

Transcription factors (TFs) bind DNA at specific sequences to regulate gene expression. This universal process is achieved via their DNA-binding domain (DBD). In mammals, the vast diversity of DBD structural conformations and the way in which they contact DNA has been used to organize TFs in the TFClass hierarchical classification. However, the numerous DBD types present in plants but absent from mammalian genomes were missing from this classification. We reviewed DBD 3D structures and models available for plant TFs to classify most of the 56 recognized plant TF types within the TFClass framework. This extended classification adds eight new classes and 37 new families corresponding to DBD structures absent in mammals. Plant-TFClass provides a unique resource for TF comparison across families and organisms.

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.

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.

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.

Generation of an engineered food-grade Lactococcus lactis strain for production of an antimicrobial peptide: in vitro and in silico evaluation

BACKGROUND

Foodborne pathogens and their biofilms are considered as one of the most serious problems in human health and food industry. Moreover, safety of foods is a main global concern because of the increasing use of chemical food additives. Ensuring food safety enhances interest in discovery of new alternative compounds such as antimicrobial peptides (AMPs), which can be used as bio-preservatives in the food industry. In this study, the most important antimicrobial peptides of camel milk lactoferrin (lactoferrampin and lactoferricin) were recombinantly expressed in the form of chimeric peptide (cLFchimera) in a food-grade L. lactis strain. P170 expression system was used to express secreted cLFchimera using pAMJ1653 expression vector which harbors a safe (non-antibiotic) selectable marker.

RESULTS

Peptide purification was carried out using Ni-NTA agarose column from culture medium with concentration of 0.13 mg/mL. The results of disk diffusion test revealed that cLFchimera had considerable antimicrobial activity against a number of major foodborne bacteria. Furthermore, this chimeric peptide showed strong and weak inhibitory effect on biofilm formation against P. aeruginosa, S. aureus E. faecalis, and E. coli, respectively. Antioxidant activity and thermal stability of the chimeric peptide was determined. The results showed that cLFchimera had antioxidant activity (IC50: 310 μ/mL) and its activity was not affected after 40 min of boiling. Finally, we evaluated the interaction of the peptide with LPS and DNA in bacteria using molecular dynamic simulation as two main intra and extra cellular targets for AMPs, respectively. Our in silico analysis showed that cLFchimera had strong affinity to both of these targets by positive charged residues after 50 ns molecular dynamic simulation.

CONCLUSIONS

Overall, the engineered food-grade L. lactis generated in the present study successfully expressed a secreted chimeric peptide with antimicrobial properties and could be considered as a promising bio-preservative in the food industry.

Genome-wide identification and expression profiling of the C2H2-type zinc finger protein genes in the silkworm Bombyx mori

Cys2-His2 zinc finger (C2H2-ZF) proteins comprise the largest class of putative eukaryotic transcription factors. The zinc finger motif array is highly divergent, indicating that most proteins will have distinctive binding sites and perform different functions. However, the binding sites and functions of the majority of C2H2-ZF proteins remain unknown. In this study, we identified 327 C2H2-ZF protein genes in the silkworm, 290 in the monarch butterfly, 243 in the fruit fly, 107 in elegans, 673 in mouse, and 1,082 in human. The C2H2-ZF protein genes of the silkworm were classified into three main grouping clades according to a phylogenetic classification, and 312 of these genes could be mapped onto 27 chromosomes. Most silkworm C2H2-ZF protein genes exhibited specific expression in larval tissues. Furthermore, several C2H2-ZF protein genes had sex-specific expression during metamorphosis. In addition, we found that some C2H2-ZF protein genes are involved in metamorphosis and female reproduction by using expression clustering and gene annotation analysis. Among them, five genes were selected, BGIBMGA002091 (CTCF), BGIBMGA006492 (fru), BGIBMGA006230 (wor), BGIBMGA004640 (lola), and BIGBMGA004569, for quantitative real-time PCR analysis from larvae to adult ovaries. The results showed that the five genes had different expression patterns in ovaries, among which BGIBMGA002091 (CTCF) gene expression level was the highest, and its expression level increased rapidly in late pupae and adult stages. These findings provide a basis for further investigation of the functions of C2H2-ZF protein genes in the silkworm, and the results offer clues for further research into the development of metamorphosis and female reproduction in the silkworm.

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.

Expression, purification and DNA-binding properties of zinc finger domains of DOF proteins from Arabidopsis thaliana

Introduction: DOF proteins are a family of plant-specific transcription factors with a conserved zinc finger (ZF) DNA-binding domain. Although several studies have demonstrated their specific DNA binding, quantitative affinity data is not available for the binding of DOF domains to their binding sites.

Methods: ZF domains of DOF2.1, DOF3.4, and DOF5.8 from Arabidopsis thaliana were expressed and purified. Their DNA binding affinities were assessed using gel retardation assays and microscale thermophoresis with two different oligonucleotide probes containing one and two copies of recognition sequence AAAG.

Results: DOF zinc finger domains (DOF-ZFs) were shown to form independently folded structures. Assessments using microscale thermophoresis demonstrated that DOF-ZFs interact more tightly (~ 100 fold) with double-motif probe than the single-motif probe. The overall K_d values for the DOF3.4-ZF and DOF5.8-ZF to the double-motif probe were ~2.3±1 and 2.5±1 µM, respectively.

Conclusion: Studied DOF-ZF domains formed stable complexes with the double-motif probe. Although DOF3.4-ZF and DOF5.8-ZF do not dimerize with an appreciable affinity in the absence of DNA (judging from size-exclusion and multiangle laser light scattering data), it is possible that these ZFs form protein-protein contacts when bound to this oligonucleotide, consistent with previous reports that DOF proteins can homo- and hetero-dimerize.

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.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4506-3) contains supplementary material, which is available to authorized users.