Lipid modification of viral proteins with fatty acids of different lengths (S-acylation) is crucial for virus pathogenesis. The reaction is catalyzed by members of the DHHC family and proceeds in two steps: the autoacylation is followed by the acyl chain transfer onto protein substrates. The crystal structure of human DHHC20 (hDHHC20), an enzyme involved in the acylation of S-protein of SARS-CoV-2, revealed that the acyl chain may be inserted into a hydrophobic cavity formed by four transmembrane (TM) α-helices. To test this model, we used molecular dynamics of membrane-embedded hDHHC20 and its mutants either in the absence or presence of various acyl-CoAs. We found that among a range of acyl chain lengths probed only C16 adopts a conformation suitable for hDHHC20 autoacylation. This specificity is altered if the small or bulky residues at the cavity's ceiling are exchanged, e.g., the V185G mutant obtains strong preferences for binding C18. Surprisingly, an unusual hydrophilic ridge was found in TM helix 4 of hDHHC20, and the responsive hydrophilic patch supposedly involved in association was found in the 3D model of the S-protein TM-domain trimer. Finally, the exchange of critical Thr and Ser residues in the spike led to a significant decrease in its S-acylation. Our data allow further development of peptide/lipid-based inhibitors of hDHHC20 that might impede replication of Corona- and other enveloped viruses.
Keywords: DHHC-PATs; S-acylation; acyl CoA; molecular dynamics; palmitoylation.