Question

在查看documentation几次并花了几个小时尝试各种想法后，我已经决定需要一只手。我正在编写Python软件，逐步完成RADA蛋白质的NAMD模拟，计算出我们感兴趣的不同值。

目前，我的代码逐步完成每个时间步骤，然后逐步执行系统中的每个原子，执行各种分析步骤。

我需要做的是将每种氨基酸合并到它自己的原子中＃34; （氨基酸的单点表示，位于残基的质心）。

我能够实例化新的Atoms，并将它们添加到MDAnalysis Universe吗？我可以使用我的氨基酸质量中心，人类可读的名称等数据填充这些新的原子吗？

类似的东西：

u = MDAnalysis.Universe(psf, coordDcd)
ag = u.selectAtoms(" the atoms in my amino acid ")
amino_acid = MDAnalysis.Atom
amino_acid.pos = ag.centerOfMass()

我知道如何阅读NAMD模拟（.dcd文件）并且所有原子都表现得很好，但最终，我需要将~20个原子转化为一个＆＃34;平均值＆＃34;原子（用于计算简化）。

Answer 1

您可以使用Universe.empty()方法创建一个包含空原子的新宇宙。在以下示例中，我使用测试中的示例文件（data.PSF和data.CD，但您可以使用自己的文件）。

import MDAnalysis as mda
from MDAnalysis.tests import datafiles as data   # only needed for the example files

u = mda.Universe(data.PSF, data.DCD)

# create the "coarse grained" universe with one Atom per Residue
cg = mda.Universe.empty(n_atoms=u.residues.n_residues, n_residues=u.residues.n_residues, 
                        atom_resindex=u.residues.ix, trajectory=True)

# add total mass for each residue
cg.add_TopologyAttr("masses", u.residues.masses)

# add total charge per residue
cg.add_TopologyAttr("charges", u.residues.charges)

# add names etc... so that cg.select_atoms("protein") works
cg.add_TopologyAttr("resnames", u.residues.resnames)
cg.add_TopologyAttr("resnums", u.residues.resnums)

# use the residue center of mass for each residue as the position of the CG particle
centers = u.residues.center_of_mass(compound="residues")
cg.atoms.positions = centers

空荡荡的宇宙是裸露的骨头；使用Universe.add_TopologyAttr()，我们可以添加诸如质量和电荷之类的拓扑属性，这对于粗颗粒而言可能是有意义的。我们使用汇总的每个残基值初始化这些属性（u.residues.masses包含每个残基的质量，而u.atoms.masses包含每个原子的质量）。添加残基名称可启用诸如

的选择

protein = cg.select_atoms("protein")

也可以在cg宇宙中工作。可以以类似方式添加其他拓扑属性。（有关可用属性的文档目前尚不完善，因此也许可以在user mailing list上提问或查看User Guide: Topology System（截至2019年11月，正在进行中））。

我们还可以为CG“原子”分配位置（因为我们用trajectory=True添加了一个空的单帧轨迹），因此典型的分析有效，例如，我们选择的蛋白质残基的回转半径：

rgyr = protein.radius_of_gyration()

以上内容应回答最初发布的问题。

但是，请注意，为轨迹中的每个时间步创建一个Universe可能效率不高。创建粗粒度Universe with a in-memory trajectory attached会很方便，但是当前的Universe.empty()方法（MDAnalysis 0.20.1）不支持此功能-主要是因为没人要求它。实施起来很简单，因此如果您需要它，请在我们的问题跟踪器https://github.com/mdanalysis/mdanalysis/issues中提出问题，并提出功能建议。

如何使用MDAnalysis实例化Atoms？

1 个答案: