Bonds are currently partly supported, only disulfide bonds are read from PDBs. The main issue is in the implementation design. Currently the bonds are saved as a list of bond type and the serial numbers of the atoms involved on the PDB struct. This removes any possibility to trace the bound state from an individual atom. The implementation is open for discussion.

Best case from a user perspective

• From an atom find its bound state (amount of bonds + type of covalent bonds (single, double, triple))
• From an atom find the bounding partners maybe something like: Atom::bonds() -> impl Iterator<Item=(Bond, &Atom)>
• Read all bond types from PDB and mmCIF
• Fill in missing bonds for common structures, like all amino acids

Implementation idea: the bonds function on atoms could use a reference to the PDB to get the reference to the other atoms.