Scoring functions and enrichment: a case study on Hsp90

Background  

The need for fast and accurate scoring functions has been driven by the increased use of in silico virtual screening twinned with high-throughput screening as a method to rapidly identify potential candidates in the early stages of drug development. We examine the ability of some the most common scoring functions (GOLD, ChemScore, DOCK, PMF, BLEEP and Consensus) to discriminate correctly and efficiently between active and non-active compounds among a library of ~3,600 diverse decoy compounds in a virtual screening experiment against heat shock protein 90 (Hsp90).     

结直肠癌肝转移原代肿瘤组织的异种移植动物模型的研究进展

肝是结直肠癌转移主要的靶器官,结直肠癌肝转移也是结直肠癌患者死亡的主要原因。近年来研究显示,结直肠癌肝转移原代肿瘤组织的异种移植(patient-derived tumor xenograft,PDTX)动物模型能较好的复制临床肿瘤患者的特征,其建模方法主要分为异种原位与异种异位种植模型。本文综述结直肠癌肝转移PDTX动物模型造模方法及其应用范围,以期为实验建模提供参考。     

测序深度对RNA编辑识别算法的影响

RNA编辑是重要的转录后修饰过程,目前已有多种算法用于识别RNA编辑,本文主要研究小鼠中测序深度对RNA编辑识别算法的影响,从而为RNA编辑的研究给出建议的方法. 本文使用STAR比对软件将小鼠的RNA-seq数据进行序列比对,然后使用GATK识别SNV,并用Separate Method、GIREMI、RNAEditor 3种方法识别出RNA编辑位点. 最后对3种方法识别RNA编辑位点的共同部分、识别效率、识别稳定性、识别与测序深度的关系进行分析. 结果发现3种方法识别的编辑位点数目差异大,共有位点较少,随着测序深度的增加,识别的RNA编辑位点数也在增加. 结果表明RNA编辑识别算法在小鼠中的识别性能与测序深度呈正相关.     

辐射诱导的染色质拓扑关联结构域层级变化及其在细胞辐射响应中的作用

基因组三维结构在基因表达调控中发挥重要作用,染色质拓扑关联结构域(topologically associated domain,TAD)是DNA复制和基因转录的基本功能单位,也是DNA损伤修复的功能单元,在辐射诱导的DNA损伤修复中发挥重要作用。近期研究表明,TAD并非是完全独立的结构单元,其内部常呈现多层级结构,对基因表达具有重要调控作用。为探究TAD多层级结构在细胞辐射响应中的作用,本研究使用TAD层级结构识别算法OnTAD对Gene expression omnibus数据库中5Gy X射线照射的淋巴细胞、成纤维细胞和毛细血管扩张性共济失调突变(ataxia telangiectasia mutated,ATM)基因缺陷的成纤维细胞,共26个样本的Hi-C(high-through chromosome conformation capture,Hi-C)数据进行分析,发现辐射后细胞的TAD层级结构出现规律性变化,高层级TAD缺失较多,低层级TAD相对保守;辐射诱导的TAD层级结构变化通过调节基因表达参与细胞辐射响应;ATM是辐射诱导TAD层级结构变化和恢复的重要因子。本研究为从TAD多层级结构角度理解基因组三维结构在细胞辐射响应中的作用提供了新思路。     

The first laminin G-like domain of protein S is essential for binding and activation of Tyro3 receptor and intracellular signalling

The homologous proteins Gas6 and protein S (ProS1) are both natural ligands for the TAM (Tyro3, Axl, MerTK) receptor tyrosine kinases. ProS1 selectively activates Tyro3; however, the precise molecular interface of the ProS1-Tyro3 contact has not been characterised. We used a set of chimeric proteins in which each of the C-terminal laminin G-like (LG) domains of ProS1 were swapped with those of Gas6, as well as a set of ProS1 mutants with novel added glycosylations within LG1. Alongside wildtype ProS1, only the chimera containing ProS1 LG1 domain stimulated Tyro3 and Erk phosphorylation in human cancer cells, as determined by Western blot. In contrast, Gas6 and chimeras containing minimally the Gas6 LG1 domain stimulated Axl and Akt phosphorylation. We performed in silico homology modelling and molecular docking analysis to construct and evaluate structural models of both ProS1-Tyro3 and Gas6-Axl ligand-receptor interactions. These analyses revealed a contact between the ProS1 LG1 domain and the first immunoglobulin domain of Tyro3, which was similar to the Gas6-Axl interaction, and involved long-range electrostatic interactions that were further stabilised by hydrophobic and polar contacts. The mutant ProS1 proteins, which had added glycosylations within LG1 but which were all outside of the modelled contact region, all activated Tyro3 in cells with no hindrance. In conclusion, we show that the LG1 domain of ProS1 is necessary for activation of the Tyro3 receptor, involving protein-protein interaction interfaces that are homologous to those of the Gas6-Axl interaction.     

Structure-based virtual ligand screening with LigandFit: pose prediction and enrichment of compound collections

Virtual ligand screening methods based on the structure of the receptor are extensively used to facilitate the discovery of lead compounds. In the present study, we investigated the LigandFit package on four different proteins (coagulation factor VIIa, estrogen receptor, thymidine kinase, and neuraminidase), a relatively large compound collection of 65,560 unique "drug-like" molecules and four focused libraries (1950 molecules each). We performed virtual screening experiments with the large database and evaluated six scoring functions available in the package (DockScore, LigScore1, LigScore2, PLP1, PLP2, and PMF). The results showed that LigandFit is an efficient program, especially when used with LigScore1. Similar computations were carried out using focused libraries. In this situation the LigScore1 scoring function outperformed the other ones on three out of the four proteins tested. Even for the difficult neuraminidase case, the LigandFit/LigScore1 combination was still reasonably successful. Assessment of docking accuracy was also performed and again, we found that LigandFit (with DockScore and the CFF parameters) was performing well. On the basis of these results and observed increased enrichments after LigandFit/Ligscore1 screening on focused libraries, we suggest that using this program as a final step of a hierarchical protocol can be very beneficial to assist lead finding.     

In comparison with nitrate nutrition,ammonium nutrition increases growth of the frostbite1 Arabidopsis mutant

Ammonium nutrition inhibits the growth of many plant species, including Arabidopsis thaliana. The toxicity of ammonium is associated with changes in the cellular redox state. The cellular oxidant/antioxidant balance is controlled by mitochondrial electron transport chain. In this study, we analysed the redox metabolism of frostbite1 (fro1) plants, which lack mitochondrial respiratory chain complex I. Surprisingly, the growth of fro1 plants increased under ammonium nutrition. Ammonium nutrition increased the reduction level of pyridine nucleotides in the leaves of wild‐type plants, but not in the leaves of fro1 mutant plants. The observed higher activities of type II NADH dehydrogenases and cytochrome c oxidase in the mitochondrial electron transport chain may improve the energy metabolism of fro1 plants grown on ammonium. Additionally, the observed changes in reactive oxygen species (ROS) metabolism in the apoplast may be important for determining the growth of fro1 under ammonium nutrition. Moreover, bioinformatic analyses showed that the gene expression changes in fro1 plants significantly overlap with the changes previously observed in plants with a modified apoplastic pH. Overall, the results suggest a pronounced connection between the mitochondrial redox system and the apoplastic pH and ROS levels, which may modify cell wall plasticity and influence growth.     

Partial glycosylation of Asn2181 in human factor V as a cause of molecular and functional heterogeneity. Modulation of glycosylation efficiency by mutagenesis of the consensus sequence for N-linked glycosylation.

Coagulation factor V (FV) circulates in two forms, FV1 and FV2, having slightly different molecular masses and phospholipid-binding properties. The aim was to determine whether this heterogeneity is due to the degree of glycosylation of Asn(2181). FVa1 and FVa2 were isolated and digested with endoglycosidase PNGase F. As judged by Western blotting, the FVa2 light chain contained two N-linked carbohydrates, whereas FVa1 contained three. Wild-type FV and three mutants, Asn(2181)Gln, Ser(2183)Thr, and Ser(2183)Ala, were expressed in COS1 cells, activated by thrombin, and analyzed by Western blotting. Wild-type FVa contained the 71 kDa-74 kDa doublet, whereas the Asn(2181)Gln and Ser(2183)Ala mutants contained only the 71 kDa light chain. In contrast, the Ser(2183)Thr mutant gave a 74 kDa light chain. This demonstrated that the third position in the Asn-X-Ser/Thr consensus affects glycosylation efficiency, Thr being associated with a higher degree of glycosylation than Ser. The Ser(2183)Thr mutant FVa was functionally indistinguishable from plasma-purified FVa1, whereas Asn(2181)Gln and Ser(2183)Ala mutants behaved like FVa2. Thus, the carbohydrate at Asn(2181) impaired the interaction between FVa and the phospholipid membrane, an interpretation consistent with a structural analysis of a three-dimensional model of the C2 domain and the position of a proposed phospholipid-binding site. In conclusion, we show that the FV1-FV2 heterogeneity is caused by differential glycosylation of Asn(2181) related to the presence of a Ser rather than a Thr at the third position in the consensus sequence of glycosylation.     

基于Methylovorus sp.J1-1基因组尺度代谢网络优化吡咯喹啉醌合成

通过构建Methylovorus sp.J1-1基因组尺度代谢网络模型(genome scale of metabolic network model,GSMM),发掘能够提升吡咯喹啉醌(pyrroloquinoline quinone,PQQ)产量的发酵策略和相关靶基因.基于其注释的基因组和生化信息,构建J1-1的G...