Checking halogen bonds

Halogen bonds (or X-bonds) are one of the less appreciated protein-ligand interactions. As we discussed in 2022, the polarized nature of a carbon-halogen bond creates a partially-positively charged “σ-hold” at the bit of the halogen furthest from the carbon, and this can make favorable interactions with lone pairs on oxygen or sulfur atoms (or nitrogen, but in most proteins this is limited to histidine residues and is rare.) Halogens can also interact with aromatic π-systems such as the side chains of phenylalanine, tryptophan, histidine, and tryptophan. Since many fragments contain halogen atoms by design, halogen bonds may occur frequently. But how do you decide whether “a halogen in proximity of a possible acceptor” actually contributes to binding? In a new (open-access) paper in Protein Science, Ida de Vries, Robbie Joosten, and colleagues at Oncode Institute and The Netherlands Cancer Institute provide a new metric.

 

The researchers examined structures of halogen-containing ligands bound to proteins in PDB-REDO, a database of carefully vetted and refined structures from the Protein Data Bank. They only included structures solved to better than 2.5 Å resolution and omitted structures where halogens had high B-factors, which may be the result of radiation damage. This led to 8423 structures in which a halogen possibly interacted with an oxygen or sulfur atom and 8096 potential halogen-π interactions, which were analyzed in detail.

 

A halogen bond to an oxygen or sulfur atom can be described by the interatom distance and two angles: θ1 (carbon-halogen-oxygen/sulfur) and θ2 (halogen-oxygen/sulfur-carbon). Halogen-π-system bonds can be defined by distance to the centroid of the π-system and θ1, the carbon-halogen-centroid angle. (The paper has a nice diagram.) These parameters were calculated and annotated for all the structures.

 

Median distances were 3.5 Å between halogen and oxygen/sulfur, regardless of the halogen. Median θ1 angles were smaller than the 150º-180º expected, particularly for fluorine atoms, while median θ2 angles were more consistent with theory, at 90º-120º.

 

For halogen-π-systems, median distances were 4.8 Å for all halogens except iodine, which came in slightly higher. But θ1 angles were still smaller than expected, mostly between 110º-140º.

 

Armed with this tranche of high-quality data, the researchers established a Halogen Bond Score, or HalBS. For any potential halogen bond in a new crystal structure or other structural model, the distance, θ1, and, if applicable, the θ2 values are calculated, and if any of these diverge too far from the median values, HalBS flags them. Importantly, the researchers acknowledge that “the current HalBS cannot be used as a direct validation metric but can provide an indication of genuine halogen bonds and ‘not so proper’ halogen bonds.”

 

With this caveat HalBS could be useful, and the researchers have made the source code available at https://github.com/PDB-REDO/HalBS (though the link doesn’t seem to work for me). As they note, more data, such as might be provided by widespread deposition of large crystallographic fragment screens, could further refine HalBS. Of course, the existence of a halogen bond exists says little about how much binding energy it contributes, but it’s a start.