being the pair interaction energy (Table ) and r ij is the distance between the beads given by the input structure. OLIVES then constructs a bond network by performing repeated rounds of maximum weight matching 28 (using NetworkX 29 v2.3) pairing hydrogen bond partners such that w ij is maximized. In the secondary network, the backbone beads are allowed to form two hydrogen bonds each, which is obtained by first performing a maximum weight matching and then removing matched pairs from the graph, followed by a second round of matching. The pattern resulting from the maximum weight matching is illustrated in Figure ,B1-4. In cases where the graph contains an odd number of vertices, some partners are not matched, which is solved by a third round on the small remaining graph. We found this final step to be closing some small gaps in the network important for the formation of bond patterns in helices. Note that maximum weight matching can be performed by forcing maximum cardinality 28 i.e. maximizing the number of edges, which we have tested but led to generally worse results in terms of protein stability. In the tertiary network, the graph does not distinguish between side chain -backbone and side chain -side chain pairs.