Search

PLCG1 gene is responsible for many T-cell lymphoma subtypes, including peripheral T-cell lymphoma (PTCL), angioimmunoblastic T-cell lymphoma (AITL), cutaneous T-cell lymphoma (CTCL), adult T-cell leukemia/lymphoma along with other diseases. Missense mutations of this gene have already been found in patients of CTCL and AITL. The non-synonymous single nucleotide polymorphisms (nsSNPs) can alter the protein structure as well as its functions. In this study, probable deleterious and disease-related nsSNPs in PLCG1 were identified using SIFT, PROVEAN, PolyPhen-2, PhD-SNP, Pmut, and SNPS&GO tools. Further, their effect on protein stability was checked along with conservation and solvent accessibility analysis by I-mutant 2.0, MUpro, Consurf, and Netsurf 2.0 server. Some SNPs were finalized for structural analysis with PyMol and BIOVIA discovery studio visualizer. Out of the 16 nsSNPs which were found to be deleterious, ten nsSNPs had an effect on protein stability, and six mutations (L411P, R355C, G493D, R1158H, A401V and L455F) were predicted to be highly conserved. Among the six highly conserved mutations, four nsSNPs (R355C, A401V, L411P and L455F) were part of the catalytic domain. L411P, L455F and G493D made significant structural change in the protein structure. Two mutations-Y210C and R1158H had post-translational modification. In the 5' and 3' untranslated region, three SNPs, rs139043247, rs543804707, and rs62621919 showed possible miRNA target sites and DNA binding sites. This in silico analysis has provided a structured dataset of PLCG1 gene for further in vivo researches. With the limitation of computational study, it can still prove to be an asset for the identification and treatment of multiple diseases associated with the target gene., Competing Interests: The authors have declared that no competing interests exist.

The twin-arginine translocation (Tat) pathway transports folded proteins across energetic membranes. Numerous Tat substrates contain co-factors that are inserted before transport with the assistance of redox enzyme maturation proteins (REMPs), which bind to the signal peptide of precursor proteins. How signal peptides are transferred from a REMP to a binding site on the Tat receptor complex remains unknown. Since the signal peptide mediates both interactions, possibilities include: i) a coordinated hand-off mechanism; or ii) a diffusional search after REMP dissociation. We investigated the binding interaction between substrates containing the TorA signal peptide (spTorA) and its cognate REMP, TorD, and the effect of TorD on the in vitro transport of such substrates. We found that Escherichia coli TorD is predominantly a monomer at low micromolar concentrations (dimerization KD > 50 μM), and this monomer binds reversibly to spTorA (KD ≈ 1 μM). While TorD binds to membranes (KD ≈ 100 nM), it has no apparent affinity for Tat translocons and it inhibits binding of a precursor substrate to the membrane. TorD has a minimal effect on substrate transport by the Tat system, being mildly inhibitory at high concentrations. These data are consistent with a model in which the REMP-bound signal peptide is shielded from recognition by the Tat translocon, and spontaneous dissociation of the REMP allows the substrate to engage the Tat machinery. Thus, the REMP does not assist with targeting to the Tat translocon, but rather temporarily shields the signal peptide., Competing Interests: The authors have declared that no competing interests exist.

Understanding how human ACE2 genetic variants differ in their recognition by SARS-CoV-2 can facilitate the leveraging of ACE2 as an axis for treating and preventing COVID-19. In this work, we experimentally interrogate thousands of ACE2 mutants to identify over one hundred human single-nucleotide variants (SNVs) that are likely to have altered recognition by the virus, and make the complementary discovery that ACE2 residues distant from the spike interface influence the ACE2-spike interaction. These findings illuminate new links between ACE2 sequence and spike recognition, and could find substantial utility in further fundamental research that augments epidemiological analyses and clinical trial design in the contexts of both existing strains of SARS-CoV-2 and novel variants that may arise in the future., Competing Interests: The authors have declared that no competing interests exist.

The spike protein receptor-binding domain (RBD) of SARS-CoV-2 is the molecular target for many vaccines and antibody-based prophylactics aimed at bringing COVID-19 under control. Such a narrow molecular focus raises the specter of viral immune evasion as a potential failure mode for these biomedical interventions. With the emergence of new strains of SARS-CoV-2 with altered transmissibility and immune evasion potential, a critical question is this: how easily can the virus escape neutralizing antibodies (nAbs) targeting the spike RBD? To answer this question, we combined an analysis of the RBD structure-function with an evolutionary modeling framework. Our structure-function analysis revealed that epitopes for RBD-targeting nAbs overlap one another substantially and can be evaded by escape mutants with ACE2 affinities comparable to the wild type, that are observed in sequence surveillance data and infect cells in vitro. This suggests that the fitness cost of nAb-evading mutations is low. We then used evolutionary modeling to predict the frequency of immune escape before and after the widespread presence of nAbs due to vaccines, passive immunization or natural immunity. Our modeling suggests that SARS-CoV-2 mutants with one or two mildly deleterious mutations are expected to exist in high numbers due to neutral genetic variation, and consequently resistance to vaccines or other prophylactics that rely on one or two antibodies for protection can develop quickly -and repeatedly- under positive selection. Predicted resistance timelines are comparable to those of the decay kinetics of nAbs raised against vaccinal or natural antigens, raising a second potential mechanism for loss of immunity in the population. Strategies for viral elimination should therefore be diversified across molecular targets and therapeutic modalities., Competing Interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: A.C., M.S., and U.T. are employees and shareholders of Fractal Therapeutics. D.V.E., A.N., B.Z., and D.J.- M. are shareholders of Fractal Therapeutics. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Cyp33 is an essential human cyclophilin prolyl isomerase that plays myriad roles in splicing and chromatin remodeling. In addition to a canonical cyclophilin (Cyp) domain, Cyp33 contains an RNA-recognition motif (RRM) domain, and RNA-binding triggers proline isomerase activity. One prominent role for Cyp33 is through a direct interaction with the mixed lineage leukemia protein 1 (MLL1, also known as KMT2A) complex, which is a histone methyltransferase that serves as a global regulator of human transcription. MLL activity is regulated by Cyp33, which isomerizes a key proline in the linker between the PHD3 and Bromo domains of MLL1, acting as a switch between gene activation and repression. The direct interaction between MLL1 and Cyp33 is critical, as deletion of the MLL1-PHD3 domain responsible for this interaction results in oncogenesis. The Cyp33 RRM is central to these activities, as it binds both the PHD3 domain and RNA. To better understand how RNA binding drives the action of Cyp33, we performed RNA-SELEX against full-length Cyp33 accompanied by deep sequencing. We have identified an enriched Cyp33 binding motif (AAUAAUAA) broadly represented in the cellular RNA pool as well as tightly binding RNA aptamers with affinities comparable and competitive with the Cyp33 MLL1-PHD3 interaction. RNA binding extends beyond the canonical RRM domain, but not to the Cyp domain, suggesting an indirect mechanism of interaction. NMR chemical shift mapping confirms an overlapping, but not identical, interface on Cyp33 for RNA and PHD3 binding. This finding suggests RNA can disrupt the gene repressive Cyp33-MLL1 complex providing another layer of regulation for chromatin remodeling by MLL1., Competing Interests: The authors have declared that no competing interests exist.

Common variants in the hepatic lipase (LIPC) gene have been shown to be associated with plasma lipid levels; however, the distribution and functional features of rare and regulatory LIPC variants contributing to the extreme lipid phenotypes are not well known. This study was aimed to catalogue LIPC variants by resequencing the entire LIPC gene in 95 non-Hispanic Whites (NHWs) and 95 African blacks (ABs) with extreme HDL-C levels followed by in silico functional analyses. A total of 412 variants, including 43 novel variants were identified; 56 were unique to NHWs and 234 were unique to ABs. Seventy-eight variants in NHWs and 89 variants in ABs were present either in high HDL-C group or low HDL-C group. Two non-synonymous variants (p.S289F, p.T405M), found in NHWs with high HDL-C group were predicted to have damaging effect on LIPC protein by SIFT, MT2 and PP2. We also found several non-coding variants that possibly reside in the circRNA and lncRNA binding sites and may have regulatory potential, as identified in rSNPbase and RegulomeDB databases. Our results shed light on the regulatory nature of rare and non-coding LIPC variants as well as suggest their important contributions in affecting the extreme HDL-C phenotypes., Competing Interests: The authors have declared that no competing interests exist.

Sex-specific transcription characterizes hundreds of genes in mouse liver, many implicated in sex-differential drug and lipid metabolism and disease susceptibility. While the regulation of liver sex differences by growth hormone-activated STAT5 is well established, little is known about autosomal genetic factors regulating the sex-specific liver transcriptome. Here we show, using genotyping and expression data from a large population of Diversity Outbred mice, that genetic factors work in tandem with growth hormone to control the individual variability of hundreds of sex-biased genes, including many long non-coding RNA genes. Significant associations between single nucleotide polymorphisms and sex-specific gene expression were identified as expression quantitative trait loci (eQTLs), many of which showed strong sex-dependent associations. Remarkably, autosomal genetic modifiers of sex-specific genes were found to account for more than 200 instances of gain or loss of sex-specificity across eight Diversity Outbred mouse founder strains. Sex-biased STAT5 binding sites and open chromatin regions with strain-specific variants were significantly enriched at eQTL regions regulating correspondingly sex-specific genes, supporting the proposed functional regulatory nature of the eQTL regions identified. Binding of the male-biased, growth hormone-regulated repressor BCL6 was most highly enriched at trans-eQTL regions controlling female-specific genes. Co-regulated gene clusters defined by overlapping eQTLs included sets of highly correlated genes from different chromosomes, further supporting trans-eQTL action. These findings elucidate how an unexpectedly large number of autosomal factors work in tandem with growth hormone signaling pathways to regulate the individual variability associated with sex differences in liver metabolism and disease., Competing Interests: The authors have declared that no competing interests exist.

The interactions between membrane receptors and their endogenous ligands are key interactions in organisms. Recently, we have shown that a high number of genes encoding human receptors appeared at the same moment as their ligand in the animal tree of life. However, a set of receptors appeared before their present ligand. Different scenarios have been proposed to explain how a receptor can be conserved if its ligand is not yet appeared. However, these scenarios have been proposed individually and have never been studied in a global way. In this study, we investigated 30 mammalian pairs of ligand/receptor for which the first ligand appeared after its receptor in the tree of life, by using common indexes of selection, and proposed different scenarios explaining the earlier appearance of a receptor relative to its ligand. Based on 3D structural studies, our indexes allowed us to classify the evolution of these partners into different scenarios: 1) a scenario where the binding interface of the receptor is already present and under purifying selection before the appearance of the ligand; 2) a scenario where the binding interface seems to have appeared progressively, and 3) a scenario where the binding site seems to have been reshuffled since its appearance. As some scenarios were confirmed by the literature, we concluded that simple indexes can give a good highlight of the evolutive history of two partners that did not appear at the same time. Based on these scenarios, we also hypothesize that the replacement of a ligand by another is a frequent phenomenon during evolution., Competing Interests: The authors have declared that no competing interests exist.

Following nerve stimulation, there are two distinct phases of Ca2+-dependent neurotransmitter release: a fast, synchronous release phase, and a prolonged, asynchronous release phase. Each of these phases is tightly regulated and mediated by distinct mechanisms. Synaptotagmin 1 is the major Ca2+ sensor that triggers fast, synchronous neurotransmitter release upon Ca2+ binding by its C2A and C2B domains. It has also been implicated in the inhibition of asynchronous neurotransmitter release, as blocking Ca2+ binding by the C2A domain of synaptotagmin 1 results in increased asynchronous release. However, the mutation used to block Ca2+ binding in the previous experiments (aspartate to asparagine mutations, sytD-N) had the unintended side effect of mimicking Ca2+ binding, raising the possibility that the increase in asynchronous release was directly caused by ostensibly constitutive Ca2+ binding. Thus, rather than modulating an asynchronous sensor, sytD-N may be mimicking one. To directly test the C2A inhibition hypothesis, we utilized an alternate C2A mutation that we designed to block Ca2+ binding without mimicking it (an aspartate to glutamate mutation, sytD-E). Analysis of both the original sytD-N mutation and our alternate sytD-E mutation at the Drosophila neuromuscular junction showed differential effects on asynchronous release, as well as on synchronous release and the frequency of spontaneous release. Importantly, we found that asynchronous release is not increased in the sytD-E mutant. Thus, our work provides new mechanistic insight into synaptotagmin 1 function during Ca2+-evoked synaptic transmission and demonstrates that Ca2+ binding by the C2A domain of synaptotagmin 1 does not inhibit asynchronous neurotransmitter release in vivo., Competing Interests: The authors have declared that no competing interests exist.

Recent evidence shows that the disruption of constitutive insulated neighbourhoods might lead to oncogene dysregulation. We present here a systematic pan-cancer characterisation of the associations between constitutive boundaries and genome alterations in cancer. Specifically, we investigate the enrichment of somatic mutation, abnormal methylation, and copy number alteration events in the proximity of CTCF bindings overlapping with topological boundaries (junctions) in 26 cancer types. Focusing on CTCF motifs that are both in-boundary (overlapping with junctions) and active (overlapping with peaks of CTCF expression), we find a significant enrichment of somatic mutations in several cancer types. Furthermore, mutated junctions are significantly conserved across cancer types, and we also observe a positive selection of transversions rather than transitions in many cancer types. We also analyzed the mutational signature found on the different classes of CTCF motifs, finding some signatures (such as SBS26) to have a higher weight within in-boundary than off-bounday motifs. Regarding methylation, we find a significant number of over-methylated active in-boundary CTCF motifs in several cancer types; similarly to somatic-mutated junctions, they also have a significant conservation across cancer types. Finally, in several cancer types we observe that copy number alterations tend to overlap with active junctions more often than in matched normal samples. While several articles have recently reported a mutational enrichment at CTCF binding sites for specific cancer types, our analysis is pan-cancer and investigates abnormal methylation and copy number alterations in addition to somatic mutations. Our method is fully replicable and suggests several follow-up tumour-specific analyses., Competing Interests: The authors have declared that no competing interests exist.

The HIV genome is rich in A but not G or U and deficient in C. This nucleotide bias controls HIV phenotype by determining the highly unusual composition of all major HIV proteins. The bias is also responsible for the high frequency of narrow DNA minor groove sites in the double-stranded HIV genome as compared to cellular protein coding sequences and the bulk of the human genome. Since drugs that bind in the DNA minor groove disrupt nucleosomes on sequences that contain closely spaced oligo-A tracts which are prevalent in HIV DNA because of its bias, it was of interest to determine if these drugs exert this selective inhibitory effect on HIV chromatin. To test this possibility, nucleosomes were reconstituted onto five double-stranded DNA fragments from the HIV-1 pol gene in the presence and in the absence of several minor groove binding drugs (MGBDs). The results demonstrated that the MGBDs inhibited the assembly of nucleosomes onto all of the HIV-1 segments in a manner that was proportional to the A-bias, but had no detectable effect on the formation of nucleosomes on control cloned fragments or genomic DNA from chicken and human. Nucleosomes preassembled onto HIV DNA were also preferentially destabilized by the drugs as evidenced by enhanced nuclease accessibility in physiological ionic strength and by the preferential loss of the histone octamer in hyper-physiological salt solutions. The drugs also selectively disrupted HIV-containing nucleosomes in yeast as revealed by enhanced nuclease accessibility of the in vivo assembled HIV chromatin and reductions in superhelical densities of plasmid chromatin containing HIV sequences. A comparison of these results to the density of A-tracts in the HIV genome indicates that a large fraction of the nucleosomes that make up HIV chromatin should be preferred in vitro targets for the MGBDs. These results show that the MGBDs preferentially disrupt HIV-1 chromatin in vitro and in vivo and raise the possibility that non-toxic derivatives of certain MGBDs might serve as a novel class of anti-HIV agents., Competing Interests: The authors have declared that no competing interests exist.

It has been reported that neuroinflammation occurs in the central nervous system (CNS) in patients with neuropathic pain, Alzheimer's disease and autism spectrum disorder. The 18-kDa translocator protein TSPO is used as an imaging target in positron emission tomography to detect neuroinflammation, and its expression is correlated with microglial activation. However, the mechanism underlying the transcriptional regulation of Tspo induced by inflammation is not clear. Here, we revealed that lipopolysaccharide (LPS) -induced Tspo expression was activated by the AP-1 complex in a mouse microglial cell line, BV-2. Knockdown of c-Fos and c-Jun, the components of AP-1, reduced LPS-induced Tspo expression. Furthermore, the enrichment of Sp1 in the proximal promoter region of Tspo was increased in the presence of LPS. In addition, the binding of histone deacetylase 1 (HDAC1) to the enhancer region, which contains the AP-1 site, was decreased by LPS treatment, but there were no significant differences in HDAC1 binding to the proximal promoter region with or without LPS. These results indicated that HDAC1 is involved not in the proximal promoter region but in the enhancer region. Our study revealed that inflammatory signals induce the recruitment of AP-1 to the enhancer region and Sp1 to the proximal promoter region of the Tspo gene and that Sp1 may regulate the basal expression of Tspo., Competing Interests: The authors have declared that no competing interests exist.

Estrogen receptor α (ERα) is the major driving transcription factor in the mammary gland development as well as breast cancer initiation and progression. However, the genomic landscape of ERα binding sites in the normal mouse mammary gland has not been completely elucidated. Here, we mapped genome-wide ERα binding events by chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) in the mouse mammary gland in response to estradiol. We identified 6237 high confidence ERα binding sites in two biological replicates and showed that many of these were located at distal enhancer regions. Furthermore, we discovered 3686 unique genes in the mouse genome that recruit ER in response to estradiol. Interrogation of ER-DNA binding sites in ER-positive luminal epithelial cells showed that the ERE, PAX2, SF1, and AP1 motifs were highly enriched at distal enhancer regions. In addition, comprehensive transcriptome analysis by RNA-seq revealed that 493 genes are differentially regulated by acute treatment with estradiol in the mouse mammary gland in vivo. Through integration of RNA-seq and ERα ChIP-seq data, we uncovered a novel ERα targetome in mouse mammary epithelial cells. Taken together, our study has identified the genomic landscape of ERα binding events in mouse mammary epithelial cells. Furthermore, our study also highlights the cis-regulatory elements and cofactors that are involved in estrogen signaling and may contribute to ductal elongation in the normal mouse mammary gland., Competing Interests: The authors have declared that no competing interests exist.

Background: Great progress has been made in the understanding of inflammatory processes in psoriasis. However, clarifying the role of genetic variability in processes regulating inflammation, including post-transcriptional regulation by microRNA (miRNA), remains a challenge., Objectives: We therefore investigated single nucleotide polymorphisms (SNPs) with a predicted change in the miRNA/mRNA interaction of genes involved in the psoriasis inflammatory processes., Methods: Studied SNPs rs2910164 C/G-miR-146a, rs4597342 T/C-ITGAM, rs1368439 G/T-IL12B, rs1468488 C/T-IL17RA were selected using a bioinformatics analysis of psoriasis inflammation-associated genes. These SNPs were then genotyped using a large cohort of women with psoriasis (n = 241) and healthy controls (n = 516)., Results: No significant association with psoriasis was observed for rs2910164, rs1368439, and rs1468488 genotypes. However, the major allele T of rs4597342 -ITGAM was associated with approximately 28% higher risk for psoriasis in comparison to the patients with the C allele (OR = 1.28, 95% CI 1.01-1.61, p = 0.037). In case of genotypes, the effect of the T allele indicates the dominant model of disease penetrance as the CT and TT genotypes increase the chance of psoriasis up to 42% in comparison to CC homozygotes of rs4597342 (OR = 1.42, 95% CI = 1.05-1.94, p = 0.025)., Conclusion: SNP rs4597342 in 3'UTR of ITGAM influencing miR-21 binding may be considered a risk factor for psoriasis development. Upregulated miR-21 in psoriasis is likely to inhibit CD11b production in the case of the rs4597342 T allele which may lead to Mac-1 dysfunction, resulting in an aberrant function of innate immune cells and leading to the production of cytokines involved in psoriasis pathogenesis., Competing Interests: The authors have declared that no competing interests exist.

C/EBPα plays a key role in specifying myeloid lineage development. HoxA9 is expressed in myeloid progenitors, with its level diminishing during myeloid maturation, and HOXA9 is over-expressed in a majority of acute myeloid leukemia cases, including those expressing NUP98-HOXD13. The objective of this study was to determine whether HoxA9 directly represses Cebpa gene expression. We find 4-fold increased HoxA9 and 5-fold reduced Cebpa in marrow common myeloid and LSK progenitors from Vav-NUP98-HOXD13 transgenic mice. Conversely, HoxA9 decreases 5-fold while Cebpa increases during granulocytic differentiation of 32Dcl3 myeloid cells. Activation of exogenous HoxA9-ER in 32Dcl3 cells reduces Cebpa mRNA even in the presence of cycloheximide, suggesting direct repression. Cebpa transcription in murine myeloid cells is regulated by a hematopoietic-specific +37 kb enhancer and by a more widely active +8 kb enhancer. ChIP-Seq analysis of primary myeloid progenitor cells expressing exogenous HoxA9 or HoxA9-ER demonstrates that HoxA9 localizes to both the +8 kb and +37 kb Cebpa enhancers. Gel shift analysis demonstrates HoxA9 binding to three consensus sites in the +8 kb enhancer, but no affinity for the single near-consensus site present in the +37 kb enhancer. Activity of a Cebpa +8 kb enhancer/promoter-luciferase reporter in 32Dcl3 or MOLM14 myeloid cells is increased ~2-fold by mutation of its three HOXA9-binding sites, suggesting that endogenous HoxA9 represses +8 kb Cebpa enhancer activity. In contrast, mutation of five C/EBPα-binding sites in the +8 kb enhancer reduces activity 3-fold. Finally, expression of a +37 kb enhancer/promoter-hCD4 transgene reporter is reduced ~2-fold in marrow common myeloid progenitors when the Vav-NUP98-HOXD13 transgene is introduced. Overall, these data support the conclusion that HoxA9 represses Cebpa expression, at least in part via inhibition of its +8 kb enhancer, potentially allowing normal myeloid progenitors to maintain immaturity and contributing to the pathogenesis of acute myeloid leukemia associated with increased HOXA9., Competing Interests: The authors have declared that no competing interests exist.

Fragile X syndrome, the most common inherited form of intellectual disability, is caused by the CGG trinucleotide expansion in the 5'-untranslated region of the Fmr1 gene on the X chromosome, which silences the expression of the fragile X mental retardation protein (FMRP). FMRP has been shown to bind to a G-rich region within the PSD-95 mRNA, which encodes for the postsynaptic density protein 95, and together with microRNA-125a to mediate the reversible inhibition of the PSD-95 mRNA translation in neurons. The miR-125a binding site within the PSD-95 mRNA 3'-untranslated region (UTR) is embedded in a G-rich region bound by FMRP, which we have previously demonstrated folds into two parallel G-quadruplex structures. The FMRP regulation of PSD-95 mRNA translation is complex, being mediated by its phosphorylation. While the requirement for FMRP in the regulation of PSD-95 mRNA translation is clearly established, the exact mechanism by which this is achieved is not known. In this study, we have shown that both unphosphorylated FMRP and its phosphomimic FMRP S500D bind to the PSD-95 mRNA G-quadruplexes with high affinity, whereas only FMRP S500D binds to miR-125a. These results point towards a mechanism by which, depending on its phosphorylation status, FMRP acts as a switch that potentially controls the stability of the complex formed by the miR-125a-guided RNA induced silencing complex (RISC) and PSD-95 mRNA., Competing Interests: The authors have declared that no competing interests exist.

In the search for new pharmaceutical leads, especially with DNA-binding molecules or genome editing methods, the issue of side and off-target effects have always been thorny in nature. A particular case is the investigation into the off-target effects of N-methylpyrrole-N-methylimidazole polyamides, a naturally inspired class of DNA binders with strong affinity to the minor-groove and sequence specificity, but at < 20 bases, their relatively short motifs also insinuate the possibility of non-unique genomic binding. Binding at non-intended loci potentially lead to the rise of off-target effects, issues that very few approaches are able to address to-date. We here report an analytical method to infer off-target binding, via expression profiling, based on probing the relative impact to various biochemical pathways; we also proposed an accompanying side effect prediction engine for the systematic screening of candidate polyamides. This method marks the first attempt in PI polyamide research to identify elements in biochemical pathways that are sensitive to the treatment of a candidate polyamide as an approach to infer possible off-target effects. Expression changes were then considered to assess possible outward phenotypic changes, manifested as side effects, should the same PI polyamide candidate be administered clinically. We validated some of these effects with a series of animal experiments, and found agreeable corroboration in certain side effects, such as changes in aspartate transaminase levels in ICR and nude mice post-administration., Competing Interests: The authors have declared that no competing interests exist.

Long-range associations between enhancers and their target gene promoters have been shown to play critical roles in executing genome function. Recent variations of chromosome capture technology have revealed a comprehensive view of intra- and interchromosomal contacts between specific genomic sites. The locus control region of the β-globin genes (β-LCR) is a super-enhancer that is capable of activating all of the β-like globin genes within the locus in cis through physical interaction by forming DNA loops. CTCF helps to mediate loop formation between LCR-HS5 and 3'HS1 in the human β-globin locus, in this way thought to contribute to the formation of a "chromatin hub". The β-globin locus is also in close physical proximity to other erythrocyte-specific genes located long distances away on the same chromosome. In this case, erythrocyte-specific genes gather together at a shared "transcription factory" for co-transcription. Theoretically, enhancers could also activate target gene promoters at the identical loci, yet on different chromosomes in trans, a phenomenon originally described as transvection in Drosophilla. Although close physical proximity has been reported for the β-LCR and the β-like globin genes when integrated at the mouse homologous loci in trans, their structural and functional interactions were found to be rare, possibly because of a lack of suitable regulatory elements that might facilitate such trans interactions. Therefore, we re-evaluated presumptive transvection-like enhancer-promoter communication by introducing CTCF binding sites and erythrocyte-specific transcription units into both LCR-enhancer and β-promoter alleles, each inserted into the mouse ROSA26 locus on separate chromosomes. Following cross-mating of mice to place the two mutant loci at the identical chromosomal position and into active chromation in trans, their transcriptional output was evaluated. The results demonstrate that there was no significant functional association between the LCR and the β-globin gene in trans even in this idealized experimental context., Competing Interests: The authors have declared that no competing interests exist.

For unknown reasons, humans appear to be particular susceptible to developing tau pathology leading to neurodegeneration. Transgenic mice are still undoubtedly the most popular and extensively used animal models for studying Alzheimer's disease and other tauopathies. While these murine models generally overexpress human tau in the mouse brain or specific brain regions, there are differences between endogenous mouse tau and human tau protein. Among them, a main difference between human and mouse tau is the presence of a short motif spanning residues 18 to 28 in the human tau protein that is missing in murine tau, and which could be at least partially responsible for that different susceptibility across species. Here we report novel data using affinity chromatography analysis indicating that the sequence containing human tau residues 18 to 28 acts a binding motif for End Binding proteins and that this interaction could facilitate tau secretion to the extracellular space., Competing Interests: The authors have declared that no competing interests exist.

Through the application of a region-focused saturation mutagenesis and randomization approach, protein engineering of the Cal-A enzyme was undertaken with the goal of conferring new triglyceride selectivity. Little is known about the mode of triglyceride binding to Cal-A. Engineering Cal-A thus requires a systemic approach. Targeted and randomized Cal-A libraries were created, recombined using the Golden Gate approach and screened to detect variants able to discriminate between long-chain (olive oil) and short-chain (tributyrin) triglyceride substrates using a high-throughput in vivo method to visualize hydrolytic activity. Discriminative variants were analyzed using an in-house script to identify predominant substitutions. This approach allowed identification of variants that exhibit strong discrimination for the hydrolysis of short-chain triglycerides and others that discriminate towards hydrolysis of long-chain triglycerides. A clear pattern emerged from the discriminative variants, identifying the 217-245 helix-loop-helix motif as being a hot-spot for triglyceride recognition. This was the consequence of introducing the entire mutational load in selected regions, without putting a strain on distal parts of the protein. Our results improve our understanding of the Cal-A lipase mode of action and selectivity. This holistic perspective to protein engineering, where parts of the gene are individually mutated and the impact evaluated in the context of the whole protein, can be applied to any protein scaffold., Competing Interests: The authors have declared that no competing interests exist.

Mammalian gene expression is a complex process regulated in part by CpG methylation. The ability to target methylation for de novo gene regulation could have therapeutic and research applications. We have previously developed a dCas9-MC/MN protein for targeting CpG methylation. dCas9-MC/MN is composed of an artificially split M.SssI methyltransferase (MC/MN), with the MC fragment fused to a nuclease-null CRISPR/Cas9 (dCas9). Guide RNAs directed dCas9-MC/MN to methylate target sites in E. coli and human cells but also caused some low-level off-target methylation. Here, in E. coli, we show that shortening the dCas9-MC linker increases methylation of CpG sites located at select distances from the dCas9 binding site. Although a shortened linker decreased methylation of other CpGs proximal to the target site, it did not reduce off-target methylation of more distant CpG sites. Instead, targeted mutagenesis of the methyltransferase's DNA binding domain, designed to reduce DNA affinity, significantly and preferentially reduced methylation of such sites., Competing Interests: TX, CDN, and MO are inventors on a patent application (US Patent Application 15/539256 Systems and methods for genome modification and regulation), which is related to the experiments found in this manuscript. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Adenosine deaminase (ADA) catalyzes the deamination of adenosine, which is important in purine metabolism. ADA is ubiquitous to almost all human tissues, and ADA abnormalities have been reported in various diseases, including rheumatoid arthritis. ADA can be divided into two conformations based on the inhibitor that it binds to: open and closed forms. Here, we chose three ligands, namely, FR117016 (FR0), FR221647 (FR2) (open form), and HDPR (PRH, closed form), to investigate the inhibition mechanism of ADA and its effect on ADA through molecular dynamics simulations. In open forms, Egap and electrostatic potential (ESP) indicated that electron transfer might occur more easily in FR0 than in FR2. Binding free energy and hydrogen bond occupation revealed that the ADA-FR0 complex had a more stable structure than ADA-FR2. The probability of residues Pro159 to Lys171 of ADA-FR0 and ADA-FR2 to form a helix moderately increased compared with that in nonligated ADA. In comparison with FR0 and FR2 PRH could maintain ADA in a closed form to inhibit the function of ADA. The α7 helix (residues Thr57 to Ala73) of ADA in the closed form was mostly unfastened because of the effect of PRH. The number of H bonds and the relative superiority of the binding free energy indicated that the binding strength of PRH to ADA was significantly lower than that of an open inhibitor, thereby supporting the comparison of the inhibitory activities of the three ligands. Alanine scanning results showed that His17, Gly184, Asp295, and Asp296 exerted the greatest effects on protein energy, suggesting that they played crucial roles in binding to inhibitors. This study served as a theoretical basis for the development of new ADA inhibitors., Competing Interests: The authors have declared that no competing interests exist.

Regulation of mRNA stability by RNA-protein interactions contributes significantly to quantitative aspects of gene expression. We have identified potential mRNA targets of the AU-rich element binding protein AUF1. Myc-tagged AUF1 p42 was induced in mouse NIH/3T3 cells and RNA-protein complexes isolated using anti-myc tag antibody beads. Bound mRNAs were analyzed with Affymetrix microarrays. We have identified 508 potential target mRNAs that were at least 3-fold enriched compared to control cells without myc-AUF1. 22.3% of the enriched mRNAs had an AU-rich cluster in the ARED Organism database, against 16.3% of non-enriched control mRNAs. The enrichment towards AU-rich elements was also visible by AREScore with an average value of 5.2 in the enriched mRNAs versus 4.2 in the control group. Yet, numerous mRNAs were enriched without a high ARE score. The enrichment of tetrameric and pentameric sequences suggests a broad AUF1 p42-binding spectrum at short U-rich sequences flanked by A or G. Still, some enriched mRNAs were highly unstable, as those of TNFSF11 (known as RANKL), KLF10, HES1, CCNT2, SMAD6, and BCL6. We have mapped some of the instability determinants. HES1 mRNA appeared to have a coding region determinant. Detailed analysis of the RANKL and BCL6 3'UTR revealed for both that full instability required two elements, which are conserved in evolution. In RANKL mRNA both elements are AU-rich and separated by 30 bases, while in BCL6 mRNA one is AU-rich and 60 bases from a non AU-rich element that potentially forms a stem-loop structure., Competing Interests: The authors have declared that no competing interests exist.

Ehrlichia chaffeensis is an obligately intracellular bacterium that establishes infection in mononuclear phagocytes through largely undefined reprogramming strategies including modulation of host gene transcription. In this study, we demonstrate that the E. chaffeensis effector TRP47 enters the host cell nucleus and binds regulatory regions of host genes relevant to infection. TRP47 was observed in the nucleus of E. chaffeensis-infected host cells, and nuclear localization was dependent on a variant MYND-binding domain. An electrophoretic mobility shift assay (EMSA) demonstrated that TRP47 directly binds host DNA via its tandem repeat domain. Utilizing chromatin immunoprecipitation followed by high-throughput DNA sequencing (ChIP-seq) with E. chaffeensis-infected cells, TRP47 was found to bind at multiple sites in the human genome (n = 2,051 at p < 10-30). Ontology analysis identified genes involved in functions such as immune response, cytoskeletal organization, and signal transduction. TRP47-bound genes included RNA-coding genes, many of these linked to cell proliferation or apoptosis. Comparison of TRP47 binding sites with those of previously-identified E. chaffeensis nucleomodulins identified multiple genes and gene functional categories in common including intracellular transport, cell signaling, and transcriptional regulation. Further, motif analysis followed by EMSA with synthetic oligonucleotides containing discovered motifs revealed a conserved TRP47 DNA-binding motif. This study reveals that TRP47 is a nucleomodulin that enters the nucleus via a MYND-binding domain and appears to play a role in host cell reprogramming by regulation of transcription., Competing Interests: The authors have declared that no competing interests exist.

Sucrose synthase (SuSy) is one key enzyme directly hydrolyzing sucrose to supply substrates for plant metabolism, and is considered to be a biomarker for plant sink strength. Improvement in plant sink strength could lead to enhanced plant growth and yield. Cultivated tomatoes are known to have a narrow genetic diversity, which hampers further breeding for novel and improved traits in new cultivars. In this study, we observed limited genetic variation in SuSy1, SuSy3 and SuSy4 in 53 accessions of cultivated tomato and landraces, but identified a wealth of genetic diversity in 32 accessions of related wild species. The variation in the deduced amino acid sequences was grouped into 23, 22, and 17 distinct haplotypes for SuSy1/3/4, respectively. Strikingly, all known substrate binding sites were highly conserved, as well as most of the phosphorylation sites except in SuSy1. Two SuSy1 and three SuSy3 protein variants were heterologously expressed to study the effect of the amino acid changes on enzyme kinetic properties, i.e. maximal sucrose hydrolyzing capacity (Vmax), affinity for sucrose (Km), and catalytic efficiency (Vmax/Km) at 25°C and 16°C. SuSy1-haplotype#3 containing phosphorylation site Ser-16 did not have an improvement in the kinetic properties compared to the reference SuSy1-haplotype#1 containing Arg-16. Meanwhile SuSy3-haplotype#9 from a wild accession, containing four amino acid changes S53A, S106I, E727D and K741E, showed an increase in Vmax/Km at 16°C compared to the reference SuSy3-haplotype#1. This study demonstrates that SuSy kinetic properties can be enhanced by exploiting natural variation, and the potential of this enzyme to improve sucrose metabolism and eventually sink strength in planta., Competing Interests: The authors have declared that no competing interests exist.

Background and Aims: Pathogenic mutations in the Low Density Lipoprotein Receptor gene (LDLR) cause Familial Hypercholesterolemia (FH), one of the most common genetic disorders with a prevalence as high as 1 in 200 in some populations. FH is an autosomal dominant disorder of lipoprotein metabolism characterized by high blood cholesterol levels, deposits of cholesterol in peripheral tissues such as tendon xanthomas and accelerated atherosclerosis. To date, 2500 LDLR variants have been identified in the LDLR gene; however, only a minority of them have been experimentally characterized and proven to be pathogenic. Here we investigated the role of Cys46 located in the first repeat of the LDL receptor binding domain in recognition of apolipoproteins., Methods: Activity of the p.(Cys46Gly) LDLR variant was assessed by immunoblotting and flow cytometry in CHO-ldlA7 expressing the receptor variant. Affinity of p.(Cys46Gly) for LDL and VLDL was determined by solid-phase immunoassays and in silico analysis was used to predict mutation effects., Results and Conclusion: Functional characterization of p.(Cys46Gly) LDLR variant showed impaired LDL and VLDL binding and uptake activity. Consistent with this, solid-phase immunoassays showed the p.(Cys46Gly) LDLR variant has decreased binding affinity for apolipoproteins. These results indicate the important role of Cys46 in LDL receptor activity and highlight the role of LR1 in LDLr activity modulation. This study reinforces the significance of in vitro functional characterization of LDL receptor activity in developing an accurate approach to FH genetic diagnosis. This is of particular importance because it enables clinicians to tailor personalized treatments for patients' mutation profile., Competing Interests: Progenika Biopharma SA has SEQPRO LIPO® in the market to detect FH causing variants and does not have any role with patents or products in development related to this work. This work will improve studied variant reporting. Since all the studied data in this study is published, no competing conflicts exist. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

P-glycoprotein (P-gp) is an ABC transporter that exports many amphipathic or hydrophobic compounds, including chemically and functionally dissimilar anticancer drugs, from cells. To understand the role of transmembrane helices (TMH) 1 and 7 in drug-binding and transport, we selected six residues from both TMH1 (V53, I59, I60, L65, M68 and F72) and TMH7 (V713, I719, I720, Q725, F728 and F732); and substituted them with alanine by gene synthesis to generate a variant termed "TMH1,7 mutant P-gp". The expression and function of TMH1,7 mutant P-gp with twelve mutations was characterized using the BacMam baculovirus-HeLa cell expression system. The expression and conformation of TMH1,7 mutant P-gp was not altered by the introduction of the twelve mutations, as confirmed by using the human P-gp-specific antibodies UIC2, MRK16 and 4E3. We tested 25 fluorescently-labeled substrates and found that only three substrates, NBD-cyclosporine A, Rhod-2-AM and X-Rhod-1-AM were transported by the TMH1,7 mutant. The basal ATPase activity of TMH1,7 mutant P-gp was lower (40-50%) compared to wild-type (WT) P-gp, despite similar level of expression. Although most of the substrates modulate ATPase activity of P-gp, the activity of TMH1,7 mutant transporter was not significantly modulated by any of the tested substrates. Docking of selected substrates in homology models showed comparable docking scores for the TMH1,7 mutant and WT P-gp, although the binding conformations were different. Both the ATPase assay and in silico docking analyses suggest that the interactions with residues in the drug-binding pocket are altered as a consequence of the mutations. We demonstrate that it is possible to generate a variant of P-gp with a loss of broad substrate specificity and propose that TMH1 and TMH7 play a critical role in the drug efflux function of this multidrug transporter., Competing Interests: The authors have declared that no competing interests exist.

The active form of vitamin B6, pyridoxal 5'-phosphate (PLP), plays an essential role in the catalytic mechanism of various proteins, including human glutamate-oxaloacetate transaminase (hGOT1), an important enzyme in amino acid metabolism. A recent molecular and genetic study showed that the E266K, R267H, and P300L substitutions in aspartate aminotransferase, the Arabidopsis analog of hGOT1, genetically suppress a developmentally arrested Arabidopsis RUS mutant. Furthermore, CD analyses suggested that the variants exist as apo proteins and implicated a possible role of PLP in the regulation of PLP homeostasis and metabolic pathways. In this work, we assessed the stability of PLP bound to hGOT1 for the three variant and wildtype (WT) proteins using a combined 6 μs of molecular dynamics (MD) simulation. For the variants and WT in the holo form, the MD simulations reproduced the "closed-open" transition needed for substrate binding. This conformational transition was associated with the rearrangement of the P15-R32 small domain loop providing substrate access to the R387/R293 binding motif. We also showed that formation of the dimer interface is essential for PLP affinity to the active site. The position of PLP in the WT binding site was stabilized by a unique hydrogen bond network of the phosphate binding cup, which placed the cofactor for formation of the covalent Schiff base linkage with K259 for catalysis. The amino acid substitutions at positions 266, 267, and 300 reduced the structural correlation between PLP and the protein active site and/or integrity of the dimer interface. Principal component analysis and energy decomposition clearly suggested dimer misalignment and dissociation for the three variants tested in our work. The low affinity of PLP in the hGOT1 variants observed in our computational work provided structural rationale for the possible role of vitamin B6 in regulating metabolic pathways., Competing Interests: The authors have declared that no competing interests exist.

Transcription of a gene can be regulated at many different levels. One such fundamental level is interaction between protein and DNA. Protein(s) binds to distinct sites on the DNA, which activate, enhance or repress transcription. Despite being such an important process, very few tools exist to identify the proteins that interact with chromosome, most of which are in vitro in nature. Here, we propose an in vivo based method for identification of DNA binding protein(s) in bacteria where the DNA-protein complex formed in vivo is crosslinked by formaldehyde. This complex is further isolated and the bound proteins are identified. The method was used to isolate promoter DNA binding proteins, which bind and regulate the dsz operon in Gordonia sp. IITR 100 responsible for biodesulfurization of organosulfurs. The promoter binding proteins were identified by MALDI ToF MS/MS and the binding was confirmed by gel shift assay. Unlike other reported in vivo methods, this improved method does not require sequence of the whole genome or a chip and can be scaled up to improve yields., Competing Interests: The authors have declared that no competing interests exist.

Herpes Simplex Virus type 2 (HSV-2) is a neurotropic human pathogen. Upon de novo infection, the viral infected cell protein 0 (ICP0) is immediately expressed and interacts with various cellular components during the viral replication cycle. ICP0 is a multifunctional regulatory protein that has been shown to be important for both efficient viral replication and virus reactivation from latency. In particular, as previously demonstrated in transfected tissue culture models, ICP0 interacts with the cellular E3 ubiquitin ligase SIAH-1, which targets ICP0 for proteasomal degradation. However, the consequence of this virus-host interaction during the establishment of HSV-2 infection in vivo has not yet been elucidated. Here we confirmed that ICP0 of HSV-2 interacts with SIAH-1 via two conserved PxAxVxP amino acid binding motifs. We also demonstrate in vitro that a SIAH-1 binding-deficient HSV-2 strain, constructed by homologous recombination technology, exhibits an attenuated growth curve and impaired DNA and protein synthesis. This attenuated phenotype was also confirmed in an in vivo ocular infection mouse model. Specifically, viral load of the SIAH-1 binding-deficient HSV-2 mutant was significantly reduced in the trigeminal ganglia and brain stem at day 5 and 7 post infection. Our findings indicate that the interplay between ICP0 and SIAH-1 is important for efficient HSV-2 replication in vivo, thereby affecting viral dissemination kinetics in newly infected organisms, and possibly revealing novel targets for antiviral therapy., Competing Interests: The commercial affiliation of an author of this study [MME] (Biomedizinische Forschungsgesellschaft mbH, Vienna, Austria) does not alter our adherence to all PLOS ONE policies on sharing data and materials.

Iron-sulphur (Fe-S) clusters are ubiquitous co-factors which require multi-protein systems for their synthesis. In Mycobacterium tuberculosis, the Rv1460-Rv1461-Rv1462-Rv1463-csd-Rv1465-Rv1466 operon (suf operon) encodes the primary Fe-S cluster biogenesis system. The first gene in this operon, Rv1460, shares homology with the cyanobacterial SufR, which functions as a transcriptional repressor of the sufBCDS operon. Rv1460's function in M. tuberculosis has however not been determined. In this study, we demonstrate that M. tuberculosis mutants lacking a functional Rv1460 protein are impaired for growth under standard culture conditions. Elevated expression of Rv1460 and Rv1461 was observed in the mutant, implicating Rv1460 in the regulation of the suf operon. Binding of an Fe-S cluster to purified recombinant Rv1460 was confirmed by UV-visible spectroscopy and circular dichroism. Furthermore, three conserved cysteine residues, C203, C216 and C244, proposed to provide ligands for the coordination of an Fe-S cluster, were shown to be required for the function of Rv1460 in M. tuberculosis. Rv1460 therefore seems to be functionally analogous to cyanobacterial SufR., Competing Interests: The authors have declared that no competing interests exist.

D-amino acid oxidase (DAAO) degrades D-amino acids to produce α-ketoacids, hydrogen peroxide and ammonia. DAAO has often been investigated and engineered for industrial and clinical applications. We combined information from literature with a detailed analysis of the structure to engineer mammalian DAAOs. The structural analysis was complemented with molecular dynamics simulations to characterize solvent accessibility and product release mechanisms. We identified non-obvious residues located on the loops on the border between the active site and the secondary binding pocket essential for pig and human DAAO substrate specificity and activity. We engineered DAAOs by mutating such critical residues and characterised the biochemical activity of the resulting variants. The results highlight the importance of the selected residues in modulating substrate specificity, product egress and enzyme activity, suggesting further steps of DAAO re-engineering towards desired clinical and industrial applications., Competing Interests: The authors have declared that no competing interests exist.

Scabin is a mono-ADP-ribosyltransferase enzyme and is a putative virulence factor produced by the plant pathogen, Streptomyces scabies. Previously, crystal structures of Scabin were solved in the presence and absence of substrate analogues and inhibitors. Herein, experimental (hydrogen-deuterium exchange), simulated (molecular dynamics), and theoretical (Gaussian Network Modeling) approaches were systematically applied to study the dynamics of apo-Scabin in the context of a Scabin·NAD+·DNA model. MD simulations revealed that the apo-Scabin solution conformation correlates well with the X-ray crystal structure, beyond the conformation of the exposed, mobile regions. In turn, the MD fluctuations correspond with the crystallographic B-factors, with the fluctuations derived from a Gaussian network model, and with the experimental H/D exchange rates. An Essential Dynamics Analysis identified the dynamic aspects of the toxin as a crab-claw-like mechanism of two topological domains, along with coupled deformations of exposed motifs. The "crab-claw" movement resembles the motion of C3-like toxins and emerges as a property of the central β scaffold of catalytic single domain toxins. The exposure and high mobility of the cis side motifs in the Scabin β-core suggest involvement in DNA substrate binding. A ternary Scabin·NAD+·DNA model was produced via an independent docking methodology, where the intermolecular interactions correspond to the region of high mobility identified by dynamics analyses and agree with binding and kinetic data reported for wild-type and Scabin variants. Based on data for the Pierisin-like toxin group, the sequence motif Rβ1-RLa-NLc-STTβ2-WPN-WARTT-(QxE)ARTT emerges as a catalytic signature involved in the enzymatic activity of these DNA-acting toxins. However, these results also show that Scabin possesses a unique DNA-binding motif within the Pierisin-like toxin group.

Long noncoding RNAs (lncRNAs) are a large class of regulatory RNAs with diverse roles in cellular processes. Thousands of lncRNAs have been discovered; however, their roles in the regulation of muscle differentiation are unclear because no comprehensive analysis of lncRNAs during this process has been performed. In the present study, by combining diverse RNA sequencing datasets obtained from public database, we discovered lncRNAs that could behave as regulators in the differentiation of smooth or skeletal muscle cells. These analyses confirmed the roles of previously reported lncRNAs in this process. Moreover, we discovered dozens of novel lncRNAs whose expression patterns suggested their possible involvement in the phenotypic switch of vascular smooth muscle cells. The comparison of lncRNA expression change suggested that many lncRNAs have common roles during the differentiation of smooth and skeletal muscles, while some lncRNAs may have opposite roles in this process. The expression change of lncRNAs was highly correlated with that of their neighboring genes, suggesting that they may function as cis-acting lncRNAs. Furthermore, within the lncRNA sequences, there were binding sites for miRNAs with expression levels inversely correlated with the expression of corresponding lncRNAs during differentiation, suggesting a possible role of these lncRNAs as competing endogenous RNAs. The lncRNAs identified in this study will be a useful resource for future studies of gene regulation during muscle differentiation.

LapA is the biggest protein in Pseudomonas putida and a key factor for biofilm formation. Its importance and posttranslational regulation is rather thoroughly studied but less is known about the transcriptional regulation. Here we give evidence that transcription of lapA in LB-grown bacteria is initiated from six promoters, three of which display moderate RpoS-dependence. The global transcription regulator Fis binds to the lapA promoter area at six positions in vitro, and Fis activates the transcription of lapA while overexpressed in cells. Two of the six Fis binding sites, Fis-A7 and Fis-A5, are necessary for the positive effect of Fis on the transcription of lapA in vivo. Our results indicate that Fis binding to the Fis-A7 site increases the level of transcription from the most distal promoter of lapA, whereas Fis binding to the Fis-A5 site could be important for modifying the promoter area topology.

The identification of metal ion binding sites is important for protein function annotation and the design of new drug molecules. This study presents an effective method of analyzing and identifying the binding residues of metal ions based solely on sequence information. Ten metal ions were extracted from the BioLip database: Zn2+, Cu2+, Fe2+, Fe3+, Ca2+, Mg2+, Mn2+, Na+, K+ and Co2+. The analysis showed that Zn2+, Cu2+, Fe2+, Fe3+, and Co2+ were sensitive to the conservation of amino acids at binding sites, and promising results can be achieved using the Position Weight Scoring Matrix algorithm, with an accuracy of over 79.9% and a Matthews correlation coefficient of over 0.6. The binding sites of other metals can also be accurately identified using the Support Vector Machine algorithm with multifeature parameters as input. In addition, we found that Ca2+ was insensitive to hydrophobicity and hydrophilicity information and Mn2+ was insensitive to polarization charge information. An online server was constructed based on the framework of the proposed method and is freely available at http://60.31.198.140:8081/metal/HomePage/HomePage.html.

Growth differentiation factor 15 (GDF15) is a strong predictor of cardiovascular events and mortality in individuals with or without cardiovascular diseases. Single nucleotide polymorphisms (SNPs) in microRNA (miRNA) target sites, also known as miRSNPs, are known to enhance or weaken miRNA-mRNA interactions and have been linked to diseases such as cardiovascular disease and cancer. In this study, we aimed to elucidate the functional significance of the miRSNP rs1054564 in regulating GDF15 levels. Two rs1054564-containing binding sites for hsa-miR-873-5p and hsa-miR-1233-3p were identified in the 3' untranslated region (UTR) of the GDF15 transcript using bioinformatics tools. Their activities were further characterized by in vitro reporter assays. Bioinformatics prediction suggested that miRNA binding sites harboring the rs1054564-G allele had lower free energies than those with the C allele and therefore were better targets with higher affinities for both hsa-miR-873-5p and hsa-miR-1233-3p. Reporter assays showed that luciferase activity was significantly decreased by rs1054564-G-containing 3' UTRs for both miRNAs (P < 0.05) and was restored by miRNA inhibitors. Comparing the fold suppression of the two miRNAs, only that of hsa-miR-1233-3p showed significant changes between the rs1054564-G- and C-containing 3' UTRs (P = 0.034). In addition, western blots showed that transfection of both miRNA mimics significantly decreased endogenous GDF15 expression in a melanoma cell line (P < 0.05). Taken together, our findings demonstrate that GDF15 is a target of hsa-miR-873-5p and hsa-miR-1233-3p and that the rs1054564-C allele partially abolishes hsa-miR-1233-3p-mediated translational suppression of GDF15. These results suggest that rs1054564 confers allele-specific translational repression of GDF15 via hsa-miR-1233-3p. Our work thus provides biological insight into the previously reported clinical association between rs1054564 and plasma GDF15 levels.

V-set and transmembrane domain-containing protein 5 (Vstm5), a newly characterized small membrane glycoprotein, can induce membrane protrusions in various cells. Vstm5 can modulate both the position and complexity of central neurons by altering their membrane morphology and dynamics. In this study, we investigated the significance of glycosylation in the expression and function of Vstm5. Four N-linked glycosylation sites (Asn43, Asn87, Asn101, and Asn108) are predicted to be located in the extracellular N-terminus of mouse Vstm5. Although all four sites were glycosylated, their functional roles may not be identical. N-glycosylation at multiple sites affects differentially the function of Vstm5. Glycosylation at individual sites not only played essential roles in surface expression of Vstm5 but also in the formation of neuronal dendritic filopodia. These results indicate that N-linked glycosylation at multiple sites plays important roles by differentially influencing the expression, targeting, and biological activity of Vstm5.

Species A rotavirus non-structural protein 3 (NSP3) is a translational regulator that inhibits or, under some conditions, enhances host cell translation. NSP3 binds to the translation initiation factor eIF4G1 and evicts poly-(A) binding protein (PABP) from eIF4G1, thus inhibiting translation of polyadenylated mRNAs, presumably by disrupting the effect of PABP bound to their 3'-ends. NSP3 has a long coiled-coil region involved in dimerization that includes a chaperone Hsp90-binding domain (HS90BD). We aimed to study the role in NSP3 dimerization of a segment of the coiled-coil region adjoining the HS90BD. We used a vaccinia virus system to express NSP3 with point mutations in conserved amino acids in the coiled-coil region and determined the effects of these mutations on translation by metabolic labeling of proteins as well as on accumulation of stable NSP3 dimers by non-dissociating Western blot, a method that separates stable NSP3 dimers from the monomer/dimerization intermediate forms of the protein. Four of five mutations reduced the total yield of NSP3 and the formation of stable dimers (W170A, K171E, R173E and R187E:K191E), whereas one mutation had the opposite effects (Y192A). Treatment with the proteasome inhibitor MG132 revealed that stable NSP3 dimers and monomers/dimerization intermediates are susceptible to proteasome degradation. Surprisingly, mutants severely impaired in the formation of stable dimers were still able to inhibit host cell translation, suggesting that NSP3 dimerization intermediates are functional. Our results demonstrate that rotavirus NSP3 acquires its function prior to stable dimer formation and remain as a proteasome target throughout dimerization.

Similarly to its homologs in other bacteria, Pseudomonas aeruginosa partitioning protein ParB facilitates segregation of newly replicated chromosomes. Lack of ParB is not lethal but results in increased frequency of anucleate cells production, longer division time, cell elongation, altered colony morphology and defective swarming and swimming motility. Unlike in other bacteria, inactivation of parB leads to major changes of the transcriptome, suggesting that, directly or indirectly, ParB plays a role in regulation of gene expression in this organism. ParB overproduction affects growth rate, cell division and motility in a similar way as ParB deficiency. To identify primary ParB targets, here we analysed the impact of a slight increase in ParB level on P. aeruginosa transcriptome. ParB excess, which does not cause changes in growth rate and chromosome segregation, significantly alters the expression of 176 loci. Most notably, the mRNA level of genes adjacent to high affinity ParB binding sites parS1-4 close to oriC is reduced. Conversely, in cells lacking either parB or functional parS sequences the orfs adjacent to parS3 and parS4 are upregulated, indicating that direct ParB- parS3/parS4 interactions repress the transcription in this region. In addition, increased ParB level brings about repression or activation of numerous genes including several transcriptional regulators involved in SOS response, virulence and adaptation. Overall, our data support the role of partitioning protein ParB as a transcriptional regulator in Pseudomonas aeruginosa.

Common wild rice (Oryza rufipogon Griff.) is an important germplasm for rice breeding, which contains many resistance genes. Re-sequencing provides an unprecedented opportunity to explore the abundant useful genes at whole genome level. Here, we identified the nucleotide-binding site leucine-rich repeat (NBS-LRR) encoding genes by re-sequencing of two wild rice lines (i.e. Huaye 1 and Huaye 2) that were developed from common wild rice. We obtained 128 to 147 million reads with approximately 32.5-fold coverage depth, and uniquely covered more than 89.6% (> = 1 fold) of reference genomes. Two wild rice lines showed high SNP (single-nucleotide polymorphisms) variation rate in 12 chromosomes against the reference genomes of Nipponbare (japonica cultivar) and 93-11 (indica cultivar). InDels (insertion/deletion polymorphisms) count-length distribution exhibited normal distribution in the two lines, and most of the InDels were ranged from -5 to 5 bp. With reference to the Nipponbare genome sequence, we detected a total of 1,209,308 SNPs, 161,117 InDels and 4,192 SVs (structural variations) in Huaye 1, and 1,387,959 SNPs, 180,226 InDels and 5,305 SVs in Huaye 2. A total of 44.9% and 46.9% genes exhibited sequence variations in two wild rice lines compared to the Nipponbare and 93-11 reference genomes, respectively. Analysis of NBS-LRR mutant candidate genes showed that they were mainly distributed on chromosome 11, and NBS domain was more conserved than LRR domain in both wild rice lines. NBS genes depicted higher levels of genetic diversity in Huaye 1 than that found in Huaye 2. Furthermore, protein-protein interaction analysis showed that NBS genes mostly interacted with the cytochrome C protein (Os05g0420600, Os01g0885000 and BGIOSGA038922), while some NBS genes interacted with heat shock protein, DNA-binding activity, Phosphoinositide 3-kinase and a coiled coil region. We explored abundant NBS-LRR encoding genes in two common wild rice lines through genome wide re-sequencing, which proved to be a useful tool to exploit elite NBS-LRR genes in wild rice. The data here provide a foundation for future work aimed at dissecting the genetic basis of disease resistance in rice, and the two wild rice lines will be useful germplasm for the molecular improvement of cultivated rice.

Kif5A is a neuronally-enriched isoform of the Kinesin-1 family of cellular transport motors. 23 separate mutations in the motor domain of Kif5A have been identified in patients with the complicated form of hereditary spastic paraplegia (HSP). We performed in vitro assays on dimeric recombinant Kif5A with HSP-causing mutations in the Switch I domain, which participates in the coordination and hydrolysis of ATP by kinesin. We observed a variety of significantly reduced catalytic and mechanical activities as a result of each mutation, with the shared phenotype from each that motility was significantly reduced. Substitution of Mn2+ for Mg2+ in our reaction buffers provides a dose-dependent rescue in both the catalytic and ensemble mechanical properties of the S203C mutant. This work provides mechanistic insight into the cause of HSP in patients with these mutations and points to future experiments to further dissect the root cause of this disease.

Nitrile hydratase (NHase) from Rhodococcus rhodochrous J1 is widely used for industrial production of acrylamide and nicotinamide. However, the two types of NHases (L-NHase and H-NHase) from R. rhodochrous J1 were only slightly expressed in E. coli by routine methods, which limits the comprehensive and systematic characterization of the enzyme properties. We successfully expressed the two types of recombinant NHases in E. coli by codon-optimization, engineering of Ribosome Binding Site (RBS) and spacer sequences. The specific activity of the purified L-NHase and H-NHase were 400 U/mg and 234 U/mg, respectively. The molecular mass of L-NHase and H-NHase was identified to be 94 kDa and 504 kDa, respectively, indicating that the quaternary structure of the two types of NHases was the same as those in R. rhodochrous J1. H-NHase exhibited higher substrate and product tolerance than L-NHase. Moreover, higher activity and shorter culture time were achieved in recombinant E. coli, and the whole cell catalyst of recombinant E. coli harboring H-NHase has equivalent efficiency in tolerance to the high-concentration product relative to that in R. rhodochrous J1. These results indicate that biotransformation of nitrile by R. rhodochrous J1 represents a potential alternative to NHase-producing E. coli.

Chronic kidney disease (CKD) has a prevalence of approximately 10% in adult populations. CKD can progress to end-stage renal disease (ESRD) and this is usually fatal unless some form of renal replacement therapy (chronic dialysis or renal transplantation) is provided. There is an inherited predisposition to CKD with several genetic risk markers now identified. The UMOD gene has been associated with CKD of varying aetiologies. An AmpliSeq next generation sequencing panel was developed to facilitate comprehensive sequencing of the UMOD gene, covering exonic and regulatory regions. SNPs and CpG sites in the genomic region encompassing UMOD were evaluated for association with CKD in two studies; the UK Wellcome Trust Case-Control 3 Renal Transplant Dysfunction Study (n = 1088) and UK-ROI GENIE GWAS (n = 1726). A technological comparison of two Ion Torrent machines revealed 100% allele call concordance between S5 XL™ and PGM™ machines. One SNP (rs183962941), located in a non-coding region of UMOD, was nominally associated with ESRD (p = 0.008). No association was identified between UMOD variants and estimated glomerular filtration rate. Analysis of methylation data for over 480,000 CpG sites revealed differential methylation patterns within UMOD, the most significant of these was cg03140788 p = 3.7 x 10-10.

In mammals, very long chain fatty acids (VLCFAs) perform pleiotropic roles in a wide range of biological processes, such as cell membrane formation, cell signal transduction, and endocrine regulation. Beef and milk are abundant of palmitic acid which can be further elongated into stearic acid for synthesizing VLCFAs. Elongase of very long chain fatty acid 6 (ELOVL6) is a rate-limiting enzyme for converting palmitic acid to stearic acid. Consequently, investigating the tissue expression patterns and transcriptional regulation of bovine ELOVL6 can provide new insights into improving the composition of beneficial fats in cattle and expanding the knowledge of transcriptional regulation mechanism among domestic animals. In the current study, we found that bovine ELOVL6 expressed ubiquitously. Dual-luciferase reporter assay identified that the core promoter region (-130/-41 bp) was located in the second CpG island. In addition, the deletion mutation of binding sites demonstrated that sterol regulatory element binding transcription factor 1 (SREBF1) and specific protein 1 (SP1) both were able to stimulate bovine ELOVL6 promoter activity independently, while resulting the similar effect. To confirm these findings, further RNA interference assays were executed in bovine mammary epithelial cells (BMECs). In summary, these data suggest that bovine ELOVL6 expressed ubiquitously and is activated by SREBF1 and SP1, via two binding sites present in the ELOVL6 promoter region between -130 bp to -41bp.

The psychrophilic and mesophilic endonucleases A (EndA) from Aliivibrio salmonicida (VsEndA) and Vibrio cholera (VcEndA) have been studied experimentally in terms of the biophysical properties related to thermal adaptation. The analyses of their static X-ray structures was no sufficient to rationalize the determinants of their adaptive traits at the molecular level. Thus, we used Molecular Dynamics (MD) simulations to compare the two proteins and unveil their structural and dynamical differences. Our simulations did not show a substantial increase in flexibility in the cold-adapted variant on the nanosecond time scale. The only exception is a more rigid C-terminal region in VcEndA, which is ascribable to a cluster of electrostatic interactions and hydrogen bonds, as also supported by MD simulations of the VsEndA mutant variant where the cluster of interactions was introduced. Moreover, we identified three additional amino acidic substitutions through multiple sequence alignment and the analyses of MD-based protein structure networks. In particular, T120V occurs in the proximity of the catalytic residue H80 and alters the interaction with the residue Y43, which belongs to the second coordination sphere of the Mg2+ ion. This makes T120V an amenable candidate for future experimental mutagenesis., Competing Interests: The authors have declared that no competing interests exist.

Coding region determinant-binding protein (CRD-BP) binds to the 3'-UTR of microphthalmia-associated transcription factor (MITF) mRNA to prevent its targeted degradation by miR-340. Here, we aim to further understand the molecular interaction between CRD-BP and MITF RNA. Using point mutation in the GXXG motif of each KH domains, we showed that all four KH domains of CRD-BP are important for their physical association with MITF RNA. We mapped the CRD-BP-binding site in the 3'-UTR of MITF RNA from nts 1330-1740 and showed that the 49-nt fragment 1621-1669 is the minimal size MITF RNA for binding. Upon deletion of nts 1621-1669 within the nts1550-1740 of MITF RNA, there was a 3-fold increase in dissociation constant Kd, which further confirms the critical role sequences within nts 1621-1669 in binding to CRD-BP. Amongst the eight antisense oligonucleotides designed against MITF RNA 1550-1740, we found MHO-1 and MHO-7 as potent inhibitors of the CRD-BP-MITF RNA interaction. Using RNase protection and fluorescence polarization assays, we showed that both MHO-1 and MHO-7 have affinity for the MITF RNA, suggesting that both antisense oligonucleotides inhibited CRD-BP-MITF RNA interaction by directly binding to MITF RNA. The new molecular insights provided in this study have important implications for understanding the oncogenic function of CRD-BP and development of specific inhibitors against CRD-BP-MITF RNA interaction., Competing Interests: The authors have declared that no competing interests exist.

In sea urchin development, structures derived from neurogenic territory control the swimming and feeding responses of the pluteus as well as the process of metamorphosis. We have previously isolated an alpha tubulin family member of Paracentrotus lividus (Pl-Tuba1a, formerly known as Pl-Talpha2) that is specifically expressed in the ciliary band and animal pole neurogenic domains of the sea urchin embryo. In order to identify cis-regulatory elements controlling its spatio-temporal expression, we conducted gene transfer experiments, transgene deletions and site specific mutagenesis. Thus, a genomic region of about 2.6 Kb of Pl-Tuba1a, containing four Interspecifically Conserved Regions (ICRs), was identified as responsible for proper gene expression. An enhancer role was ascribed to ICR1 and ICR2, while ICR3 exerted a pivotal role in basal expression, restricting Tuba1a expression to the proper territories of the embryo. Additionally, the mutation of the forkhead box consensus sequence binding site in ICR3 prevented Pl-Tuba1a expression., Competing Interests: The authors have declared that no competing interests exist.

Voltage-gated sodium channels (VGSC) regulate neuronal excitability by governing action potential (AP) generation and propagation. Recent studies have revealed that AMP-activated protein kinase (AMPK) activators decrease sensory neuron excitability, potentially by preventing sodium (Na+) channel phosphorylation by kinases such as ERK or via modulation of translation regulation pathways. The direct positive allosteric modulator A769662 displays substantially greater efficacy than other AMPK activators in decreasing sensory neuron excitability suggesting additional mechanisms of action. Here, we show that A769662 acutely inhibits AP firing stimulated by ramp current injection in rat trigeminal ganglion (TG) neurons. PT1, a structurally dissimilar AMPK activator that reduces nerve growth factor (NGF) -induced hyperexcitability, has no influence on AP firing in TG neurons upon acute application. In voltage-clamp recordings, application of A769662 reduces VGSC current amplitudes. These findings, based on acute A769662 application, suggest a direct channel blocking effect. Indeed, A769662 dose-dependently blocks VGSC in rat TG neurons and in Nav1.7-transfected cells with an IC50 of ~ 10 μM. A769662 neither displayed use-dependent inhibition nor interacted with the local anesthetic (LA) binding site. Popliteal fossa administration of A769662 decreased noxious thermal responses with a peak effect at 5 mins demonstrating an analgesic effect. These data indicate that in addition to AMPK activation, A769662 acts as a direct blocker/modulator of VGSCs, a potential mechanism enhancing the analgesic property of this compound., Competing Interests: The authors have declared that no competing interests exist.

Expression of the tissue-type plasminogen activator gene (t-PA; gene name PLAT) is regulated, in part, by epigenetic mechanisms. We investigated the relationship between PLAT methylation and PLAT expression in five primary human cell types and six transformed cell lines. CpG methylation was analyzed in the proximal PLAT gene promoter and near the multihormone responsive enhancer (MHRE) -7.3 kilobase pairs upstream of the PLAT transcriptional start site (TSS, -7.3 kb). In Bowes melanoma cells, the PLAT promoter and the MHRE were fully unmethylated and t-PA secretion was extremely high. In other cell types the region from -647 to -366 was fully methylated, whereas an unmethylated stretch of DNA from -121 to +94 was required but not sufficient for detectable t-PA mRNA and t-PA secretion. DNA methylation near the MHRE was not correlated with t-PA secretion. Specific methylation of the PLAT promoter region -151 to +151, inserted into a firefly luciferase reporter gene, abolished reporter gene activity. The region -121 to + 94 contains two well-described regulatory elements, a PMA-responsive element (CRE) near -106 and a GC-rich region containing an Sp1 binding site near +59. Methylation of double-stranded DNA oligonucleotides containing the CRE or the GC-rich region had little or no effect on transcription factor binding. Methylated CpGs may attract co-repressor complexes that contain histone deacetylases (HDAC). However, reporter gene activity of methylated plasmids was not restored by the HDAC inhibitor trichostatin. In conclusion, efficient PLAT gene expression requires a short stretch of unmethylated CpG sites in the proximal promoter., Competing Interests: The authors have declared that no competing interests exist.