MoleculeAlign

MoleculeAlign[ref,mol]

returns a version of mol that is aligned with reference molecule ref.

MoleculeAlign[ref,mol,{r1m1,r2m2,}]

aligns atoms in mol with index mi to the atoms in ref with index ri.

MoleculeAlign[ref,mol,patt]

uses the molecule pattern patt to find an atom mapping between ref and mol.

MoleculeAlign[ref,{mol1,mol2,},patt]

aligns each of the moli with ref.

MoleculeAlign[ref,mols,patt,prop]

aligns the molecules and returns the property prop of the alignment.

Details and Options

  • MoleculeAlign returns a Molecule object with the AtomCoordinates option explicitly set.
  • The atom mapping {r1m1,r2m2,} can indicate a subset of the atoms, which will be aligned.
  • When at least three non-colinear atoms are included in the mapping, the alignment will be unique. When colinear atoms are provided the orientation about that axis is undefined.
  • When no atom mapping is provided, ref must be present as a substructure in mol. If ref occurs multiple times, the first substructure match as found by FindMoleculeSubstructure is used.
  • When multiple molecules are aligned, the atom mapping can be provided in the form {r1{},r2{},} where the right-hand sides are lists of atom indices for each molecule to be aligned.
  • Multiple molecules can be aligned either pairwise to the reference molecule with MethodAutomatic, or jointly to each other and the reference molecule with Method"Joint".
  • If ref does not contain 3D coordinates, they are generated automatically prior to alignment.
  • The following options can be specified:
  • WeightsAutomaticweights for the atomic positions
    MethodAutomaticthe method to use
  • prop can be one of the following:
  • "Molecule"the aligned molecules
    "Transformation"the transformation functions used
    "Error"the alignment root mean squared error

Examples

open allclose all

Basic Examples  (5)

Align one conformation to another:

Align one conformation of water to another:

The coordinates have been rotated and translated into alignment with the reference:

Align two molecules:

Align two molecules with an explicit atom mapping:

Align two molecules using a MoleculePattern to determine the atom mapping:

Scope  (3)

MoleculeAlign can be used to align several conformers. Start by generating 10 conformations of adenosine with randomly oriented conformations:

Align the molecules and display them together:

Provide a pattern to specify which parts of the conformers should be aligned:

Menthol, a monoterpene, is a substructure of jujenol, a sesquiterpene. Here are the molecules:

Align jujenol to menthol:

Show the aligned result:

Align a resveratrol, an endocrine disruptor, to estradiol, a hormone, using an explicitly specified atom mapping. Here are the molecules:

Align the phenol ring carbon atoms and two of the hydroxyl oxygen atoms of resveratrol to the A-ring carbon atoms and the two hydroxyl oxygen atoms of estradiol:

Show the aligned molecules:

Compute the RMS alignment error:

Options  (2)

Method  (1)

When aligning multiple molecules, the default method aligns each of the moving molecules pairwise with the reference molecule. Here are three conformations of butane:

The first has the anti conformation and will serve as the reference molecule. The other two are the gauche conformations. First align with the default pairwise method:

Here are the RMS deviations from the reference molecule:

Now, align the molecules simultaneously by minimizing the joint sum of squared deviations:

Note that the second and third molecules are now more closely aligned to each other, as well as to the reference molecule, than in the pairwise alignment.

Here is the joint RMS deviation:

Weights  (1)

Atom weights allow specified atoms to have more influence on the alignment. Here are three molecules to demonstrate this:

Here is the alignment with the default atom weights:

Increase the weight on the second atom:

Wolfram Research (2021), MoleculeAlign, Wolfram Language function, https://reference.wolfram.com/language/ref/MoleculeAlign.html.

Text

Wolfram Research (2021), MoleculeAlign, Wolfram Language function, https://reference.wolfram.com/language/ref/MoleculeAlign.html.

BibTeX

@misc{reference.wolfram_2021_moleculealign, author="Wolfram Research", title="{MoleculeAlign}", year="2021", howpublished="\url{https://reference.wolfram.com/language/ref/MoleculeAlign.html}", note=[Accessed: 26-September-2021 ]}

BibLaTeX

@online{reference.wolfram_2021_moleculealign, organization={Wolfram Research}, title={MoleculeAlign}, year={2021}, url={https://reference.wolfram.com/language/ref/MoleculeAlign.html}, note=[Accessed: 26-September-2021 ]}

CMS

Wolfram Language. 2021. "MoleculeAlign." Wolfram Language & System Documentation Center. Wolfram Research. https://reference.wolfram.com/language/ref/MoleculeAlign.html.

APA

Wolfram Language. (2021). MoleculeAlign. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/MoleculeAlign.html