Comparison of structural networks across homologous proteins.

Protein sequence determines its structure and function. The indirect relationship between protein function and structure lies deep-rooted in the structural topology that has evolved into performing optimal function. The evolution of structure and its interconnectivity has been conventionally studied by comparing the root means square deviation between protein structures at the backbone level. Two factors that are necessary for the quantitative comparison of non-covalent interactions are (a) explicit inclusion of the coordinates of side-chain atoms and (b) consideration of multiple structures from the conformational landscape to account for structural variability. We have recently addressed these fundamental issues by investigating the alteration of inter-residue interactions across an ensemble of protein structure networks through a graph spectral approach. In this study, we have developed a rigorous method to compare the structure networks of homologous proteins, with a wide range of sequence identity percentages. A range of dissimilarity measures that show the extent of change in the network across homologous structures are generated, which also includes the comparison of the protein structure variability. We discuss in detail, scenarios where the variation of structure is not accompanied by loss or gain of the overall network and its vice versa. The sequence-based phylogeny among the homologs is also compared with the lineage obtained from information from such a robust structure comparison. In summary, we can obtain a quantitative comparison score for the structure networks of homologous proteins, which also enables us to study the evolution of protein function based on the variation of their topologies.

Belinostat, at Its Clinically Relevant Concentrations, Inhibits Rifampicin-Induced CYP3A4 and MDR1 Gene Expression.

Activation of human pregnane X receptor (hPXR) has been associated with induction of chemoresistance. It has been proposed that such chemoresistance via cytochrome P450/drug transporters can be reversed with the use of antagonists that specifically abrogate agonist-mediated hPXR activation. Unfortunately, proposed antagonists lack the specificity and appropriate pharmacological characteristics that allow these features to be active in the clinic. We propose that, ideally, an hPXR antagonist would be a cancer drug itself that is part of a "cancer drug cocktail" and effective as an hPXR antagonist at therapeutic concentrations. Belinostat (BEL), a histone deacetylase inhibitor approved for the treatment of relapsed/refractory peripheral T-cell lymphoma, and often used in combination with chemotherapy, is an attractive candidate based on its hPXR ligand-like features. We sought to determine whether these features of BEL might allow it to behave as an antagonist in combination chemotherapy regimens that include hPXR activators. BEL represses agonist-activated hPXR target gene expression at its therapeutic concentrations in human primary hepatocytes and LS174T human colon cancer cells. BEL repressed rifampicin-induced gene expression of CYP3A4 and multidrug resistance protein 1, as well as their respective protein activities. BEL decreased rifampicin-induced resistance to SN-38, the active metabolite of irinotecan, in LS174T cells. This finding indicates that BEL could suppress hPXR agonist-induced chemoresistance. BEL attenuated the agonist-induced steroid receptor coactivator-1 interaction with hPXR, and, together with molecular docking studies, the study suggests that BEL directly interacts with multiple sites on hPXR. Taken together, our results suggest that BEL, at its clinically relevant therapeutic concentration, can antagonize hPXR agonist-induced gene expression and chemoresistance.

From Quantum Chemistry to Networks in Biology: A Graph Spectral Approach to Protein Structure Analyses.

In this perspective article, we present a multidisciplinary approach for characterizing protein structure networks. We first place our approach in its historical context and describe the manner in which it synthesizes concepts from quantum chemistry, biology of polymer conformations, matrix mathematics, and percolation theory. We then explicitly provide the method for constructing the protein structure network in terms of noncovalently interacting amino acid side chains and show how a mine of information can be obtained from the graph spectra of these networks. Employing suitable mathematical approaches, such as the use of a weighted, Laplacian matrix to generate the spectra, enables us to develop rigorous methods for network comparison and to identify crucial nodes responsible for the network integrity through a perturbation approach. Our scoring methods have several applications in structural biology that are elusive to conventional methods of analyses. Here, we discuss the instances of (a) protein structure comparison that includes the details of side chain connectivity, (b) contribution to node clustering as a function of bound ligand, explaining the global effect of local changes in phenomena such as allostery, and (c) identification of crucial amino acids for structural integrity, derived purely from the spectra of the graph. We demonstrate how our method enables us to obtain valuable information on key proteins involved in cellular functions and diseases such as GPCR and HIV protease and discuss the biological implications. We then briefly describe how concepts from percolation theory further augment our analyses. In our concluding perspective for future developments, we suggest a further unifying approach to protein structure analyses and a judicious choice of questions to employ our methods for larger, more complex networks, such as metabolic and disease networks.

Validation of protein structure models using network similarity score.

Accurate structural validation of proteins is of extreme importance in studies like protein structure prediction, analysis of molecular dynamic simulation trajectories and finding subtle changes in very similar structures. The benchmarks for today's structure validation are scoring methods like global distance test-total structure (GDT-TS), TM-score and root mean square deviations (RMSD). However, there is a lack of methods that look at both the protein backbone and side-chain structures at the global connectivity level and provide information about the differences in connectivity. To address this gap, a graph spectral based method (NSS-network similarity score) which has been recently developed to rigorously compare networks in diverse fields, is adopted to compare protein structures both at the backbone and at the side-chain noncovalent connectivity levels. In this study, we validate the performance of NSS by investigating protein structures from X-ray structures, modeling (including CASP models), and molecular dynamics simulations. Further, we systematically identify the local and the global regions of the structures contributing to the difference in NSS, through the components of the score, a feature unique to this spectral based scoring scheme. It is demonstrated that the method can quantify subtle differences in connectivity compared to a reference protein structure and can form a robust basis for protein structure comparison. Additionally, we have also introduced a network-based method to analyze fluctuations in side chain interactions (edge-weights) in an ensemble of structures, which can be an useful tool for the analysis of MD trajectories.

Mechanism of Iron-Dependent Repressor (IdeR) Activation and DNA Binding: A Molecular Dynamics and Protein Structure Network Study.

Metalloproteins form a major class of enzymes in the living system that are involved in crucial biological functions such as catalysis, redox reactions and as 'switches' in signal transductions. Iron dependent repressor (IdeR) is a metal-sensing transcription factor that regulates free iron concentration in Mycobacterium tuberculosis. IdeR is also known to promote bacterial virulence, making it an important target in the field of therapeutics. Mechanistic details of how iron ions modulate IdeR such that it dimerizes and binds to DNA is not understood clearly. In this study, we have performed molecular dynamic simulations and integrated it with protein structure networks to study the influence of iron on IdeR structure and function. A significant structural variation between the metallated and the non-metallated system is observed. Our simulations clearly indicate the importance of iron in stabilizing its monomeric subunit, which in turn promotes dimerization. However, the most striking results are obtained from the simulations of IdeR-DNA complex in the absence of metals, where at the end of 100ns simulations, the protein subunits are seen to rapidly dissociate away from the DNA, thereby forming an excellent resource to investigate the mechanism of DNA binding. We have also investigated the role of iron as an allosteric regulator of IdeR that positively induces IdeR-DNA complex formation. Based on this study, a mechanistic model of IdeR activation and DNA binding has been proposed.

A histidine aspartate ionic lock gates the iron passage in miniferritins from Mycobacterium smegmatis.

Dps (DNA-binding protein from starved cells) are dodecameric assemblies belonging to the ferritin family that can bind DNA, carry out ferroxidation, and store iron in their shells. The ferritin-like trimeric pore harbors the channel for the entry and exit of iron. By representing the structure of Dps as a network we have identified a charge-driven interface formed by a histidine aspartate cluster at the pore interface unique to Mycobacterium smegmatis Dps protein, MsDps2. Site-directed mutagenesis was employed to generate mutants to disrupt the charged interactions. Kinetics of iron uptake/release of the wild type and mutants were compared. Crystal structures were solved at a resolution of 1.8-2.2 Å for the various mutants to compare structural alterations vis à vis the wild type protein. The substitutions at the pore interface resulted in alterations in the side chain conformations leading to an overall weakening of the interface network, especially in cases of substitutions that alter the charge at the pore interface. Contrary to earlier findings where conserved aspartate residues were found crucial for iron release, we propose here that in the case of MsDps2, it is the interplay of negative-positive potentials at the pore that enables proper functioning of the protein. In similar studies in ferritins, negative and positive patches near the iron exit pore were found to be important in iron uptake/release kinetics. The unique ionic cluster in MsDps2 makes it a suitable candidate to act as nano-delivery vehicle, as these gated pores can be manipulated to exhibit conformations allowing for slow or fast rates of iron release.

Distinctive contributions of the ribosomal P-site elements m2G966, m5C967 and the C-terminal tail of the S9 protein in the fidelity of initiation of translation in Escherichia coli.

The accuracy of pairing of the anticodon of the initiator tRNA (tRNA(fMet)) and the initiation codon of an mRNA, in the ribosomal P-site, is crucial for determining the translational reading frame. However, a direct role of any ribosomal element(s) in scrutinizing this pairing is unknown. The P-site elements, m(2)G966 (methylated by RsmD), m(5)C967 (methylated by RsmB) and the C-terminal tail of the protein S9 lie in the vicinity of tRNA(fMet). We investigated the role of these elements in initiation from various codons, namely, AUG, GUG, UUG, CUG, AUA, AUU, AUC and ACG with tRNA(fMet(CAU) (tRNA(fMet) with CAU anticodon); CAC and CAU with tRNA(fMet(GUG); UAG with tRNA(fMet(CAU) ; UAC with tRNA(fMet(GUG) ; and AUC with tRNA(fMet(GUG) using in vivo and computational methods. Although RsmB deficiency did not impact initiation from most codons, RsmD deficiency increased initiation from AUA, CAC and CAU (2- to 3.6-fold). Deletion of the S9 C-terminal tail resulted in poorer initiation from UUG, GUG and CUG, but in increased initiation from CAC, CAU and UAC codons (up to 4-fold). Also, the S9 tail suppressed initiation with tRNA(fMet(CAU) lacking the 3GC base pairs in the anticodon stem. These observations suggest distinctive roles of 966/967 methylations and the S9 tail in initiation.

GraSp-PSN: A web server for graph spectra based analysis of protein structure networks.

The function of a protein is most of the time achieved due to minute conformational changes in its structure due to ligand binding or environmental changes or other interactions. Hence the analysis of structure of proteins should go beyond the analysis of mere atom contacts and should include the emergent global structure as a whole. This can be achieved by graph spectra based analysis of protein structure networks. GraSp-PSN is a web server that can assist in (1) acquiring weighted protein structure network (PSN) and network parameters ranging from atomic level to global connectivity from the three dimensional coordinates of a protein, (2) generating scores for comparison of a pair of protein structures with detailed information of local to global connectivity, and (3) assigning perturbation scores to the residues and their interactions, that can prioritise them in terms of residue clusters. The methods implemented in the server are generic in nature and can be used for comparing networks in any discipline by uploading adjacency matrices in the server. The webserver can be accessed using the following link: https://pople.mbu.iisc.ac.in/.

Insights into the fold organization of TIM barrel from interaction energy based structure networks.

There are many well-known examples of proteins with low sequence similarity, adopting the same structural fold. This aspect of sequence-structure relationship has been extensively studied both experimentally and theoretically, however with limited success. Most of the studies consider remote homology or "sequence conservation" as the basis for their understanding. Recently "interaction energy" based network formalism (Protein Energy Networks (PENs)) was developed to understand the determinants of protein structures. In this paper we have used these PENs to investigate the common non-covalent interactions and their collective features which stabilize the TIM barrel fold. We have also developed a method of aligning PENs in order to understand the spatial conservation of interactions in the fold. We have identified key common interactions responsible for the conservation of the TIM fold, despite high sequence dissimilarity. For instance, the central beta barrel of the TIM fold is stabilized by long-range high energy electrostatic interactions and low-energy contiguous vdW interactions in certain families. The other interfaces like the helix-sheet or the helix-helix seem to be devoid of any high energy conserved interactions. Conserved interactions in the loop regions around the catalytic site of the TIM fold have also been identified, pointing out their significance in both structural and functional evolution. Based on these investigations, we have developed a novel network based phylogenetic analysis for remote homologues, which can perform better than sequence based phylogeny. Such an analysis is more meaningful from both structural and functional evolutionary perspective. We believe that the information obtained through the "interaction conservation" viewpoint and the subsequently developed method of structure network alignment, can shed new light in the fields of fold organization and de novo computational protein design.

Allosteric communication in cysteinyl tRNA synthetase: a network of direct and indirect readout.

Protein structure networks are constructed for the identification of long-range signaling pathways in cysteinyl tRNA synthetase (CysRS). Molecular dynamics simulation trajectory of CysRS-ligand complexes were used to determine conformational ensembles in order to gain insight into the allosteric signaling paths. Communication paths between the anticodon binding region and the aminoacylation region have been identified. Extensive interaction between the helix bundle domain and the anticodon binding domain, resulting in structural rigidity in the presence of tRNA, has been detected. Based on the predicted model, six residues along the communication paths have been examined by mutations (single and double) and shown to mediate a coordinated coupling between anticodon recognition and activation of amino acid at the active site. This study on CysRS clearly shows that specific key residues, which are involved in communication between distal sites in allosteric proteins but may be elusive in direct structure analysis, can be identified from dynamics of protein structure networks.

Understanding structural variability in proteins using protein structural networks.

Proteins perform their function by accessing a suitable conformer from the ensemble of available conformations. The conformational diversity of a chosen protein structure can be obtained by experimental methods under different conditions. A key issue is the accurate comparison of different conformations. A gold standard used for such a comparison is the root mean square deviation (RMSD) between the two structures. While extensive refinements of RMSD evaluation at the backbone level are available, a comprehensive framework including the side chain interaction is not well understood. Here we employ protein structure network (PSN) formalism, with the non-covalent interactions of side chain, explicitly treated. The PSNs thus constructed are compared through graph spectral method, which provides a comparison at the local and at the global structural level. In this work, PSNs of multiple crystal conformers of single-chain, single-domain proteins, are subject to pair-wise analysis to examine the dissimilarity in their network topologies and in order to determine the conformational diversity of their native structures. This information is utilized to classify the structural domains of proteins into different categories. It is observed that proteins typically tend to retain structure and interactions at the backbone level. However, some of them also depict variability in either their overall structure or only in their inter-residue connectivity at the sidechain level, or both. Variability of sub-networks based on solvent accessibility and secondary structure is studied. The types of specific interactions are found to contribute differently to structure variability. An ensemble analysis by computing the mathematical variance of edge-weights across multiple conformers provided information on the contribution to overall variability from each edge of the PSN. Interactions that are highly variable are identified and their impact on structure variability has been discussed with the help of a case study. The classification based on the present side-chain network-based studies provides a framework to correlate the structure-function relationships in protein structures.

Gefitinib Inhibits Rifampicin-Induced CYP3A4 Gene Expression in Human Hepatocytes.

During multidrug combination chemotherapy, activation of the nuclear receptor and the transcription factor human pregnane xenobiotic receptor (hPXR) has been shown to play a role in the development of chemoresistance. Mechanistically, this could occur due to the cancer drug activation of hPXR and the subsequent upregulation of hPXR target genes such as the drug metabolism enzyme, cytochrome P450 3A4 (CYP3A4). In the context of hPXR-mediated drug resistance, hPXR antagonists would be useful adjuncts to PXR-activating chemotherapy. However, there are currently no clinically approved hPXR antagonists in the market. Gefitinib (GEF), a tyrosine kinase inhibitor used for the treatment of advanced non-small-cell lung cancer and effectively used in combinational chemotherapy treatments, is a promising candidate owing to its hPXR ligand-like features. We, therefore, investigated whether GEF would act as an hPXR antagonist when combined with a known hPXR agonist, rifampicin (RIF). At therapeutically relevant concentrations, GEF successfully inhibited the RIF-induced upregulation of endogenous CYP3A4 gene expression in human primary hepatocytes and human hepatocells. Additionally, GEF inhibited the RIF induction of hPXR-mediated CYP3A4 promoter activity in HepG2 human liver carcinoma cells. The computational modeling of molecular docking predicted that GEF could bind to multiple sites on hPXR including the ligand-binding pocket, allowing for potential as a direct antagonist as well as an allosteric inhibitor. Indeed, GEF bound to the ligand-binding domain of the hPXR in cell-free assays, suggesting that GEF directly interacts with the hPXR. Taken together, our results suggest that GEF, at its clinically relevant therapeutic concentration, can antagonize the hPXR agonist-induced CYP3A4 gene expression in human hepatocytes. Thus, GEF could be a potential candidate for use in combinational chemotherapies to combat hPXR agonist-induced chemoresistance. Further studies are warranted to determine whether GEF has sufficient hPXR inhibitor abilities to overcome the hPXR agonist-induced chemoresistance.

Probing the allosteric mechanism in pyrrolysyl-tRNA synthetase using energy-weighted network formalism.

Pyrrolysyl-tRNA synthetase (PylRS) is an atypical enzyme responsible for charging tRNA(Pyl) with pyrrolysine, despite lacking precise tRNA anticodon recognition. This dimeric protein exhibits allosteric regulation of function, like any other tRNA synthetases. In this study we examine the paths of allosteric communication at the atomic level, through energy-weighted networks of Desulfitobacterium hafniense PylRS (DhPylRS) and its complexes with tRNA(Pyl) and activated pyrrolysine. We performed molecular dynamics simulations of the structures of these complexes to obtain an ensemble conformation-population perspective. Weighted graph parameters relevant to identifying key players and ties in the context of social networks such as edge/node betweenness, closeness index, and the concept of funneling are explored in identifying key residues and interactions leading to shortest paths of communication in the structure networks of DhPylRS. Further, the changes in the status of important residues and connections and the costs of communication due to ligand induced perturbations are evaluated. The optimal, suboptimal, and preexisting paths are also investigated. Many of these parameters have exhibited an enhanced asymmetry between the two subunits of the dimeric protein, especially in the pretransfer complex, leading us to conclude that encoding of function goes beyond the sequence/structure of proteins. The local and global perturbations mediated by appropriate ligands and their influence on the equilibrium ensemble of conformations also have a significant role to play in the functioning of proteins. Taking a comprehensive view of these observations, we propose that the origin of many functional aspects (allostery and half-sites reactivity in the case of DhPylRS) lies in subtle rearrangements of interactions and dynamics at a global level.

Network approach for capturing ligand-induced subtle global changes in protein structures.

Ligand-induced conformational changes in proteins are of immense functional relevance. It is a major challenge to elucidate the network of amino acids that are responsible for the percolation of ligand-induced conformational changes to distal regions in the protein from a global perspective. Functionally important subtle conformational changes (at the level of side-chain noncovalent interactions) upon ligand binding or as a result of environmental variations are also elusive in conventional studies such as those using root-mean-square deviations (r.m.s.d.s). In this article, the network representation of protein structures and their analyses provides an efficient tool to capture these variations (both drastic and subtle) in atomistic detail in a global milieu. A generalized graph theoretical metric, using network parameters such as cliques and/or communities, is used to determine similarities or differences between structures in a rigorous manner. The ligand-induced global rewiring in the protein structures is also quantified in terms of network parameters. Thus, a judicious use of graph theory in the context of protein structures can provide meaningful insights into global structural reorganizations upon perturbation and can also be helpful for rigorous structural comparison. Data sets for the present study include high-resolution crystal structures of serine proteases from the S1A family and are probed to quantify the ligand-induced subtle structural variations.

Network Re-Wiring During Allostery and Protein-Protein Interactions: A Graph Spectral Approach.

The process of allostery is often guided by subtle changes in the non-covalent interactions between residues of a protein. These changes may be brought about by minor perturbations by natural processes like binding of a ligand or protein-protein interaction. The challenge lies in capturing minute changes at the residue interaction level and following their propagation at local as well as global distances. While macromolecular effects of the phenomenon of allostery are inferred from experiments, a computational microscope can elucidate atomistic-level details leading to such macromolecular effects. Network formalism has served as an attractive means to follow this path and has been pursued further for the past couple of decades. In this chapter some concepts and methods are summarized, and recent advances are discussed. Specifically, the changes in strength of interactions (edge weight) and their repercussion on the overall protein organization (residue clustering) are highlighted. In this review, we adopt a graph spectral method to probe these subtle changes in a quantitative manner. Further, the power of this method is demonstrated for capturing re-ordering of side-chain interactions in response to ligand binding, which culminates into formation of a protein-protein complex in ß2-adrenergic receptors.

Interaction energy based protein structure networks.

The three-dimensional structure of a protein is formed and maintained by the noncovalent interactions among the amino-acid residues of the polypeptide chain. These interactions can be represented collectively in the form of a network. So far, such networks have been investigated by considering the connections based on distances between the amino-acid residues. Here we present a method of constructing the structure network based on interaction energies among the amino-acid residues in the protein. We have investigated the properties of such protein energy-based networks (PENs) and have shown correlations to protein structural features such as the clusters of residues involved in stability, formation of secondary and super-secondary structural units. Further we demonstrate that the analysis of PENs in terms of parameters such as hubs and shortest paths can provide a variety of biologically important information, such as the residues crucial for stabilizing the folded units and the paths of communication between distal residues in the protein. Finally, the energy regimes for different levels of stabilization in the protein structure have clearly emerged from the PEN analysis.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks.

BACKGROUND: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. METHODS: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. RESULTS: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. CONCLUSION: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Allostery and conformational free energy changes in human tryptophanyl-tRNA synthetase from essential dynamics and structure networks.

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of allostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side-chains and then to compute the ligand-induced population shift. Finally, we obtain the free-energy landscape of the protein in equilibrium, characterizing the free-energy minima accessed by the protein complexes. We have chosen human tryptophanyl-tRNA synthetase (hTrpRS), a protein responsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure.

Elucidation of the conformational free energy landscape in H.pylori LuxS and its implications to catalysis.

BACKGROUND: One of the major challenges in understanding enzyme catalysis is to identify the different conformations and their populations at detailed molecular level in response to ligand binding/environment. A detail description of the ligand induced conformational changes provides meaningful insights into the mechanism of action of enzymes and thus its function. RESULTS: In this study, we have explored the ligand induced conformational changes in H.pylori LuxS and the associated mechanistic features. LuxS, a dimeric protein, produces the precursor (4,5-dihydroxy-2,3-pentanedione) for autoinducer-2 production which is a signalling molecule for bacterial quorum sensing. We have performed molecular dynamics simulations on H.pylori LuxS in its various ligand bound forms and analyzed the simulation trajectories using various techniques including the structure network analysis, free energy evaluation and water dynamics at the active site. The results bring out the mechanistic details such as co-operativity and asymmetry between the two subunits, subtle changes in the conformation as a response to the binding of active and inactive forms of ligands and the population distribution of different conformations in equilibrium. These investigations have enabled us to probe the free energy landscape and identify the corresponding conformations in terms of network parameters. In addition, we have also elucidated the variations in the dynamics of water co-ordination to the Zn2+ ion in LuxS and its relation to the rigidity at the active sites. CONCLUSIONS: In this article, we provide details of a novel method for the identification of conformational changes in the different ligand bound states of the protein, evaluation of ligand-induced free energy changes and the biological relevance of our results in the context of LuxS structure-function. The methodology outlined here is highly generalized to illuminate the linkage between structure and function in any protein of known structure.

Surveying the Side-Chain Network Approach to Protein Structure and Dynamics: The SARS-CoV-2 Spike Protein as an Illustrative Case.

Network theory-based approaches provide valuable insights into the variations in global structural connectivity between different dynamical states of proteins. Our objective is to review network-based analyses to elucidate such variations, especially in the context of subtle conformational changes. We present technical details of the construction and analyses of protein structure networks, encompassing both the non-covalent connectivity and dynamics. We examine the selection of optimal criteria for connectivity based on the physical concept of percolation. We highlight the advantages of using side-chain-based network metrics in contrast to backbone measurements. As an illustrative example, we apply the described network approach to investigate the global conformational changes between the closed and partially open states of the SARS-CoV-2 spike protein. These conformational changes in the spike protein is crucial for coronavirus entry and fusion into human cells. Our analysis reveals global structural reorientations between the two states of the spike protein despite small changes between the two states at the backbone level. We also observe some differences at strategic locations in the structures, correlating with their functions, asserting the advantages of the side-chain network analysis. Finally, we present a view of allostery as a subtle synergistic-global change between the ligand and the receptor, the incorporation of which would enhance drug design strategies.