code: def show_wbridges(self):
          
          grp = self.getPseudoBondGroup("Water Bridges-%i" % self.tid, associateWith=[self.model])
          grp.lineWidth = 3
          for i, wbridge in enumerate(self.plcomplex.waterbridges):
              c = grp.newPseudoBond(self.atoms[wbridge.water_id], self.atoms[wbridge.acc_id])
              c.color = self.colorbyname('cornflower blue')
              self.water_ids.append(wbridge.water_id)
              b = grp.newPseudoBond(self.atoms[wbridge.don_id], self.atoms[wbridge.water_id])
              b.color = self.colorbyname('cornflower blue')
              self.water_ids.append(wbridge.water_id)
              if wbridge.protisdon:
                  self.bs_res_ids.append(wbridge.don_id)
              else:
                  self.bs_res_ids.append(wbridge.acc_id)
docs: Visualizes water bridges