Naganathan Athi N.

Identifying the series of molecular events that take place when a protein molecule folds or unfolds has confounded and delighted researchers alike while challenging the limits of experimental approaches and the applicability of theoretical models and simulations. This chapter discusses the vast conformational space accessible to a protein chain, the diversity of weak noncovalent interactions it makes, the role of solvent compensating energetic-entropic terms, and the cellular environment, all of which weave an intricate fabric of complexity during the folding of a protein. Emphasis is placed on how generating conformational landscapes in a quantitative fashion can provide an unparalleled view of competing substates, which can either aid in or hinder folding, thus playing a role in function and disease. The role and effect of mutations, the drivers of evolution, are discussed in detail along with the current high pedestal of disordered proteins that has overturned the conventional structure-function paradigm.

Aggregation of Cu–Zn superoxide dismutase (SOD1) is implicated in the motor neuron disease, amyotrophic lateral sclerosis (ALS). Although more than 140 disease mutations of SOD1 are available, their stability or aggregation behaviors in membrane environment are not correlated with disease pathophysiology. Here, we use multiple mutational variants of SOD1 to show that the absence of Zn, and not Cu, significantly impacts membrane attachment of SOD1 through two loop regions facilitating aggregation driven by lipid-induced conformational changes. These loop regions influence both the primary (through Cu intake) and the gain of function (through aggregation) of SOD1 presumably through a shared conformational landscape. Combining experimental and theoretical frameworks using representative ALS disease mutants, we develop a ‘co-factor derived membrane association model’ wherein mutational stress closer to the Zn (but not to the Cu) pocket is responsible for membrane association-mediated toxic aggregation and survival time scale after ALS diagnosis.

Phosphoglycerate kinase has been a model for the stability, folding cooperativity and catalysis of a two-domain protein. The human isoform 1 (hPGK1) is associated with cancer development and rare genetic diseases that affect several of its features. To investigate how mutations affect hPGK1 folding landscape and interaction networks, we have introduced mutations at a buried site in the N-terminal domain (F25 mutants) that either created cavities (F25L, F25V, F25A), enhanced conformational entropy (F25G) or introduced structural strain (F25W) and evaluated their effects using biophysical experimental and theoretical methods. All F25 mutants folded well, but showed reduced unfolding cooperativity, kinetic stability and altered activation energetics according to the results from thermal and chemical denaturation analyses. These alterations correlated well with the structural perturbation caused by mutations in the N-terminal domain and the destabilization caused in the interdomain interface as revealed by H/D exchange under native conditions. Importantly, experimental and theoretical analyses showed that these effects are significant even when the perturbation is mild and local. Our approach will be useful to establish the molecular basis of hPGK1 genotype–phenotype correlations due to phosphorylation events and single amino acid substitutions associated with disease.

The Publisher regrets that the “Conflict of Interest” statement was not included in the published article at the time of publication. The Authors confirm that they do not have any Conflict of interest to report for this article. The Publisher would like to apologise for any inconvenience caused.

Over 40% of eukaryotic proteomes and 15% of bacterial proteomes are predicted to be intrinsically disordered based on their amino acid sequence. Intrinsically disordered proteins (IDPs) exist as heterogeneous ensembles of interconverting conformations and pose a challenge to the structure–function paradigm by apparently functioning without possessing stable structural elements. IDPs play a prominent role in biological processes involving extensive intermolecular interaction networks and their inherently dynamic nature facilitates their promiscuous interaction with multiple structurally diverse partner molecules. NMR spectroscopy has made pivotal contributions to our understanding of IDPs because of its unique ability to characterize heterogeneity at atomic resolution. NMR methods such as Chemical Exchange Saturation Transfer (CEST) and relaxation dispersion have enabled the detection of ‘invisible’ excited states in biomolecules which are transiently and sparsely populated, yet central for function. Here, we develop a 1Hα CEST pulse sequence which overcomes the resonance overlap problem in the 1Hα-13Cα plane of IDPs by taking advantage of the superior resolution in the 1H-15N correlation spectrum. In this sequence, magnetization is transferred after 1H CEST using a triple resonance coherence transfer pathway from 1Hα (i) to 1HN(i + 1) during which the 15N(t1) and 1HN(t2) are frequency labelled. This approach is integrated with spin state-selective CEST for eliminating spurious dips in CEST profiles resulting from dipolar cross-relaxation. We apply this sequence to determine the excited state 1Hα chemical shifts of the intrinsically disordered DNA binding domain (CytRN) of the bacterial cytidine repressor (CytR), which transiently acquires a functional globally folded conformation. The structure of the excited state, calculated using 1Hα chemical shifts in conjunction with other excited state NMR restraints, is a three-helix bundle incorporating a helix-turn-helix motif that is vital for binding DNA.

The functioning of proteins is intimately tied to their fluctuations in the native ensemble. The structural-energetic features that determine fluctuation amplitudes and hence the shape of the underlying landscape, which in turn determine the magnitude of the functional output, are often confounded by multiple variables. Here, we employ the FF1 domain from human p190A RhoGAP protein as a model system to uncover the molecular basis for phosphorylation of a buried tyrosine, which is crucial to the transcriptional activity associated with transcription factor TFII-I. Combining spectroscopy, calorimetry, statistical-mechanical modeling, molecular simulations, and in vitro phosphorylation assays, we show that the FF1 domain samples a diverse array of conformations in its native ensemble, some of which are phosphorylation-competent. Upon eliminating unfavorable charge-charge interactions through a single charge-reversal (K53E) or charge-neutralizing (K53Q) mutation, we observe proportionately lower phosphorylation extents due to the altered structural coupling, damped equilibrium fluctuations, and a more compact native ensemble. We thus establish a conformational selection mechanism for phosphorylation in the FF1 domain with K53 acting as a “gatekeeper”, modulating the solvent exposure of the buried tyrosine. Our work demonstrates the role of unfavorable charge-charge interactions in governing functional events through the modulation of native ensemble characteristics, a feature that could be prevalent in ordered protein domains.

The strength of intraprotein interactions or contact network is one of the dominant factors determining the thermodynamic stabilities of proteins. The nature and the extent of connectivity of this network also play a role in allosteric signal propagation characteristics upon ligand binding to a protein domain. Here, we develop a server for rapid quantification of the strength of an interaction network by employing an experimentally consistent perturbation approach previously validated against a large data set of 375 mutations in 19 different proteins. The web server can be employed to predict the extent of destabilization of proteins arising from mutations in the protein interior in experimentally relevant units. Moreover, coupling distances - a measure of the extent of percolation on perturbation - and overall perturbation magnitudes are predicted in a residue-specific manner, enabling a first look at the distribution of energetic couplings in a protein or its changes upon ligand binding. We show specific examples of how the server can be employed to probe for the distribution of local stabilities in a protein, to examine changes in side chain orientations or packing before and after ligand binding, and to predict changes in stabilities of proteins upon mutations of buried residues. The web server is freely available at http://pbl.biotech.iitm.ac.in/pPerturb and supports recent versions of all major browsers.

A longstanding goal in the field of intrinsically disordered proteins (IDPs) is to characterize their structural heterogeneity and pinpoint the role of this heterogeneity in IDP function. Here, we use multinuclear chemical exchange saturation (CEST) nuclear magnetic resonance to determine the structure of a thermally accessible globally folded excited state in equilibrium with the intrinsically disordered native ensemble of a bacterial transcriptional regulator CytR. We further provide evidence from double resonance CEST experiments that the excited state, which structurally resembles the DNA-bound form of cytidine repressor (CytR), recognizes DNA by means of a “folding-before-binding” conformational selection pathway. The disorder-to-order regulatory switch in DNA recognition by natively disordered CytR therefore operates through a dynamical variant of the lock-and-key mechanism where the structurally complementary conformation is transiently accessed via thermal fluctuations.

Mutational effects in globular proteins exhibit an exponential-like decreasing dependence on distance from the mutated site, suggestive of long-range modulation of structural-thermodynamic features. Here, we extract the physical origins of this pattern by employing a statistical-mechanical model to construct conformational ensembles of three archetypal proteins. Through large-scale in silico alanine-scanning mutagenesis, we show that inter-residue differential coupling free energies, which are characteristic ensemble thermodynamic properties, follow a similar exponential distance dependence with the effects felt until ∼15–20 Å from the mutated site. From the perspective of an ensemble-averaged structure, this feature arises via long-range reorganization of the interaction network on mutations which is more significant for charged residues compared to hydrophobic residues. Our work highlights how subtle alterations in the microscopic distribution of states manifest as a macroscopic distance dependence, the physical origins of mutation-induced dynamic allostery, and the necessity to consider the global intra-protein interaction network to understand mutational outcomes.

The multifunctional nature of human flavoproteins is critically linked to their ability to populate multiple conformational states. Ligand binding, post-translational modifications and disease-associated mutations can reshape this functional landscape, although the structure-function relationships of these effects are not well understood. Herein, we characterized the structural and functional consequences of two mutations (the cancer-associated P187S and the phosphomimetic S82D) on different ligation states which are relevant to flavin binding, intracellular stability and catalysis of the disease-associated NQO1 flavoprotein. We found that these mutations affected the stability locally and their effects propagated differently through the protein structure depending both on the nature of the mutation and the ligand bound, showing directional preference from the mutated site and leading to specific phenotypic manifestations in different functional traits (FAD binding, catalysis and inhibition, intracellular stability and pharmacological response to ligands). Our study thus supports that pleitropic effects of disease-causing mutations and phosphorylation events on human flavoproteins may be caused by long-range structural propagation of stability effects to different functional sites that depend on the ligation-state and site-specific perturbations. Our approach can be of general application to investigate these pleiotropic effects at the flavoproteome scale in the absence of high-resolution structural models.