Molecular Determinants for Guanine Binding in GTP-Binding Proteins: A Data Mining and Quantum Chemical Study

GTP-binding proteins are essential molecular switches that regulate a wide range of cellular processes. Their function relies on the specific recognition and binding of guanine within their binding pockets. This study aims to elucidate the molecular determinants underlying this recognition. A large-scale data mining of the Protein Data Bank yielded 298 GTP-binding protein complexes, which provided a structural foundation for a systematic analysis of the intermolecular interactions that are responsible for the molecular recognition of guanine in proteins. It was found that multiple modes of non-bonded interactions including hydrogen bonding, cation-π interactions, and π-π stacking interactions are employed by GTP-binding proteins for binding. Subsequently, the strengths of non-bonded interaction energies between guanine and its surrounding protein residues were quantified by means of the double-hybrid DFT method B2PLYP-D3/cc-pVDZ. Hydrogen bonds, particularly those involving the N2 and O6 atoms of guanine, confer specificity to guanine recognition. Cation-π interactions between the guanine ring and basic residues (Lys and Arg) provide significant electrostatic stabilization. π-π stacking interactions with aromatic residues (Phe, Tyr, and Trp) further contribute to the overall binding affinity. This synergistic interplay of multiple interaction modes enables GTP-binding proteins to achieve high specificity and stability in guanine recognition, ultimately underpinning their crucial roles in cellular signaling and regulation. Notably, the NKXD motif, while historically considered crucial for guanine binding in GTP-binding proteins, is not universally required. Our study revealed significant variability in hydrogen bonding patterns, with many proteins lacking the NKXD motif but still effectively binding guanine through alternative arrangements of interacting residues.