Amber中的信息传递——章节1.1-第二部分

2023-10-08 19:20

本文主要是介绍Amber中的信息传递——章节1.1-第二部分,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

Amber中的信息传递在实操中共分为预备程序、模拟程序和分析程序三个部分,具体相关文件如下:

1. 预备程序
**LEaP:**是在 Amber 中创建新系统或修改现有系统的主要程序。 它有命令行程序 tleap 和图形用户界面 xleap 两种形式。它结合了 Amber 早期版本中的 prep、link、edit 和parm 功能。
**pdb4amber:**通常帮助准备来自其他地方(如 rcsb.org)的 pdb 格式文件,使其与 LEaP 兼容。
prepareforleap:并不是一个程序,而是 cpptraj 中的一个操作,也可以帮助制作与 LEaP 兼容的 pdb 格式文件。它对碳水化合物特别有用。
parmed:提供了一种提取参数拓扑文件中定义的参数信息的简单方法。它还可
用于检查参数拓扑文件对复杂系统是否有效
(参见 checkValidity 命令),并可对该文件进行简单修改。
antechamber:是使用通用 Amber 力场 (general Amber force field,GAFF) 版本为小分子有机物(如药物、修饰氨基酸)开发力场的主要程序。这些力场可直接用于 LEaP,也可作为进一步开发参数的起点。
MCPB.py: 为金属蛋白和有机金属化合物 MM 模型的构建、原型设计和验证提供了一种方法。 它使用键合加静电模型来扩展现有的成对相加力场。它是 MCPB 在 Python 中的重新实现,具有更高效的工作流程,并自动集成了以前版本中的许多建模过程。
IPMach.py: 为离子非键模型(12-6 LJ 模型和 12-6-4 LJ 型模型)的参数化提供了便利。
mdgx: 允许通过量子数据拟合生成任何分子的键合力场参数。
packmol-memgen: 提供了一种生成膜系统的简单方法,无论是否含有蛋白质,只要用 Memembed 确定输入蛋白质的方向,并使用 Packmol 作为堆积引擎即可。 它可以处理复杂的脂质混合物以及多层系统。 输出结果通过 charmm- lipid2amber.py 与 Amber 兼容。

2. 模拟程序
sander(AmberTools 的一部分): 是基本的能量最小化和分子动力学程序。 该程序通过向下迭代移动原子的能量梯度来松弛结构,直到获得足够低的平均梯度。分子动力学部分通过整合牛顿运动方程生成系统的构型。 与最小化相比,分子动力学将采样更多的构型空间,并允许结构跨越较小的势能障碍。 可以在模拟过程中定时保存构型,以便日后进行分析,还可以利用热力学积分进行基本的自由能计算。 还可以使用 sander 模块进行更复杂的构象搜索和建模 MD 研究。 该模块允许在基本力场中添加各种约束条件,是专为核磁共振、X 射线或低温电子显微镜结构完善所涉及的计算类型而设计的。

pmemd(Amber 的一部分): 是针对速度和并行扩展进行了优化的 sander 版本;pmemd.cuda 变体可在 GPU 上运行。其名称代表 “粒子网格埃沃德分子动力学”,但该代码现在也可以进行广义玻恩模拟。输入和输出与 sander 相比只有一些变化。

gem.pmemd(AmberTools 的一部分): 是 pmemd 程序的(CPU 专用)变体,专为使用 AMOEBA[16] 和 GEM 等 "高级 "力场进行计算而设计。

3. 分析程序
mdout_analyzer.py: 是一个运行简单的 Python 脚本,可提供 sander 或 pmemd 输出文件中的信息摘要
cpptraj:主要的轨迹分析工具(用 C++ 编写),用于叠加、提取坐标、计算键值/角值/二面体值、原子位置波动、相关函数、氢键分析等。更多信息请参见第 35 章。
pytraj: 是 cpptraj 的 Python 封装器。 它通过与 Python 丰富的生态系统(如 numpy、scipy 和 ipython-notebook 等)相结合,为数据分析引入了更多灵活性
pbsa: 是一个分析程序,用于分析以溶剂为媒介的生物分子能量学。 pbsa.cuda 变体可在 GPU 上运行。 它可以利用分子动力学模拟和其他来源(pqr 格式)的输入坐标文件执行静电和非静电连续溶解计算。 它还支持在各种可视化程序中将溶剂介导的静电位可视化。更多信息请参见第 6 章。
MMPBSA.py: 是一个 python 脚本,利用连续介质模型产生的思想自动对分子动力学模拟的快照进行能量分析。(还有一个更早的 perl 脚本,名为 mm_pbsa.pl,也具有类似功能)。
FEW:(自由能工作流)使用 TI、MM/PBSA 类型或 LIE 计算方法自动计算蛋白质与配体结合的自由能。

这篇关于Amber中的信息传递——章节1.1-第二部分的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/167472

相关文章

闲置电脑也能活出第二春?鲁大师AiNAS让你动动手指就能轻松部署

对于大多数人而言,在这个“数据爆炸”的时代或多或少都遇到过存储告急的情况,这使得“存储焦虑”不再是个别现象,而将会是随着软件的不断臃肿而越来越普遍的情况。从不少手机厂商都开始将存储上限提升至1TB可以见得,我们似乎正处在互联网信息飞速增长的阶段,对于存储的需求也将会不断扩大。对于苹果用户而言,这一问题愈发严峻,毕竟512GB和1TB版本的iPhone可不是人人都消费得起的,因此成熟的外置存储方案开

如何在页面调用utility bar并传递参数至lwc组件

1.在app的utility item中添加lwc组件: 2.调用utility bar api的方式有两种: 方法一,通过lwc调用: import {LightningElement,api ,wire } from 'lwc';import { publish, MessageContext } from 'lightning/messageService';import Ca

usaco 1.1 Broken Necklace(DP)

直接上代码 接触的第一道dp ps.大概的思路就是 先从左往右用一个数组在每个点记下蓝或黑的个数 再从右到左算一遍 最后取出最大的即可 核心语句在于: 如果 str[i] = 'r'  ,   rl[i]=rl[i-1]+1, bl[i]=0 如果 str[i] = 'b' ,  bl[i]=bl[i-1]+1, rl[i]=0 如果 str[i] = 'w',  bl[i]=b

业务中14个需要进行A/B测试的时刻[信息图]

在本指南中,我们将全面了解有关 A/B测试 的所有内容。 我们将介绍不同类型的A/B测试,如何有效地规划和启动测试,如何评估测试是否成功,您应该关注哪些指标,多年来我们发现的常见错误等等。 什么是A/B测试? A/B测试(有时称为“分割测试”)是一种实验类型,其中您创建两种或多种内容变体——如登录页面、电子邮件或广告——并将它们显示给不同的受众群体,以查看哪一种效果最好。 本质上,A/B测

【北交大信息所AI-Max2】使用方法

BJTU信息所集群AI_MAX2使用方法 使用的前提是预约到相应的算力卡,拥有登录权限的账号密码,一般为导师组共用一个。 有浏览器、ssh工具就可以。 1.新建集群Terminal 浏览器登陆10.126.62.75 (如果是1集群把75改成66) 交互式开发 执行器选Terminal 密码随便设一个(需记住) 工作空间:私有数据、全部文件 加速器选GeForce_RTX_2080_Ti

【WebGPU Unleashed】1.1 绘制三角形

一部2024新的WebGPU教程,作者Shi Yan。内容很好,翻译过来与大家共享,内容上会有改动,加上自己的理解。更多精彩内容尽在 dt.sim3d.cn ,关注公众号【sky的数孪技术】,技术交流、源码下载请添加微信号:digital_twin123 在 3D 渲染领域,三角形是最基本的绘制元素。在这里,我们将学习如何绘制单个三角形。接下来我们将制作一个简单的着色器来定义三角形内的像素

poj 2976 分数规划二分贪心(部分对总体的贡献度) poj 3111

poj 2976: 题意: 在n场考试中,每场考试共有b题,答对的题目有a题。 允许去掉k场考试,求能达到的最高正确率是多少。 解析: 假设已知准确率为x,则每场考试对于准确率的贡献值为: a - b * x,将贡献值大的排序排在前面舍弃掉后k个。 然后二分x就行了。 代码: #include <iostream>#include <cstdio>#incl

《数据结构(C语言版)第二版》第八章-排序(8.3-交换排序、8.4-选择排序)

8.3 交换排序 8.3.1 冒泡排序 【算法特点】 (1) 稳定排序。 (2) 可用于链式存储结构。 (3) 移动记录次数较多,算法平均时间性能比直接插入排序差。当初始记录无序,n较大时, 此算法不宜采用。 #include <stdio.h>#include <stdlib.h>#define MAXSIZE 26typedef int KeyType;typedef char In

CSP 2023 提高级第一轮 CSP-S 2023初试题 完善程序第二题解析 未完

一、题目阅读 (最大值之和)给定整数序列 a0,⋯,an−1,求该序列所有非空连续子序列的最大值之和。上述参数满足 1≤n≤105 和 1≤ai≤108。 一个序列的非空连续子序列可以用两个下标 ll 和 rr(其中0≤l≤r<n0≤l≤r<n)表示,对应的序列为 al,al+1,⋯,ar​。两个非空连续子序列不同,当且仅当下标不同。 例如,当原序列为 [1,2,1,2] 时,要计算子序列 [

笔记整理—内核!启动!—kernel部分(2)从汇编阶段到start_kernel

kernel起始与ENTRY(stext),和uboot一样,都是从汇编阶段开始的,因为对于kernel而言,还没进行栈的维护,所以无法使用c语言。_HEAD定义了后面代码属于段名为.head .text的段。         内核起始部分代码被解压代码调用,前面关于uboot的文章中有提到过(eg:zImage)。uboot启动是无条件的,只要代码的位置对,上电就工作,kern