The SARS-CoV-2 spike protein contains a fatty acid binding site, also found in some other coronaviruses (e.g. SARS-CoV), which binds linoleic acid and is functionally important. When occupied by linoleic acid, it reduces infectivity, by ‘locking’ the spike in a less infectious conformation. Here, we use dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to compare the response of spike variants to linoleic acid removal. These simulations show that the fatty acid site is coupled to functional regions of the protein, some of them far from the site (e.g. in the receptor-binding motif, N-terminal domain, the furin cleavage site located in position 679-685 and the fusion peptide-surrounding regions) and identify the allosteric networks involved in these connections. Comparison of the response of the original (‘Wuhan’) spike with four variants: Alpha, Delta, Delta plus and Omicron BA.1 show that the variants differ significantly in their response to linoleic acid removal. The allosteric connections to the fatty acid site on Alpha are generally similar to the original protein, except for the receptor-binding motif and S71-R78 region which show a weaker link to the FA site. In contrast, Omicron is the most affected variant exhibiting significant differences in the receptor-binding motif, N-terminal domain, V622-L629 and the furin cleavage site. These differences in allosteric modulation may be of functional relevance, e.g. in differences in transmissibility and virulence. Experimental comparison of the effects of linoleic acid on different variants is warranted.