Subsets of Slow Dynamic Modes Reveal Global Information Sources as Allosteric Sites

Exploring Druggable Binding Sites on the Class A GPCRs Using the Residue Interaction Network and Site Identification by Ligand Competitive Saturation

G protein-coupled receptors (GPCRs) play a central role in cellular signaling and are linked to many diseases. Accordingly, computational methods to explore potential allosteric sites for this class of proteins to facilitate the identification of potential modulators are needed. Importantly, the availability of rich structural data providing the locations of the orthosteric ligands and allosteric modulators targeting different GPCRs allows for the validation of approaches to identify new allosteric binding sites. Here, we validate the combination of two computational techniques, the residue interaction network (RIN) model and the site identification by ligand competitive saturation (SILCS) method, to predict putative allosteric binding sites of class A GPCRs. RIN analysis identifies hub residues that mediate allosteric signaling within a receptor and have a high capacity to alter receptor dynamics upon ligand binding. The known orthosteric (and allosteric) binding sites of 18 distinct class A GPCRs were successfully predicted by RIN through a dataset of 105 crystal structures (91 ligand-bound, 14 unbound) with up to 77.8% (76.9%) sensitivity, 92.5% (95.3%) specificity, 51.9% (50%) precision, and 86.2% (92.4%) accuracy based on the experimental and theoretical binding site data. Moreover, graph spectral analysis of the residue networks revealed that the proposed sites were located at the interfaces of highly interconnected residue clusters with a high ability to coordinate the functional dynamics. Then, we employed the SILCS-Hotspots method to assess the druggability of the novel sites predicted for 7 distinct class A GPCRs that are critical for a variety of diseases. While the known orthosteric and allosteric binding sites are successfully explored by our approach, numerous putative allosteric sites with the potential to bind drug-like molecules are proposed. The computational approach presented here promises to be a highly effective tool to predict putative allosteric sites of GPCRs to facilitate the design of effective modulators.

Computational analysis of long-range allosteric communications in CFTR

Malfunction of the CFTR protein results in cystic fibrosis, one of the most common hereditary diseases. CFTR functions as an anion channel, the gating of which is controlled by long-range allosteric communications. Allostery also has direct bearings on CF treatment: the most effective CFTR drugs modulate its activity allosterically. Herein, we integrated Gaussian network model, transfer entropy, and anisotropic normal mode-Langevin dynamics and investigated the allosteric communications network of CFTR. The results are in remarkable agreement with experimental observations and mutational analysis and provide extensive novel insight. We identified residues that serve as pivotal allosteric sources and transducers, many of which correspond to disease-causing mutations. We find that in the ATP-free form, dynamic fluctuations of the residues that comprise the ATP-binding sites facilitate the initial binding of the nucleotide. Subsequent binding of ATP then brings to the fore and focuses on dynamic fluctuations that were present in a latent and diffuse form in the absence of ATP. We demonstrate that drugs that potentiate CFTR's conductance do so not by directly acting on the gating residues, but rather by mimicking the allosteric signal sent by the ATP-binding sites. We have also uncovered a previously undiscovered allosteric 'hotspot' located proximal to the docking site of the phosphorylated regulatory (R) domain, thereby establishing a molecular foundation for its phosphorylation-dependent excitatory role. This study unveils the molecular underpinnings of allosteric connectivity within CFTR and highlights a novel allosteric 'hotspot' that could serve as a promising target for the development of novel therapeutic interventions.

Structural Basis of Sequential and Concerted Cooperativity

Allostery is a property of biological macromolecules featuring cooperative ligand binding and regulation of ligand affinity by effectors. The definition was introduced by Monod and Jacob in 1963, and formally developed as the "concerted model" by Monod, Wyman, and Changeux in 1965. Since its inception, this model of cooperativity was seen as distinct from and not reducible to the "sequential model" originally formulated by Pauling in 1935, which was developed further by Koshland, Nemethy, and Filmer in 1966. However, it is difficult to decide which model is more appropriate from equilibrium or kinetics measurements alone. In this paper, we examine several cooperative proteins whose functional behavior, whether sequential or concerted, is established, and offer a combined approach based on functional and structural analysis. We find that isologous, mostly helical interfaces are common in cooperative proteins regardless of their mechanism. On the other hand, the relative contribution of tertiary and quaternary structural changes, as well as the asymmetry in the liganded state, may help distinguish between the two mechanisms.