Integrative modeling of guanylate binding protein dimers

Guanylate binding proteins (GBPs) are interferon-γ-activated large GTPases, effective against intracellular pathogens like Toxoplasma gondii. Their host-protective functions require oligomerization, however, the oligomer structures have not been completely resolved yet. Here, we provide dimer models for hGBP1 and the murine GBPs 2 and 7 (mGBP2 and mGBP7) based on integrative modeling that involves the crystal structure of the G domain dimer of hGBP1, cross-linking mass spectrometry (XL-MS), small angle X-ray scattering (SAXS), protein-protein docking, and molecular dynamics (MD) simulations of hGBP1, mGBP2 and mGBP7. We first compare the sequences and protein dynamics of the monomeric hGBP1, mGBP2, and mGBP7, finding that the M/E domain of all three proteins is highly mobile featuring a hinge movement, yet this motion is less pronounced in mGBP7 while its GTPase (G) domain is more flexible. These differences can be explained by the variations in the sequences between mGBP7 and hGBP1/mGBP2 and extend to their dimers. While hGBP1 and its close orthologue mGBP2 dimerize via their G domains, mGBP7 shows a variety of possible dimer structures, among them parallel and crossed-stalk conformations. The G domain is only partly involved in mGBP7 dimerization, which provides a rational why mGBP7, unlike hGBP1 and mGBP2, can dimerize in the absence of GTP. The different GBP dimer structures, which still exhibit hinge movements to certain degrees, are expected to encode diverging functions, such as a destabilization of pathogenic membranes or fusion of the parasitophorous vacuole membrane with the autophagic machinery.