Pawan Bhatta

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, continually undergoes mutation, leading to variants with altered pathogenicity and transmissibility. The Omicron variant (B.1.1.529), first identified in South Africa in 2021, has become the dominant strain worldwide. It harbors approximately 50 mutations compared to the original strain, with 15 located in the receptor-binding domain (RBD) of the spike protein that facilitates viral entry via binding to the human angiotensin-converting enzyme 2 (ACE2) receptor. How do these mutated residues modulate the intermolecular interactions and binding affinity between the RBD and ACE2? This is a question of great theoretical importance and practical implication. In this study, we employed quantum chemical calculations at the B2PLYP-D3/def2-TZVP level of theory to investigate the molecular determinants governing Omicron’s ACE2 interaction. Comparative analysis of the Omicron and wild-type RBD–ACE2 interfaces revealed that mutations including S477N, Q493R, Q498R, and N501Y enhance binding through the formation of bifurcated hydrogen bonds, π–π stacking, and cation–π interactions. These favorable interactions counterbalance such destabilizing mutations as K417N, G446S, G496S, and Y505H, which disrupt salt bridges and hydrogen bonds. Additionally, allosteric effects improve the contributions of non-mutated residues (notably A475, Y453, and F486) via structural realignment and novel hydrogen bonding with ACE2 residues such as S19, leading to an overall increase in the electrostatic and π-system interaction energy. In conclusion, our findings provide a mechanistic basis for Omicron’s increased infectivity and offer valuable insights for the development of targeted antiviral therapies. Moreover, from a methodological perspective, we directly calculated mutation-induced binding energy changes at the residue level using advanced quantum chemical methods rather than relying on the indirect decomposition schemes typical of molecular dynamics-based free energy analyses. The strong correlation between calculated energy differences and experimental deep mutational scanning (DMS) data underscores the robustness of the theoretical framework in predicting the effects of RBD mutations on ACE2 binding affinity. This demonstrates the potential of quantum chemical methods as predictive tools for studying mutation-induced changes in protein–protein interactions and guiding rational therapeutic design.

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.