Guanylate kinase, induced fit, and the allosteric spring probe

Enzyme Engineering by Force: DNA Springs for the Modulation of Biocatalytic Trajectories

The engineering of enzymatic activity generally involves alteration of the protein primary sequences, which introduce structural changes that give rise to functional improvements. Mechanical forces have been used to interrogate protein biophysics, leading to deep mechanistic insights in single-molecule studies. Here, we use simple DNA springs to apply small pulling forces to perturb the active site of a thermostable alcohol dehydrogenase. Methods were developed to enable the study of different spring lengths and spring orientations under bulk catalysis conditions. Tension applied across the active site expanded the binding pocket volume and shifted the preference of the enzyme for longer chain-length substrates, which could be tuned by altering the spring length and the resultant applied force. The substrate specificity changes did not occur when the DNA spring was either severed or rotated by ∼90°. These findings demonstrate an alternative approach in protein engineering, where active site architectures can be dynamically and reversibly remodeled using applied mechanical forces.

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.

Biotechnology applications of proteins functionalized with DNA oligonucleotides

The functionalization of proteins with DNA through the formation of covalent bonds enables a wide range of biotechnology advancements. For example, single-molecule analytical methods rely on bioconjugated DNA as elastic biolinkers for protein immobilization. Labeling proteins with DNA enables facile protein identification, as well as spatial and temporal organization and control of protein within DNA-protein networks. Bioconjugation reactions can target native, engineered, and non-canonical amino acids (NCAAs) within proteins. In addition, further protein engineering via the incorporation of peptide tags and self-labeling proteins can also be used for conjugation reactions. The selection of techniques will depend on application requirements such as yield, selectivity, conjugation position, potential for steric hindrance, cost, commercial availability, and potential impact on protein function.

Solution structure and functional investigation of human guanylate kinase reveals allosteric networking and a crucial role for the enzyme in cancer

Human guanylate kinase (hGMPK) is the only known enzyme responsible for cellular GDP production, making it essential for cellular viability and proliferation. Moreover, hGMPK has been assigned a critical role in metabolic activation of antiviral and antineoplastic nucleoside-analog prodrugs. Given that hGMPK is indispensable for producing the nucleotide building blocks of DNA, RNA, and cGMP and that cancer cells possess elevated GTP levels, it is surprising that a detailed structural and functional characterization of hGMPK is lacking. Here, we present the first high-resolution structure of hGMPK in the apo form, determined with NMR spectroscopy. The structure revealed that hGMPK consists of three distinct regions designated as the LID, GMP-binding (GMP-BD), and CORE domains and is in an open configuration that is nucleotide binding-competent. We also demonstrate that nonsynonymous single-nucleotide variants (nsSNVs) of the hGMPK CORE domain distant from the nucleotide-binding site of this domain modulate enzymatic activity without significantly affecting hGMPK's structure. Finally, we show that knocking down the hGMPK gene in lung adenocarcinoma cell lines decreases cellular viability, proliferation, and clonogenic potential while not altering the proliferation of immortalized, noncancerous human peripheral airway cells. Taken together, our results provide an important step toward establishing hGMPK as a potential biomolecular target, from both an orthosteric (ligand-binding sites) and allosteric (location of CORE domain-located nsSNVs) standpoint.

Photo-electrochemical Bioanalysis of Guanosine Monophosphate Using Coupled Enzymatic Reactions at a CdS/ZnS Quantum Dot Electrode

A photo-electrochemical sensor for the specific detection of guanosine monophosphate (GMP) is demonstrated, based on three enzymes combined in a coupled reaction assay. The first reaction involves the adenosine triphosphate (ATP)-dependent conversion of GMP to guanosine diphosphate (GDP) by guanylate kinase, which warrants substrate specificity. The reaction products ADP and GDPare co-substrates for the enzymatic conversion of phosphoenolpyruvate to pyruvate in a second reaction mediated by pyruvate kinase. Pyruvate in turn is the co-substrate for lactate dehydrogenase that generates lactate via oxidation of nicotinamide adenine dinucleotide (reduced form) NADH to NAD(+). This third enzymatic reaction is electrochemically detected. For this purpose a CdS/ZnS quantum dot (QD) electrode is illuminated and the photocurrent response under fixed potential conditions is evaluated. The sequential enzyme reactions are first evaluated in solution. Subsequently, a sensor for GMP is constructed using polyelectrolytes for enzyme immobilization.

Insights into open/closed conformations of the catalytically active human guanylate kinase as investigated by small-angle X-ray scattering

Abstract

Bio-catalysis is the outcome of a subtle interplay between internal motions in enzymes and chemical kinetics. Small-angle X-ray scattering (SAXS) investigation of an enzyme’s internal motions during catalysis offers an integral view of the protein’s structural plasticity, dynamics, and function, which is useful for understanding allosteric effects and developing novel medicines. Guanylate kinase (GMPK) is an essential enzyme involved in the guanine nucleotide metabolism of unicellular and multicellular organisms. It is also required for the intracellular activation of numerous antiviral and anticancer purine nucleoside analog prodrugs. Catalytically active recombinant human GMPK (hGMPK) was purified for the first time and changes in the size and shape of open/closed hGMPK were tracked by SAXS. The binding of substrates (GMP + AMPPNP or Ap5G or GMP + ADP) resulted in the compaction of size and shape of hGMPK. The structural changes between open and completely closed hGMPK conformation were confirmed by observing differences in the hGMPK secondary structures with circular dichroism spectroscopy.

Graphical abstract

Electronic supplementary material

The online version of this article (doi:10.1007/s00249-015-1079-9) contains supplementary material, which is available to authorized users.

Mechanistic insight into the functional transition of the enzyme guanylate kinase induced by a single mutation

Dramatic functional changes of enzyme usually require scores of alterations in amino acid sequence. However, in the case of guanylate kinase (GK), the functional novelty is induced by a single (S→P) mutation, leading to the functional transition of the enzyme from a phosphoryl transfer kinase into a phosphorprotein interaction domain. Here, by using molecular dynamic (MD) and metadynamics simulations, we provide a comprehensive description of the conformational transitions of the enzyme after mutating serine to proline. Our results suggest that the serine plays a crucial role in maintaining the closed conformation of wild-type GK and the GMP recognition. On the contrary, the S→P mutant exhibits a stable open conformation and loses the ability of ligand binding, which explains its functional transition from the GK enzyme to the GK domain. Furthermore, the free energy profiles (FEPs) obtained by metadymanics clearly demonstrate that the open-closed conformational transition in WT GK is positive correlated with the process of GMP binding, indicating the GMP-induced closing motion of GK enzyme, which is not observed in the mutant. In addition, the FEPs show that the S→P mutation can also leads to the mis-recognition of GMP, explaining the vanishing of catalytic activity of the mutant.

Dissipative dynamics of enzymes

We explore enzyme conformational dynamics at sub-Å resolution, specifically, temperature effects. The ensemble-averaged mechanical response of the folded enzyme is viscoelastic in the whole temperature range between the warm and cold denaturation transitions. The dissipation parameter γ of the viscoelastic description decreases by a factor of 2 as the temperature is raised from 10 to 45 °C; the elastic parameter K shows a similar decrease. Thus, when probed dynamically, the enzyme softens for increasing temperature. Equilibrium mechanical experiments with the DNA spring (and a different enzyme) also show, qualitatively, a small softening for increasing temperature.

Protein mechanics: how force regulates molecular function

BACKGROUND

Regulation of proteins is ubiquitous and vital for any organism. Protein activity can be altered chemically, by covalent modifications or non-covalent binding of co-factors. Mechanical forces are emerging as an additional way of regulating proteins, by inducing a conformational change or by partial unfolding.

SCOPE

We review some advances in experimental and theoretical techniques to study protein allostery driven by mechanical forces, as opposed to the more conventional ligand driven allostery. In this respect, we discuss recent single molecule pulling experiments as they have substantially augmented our view on the protein allostery by mechanical signals in recent years. Finally, we present a computational analysis technique, Force Distribution Analysis, that we developed to reveal allosteric pathways in proteins.

MAJOR CONCLUSIONS

Any kind of external perturbation, being it ligand binding or mechanical stretching, can be viewed as an external force acting on the macromolecule, rendering force-based experimental or computational techniques, a very general approach to the mechanics involved in protein allostery.

GENERAL SIGNIFICANCE

This unifying view might aid to decipher how complex allosteric protein machineries are regulated on the single molecular level.

Viscoelastic transition and yield strain of the folded protein

For proteins, the mechanical properties of the folded state are directly related to function, which generally entails conformational motion. Through sub-Angstrom resolution measurements of the AC mechanical susceptibility of a globular protein we describe a new fundamental materials property of the folded state. For increasing amplitude of the forcing, there is a reversible transition from elastic to viscoelastic response. At fixed frequency, the amplitude of the deformation is piecewise linear in the force, with different slopes in the elastic and viscoelastic regimes. Effectively, the protein softens beyond a yield point defined by this transition. We propose that ligand induced conformational changes generally operate in this viscoelastic regime, and that this is a universal property of the folded state.

Enzyme closure and nucleotide binding structurally lock guanylate kinase

We investigate the conformational dynamics and mechanical properties of guanylate kinase (GK) using a multiscale approach combining high-resolution atomistic molecular dynamics and low-resolution Brownian dynamics simulations. The GK enzyme is subject to large conformational changes, leading from an open to a closed form, which are further influenced by the presence of nucleotides. As suggested by recent work on simple coarse-grained models of apo-GK, we primarily focus on GK's closure mechanism with the aim to establish a detailed picture of the hierarchy and chronology of structural events essential for the enzymatic reaction. We have investigated open-versus-closed, apo-versus-holo, and substrate-versus-product-loaded forms of the GK enzyme. Bound ligands significantly modulate the mechanical and dynamical properties of GK and rigidity profiles of open and closed states hint at functionally important differences. Our data emphasizes the role of magnesium, highlights a water channel permitting active site hydration, and reveals a structural lock that stabilizes the closed form of the enzyme.

Functional modes and residue flexibility control the anisotropic response of guanylate kinase to mechanical stress

The coupling between the mechanical properties of enzymes and their biological activity is a well-established feature that has been the object of numerous experimental and theoretical works. In particular, recent experiments show that enzymatic function can be modulated anisotropically by mechanical stress. We study such phenomena using a method for investigating local flexibility on the residue scale that combines a reduced protein representation with Brownian dynamics simulations. We performed calculations on the enzyme guanylate kinase to study its mechanical response when submitted to anisotropic deformations. The resulting modifications of the protein's rigidity profile can be related to the changes in substrate binding affinity observed experimentally. Further analysis of the principal components of motion of the trajectories shows how the application of a mechanical constraint on the protein can disrupt its dynamics, thus leading to a decrease of the enzyme's catalytic rate. Eventually, a systematic probe of the protein surface led to the prediction of potential hotspots where the application of an external constraint would produce a large functional response both from the mechanical and dynamical points of view. Such enzyme-engineering approaches open the possibility to tune catalytic function by varying selected external forces.

Regulation of catch binding by allosteric transitions

An allosteric model is used to describe changes in lifetimes of biological receptor-ligand bonds subjected to an external force. Force-induced transitions between the two states of the allosteric site lead to changes in the receptor conformation. The ligand bound to the receptor fluctuates between two different potentials formed by the two receptor conformations. The effect of the force on the receptor-ligand interaction potential is described by the Bell mechanism. The probability of detecting the ligand in the bound state is found to depend on the relaxation times of both ligand and allosteric sites. An analytic expression for the bond lifetime is derived as a function of force. The formal theoretical results are used to explain the anomalous force and time dependences of the integrin-fibronectin bond lifetimes measured by atomic force microscopy (Kong, F.; et al J. Cell Biol. 2009, 185, 1275-1284). The analytic expression and model parameters describe very well all anomalous dependences identified in the experiments.

Complex molecular assemblies at hand via interactive simulations

Studying complex molecular assemblies interactively is becoming an increasingly appealing approach to molecular modeling. Here we focus on interactive molecular dynamics (IMD) as a textbook example for interactive simulation methods. Such simulations can be useful in exploring and generating hypotheses about the structural and mechanical aspects of biomolecular interactions. For the first time, we carry out low-resolution coarse-grain IMD simulations. Such simplified modeling methods currently appear to be more suitable for interactive experiments and represent a well-balanced compromise between an important gain in computational speed versus a moderate loss in modeling accuracy compared to higher resolution all-atom simulations. This is particularly useful for initial exploration and hypothesis development for rare molecular interaction events. We evaluate which applications are currently feasible using molecular assemblies from 1900 to over 300,000 particles. Three biochemical systems are discussed: the guanylate kinase (GK) enzyme, the outer membrane protease T and the soluble N-ethylmaleimide-sensitive factor attachment protein receptors complex involved in membrane fusion. We induce large conformational changes, carry out interactive docking experiments, probe lipid-protein interactions and are able to sense the mechanical properties of a molecular model. Furthermore, such interactive simulations facilitate exploration of modeling parameters for method improvement. For the purpose of these simulations, we have developed a freely available software library called MDDriver. It uses the IMD protocol from NAMD and facilitates the implementation and application of interactive simulations. With MDDriver it becomes very easy to render any particle-based molecular simulation engine interactive. Here we use its implementation in the Gromacs software as an example.

Elastic energy of protein-DNA chimeras

We present experimental measurements of the equilibrium elastic energy of protein-DNA chimeras, for two different sets of attachment points of the DNA "molecular spring" on the surface of the protein. Combining these with measurements of the enzyme's activity under stress and a mechanical model of the system, we determine how the elastic energy is partitioned between the DNA and the protein. The analysis shows that the protein is mechanically stiffer than the DNA spring.

Protein-DNA chimeras: synthesis of two-arm chimeras and non-mechanical effects of the DNA spring

DNA molecular springs have recently been used to control the activity of enzymes and ribozymes. In this approach, the mechanical stress exerted by the molecular spring alters the enzyme's conformation and thus the enzymatic activity. Here we describe a method alternative to our previous one to attach DNA molecular springs to proteins, where two separate DNA 'arms' are coupled to the protein and subsequently ligated. We report certain non-mechanical effects associated with the DNA spring observed in some chimeras with specific DNA sequences and the nucleotide binding enzyme guanylate kinase. If a ssDNA 'arm' is attached to the protein by one end only, we find that in some cases (depending on the DNA sequence and attachment point on the protein's surface) the unhybridized DNA arm inhibits the enzyme, while hybridization of the DNA arm leads to an apparent activation of the enzyme. One interpretation is that, in these cases, hybridization of the DNA arm removes it from the vicinity of the active site of the enzyme. We show how mechanical and non-mechanical effects of the DNA spring can be distinguished. This is important if one wants to use the protein-DNA chimeras to quantitatively study the response of the enzyme to mechanical perturbations.

Controlling proteins through molecular springs

We argue that the mechanical control of proteins-the notion of controlling chemical reactions and processes by mechanics-is conceptually interesting. We give a brief review of the main accomplishments so far, leading to our present approach of using DNA molecular springs to exert controlled stresses on proteins. Our focus is on the physical principles that underlie both artificial mechanochemical devices and natural mechanisms of allostery.

Elastic energy driven polymerization

We present a molecular system where polymerization is controlled externally by tuning the elastic energy of the monomers. The elastic energy, provided by a DNA molecular spring, destabilizes the monomer state through a process analogous to domain swapping. This energy can be large (of approximately 10 kT) and thus drive polymerization at relatively low monomer concentrations. The monomer-dimer equilibrium provides a measurement of the elastic energy of the monomer, which in this construction appears limited by kink formation in the DNA molecular spring, in accord with previous theoretical and experimental investigations of the elasticity of sharply bent DNA.

DNA and RNA-controlled switching of protein kinase activity

Protein switches use the binding energy gained upon recognition of ligands to modulate the conformation and binding properties of protein segments. We explored whether the programmable nucleic acid mediated recognition might be used to design or mimic constraints that limit the conformational freedom of peptide segments. The aim was to design nucleic acid-peptide conjugates in which the peptide portion of the conjugate would change the affinity for a protein target upon hybridization. This approach was used to control the affinity of a PNA-phosphopeptide conjugate for the signal transduction protein Src kinase, which binds the cognate phosphopeptides in a linear conformation. Peptide-nucleic acid arms were attached to known peptide binders. The chimeric molecules were studied in three modes: 1) as single strands, 2) constrained by intermolecular hybridization (duplex formation) and 3) constrained by intramolecular hybridization (hairpin formation). Of note, duplexes that were designed to accommodate bulged peptide structures (for example, in hairpins or bulges) had lower binding affinities than duplexes in which the peptide was allowed to adopt a more relaxed conformation. Greater than 90-fold differences in binding affinities were observed. It was, thus, feasible to make use of DNA hybridization to reversibly switch from no to almost complete inhibition of Src-SH2-peptide binding, and vice versa. A series of DNA and PNA-based hybridization experiments revealed the importance of charges and conformational effects. Nucleic acid mediated switching was extended to the use of RNA; this enabled a regulation of the enzymatic activity of the Src kinase. The proof-of-principle results demonstrate for the first time that PNA-peptide chimeras can transduce changes of the concentration of a given RNA molecule to changes of the activity of a signal transduction enzyme.

Allosteric role of the large-scale domain opening in biological catch-binding

The proposed model demonstrates the allosteric role of the two-domain region of the receptor protein in the increased lifetimes of biological receptor/ligand bonds subjected to an external force. The interaction between the domains is represented by a bounded potential, containing two minima corresponding to the attached and separated conformations of the two protein domains. The dissociative potential with a single minimum describing receptor/ligand binding fluctuates between deep and shallow states, depending on whether the domains are attached or separated. A number of valuable analytic expressions are derived and are used to interpret experimental data for two catch bonds. The P-selectin/P-selectin-glycoprotein-ligand-1 (PSGL-1) bond is controlled by the interface between the epidermal growth factor (EGF) and lectin domains of P-selectin, and the type 1 fimbrial adhesive protein (FimH)/mannose bond is governed by the interface between the lectin and pilin domains of FimH. Catch-binding occurs in these systems when the external force stretches the receptor proteins and increases the interdomain distance. The allosteric effect is supported by independent measurements, in which the domains are kept separated by attachment of another ligand. The proposed model accurately describes the experimentally observed anomalous behavior of the lifetimes of the P-selectin/PSGL-1 and FimH/mannose complexes as a function of applied force and provides valuable insights into the mechanism of catch-binding.

Efficient control of group I intron ribozyme catalysis by DNA constraints

Double-stranded DNA constraints enable efficient control of catalysis by a large multi-domain group I intron ribozyme.