OECapNTerminiΒΆ

Attention

PRELIMINARY-IMAGE This is a preliminary API until 2020.Apr and may be improved based on user feedback. It is currently available in C++ and Python.

bool OECapNTermini(OEChem::OEMolBase &mol,
                   const OESystem::OEUnaryPredicate<OEChem::OEAtomBase> &exclude=OEChem::OEIsFalseAtom(),
                   const bool deleteClashingSolvent=true)

This function caps all the N-termini in a protein with an acetyl (OEResidueIndex_ACE) group. The acetyl is built assuming ideal backbone geometry, the \(\phi\) value of the residue being capped is set to match the \(\phi\) value of the previous residue (if it’s a single residue in space, a value of \(-90.0^{\circ}\) degrees, is set which is in the allowed range for both an alpha helix and beta sheet secondary structure).

A predicate can be passed to exclude certain termini from being capped, such as the true biological N-terminal residue since this should be charged to match the biophysical experiments.

The function by default can delete a solvent molecule if it clashes with where the cap is being built, but this behavior can be turned off. A solvent molecule blocking the protein backbone is not very realistic and the solvent molecule is likely placed to account for density that should have belonged to the residue not resolved in the crystal structure.

The cap will be be numbered with -1 from the residue number of the terminal residue.

Note

If the backbone of the residue being capped is incomplete, or a cap can not be built because of a clash, the residue being capped will be turned into the cap by removing atoms in the residue. A clash can either be into a non-solvent or into a solvent with the deleteClashingSolvent option set to false.

See also