BRAF蛋白V600E突变的影响

最新推荐文章于 2021-02-20 10:59:19 发布
最新推荐文章于 2021-02-20 10:59:19 发布
阅读量8.1k 收藏 13
点赞数
分类专栏： 分析流程 文章标签： 蛋白质 突变
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_29300341/article/details/80586459
版权

BRAF蛋白V600E突变的影响

目的：使用Cosmic、Clinvar、dbsnp等多种数据库，注释错义突变对蛋白质功能、稳定性和相互作用的影响。

步骤：

  1. 查找突变数据：选择癌症相关基因BARF，从COSMIC数据库下载BARF错义突变数据。
  2. 查找结构数据：进入PDB，搜索BRAF，下载BRAF蛋白质单体的结构pdb文件、以及BRAF与其他蛋白结合的复合物pdb文件。
  3. 突变数据分析：使用R语言，统计突变相关位点，筛选出高频率突变位点。
  4. 使用cBioPortal工具对突变进行可视化。
  5. 使用dbsnp、ClinVar、pdbsum等数据库对BRAF相关突变信息调研。
  6. 使用IBIS数据库对BRAF结合位点进行调研。
  7. NCBI Biosystem、uniprot、pubmed等数据库对BRAF相关通路信息调研。

一、COSMIC数据下载

进入COSMIC，搜索BRAF基因，查找突变信息，可以看出来，Mutation中99.11%为Missense substitution，选择Missense substitution，下载全部错义突变数据，共计50046条记录。
在这里插入图片描述
图表 1 突变类型分布
接下来，使用R语言，统计错义突变的分布，由饼状图可以发现，在93.8%的样本中，第600位缬氨酸(V)突变为谷氨酸(E)，这暗示着V600E在全部的Missense mutation中有着很强的统计学意义。并且可以发现： D594，F595，G596，L597，A598，T599，V600，K601这一小段的氨基酸序列突变频率远高出其他区域，暗示着这段氨基酸序列对整体结构有着极为重要的意义。
在这里插入图片描述
图表 2错义突变类型分布

二、pdb数据库调研

进入pdb数据库，搜索BRAF蛋白，选取一个“BRAF蛋白+其他蛋白复合物”和一个BRAF蛋白质单体。蛋白质复合物需要完整蛋白，不能是肽段，而且均需要来源于人类。经过一系列仔细审查，最终选择编号“4MNE的BRAF+MEK1蛋白质复合物”以及“4WO5的BRAF蛋白质单体”。

三、cBioPortal可视化

由于记录总数达到50046条，因此选择至少在两个样本中出现的突变，共计78种突变类型，转化成bed文件格式，放入cBioPortal进行可视化。可以发现在“酪氨酸激酶功能域”上突变有明显地富集，特别是第600位氨基酸附近，突变种类十分丰富。
在这里插入图片描述
图表 3cBioPortal突变可视化

接下来，显示蛋白质的3D structure，可以发现这一小段氨基酸位于蛋白质的环状区域。
在这里插入图片描述
图表 4 cBioPortal chain B结构可视化
在这里插入图片描述
图表 5 cBioPortal chain C结构可视化
在这里插入图片描述
图表 6 cBioPortal chain F结构可视化
在这里插入图片描述
图表 7cBioPortal chain G结构可视化

四、NCBI数据库调研

将COSMIC下载的数据进行格式转换，处理成bed文件格式，提交至NCBI variation viewer，在ClinVar一栏有着许多紫色与绿色盒子，紫色盒子表明致病变体，绿色盒子表明良性变体，发现在chr7:140753335附近处，ClinVar的致病突变明显富集，加号紫色盒子里面有着33个突变样本，而且COSMIC里的突变数据也富集在此处。
在这里插入图片描述
图表 8 NCBI variation viewer
选择V600E疾病样本点进ClinVar数据库，发现该突变可以引发很多种疾病，提交次数多达25次，会导致以下疾病：Gastrointestinal stromal tumor、Lung cancer、Malignant melanoma of skin、Multiple myeloma、Squamous cell carcinoma of the head and neck、Papillary thyroid carcinoma、Adenocarcinoma of lung、Non-small cell lung cancer、Carcinoma of colon、Astrocytoma, low-grade, somatic、Germ cell tumor, nonseminomatous、NONSEMINOMATOUS GERM CELL TUMORS, SOMATIC、VARIANT OF UNKNOWN SIGNIFICANCE、Neoplasm of brain、Malignant melanoma、Glioblastoma、Brainstem glioma、Neoplasm、Cardio-facio-cutaneous syndrome、Ovarian Neoplasms、Colorectal Neoplasms、Papillary renal cell carcinoma, sporadic、Colonic Neoplasms。
在这里插入图片描述
图表 9 ClinVar数据库调研
接着，选择同一突变，点进dbSNP数据库，发现V600的各种missense mutation同样也在dbSNP数据库存储。
在这里插入图片描述
图表 10dbSNP数据库调研
然后，链接进入NCBI Structure数据库，使用Cn3D对突变进行可视化，发现V600处于一个环状结构的转角处。
在这里插入图片描述
图表 11 Cn3D可视化

五、其他数据库调研

Uniprot

进入Uniprot搜索BRAF蛋白，可以得到关于BRAF蛋白更多更全面的信息。首先，该蛋白与RGS14和RAF1(Ras-Raf-MAPK通路中重要蛋白)共同存在于“细胞质”和“细胞膜”上，在“细胞核”内也有存在。
在这里插入图片描述
图表 12uniprot BRAF定位
其次，在OMIM中有报道，BRAF蛋白质的损坏会导致包括Colorectal cancer 、Lung cancer、Familial non-Hodgkin lymphoma、Cardiofaciocutaneous syndrome 1、Noonan syndrome 7、LEOPARD syndrome 3在内的多种疾病。
并且，BRAF的V600E是最常见的体细胞突变，能有效提高BRAF的激酶活性，会导致包括sarcoma、Colorectal cancer、metastatic melanoma、ovarian serous carcinoma、pilocytic astrocytoma在内的多种疾病。

Pdbsum

进入Pdbsum数据库搜索4MNE，选择BRAF蛋白质，发现该蛋白有着“Raf-like Ras-binding domain”、“C1_1”和“Protein tyrosine kinase”三个保守结构域，而V600位于“酪氨酸激酶”结构域，该蛋白催化着ATP + a protein = ADP + a phosphoprotein的反应
在这里插入图片描述
图表 13Pdbsum催化反应

六、多序列比对

下载uniprot中源自不同物种BRAF蛋白质序列，分别来自“绿脓杆菌”、“老鼠”、“人类”、 “鸡”和“鹌鹑”。使用cluster Omega进行多序列比对，并且绘制进化树，从进化树可以看出来多序列比多没有问题，“绿脓杆菌”作为外类，和其他物种距离较远。
在这里插入图片描述
图表 14进化树
而且，仔细观察多序列比对结果，发现G596、L597、T599、V600有着很高的保守性。
在这里插入图片描述
图表 15多序列比对结果

七、IBIS蛋白质相互作用

打开IBIS网页，并搜索4MNE结构链B (4MNE的B, C, F, G链为BRAF蛋白)，直接可视化所有与4MNE结构中BRAF蛋白有着相互作用的结构，可以观察到MEK1，Cyclin这些信号通路中关键蛋白(后续pubmed文献中也证实这些蛋白与BRAF有着结合)。
在这里插入图片描述
图表 16 IBIS Protein-protein interactions
接着，转到“protein-protein”选项卡，可以看到4MNE自身在“STKc_Raf”和“PKc_MEK1”域的标签栏目中被检索到，显示BRAF与“其他蛋白的STKc_Raf和PKc_MEK1结构域”存在相互作用，在BRAF序列上的三角形表示结合位点，红色表示高度保守的结合位点。

Interaction PartnerRanking ScoreNumber of Cluster MembersAverage %Identity to QueryNumber of Binding Site ResiduesTaxonomic Diversity
STKc_Raf2.8468840Eukaryota
PKc_MEK1singleton210029Homo sapiens

图表 17 List of protein interaction partners

八、Biosystems信号通路查询

进入NCBI biosystems，搜索“BRAF V600E”，可以查到三条记录“Signaling by high-kinase activity BRAF mutants”、“Signaling by moderate kinase activity BRAF mutants”和“Paradoxical activation of RAF signaling by kinase inactive BRAF”，阐述了BRAF在8%人类癌症中突变，而且大多数BRAF突变落在“激酶区”或者其相邻的“激活环区”，通过“模拟活化环的磷酸化作用”或者“促进酶的活性构象”来增加激酶活性。而且，大约90%的BRAF突变体都是单次错义突变V600E。

相关链接：
https://www.ncbi.nlm.nih.gov/biosystems/1457783
https://www.ncbi.nlm.nih.gov/biosystems/1457784
https://www.ncbi.nlm.nih.gov/biosystems/1457785

九、Pubmed文献查询

BRAF的研究：

其实，早在2002年，Nature上就有关于BRAF V600E突变的研究，第1799位核苷酸T突变为A(T1799A2)，导致其编码蛋白第600位谷氨酸取代了缬氨酸，通过消除激酶域与甘氨酸P环间的相互作用，模仿T599和S602两个位点同时被磷酸化的效果，突变后的BRAF蛋白的激酶活性会较野生型提高500倍。
后来，2010年，Cell上又发表封面文章，研究发现致死性突变的BRAF V600E需要与致癌基因RAS共同作用，一起促进肿瘤的生长发育。
事实上，BRAF V600E突变位点已经应用于临床检测。2011 年,首个BRAFV600E 靶向抑制剂威罗菲尼被FDA 批准上市，用于治疗BRAFV600E 突变的晚期黑色素瘤患者，有效延长了患者无进展生存期及总生存期，取得了突破性的治疗效果，也是典型的基于基因诊断选择用药的靶向治疗药物。

相关文献：

Hubbard SR,Oncogenic mutations in B-Raf: some losses yield gains. Cell, 116(6):764-6.
Wan PT,Mechanism of activation of the RAF-ERK signaling pathway by oncogenic mutations of B-RAF. Cell, 2004.116(6):855-67. PMID: 15035987
Davies H,Mutations of the BRAF gene in human cancer. Nature,2002 Jun 27;417(6892):949-54. PMID: 12068308

wangyunpeng_bio
关注
专栏目录
BRAF蛋白F595S G615R突变的影响
生物医学信息学博客
06-05 2440
BRAF蛋白F595S G615R突变的影响 目的：使用在线Webserver，分析突变对蛋白质表型、稳定性、结合亲和力的影响。 步骤： 1. 查找突变数据：选择癌症相关基因BRAF，从COSMIC数据库下载BARF错义突变数据。 2. 查找结构数据：进入PDB，搜索BRAF，下载BRAF蛋白质单体的结构pdb文件、以及BRAF与其他蛋白结合的复合物pdb文件。 3. 突变数据分析...
生物信息学数据库大全
qiuxue198598的专栏
01-09 8774
 综合数据库★ INSD,国际核酸序列数据库（International Nucleotide Sequence Databank）。由日本的DDBJ、欧洲的EMBL和美国的GenBank三家各自建立和共同维护。★ EMBL库，欧洲分子生物学实验室的DNA和RNA 序列库。 http://www.ebi.ac.uk/embl.html ★ GenBank ，美国国家生物技术
转录因子详细介绍(motif)
wangprince2017
02-20 2万+
转录因子详细介绍(motif) TF: transcription factor转录因子TFBS: transcription factor binding site转录因子结合位点TFBS是序列内的location，TF特异结合在这里，这个site有这种特点 1 和一些参考相关的一个位置（开始，结束，strand），这些reference可以是染色体开始，geneTSS。也可以是一段sequence 2 A SITE可以是实验证实的（已知的），也可以是一些算法（预测的） 3 例子，下面这个图是酵母
SARS2-Variation-Viewer:GenBank中SARS-CoV-2序列的变异分析和可视化
05-28
GenBank中SARS-CoV-2序列的变异分析和可视化 NCBI数据集代码马拉松小组4 该项目的目标是识别和可视化SARS-CoV-2基因组中的变异（SNP） 背景 了解跨SARS-CoV-2基因组序列的变异对流行病学和分子生物学具有重要意义。 现有工具可让研究人员可视化SARS-CoV-2变异，但访问底层数据（包括基因组序列和有价值的元数据，例如在世界上某个特定基因组的位置受到限制）的访问受到限制。 SARS-CoV-2变异查看器使研究人员可以轻松地可视化映射到带注释的武汉-Hu-1参考基因组（NC_045512.2）的SARS-CoV-2 SNP，并带有可搜索的表格视图，显示氨基酸的变化，收集位置，宿主物种，并链接到iCn3D中的蛋白质结构视图。 工作流程 使用从GenBank下载所有SARS-CoV-2基因组 识别这些SARS-CoV-2基因组中的SNP 将SNP与SAR
NCBI上查看SNP位点在哪个基因座上（locus）
weixin_33725722的博客
04-12 3364
首先，进入NCBI的主页网站：https://www.ncbi.nlm.nih.gov/variation/view/ 进入后，在下图红色框框位置输入目的SNP，比如rs608139 输完后，出现如下结果，箭头指向的是该SNP当前所在基因座的位置，即p21，由于该SNP在2号染色体上，因为我们也写为2p21 转载于:https://www.cnblogs.com/chenwenyan...
论文研究 - 乳头状甲状腺癌BRAF V600E突变与临床病理特征的关系
05-13
在这里，我们通过执行qPCR和免疫组织化学（IHC）检查了BRAF V600E突变，并分析了总共188例PTC患者的BRAF V600E突变状态与多种临床病理特征之间的关系。 在患者性别（χ2= 1.252，P = 0.263），种族（χ2= 0.756，P ...
BRAF突变对成人胶质瘤的研究有什么意义
最新发布
05-29
BRAF基因突变在许多类型的肿瘤中都是重要的致病因素，包括成人胶质瘤。近年来，随着对BRAF突变的深入研究，科学家们发现了针对BRAF突变的药物，这些药物可以抑制肿瘤细胞的生长和扩散。这些药物可以用于治疗那些携带...
EGFR突变NSCLC治疗新进展.pptx
09-14
非小细胞肺癌（NSCLC）是肺癌中最常见的一种类型，其中EGFR（表皮生长因子受体）突变是重要的驱动基因，影响着疾病的发展和治疗策略。这篇内容主要讨论了EGFR突变在NSCLC中的角色、EGFR-TKI（酪氨酸激酶抑制剂）的...
REGy，p53，NOTCHl和BRAF在甲状腺癌中的表达及其临床意义 (2012年)
05-09
分析REGγ与p53，NOTCHl，BRAF之间的关系及其与临床病理因素的相关性．REGr在乳头状甲状腺癌组织中表达明显增高，且其阳性率为87．5％(35／40)．REGγ表达与患者年龄、淋巴结转移情况有关(P=0．015 1<0 uff0e05)，...
在python平台上利用pymol来查找PDB文件中蛋白质的相互作用位点
weixin_45832250的博客
08-23 1万+
在python平台上利用pymol来查找PDB文件中蛋白质的相互作用位点 最近要做一些关于蛋白质结合位点的实验，查阅了很多篇文献，大多数都是关于如何定义结合位点的，却没有说明该如何去找结合位点。大多数都是说从PDB库中下载蛋白质结构数据，但是却没有说下载过数据之后该怎么做。偶然间从论坛上看到利用pymol软件可以解析PDB文件中的原子坐标，也可以查找选中片段范围内的氨基酸，比如3.5A内的氨基酸。经过半个多月的摸索，发现可以直接在python上利用pymol包来找结合位点。 下图是利用pymol来找结合位点
分子对接（第一步：相关蛋白的信息查询）
qq_38016547的博客
03-08 1976
https://www.uniprot.org/uniprot/P05804 此界面可以直接转入PDB官网、Protein Data Bank in Europe、https://pdbj.org/、https://www.ebi.ac.uk/pdbsum（可能用于同源模建，包含具体的蛋白结构区域） 人之相识，贵在相知，人之相知，贵在知心。----------孟子 ...
检测蛋白质构象是否合理————Ramachandran plot（拉氏图）
MCANDML的博客
06-12 4万+
Ramachandran plot（拉氏图）是由G. N. Ramachandran等人[1]于1963年开发的，用来描述蛋白质结构中氨基酸残基二面角ψ和φ是否在合理区域的一种可视化方法。同时也可以反映出该蛋白质的构象是否合理。 在介绍Ramachandran plot前有必要简单说明一下二面角。图1 蛋白质链上肽平面之间形成的二面角以及 两个残基脱水缩合后的 ψ，φ和W（图片来自...
正常健康人组织也有癌基因和抑癌基因突变
weixin_44054863的博客
07-25 1398
大型癌症基因组计划如TCGA研究分析了各种类型癌症的基因组变异。很少的研究人员会关注正常组织中的基因突变。在今年6月份的Science杂志上发表了一项研究，科学家开发了一种新的方法，能够从转录组测序（RNA-seq）数据中更加精准的研究突变。这种方法被称为RNA-MuTect，这种方法能够用利用组织来源的RNA及配对的正常DNA来确定组织中的体细胞突变。 采用这种方法，他们分析了GTEx项目中正...
蛋白质结构预测(swiss-model,modeller )
热门推荐
生物医学信息学博客
04-07 8万+
用swiss-model (3 models)和modeller (6 models)分别预测给定序列的结构，并用PROCHECK，Molprobity，Errat, Verify_3D分别评价所得预测模型，详细注释所得结果 依据打分情况，选出最佳的模型，说明理由。 建模注意事项： 了解本序列的结构组成（有几个结构域，分别是什么结构域……） 有针对性选择合适的模板进行建模 最后选出的最优模型用VMD显示三维结构，并注明每个结构域的位置
蛋白质结构模型和功能预测：Swiss-model工具的使用
weixin_33853794的博客
12-01 8647
Swiss-model也是一款预测蛋白质结构模型的工具。网页地址：https://swissmodel.expasy.org/ 首先，进行常规的注册后，点击start modelling 以搜索BRCA1蛋白质为例： 进入NCBI的protein网页：https://www.ncbi.nlm.nih.gov/protein/，输入要查询的BRCA1 点击搜索后，选择最后的氨基酸序列： ...
Uniprot数据库
weixin_30847939的博客
04-13 7809
Uniprot数据库是Universal Protein的英文缩写，是信息最丰富、资源最广的蛋白质数据库。 UniprotKB由两部分组成： UniProtKB/Swiss-Prot 高质量的、手工注释的、非冗余的数据集，这些数据都是由质量保证的。 UniProtKB/TrEMBL 该数据集包含高质量的计算分析结果，需要我们手工注释。 当我们搜索蛋白质时，就会如下...
COSMIC数据库简介
庐州月光的博客
06-03 4526
欢迎关注"生信修炼手册"！COSMIC是癌症相关体细胞突变位点的最大的数据库之一，网址如下：https://cancer.sanger.ac.uk/cosmic/最新版本...
光学镊子测量RAS蛋白与病理BRAF突变体间的结合力
本文研究了通过光学镊子测量RAS蛋白与病理BRAF突变体之间的结合力，这两种蛋白在人类肿瘤和黑色素瘤中极为常见，占至少三分之一的病例。RAS和BRAF是细胞信号传导中的关键分子，它们的激活突变通常与癌症的发生密切...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值