您的位置 首页 工具使用

Rosetta Cartesian_ddG: 更快更准确的蛋白质稳定性预测方法

前言

计算单点突变自由能是分子稳定性改造的”老套路”。FoldX以及ddg_monomer application已经久经实战考验(具体用法参考前人总结:http://bioengx.hk1.91site.net/rosetta-ddg_monomer/)。今天我们就来说说Rosetta新一代的预测方法: Cartesian_ddG。

1 Cartesian_ddG方法概述

参考: Hahnbeom Park, Philip Bradley, Per Greisen Jr., Yuan Liu, Vikram Khipple Mulligan, David E Kim, David Baker, and Frank DiMaio (2016) “Simultaneous optimization of biomolecular energy function on features from small molecules and macromolecules”, JCTC.

Cartesian_ddG是一种新的采样方法,与ddg_monomer不同,对于骨架柔性的计算并不需要大量重复的强限制约束,而是采取了卡迪尔空间的优化来允许小幅度的骨架运动。

在计算∆∆G时,Cartesian_ddG方法使用两步Relax的方法:

-预先使用FastRelax进行优化,确定野生型以及点突变氨基酸最佳的侧链排布方式
-再使用卡迪尔空间版本的FastRelax对局部区域(点突变位置6埃范围内)进行彻底放松。
-评估野生型以及突变型蛋白的能量差,并使用能量函数特异性的校正因子对∆∆G值进一步矫正,以对应实验测定值单位kcal/mole。

效果:更快、更准。
Cartesian_ddG(使用opt-nov15函数)比ddg_monomer最佳Protocol(row16)要快上10倍有余。并且Cartesian_ddG预测∆∆G与实验∆∆G的皮尔森相关性从0.703升至0.743。

2 Cartesian_ddG的使用方法

Cartesian_ddG的使用方法十分简单。

-预优化野生型蛋白结构
-准备mutfile文件
-运行Cartesian_ddG
-分析∆∆G值

2.1 预优化野生型蛋白结构

创建relax_flag文件,文件内容如下:

-use_input_sc
-constrain_relax_to_start_coords
-ignore_unrecognized_res
-nstruct 20
-relax:coord_constrain_sidechains
-relax:cartesian-score:weights ref2015_cart # 务必设为ref2015_cart
-relax:min_type lbfgs_armijo_nonmonotone # 能量最小化的算法
-relax:script cart2.script
创建cart2.script文件,文件内容如下:(无需额外更改参数)

switch:cartesian
repeat 2
ramp_repack_min 0.02 0.01 1.0 50
ramp_repack_min 0.250 0.01 0.5 50
ramp_repack_min 0.550 0.01 0.0 100
ramp_repack_min 1 0.00001 0.0 200
accept_to_best
endrepeat
进行relax优化,打开终端并输入:

relax.mpi.macosclangrelease -s $pdb @relax_flag
打开relax的打分文件,查看最低能量的结构pdb编号,并以此作为出发的输入结构。

2.2 准备mutfile文件
单点突变的mutfile文件的格式如下:

注意!此处的89为pose序号,而非PDB中的残基号,设置错误一定会报错。
total 1 # mutfile中总突变数量
1 # 本轮突变氨基酸个数,单点突变设置为1.
T 89 A # 将89位的T突变为A
将文件保存为T89A.mutfile

2.3 运行Cartesian_ddG
Cartesian_ddG有两种运行模式: 蛋白稳定性预测模式(Protein stability mode)和 相互作用界面模式(Interface mode)。

1)蛋白稳定性预测模式(Protein stability mode)
创建cartddg_flag文件, 内容如下:

-ddg:iterations 5 # 默认为3,可以根据自身计算资源调整。
-ddg::cartesian
-ddg::dump_pdbs False # 是否输出突变后的蛋白PDB文件,推荐False,否则文件夹会很乱
-bbnbrs 1 # 骨架自由度,额外考虑线性邻居氨基酸的数目,1代表 邻近的3个氨基酸骨架自由度被考虑在计算中。
-fa_max_dis 9.0 # 控制范德华和溶剂化能量计算的范围,默认为6埃以内的氨基酸被考虑。
-score:weights ref2015_cart # 务必设为ref2015_cart
运行ddg计算,打开终端并输入,这么就会开始计算T89A的∆∆G。

cartesian_ddg.linuxgccrelease -s [inputpdb] @cartddg_flag -ddg:mut_file T89A.mutfile

2 )相互作用界面模式(Interface mode)
该模式开发者并未做详细测试,斟酌使用,此处笔者未做测试。如有需要请参考: https://www.rosettacommons.org/docs/latest/cartesian-ddG

3 结果分析

运行完毕后,文件夹中会产生T89A.ddg这个文件。打开后,文件内容如下:

COMPLEX: Round1: WT: -333.782 fa_atr: -840.140 fa_rep: ………….
COMPLEX: Round2: WT: -333.782 fa_atr: -840.140 fa_rep: ………….
COMPLEX: Round3: WT: -333.782 fa_atr: -840.140 fa_rep: ………….
COMPLEX: Round4: WT: -333.782 fa_atr: -840.140 fa_rep: ………….
COMPLEX: Round5: WT: -333.782 fa_atr: -840.140 fa_rep: ………….
COMPLEX: Round1: MUT_89ALA: -330.352 fa_atr: -834.516 ………….
COMPLEX: Round2: MUT_89ALA: -330.352 fa_atr: -834.516 ………….
COMPLEX: Round3: MUT_89ALA: -330.352 fa_atr: -834.516 ………….
COMPLEX: Round4: MUT_89ALA: -330.352 fa_atr: -834.516 ………….
COMPLEX: Round5: MUT_89ALA: -330.352 fa_atr: -834.516 ………….
此处Round代表第几次迭代,根据-ddg:iterations设置定义。

WT:代表本轮能量计算,野生型的能量为-333.782。

MUT_89ALA:代表本轮能量计算,突变型的能量为-330.352。

我们需要将WT和MUT_89ALA平均后,用MUT能量平均值减去WT能量平均值就可以得到T89A的∆∆G值,为3.43。

稳定性变化评判标准:(务必根据自身需要调整)

∆∆G < -1.0 代表点突变为蛋白带来稳定效应; ∆∆G > 1.0 代表点突变为蛋白带来不稳定效应;

备注: 如果需要批量分析和点突变扫描,请使用脚本生成多个mutfile,然后按照以上流程进行计算。



发表回复

您的电子邮箱地址不会被公开。 必填项已用*标注

评论列表(266)

  1. This design is steller! You obviously know how to keep a reader amused.
    Between your wit and your videos, I was almost moved
    to start my own blog (well, almost…HaHa!) Great job.

    I really enjoyed what you had to say, and more than that, how you presented it.
    Too cool!

    My site JUDI ONLINE

  2. “Was it because of that?” His dad asks as he points a finger towards him. And his midsection.

  3. I am curious to find out what blog platform you are working
    with? I’m experiencing some minor security issues with my latest blog
    and I’d like to find something more safe. Do you have any recommendations?

    Review my website :: slot gacor bonus deposit

  4. This is really interesting, You’re a very skilled blogger.
    I have joined your feed and look forward to seeking more of your excellent post.
    Also, I have shared your web site in my social networks!

    Here is my site JUDI ONLINE

  5. “I look forward to the day when you are a fully grown man, my son.” His father says. “I am anxious to see the fruit of my loins reach his fully-adorned manhood. Adulthood.” “That’s it! That’s it! That’s it!” His father bellows.

  6. “Once you plant that cock of yours in some squirmy hole and empty those warm balls of yours into a moist wet hole. Then you are on your way to becoming a man, and only then will you, come-of-age. But it is only a step onto the winding pathway towards manhood. It is my job to teach you what it means to be a man.” His father stands as he finishes his sentence. Loosening his belt and pulling his shirt out from the tucked confines of his pants. “I ain’t a-talking ‘bout how big ya cock is, boy. I am talking ‘bout ya, weight. Your height.” His dad says. “I kinda figured on how big you are down there. That’s obvious. It makes me proud of you. My son. Of what I created, along with ya mother.”

  7. Thank you for sharing your info. I really appreciate your efforts and I will be waiting for your further write
    ups thanks once again.

    Look at my web blog … JUDI ONLINE

  8. Hi there, just became aware of your blog through Google,
    and found that it’s really informative. I am gonna watch out for brussels.
    I will be grateful if you continue this in future.
    Lots of people will be benefited from your writing.
    Cheers!

    Take a look at my web site … slot gacor china

  9. Well voiced really. !
    casino aladdin online [url=https://usagamblinghub.com/#]us online casino[/url] billionaire casino online

  10. Very good write ups, Cheers!
    free online casino slots with bonus rounds [url=https://usacasinomaster.com/#]free casino online[/url] bandar taruhan casino joker123 online

  11. His father’s words are what he hears when he erupts. His cum streaming like liquid threads from the pee-hole of his rigid cock. gay sex porn “Yes, sir.” He says as he straightens up, standing erect as his member pulses to life between his legs, and in his father’s right hand, which are firmly locked on his balls. “Yes, sir, I do.”

  12. Wonderful write ups. Thank you!
    aztec riches casino online [url=https://luckyusaplay.com/#]best online usa casino[/url] play online casino free win real money

  13. Truly a good deal of great tips.
    william hill casino online [url=https://luckyusaplay.com/#]casino bonus online[/url] encore online casino

  14. I am really inspired with your writing talents as well as with the structure
    for your weblog. Is this a paid subject matter or
    did you customize it yourself? Either way stay up the nice quality writing, it is
    uncommon to peer a nice weblog like this one nowadays..

    Also visit my blog – JUDI ONLINE

  15. Hey just wanted to give you a quick heads up. The
    words in your article seem to be running off the screen in Chrome.

    I’m not sure if this is a format issue or something to do with browser compatibility but I figured I’d
    post to let you know. The design look great though!

    Hope you get the problem fixed soon. Thanks

    Stop by my web-site :: JUDI ONLINE

  16. Hello friends, its fantastic post concerning teachingand entirely explained, keep it up all the time.

    Feel free to visit my web site; JUDI ONLINE

  17. You have made your point!
    online casino app philippines [url=https://uscasinoguides.com/#]online casino gaming[/url] how to report online casino

  18. Thank you! A lot of posts.
    top dragonfish online casinos [url=https://usagamblinghub.com/#]free casino online[/url] 1 win casino online

  19. Great info, Kudos!
    list of the best online casinos [url=https://usacasinomaster.com/#]best online casinos that payout[/url] como ganar en casino online

  20. Nicely spoken without a doubt! !
    online casino schleswig holstein [url=https://luckyusaplay.com/#]play casino online[/url] peermont online casino

  21. Thanks a lot! Ample data.
    win big online casino [url=https://uscasinoguides.com/#]live dealer online casino[/url] online kostenlos casino

  22. Amazing write ups. Thanks.
    poker star casino online [url=https://usaplayerscasino.com/#]casino gambling online[/url] casino bet online real money

  23. Superb knowledge. Appreciate it!
    online casino legit [url=https://usagamblinghub.com/#]play online casino[/url] el royale casino online reviews

  24. Many thanks! A lot of posts!
    best usa online casinos with no deposit bonus [url=https://usacasinomaster.com/#]casino online no deposit bonus[/url] horseshoe casino application online

  25. Regards. I appreciate this.
    online casino paddy power [url=https://luckyusaplay.com/#]best online casino bonus[/url] tips casino online

  26. Hello There. I found your blog using msn. This is
    a very well written article. I will make sure to
    bookmark it and return to read more of your useful info. Thanks for the post.
    I’ll certainly comeback.

    My blog post; betxhub 247

  27. Hi there, You have done a fantastic job. I’ll certainly digg it and personally suggest to my friends.
    I’m sure they’ll be benefited from this site.

    Feel free to surf to my web page :: JUDI ONLINE

  28. You said this fantastically!
    online casino lastschrift [url=https://uscasinoguides.com/#]best online usa casino[/url] free casino online slot games

  29. What’s up all, here every one is sharing these kinds of experience, thus it’s nice to
    read this website, and I used to visit this webpage every day.

    Look at my web page … ro top serv

  30. Superb tips. Appreciate it.
    leovegas casino online [url=https://usagamblinghub.com/#]best usa online casino[/url] are online casinos legal in india

  31. Beneficial knowledge. Kudos!
    gratis geld online casino [url=https://usacasinomaster.com/#]casino games free online[/url] new mexico online casinos

  32. Thanks! I like it.
    10 online casino [url=https://luckyusaplay.com/#]no deposit online casino[/url] online casino pa real money

  33. You actually said it well.
    golden lady casino online [url=https://uscasinoguides.com/#]online casino games free[/url] white lion online casino

  34. Incredible a lot of helpful knowledge!
    usa friendly online casinos with free daily slots tournaments [url=https://usagamblinghub.com/#]casino games free online[/url] best online casino ios

  35. Incredible a good deal of helpful knowledge!
    vt online casino bonus [url=https://usacasinomaster.com/#]online casinos with no deposit bonus[/url] 21 grand casino online

  36. What a information of un-ambiguity and preserveness of precious familiarity concerning unpredicted feelings.

    Feel free to visit my web blog bollyflix

  37. Superb posts. Thanks!
    agen betting maxbet casino online [url=https://uscasinoguides.com/#]new online casino usa ok[/url] jack’s casino online review

  38. Perfectly spoken genuinely! .
    casino online che accettano american express [url=https://usaplayerscasino.com/#]deposit bonus casino online[/url] bandar judi osg777 casino online

  39. Incredible a good deal of great data!
    do online casinos cheat blackjack [url=https://usagamblinghub.com/#]best online casinos that payout[/url] alle neuen online casinos

  40. You explained it really well!
    malaysia online casino website [url=https://usacasinomaster.com/#]casino online casino[/url] online casino rollover

  41. Kudos, Terrific information!
    online casino real money india [url=https://luckyusaplay.com/#]live online casino[/url] rtg us online casinos no deposit

  42. You really make it seem really easy along with your presentation however I to find this matter to be
    really something that I believe I would by no
    means understand. It kind of feels too complex and extremely wide for me.
    I am having a look forward to your subsequent publish,
    I’ll attempt to get the hold of it!

    Here is my web blog :: Free Link Shortener

  43. Whats up this is somewhat of off topic but I was wanting
    to know if blogs use WYSIWYG editors or if you have to manually code with HTML.
    I’m starting a blog soon but have no coding expertise so I wanted to get guidance from someone with experience.
    Any help would be greatly appreciated!

    Also visit my webpage :: URL shortener algorithm

  44. You said it perfectly..
    ignition casino online poker [url=https://usagamblinghub.com/#]best online casinos that payout[/url] no deposit bonus codes online casino

  45. If some one needs expert view on the topic of running a blog then i
    suggest him/her to visit this web site, Keep up the pleasant job.

    Feel free to visit my site: אידוי קנאביס

  46. Regards! Loads of facts!
    river rock casino online [url=https://usacasinomaster.com/#]best online casinos[/url] best international online casino

  47. You actually explained this well.
    new online casino usa real money [url=https://luckyusaplay.com/#]casino games free online[/url] analisis casinos online colombia

  48. Regards, I like it!
    2022 online casinos real money [url=https://uscasinoguides.com/#]legit online casino[/url] stake online casino location

  49. You stated it well.
    casino cyprus online [url=https://luckyusaplay.com/#]no deposit online casino[/url] online casinos portugal

  50. This is nicely said! .
    casino online argentina dinero real [url=https://uscasinoguides.com/#]casino online usa[/url] real money online casino no deposit bonus usa

  51. Thank you. Loads of stuff.
    best online casinos india [url=https://usaplayerscasino.com/#]casino gambling online[/url] river sweepstakes online casino

  52. Valuable facts. Thanks a lot.
    caesars online casino review [url=https://usagamblinghub.com/#]online casino with no deposit bonus[/url] online casino mastercard withdrawal

  53. Many thanks. Lots of stuff!
    south african online casinos with no deposit coupons [url=https://usacasinomaster.com/#]best online casinos[/url] online casino foxwoods

  54. Many thanks! An abundance of postings.
    how to win online casino jackpot [url=https://luckyusaplay.com/#]online casino top[/url] online casinos mit sofortГјberweisung

  55. It’s truly a great and useful piece of information. I’m happy that you shared this
    helpful information with us. Please keep us
    up to date like this. Thank you for sharing.

    My page … get visibility with seo and markerting (Donte)

  56. It’s a shame you don’t have a donate button! I’d definitely donate to this outstanding blog! I guess for now i’ll settle for book-marking and adding your RSS feed to my Google account. I look forward to brand new updates and will share this website with my Facebook group. Chat soon!

    https://youtu.be/bqeueF7oPzE

  57. Attractive section of content. I just stumbled upon your site and in accession capital to say that I acquire in fact loved account your weblog posts. Anyway I’ll be subscribing in your feeds and even I achievement you get admission to constantly rapidly.

    https://youtu.be/B16dCq5p6FI

  58. Greetings! This is my first visit to your blog! We are a collection of volunteers and starting a new initiative in a community in the same niche. Your blog provided us beneficial information to work on. You have done a marvellous job!

    https://rj-rajaslot88.com

  59. Virtually all of whatever you mention happens to be astonishingly accurate and it makes me wonder the reason why I hadn’t looked at this with this light previously. This particular piece truly did switch the light on for me as far as this specific subject goes. However at this time there is 1 factor I am not really too comfy with and while I attempt to reconcile that with the actual core theme of your position, let me see what all the rest of your visitors have to point out.Nicely done.

    https://youtu.be/2EAExRR5Du0

  60. I¦ll immediately grasp your rss as I can not in finding your e-mail subscription link or newsletter service. Do you’ve any? Kindly allow me understand so that I may subscribe. Thanks.

    https://youtu.be/FgevmHDERA0

  61. Somebody essentially help to make seriously posts I would state. This is the very first time I frequented your web page and thus far? I amazed with the research you made to create this particular publish amazing. Wonderful job!

    https://youtu.be/Iwo1I4DrVhM

  62. Hiya, I am really glad I have found this information. Nowadays bloggers publish just about gossips and internet and this is actually irritating. A good web site with interesting content, that is what I need. Thanks for keeping this website, I will be visiting it. Do you do newsletters? Can not find it.

    https://youtu.be/n3ZJVtprREU

  63. What is Lottery Defeater Software? Lottery Defeater Software is a plug-and-play Lottery Winning Software that is fully automated. Kenneth created the Lottery Defeater software. Every time someone plays the lottery, it increases their odds of winning by around 98.

    https://youtu.be/Wb3pO9Kra0M

  64. I have been surfing online greater than 3 hours lately, yet I by no means found any attention-grabbing article like yours. It¦s beautiful price sufficient for me. Personally, if all webmasters and bloggers made good content as you probably did, the web might be much more helpful than ever before.

    https://youtu.be/VEXFdh00ODk

  65. I was just seeking this information for some time. After six hours of continuous Googleing, finally I got it in your web site. I wonder what is the lack of Google strategy that do not rank this kind of informative web sites in top of the list. Generally the top websites are full of garbage.

    https://youtu.be/RZj5HVtZ4pQ

  66. This is very interesting, You are a very skilled blogger. I have joined your rss feed and look forward to seeking more of your fantastic post. Also, I have shared your site in my social networks!

    http://www.tlovertonet.com/

  67. Howdy! I know this is kind of off topic but I was
    wondering which blog platform are you using for this site?
    I’m getting fed up of WordPress because I’ve had issues with hackers and I’m looking at options for another platform.

    I would be great if you could point me in the direction of a
    good platform.

    Also visit my homepage eharmony special coupon code 2024

  68. Hi,

    I just visited bioengx.com and wondered if you’d ever thought about having an engaging video to explain what you do?

    I can show you some previous videos we’ve done if you want me to send some over. Let me know if you’re interested in seeing samples of our previous work.

    Regards,
    Georgina

    Unsubscribe: https://removeme.click/ev/unsubscribe.php?d=bioengx.com

  69. You actually make it appear really easy along with your presentation however I to find this topic to be actually
    one thing that I feel I’d by no means understand.
    It seems too complicated and extremely large for me. I’m looking forward for your subsequent post, I’ll try
    to get the hold of it!

    Feel free to surf to my web-site vpn coupon code 2024

  70. Hi there! Would you mind if I share your blog with my zynga group?
    There’s a lot of people that I think would really enjoy your content.
    Please let me know. Thanks

    Feel free to surf to my web site; vpn special code

  71. Hi there,

    We run a Youtube, Instagram and Twitter/X service, where we can increase your subscribers/followers.

    This is all done manually and so all subscribers and followers are real people.

    Our prices start from just $60/month

    If you are interested, please let us know, which of your social media profiles are you wishing to grow?

    If you are not interested, no problem, just ignore and delete this email.

    Kind Regards,
    Sarah

  72. 创建relax_flag文件,文件内容如下:

    -use_input_sc
    -constrain_relax_to_start_coords
    -ignore_unrecognized_res
    -nstruct 20
    -relax:coord_constrain_sidechains
    -relax:cartesian-score:weights ref2015_cart # 务必设为ref2015_cart
    -relax:min_type lbfgs_armijo_nonmonotone # 能量最小化的算法
    -relax:script cart2.script

    should be

    -use_input_sc
    -constrain_relax_to_start_coords
    -ignore_unrecognized_res
    -nstruct 20
    -relax:coord_constrain_sidechains
    -relax:cartesian
    -score:weights ref2015_cart # 务必设为ref2015_cart
    -relax:min_type lbfgs_armijo_nonmonotone # 能量最小化的算法
    -relax:script cart2.script

    1. Hi there,

      We run a YouTube growth service, which increases your number of subscribers both safely and practically.

      – We guarantee to gain you 700-1500 new subscribers per month.
      – People subscribe because they are interested in your videos/channel, increasing video likes, comments and interaction.
      – All actions are made manually by our team. We do not use any ‘bots’.

      The price is just $60 (USD) per month, and we can start immediately.

      If you’d like to see some of our previous work, let me know, and we can discuss it further.

      Kind Regards,
      Emily

联系我们

联系我们

(44)07934433023

在线咨询: QQ交谈

邮箱: info@bioengx.org

关注微信
微信扫一扫关注我们

微信扫一扫关注我们

关注微博
返回顶部