code: def show_metal(self):
          
          grp = self.getPseudoBondGroup("Metal Coordination-%i" % self.tid, associateWith=[self.model])
          grp.lineWidth = 3
          for i, metal in enumerate(self.plcomplex.metal_complexes):
              c = grp.newPseudoBond(self.atoms[metal.metal_id], self.atoms[metal.target_id])
              c.color = self.colorbyname('magenta')
  
              if metal.location == 'water':
                  self.water_ids.append(metal.target_id)
  
              if metal.location.startswith('protein'):
                  self.bs_res_ids.append(metal.target_id)
docs: Visualizes metal coordination.