Conformational heterogeneity within the LID domain mediates substrate binding to Escherichia coli adenylate kinase: function follows fluctuations

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.

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.

Alkaloid binding to opium poppy major latex proteins triggers structural modification and functional aggregation

Opium poppy accumulates copious amounts of several benzylisoquinoline alkaloids including morphine, noscapine, and papaverine, in the specialized cytoplasm of laticifers, which compose an internal secretory system associated with phloem throughout the plant. The contiguous latex includes an abundance of related proteins belonging to the pathogenesis-related (PR)10 family known collectively as major latex proteins (MLPs) and representing at least 35% of the total cellular protein content. Two latex MLP/PR10 proteins, thebaine synthase and neopione isomerase, have recently been shown to catalyze late steps in morphine biosynthesis previously assigned as spontaneous reactions. Using a combination of sucrose density-gradient fractionation-coupled proteomics, differential scanning fluorimetry, isothermal titration calorimetry, and X-ray crystallography, we show that the major latex proteins are a family of alkaloid-binding proteins that display altered conformation in the presence of certain ligands. Addition of MLP/PR10 proteins to yeast strains engineered with morphine biosynthetic genes from the plant significantly enhanced the conversion of salutaridine to morphinan alkaloids.

Explore Protein Conformational Space With Variational Autoencoder

Molecular dynamics (MD) simulations have been actively used in the study of protein structure and function. However, extensive sampling in the protein conformational space requires large computational resources and takes a prohibitive amount of time. In this study, we demonstrated that variational autoencoders (VAEs), a type of deep learning model, can be employed to explore the conformational space of a protein through MD simulations. VAEs are shown to be superior to autoencoders (AEs) through a benchmark study, with low deviation between the training and decoded conformations. Moreover, we show that the learned latent space in the VAE can be used to generate unsampled protein conformations. Additional simulations starting from these generated conformations accelerated the sampling process and explored hidden spaces in the conformational landscape.

Integrated structural and evolutionary analysis reveals common mechanisms underlying adaptive evolution in mammals

Understanding the molecular basis of adaptation to the environment is a central question in evolutionary biology, yet linking detected signatures of positive selection to molecular mechanisms remains challenging. Here we demonstrate that combining sequence-based phylogenetic methods with structural information assists in making such mechanistic interpretations on a genomic scale. Our integrative analysis shows that positively selected sites tend to colocalize on protein structures and that positively selected clusters are found in functionally important regions of proteins, indicating that positive selection can contravene the well-known principle of evolutionary conservation of functionally important regions. This unexpected finding, along with our discovery that positive selection acts on structural clusters, opens previously unexplored strategies for the development of better models of protein evolution. Remarkably, proteins where we detect the strongest evidence of clustering belong to just two functional groups: Components of immune response and metabolic enzymes. This gives a coherent picture of pathogens and xenobiotics as important drivers of adaptive evolution of mammals.

Stabilization of Active Site Dynamics Leads to Increased Activity with 3'-Azido-3'-deoxythymidine Monophosphate for F105Y Mutant Human Thymidylate Kinase

Thymidylate kinases are essential enzymes with roles in DNA synthesis and repair and have been the target of drug development for antimalarials, antifungals, HIV treatment, and cancer therapeutics. Human thymidylate kinase (hTMPK) conversion of the anti-HIV prodrug 3'-azido-3'-deoxythymidine (AZT or zidovudine) monophosphate to diphosphate is the rate-limiting step in the activation of AZT. A point mutant (F105Y) has been previously reported with significantly increased activity for the monophosphate form of the drug [3'-azidothymidine-5'-monophosphate (AZTMP)]. Using solution nuclear magnetic resonance (NMR) techniques, we show that while the wild-type (WT) and F105Y hTMPK adopt the same structure in solution, significant changes in dynamics may explain their different activities toward TMP and AZTMP. ¹³C spin-relaxation measurements show that there is little change in dynamics on the ps to ns time scale. In contrast, methyl ¹H relaxation dispersion shows that AZTMP alters adenosine nucleotide handling in the WT protein but not in the mutant. Additionally, the F105Y mutant has reduced conformational flexibility, leading to an increase in affinity for the product ADP and a slower rate of phosphorylation of TMP. The dynamics at the catalytic center for F105Y bound to AZTMP are tuned to the same frequency as WT bound to TMP, which may explain the mutant's catalytic efficiency toward the prodrug.

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.

Structural analyses of adenylate kinases from Antarctic and tropical fishes for understanding cold adaptation of enzymes

Psychrophiles are extremophilic organisms capable of thriving in cold environments. Proteins from these cold-adapted organisms can remain physiologically functional at low temperatures, but are structurally unstable even at moderate temperatures. Here, we report the crystal structure of adenylate kinase (AK) from the Antarctic fish Notothenia coriiceps, and identify the structural basis of cold adaptation by comparison with homologues from tropical fishes including Danio rerio. The structure of N. coriiceps AK (AKNc) revealed suboptimal hydrophobic packing around three Val residues in its central CORE domain, which are replaced with Ile residues in D. rerio AK (AKDr). The Val-to-Ile mutations that improve hydrophobic CORE packing in AKNc increased stability at high temperatures but decreased activity at low temperatures, suggesting that the suboptimal hydrophobic CORE packing is important for cold adaptation. Such linkage between stability and activity was also observed in AKDr. Ile-to-Val mutations that destabilized the tropical AK resulted in increased activity at low temperatures. Our results provide the structural basis of cold adaptation of a psychrophilic enzyme from a multicellular, eukaryotic organism, and highlight the similarities and differences in the structural adjustment of vertebrate and bacterial psychrophilic AKs during cold adaptation.

The phosphocarrier protein HPr of the bacterial phosphotransferase system globally regulates energy metabolism by directly interacting with multiple enzymes in Escherichia coli

The histidine-phosphorylatable phosphocarrier protein (HPr) is an essential component of the sugar-transporting phosphotransferase system (PTS) in many bacteria. Recent interactome findings suggested that HPr interacts with several carbohydrate-metabolizing enzymes, but whether HPr plays a regulatory role was unclear. Here, we provide evidence that HPr interacts with a large number of proteins in Escherichia coli We demonstrate HPr-dependent allosteric regulation of the activities of pyruvate kinase (PykF, but not PykA), phosphofructokinase (PfkB, but not PfkA), glucosamine-6-phosphate deaminase (NagB), and adenylate kinase (Adk). HPr is either phosphorylated on a histidyl residue (HPr-P) or non-phosphorylated (HPr). PykF is activated only by non-phosphorylated HPr, which decreases the PykF K_half for phosphoenolpyruvate by 10-fold (from 3.5 to 0.36 mm), thus influencing glycolysis. PfkB activation by HPr, but not by HPr-P, resulted from a decrease in the K_half for fructose-6-P, which likely influences both gluconeogenesis and glycolysis. Moreover, NagB activation by HPr was important for the utilization of amino sugars, and allosteric inhibition of Adk activity by HPr-P, but not by HPr, allows HPr to regulate the cellular energy charge coordinately with glycolysis. These observations suggest that HPr serves as a directly interacting global regulator of carbon and energy metabolism and probably of other physiological processes in enteric bacteria.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

WEBnm@ v2.0: Web server and services for comparing protein flexibility

BACKGROUND

Normal mode analysis (NMA) using elastic network models is a reliable and cost-effective computational method to characterise protein flexibility and by extension, their dynamics. Further insight into the dynamics-function relationship can be gained by comparing protein motions between protein homologs and functional classifications. This can be achieved by comparing normal modes obtained from sets of evolutionary related proteins.

RESULTS

We have developed an automated tool for comparative NMA of a set of pre-aligned protein structures. The user can submit a sequence alignment in the FASTA format and the corresponding coordinate files in the Protein Data Bank (PDB) format. The computed normalised squared atomic fluctuations and atomic deformation energies of the submitted structures can be easily compared on graphs provided by the web user interface. The web server provides pairwise comparison of the dynamics of all proteins included in the submitted set using two measures: the Root Mean Squared Inner Product and the Bhattacharyya Coefficient. The Comparative Analysis has been implemented on our web server for NMA, WEBnm@, which also provides recently upgraded functionality for NMA of single protein structures. This includes new visualisations of protein motion, visualisation of inter-residue correlations and the analysis of conformational change using the overlap analysis. In addition, programmatic access to WEBnm@ is now available through a SOAP-based web service. Webnm@ is available at http://apps.cbu.uib.no/webnma .

CONCLUSION

WEBnm@ v2.0 is an online tool offering unique capability for comparative NMA on multiple protein structures. Along with a convenient web interface, powerful computing resources, and several methods for mode analyses, WEBnm@ facilitates the assessment of protein flexibility within protein families and superfamilies. These analyses can give a good view of how the structures move and how the flexibility is conserved over the different structures.