biotite 1.0.1__cp310-cp310-win_amd64.whl → 1.2.0__cp310-cp310-win_amd64.whl

biotite/structure/io/general.py

@@ -34,7 +34,7 @@ def load_structure(file_path, template=None, **kwargs):
         The path to structure file.
     template : AtomArray or AtomArrayStack or file-like object or str, optional
         Only required when reading a trajectory file.
-    kwargs
+    **kwargs
         Additional parameters will be passed to either the
         :func:`get_structure()` or :func:`read()` method of the file
         object.
@@ -146,7 +146,7 @@ def save_structure(file_path, array, **kwargs):
         The path to structure file.
     array : AtomArray or AtomArrayStack
         The structure to be saved.
-    kwargs
+    **kwargs
         Additional parameters will be passed to the respective `set_structure`
         method.
 

biotite/structure/io/gro/file.py

@@ -48,7 +48,6 @@ class GROFile(TextFile):
     >>> file = GROFile()
     >>> file.set_structure(array_stack_mod)
     >>> file.write(os.path.join(path_to_directory, "1l2y_mod.gro"))
-
     """
 
     def get_model_count(self):
@@ -89,13 +88,13 @@ class GROFile(TextFile):
             The return type depends on the `model` parameter.
         """
 
-        def get_atom_line_i(model_start_i, model_atom_counts):
+        def _get_atom_line_i(model_start_i, model_atom_counts):
             """
             Helper function to get the indices of all atoms for a model
             """
             return np.arange(model_start_i + 1, model_start_i + 1 + model_atom_counts)
 
-        def set_box_dimen(box_param):
+        def _set_box_dimen(box_param):
             """
             Helper function to create the box vectors from the values
             in the GRO file
@@ -145,7 +144,7 @@ class GROFile(TextFile):
 
             # Line indices for annotation determination is determined
             # from model 1
-            annot_i = get_atom_line_i(model_start_i[0], length)
+            annot_i = _get_atom_line_i(model_start_i[0], length)
         else:
             if model == 0:
                 raise ValueError("The model index must not be 0")
@@ -160,7 +159,7 @@ class GROFile(TextFile):
             length = model_atom_counts[model - 1]
             array = AtomArray(length)
 
-            annot_i = get_atom_line_i(model_start_i[model - 1], length)
+            annot_i = _get_atom_line_i(model_start_i[model - 1], length)
 
         # Replace empty strings for elements with guessed types
         # i is index in array, line_i is line index
@@ -183,11 +182,11 @@ class GROFile(TextFile):
             # Box is stored in last line (after coordinates)
             box_i = atom_i[-1] + 1
             box_param = [float(e) * 10 for e in self.lines[box_i].split()]
-            array.box = set_box_dimen(box_param)
+            array.box = _set_box_dimen(box_param)
 
         elif isinstance(array, AtomArrayStack):
             for m in range(len(model_start_i)):
-                atom_i = get_atom_line_i(model_start_i[m], model_atom_counts[m])
+                atom_i = _get_atom_line_i(model_start_i[m], model_atom_counts[m])
                 for i, line_i in enumerate(atom_i):
                     line = self.lines[line_i]
                     array.coord[m, i, 0] = float(line[20:28]) * 10
@@ -196,7 +195,7 @@ class GROFile(TextFile):
                 # Box is stored in last line (after coordinates)
                 box_i = atom_i[-1] + 1
                 box_param = [float(e) * 10 for e in self.lines[box_i].split()]
-                box = set_box_dimen(box_param)
+                box = _set_box_dimen(box_param)
                 # Create a box in the stack if not already existing
                 # and the box is not a dummy
                 if box is not None:
@@ -311,7 +310,7 @@ class GROFile(TextFile):
 
             for i in range(array.stack_depth()):
                 self.lines.append(
-                    f"Generated by Biotite at {datetime.now()}, model={i+1}"
+                    f"Generated by Biotite at {datetime.now()}, model={i + 1}"
                 )
                 self.lines.append(str(array.array_length()))
 

biotite/structure/io/mol/convert.py

@@ -52,7 +52,7 @@ def set_structure(
     ----------
     mol_file : MOLFile
         The MOL file.
-    array : AtomArray
+    atoms : AtomArray
         The array to be saved into this file.
         Must have an associated :class:`BondList`.
         Bond type fallback for the *Bond block*, if a

biotite/structure/io/mol/ctab.py

@@ -19,24 +19,19 @@ from biotite.file import InvalidFileError
 from biotite.structure.atoms import AtomArray, AtomArrayStack
 from biotite.structure.bonds import BondList, BondType
 from biotite.structure.error import BadStructureError
+from biotite.structure.io.util import number_of_integer_digits
 
 BOND_TYPE_MAPPING = {
     1: BondType.SINGLE,
     2: BondType.DOUBLE,
     3: BondType.TRIPLE,
+    4: BondType.AROMATIC,
     5: BondType.ANY,
-    6: BondType.SINGLE,
-    7: BondType.DOUBLE,
+    6: BondType.AROMATIC_SINGLE,
+    7: BondType.AROMATIC_DOUBLE,
     8: BondType.ANY,
 }
-BOND_TYPE_MAPPING_REV = {
-    BondType.SINGLE: 1,
-    BondType.DOUBLE: 2,
-    BondType.TRIPLE: 3,
-    BondType.AROMATIC_SINGLE: 1,
-    BondType.AROMATIC_DOUBLE: 2,
-    BondType.ANY: 8,
-}
+BOND_TYPE_MAPPING_REV = {v: k for k, v in BOND_TYPE_MAPPING.items()}
 
 CHARGE_MAPPING = {0: 0, 1: 3, 2: 2, 3: 1, 5: -1, 6: -2, 7: -3}
 CHARGE_MAPPING_REV = {val: key for key, val in CHARGE_MAPPING.items()}
@@ -56,7 +51,7 @@ def read_structure_from_ctab(ctab_lines):
     ----------
     ctab_lines : lines of str
         The lines containing the *ctab*.
-        Must begin with the *counts* line and end with the `M END` line
+        Must begin with the *counts* line and end with the `M END` line.
 
     Returns
     -------
@@ -71,7 +66,6 @@ def read_structure_from_ctab(ctab_lines):
     `<https://discover.3ds.com/sites/default/files/2020-08/biovia_ctfileformats_2020.pdf>`_.
 
     .. footbibliography::
-
     """
     match _get_version(ctab_lines[0]):
         case "V2000":
@@ -119,11 +113,10 @@ def write_structure_to_ctab(atoms, default_bond_type=BondType.ANY, version=None)
     `<https://discover.3ds.com/sites/default/files/2020-08/biovia_ctfileformats_2020.pdf>`_.
 
     .. footbibliography::
-
     """
     if isinstance(atoms, AtomArrayStack):
         raise TypeError(
-            "An 'AtomArrayStack' was given, " "but only a single model can be written"
+            "An 'AtomArrayStack' was given, but only a single model can be written"
         )
     if atoms.bonds is None:
         raise BadStructureError("Input AtomArray has no associated BondList")
@@ -141,8 +134,7 @@ def write_structure_to_ctab(atoms, default_bond_type=BondType.ANY, version=None)
                 atoms.array_length(), atoms.bonds.get_bond_count()
             ):
                 raise ValueError(
-                    "The given number of atoms or bonds is too large "
-                    "for V2000 format"
+                    "The given number of atoms or bonds is too large for V2000 format"
                 )
             return _write_structure_to_ctab_v2000(atoms, default_bond_type)
         case "V3000":
@@ -174,7 +166,7 @@ def _read_structure_from_ctab_v2000(ctab_lines):
             charge = CHARGE_MAPPING.get(int(line[36:39]))
             if charge is None:
                 warnings.warn(
-                    f"Cannot handle MDL charge type {int(line[36 : 39])}, "
+                    f"Cannot handle MDL charge type {int(line[36:39])}, "
                     f"0 is used instead"
                 )
                 charge = 0
@@ -194,7 +186,7 @@ def _read_structure_from_ctab_v2000(ctab_lines):
         bond_type = BOND_TYPE_MAPPING.get(int(line[6:9]))
         if bond_type is None:
             warnings.warn(
-                f"Cannot handle MDL bond type {int(line[6 : 9])}, "
+                f"Cannot handle MDL bond type {int(line[6:9])}, "
                 f"BondType.ANY is used instead"
             )
             bond_type = BondType.ANY
@@ -247,8 +239,7 @@ def _read_structure_from_ctab_v3000(ctab_lines):
         bond_type = BOND_TYPE_MAPPING.get(v30_type)
         if bond_type is None:
             warnings.warn(
-                f"Cannot handle MDL bond type {v30_type}, "
-                f"BondType.ANY is used instead"
+                f"Cannot handle MDL bond type {v30_type}, BondType.ANY is used instead"
             )
             bond_type = BondType.ANY
         bond_array[i, 0] = v30_atom_indices[v30_atom_index_1]
@@ -307,10 +298,17 @@ def _write_structure_to_ctab_v2000(atoms, default_bond_type):
         "  0     0  0  0  0  0  0  1 V2000"
     )
 
+    for i, coord_name in enumerate(["x", "y", "z"]):
+        n_coord_digits = number_of_integer_digits(atoms.coord[:, i])
+        if n_coord_digits > 5:
+            raise BadStructureError(
+                f"5 pre-decimal columns for {coord_name}-coordinates are "
+                f"available, but array would require {n_coord_digits}"
+            )
     atom_lines = [
-        f"{atoms.coord[i,0]:>10.4f}"
-        f"{atoms.coord[i,1]:>10.4f}"
-        f"{atoms.coord[i,2]:>10.4f}"
+        f"{atoms.coord[i, 0]:>10.4f}"
+        f"{atoms.coord[i, 1]:>10.4f}"
+        f"{atoms.coord[i, 2]:>10.4f}"
         f" {atoms.element[i].capitalize():3}"
         f"{0:>2}"  # Mass difference -> unused
         f"{CHARGE_MAPPING_REV.get(charge[i], 0):>3d}"
@@ -321,7 +319,7 @@ def _write_structure_to_ctab_v2000(atoms, default_bond_type):
 
     default_bond_value = BOND_TYPE_MAPPING_REV[default_bond_type]
     bond_lines = [
-        f"{i+1:>3d}{j+1:>3d}"
+        f"{i + 1:>3d}{j + 1:>3d}"
         f"{BOND_TYPE_MAPPING_REV.get(bond_type, default_bond_value):>3d}"
         + f"{0:>3d}"
         * 4
@@ -337,7 +335,7 @@ def _write_structure_to_ctab_v2000(atoms, default_bond_type):
     ):
         charge_lines.append(
             f"M  CHG{len(batch):>3d}"
-            + "".join(f" {atom_i+1:>3d} {c:>3d}" for atom_i, c in batch)
+            + "".join(f" {atom_i + 1:>3d} {c:>3d}" for atom_i, c in batch)
         )
 
     return [counts_line] + atom_lines + bond_lines + charge_lines + ["M  END"]
@@ -351,12 +349,19 @@ def _write_structure_to_ctab_v3000(atoms, default_bond_type):
 
     counts_line = f"COUNTS {atoms.array_length()} {atoms.bonds.get_bond_count()} 0 0 0"
 
+    for i, coord_name in enumerate(["x", "y", "z"]):
+        n_coord_digits = number_of_integer_digits(atoms.coord[:, i])
+        if n_coord_digits > 5:
+            raise BadStructureError(
+                f"5 pre-decimal columns for {coord_name}-coordinates are "
+                f"available, but array would require {n_coord_digits}"
+            )
     atom_lines = [
         f"{i + 1}"
         f" {_quote(atoms.element[i].capitalize())}"
-        f" {atoms.coord[i,0]:.4f}"
-        f" {atoms.coord[i,1]:.4f}"
-        f" {atoms.coord[i,2]:.4f}"
+        f" {atoms.coord[i, 0]:.4f}"
+        f" {atoms.coord[i, 1]:.4f}"
+        f" {atoms.coord[i, 2]:.4f}"
         # 'aamap' is unused
         f" 0"
         f" {_to_property(charges[i])}"

biotite/structure/io/mol/mol.py

@@ -122,7 +122,7 @@ class MOLFile(TextFile):
 
         Parameters
         ----------
-        array : AtomArray
+        atoms : AtomArray
             The array to be saved into this file.
             Must have an associated :class:`BondList`.
         default_bond_type : BondType, optional

biotite/structure/io/mol/sdf.py

@@ -78,7 +78,6 @@ class Metadata(MutableMapping):
     >>> print(metadata[Metadata.Key(number=42, name="bar")])
     dolor sit amet,
     consectetur
-
     """
 
     @dataclass(frozen=True, kw_only=True)
@@ -144,12 +143,18 @@ class Metadata(MutableMapping):
         @staticmethod
         def deserialize(text):
             """
-            Create an object by deserializing the given text content.
+            Create a :class:`Metadata.Key` object by deserializing the given text
+            content.
 
             Parameters
             ----------
-            content : str
+            text : str
                 The content to be deserialized.
+
+            Returns
+            -------
+            key : Metadata.Key
+                The parsed key.
             """
             # Omit the leading '>'
             key_components = text[1:].split()
@@ -207,12 +212,17 @@ class Metadata(MutableMapping):
     @staticmethod
     def deserialize(text):
         """
-        Create an object by deserializing the given text content.
+        Create a :class:`Metadata` objtect by deserializing the given text content.
 
         Parameters
         ----------
-        content : str
+        text : str
             The content to be deserialized.
+
+        Returns
+        -------
+        metadata : Metadata
+            The parsed metadata.
         """
         metadata = {}
         current_key = None
@@ -314,19 +324,19 @@ class SDRecord:
     >>> record = SDRecord(header=Header(mol_name="ALA", dimensions="3D"))
     >>> record.set_structure(atoms)
     >>> print(record.get_structure())
-            0             N        -0.970    0.490    1.500
-            0             C         0.260    0.420    0.690
-            0             C        -0.090    0.020   -0.720
-            0             O        -1.060   -0.680   -0.920
-            0             C         1.200   -0.620    1.300
-            0             O         0.660    0.440   -1.740
-            0             H        -1.380   -0.420    1.480
-            0             H        -0.680    0.660    2.450
-            0             H         0.750    1.390    0.680
-            0             H         1.460   -0.330    2.320
-            0             H         0.720   -1.590    1.310
-            0             H         2.110   -0.680    0.700
-            0             H         0.440    0.180   -2.650
+                0             N        -0.966    0.493    1.500
+                0             C         0.257    0.418    0.692
+                0             C        -0.094    0.017   -0.716
+                0             O        -1.056   -0.682   -0.923
+                0             C         1.204   -0.620    1.296
+                0             O         0.661    0.439   -1.742
+                0             H        -1.383   -0.425    1.482
+                0             H        -0.676    0.661    2.452
+                0             H         0.746    1.392    0.682
+                0             H         1.459   -0.330    2.316
+                0             H         0.715   -1.594    1.307
+                0             H         2.113   -0.676    0.697
+                0             H         0.435    0.182   -2.647
     >>> # Add the record to an SD file
     >>> file = SDFile()
     >>> file["ALA"] = record
@@ -335,19 +345,19 @@ class SDRecord:
                         3D
     <BLANKLINE>
      13 12  0     0  0  0  0  0  0  1 V2000
-       -0.9700    0.4900    1.5000 N   0  0  0  0  0  0  0  0  0  0  0  0
-        0.2600    0.4200    0.6900 C   0  0  0  0  0  0  0  0  0  0  0  0
-       -0.0900    0.0200   -0.7200 C   0  0  0  0  0  0  0  0  0  0  0  0
-       -1.0600   -0.6800   -0.9200 O   0  0  0  0  0  0  0  0  0  0  0  0
-        1.2000   -0.6200    1.3000 C   0  0  0  0  0  0  0  0  0  0  0  0
-        0.6600    0.4400   -1.7400 O   0  0  0  0  0  0  0  0  0  0  0  0
-       -1.3800   -0.4200    1.4800 H   0  0  0  0  0  0  0  0  0  0  0  0
-       -0.6800    0.6600    2.4500 H   0  0  0  0  0  0  0  0  0  0  0  0
-        0.7500    1.3900    0.6800 H   0  0  0  0  0  0  0  0  0  0  0  0
-        1.4600   -0.3300    2.3200 H   0  0  0  0  0  0  0  0  0  0  0  0
-        0.7200   -1.5900    1.3100 H   0  0  0  0  0  0  0  0  0  0  0  0
-        2.1100   -0.6800    0.7000 H   0  0  0  0  0  0  0  0  0  0  0  0
-        0.4400    0.1800   -2.6500 H   0  0  0  0  0  0  0  0  0  0  0  0
+       -0.9660    0.4930    1.5000 N   0  0  0  0  0  0  0  0  0  0  0  0
+        0.2570    0.4180    0.6920 C   0  0  0  0  0  0  0  0  0  0  0  0
+       -0.0940    0.0170   -0.7160 C   0  0  0  0  0  0  0  0  0  0  0  0
+       -1.0560   -0.6820   -0.9230 O   0  0  0  0  0  0  0  0  0  0  0  0
+        1.2040   -0.6200    1.2960 C   0  0  0  0  0  0  0  0  0  0  0  0
+        0.6610    0.4390   -1.7420 O   0  0  0  0  0  0  0  0  0  0  0  0
+       -1.3830   -0.4250    1.4820 H   0  0  0  0  0  0  0  0  0  0  0  0
+       -0.6760    0.6610    2.4520 H   0  0  0  0  0  0  0  0  0  0  0  0
+        0.7460    1.3920    0.6820 H   0  0  0  0  0  0  0  0  0  0  0  0
+        1.4590   -0.3300    2.3160 H   0  0  0  0  0  0  0  0  0  0  0  0
+        0.7150   -1.5940    1.3070 H   0  0  0  0  0  0  0  0  0  0  0  0
+        2.1130   -0.6760    0.6970 H   0  0  0  0  0  0  0  0  0  0  0  0
+        0.4350    0.1820   -2.6470 H   0  0  0  0  0  0  0  0  0  0  0  0
       1  2  1  0  0  0  0
       1  7  1  0  0  0  0
       1  8  1  0  0  0  0
@@ -362,6 +372,7 @@ class SDRecord:
       6 13  1  0  0  0  0
     M  END
     $$$$
+    <BLANKLINE>
     """
 
     def __init__(self, header=None, ctab=None, metadata=None):
@@ -422,19 +433,23 @@ class SDRecord:
             self._metadata = Metadata(metadata)
         else:
             raise TypeError(
-                "Expected 'Metadata' or Mapping, "
-                f"but got '{type(metadata).__name__}'"
+                f"Expected 'Metadata' or Mapping, but got '{type(metadata).__name__}'"
             )
 
     @staticmethod
     def deserialize(text):
         """
-        Create an object by deserializing the given text content.
+        Create an :class:`SDRecord` by deserializing the given text content.
 
         Parameters
         ----------
-        content : str
+        text : str
             The content to be deserialized.
+
+        Returns
+        -------
+        record : SDRecord
+            The parsed record.
         """
         lines = text.splitlines()
         ctab_end = _get_ctab_stop(lines)
@@ -493,7 +508,7 @@ class SDRecord:
 
         Parameters
         ----------
-        array : AtomArray
+        atoms : AtomArray
             The array to be saved into this file.
             Must have an associated :class:`BondList`.
         default_bond_type : BondType, optional
@@ -538,11 +553,18 @@ class SDFile(File, MutableMapping):
     :class:`SDRecord` object via :func:`get_structure()` or
     :func:`set_structure()`, respectively.
 
+    Parameters
+    ----------
+    records : dict (str -> SDRecord), optional
+        The initial records of the file.
+        Maps the record names to the corresponding :class:`SDRecord` objects.
+        By default no initial records are added.
+
     Attributes
     ----------
     record : CIFBlock
         The sole record of the file.
-        If the file contains multiple blocks, an exception is raised.
+        If the file contains multiple records, an exception is raised.
 
     Examples
     --------
@@ -678,19 +700,19 @@ class SDFile(File, MutableMapping):
     <BLANKLINE>
     <BLANKLINE>
      13 12  0     0  0  0  0  0  0  1 V2000
-       -0.9700    0.4900    1.5000 N   0  0  0  0  0  0  0  0  0  0  0  0
-        0.2600    0.4200    0.6900 C   0  0  0  0  0  0  0  0  0  0  0  0
-       -0.0900    0.0200   -0.7200 C   0  0  0  0  0  0  0  0  0  0  0  0
-       -1.0600   -0.6800   -0.9200 O   0  0  0  0  0  0  0  0  0  0  0  0
-        1.2000   -0.6200    1.3000 C   0  0  0  0  0  0  0  0  0  0  0  0
-        0.6600    0.4400   -1.7400 O   0  0  0  0  0  0  0  0  0  0  0  0
-       -1.3800   -0.4200    1.4800 H   0  0  0  0  0  0  0  0  0  0  0  0
-       -0.6800    0.6600    2.4500 H   0  0  0  0  0  0  0  0  0  0  0  0
-        0.7500    1.3900    0.6800 H   0  0  0  0  0  0  0  0  0  0  0  0
-        1.4600   -0.3300    2.3200 H   0  0  0  0  0  0  0  0  0  0  0  0
-        0.7200   -1.5900    1.3100 H   0  0  0  0  0  0  0  0  0  0  0  0
-        2.1100   -0.6800    0.7000 H   0  0  0  0  0  0  0  0  0  0  0  0
-        0.4400    0.1800   -2.6500 H   0  0  0  0  0  0  0  0  0  0  0  0
+       -0.9660    0.4930    1.5000 N   0  0  0  0  0  0  0  0  0  0  0  0
+        0.2570    0.4180    0.6920 C   0  0  0  0  0  0  0  0  0  0  0  0
+       -0.0940    0.0170   -0.7160 C   0  0  0  0  0  0  0  0  0  0  0  0
+       -1.0560   -0.6820   -0.9230 O   0  0  0  0  0  0  0  0  0  0  0  0
+        1.2040   -0.6200    1.2960 C   0  0  0  0  0  0  0  0  0  0  0  0
+        0.6610    0.4390   -1.7420 O   0  0  0  0  0  0  0  0  0  0  0  0
+       -1.3830   -0.4250    1.4820 H   0  0  0  0  0  0  0  0  0  0  0  0
+       -0.6760    0.6610    2.4520 H   0  0  0  0  0  0  0  0  0  0  0  0
+        0.7460    1.3920    0.6820 H   0  0  0  0  0  0  0  0  0  0  0  0
+        1.4590   -0.3300    2.3160 H   0  0  0  0  0  0  0  0  0  0  0  0
+        0.7150   -1.5940    1.3070 H   0  0  0  0  0  0  0  0  0  0  0  0
+        2.1130   -0.6760    0.6970 H   0  0  0  0  0  0  0  0  0  0  0  0
+        0.4350    0.1820   -2.6470 H   0  0  0  0  0  0  0  0  0  0  0  0
       1  2  1  0  0  0  0
       1  7  1  0  0  0  0
       1  8  1  0  0  0  0
@@ -731,12 +753,17 @@ class SDFile(File, MutableMapping):
     @staticmethod
     def deserialize(text):
         """
-        Create an object by deserializing the given text content.
+        Create an :class:`SDFile` by deserializing the given text content.
 
         Parameters
         ----------
-        content : str
+        text : str
             The content to be deserialized.
+
+        Returns
+        -------
+        file_object : SDFile
+            The parsed file.
         """
         lines = text.splitlines()
         record_ends = np.array(
@@ -895,7 +922,7 @@ def _to_metadata_key(key):
         return Metadata.Key(name=key)
     else:
         raise TypeError(
-            "Expected 'Metadata.Key' or str, " f"but got '{type(key).__name__}'"
+            f"Expected 'Metadata.Key' or str, but got '{type(key).__name__}'"
         )
 
 

biotite/structure/io/pdb/convert.py

@@ -89,14 +89,13 @@ def get_structure(
     -------
     array : AtomArray or AtomArrayStack
         The return type depends on the `model` parameter.
-
     """
     return pdb_file.get_structure(model, altloc, extra_fields, include_bonds)
 
 
 def set_structure(pdb_file, array, hybrid36=False):
     """
-    write an :class:`AtomArray` or :class:`AtomArrayStack` into a
+    Write an :class:`AtomArray` or :class:`AtomArrayStack` into a
     :class:`PDBFile`.
 
     This function is a thin wrapper around the :class:`PDBFile` method
@@ -114,7 +113,7 @@ def set_structure(pdb_file, array, hybrid36=False):
     array : AtomArray or AtomArrayStack
         The structure to be written. If a stack is given, each array in
         the stack will be in a separate model.
-    hybrid36: boolean, optional
+    hybrid36 : boolean, optional
         Defines wether the file should be written in hybrid-36 format.
 
     Notes
@@ -218,6 +217,8 @@ def get_assembly(
     assembly : AtomArray or AtomArrayStack
         The assembly.
         The return type depends on the `model` parameter.
+        Contains the `sym_id` annotation, which enumerates the copies of the asymmetric
+        unit in the assembly.
 
     Examples
     --------