In the multi-domain protein adenylate kinase, domain insertion facilitates cooperative folding while accommodating function at domain interfaces

The Biochemical Impact of Extracting an Embedded Adenylate Kinase Domain Using Circular Permutation

Adenylate kinases (AKs) have evolved AMP-binding and lid domains that are encoded as continuous polypeptides embedded at different locations within the discontinuous polypeptide encoding the core domain. A prior study showed that AK homologues of different stabilities consistently retain cellular activity following circular permutation that splits a region with high energetic frustration within the AMP-binding domain into discontinuous fragments. Herein, we show that mesophilic and thermophilic AKs having this topological restructuring retain activity and substrate-binding characteristics of the parental AK. While permutation decreased the activity of both AK homologues at physiological temperatures, the catalytic activity of the thermophilic AK increased upon permutation when assayed >30 °C below the melting temperature of the native AK. The thermostabilities of the permuted AKs were uniformly lower than those of native AKs, and they exhibited multiphasic unfolding transitions, unlike the native AKs, which presented cooperative thermal unfolding. In addition, proteolytic digestion revealed that permutation destabilized each AK in differing manners, and mass spectrometry suggested that the new termini within the AMP-binding domain were responsible for the increased proteolysis sensitivity. These findings illustrate how changes in contact order can be used to tune enzyme activity and alter folding dynamics in multidomain enzymes.

Cooperativity and Folding Kinetics in a Multidomain Protein with Interwoven Chain Topology

Although a large percentage of eukaryotic proteomes consist of proteins with multiple domains, not much is known about their assembly mechanism, especially those with intricate native state architectures. Some have a complex topology in which the structural elements along the sequence are interwoven in such a manner that the domains cannot be separated by cutting at any location along the sequence. Such proteins are multiply connected multidomain proteins (MMPs) with the three-domain (NMP, LID, and CORE) phosphotransferase enzyme adenylate kinase (ADK) being an example. We devised a coarse-grained model to simulate ADK folding initiated by changing either the temperature or guanidinium chloride (GdmCl) concentration. The simulations reproduce the experimentally measured melting temperatures (associated with two equilibrium transitions), FRET efficiency as a function of GdmCl concentration, and the folding times quantitatively. Although the NMP domain orders independently, cooperative interactions between the LID and the CORE domains are required for complete assembly of the enzyme. Kinetic simulations show that, on the collapse time scale, multiple interconnected metastable states are populated, attesting to the folding heterogeneity. The network of kinetically connected states reveals that the CORE domain folds only after the NMP and LID domains, reflecting the interwoven nature of the chain topology.

Investigating the trade-off between folding and function in a multidomain Y-family DNA polymerase

The way in which multidomain proteins fold has been a puzzling question for decades. Until now, the mechanisms and functions of domain interactions involved in multidomain protein folding have been obscure. Here, we develop structure-based models to investigate the folding and DNA-binding processes of the multidomain Y-family DNA polymerase IV (DPO4). We uncover shifts in the folding mechanism among ordered domain-wise folding, backtracking folding, and cooperative folding, modulated by interdomain interactions. These lead to ‘U-shaped’ DPO4 folding kinetics. We characterize the effects of interdomain flexibility on the promotion of DPO4–DNA (un)binding, which probably contributes to the ability of DPO4 to bypass DNA lesions, which is a known biological role of Y-family polymerases. We suggest that the native topology of DPO4 leads to a trade-off between fast, stable folding and tight functional DNA binding. Our approach provides an effective way to quantitatively correlate the roles of protein interactions in conformational dynamics at the multidomain level.

A Change in Domain Cooperativity Drives the Function of Calnuc

With the increasing incidence of neurodegenerative disorders, there is an urgent need to understand the protein folding process. Examining the folding process of multidomain proteins remains a prime challenge, as their complex conformational dynamics make them highly susceptible to misfolding and/or aggregation. The presence of multiple domains in a protein can lead to interaction between the partially folded domains, thereby driving misfolding and/or aggregation. Calnuc is one such multidomain protein for which Ca²⁺ binding plays a pivotal role in governing its structural dynamics and stability and, presumably, in directing its interactions with other proteins. We demonstrate differential structural dynamics between the Ca²⁺-free and Ca²⁺-bound forms of calnuc. In the absence of Ca²⁺, full-length calnuc displays equilibrium structural transitions with four intermediate states, reporting a sum of the behavioral properties of its individual domains. Fragment-based studies illustrate the sequential events of structure adoption proceeding in the following order: EF domain followed by the NT and LZ domains in the apo state. On the other hand, Ca²⁺ binding increases domain cooperativity and enables the protein to fold as a single unit. Single-tryptophan mutant proteins, designed in a domain-dependent manner, confirm an increase in the number of interdomain interactions in the Ca²⁺-bound form as compared to the Ca²⁺-free state of the protein, thereby providing insight into its folding process. The attenuated domain crosstalk in apo-calnuc is likely to influence and regulate its physiologically important intermolecular interactions.

Successes and challenges in simulating the folding of large proteins

Computational simulations of protein folding can be used to interpret experimental folding results, to design new folding experiments, and to test the effects of mutations and small molecules on folding. However, whereas major experimental and computational progress has been made in understanding how small proteins fold, research on larger, multidomain proteins, which comprise the majority of proteins, is less advanced. Specifically, large proteins often fold via long-lived partially folded intermediates, whose structures, potentially toxic oligomerization, and interactions with cellular chaperones remain poorly understood. Molecular dynamics based folding simulations that rely on knowledge of the native structure can provide critical, detailed information on folding free energy landscapes, intermediates, and pathways. Further, increases in computational power and methodological advances have made folding simulations of large proteins practical and valuable. Here, using serpins that inhibit proteases as an example, we review native-centric methods for simulating the folding of large proteins. These synergistic approaches range from Gō and related structure-based models that can predict the effects of the native structure on folding to all-atom-based methods that include side-chain chemistry and can predict how disease-associated mutations may impact folding. The application of these computational approaches to serpins and other large proteins highlights the successes and limitations of current computational methods and underscores how computational results can be used to inform experiments. These powerful simulation approaches in combination with experiments can provide unique insights into how large proteins fold and misfold, expanding our ability to predict and manipulate protein folding.

Nature's Shortcut to Protein Folding

This Feature Article presents a view of the protein folding transition based on the hypothesis that Nature has built features within the sequences that enable a Shortcut to efficient folding. Nature's Shortcut is proposed to be the early establishment of a set of nonlocal weak contacts, constituting protein loops that significantly constrain regions of the collapsed disordered protein into a native-like low-resolution fluctuating topology of major sections of the backbone. Nature's establishment of this scaffold of nonlocal contacts is claimed to bypass what would otherwise be a nearly hopeless unaided search for the final three-dimensional structure in proteins longer than ∼100 amino acids. To support this main contention of the Feature Article, the loop hypothesis (LH) description of early folding events is experimentally tested with time-resolved Förster resonance energy transfer techniques for adenylate kinase, and the data are shown to be consistent with theoretical predictions from the sequential collapse model (SCM). The experimentally based LH and the theoretically founded SCM are argued to provide a unified picture of the role of nonlocal contacts as constituting Nature's Shortcut to protein folding. Importantly, the SCM is shown to reliably predict key nonlocal contacts utilizing only primary sequence information. This view on Nature's Shortcut is open to the protein community for further detailed assessment, including its practical consequences, by suitable application of advanced experimental and computational techniques.

Effect of ligand binding on a protein with a complex folding landscape

Ligand binding to a protein can stabilize it significantly against unfolding. The variation of the folding free energy, ΔΔG⁰, due to ligand binding can be derived from a simple reaction scheme involving exclusive binding to the native state. One obtains the following expression: , where K_d is the ligand dissociation constant and L is its concentration, R is the universal gas constant and T is the temperature. This expression has been shown to correctly describe experimental results on multiple proteins. In the current work we studied the effect of ligand binding on the stability of the multi-domain protein adenylate kinase from E. coli (AKE). Unfolding experiments were conducted using single-molecule FRET spectroscopy, which allowed us to directly obtain the fraction of unfolded protein in a model-free way from FRET efficiency histograms. Surprisingly, it was found that the effect of two inhibitors (Ap₅A and AMPPNP) and a substrate (AMP) on the stability of AKE was much smaller than expected based on K_d values obtained independently using microscale thermophoresis. To shed light on this issue, we measured the K_d for Ap₅A over a range of chemical denaturant concentrations where the protein is still folded. It was found that K_d increases dramatically over this range, likely due to the population of folding intermediates, whose binding to the ligand is much weaker than that of the native state. We propose that binding to folding intermediates may dominate the effect of ligands on the stability of multi-domain proteins, and could therefore have a strong impact on protein homeostasis in vivo.

A computational method for the design of nested proteins by loop-directed domain insertion

The computational design of novel nested proteins-in which the primary structure of one protein domain (insert) is flanked by the primary structure segments of another (parent)-would enable the generation of multifunctional proteins. Here we present a new algorithm, called Loop-Directed Domain Insertion (LooDo), implemented within the Rosetta software suite, for the purpose of designing nested protein domain combinations connected by flexible linker regions. Conformational space for the insert domain is sampled using large libraries of linker fragments for linker-to-parent domain superimposition followed by insert-to-linker superimposition. The relative positioning of the two domains (treated as rigid bodies) is sampled efficiently by a grid-based, mutual placement compatibility search. The conformations of the loop residues, and the identities of loop as well as interface residues, are simultaneously optimized using a generalized kinematic loop closure algorithm and Rosetta EnzymeDesign, respectively, to minimize interface energy. The algorithm was found to consistently sample near-native conformations and interface sequences for a benchmark set of structurally similar but functionally divergent domain-inserted enzymes from the α/β hydrolase superfamily, and discriminates well between native and nonnative conformations and sequences, although loop conformations tended to deviate from the native conformations. Furthermore, in cross-domain placement tests, native insert-parent domain combinations were ranked as the best-scoring structures compared to nonnative domain combinations. This algorithm should be broadly applicable to the design of multi-domain protein complexes with any combination of inserted or tandem domain connections.

Tracking the Catalytic Cycle of Adenylate Kinase by Ultraviolet Photodissociation Mass Spectrometry

The complex interplay of dynamic protein plasticity and specific side-chain interactions with substrate molecules that allows enzymes to catalyze reactions has yet to be fully unraveled. Top-down ultraviolet photodissociation (UVPD) mass spectrometry is used to track snapshots of conformational fluctuations in the phosphotransferase adenylate kinase (AK) throughout its active reaction cycle by characterization of complexes containing AK and each of four different adenosine phosphate ligands. Variations in efficiencies of UVPD backbone cleavages were consistently observed for three α-helices and the adenosine binding regions for AK complexes representing different steps of the catalytic cycle, implying that these stretches of the protein sample various structural microstates as the enzyme undergoes global open-to-closed transitions. Focusing on the conformational impact of recruiting or releasing the Mg2+ cofactor highlights two loop regions for which fragmentation increases upon UVPD, signaling an increase in loop flexibility as the metal cation disrupts the loop interactions with the substrate ligands. Additionally, the observation of holo ions and variations in UVPD backbone cleavage efficiency at R138 implicate this conserved active site residue in stabilizing the donor phosphoryl group during catalysis. This study showcases the utility of UVPD-MS to provide insight into conformational fluctuations of single residues for active enzymes.

Power of protein/tRNA functional assembly against aberrant aggregation

Understanding the mechanisms of protein oligomerization and aggregation is a major concern for biotechnology and medical purposes. However, significant challenges remain in determining the mechanism of formation of these superstructures and the environmental factors that can precisely modulate them. Notably the role that a functional ligand plays in the process of protein aggregation is largely unexplored. We herein address these issues with an original flavin-dependent RNA methyltransferase (TrmFO) used as a protein model since this protein employs a complex set of cofactors and ligands for catalysis. Here, we show that TrmFO carries an unstable protein structure that can partially mis-unfold leading to either formation of irregular and nonfunctional soluble oligomers endowed with hyper-thermal stability or large amorphous aggregates in the presence of salts. Mutagenesis confirmed that this peculiarity is an intrinsic property of a polypeptide and it is independent of the flavin coenzyme. Structural characterization and kinetic studies identified several regions of the protein that enjoy conformational changes and more particularly pinpointed the N-terminal subdomain as being a key element in the mechanisms of oligomerization and aggregation. Only stabilization of this region via tRNA suppresses these aberrant protein states. Although protein chaperones emerged as major actors against aggregation, our study emphasizes that other powerful mechanisms exist such as the stabilizing effect of functional assemblies that provide an additional layer of protection against the instability of the proteome.

Temperature stability of proteins: Analysis of irreversible denaturation using isothermal calorimetry

The structural stability of proteins has been traditionally studied under conditions in which the folding/unfolding reaction is reversible, since thermodynamic parameters can only be determined under these conditions. Achieving reversibility conditions in temperature stability experiments has often required performing the experiments at acidic pH or other nonphysiological solvent conditions. With the rapid development of protein drugs, the fastest growing segment in the pharmaceutical industry, the need to evaluate protein stability under formulation conditions has acquired renewed urgency. Under formulation conditions and the required high protein concentration (∼100 mg/mL), protein denaturation is irreversible and frequently coupled to aggregation and precipitation. In this article, we examine the thermal denaturation of hen egg white lysozyme (HEWL) under irreversible conditions and concentrations up to 100 mg/mL using several techniques, especially isothermal calorimetry which has been used to measure the enthalpy and kinetics of the unfolding and aggregation/precipitation at 12°C below the transition temperature measured by DSC. At those temperatures the rate of irreversible protein denaturation and aggregation of HEWL is measured to be on the order of 1 day^-1 . Isothermal calorimetry appears a suitable technique to identify buffer formulation conditions that maximize the long term stability of protein drugs.

Molecular simulations of substrate release and coupled conformational motions in adenylate kinase

Conformational opening coupled substrate release is believed to be related to the rate limiting step in the catalysis cycle of the adenylate kinase. However, it is still unclear how the substrate dissociates from its active site and how the substrate release is coupled to conformational changes of the kinase. In this work, by using metadynamics simulations, we investigated the ADP release process and the coupled protein dynamics. We found that the ADP release involves overcoming a high free energy barrier, and protonation of the [Formula: see text]-phosphate of the ADP molecules can drastically reduce the barrier height, therefore, promote the ADP release. We identified several key residues contributing to the high free energy barrier. We also showed that the ADP attached to LID domain leaves the binding pocket earlier than the one attached to the NMP domain. We further observed that the ADP release is accompanied by almost fully opening of the LID domain and partially opening of the NMP domain. Our results provide insight into the molecular mechanism of the substrate release of adenylate kinase and the coupled conformational motions.

Using the folding landscapes of proteins to understand protein function

Proteins fold on a biologically-relevant timescale because of a funnel-shaped energy landscape. This landscape is sculpted through evolution by selecting amino-acid sequences that stabilize native interactions while suppressing stable non-native interactions that occur during folding. However, there is strong evolutionary selection for functional residues and these cannot be chosen to optimize folding. Their presence impacts the folding energy landscape in a variety of ways. Here, we survey the effects of functional residues on folding by providing several examples. We then review how such effects can be detected computationally and be used as assays for protein function. Overall, an understanding of how functional residues modulate folding should provide insights into the design of natural proteins and their homeostasis.

Effects of heme binding on myoglobin folding: Coarse grained molecular simulations

Many proteins contain cofactors, such as heme, ATP and metal ions. Binding of cofactors is not only essential for their biological functions, but also can reshape the intrinsic energy landscape of protein molecules and modulate the folding and stability. However, the molecular mechanism of cofactor coupled protein folding is not well understood. In this work, we study the cofactor coupled folding of myoglobin, which is a typical cofactor (heme) containing protein, by performing molecular dynamics simulations with a structure-based protein model developed based on the energy landscape theory. We showed that the heme binding increases the stability of the myoglobin. More importantly, the heme binding tends to increase the protein folding cooperativity, and switch the folding process from a “three-state” mechanism to a “two-state” mechanism. We also showed that the folding pathways of the myoglobin can be modulated by the heme binding. By performing comparative simulations, we revealed that the above effects of heme binding are resulted from the heme induced folding of F-helix, which is otherwise unstructured at apo state, and the heme mediated contacting interactions around the heme binding site. The simulation results are consistent with available experimental data, and provide insights into the molecular mechanism of the effects of cofactor binding on the protein folding and stability.

Protein-Protein Interactions Mediated by Helical Tertiary Structure Motifs

The modulation of protein-protein interactions (PPIs) by means of creating or stabilizing secondary structure conformations is a rapidly growing area of research. Recent success in the inhibition of difficult PPIs by secondary structure mimetics also points to potential limitations, because often, specific cases require tertiary structure mimetics. To streamline protein structure-based inhibitor design, we have previously described the examination of protein complexes in the Protein Data Bank where α-helices or β-strands form critical contacts. Here, we examined coiled coils and helix bundles that mediate complex formation to create a platform for the discovery of potential tertiary structure mimetics. Though there has been extensive analysis of coiled coil motifs, the interactions between pre-formed coiled coils and globular proteins have not been systematically analyzed. This article identifies critical features of these helical interfaces with respect to coiled coil and other helical PPIs. We expect the analysis to prove useful for the rational design of modulators of this fundamental class of protein assemblies.

Nonspecific yet decisive: Ubiquitination can affect the native-state dynamics of the modified protein

Ubiquitination is one of the most common post-translational modifications of proteins, and mediates regulated protein degradation among other cellular processes. A fundamental question regarding the mechanism of protein ubiquitination is whether and how ubiquitin affects the biophysical nature of the modified protein. For some systems, it was shown that the position of ubiquitin within the attachment site is quite flexible and ubiquitin does not specifically interact with its substrate. Nevertheless, it was revealed that polyubiquitination can decrease the thermal stability of the modified protein in a site-specific manner because of alterations of the thermodynamic properties of the folded and unfolded states. In this study, we used detailed atomistic simulations to focus on the molecular effects of ubiquitination on the native structure of the modified protein. As a model, we used Ubc7, which is an E2 enzyme whose in vivo ubiquitination process is well characterized and known to lead to degradation. We found that, despite the lack of specific direct interactions between the ubiquitin moiety and Ubc7, ubiquitination decreases the conformational flexibility of certain regions of the substrate Ubc7 protein, which reduces its entropy and thus destabilizes it. The strongest destabilizing effect was observed for systems in which Lys48-linked tetra-ubiquitin was attached to sites used for in vivo degradation. These results reveal how changes in the configurational entropy of the folded state may modulate the stability of the protein's native state. Overall, our results imply that ubiquitination can modify the biophysical properties of the attached protein in the folded state and that, in some proteins, different ubiquitination sites will lead to different biophysical outcomes. We propose that this destabilizing effect of polyubiquitin on the substrate is linked to the functions carried out by the modification, and in particular, regulatory control of protein half-life through proteasomal degradation.