An exploration of the SARS-CoV-2 spike receptor binding domain (RBD), a complex palette of evolutionary and structural features

SARS-CoV-2 spike protein (S) is associated with the entry of virus inside the host cell by recruiting its loop dominant receptor binding domain (RBD) and interacting with the host ACE2 receptor. Our study deploying a two-tier approach encompassing evolutionary and structural analysis provides a comprehensive picture of the RBD, which could be of potential use for better understanding the RBD and address its druggability issues. Resorting to an ensemble of sequence space exploratory tools including co-evolutionary analysis and deep mutational scans we provide a quantitative insight into the evolutionarily constrained subspace of the RBD sequence space. Guided by structure network analysis and Monte Carlo simulation we highlight regions inside the RBD, which are critical for providing structural integrity and conformational flexibility of the binding cleft. We further deployed fuzzy C-means clustering by plugging the evolutionary and structural features of discrete structure blocks of RBD to understand which structure blocks share maximum overlap based on their evolutionary and structural features. Deploying this multi-tier interlinked approach, which essentially distilled the evolutionary and structural features of RBD, we highlight discrete region, which could be a potential druggable pocket thereby destabilizing the structure and addressing evolutionary routes.