Nucleosomal positioning and genetic divergence study based on DNA flexibility map

Based on worm like chain model, DNA structural parameters--tilt, roll and rise, derived from crystallographic database have been used to determine the flexibility of DNA that regulates the nucleosomal translational positioning. Theoretically derived data has been compared to the experimental values available in loshikhes and Trifonov's database. The methodology has been extended to determine the flexibility of 18S rRNA genome in eukarya, where yeast shows a distinct difference when compared with mammals like human, mouse and rabbit.     

牛痘病毒感染人巨噬细胞基因表达谱的RNA-Seq分析

[目的]利用RNA-Seq,分析人巨噬细胞在牛痘病毒(VACV)感染前后基因表达的变化,探索牛痘病毒与宿主细胞相互作用的机制.[方法]用牛痘病毒感染人巨噬细胞,采用RNA-Seq比较感染组和对照组的差异表达基因,并进行KEGG、GO以及STRING网络分析.[结果]感染组与对照组相比,筛选出显著性差异表达基因4796个...     

Ab initio Investigation of the Molecular Structure of Methyl Methoxymethyl Phosphonate,a Promising Nuclease-resistant Alternative of the Phosphodiester Linkage

Abstract

Conformational flexibility of the methyl methoxymethyl phosphonate anion (CH₃-O-PO₂- CH₂-O-CH₃)^?, a nuclease resistant alternative to the phosphodiester linkage in DNA, have been investigated by ab initio quantum mechanical calculations. The potential of backbone torsional degrees of freedom of methyl methoxymethyl phosphonate anion (MMP) was determined at the Hartree-Fock (HF) 3–21G* level using the adiabatic mapping technique. Energies, geometries, and effective atomic charges of different conformers were calculated at HF/6–31G* and MP2/6–31G* levels of theory. These were compared to the results obtained for dimethyl phosphate calculated at the same level. The impact on DNA structure from inserting a methylene group between phosphorus and oxygen of the nucleoside sugar moiety was examined via distance- and angle-constrained geometry optimizations. Due to its high flexibility, MMP has been shown to be compatible with both A and B forms of DNA.     

Identification of differentially expressed genes and biological characteristics of colorectal cancer by integrated bioinformatics analysis

Colorectal cancer (CRC) ranks as one of the most common malignant tumors worldwide. Its mortality rate has remained high in recent years. Therefore, the aim of this study was to identify significant differentially expressed genes (DEGs) involved in its pathogenesis, which may be used as novel biomarkers or potential therapeutic targets for CRC. The gene expression profiles of GSE21510, GSE32323, GSE89076, and GSE113513 were downloaded from the Gene Expression Omnibus (GEO) database. After screening DEGs in each GEO data set, we further used the robust rank aggregation method to identify 494 significant DEGs including 212 upregulated and 282 downregulated genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed by DAVID and the KOBAS online database, respectively. These DEGs were shown to be significantly enriched in different cancer-related functions and pathways. Then, the STRING database was used to construct the protein–protein interaction network. The module analysis was performed by the MCODE plug-in of Cytoscape based on the whole network. We finally filtered out seven hub genes by the cytoHubba plug-in, including PPBP, CCL28, CXCL12, INSL5, CXCL3, CXCL10, and CXCL11. The expression validation and survival analysis of these hub genes were analyzed based on The Cancer Genome Atlas database. In conclusion, the robust DEGs associated with the carcinogenesis of CRC were screened through the GEO database, and integrated bioinformatics analysis was conducted. Our study provides reliable molecular biomarkers for screening and diagnosis, prognosis as well as novel therapeutic targets for CRC.     

Identification of differentially expressed genes related to metabolic syndrome induced with high-fat diet in E3 rats

Understanding the genes differentially expressing in aberrant organs of metabolic syndrome (MetS) facilitates the uncovering of molecular mechanisms and the identification of novel therapeutic targets for the disease. This study aimed to identify differentially expressed genes related to MetS in livers of E3 rats with high-fat-diet-induced metabolic syndrome (HFD-MetS). E3 rats were fed with high-fat diet for 24 weeks to induce MetS. Then, suppression subtractive hybridization (SSH) technology was used to identify the genes differentially expressed between HFD-MetS and control E3 rat livers. Twenty positive recombinant clones were chosen randomly from forward subtractive library and sent to sequence. BLAST analysis in GenBank database was used to determine the property of each cDNA fragment. In total, 11 annotated genes, 3 ESTs, and 2 novel gene fragments were identified by SSH technology. The expression of four genes (Alb, Pip4k2a, Scd1, and Tf) known to be associated with MetS and other five genes (Eif1, Rnase4, Rps12, Rup2, and Tmsb4) unknown to be relevant to MetS was significantly up-regulated in the livers of HFD-MetS E3 rats compared with control rats using real-time quantitative PCR (RT-qPCR). By analyzing the correlations between the expression of these nine genes and serum concentrations of TG, Tch, HDL-C, and LDL-C, we found that there were significant positive correlations between TG and the expression of five genes (Alb, Eif1, Pip4k2a, Rps12, and Tmsb4x), Tch and three genes (Rnase4, Scd1, and Tmsb4x), and LDL-C and two genes (Rnase4 and Scd1), as well there were significant negative correlations between HDL-C and the expression of three genes (Rup2, Scd1, and Tf). This study provides important clues for unraveling the molecular mechanisms of MetS.     

Codon usage and gene expression pattern of Stenotrophomonas maltophilia R551-3 for pathogenic mode of living

Stenotrophomonas maltophilia strain R551-3 is a multiple-antibiotic-resistant opportunistic human pathogen involved in nosocomial infections. It has a widely distributed GC-rich (>66%) genome. Analysis of differential expression of the genes of this genome reveals that majority of genes belonging to highly expressed category are mostly present on lagging strand without showing any strand specific codon usage bias. Relatively small number of lowly expressed genes is equally distributed on both leading and lagging strands with a difference in codon usage pattern between them. Among several multi drug resistance genes of S. maltophilia involving lowly expressed category some are predicted as horizontally transferred. It can be inferred that horizontally transferred genes may have been imported into this genome for their pathogenic mode of living. Our study may help to modify the expression level of the target genes of this human pathogen in order to control its infection.     

The relative flexibility of B-DNA and A-RNA duplexes: database analysis

An extensive analysis of structural databases is carried out to investigate the relative flexibility of B-DNA and A-RNA duplexes in crystal form. Our results show that the general anisotropic concept of flexibility is not very useful to compare the deformability of B-DNA and A-RNA duplexes, since the flexibility patterns of B-DNA and A-RNA are quite different. In other words, ‘flexibility’ is a dangerous word for describing macromolecules, unless it is clearly defined. A few soft essential movements explain most of the natural flexibility of A-RNA, whereas many are necessary for B-DNA. Essential movements occurring in naked B-DNAs are identical to those necessary to deform DNA in DNA–protein complexes, which suggest that evolution has designed DNA–protein complexes so that B-DNA is deformed according to its natural tendency. DNA is generally more flexible, but for some distortions A-RNA is easier to deform. Local stiffness constants obtained for naked B-DNAs and DNA complexes are very close, demonstrating that global distortions in DNA necessary for binding to proteins are the result of the addition of small concerted deformations at the base-pair level. Finally, it is worth noting that in general the picture of the relative deformability of A-RNA and DNA derived from database analysis agrees very well with that derived from state of the art molecular dynamics (MD) simulations.     

The thermodynamic signature of ligand binding to histone deacetylase-like amidohydrolases is most sensitive to the flexibility in the L2-loop lining the active site pocket

BackgroundThe analysis of the thermodynamic driving forces of ligand-protein binding has been suggested to be a key component for the selection and optimization of active compounds into drug candidates. The binding enthalpy as deduced from isothermal titration calorimetry (ITC) is usually interpreted assuming single-step binding of a ligand to one conformation of the target protein. Although successful in many cases, these assumptions are oversimplified approximations of the reality with flexible proteins and complicated binding mechanism in many if not most cases. The relationship between protein flexibility and thermodynamic signature of ligand binding is largely understudied.MethodsDirected mutagenesis, X-ray crystallography, enzyme kinetics and ITC methods were combined to dissect the influence of loop flexibility on the thermodynamics and mechanism of ligand binding to histone deacetylase (HDAC)-like amidohydrolases.ResultsThe general ligand-protein binding mechanism comprises an energetically demanding gate opening step followed by physical binding. Increased flexibility of the L2-loop in HDAC-like amidohydrolases facilitates access of ligands to the binding pocket resulting in predominantly enthalpy-driven complex formation.ConclusionsThe study provides evidence for the great importance of flexibility adjacent to the active site channel for the mechanism and observed thermodynamic driving forces of molecular recognition in HDAC like enzymes.General significanceThe flexibility or malleability in regions adjacent to binding pockets should be given more attention when designing better drug candidates. The presented case study also suggests that the observed binding enthalpy of protein-ligand systems should be interpreted with caution, since more complicated binding mechanisms may obscure the significance regarding potential drug likeness.     

Partial Genome Scale Analysis of Gene Expression in Human Adipose Tissue Using DNA Array

Objective: Large scale analysis of gene expression in adipose tissue provides a basis for the identification of novel candidate genes involved in the pathophysiology of obesity. Our goal was to explore gene expression in human adipose tissue at a partial genome scale using DNA array. Research Methods and Procedures: Labeled cDNA, derived from human adipose tissue poly(A⁺) RNA, was hybridized to a DNA array containing over 18,000 human expressed sequence‐tagged (EST) clones. The results were analyzed by database searches. Results: Homology searches of the 300 EST clones with highest hybridization signals revealed that 145 contained DNA sequences identical to known genes and 79 could be linked to UniGene clusters. Of the 145 identified genes, 136 were nonredundant and subsequently characterized with respect to function and chromosomal localization by searching MEDLINE, UniGene, GeneMap, OMIM, SWISS‐PROT, the Genome Database, and the Location Data Base. The identified genes were grouped according to their putative functions; cell/organism defense (9.6%), cell division (5.1%), cell signaling/communication (19.8%), cell structure/motility (12.5%), gene/protein expression (16.9%), metabolism (16.2%), and unclassified (19.8%). Less than 50% of these genes have previously been reported to be expressed in adipose tissue. The chromosomal localization of 268 genes strongly expressed in adipose tissue showed that their relative abundance was significantly increased on chromosomes 11, 19, and 22 compared to the expected distribution of the same number of random genes. Discussion: Our study resulted in the identification of numerous genes previously not reported to be expressed in adipose tissue. These results suggest that DNA array is a powerful tool in the search for novel regulatory pathways within adipose tissue on a scale that is not possible using conventional methods.     

Structural Property of DNA That Migrates Faster in Gel Electrophoresis,as Deduced by CD Spectroscopy

Bent DNAs are known to migrate slower than ordinary DNA in non-denaturing polyacrylamide gel electrophoresis. In contrast, several satellite DNAs have been shown to migrate fast. The structural property that causes the fast migration, however, is not clarified so far on molecular basis. We have investigated the structural property of a satellite DNA, which contains consecutive purine sequences and migrates faster in gel, by CD spectroscopy. Partial formation of an A-form–like structure has been suggested. Reduction in DNA length due to the formation of the A-form–like structure may be responsible for the fast migration. The pronounced rigidity of DNA may also contribute to the behavior.     

Binding of Non-Intercalating Antibiotics to B-DNA: A Theoretical Study Taking Into Account Nucleic Acid Flexibility

Abstract

A detailed theoretical study has been made for five antibiotics which all bind selectively to AT sequences in the minor groove of B-DNA: SN-18071, NSCT-101327, distamycin-2, distamycin-3 and netropsin. The optimal complexes were found for systems in which the flexibility of DNA, as well as that of the antibiotics, was taken into account. Explicit, mobile counterions and a dielectric function modelling aqueous solution were also included. The binding geometries of the most strongly interacting antibiotics, distamycin-3 and netropsin, are compared in considerable detail and it is shown that notable differences exist between them. The results for netropsin are also discussed in the light of recent disagreements concerning its exact binding location within DNA.     

How many differentially expressed genes: A perspective from the comparison of genotypic and phenotypic distances

Identifying differentially expressed genes is critical in microarray data analysis. Many methods have been developed by combining p-value, fold-change, and various statistical models to determine these genes. When using these methods, it is necessary to set up various pre-determined cutoff values. However, many of these cutoff values are somewhat arbitrary and may not have clear connections to biology. In this study, a genetic distance method based on gene expression level was developed to analyze eight sets of microarray data extracted from the GEO database. Since the genes used in distance calculation have been ranked by fold-change, the genetic distance becomes more stable when adding more genes in the calculation, indicating there is an optimal set of genes which are sufficient to characterize the stable difference between samples. This set of genes is differentially expressed genes representing both the genotypic and phenotypic differences between samples.     

基于高通量测序的花马盐湖原核微生物及其耐盐基因

【目的】揭示陕北花马盐湖沉积物原核微生物群落组成,并分析其潜在的耐盐功能基因。【方法】构建盐湖沉积物宏基因组16S r RNA文库和fosmid文库,利用Illumina HiSeq高通量测序及生物信息技术分析细菌古菌群落组成和耐盐菌株(5-5)外源宏基因组的潜在耐盐基因。【结果】获得18978条有效的16Sr RNA序列,共5221个OTUs,包括23个门,155个属,其中广古菌门(Euryarchaeota)和变形菌门(Proteobacteria)为优势菌门,盐杆状菌属(Halorhabdus)、盐红菌属(Halorubrum)及假单胞菌属(Pseudomonas)等16个属为优势属,以及嗜盐单胞菌属(Halomonas)、冷弯菌属(Psychroflexus)及不动细菌属(Acinetobacter)等139个属为非优势属。从4126个fosmid文库菌株中筛选出37株耐盐菌株,其中菌株5-5、2E4和2F4对不同浓度的NaCl、CuSO_4、ZnSO_4及CdSO_4具有耐受性,从5-5的外源宏基因组序列中获得61个Unigene,其中12个Unigene的同源基因编码的蛋白质如无机焦磷酸酶、转座酶、亚碲酸钾抗性蛋白及钙调蛋白等广泛参与其他生物的耐盐逆境。【结论】盐湖沉积物中蕴藏着丰富多样的细菌古菌类群以及潜在耐盐功能基因资源。     

A statistical mechanical model for predicting B-DNA curvature and flexibility

A statistical mechanical model taking into account the symmetric twisting, tilting, sliding fluctuations and asymmetric rolling fluctuations has been proposed to predict the macroscopic curvature and flexibility of B-DNA. Based on the statistical data of structural parameters of double helix in nucleic acid database and the related theoretical analysis, the equilibrium angular parameters (Omega, rho and tau) describing the orientation of successive base-pair planes, the translation parameters (D(y)) along the long axis of neighboring base-pair step and the corresponding force constants are arranged for ten dimers appropriately. Under the assumption of independent angular parameters, independent base-pair steps and a simple energy function, we can calculate the macroscopic curvature and the flexibility of DNA sequences through the transformation matrix and the Boltzmann ensemble average. The predictions on curvature and flexibility of DNA have been compared with the corresponding experimental data. The agreement is remarkably good. It is demonstrated that the lowering of the temperature does increase the DNA curvature.     

Screening of hub genes and prediction of putative drugs in arsenic-related bladder carcinoma: An in silico study

BackgroundEvidence showed that inorganic arsenic (iAs) can trigger malignant transformation in cells with complex mechanisms. Thus, we aimed to investigate the possible molecules, pathways and therapeutic drugs for iAs-induced bladder cancer (BC) by using bioinformatics approaches.MethodsMicroarray-based data were analyzed to screen the differentially expressed genes (DEGs) between iAs-related BC cells and controls. Then, the roles of DEGs were annotated and the hub genes were screened out by protein-protein interaction network. The key genes were further selected from the hub genes through an assessment of the prognostic values. Afterward, potential drugs were predicted by using CMAP analysis.ResultsAnalysis of a dataset (GSE90023) generated 21 upregulated and 47 downregulated DEGs, which were enriched in various signaling pathways. Among the DEGs, four hub genes including WNT7B, SFRP1, DNAJB2, and ATF3, were identified as the key genes because they might predict poor prognosis in BC patients. Lastly, Cantharidin was predicted to be a potential drug reversing iAs-induced malignant transformation in urinary epithelium cells.ConclusionThe present study found several hub genes involved in iAs-induced malignant transformation in urinary epithelium cells, and predicted several small agents for iAs toxicity prevention or therapy.