๐Ÿ†• Theory๏ƒ

This chapter describes HELM support in OEChem TK, which currently enables the representation of linear and macrocyclic peptides through their amino acid sequences. To become familiar with the notation, it is highly recommended to read the original publication [Zhang-2012].

Warning

HELM representation is useful for encoding peptides with complex structures, however we DO NOT RECOMMEND using it as a unique identifier or as a primary format for storing peptide information.

The reason for the former is that a single peptide can be represented by multiple valid HELM strings. While the algorithm implemented in OEChem TK for generating HELM strings has been extensively tested for consistency, it does not guarantee that the same peptide will always be represented by the same HELM string. Furthermore, the capability to modify and extend the monomer set by the user can lead to changes in the generated HELM strings for the same peptide even with the same algorithm.

The reason for the latter is that HELM is not a self contained format. It relies on external monomer definitions which can be modified by the user at will. Storing peptide purely as HELM strings without knowing which monomer set was used to generate them can lead to losing the ability to reconstruct the original atomic-level molecular structure. We recommend using other representations such as SMILES as your primary format and generating HELM strings on demand for specific use cases.

Monomers๏ƒ

Monomers are the fundamental building blocks of peptides. To define a monomer for use in OEChem TK, the following information must be provided:

  • the monomerโ€™s structure as SMILES with explicitly defined connection points

  • the polymer type of the monomer (Peptide or Chem) โ€” OEChem TK does not currently support RNA monomers

  • the monomer type (Backbone or Terminal)

  • a unique mnemonic code that unambiguously identifies the monomer

Assigning a descriptive name (for example, an IUPAC name) to a monomer is optional.

Monomers with Peptide Polymer Type๏ƒ

Monomers are characterized as Peptide if they can appear in the PEPTIDE polymer chain.

Most of the monomers in this category are backbone amino-acids with the following constraints:

  • R1 (N-terminus) defines connection point on the nitrogen atom of the amino group

  • R2 (C-terminus) defines connection point on the carbon atom of the carboxyl group

  • R3 optionally can define side-chain connections

  • having ฮฑ, ฮฒ, ฮณ or ฮด backbone

Table 1. Backbone peptide monomers๏ƒ
../_images/NMePhe.svg ../_images/Abz.svg ../_images/App.svg
Table 2. Backbone peptide monomers with side-chain connections๏ƒ
../_images/C.svg ../_images/Asu.svg ../_images/Orn.svg

However, non-peptidic monomers can also belong to the Peptide category if

  • it has only one R1 amino group connection point โ€” used at the C-terminus

  • it has only one R2 carboxyl group connection point โ€” used at the N-terminus

Table 3. Terminal peptide monomers๏ƒ
../_images/NMe2.svg ../_images/Bua.svg

C-terminal

N-terminal

Monomers with Chem Polymer Type๏ƒ

Non-peptidic monomers that do not fit into the Peptide category, with arbitrary connection points.

Monomers with Unspecified Stereochemistry๏ƒ

The current implementation of OEChem TK does not support monomers with unspecified stereochemistry. Atom or bond stereo must be explicitly defined in the monomer structure for it to be successfully mapped to a molecule when generating a HELM string.

See also: Handling Molecules with Unspecified Stereochemistry section below for handling molecules with unspecified stereochemistry when generating HELM strings.

Defining Monomer Connection Points๏ƒ

When adding a monomer to a monomer set, an equivalency test is performed based on the canonical SMILES of the monomer without considering the R-groups. The table below shows the same cysteine monomer defined with different connection points. Only one of these can be added to the monomer set, since they are considered equivalent. It is recommended to define monomers with as many connection points as possible to maximize their versatility in peptide design. For example, the cysteine monomer defined with three connection points (see Table 4, example D) can be used in a wider variety of contexts.

Table 4. Cysteine monomer defined with various connection points.๏ƒ

A

B

C

D - RECOMMENDED

../_images/cysteine-R1.svg ../_images/cysteine-R2.svg ../_images/cysteine-R12.svg ../_images/cysteine-R123.svg

can be used only as C-terminal

can be used only as N-terminal

can be used as any terminal or backbone

can be used as terminal, backbone or cross-linker

OEChem TK Built-in Monomer Sets๏ƒ

There are two built-in monomer sets supported by OEChem TK:

  • โ€œStandardโ€ โ€” contains 20 standard amino acids in the code-set named Standard

  • โ€œOpenEyeโ€ โ€” contains about 260 diverse monomers defined in five code-sets named: OpenEye, Standard, PDB, and ChEMBL.

Warning

These built-in monomer sets are provided to test and demonstrate support for diverse peptide structures and to serve as guidance for customers building their own private monomer sets. We strongly discourage users from relying solely on the OpenEye-provided monomer sets, as they are subject to change. No guarantee is made that monomers will not be added, removed, or modified within these sets. Instead, users should maintain their own version, allowing them to add the monomers needed to support their research and to adopt the naming conventions preferred by their organization.

See also

HELM๏ƒ

The order of monomers in a HELM peptide chain is directional. Each chain lists its monomers from the N-terminus to the C-terminus. In the example below, the peptide begins with proline (P) at the N-terminus and ends with glutamic acid (E) at the C-terminus.

../_images/PEPTIDE.svg

PEPTIDE1{P.E.P.T.I.D.E}$$$$๏ƒ

Each dot in the chain represents the amide bond formed between adjacent amino acids. It indicates that the carboxyl connection point (R2) of the monomer on the left is connected to the amino (R1) connection point of the monomer on the right.

Additional bonds between the monomers are defined in the second component of the HELM string. In the oxytocin example below, the first and sixth monomers of the same PEPTIDE1 polymer chain are connected. In this case, the connection is a disulfide bond formed through the R3 side-chain connection points of the two cystine monomers. This connection is encoded as PEPTIDE1,PEPTIDE1,1:R3-6:R3

oxytocin

PEPTIDE1{C.Y.I.Q.N.C.P.L.G.[am]}$PEPTIDE1,PEPTIDE1,1:R3-6:R3$$$$$

../_images/Oxytocin.svg

Peptides can also be cyclized by a lactam bridge commonly connecting side chains (e.g., Lys-Asp or Lys-Glu).

Multiple connections are separated by |. Conotoxins such as ฮฑ-conotoxin PnIA, depicted below, generally feature several disulfide bonds that stabilize their three-dimensional structure.

ฮฑ-conotoxin PnIA

PEPTIDE1{G.C.C.S.L.P.P.C.A.A.N.N.P.D.Y.C.[am]}$PEPTIDE1,PEPTIDE1,2:R3-8:R3|PEPTIDE1,PEPTIDE1,3:R3-16:R3$$$

../_images/alpha_conotoxin_PnIA.svg

HELM Parsing๏ƒ

In OEChem, a HELM string can be parsed using the OEHelmToMol function.

Embedded SMILES๏ƒ

OEChem TK supports reading of HELM with embedded SMILES in various styles. The usage of embedded SMILES is necessary for defining a peptide that contains an unnamed monomer.

The following three formats define the same eledoisin structure

  1. The most compact representation with R-group style:

    PEPTIDE1{[N1[C@@H](CCC1=O)C([R2])=O].P.S.K.D.A.F.I.G.L.M.[am]}$$$$

  2. Representation using atom mapping to mark attachment points

    PEPTIDE1{[N1[C@@H](CCC1=O)C([*:2])=O].P.S.K.D.A.F.I.G.L.M.[am]}$$$$

  3. CXSMILES style representation

    PEPTIDE1{[N1[C@@H](CCC1=O)C(*)=O |$;;;;;;;_R2;$|].P.S.K.D.A.F.I.G.L.M.[am]}$$$$

../_images/eledoisin.svg

Eledoisin๏ƒ

In the case of an embedded SMILES, the attachment-point mapping must follow the same rules as it would if the structure were defined as a named peptide monomer.

../_images/embedded-monomer-CyclosporinA.svg

[R2]C(=O)[C@H]([C@@H]([C@H](C)C/C=C/C)O)N([R1])C๏ƒ

  • R1 (N-terminus) in the embedded SMILES is defined as the connection point on the nitrogen atom of the amino group

  • R2 (C-terminus) in the embedded SMILES is defined as the connection point on the carbon atom of the carboxyl group

See the unnamed fragment depicted on the left and how it is used to define CyclosporinA below.

Cyclosporin A

PEPTIDE1{[Abu].[Sar].[MeL].V.[MeL].A.a.[MeL].[MeL].[MeV].[[R2]C(=O)[C@H]([C@@H]([C@H](C)C/C=C/C)O)N([R1])C]}$PEPTIDE1,PEPTIDE1,11:R2-1:R1$$$

../_images/CyclosporinA.svg

Helm Generation๏ƒ

In OEChem, a HELM string can be generated using the OEMolToHelm function.

Tautomer Handling๏ƒ

When mapping monomers onto a molecule to determine its sequence, bond orders in aromatic rings and certain functional groups (such as guanidine) are handled in a way that enables successful generation of HELM strings across different tautomeric forms.

Table 4. Tautomer handling in HELM generation.๏ƒ
../_images/helm-generation-Tautomer-A.svg ../_images/helm-generation-Tautomer-B.svg

input SMILES

c1c([nH]cn1)C[C@@H](C(=O)N[C@@H](CCCNC(=N)N)C(=O)O)N

c1c(nc[nH]1)C[C@@H](C(=O)N[C@@H](CCCN=C(N)N)C(=O)O)N

output HELM

PEPTIDE1{H.R}$$$$

PEPTIDE1{H.R}$$$$

Handling Charged Molecules๏ƒ

It is strongly recommended to define monomers in their neutral form, enabling matching of both charged and uncharged peptides when generating HELM representations.

Table 5. Charge handling in HELM generation.๏ƒ

Structure A

Structure B

Structure C

../_images/helm-generation-charge-A.svg ../_images/helm-generation-charge-B.svg ../_images/helm-generation-charge-C.svg

PEPTIDE1{D.A.R.K}$$$$

PEPTIDE1{D.A.R.K}$$$$

PEPTIDE1{D.A.R.K}$$$$

If a monomer contains a positively or negatively charged atom, that atom can only be matched to an atom in a peptide that has the same formal charge. For example, if protonated lysine K_protonated ([NH3+]CCCC[C@@H](C(=O)[OH:2])N[H:1]) is added to the โ€œStandardโ€ monomer set with the code K+, this monomer will match both structure B and C in the example below.

Table 6. Charge handling in HELM generation.๏ƒ

Structure A

Structure B

Structure C

../_images/helm-generation-charge-A-K%2B.svg ../_images/helm-generation-charge-B-K%2B.svg ../_images/helm-generation-charge-C-K%2B.svg

PEPTIDE1{D.A.R.K}$$$$

PEPTIDE1{D.A.R.[K+]}$$$$

PEPTIDE1{D.A.R.[K+]}$$$$

Handling Molecules with Unspecified Stereochemistry๏ƒ

Monomer detection is stereochemistry-sensitive when generating a peptide sequence. Because the current implementation does not support monomers with unspecified stereochemistry, HELM generation automatically fails for any structure with undefined stereo. However, it is still possible to generate a valid HELM string by disabling the stereochemistry check and allowing monomers with unspecified stereo to be encoded as embedded SMILES in the generated HELM.

See also Helm Generation Options section of smiles2helm for more explanation and examples.