Overlap between folding and functional energy landscapes for adenylate kinase conformational change

Allostery can convert binding free energies into concerted domain motions in enzymes

Enzymatic mechanisms are typically inferred from structural data. However, understanding enzymes require unravelling the intricate dynamic interplay between dynamics, conformational substates, and multiple protein structures. Here, we use single-molecule nanopore analysis to investigate the catalytic conformational changes of adenylate kinase (AK), an enzyme that catalyzes the interconversion of various adenosine phosphates (ATP, ADP, and AMP). Kinetic analysis validated by hidden Markov models unravels the details of domain motions during catalysis. Our findings reveal that allosteric interactions between ligands and cofactor enable converting binding energies into directional conformational changes of the two catalytic domains of AK. These coordinated motions emerged to control the exact sequence of ligand binding and the affinity for the three different substrates, thereby guiding the reactants along the reaction coordinates. Interestingly, we find that about 10% of enzymes show altered allosteric regulation and ligand affinities, indicating that a subset of enzymes folds in alternative catalytically active forms. Since molecules or proteins might be able to selectively stabilize one of the folds, this observation suggests an evolutionary path for allostery in enzymes. In AK, this complex catalytic framework has likely emerged to prevent futile ATP/ADP hydrolysis and to regulate the enzyme for different energy needs of the cell.

Evolutionary-Scale Enzymology Enables Biochemical Constant Prediction Across a Multi-Peaked Catalytic Landscape

Quantitatively mapping enzyme sequence-catalysis landscapes remains a critical challenge in understanding enzyme function, evolution, and design. Here, we expand an emerging microfluidic platform to measure catalytic constants-k cat and K M-for hundreds of diverse naturally occurring sequences and mutants of the model enzyme Adenylate Kinase (ADK). This enables us to dissect the sequence-catalysis landscape's topology, navigability, and mechanistic underpinnings, revealing distinct catalytic peaks organized by structural motifs. These results challenge long-standing hypotheses in enzyme adaptation, demonstrating that thermophilic enzymes are not slower than their mesophilic counterparts. Combining the rich representations of protein sequences provided by deep-learning models with our custom high-throughput kinetic data yields semi-supervised models that significantly outperform existing models at predicting catalytic parameters of naturally occurring ADK sequences. Our work demonstrates a promising strategy for dissecting sequence-catalysis landscapes across enzymatic evolution and building family-specific models capable of accurately predicting catalytic constants, opening new avenues for enzyme engineering and functional prediction.

Magnesium induced structural reorganization in the active site of adenylate kinase

Phosphoryl transfer is a fundamental reaction in cellular signaling and metabolism that requires Mg2+ as an essential cofactor. While the primary function of Mg2+ is electrostatic activation of substrates, such as ATP, the full spectrum of catalytic mechanisms exerted by Mg2+ is not known. In this study, we integrate structural biology methods, molecular dynamic (MD) simulations, phylogeny, and enzymology assays to provide molecular insights into Mg2+-dependent structural reorganization in the active site of the metabolic enzyme adenylate kinase. Our results demonstrate that Mg2+ induces a conformational rearrangement of the substrates (ATP and ADP), resulting in a 30° adjustment of the angle essential for reversible phosphoryl transfer, thereby optimizing it for catalysis. MD simulations revealed transitions between conformational substates that link the fluctuation of the angle to large-scale enzyme dynamics. The findings contribute detailed insight into Mg2+ activation of enzymes and may be relevant for reversible and irreversible phosphoryl transfer reactions.

Towards the understanding of molecular motors and its relationship with local unfolding

Molecular motors are machines essential for life since they convert chemical energy into mechanical work. However, the precise mechanism by which nucleotide binding, catalysis, or release of products is coupled to the work performed by the molecular motor is still not entirely clear. This is due, in part, to a lack of understanding of the role of force in the mechanical-structural processes involved in enzyme catalysis. From a mechanical perspective, one promising hypothesis is the Haldane-Pauling hypothesis which considers the idea that part of the enzymatic catalysis is strain-induced. It suggests that enzymes cannot be efficient catalysts if they are fully complementary to the substrates. Instead, they must exert strain on the substrate upon binding, using enzyme-substrate energy interaction (binding energy) to accelerate the reaction rate. A novel idea suggests that during catalysis, significant strain energy is built up, which is then released by a local unfolding/refolding event known as 'cracking'. Recent evidence has also shown that in catalytic reactions involving conformational changes, part of the heat released results in a center-of-mass acceleration of the enzyme, raising the possibility that the heat released by the reaction itself could affect the enzyme's integrity. Thus, it has been suggested that this released heat could promote or be linked to the cracking seen in proteins such as adenylate kinase (AK). We propose that the energy released as a consequence of ligand binding/catalysis is associated with the local unfolding/refolding events (cracking), and that this energy is capable of driving the mechanical work.

Associating protein sequence positions with the modulation of quantitative phenotypes

Although protein sequences encode the information for folding and function, understanding their link is not an easy task. Unluckily, the prediction of how specific amino acids contribute to these features is still considerably impaired. Here, we developed a simple algorithm that finds positions in a protein sequence with potential to modulate the studied quantitative phenotypes. From a few hundred protein sequences, we perform multiple sequence alignments, obtain the per-position pairwise differences for both the sequence and the observed phenotypes, and calculate the correlation between these last two quantities. We tested our methodology with four cases: archaeal Adenylate Kinases and the organisms optimal growth temperatures, microbial rhodopsins and their maximal absorption wavelengths, mammalian myoglobins and their muscular concentration, and inhibition of HIV protease clinical isolates by two different molecules. We found from 3 to 10 positions tightly associated with those phenotypes, depending on the studied case. We showed that these correlations appear using individual positions but an improvement is achieved when the most correlated positions are jointly analyzed. Noteworthy, we performed phenotype predictions using a simple linear model that links per-position divergences and differences in the observed phenotypes. Predictions are comparable to the state-of-art methodologies which, in most of the cases, are far more complex. All of the calculations are obtained at a very low information cost since the only input needed is a multiple sequence alignment of protein sequences with their associated quantitative phenotypes. The diversity of the explored systems makes our work a valuable tool to find sequence determinants of biological activity modulation and to predict various functional features for uncharacterized members of a protein family.

Elucidating Dynamics of Adenylate Kinase from Enzyme Opening to Ligand Release

This study explores ligand-driven conformational changes in adenylate kinase (AK), which is known for its open-to-close conformational transitions upon ligand binding and release. By utilizing string free energy simulations, we determine the free energy profiles for both enzyme opening and ligand release and compare them with profiles from the apoenzyme. Results reveal a three-step ligand release process, which initiates with the opening of the adenosine triphosphate-binding subdomain (ATP lid), followed by ligand release and concomitant opening of the adenosine monophosphate-binding subdomain (AMP lid). The ligands then transition to nonspecific positions before complete dissociation. In these processes, the first step is energetically driven by ATP lid opening, whereas the second step is driven by ATP release. In contrast, the AMP lid opening and its ligand release make minor contributions to the total free energy for enzyme opening. Regarding the ligand binding mechanism, our results suggest that AMP lid closure occurs via an induced-fit mechanism triggered by AMP binding, whereas ATP lid closure follows conformational selection. This difference in the closure mechanisms provides an explanation with implications for the debate on ligand-driven conformational changes of AK. Additionally, we determine an X-ray structure of an AK variant that exhibits significant rearrangements in the stacking of catalytic arginines, explaining its reduced catalytic activity. In the context of apoenzyme opening, the sequence of events is different. Here, the AMP lid opens first while the ATP lid remains closed, and the free energy associated with ATP lid opening varies with orientation, aligning with the reported AK opening and closing rate heterogeneity. Finally, this study, in conjunction with our previous research, provides a comprehensive view of the intricate interplay between various structural elements, ligands, and catalytic residues that collectively contribute to the robust catalytic power of the enzyme.

Conformational dynamics of adenylate kinase in crystals

Adenylate kinase is a ubiquitous enzyme in living systems and undergoes dramatic conformational changes during its catalytic cycle. For these reasons, it is widely studied by genetic, biochemical, and biophysical methods, both experimental and theoretical. We have determined the basic crystal structures of three differently liganded states of adenylate kinase from Methanotorrus igneus, a hyperthermophilic organism whose adenylate kinase is a homotrimeric oligomer. The multiple copies of each protomer in the asymmetric unit of the crystal provide a unique opportunity to study the variation in the structure and were further analyzed using advanced crystallographic refinement methods and analysis tools to reveal conformational heterogeneity and, thus, implied dynamic behaviors in the catalytic cycle.

Temperature as a modulator of allosteric motions and crosstalk in mesophilic and thermophilic enzymes

Mesophilic and thermophilic enzyme counterparts are often studied to understand how proteins function under harsh conditions. To function well outside of standard temperature ranges, thermophiles often tightly regulate their structural ensemble through intra-protein communication (via allostery) and altered interactions with ligands. It has also become apparent in recent years that the enhancement or diminution of allosteric crosstalk can be temperature-dependent and distinguish thermophilic enzymes from their mesophilic paralogs. Since most studies of allostery utilize chemical modifications from pH, mutations, or ligands, the impact of temperature on allosteric function is comparatively understudied. Here, we discuss the biophysical methods, as well as critical case studies, that dissect temperature-dependent function of mesophilic-thermophilic enzyme pairs and their allosteric regulation across a range of temperatures.

Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge-driven or artificial intelligence-driven and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary-structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models that contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly while refining them against the density maps. We introduce such an adaptive decision-making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two prototypical protein systems, adenylate kinase and carbon monoxide dehydrogenase. For these test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2-3 Å map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that are inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from the chimpanzee adenovirus. The overhead for decision making remains low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install R-MDFF on different computing environments, from local Linux-based workstations to high-performance computing environments.

Thermostable Human Basic Fibroblast Growth Factor (TS-bFGF) Engineered with a Disulfide Bond Demonstrates Superior Culture Outcomes in Human Pluripotent Stem Cell

Human embryonic stem cells (hESCs) and induced pluripotent stem cells (iPSCs) can differentiate into various tissues and are an essential source of various disease models and therapeutics. Various growth factors are required in order to culture pluripotent stem cells, among which basic fibroblast growth factor (bFGF) is essential for maintaining stem cell ability. However, bFGF has a short half-life (8 h) under normal mammalian cell culture conditions, and its activity decreases after 72 h, posing a serious problem in the production of high-quality stem cells. Here, we evaluated the various functions of pluripotent stem cells (PSCs) by utilizing an engineered thermostable bFGF (TS-bFGF) that is thermally stable and maintains activity longer under mammalian culture conditions. PSCs cultured with TS-bFGF showed better proliferation, stemness, morphology, and differentiation than cells cultured with wild-type bFGF. In light of the importance of stem cells in a wide range of applications in the medical and biotechnology fields, we anticipate that TS-bFGF, as a thermostable and long-acting bFGF, can play a key role in securing high-quality stem cells through various sets of stem cell culture processes.

Allosteric communication between ligand binding domains modulates substrate inhibition in adenylate kinase

Enzymes play a vital role in life processes; they control chemical reactions and allow functional cycles to be synchronized. Many enzymes harness large-scale motions of their domains to achieve tremendous catalytic prowess and high selectivity for specific substrates. One outstanding example is provided by the three-domain enzyme adenylate kinase (AK), which catalyzes phosphotransfer between ATP to AMP. Here we study the phenomenon of substrate inhibition by AMP and its correlation with domain motions. Using single-molecule FRET spectroscopy, we show that AMP does not block access to the ATP binding site, neither by competitive binding to the ATP cognate site nor by directly closing the LID domain. Instead, inhibitory concentrations of AMP lead to a faster and more cooperative domain closure by ATP, leading in turn to an increased population of the closed state. The effect of AMP binding can be modulated through mutations throughout the structure of the enzyme, as shown by the screening of an extensive AK mutant library. The mutation of multiple conserved residues reduces substrate inhibition, suggesting that substrate inhibition is an evolutionary well conserved feature in AK. Combining these insights, we developed a model that explains the complex activity of AK, particularly substrate inhibition, based on the experimentally observed opening and closing rates. Notably, the model indicates that the catalytic power is affected by the microsecond balance between the open and closed states of the enzyme. Our findings highlight the crucial role of protein motions in enzymatic activity.

Mechanistic Basis for a Connection between the Catalytic Step and Slow Opening Dynamics of Adenylate Kinase

Escherichia coli adenylate kinase (AdK) is a small, monomeric enzyme that synchronizes the catalytic step with the enzyme's conformational dynamics to optimize a phosphoryl transfer reaction and the subsequent release of the product. Guided by experimental measurements of low catalytic activity in seven single-point mutation AdK variants (K13Q, R36A, R88A, R123A, R156K, R167A, and D158A), we utilized classical mechanical simulations to probe mutant dynamics linked to product release, and quantum mechanical and molecular mechanical calculations to compute a free energy barrier for the catalytic event. The goal was to establish a mechanistic connection between the two activities. Our calculations of the free energy barriers in AdK variants were in line with those from experiments, and conformational dynamics consistently demonstrated an enhanced tendency toward enzyme opening. This indicates that the catalytic residues in the wild-type AdK serve a dual role in this enzyme's function─one to lower the energy barrier for the phosphoryl transfer reaction and another to delay enzyme opening, maintaining it in a catalytically active, closed conformation for long enough to enable the subsequent chemical step. Our study also discovers that while each catalytic residue individually contributes to facilitating the catalysis, R36, R123, R156, R167, and D158 are organized in a tightly coordinated interaction network and collectively modulate AdK's conformational transitions. Unlike the existing notion of product release being rate-limiting, our results suggest a mechanistic interconnection between the chemical step and the enzyme's conformational dynamics acting as the bottleneck of the catalytic process. Our results also suggest that the enzyme's active site has evolved to optimize the chemical reaction step while slowing down the overall opening dynamics of the enzyme.

Protein dynamics: The future is bright and complicated!

Biological life depends on motion, and this manifests itself in proteins that display motion over a formidable range of time scales spanning from femtoseconds vibrations of atoms at enzymatic transition states, all the way to slow domain motions occurring on micro to milliseconds. An outstanding challenge in contemporary biophysics and structural biology is a quantitative understanding of the linkages among protein structure, dynamics, and function. These linkages are becoming increasingly explorable due to conceptual and methodological advances. In this Perspective article, we will point toward future directions of the field of protein dynamics with an emphasis on enzymes. Research questions in the field are becoming increasingly complex such as the mechanistic understanding of high-order interaction networks in allosteric signal propagation through a protein matrix, or the connection between local and collective motions. In analogy to the solution to the "protein folding problem," we argue that the way forward to understanding these and other important questions lies in the successful integration of experiment and computation, while utilizing the present rapid expansion of sequence and structure space. Looking forward, the future is bright, and we are in a period where we are on the doorstep to, at least in part, comprehend the importance of dynamics for biological function.

Insights into the evolution of enzymatic specificity and catalysis: From Asgard archaea to human adenylate kinases

Enzymatic catalysis is critically dependent on selectivity, active site architecture, and dynamics. To contribute insights into the interplay of these properties, we established an approach with NMR, crystallography, and MD simulations focused on the ubiquitous phosphotransferase adenylate kinase (AK) isolated from Odinarchaeota (OdinAK). Odinarchaeota belongs to the Asgard archaeal phylum that is believed to be the closest known ancestor to eukaryotes. We show that OdinAK is a hyperthermophilic trimer that, contrary to other AK family members, can use all NTPs for its phosphorylation reaction. Crystallographic structures of OdinAK-NTP complexes revealed a universal NTP-binding motif, while 19F NMR experiments uncovered a conserved and rate-limiting dynamic signature. As a consequence of trimerization, the active site of OdinAK was found to be lacking a critical catalytic residue and is therefore considered to be "atypical." On the basis of discovered relationships with human monomeric homologs, our findings are discussed in terms of evolution of enzymatic substrate specificity and cold adaptation.

Role of Repeated Conformational Transitions in Substrate Binding of Adenylate Kinase

The catalytic cycle of the enzyme adenylate kinase involves large conformational motions between open and closed states. A previous single-molecule experiment showed that substrate binding tends to accelerate both the opening and the closing rates and that a single turnover event often involves multiple rounds of conformational switching. In this work, we showed that the repeated conformational transitions of adenylate kinase are essential for the relaxation of incorrectly bound substrates into the catalytically competent conformation by combining all-atom and coarse-grained molecular simulations. In addition, free energy calculations based on all-atom and coarse-grained models demonstrated that the enzyme with incorrectly bound substrates has much a lower free energy barrier for domain opening compared to that with the correct substrate conformation, which may explain the the acceleration of the domain opening rate by substrate binding. The results of this work provide mechanistic understanding to previous experimental observations and shed light onto the interplay between conformational dynamics and enzyme catalysis.

Machine learning/molecular dynamic protein structure prediction approach to investigate the protein conformational ensemble

Proteins exist in several different conformations. These structural changes are often associated with fluctuations at the residue level. Recent findings show that co-evolutionary analysis coupled with machine-learning techniques improves the precision by providing quantitative distance predictions between pairs of residues. The predicted statistical distance distribution from Multi Sequence Analysis reveals the presence of different local maxima suggesting the flexibility of key residue pairs. Here we investigate the ability of the residue-residue distance prediction to provide insights into the protein conformational ensemble. We combine deep learning approaches with mechanistic modeling to a set of proteins that experimentally showed conformational changes. The predicted protein models were filtered based on energy scores, RMSD clustering, and the centroids selected as the lowest energy structure per cluster. These models were compared to the experimental-Molecular Dynamics (MD) relaxed structure by analyzing the backbone residue torsional distribution and the sidechain orientations. Our pipeline allows to retrieve the experimental structural dynamics experimentally represented by different X-ray conformations for the same sequence as well the conformational space observed with the MD simulations. We show the potential correlation between the experimental structure dynamics and the predicted model ensemble demonstrating the susceptibility of the current state-of-the-art methods in protein folding and dynamics prediction and pointing out the areas of improvement.

AK2 is an AMP-sensing negative regulator of BRAF in tumorigenesis

The RAS-BRAF signaling is a major pathway of cell proliferation and their mutations are frequently found in human cancers. Adenylate kinase 2 (AK2), which modulates balance of adenine nucleotide pool, has been implicated in cell death and cell proliferation independently of its enzyme activity. Recently, the role of AK2 in tumorigenesis was in part elucidated in some cancer types including lung adenocarcinoma and breast cancer, but the underlying mechanism is not clear. Here, we show that AK2 is a BRAF-suppressor. In in vitro assays and cell model, AK2 interacted with BRAF and inhibited BRAF activity and downstream ERK phosphorylation. Energy-deprived conditions in cell model and the addition of AMP to cell lysates strengthened the AK2-BRAF interaction, suggesting that AK2 is involved in the regulation of BRAF activity in response to cell metabolic state. AMP facilitated the AK2-BRAF complex formation through binding to AK2. In a panel of HCC cell lines, AK2 expression was inversely correlated with ERK/MAPK activation, and AK2-knockdown or -knockout increased BRAF activity and promoted cell proliferation. Tumors from HCC patients showed low-AK2 protein expression and increased ERK activation compared to non-tumor tissues and the downregulation of AK2 was also verified by two microarray datasets (TCGA-LIHC and GSE14520). Moreover, AK2/BRAF interaction was abrogated by RAS activation in in vitro assay and cell model and in a mouse model of HRAS^G12V-driven HCC, and AK2 ablation promoted tumor growth and BRAF activity. AK2 also bound to BRAF inhibitor-insensitive BRAF mutants and attenuated their activities. These findings indicate that AK2 monitoring cellular AMP levels is indeed a negative regulator of BRAF, linking the metabolic status to tumor growth.

Moderate molecular recognitions on ZnO m-plane and their selective capture/release of bio-related phosphoric acids

Herein, we explore the hidden molecular recognition abilities of ZnO nanowires uniformly grown on the inner surface of an open tubular fused silica capillary via liquid chromatography. Chromatographic evaluation revealed that ZnO nanowires showed a stronger intermolecular interaction with phenylphosphoric acid than any other monosubstituted benzene. Furthermore, ZnO nanowires specifically recognized the phosphate groups present in nucleotides even in the aqueous mobile phase, and the intermolecular interaction increased with the number of phosphate groups. This discrimination of phosphate groups in nucleotides was unique to the rich (101̄0) m-plane of ZnO nanowires with a moderate hydrophilicity and negative charge. The discrimination could be evidenced by the changes in the infrared bands of the phosphate groups on nucleotides on ZnO nanowires. Finally, as an application of the molecular recognition, nucleotides were separated by the number of phosphate groups, utilizing optimized gradient elution on ZnO nanowire column. Thus, the present results elucidate the unique and versatile molecular selectivity of well-known ZnO nanostructures for the capture and separation of biomolecules.

Tracking the ATP-binding response in adenylate kinase in real time

The biological function of proteins is critically dependent on dynamics inherent to the native structure. Such structural dynamics obey a predefined order and temporal timing to execute the specific reaction. Determination of the cooperativity of key structural rearrangements requires monitoring protein reactions in real time. In this work, we used time-resolved x-ray solution scattering (TR-XSS) to visualize structural changes in the Escherichia coli adenylate kinase (AdK) enzyme upon laser-induced activation of a protected ATP substrate. A 4.3-ms transient intermediate showed partial closing of both the ATP- and AMP-binding domains, which indicates a cooperative closing mechanism. The ATP-binding domain also showed local unfolding and breaking of an Arg¹³¹-Asp¹⁴⁶ salt bridge. Nuclear magnetic resonance spectroscopy data identified similar unfolding in an Arg¹³¹Ala AdK mutant, which refolded in a closed, substrate-binding conformation. The observed structural dynamics agree with a “cracking mechanism” proposed to underlie global structural transformation, such as allostery, in proteins.

The lipopolysaccharide-transporter complex LptB2FG also displays adenylate kinase activity in vitro dependent on the binding partners LptC/LptA

Lipopolysaccharide (LPS) is an essential glycolipid that covers the surface of gram-negative bacteria. The transport of LPS involves a dedicated seven-protein transporter system called the lipopolysaccharide transport system (Lpt) machinery that physically spans the entire cell envelope. The LptB₂FG complex is an ABC transporter that hydrolyzes ATP to extract LPS from the inner membrane for transport to the outer membrane. Here, we extracted LptB₂FG directly from the inner membrane with its original lipid environment using styrene-maleic acid polymers. We found that styrene-maleic acid polymers–LptB₂FG in nanodiscs display not only ATPase activity but also a previously uncharacterized adenylate kinase (AK) activity, as it catalyzed phosphotransfer between two ADP molecules to generate ATP and AMP. The ATPase and AK activities of LptB₂FG were both stimulated by the interaction on the periplasmic side with the periplasmic LPS transport proteins LptC and LptA and inhibited by the presence of the LptC transmembrane helix. We determined that the isolated ATPase module (LptB) had weak AK activity in the absence of transmembrane proteins LptF and LptG, and one mutation in LptB that weakens its affinity for ADP led to AK activity similar to that of fully assembled complex. Thus, we conclude that LptB₂FG is capable of producing ATP from ADP, depending on the assembly of the Lpt bridge, and that this AK activity might be important to ensure efficient LPS transport in the fully assembled Lpt system.

Milligram scale expression, refolding, and purification of Bombyx mori cocoonase using a recombinant E. coli system

Silk is one of the most versatile biomaterials with signature properties of outstanding mechanical strength and flexibility. A potential avenue for developing more environmentally friendly silk production is to make use of the silk moth (Bombyx mori) cocoonase, this will at the same time increase the possibility for using the byproduct, sericin, as a raw material for other applications. Cocoonase is a serine protease utilized by the silk moth to soften the cocoon to enable its escape after completed metamorphosis. Cocoonase selectively degrades the glue protein of the cocoon, sericin, without affecting the silk-fiber made of the protein fibroin. Cocoonase can be recombinantly produced in E. coli, however, it is exclusively found as insoluble inclusion bodies. To solve this problem and to be able to utilize the benefits associated with an E. coli based expression system, we have developed a protocol that enables the production of soluble and functional protease in the milligram/liter scale. The core of the protocol is refolding of the protein in a buffer with a redox potential that is optimized for formation of native and intramolecular di-sulfide bridges. The redox potential was balanced with defined concentrations of reduced and oxidized glutathione. This E.coli based production protocol will, in addition to structure determination, also enable modification of cocoonase both in terms of catalytic function and stability. These factors will be valuable components in the development of alternate silk production methodology.

Moving beyond static snapshots: Protein dynamics and the Protein Data Bank

Proteins are the molecular machines of living systems. Their dynamics are an intrinsic part of their evolutionary selection in carrying out their biological functions. Although the dynamics are more difficult to observe than a static, average structure, we are beginning to observe these dynamics and form sound mechanistic connections between structure, dynamics, and function. This progress is highlighted in case studies from myoglobin and adenylate kinase to the ribosome and molecular motors where these molecules are being probed with a multitude of techniques across many timescales. New approaches to time-resolved crystallography are allowing simple “movies” to be taken of proteins in action, and new methods of mapping the variations in cryo-electron microscopy are emerging to reveal a more complete description of life’s machines. The results of these new methods are aided in their dissemination by continual improvements in curation and distribution by the Protein Data Bank and their partners around the world.

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. Here, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

The change of conditions does not affect Ros87 downhill folding mechanism

Downhill folding has been defined as a unique thermodynamic process involving a conformations ensemble that progressively loses structure with the decrease of protein stability. Downhill folders are estimated to be rather rare in nature as they miss an energetically substantial folding barrier that can protect against aggregation and proteolysis. We have previously demonstrated that the prokaryotic zinc finger protein Ros87 shows a bipartite folding/unfolding process in which a metal binding intermediate converts to the native structure through a delicate barrier-less downhill transition. Significant variation in folding scenarios can be detected within protein families with high sequence identity and very similar folds and for the same sequence by varying conditions. For this reason, we here show, by means of DSC, CD and NMR, that also in different pH and ionic strength conditions Ros87 retains its partly downhill folding scenario demonstrating that, at least in metallo-proteins, the downhill mechanism can be found under a much wider range of conditions and coupled to other different transitions. We also show that mutations of Ros87 zinc coordination sphere produces a different folding scenario demonstrating that the organization of the metal ion core is determinant in the folding process of this family of proteins.

A Two-Ended Data-Driven Accelerated Sampling Method for Exploring the Transition Pathways between Two Known States of Protein

Conformational transitions of protein between different states are often associated with their biological functions. These dynamic processes, however, are usually not easy to be well characterized by experimental measurements, mainly because of inadequate temporal and spatial resolution. Meantime, sampling of configuration space with molecular dynamics (MD) simulations is still a challenge. Here we proposed a robust two-ended data-driven accelerated (teDA2) conformational sampling method, which drives the structural change in an adaptively updated feature space without introducing a bias potential. teDA2 was applied to explore adenylate kinase (ADK), a model with well characterized "open" and "closed" states. A single conformational transition event of ADK could be achieved within only a few or tens of nanoseconds sampled with teDA2. By analyzing hundreds of transition events, we reproduced different mechanisms and the associated pathways for domain motion of ADK reported in the literature. The multiroute characteristic of ADK was confirmed by the fact that some metastable states identified with teDA2 resemble available crystal structures determined at different conditions. This feature was further validated with Markov state modeling with independent MD simulations. Therefore, our work provides strong evidence for the conformational plasticity of protein, which is mainly due to the inherent degree of flexibility. As a reliable and efficient enhanced sampling protocol, teDA2 could be used to study the dynamics between functional states of various biomolecular machines.

Conformational Flexibility Drives Cold Adaptation in Pseudoalteromonas haloplanktis TAC125 Globins

Significance: Temperature is one of the most important drivers in shaping protein adaptations. Many biochemical and physiological processes are influenced by temperature. Proteins and enzymes from organisms living at low temperature are less stable in comparison to high-temperature adapted proteins. The lower stability is generally due to greater conformational flexibility. Recent Advances: Adaptive changes in the structure of cold-adapted proteins may occur at subunit interfaces, distant from the active site, thus producing energy changes associated with conformational transitions transmitted to the active site by allosteric modulation, valid also for monomeric proteins in which tertiary structural changes may play an essential role. Critical Issues: Despite efforts, the current experimental and computational methods still fail to produce general principles on protein evolution, since many changes are protein and species dependent. Environmental constraints or other biological cellular signals may override the ancestral information included in the structure of the protein, thus introducing inaccuracy in estimates and predictions on the evolutionary adaptations of proteins in response to cold adaptation. Future Directions: In this review, we describe the studies and approaches used to investigate stability and flexibility in the cold-adapted globins of the Antarctic marine bacterium Pseudoalteromonas haloplanktis TAC125. In fact, future research directions will be prescient on more detailed investigation of cold-adapted proteins and the role of fluctuations between different conformational states.

Protein tolerance to random circular permutation correlates with thermostability and local energetics of residue-residue contacts

Adenylate kinase (AK) orthologs with a range of thermostabilities were subjected to random circular permutation, and deep mutational scanning was used to evaluate where new protein termini were nondisruptive to activity. The fraction of circularly permuted variants that retained function in each library correlated with AK thermostability. In addition, analysis of the positional tolerance to new termini, which increase local conformational flexibility, showed that bonds were either functionally sensitive to cleavage across all homologs, differentially sensitive, or uniformly tolerant. The mobile AMP-binding domain, which displays the highest calculated contact energies, presented the greatest tolerance to new termini across all AKs. In contrast, retention of function in the lid and core domains was more dependent upon AK melting temperature. These results show that family permutation profiling identifies primary structure that has been selected by evolution for dynamics that are critical to activity within an enzyme family. These findings also illustrate how deep mutational scanning can be applied to protein homologs in parallel to differentiate how topology, stability, and local energetics govern mutational tolerance.

Folding of copper proteins: role of the metal?

Copper is a redox-active transition metal ion required for the function of many essential human proteins. For biosynthesis of proteins coordinating copper, the metal may bind before, during or after folding of the polypeptide. If the metal binds to unfolded or partially folded structures of the protein, such coordination may modulate the folding reaction. The molecular understanding of how copper is incorporated into proteins requires descriptions of chemical, thermodynamic, kinetic and structural parameters involved in the formation of protein-metal complexes. Because free copper ions are toxic, living systems have elaborate copper-transport systems that include particular proteins that facilitate efficient and specific delivery of copper ions to target proteins. Therefore, these pathways become an integral part of copper protein folding in vivo. This review summarizes biophysical-molecular in vitro work assessing the role of copper in folding and stability of copper-binding proteins as well as protein-protein copper exchange reactions between human copper transport proteins. We also describe some recent findings about the participation of copper ions and copper proteins in protein misfolding and aggregation reactions in vitro.

Nucleation of an Activating Conformational Change by a Cation-π Interaction

As a key molecule in biology, adenosine triphosphate (ATP) has numerous crucial functions in, for instance, energetics, post-translational modifications, nucleotide biosynthesis, and cofactor metabolism. Here, we have discovered an intricate interplay between the enzyme adenylate kinase and its substrate ATP. The side chain of an arginine residue was found to be an efficient sensor of the aromatic moiety of ATP through the formation of a strong cation-π interaction. In addition to recognition, the interaction was found to have dual functionality. First, it nucleates the activating conformational transition of the ATP binding domain and also affects the specificity in the distant AMP binding domain. In light of the functional consequences resulting from the cation-π interaction, it is possible that the mode of ATP recognition may be a useful tool in enzyme design.

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.

Biophysical and flavonoid-binding studies of the G protein ectodomain of group A human respiratory syncytial virus

The human Respiratory Syncytial Virus (hRSV) is the major causative agent of lower respiratory tract diseases in infants, young children and elderly. The membrane protein G is embedded in the viral lipid envelope and plays an adhesion function of the virus to host cells. The present study reports the production of the group A hRSV recombinant G protein ectodomain (edG) and its characterization of secondary structure and thermal unfolding by circular dichroism (CD), as well as the binding investigation of flavonoids quercetin and morin to this protein by fluorescent quenching. CD data reveal that edG is composed mostly of β-structure and its melting temperature is of 325 K. Fluorescence quenching experiments of hRSV edG show that the dissociation constants for the flavonoids binding are micromolar and the binding affinity for the edG/quercetin complex is inversely dependent on rising temperature while is directly dependent for the edG/morin interaction. The thermodynamic parameters suggest that hydrophobic contacts are important for the edG/morin association while van der Waals forces and hydrogen bonds contribute to the stabilization of the edG/quercetin complex. Thus, data reported herein may contribute to the development of new treatment strategies that prevent the viral infection by hRSV.

Structural and mutational analyses of psychrophilic and mesophilic adenylate kinases highlight the role of hydrophobic interactions in protein thermal stability

Protein thermal stability is an important field since thermally stable proteins are desirable in many academic and industrial settings. Information on protein thermal stabilization can be obtained by comparing homologous proteins from organisms living at distinct temperatures. Here, we report structural and mutational analyses of adenylate kinases (AKs) from psychrophilic Bacillus globisporus (AKp) and mesophilic Bacillus subtilis (AKm). Sequence and structural comparison showed suboptimal hydrophobic packing around Thr26 in the CORE domain of AKp, which was replaced with an Ile residue in AKm. Mutations that improved hydrophobicity of the Thr residue increased the thermal stability of the psychrophilic AKp, and the largest stabilization was observed for a Thr-to-Ile substitution. Furthermore, a reverse Ile-to-Thr mutation in the mesophilic AKm significantly decreased thermal stability. We determined the crystal structures of mutant AKs to confirm the impact of the residue substitutions on the overall stability. Taken together, our results provide a structural basis for the stability difference between psychrophilic and mesophilic AK homologues and highlight the role of hydrophobic interactions in protein thermal stability.

Engineering Allostery into Proteins

Our ability to engineer protein structure and function has grown dramatically over recent years. Perhaps the next level in protein design is to develop proteins whose function can be regulated in response to various stimuli, including ligand binding, pH changes, and light. Endeavors toward these goals have tested and expanded on our understanding of protein function and allosteric regulation. In this chapter, we provide examples from different methods for developing new allosterically regulated proteins. These methods range from whole insertion of regulatory domains into new host proteins, to covalent attachment of photoswitches to generate light-responsive proteins, and to targeted changes to specific amino acid residues, especially to residues identified to be important for relaying allosteric information across the protein framework. Many of the examples we discuss have already found practical use in medical and biotechnology applications.

Conformational dynamics and enzyme evolution

Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.

Dynamic allostery can drive cold adaptation in enzymes

Adaptation of organisms to environmental niches is a hallmark of evolution. One prevalent example is that of thermal adaptation, wherein two descendants evolve at different temperature extremes^1,2. Underlying the physiological differences between such organisms are changes in enzymes catalyzing essential reactions³, with orthologues from each organism undergoing adaptive mutations that preserve similar catalytic rates at their respective physiological temperatures ^4,5. The sequence changes responsible for these adaptive differences, however, are often at surface exposed sites distant from the substrate binding site, leaving the active site of the enzyme structurally unperturbed^6,7. How such changes are allosterically propagated to the active site, to modulate activity, is not known. Here we show that entropy-tuning changes can be engineered into distal sites of Escherichia coli adenylate kinase (AK) to quantitatively assess the role of dynamics in determining affinity, turnover, and the role in driving adaptation. The results not only reveal a dynamics-based allosteric tuning mechanism, but also uncover a spatial separation of the control of key enzymatic parameters. Fluctuations in one mobile domain (i.e. the LID) control substrate affinity, while dynamic attenuation in the other (i.e. the AMPbd) affects rate-limiting conformational changes governing enzyme turnover. Dynamics-based regulation may thus represent an elegant, widespread, and previously unrealized evolutionary adaptation mechanism that fine-tunes biological function without altering the ground state structure. Furthermore, because rigid-body conformational changes in both domains were thought to be rate limiting for turnover^8,9, these adaptation studies reveal a new paradigm for understanding the relationship between dynamics and turnover in AK.

The Molecular Origin of Enthalpy/Entropy Compensation in Biomolecular Recognition

Biomolecular recognition can be stubborn; changes in the structures of associating molecules, or the environments in which they associate, often yield compensating changes in enthalpies and entropies of binding and no net change in affinities. This phenomenon—termed enthalpy/entropy (H/S) compensation—hinders efforts in biomolecular design, and its incidence—often a surprise to experimentalists—makes interactions between biomolecules difficult to predict. Although characterizing H/S compensation requires experimental care, it is unquestionably a real phenomenon that has, from an engineering perspective, useful physical origins. Studying H/S compensation can help illuminate the still-murky roles of water and dynamics in biomolecular recognition and self-assembly. This review summarizes known sources of H/ S compensation (real and perceived) and lays out a conceptual framework for understanding and dissecting—and, perhaps, avoiding or exploiting—this phenomenon in biophysical systems.

Nonnative Energetic Frustrations in Protein Folding at Residual Level: A Simulation Study of Homologous Immunoglobulin-like β-Sandwich Proteins

Nonnative interactions cause energetic frustrations in protein folding and were found to dominate key events in folding intermediates. However, systematically characterizing energetic frustrations that are caused by nonnative intra-residue interactions at residual resolution is still lacking. Recently, we studied the folding of a set of homologous all-α proteins and found that nonnative-contact-based energetic frustrations are highly correlated to topology of the protein native-contact network. Here, we studied the folding of nine homologous immunoglobulin-like (Ig-like) β-sandwich proteins, and examined nonnative-contact-based energetic frustrations Gō-like model. Our calculations showed that nonnative-interaction-based energetic frustrations in β-sandwich proteins are much more complicated than those in all-α proteins, and they exhibit highly heterogeneous effects on the folding of secondary structures. Further, the nonnative interactions introduced distinct correlations in the folding of different folding-patches of β-sandwich proteins. Taken together, a strong interplay might exist between nonnative-interaction energetic frustrations and the protein native-contact networks, which ensures that β-sandwich domains adopt a common folding mechanism.

Molecular mechanism of ATP versus GTP selectivity of adenylate kinase

Enzymatic substrate selectivity is critical for the precise control of metabolic pathways. In cases where chemically related substrates are present inside cells, robust mechanisms of substrate selectivity are required. Here, we report the mechanism utilized for catalytic ATP versus GTP selectivity during adenylate kinase (Adk) -mediated phosphorylation of AMP. Using NMR spectroscopy we found that while Adk adopts a catalytically competent and closed structural state in complex with ATP, the enzyme is arrested in a catalytically inhibited and open state in complex with GTP. X-ray crystallography experiments revealed that the interaction interfaces supporting ATP and GTP recognition, in part, are mediated by coinciding residues. The mechanism provides an atomic view on how the cellular GTP pool is protected from Adk turnover, which is important because GTP has many specialized cellular functions. In further support of this mechanism, a structure-function analysis enabled by synthesis of ATP analogs suggests that a hydrogen bond between the adenine moiety and the backbone of the enzyme is vital for ATP selectivity. The importance of the hydrogen bond for substrate selectivity is likely general given the conservation of its location and orientation across the family of eukaryotic protein kinases.

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.

Multiple Pathways and Time Scales for Conformational Transitions in apo-Adenylate Kinase

The open/close transition in adenylate kinase (AK) is regarded as a representative example for large-scale conformational transition in proteins, yet its mechanism remains unclear despite numerous experimental and computational studies. Using extensive (∼50 μs) explicit solvent atomistic simulations and Markov state analysis, we shed new lights on the mechanism of this transition in the apo form of AK. The closed basin of apo AK features an open NMP domain while the LID domain closes and rotates toward it. Therefore, although the computed structural properties of the closed ensemble are consistent with previously reported FRET and PRE measurements, our simulations suggest that NMP closure is likely to follow AMP binding, in contrast to the previous interpretation of FRET and PRE data that the apo state was able to sample the fully closed conformation for "ligand selection". The closed state ensemble is found to be kinetically heterogeneous; multiple pathways and time scales are associated with the open/close transition, providing new clues to the disparate time scales observed in different experiments. Besides interdomain interactions, a novel mutual information analysis identifies specific intradomain interactions that correlate strongly to transition kinetics, supporting observations from previous chimera experiments. While our results underscore the role of internal domain properties in determining the kinetics of open/close transition in apo AK, no evidence is observed for any significant degree of local unfolding during the transition. These observations about AK have general implications to our view of conformational states, transition pathways, and time scales of conformational changes in proteins. The key features and time scales of observed transition pathways are robust and similar from simulations using two popular fixed charge force fields.

Effects of structurally stabilized EGF and bFGF on wound healing in type I and type II diabetic mice

Diabetes mellitus comprises a multiple metabolic disorder that affects millions of people worldwide and consequentially poses challenges for clinical treatment. Among the various complications, diabetic ulcer constitutes the most prevalent associated disorder and leads to delayed wound healing. To enhance wound healing capacity, we developed structurally stabilized epidermal growth factor (ST-EGF) and basic fibroblast growth factor (ST-bFGF) to overcome limitations of commercially available EGF (CA-EGF) and bFGF (CA-bFGF), such as short half-life and loss of activity after loading onto a matrix. Neither ST-EGF nor ST-bFGF was toxic, and both were more stable at higher temperatures than CA-EGF and CA-bFGF. We loaded ST-EGF and ST-bFGF onto a hyaluronate-collagen dressing (HCD) matrix, a biocompatible carrier, and tested the effectiveness of this system in promoting wound healing in a mouse model of diabetes. Wounds treated with HCD matrix loaded with 0.3 μg/cm² ST-EGF or 1 μg/cm² ST-bFGF showed a more rapid rate of tissue repair as compared to the control in type I and II diabetes models. Our results indicate that an HDC matrix loaded with 0.3 μg/cm² ST-EGF or 1 μg/cm² ST-bFGF can promote wound healing in diabetic ulcers and are suitable for use in wound dressings owing to their stability for long periods at room temperature.

STATEMENT OF SIGNIFICANCE

Various types of dressing materials loaded with growth factors, such as VEGF, EGF, and bFGF, are widely used to effect wound repair. However, such growth factor-loaded materials have several limitations for use as therapeutic agents in healing-impaired diabetic wounds. To overcome these limitations, we have developed new materials containing structurally stabilized EGF (ST-EGF) and bFGF (ST-bFGF). To confirm the wound healing capacity of newly developed materials (ST-EGF and ST-bFGF-loaded hyaluronate-collagen dressing [HCD] matrix), we applied these matrices in type I and type II diabetic wounds. Notably, these matrices were able to accelerate wound healing including re-epithelialization, neovascularization, and collagen deposition. Consequentially, these ST-EGF and ST-bFGF-loaded HCD matrix may be used as future therapeutic agents in patients with diabetic foot ulcers.

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.

The relationship between folding and activity in UreG, an intrinsically disordered enzyme

A growing body of literature on intrinsically disordered proteins (IDPs) led scientists to rethink the structure-function paradigm of protein folding. Enzymes are often considered an exception to the rule of intrinsic disorder (ID), believed to require a unique structure for catalysis. However, recent studies revealed the presence of disorder in several functional native enzymes. In the present work, we address the importance of dynamics for catalysis, by investigating the relationship between folding and activity in Sporosarcina pasteurii UreG (SpUreG), a P-loop GTPase and the first discovered native ID enzyme, involved in the maturation of the nickel-containing urease. The effect of denaturants and osmolytes on protein structure and activity was analyzed using circular dichroism (CD), Site-Directed Spin Labeling (SDSL) coupled to EPR spectroscopy, and enzymatic assays. Our data show that SpUreG needs a "flexibility window" to be catalytically competent, with both too low and too high mobility being detrimental for its activity.