Substrate Binding Specifically Modulates Domain Arrangements in Adenylate Kinase

Predicting binding events in very flexible, allosteric, multi-domain proteins

Knowledge of the structures formed by proteins and small ligands is of fundamental importance for understanding molecular principles of chemotherapy and for designing new and more effective drugs. Due to the still high costs and to the several limitations of experimental techniques, it is most often desirable to predict these ligand-protein complexes in silico, particularly when screening for new putative drugs from databases of millions of compounds. While virtual screening based on molecular docking is widely used for this purpose, it generally fails in mimicking binding events associated with large conformational changes in the protein, particularly when the latter involve multiple domains. In this work, we describe a new methodology aimed at generating bound-like conformations of very flexible and allosteric proteins bearing multiple binding sites. Validation was performed on the enzyme adenylate kinase (ADK), a paradigmatic example of proteins that undergo very large conformational changes upon ligand binding. By only exploiting the unbound structure and the putative binding sites of the protein, we generated a significant fraction of bound-like structures, which employed in ensemble-docking calculations allowed to find native-like poses of substrates, inhibitors, and catalytically incompetent binders. Our protocol provides a general framework for the generation of bound-like conformations of flexible proteins that are suitable to host different ligands, demonstrating high sensitivity to the fine chemical details that regulate protein's activity. We foresee applications in virtual screening for difficult targets, prediction of the impact of amino acid mutations on structure and dynamics, and protein engineering.

Tackling Hysteresis in Conformational Sampling: How to Be Forgetful with MEMENTO

The structure of proteins has long been recognized to hold the key to understanding and engineering their function, and rapid advances in structural biology and protein structure prediction are now supplying researchers with an ever-increasing wealth of structural information. Most of the time, however, structures can only be determined in free energy minima, one at a time. While conformational flexibility may thus be inferred from static end-state structures, their interconversion mechanisms─a central ambition of structural biology─are often beyond the scope of direct experimentation. Given the dynamical nature of the processes in question, many studies have attempted to explore conformational transitions using molecular dynamics (MD). However, ensuring proper convergence and reversibility in the predicted transitions is extremely challenging. In particular, a commonly used technique to map out a path from a starting to a target conformation called steered MD (SMD) can suffer from starting-state dependence (hysteresis) when combined with techniques such as umbrella sampling (US) to compute the free energy profile of a transition. Here, we study this problem in detail on conformational changes of increasing complexity. We also present a new, history-independent approach that we term "MEMENTO" (Morphing End states by Modelling Ensembles with iNdependent TOpologies) to generate paths that alleviate hysteresis in the construction of conformational free energy profiles. MEMENTO utilizes template-based structure modelling to restore physically reasonable protein conformations based on coordinate interpolation (morphing) as an ensemble of plausible intermediates, from which a smooth path is picked. We compare SMD and MEMENTO on well-characterized test cases (the toy peptide deca-alanine and the enzyme adenylate kinase) before discussing its use in more complicated systems (the kinase P38α and the bacterial leucine transporter LeuT). Our work shows that for all but the simplest systems SMD paths should not in general be used to seed umbrella sampling or related techniques, unless the paths are validated by consistent results from biased runs in opposite directions. MEMENTO, on the other hand, performs well as a flexible tool to generate intermediate structures for umbrella sampling. We also demonstrate that extended end-state sampling combined with MEMENTO can aid the discovery of collective variables on a case-by-case basis.

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.

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.

Structural dynamics support electrostatic interactions in the active site of Adenylate Kinase

Electrostatic preorganization as well as structural and dynamic heterogeneity are often used to rationalize the remarkable catalytic efficiency of enzymes. However, they are often presented as incompatible because the generation of permanent electrostatic effects implies that the protein structure remains rigid. Here, we use a metric, electric fields, that can treat electrostatic contributions and dynamics effects on equal footing, for a unique perspective on enzymatic catalysis. We find that the residues that contribute the most to electrostatic interactions with the substrate in the active site of Adenylate Kinase (our working example) are also the most flexible residues. Further, entropy-tuning mutations raise flexibility at the picosecond timescale where more conformations can be visited on short time periods, thereby softening the sharp heterogeneity normally visible at the microsecond timescale.

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.

Adenylate kinases of thermophiles Aquifex aeolicus and Geobacillus stearothermophilus: biochemical and kinetic studies

Adenylate kinases (AK) play a pivotal role in the regulation of cellular energy. The aim of our work was to achieve the overproduction and purification of AKs from two groups of bacteria and to determine, for the first time, the comprehensive biochemical and kinetic properties of adenylate kinase from Gram-negative Aquifex aeolicus (AK_aq) and Gram-positive Geobacillus stearothermophilus (AK_st). Therefore we determined K_M and V_max values, and the effects of temperature, pH, metal ions, donors of the phosphate groups and inhibitor Ap₅A for both thermophilic AKs. The kinetic studies indicate that both AKs exhibit significantly higher affinity for substrates with the pyrophosphate group than for adenosine monophosphate. AK activation by Mg²⁺ and Mn²⁺ revealed that both ions are efficient in the synthesis of adenosine diphosphate and adenosine triphosphate; however, Mn²⁺ ions at 0.2-2.0 mmol/L concentration were more efficient in the activation of the ATP synthesis than Mg²⁺ ions. Our research demonstrates that zinc ions inhibit the activity of enzymes in both directions, while Ap₅A at a concentration of 10 µmol/L and 50 µmol/L inhibited both enzymes with a different efficiency. Sigmoid-like kinetics were detected at high ATP concentrations not balanced by Mg²⁺, suggesting the allosteric effect of ATP for both bacterial AKs.

Understanding enzyme behavior in a crowded scenario through modulation in activity, conformation and dynamics

Macromolecular crowding, inside the physiological interior, modulates the energy landscape of biological macromolecules in multiple ways. Amongst these, enzymes occupy a special place and hence understanding the function of the same in the crowded interior is of utmost importance. In this study, we have investigated the manner in which the multidomain enzyme, AK3L1 (PDB ID: 1ZD8), an isoform of adenylate kinase, has its features affected in presence of commonly used crowders (PEG 8, Dextran 40, Dextran 70, and Ficoll 70). Michaelis Menten plots reveal that the crowders in general enhance the activity of the enzyme, with the K_m and V_max values showing significant variations. Ficoll 70, induced the maximum activity for AK3L1 at 100 g/L, beyond which the activity reduced. Ensemble FRET studies were performed to provide insights into the relative domain (LID and CORE) displacements in presence of the crowders. Solvation studies reveal that the protein matrix surrounding the probe CPM (7-diethylamino-3-(4-maleimido-phenyl)-4-methylcoumarin) gets restricted in presence of the crowders, with Ficoll 70 providing the maximum rigidity, the same being linked to the decrease in the activity of the enzyme. Through our multipronged approach, we have observed a distinct correlation between domain displacement, enzyme activity and associated dynamics. Thus, keeping in mind the complex nature of enzyme activity and the surrounding bath of dense soup that the biological entity remains immersed in, indeed more such approaches need to be undertaken to have a better grasp of the "enzymes in the crowd".

Dynamic Connection between Enzymatic Catalysis and Collective Protein Motions

Enzymes employ a wide range of protein motions to achieve efficient catalysis of chemical reactions. While the role of collective protein motions in substrate binding, product release, and regulation of enzymatic activity is generally understood, their roles in catalytic steps per se remain uncertain. Here, molecular dynamics simulations, enzyme kinetics, X-ray crystallography, and nuclear magnetic resonance spectroscopy are combined to elucidate the catalytic mechanism of adenylate kinase and to delineate the roles of catalytic residues in catalysis and the conformational change in the enzyme. This study reveals that the motions in the active site, which occur on a time scale of picoseconds to nanoseconds, link the catalytic reaction to the slow conformational dynamics of the enzyme by modulating the free energy landscapes of subdomain motions. In particular, substantial conformational rearrangement occurs in the active site following the catalytic reaction. This rearrangement not only affects the reaction barrier but also promotes a more open conformation of the enzyme after the reaction, which then results in an accelerated opening of the enzyme compared to that of the reactant state. The results illustrate a linkage between enzymatic catalysis and collective protein motions, whereby the disparate time scales between the two processes are bridged by a cascade of intermediate-scale motion of catalytic residues modulating the free energy landscapes of the catalytic and conformational change processes.

Assessing the Interactions of Statins with Human Adenylate Kinase Isoenzyme 1: Fluorescence and Enzyme Kinetic Studies

Statins are the most effective cholesterol-lowering drugs. They also exert many pleiotropic effects, including anti-cancer and cardio- and neuro-protective. Numerous nano-sized drug delivery systems were developed to enhance the therapeutic potential of statins. Studies on possible interactions between statins and human proteins could provide a deeper insight into the pleiotropic and adverse effects of these drugs. Adenylate kinase (AK) was found to regulate HDL endocytosis, cellular metabolism, cardiovascular function and neurodegeneration. In this work, we investigated interactions between human adenylate kinase isoenzyme 1 (hAK1) and atorvastatin (AVS), fluvastatin (FVS), pravastatin (PVS), rosuvastatin (RVS) and simvastatin (SVS) with fluorescence spectroscopy. The tested statins quenched the intrinsic fluorescence of hAK1 by creating stable hAK1-statin complexes with the binding constants of the order of 10⁴ M⁻¹. The enzyme kinetic studies revealed that statins inhibited hAK1 with significantly different efficiencies, in a noncompetitive manner. Simvastatin inhibited hAK1 with the highest yield comparable to that reported for diadenosine pentaphosphate, the only known hAK1 inhibitor. The determined AK sensitivity to statins differed markedly between short and long type AKs, suggesting an essential role of the LID domain in the AK inhibition. Our studies might open new horizons for the development of new modulators of short type AKs.

Coarse-Grained Modeling of Multiple Pathways in Conformational Transitions of Multi-Domain Proteins

Large-scale conformational transitions in multi-domain proteins are often essential for their functions. To investigate the transitions, it is necessary to explore multiple potential pathways, which involve different intermediate structures. Here, we present a multi-basin (MB) coarse-grained (CG) structure-based Go̅ model for describing transitions in proteins with more than two moving domains. This model is an extension of our dual-basin Go̅ model in which system-dependent parameters are determined systematically using the multistate Bennett acceptance ratio method. In the MB Go̅ model for multi-domain proteins, we assume that intermediate structures may have partial inter-domain native contacts. This approach allows us to search multiple transition pathways that involve distinct intermediate structures using the CG molecular dynamics (MD) simulations. We apply this scheme to an enzyme, adenylate kinase (AdK), which has three major domains and can move along two different pathways. Using the optimized mixing parameters for each pathway, AdK shows frequent transitions between the Open, Closed, and the intermediate basins and samples a wide variety of conformations within each basin. The explored multiple transition pathways could be compared with experimental data and examined in more detail by atomistic MD simulations.

Increasing the Sampling Efficiency of Protein Conformational Change by Combining a Modified Replica Exchange Molecular Dynamics and Normal Mode Analysis

Understanding conformational change at an atomic level is significant when determining a protein functional mechanism. Replica exchange molecular dynamics (REMD) is a widely used enhanced sampling method to explore protein conformational space. However, REMD with an explicit solvent model requires huge computational resources, immensely limiting its application. In this study, a variation of parallel tempering metadynamics (PTMetaD) with the omission of solvent-solvent interactions in exchange attempts and the use of low-frequency modes calculated by normal-mode analysis (NMA) as collective variables (CVs), namely ossPTMetaD, is proposed with the aim to accelerate MD simulations simultaneously in temperature and geometrical spaces. For testing the performance of ossPTMetaD, five protein systems with diverse biological functions and motion patterns were selected, including large-scale domain motion (AdK), flap movement (HIV-1 protease and BACE1), and DFG-motif flip in kinases (p38α and c-Abl). The simulation results showed that ossPTMetaD requires much fewer numbers of replicas than temperature REMD (T-REMD) with a reduction of ∼70% to achieve a similar exchange ratio. Although it does not obey the detailed balance condition, ossPTMetaD provides consistent results with T-REMD and experimental data. The high accessibility of the large conformational change of protein systems by ossPTMetaD, especially in simulating the very challenging DFG-motif flip of protein kinases, demonstrated its high efficiency and robustness in the characterization of the large-scale protein conformational change pathway and associated free energy profile.

Hypoxanthine-Guanine Phosphoribosyltransferase/adenylate Kinase From Zobellia galactanivorans: A Bifunctional Catalyst for the Synthesis of Nucleoside-5'-Mono-, Di- and Triphosphates

In our search for novel biocatalysts for the synthesis of nucleic acid derivatives, we found a good candidate in a putative dual-domain hypoxanthine-guanine phosphoribosyltransferase (HGPRT)/adenylate kinase (AMPK) from Zobellia galactanivorans (ZgHGPRT/AMPK). In this respect, we report for the first time the recombinant expression, production, and characterization of a bifunctional HGPRT/AMPK. Biochemical characterization of the recombinant protein indicates that the enzyme is a homodimer, with high activity in the pH range 6-7 and in a temperature interval from 30 to 80°C. Thermal denaturation experiments revealed that ZgHGPRT/AMPK exhibits an apparent unfolding temperature (Tm) of 45°C and a retained activity of around 80% when incubated at 40°C for 240 min. This bifunctional enzyme shows a dependence on divalent cations, with a remarkable preference for Mg²⁺ and Co²⁺ as cofactors. More interestingly, substrate specificity studies revealed ZgHGPRT/AMPK as a bifunctional enzyme, which acts as phosphoribosyltransferase or adenylate kinase depending upon the nature of the substrate. Finally, to assess the potential of ZgHGPRT/AMPK as biocatalyst for the synthesis of nucleoside-5′-mono, di- and triphosphates, the kinetic analysis of both activities (phosphoribosyltransferase and adenylate kinase) and the effect of water-miscible solvents on enzyme activity were studied.

Exploring Conformational Change of Adenylate Kinase by Replica Exchange Molecular Dynamic Simulation

Replica exchange molecular dynamics (REMD) simulation is a popular enhanced sampling method that is widely used for exploring the atomic mechanism of protein conformational change. However, the requirement of huge computational resources for REMD, especially with the explicit solvent model, largely limits its application. In this study, the availability and efficiency of a variant of velocity-scaling REMD (vsREMD) was assessed with adenylate kinase as an example. Although vsREMD achieved results consistent with those from conventional REMD and experimental studies, the number of replicas required for vsREMD (30) was much less than that for conventional REMD (80) to achieve a similar acceptance rate (∼0.2), demonstrating high efficiency of vsREMD to characterize the protein conformational change and associated free-energy profile. Thus, vsREMD is a highly efficient approach for studying the large-scale conformational change of protein systems.

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.