Deciphering the molecular organization of GET pathway chaperones through native mass spectrometry.

Get3/4/5 chaperone complex is responsible for targeting C-terminal tail-anchored membrane proteins to the endoplasmic reticulum. Despite the availability of several crystal structures of independent proteins and partial structures of subcomplexes, different models of oligomeric states and structural organization have been proposed for the protein complexes involved. Here, using native mass spectrometry (Native-MS), coupled with intact dissociation, we show that Get4/5 exclusively forms a tetramer using both Get5/5 and a novel Get4/4 dimerization interface. Addition of Get3 to this leads to a hexameric (Get3)2-(Get4)2-(Get5)2 complex with closed-ring cyclic architecture. We further validate our claims through molecular modeling and mutational abrogation of the proposed interfaces. Native-MS has become a principal tool to determine the state of oligomeric organization of proteins. The work demonstrates that for multiprotein complexes, native-MS, coupled with molecular modeling and mutational perturbation, can provide an alternative route to render a detailed view of both the oligomeric states as well as the molecular interfaces involved. This is especially useful for large multiprotein complexes with large unstructured domains that make it recalcitrant to conventional structure determination approaches.

The flexibility and dynamics of protein disulfide isomerase.

We have studied the mobility of the multidomain folding catalyst, protein disulfide isomerase (PDI), by a coarse-graining approach based on flexibility. We analyze our simulations of yeast PDI (yPDI) using measures of backbone movement, relative positions and orientations of domains, and distances between functional sites. We find that there is interdomain flexibility at every interdomain junction but these show very different characteristics. The extent of interdomain flexibility is such that yPDI's two active sites can approach much more closely than is found in crystal structures-and indeed hinge motion to bring these sites into proximity is the lowest energy normal mode of motion of the protein. The flexibility predicted for yPDI (based on one structure) includes the other known conformation of yPDI and is consistent with (i) the mobility observed experimentally for mammalian PDI and (ii) molecular dynamics. We also observe intradomain flexibility and clear differences between the domains in their propensity for internal motion. Our results suggest that PDI flexibility enables it to interact with many different partner molecules of widely different sizes and shapes, and highlights considerable similarities of yPDI and mammalian PDI. Proteins 2016; 84:1776-1785. © 2016 Wiley Periodicals, Inc.

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.

Cyclic AMP-dependent protein lysine acylation in mycobacteria regulates fatty acid and propionate metabolism.

Acetylation of lysine residues is a posttranslational modification that is used by both eukaryotes and prokaryotes to regulate a variety of biological processes. Here we identify multiple substrates for the cAMP-dependent protein lysine acetyltransferase from Mycobacterium tuberculosis (KATmt). We demonstrate that a catalytically important lysine residue in a number of FadD (fatty acyl CoA synthetase) enzymes is acetylated by KATmt in a cAMP-dependent manner and that acetylation inhibits the activity of FadD enzymes. A sirtuin-like enzyme can deacetylate multiple FadDs, thus completing the regulatory cycle. Using a strain deleted for the KATmt ortholog in Mycobacterium bovis Bacillus Calmette-Guérin (BCG), we show for the first time that acetylation is dependent on intracellular cAMP levels. KATmt can utilize propionyl CoA as a substrate and, therefore, plays a critical role in alleviating propionyl CoA toxicity in mycobacteria by inactivating acyl CoA synthetase (ACS). The precision by which mycobacteria can regulate the metabolism of fatty acids in a cAMP-dependent manner appears to be unparalleled in other biological organisms and is ideally suited to adapt to the complex environment that pathogenic mycobacteria experience in the host.

Oligomeric organization of membrane proteins from native membranes at nanoscale spatial and single-molecule resolution.

The oligomeric organization of membrane proteins in native cell membranes is a critical regulator of their function. High-resolution quantitative measurements of oligomeric assemblies and how they change under different conditions are indispensable to understanding membrane protein biology. We report Native-nanoBleach, a total internal reflection fluorescence microscopy-based single-molecule photobleaching step analysis technique to determine the oligomeric distribution of membrane proteins directly from native membranes at an effective spatial resolution of ~10 nm. We achieved this by capturing target membrane proteins in native nanodiscs with their proximal native membrane environment using amphipathic copolymers. We applied Native-nanoBleach to quantify the oligomerization status of structurally and functionally diverse membrane proteins, including a receptor tyrosine kinase (TrkA) and a small GTPase (KRas) under growth-factor binding and oncogenic mutations, respectively. Our data suggest that Native-nanoBleach provides a sensitive, single-molecule platform to quantify membrane protein oligomeric distributions in native membranes under physiologically and clinically relevant conditions.

A proteome-wide quantitative platform for nanoscale spatially resolved extraction of membrane proteins into native nanodiscs.

The intricate molecular environment of the native membrane profoundly influences every aspect of membrane protein (MP) biology. Despite this, the most prevalent method of studying MPs uses detergent-like molecules that disrupt and remove this vital local membrane context. This severely impedes our ability to quantitatively decipher the local molecular context and comprehend its regulatory role in the structure, function, and biogenesis of MPs. Using a library of membrane-active polymers we have developed a platform for the high-throughput analysis of the membrane proteome. The platform enables near-complete spatially resolved extraction of target MPs directly from their endogenous membranes into native nanodiscs that maintain the local membrane context. We accompany this advancement with an open-access quantitative database that provides the most efficient extraction conditions of 2065 unique mammalian MPs. Our method enables rapid and near-complete extraction and purification of target MPs directly from their endogenous organellar membranes at physiological expression levels while maintaining the nanoscale local membrane environment. Going beyond the plasma membrane proteome, our platform enables extraction from any target organellar membrane including the endoplasmic reticulum, mitochondria, lysosome, Golgi, and even transient organelles such as the autophagosome. To further validate this platform we took several independent MPs and demonstrated how our resource can enable rapid extraction and purification of target MPs from different organellar membranes with high efficiency and purity. Further, taking two synaptic vesicle MPs, we show how the database can be extended to capture multiprotein complexes between overexpressed MPs. We expect these publicly available resources to empower researchers across disciplines to capture membrane 'nano-scoops' containing a target MP efficiently and interface with structural, functional, and other bioanalytical approaches. We demonstrate an example of this by combining our extraction platform with single-molecule TIRF imaging to demonstrate how it can enable rapid determination of homo-oligomeric states of target MPs in native cell membranes.

Determination of oligomeric organization of membrane proteins from native membranes at nanoscale-spatial and single-molecule resolution.

The oligomeric organization of membrane proteins in native cell membranes is a critical regulator of their function. High-resolution quantitative measurements of oligomeric assemblies and how they change under different conditions are indispensable to the understanding of membrane protein biology. We report a single-molecule imaging technique (Native-nanoBleach) to determine the oligomeric distribution of membrane proteins directly from native membranes at an effective spatial resolution of ∼10 nm. We achieved this by capturing target membrane proteins in "native nanodiscs" with their proximal native membrane environment using amphipathic copolymers. We established this method using structurally and functionally diverse membrane proteins with well-established stoichiometries. We then applied Native-nanoBleach to quantify the oligomerization status of a receptor tyrosine kinase (TrkA) and a small GTPase (KRas) under conditions of growth-factor binding or oncogenic mutations, respectively. Native-nanoBleach provides a sensitive, single-molecule platform to quantify membrane protein oligomeric distributions in native membranes at an unprecedented spatial resolution.

Gene body DNA hydroxymethylation restricts the magnitude of transcriptional changes during aging.

DNA hydroxymethylation (5hmC) is the most abundant oxidative derivative of DNA methylation (5mC) and is typically enriched at enhancers and gene bodies of transcriptionally active and tissue-specific genes. Although aberrant genomic 5hmC has been implicated in many age-related diseases, the functional role of the modification in aging remains largely unknown. Here, we report that 5hmC is stably enriched in multiple aged organs. Using the liver and cerebellum as model organs, we show that 5hmC accumulates in gene bodies associated with tissue-specific function and thereby restricts the magnitude of gene expression changes during aging. Mechanistically, we found that 5hmC decreases binding affinity of splicing factors compared to unmodified cytosine and 5mC, and is correlated with age-related alternative splicing events, suggesting RNA splicing as a potential mediator of 5hmC's transcriptionally restrictive function. Furthermore, we show that various age-related contexts, such as prolonged quiescence and senescence, are partially responsible for driving the accumulation of 5hmC with age. We provide evidence that this age-related function is conserved in mouse and human tissues, and further show that the modification is altered by regimens known to modulate lifespan. Our findings reveal that 5hmC is a regulator of tissue-specific function and may play a role in regulating longevity.

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.

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.

Structural Insights into the Regulation of Ca2+/Calmodulin-Dependent Protein Kinase II (CaMKII).

Ca2+/calmodulin-dependent protein kinase II (CaMKII) is a highly conserved serine/threonine kinase that is ubiquitously expressed throughout the human body. Specialized isoforms of CaMKII play key roles in neuronal and cardiac signaling. The distinctive holoenzyme architecture of CaMKII, with 12-14 kinase domains attached by flexible linkers to a central hub, poses formidable challenges for structural characterization. Nevertheless, progress in determining the structural mechanisms underlying CaMKII functions has come from studying the kinase domain and the hub separately, as well as from a recent electron microscopic investigation of the intact holoenzyme. In this review, we discuss our current understanding of the structure of CaMKII. We also discuss the intriguing finding that the CaMKII holoenzyme can undergo activation-triggered subunit exchange, a process that has implications for the potentiation and perpetuation of CaMKII activity.

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.

Flexible linkers in CaMKII control the balance between activating and inhibitory autophosphorylation.

The many variants of human Ca2+/calmodulin-dependent protein kinase II (CaMKII) differ in the lengths and sequences of disordered linkers connecting the kinase domains to the oligomeric hubs of the holoenzyme. CaMKII activity depends on the balance between activating and inhibitory autophosphorylation (on Thr 286 and Thr 305/306, respectively, in the human α isoform). Variation in the linkers could alter transphosphorylation rates within a holoenzyme and the balance of autophosphorylation outcomes. We show, using mammalian cell expression and a single-molecule assay, that the balance of autophosphorylation is flipped between CaMKII variants with longer and shorter linkers. For the principal isoforms in the brain, CaMKII-α, with a ~30 residue linker, readily acquires activating autophosphorylation, while CaMKII-ß, with a ~200 residue linker, is biased towards inhibitory autophosphorylation. Our results show how the responsiveness of CaMKII holoenzymes to calcium signals can be tuned by varying the relative levels of isoforms with long and short linkers.

Breakage of the oligomeric CaMKII hub by the regulatory segment of the kinase.

Ca2+/calmodulin-dependent protein kinase II (CaMKII) is an oligomeric enzyme with crucial roles in neuronal signaling and cardiac function. Previously, we showed that activation of CaMKII triggers the exchange of subunits between holoenzymes, potentially increasing the spread of the active state (Stratton et al., 2014; Bhattacharyya et al., 2016). Using mass spectrometry, we show now that unphosphorylated and phosphorylated peptides derived from the CaMKII-α regulatory segment bind to the CaMKII-α hub and break it into smaller oligomers. Molecular dynamics simulations show that the regulatory segments dock spontaneously at the interface between hub subunits, trapping large fluctuations in hub structure. Single-molecule fluorescence intensity analysis of CaMKII-α expressed in mammalian cells shows that activation of CaMKII-α results in the destabilization of the holoenzyme. Our results suggest that release of the regulatory segment by activation and phosphorylation allows it to destabilize the hub, producing smaller assemblies that might reassemble to form new holoenzymes.

Functional correlation of bacterial LuxS with their quaternary associations: interface analysis of the structure networks.

BACKGROUND: The genome of a wide variety of prokaryotes contains the luxS gene homologue, which encodes for the protein S-ribosylhomocysteinelyase (LuxS). This protein is responsible for the production of the quorum sensing molecule, AI-2 and has been implicated in a variety of functions such as flagellar motility, metabolic regulation, toxin production and even in pathogenicity. A high structural similarity is present in the LuxS structures determined from a few species. In this study, we have modelled the structures from several other species and have investigated their dimer interfaces. We have attempted to correlate the interface features of LuxS with the phenotypic nature of the organisms. RESULTS: The protein structure networks (PSN) are constructed and graph theoretical analysis is performed on the structures obtained from X-ray crystallography and on the modelled ones. The interfaces, which are known to contain the active site, are characterized from the PSNs of these homodimeric proteins. The key features presented by the protein interfaces are investigated for the classification of the proteins in relation to their function. From our analysis, structural interface motifs are identified for each class in our dataset, which showed distinctly different pattern at the interface of LuxS for the probiotics and some extremophiles. Our analysis also reveals potential sites of mutation and geometric patterns at the interface that was not evident from conventional sequence alignment studies. CONCLUSION: The structure network approach employed in this study for the analysis of dimeric interfaces in LuxS has brought out certain structural details at the side-chain interaction level, which were elusive from the conventional structure comparison methods. The results from this study provide a better understanding of the relation between the luxS gene and its functional role in the prokaryotes. This study also makes it possible to explore the potential direction towards the design of inhibitors of LuxS and thus towards a wide range of antimicrobials.

Hierarchical Organization Endows the Kinase Domain with Regulatory Plasticity.

The functional diversity of kinases enables specificity in cellular signal transduction. Yet how more than 500 members of the human kinome specifically receive regulatory inputs and convey information to appropriate substrates-all while using the common signaling output of phosphorylation-remains enigmatic. Here, we perform statistical co-evolution analysis, mutational scanning, and quantitative live-cell assays to reveal a hierarchical organization of the kinase domain that facilitates the orthogonal evolution of regulatory inputs and substrate outputs while maintaining catalytic function. We find that three quasi-independent "sectors"-groups of evolutionarily coupled residues-represent functional units in the kinase domain that encode for catalytic activity, substrate specificity, and regulation. Sector positions impact both disease and pharmacology: the catalytic sector is significantly enriched for somatic cancer mutations, and residues in the regulatory sector interact with allosteric kinase inhibitors. We propose that this functional architecture endows the kinase domain with inherent regulatory plasticity.

Deconstruction of the Ras switching cycle through saturation mutagenesis.

Ras proteins are highly conserved signaling molecules that exhibit regulated, nucleotide-dependent switching between active and inactive states. The high conservation of Ras requires mechanistic explanation, especially given the general mutational tolerance of proteins. Here, we use deep mutational scanning, biochemical analysis and molecular simulations to understand constraints on Ras sequence. Ras exhibits global sensitivity to mutation when regulated by a GTPase activating protein and a nucleotide exchange factor. Removing the regulators shifts the distribution of mutational effects to be largely neutral, and reveals hotspots of activating mutations in residues that restrain Ras dynamics and promote the inactive state. Evolutionary analysis, combined with structural and mutational data, argue that Ras has co-evolved with its regulators in the vertebrate lineage. Overall, our results show that sequence conservation in Ras depends strongly on the biochemical network in which it operates, providing a framework for understanding the origin of global selection pressures on proteins.