Aligning a structure to another¶
align.alignto to align a structure to another.
Last executed: Feb 06, 2020 with MDAnalysis 0.20.2-dev0
Last updated: February 2020
Minimum version of MDAnalysis: 0.17.0
Optional packages for molecular visualisation:
MDAnalysis implements RMSD calculation using the fast QCP algorithm ([The05]) and a rotation matrix R that minimises the RMSD ([LAT09]). Please cite ([The05]) and ([LAT09]) when using the
MDAnalysis.analysis.align module in published work.
import MDAnalysis as mda from MDAnalysis.analysis import align from MDAnalysis.tests.datafiles import CRD, PSF, DCD, DCD2 import nglview as nv
The test files we will be working with here are trajectories of a adenylate kinase (AdK), a phosophotransferase enzyme. ([BDPW09]) The trajectories sample a transition from a closed to an open conformation.
adk_open = mda.Universe(CRD, DCD2) nv.show_mdanalysis(adk_open)
adk_closed = mda.Universe(PSF, DCD) nv.show_mdanalysis(adk_closed.atoms)
Currently, the proteins are not aligned to each other. The difference becomes even more obvious when the closed conformation is compared to the open. Below, we set
adk_open to the last frame and see the relative positions of each protein in a merged Universe.
adk_open.trajectory[-1] # last frame merged = mda.Merge(adk_open.atoms, adk_closed.atoms) nv.show_mdanalysis(merged)
Aligning a structure with align.alignto¶
alignto aligns the mobile AtomGroup to the target AtomGroup by minimising the RMSD. It returns (old_rmsd, new_rmsd).
align.alignto(adk_open, # mobile adk_closed, # reference select='name CA') # selection to operate on
However, positions are set temporarily. If we flip to the first frame of
adk_open and back to the last frame, we can see that it has returned to its original location.
adk_open.trajectory # set to first frame adk_open.trajectory[-1] # set to last frame nv.show_mdanalysis(mda.Merge(adk_open.atoms, adk_closed.atoms))
You can save the aligned positions by writing them out to a PDB file and creating a new Universe.
align.alignto(adk_open, adk_closed, select='name CA') adk_open.atoms.write('aligned.pdb') nv.show_mdanalysis(mda.Universe('aligned.pdb'))
 Oliver Beckstein, Elizabeth J. Denning, Juan R. Perilla, and Thomas B. Woolf. Zipping and Unzipping of Adenylate Kinase: Atomistic Insights into the Ensemble of Open↔Closed Transitions. Journal of Molecular Biology, 394(1):160–176, November 2009. 00107. URL: https://linkinghub.elsevier.com/retrieve/pii/S0022283609011164, doi:10.1016/j.jmb.2009.09.009.
 Richard J. Gowers, Max Linke, Jonathan Barnoud, Tyler J. E. Reddy, Manuel N. Melo, Sean L. Seyler, Jan Domański, David L. Dotson, Sébastien Buchoux, Ian M. Kenney, and Oliver Beckstein. MDAnalysis: A Python Package for the Rapid Analysis of Molecular Dynamics Simulations. Proceedings of the 15th Python in Science Conference, pages 98–105, 2016. 00152. URL: https://conference.scipy.org/proceedings/scipy2016/oliver_beckstein.html, doi:10.25080/Majora-629e541a-00e.
 Pu Liu, Dimitris K. Agrafiotis, and Douglas L. Theobald. Fast determination of the optimal rotational matrix for macromolecular superpositions. Journal of Computational Chemistry, pages n/a–n/a, 2009. URL: http://doi.wiley.com/10.1002/jcc.21439, doi:10.1002/jcc.21439.
 Naveen Michaud-Agrawal, Elizabeth J. Denning, Thomas B. Woolf, and Oliver Beckstein. MDAnalysis: A toolkit for the analysis of molecular dynamics simulations. Journal of Computational Chemistry, 32(10):2319–2327, July 2011. 00778. URL: http://doi.wiley.com/10.1002/jcc.21787, doi:10.1002/jcc.21787.
 Hai Nguyen, David A Case, and Alexander S Rose. NGLview–interactive molecular graphics for Jupyter notebooks. Bioinformatics, 34(7):1241–1242, April 2018. 00024. URL: https://academic.oup.com/bioinformatics/article/34/7/1241/4721781, doi:10.1093/bioinformatics/btx789.
 Douglas L. Theobald. Rapid calculation of RMSDs using a quaternion-based characteristic polynomial. Acta Crystallographica Section A Foundations of Crystallography, 61(4):478–480, July 2005. 00127. URL: http://scripts.iucr.org/cgi-bin/paper?S0108767305015266, doi:10.1107/S0108767305015266.