Structural basis for ligand binding to an enzyme by a conformational selection pathway

Protecting Lyophilized Escherichia coli Adenylate Kinase

Drying protein-based drugs, usually via lyophilization, can facilitate storage at ambient temperature and improve accessibility but many proteins cannot withstand drying and must be formulated with protective additives called excipients. However, mechanisms of protection are poorly understood, precluding rational formulation design. To better understand dry proteins and their protection, we examine Escherichia coli adenylate kinase (AdK) lyophilized alone and with the additives trehalose, maltose, bovine serum albumin, cytosolic abundant heat soluble protein D, histidine, and arginine. We apply liquid-observed vapor exchange NMR to interrogate the residue-level structure in the presence and absence of additives. We pair these observations with differential scanning calorimetry data of lyophilized samples and AdK activity assays with and without heating. We show that the amino acids do not preserve the native structure as well as sugars or proteins and that after heating the most stable additives protect activity best.

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.

CryoEM-sampling of metastable conformations appearing in cofactor-ligand association and catalysis of glutamate dehydrogenase

Kinetic aspects of enzymatic reactions are described by equations based on the Michaelis-Menten theory for the initial stage. However, the kinetic parameters provide little information on the atomic mechanism of the reaction. In this study, we analyzed structures of glutamate dehydrogenase in the initial and steady stages of the reaction using cryoEM at near-atomic resolution. In the initial stage, four metastable conformations displayed different domain motions and cofactor/ligand association modes. The most striking finding was that the enzyme-cofactor-substrate complex, treated as a single state in the enzyme kinetic theory, comprised at least three different metastable conformations. In the steady stage, seven conformations, including derivatives from the four conformations in the initial stage, made the reaction pathway complicated. Based on the visualized conformations, we discussed stage-dependent pathways to illustrate the dynamics of the enzyme in action.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

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.

Coenzyme-binding pathway on glutamate dehydrogenase suggested from multiple-binding sites visualized by cryo-electron microscopy

The structure of hexameric glutamate dehydrogenase (GDH) in the presence of the coenzyme nicotinamide adenine dinucleotide phosphate (NADP) was visualized using cryogenic transmission electron microscopy to investigate the ligand-binding pathways to the active site of the enzyme. Each subunit of GDH comprises one hexamer-forming core domain and one nucleotide-binding domain (NAD domain), which spontaneously opens and closes the active-site cleft situated between the two domains. In the presence of NADP, the potential map of GDH hexamer, assuming D3 symmetry, was determined at a resolution of 2.4 Å, but the NAD domain was blurred due to the conformational variety. After focused classification with respect to the NAD domain, the potential maps interpreted as NADP molecules appeared at five different sites in the active-site cleft. The subunits associated with NADP molecules were close to one of the four metastable conformations in the unliganded state. Three of the five binding sites suggested a pathway of NADP molecules to approach the active-site cleft for initiating the enzymatic reaction. The other two binding modes may rarely appear in the presence of glutamate, as demonstrated by the reaction kinetics. Based on the visualized structures and the results from the enzymatic kinetics, we discussed the binding modes of NADP to GDH in the absence and presence of glutamate.

Insights into Enzymatic Catalysis from Binding and Hydrolysis of Diadenosine Tetraphosphate by E. coli Adenylate Kinase

Adenylate kinases play a crucial role in cellular energy homeostasis through the interconversion of ATP, AMP, and ADP in all living organisms. Here, we explore how adenylate kinase (AdK) from Escherichia coli interacts with diadenosine tetraphosphate (AP4A), a putative alarmone associated with transcriptional regulation, stress, and DNA damage response. From a combination of EPR and NMR spectroscopy together with X-ray crystallography, we found that AdK interacts with AP4A with two distinct modes that occur on disparate time scales. First, AdK dynamically interconverts between open and closed states with equal weights in the presence of AP4A. On a much slower time scale, AdK hydrolyses AP4A, and we suggest that the dynamically accessed substrate-bound open AdK conformation enables this hydrolytic activity. The partitioning of the enzyme into open and closed states is discussed in relation to a recently proposed linkage between active site dynamics and collective conformational dynamics.

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.

MCL-1 promiscuity and the structural resilience of its binding partners

The allosteric protein MCL-1 and its natural inhibitors, the BH3-only proteins PUMA, BIM, and NOXA regulate apoptosis by interacting promiscuously within an entangled binding network. Little is known about the transient processes and dynamic conformational fluctuations that are the basis for the formation and stability of the MCL-1/BH3-only complex. In this study, we designed photoswitchable versions of MCL-1/PUMA and MCL-1/NOXA, and investigated the protein response after an ultrafast photo-perturbation with transient infrared spectroscopy. We observed partial α-helical unfolding in all cases, albeit on strongly varying timescales (1.6 ns for PUMA, 9.7 ns for the previously studied BIM, and 85 ns for NOXA). These differences are interpreted as a BH3-only-specific "structural resilience" to defy the perturbation while remaining in MCL-1's binding pocket. Thus, the presented insights could help to better understand the differences between PUMA, BIM, and NOXA, the promiscuity of MCL-1, in general, and the role of the proteins in the apoptotic network.

Conformational exchange divergence along the evolutionary pathway of eosinophil-associated ribonucleases

The evolutionary role of conformational exchange in the emergence and preservation of function within structural homologs remains elusive. While protein engineering has revealed the importance of flexibility in function, productive modulation of atomic-scale dynamics has only been achieved on a finite number of distinct folds. Allosteric control of unique members within dynamically diverse structural families requires a better appreciation of exchange phenomena. Here, we examined the functional and structural role of conformational exchange within eosinophil-associated ribonucleases. Biological and catalytic activity of various EARs was performed in parallel to mapping their conformational behavior on multiple timescales using NMR and computational analyses. Despite functional conservation and conformational seclusion to a specific domain, we show that EARs can display similar or distinct motional profiles, implying divergence rather than conservation of flexibility. Comparing progressively more distant enzymes should unravel how this subfamily has evolved new functions and/or altered their behavior at the molecular level.

Simple mechanisms for the evolution of protein complexity

Proteins are tiny models of biological complexity: specific interactions among their many amino acids cause proteins to fold into elaborate structures, assemble with other proteins into higher-order complexes, and change their functions and structures upon binding other molecules. These complex features are classically thought to evolve via long and gradual trajectories driven by persistent natural selection. But a growing body of evidence from biochemistry, protein engineering, and molecular evolution shows that naturally occurring proteins often exist at or near the genetic edge of multimerization, allostery, and even new folds, so just one or a few mutations can trigger acquisition of these properties. These sudden transitions can occur because many of the physical properties that underlie these features are present in simpler proteins as fortuitous by-products of their architecture. Moreover, complex features of proteins can be encoded by huge arrays of sequences, so they are accessible from many different starting points via many possible paths. Because the bridges to these features are both short and numerous, random chance can join selection as a key factor in explaining the evolution of molecular complexity.

Frustration and the Kinetic Repartitioning Mechanism of Substrate Inhibition in Enzyme Catalysis

Substrate inhibition, whereby enzymatic activity decreases with excess substrate after reaching a maximum turnover rate, is among the most elusive phenomena in enzymatic catalysis. Here, based on a dynamic energy landscape model, we investigate the underlying mechanism by performing molecular simulations and frustration analysis for a model enzyme adenylate kinase (AdK), which catalyzes the phosphoryl transfer reaction ATP + AMP ⇋ ADP + ADP. Intriguingly, these reveal a kinetic repartitioning mechanism of substrate inhibition, whereby excess substrate AMP suppresses the population of an energetically frustrated, but kinetically activated, catalytic pathway going through a substrate (ATP)-product (ADP) cobound complex with steric incompatibility. Such a frustrated pathway plays a crucial role in facilitating the bottleneck product ADP release, and its suppression by excess substrate AMP leads to a slow down of product release and overall turnover. The simulation results directly demonstrate that substrate inhibition arises from the rate-limiting product-release step, instead of the steps for populating the catalytically competent complex as often suggested in previous works. Furthermore, there is a tight interplay between the enzyme conformational equilibrium and the extent of substrate inhibition. Mutations biasing to more closed conformations tend to enhance substrate inhibition. We also characterized the key features of single-molecule enzyme kinetics with substrate inhibition effect. We propose that the above molecular mechanism of substrate inhibition may be relevant to other multisubstrate enzymes in which product release is the bottleneck step.

Pseudo-Isolated α-Helix Platform for the Recognition of Deep and Narrow Targets

Although interest in stabilized α-helical peptides as next-generation therapeutics for modulating biomolecular interfaces is increasing, peptides have limited functionality and stability due to their small size. In comparison, α-helical ligands based on proteins can make steric clash with targets due to their large size. Here, we report the design of a monomeric pseudo-isolated α-helix (mPIH) system in which proteins behave as if they are peptides. The designed proteins contain α-helix ligands that do not require any covalent chemical modification, do not have frayed ends, and importantly can make sterically favorable interactions similar to isolated peptides. An optimal mPIH showed a more than 100-fold increase in target selectivity, which might be related to the advantages in conformational selection due to the absence of frayed ends. The α-helical ligand in the mPIH displayed high thermal stability well above human body temperature and showed reversible and rapid folding/unfolding transitions. Thus, mPIH can become a promising protein-based platform for developing stabilized α-helix pharmaceuticals.

Elucidating the mechanisms underlying protein conformational switching using NMR spectroscopy

How proteins switch between various ligand-free and ligand-bound structures has been a key biophysical question ever since the postulation of the Monod-Wyman-Changeux and Koshland-Nemethy-Filmer models over six decades ago. The ability of NMR spectroscopy to provide structural and kinetic information on biomolecular conformational exchange places it in a unique position as an analytical tool to interrogate the mechanisms of biological processes such as protein folding and biomolecular complex formation. In addition, recent methodological developments in the areas of saturation transfer and relaxation dispersion have expanded the scope of NMR for probing the mechanics of transitions in systems where one or more states constituting the exchange process are sparsely populated and 'invisible' in NMR spectra. In this review, we highlight some of the strategies available from NMR spectroscopy for examining the nature of multi-site conformational exchange, using five case studies that have employed NMR, either in isolation, or in conjunction with other biophysical tools.

Time Evolution of the Millisecond Allosteric Activation of Imidazole Glycerol Phosphate Synthase

Deciphering the molecular mechanisms of enzymatic allosteric regulation requires the structural characterization of functional states and also their time evolution toward the formation of the allosterically activated ternary complex. The transient nature and usually slow millisecond time scale interconversion between these functional states hamper their experimental and computational characterization. Here, we combine extensive molecular dynamics simulations, enhanced sampling techniques, and dynamical networks to describe the allosteric activation of imidazole glycerol phosphate synthase (IGPS) from the substrate-free form to the active ternary complex. IGPS is a heterodimeric bienzyme complex whose HisH subunit is responsible for hydrolyzing glutamine and delivering ammonia for the cyclase activity in HisF. Despite significant advances in understanding the underlying allosteric mechanism, essential molecular details of the long-range millisecond allosteric activation of IGPS remain hidden. Without using a priori information of the active state, our simulations uncover how IGPS, with the allosteric effector bound in HisF, spontaneously captures glutamine in a catalytically inactive HisH conformation, subsequently attains a closed HisF:HisH interface, and finally forms the oxyanion hole in HisH for efficient glutamine hydrolysis. We show that the combined effector and substrate binding dramatically decreases the conformational barrier associated with oxyanion hole formation, in line with the experimentally observed 4500-fold activity increase in glutamine hydrolysis. The allosteric activation is controlled by correlated time-evolving dynamic networks connecting the effector and substrate binding sites. This computational strategy tailored to describe millisecond events can be used to rationalize the effect of mutations on the allosteric regulation and guide IGPS engineering efforts.

Structure determination of high-energy states in a dynamic protein ensemble

Macromolecular function frequently requires that proteins change conformation into high-energy states^1-4. However, methods for solving the structures of these functionally essential, lowly populated states are lacking. Here we develop a method for high-resolution structure determination of minorly populated states by coupling NMR spectroscopy-derived pseudocontact shifts⁵ (PCSs) with Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion⁶ (PCS-CPMG). Our approach additionally defines the corresponding kinetics and thermodynamics of high-energy excursions, thereby characterizing the entire free-energy landscape. Using a large set of simulated data for adenylate kinase (Adk), calmodulin and Src kinase, we find that high-energy PCSs accurately determine high-energy structures (with a root mean squared deviation of less than 3.5 angström). Applying our methodology to Adk during catalysis, we find that the high-energy excursion involves surprisingly small openings of the AMP and ATP lids. This previously unresolved high-energy structure solves a longstanding controversy about conformational interconversions that are rate-limiting for catalysis. Primed for either substrate binding or product release, the high-energy structure of Adk suggests a two-step mechanism combining conformational selection to this state, followed by an induced-fit step into a fully closed state for catalysis of the phosphoryl-transfer reaction. Unlike other methods for resolving high-energy states, such as cryo-electron microscopy and X-ray crystallography, our solution PCS-CPMG approach excels in cases involving domain rearrangements of smaller systems (less than 60 kDa) and populations as low as 0.5%, and enables the simultaneous determination of protein structure, kinetics and thermodynamics while proteins perform their function.

A thermodynamic atlas of proteomes reveals energetic innovation across the tree of life

Protein stability is a fundamental molecular property enabling organisms to adapt to their biological niches. How this is facilitated and whether there are kingdom specific or more general universal strategies is not known. A principal obstacle to addressing this issue is that the vast majority of proteins lack annotation, specifically thermodynamic annotation, beyond the amino acid and chromosome information derived from genome sequencing. To address this gap and facilitate future investigation into large-scale patterns of protein stability and dynamics within and between organisms, we applied a unique ensemble-based thermodynamic characterization of protein folds to a substantial portion of extant sequenced genomes. Using this approach, we compiled a database resource focused on the position-specific variation in protein stability. Interrogation of the database reveals; 1) domains of life exhibit distinguishing thermodynamic features, with eukaryotes particularly different from both archaea and bacteria, 2) the optimal growth temperature of an organism is proportional to the average apolar enthalpy of its proteome, 3) intrinsic disorder content is also proportional to the apolar enthalpy (but unexpectedly not the predicted stability at 25 °C), and 4) secondary structure and global stability information of individual proteins is extractable. We hypothesize that wider access to residue-specific thermodynamic information of proteomes will result in deeper understanding of mechanisms driving functional adaptation and protein evolution. Our database is free for download at https://afc-science.github.io/thermo-env-atlas/.

Antibodies targeting enzyme inhibition as potential tools for research and drug development

Antibodies have transformed biomedical research and are now being used for different experimental applications. Generally, the interaction of enzymes with their specific antibodies can lead to a reduction in their enzymatic activity. The effect of the antibody is dependent on its narrow i.e. the regions of the enzyme to which it is directed. The mechanism of this inhibition is rarely a direct combination of the antibodies with the catalytic site, but is rather due to steric hindrance, barring the substrate access to the active site. In several systems, however, the interaction with the antibody induces conformational changes on the enzyme that can either inhibit or enhance its catalytic activity. The extent of enzyme inhibition or enhancement is, therefore, a reflection of the nature and distribution of the various antigenic determinants on the enzyme molecule. Currently, the mode of action of many enzymes has been elucidated at the molecular level. We here review the molecular mechanisms and recent trends by which antibodies inhibit the catalytic activity of enzymes and provide examples of how specific antibodies can be useful for the neutralization of biologically active molecules.

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.

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.

Understanding the Binding Transition State After the Conformational Selection Step: The Second Half of the Molecular Recognition Process Between NS1 of the 1918 Influenza Virus and Host p85β

Biomolecular recognition often involves conformational changes as a prerequisite for binding (i.e., conformational selection) or concurrently with binding (i.e., induced-fit). Recent advances in structural and kinetic approaches have enabled the detailed characterization of protein motions at atomic resolution. However, to fully understand the role of the conformational dynamics in molecular recognition, studies on the binding transition state are needed. Here, we investigate the binding transition state between nonstructural protein 1 (NS1) of the pandemic 1918 influenza A virus and the human p85β subunit of PI3K. 1918 NS1 binds to p85β via conformational selection. We present the free-energy mapping of the transition and bound states of the 1918 NS1:p85β interaction using linear free energy relationship and ϕ-value analyses. We find that the binding transition state of 1918 NS1 and p85β is structurally similar to the bound state with well-defined binding orientation and hydrophobic interactions. Our finding provides a detailed view of how protein motion contributes to the development of intermolecular interactions along the binding reaction coordinate.

Making the invisible visible: Toward structural characterization of allosteric states, interaction networks, and allosteric regulatory mechanisms in protein kinases

Despite the established view of protein kinases as dynamic and versatile allosteric regulatory machines, our knowledge of allosteric functional states, allosteric interaction networks, and the intrinsic folding energy landscapes is surprisingly limited. We discuss the latest developments in structural characterization of allosteric molecular events underlying protein kinase dynamics and functions using structural, biophysical, and computational biology approaches. The recent studies highlighted progress in making the invisible aspects of protein kinase 'life' visible, including the determination of hidden allosteric states and mapping of allosteric energy landscapes, discovery of new mechanisms underlying ligand-induced modulation of allosteric activity, evolutionary adaptation of kinase allostery, and characterization of allosteric interaction networks as the intrinsic driver of kinase adaptability and signal transmission in the regulatory assemblies.

Sequence of Events during Peptide Unbinding from RNase S: A Complete Experimental Description

The phototriggered unbinding of the intrinsically disordered S-peptide from the RNase S complex is studied with the help of transient IR spectroscopy, covering a wide range of time scales from 100 ps to 10 ms. To that end, an azobenzene moiety has been linked to the S-peptide in a way that its helicity is disrupted by light, thereby initiating its complete unbinding. The full sequence of events is observed, starting from unfolding of the helical structure of the S-peptide on a 20 ns time scale while still being in the binding pocket of the S-protein, S-peptide unbinding after 300 μs, and the structural response of the S-protein after 3 ms. With regard to the S-peptide dynamics, the binding mechanism can be classified as an induced fit, while the structural response of the S-protein is better described as conformational selection.

Drug discovery for enzymes

Enzymes are essential, physiological catalysts involved in all processes of life, including metabolism, cellular signaling and motility, as well as cell growth and division. They are attractive drug targets because of the presence of defined substrate-binding pockets, which can be exploited as binding sites for pharmaceutical enzyme inhibitors. Understanding the reaction mechanisms of enzymes and the molecular mode of action of enzyme inhibitors is indispensable for the discovery and development of potent, efficacious, and safe novel drugs. The combination of classical concepts of enzymology with new experimental and data analysis methods opens new routes for drug discovery.

In vitro effects of synthetic muscimol and an extract from Amanita muscaria on human recombinant MAOB enzyme

The effects of synthetic muscimol and an extract from Amanita muscaria, containing this compound on the activity of human recombinant MAOB enzyme (hMAOB) were studied. Muscimol had statistically significant inducing effect on hMAOB at concentrations 0.25–5 μM, while A. muscaria extract did not influence the enzyme activity at all.

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.

Mechanisms of ligand binding

Many processes in chemistry and biology involve interactions of a ligand with its molecular target. Interest in the mechanism governing such interactions has dominated theoretical and experimental analysis for over a century. The interpretation of molecular recognition has evolved from a simple rigid body association of the ligand with its target to appreciation of the key role played by conformational transitions. Two conceptually distinct descriptions have had a profound impact on our understanding of mechanisms of ligand binding. The first description, referred to as induced fit, assumes that conformational changes follow the initial binding step to optimize the complex between the ligand and its target. The second description, referred to as conformational selection, assumes that the free target exists in multiple conformations in equilibrium and that the ligand selects the optimal one for binding. Both descriptions can be merged into more complex reaction schemes that better describe the functional repertoire of macromolecular systems. This review deals with basic mechanisms of ligand binding, with special emphasis on induced fit, conformational selection, and their mathematical foundations to provide rigorous context for the analysis and interpretation of experimental data. We show that conformational selection is a surprisingly versatile mechanism that includes induced fit as a mathematical special case and even captures kinetic properties of more complex reaction schemes. These features make conformational selection a dominant mechanism of molecular recognition in biology, consistent with the rich conformational landscape accessible to biological macromolecules being unraveled by structural biology.

Protein Conformational Dynamics upon Association with the Surfaces of Lipid Membranes and Engineered Nanoparticles: Insights from Electron Paramagnetic Resonance Spectroscopy

Detailed study of conformational rearrangements and dynamics of proteins is central to our understanding of their physiological functions and the loss of function. This review outlines the applications of the electron paramagnetic resonance (EPR) technique to study the structural aspects of proteins transitioning from a solution environment to the states in which they are associated with the surfaces of biological membranes or engineered nanoobjects. In the former case these structural transitions generally underlie functional protein states. The latter case is mostly relevant to the application of protein immobilization in biotechnological industries, developing methods for protein purification, etc. Therefore, evaluating the stability of the protein functional state is particularly important. EPR spectroscopy in the form of continuous-wave EPR or pulse EPR distance measurements in conjunction with protein spin labeling provides highly versatile and sensitive tools to characterize the changes in protein local dynamics as well as large conformational rearrangements. The technique can be widely utilized in studies of both protein-membrane and engineered nanoobject-protein complexes.

The crystal structures of Thermus thermophilus CMP kinase complexed with a phosphoryl group acceptor and donor

Nucleoside monophosphate kinases play crucial roles in biosynthesis and regeneration of nucleotides. These are bi-substrate enzymes that catalyze reversible transfers of a phosphoryl group between ATP and nucleoside monophosphate. These enzymes are comprised of the CORE domain, the NMP-binding domain, and the LID domain. Large conformational rearrangement of the three domains occurs during the catalytic cycle. Although many structures of CMP kinase have been determined, only limited structural information has been available on the conformational changes along the reaction pathway. We determined five crystal structures of CMP kinase of Thermus thermophilus HB8 in ligand-free form and the CMP "open", CMP "closed", ADP-CDP-Gd3+-, and CDP-bound forms at resolutions of 1.7, 2.2, 1.5, 1.6, and 1.7 Å, respectively. The ligand-free form was in an open conformation, whereas the structures of the CMP "closed", ADP-CDP-Gd3+-, and CDP-bound forms were in a closed conformation, in which the shift of the NMP-binding domain and LID domain caused closure of the substrate-binding cleft. Interestingly, the CMP "open" form was in an open conformation even with CMP bound, implying intrinsic conformational fluctuation. The structure of the ADP-CDP complex is the first structure of CMP kinase with a phosphoryl group donor and an acceptor. Upon simultaneous binding of ADP and CDP, the side chains of several residues in the LID domain moved toward the nucleotides without global open-closed conformational changes compared to those in the CMP "closed" and CDP complexes. These global and local conformational changes may be crucial for the substrate recognition and catalysis. The terminal phosphate groups of ADP and CDP had similar geometry to those of two ADP in AMP kinase, suggesting common catalytic mechanisms to other nucleoside monophosphate kinases. Our findings are expected to contribute to detailed understanding of the reaction mechanism of CMP kinase.

Role of substrate-product frustration on enzyme functional dynamics

Natural enzymes often have enormous catalytic power developed by evolution. Revealing the underlying physical strategy used by enzymes to achieve high catalysis efficiency is one of the central focuses in the field of biological physics. Our recent work demonstrated that multisubstrate enzymes can utilize steric frustration encountered in the substrate-product cobound complex to overcome the bottleneck of the enzymatic cycle [W. Li et al., Phys. Rev. Lett. 122, 238102 (2019)10.1103/PhysRevLett.122.238102]. However, the key atomic-level interactions by which the steric frustration contributes to the enzymatic cycle remain elusive. In this work we study the microscopic mechanism for the role of the substrate-product frustration on the key physical steps in the enzymatic cycle of adenylate kinase (AdK), a multisubstrate enzyme catalyzing the reversible phosphoryl transfer reaction ATP+AMP⇋ADP+ADP. By using atomistic molecular dynamics simulations with enhanced sampling, we showed that the competitive interactions from the phosphate groups of the substrate ATP and product ADP in the ATP-ADP cobound complex of the AdK lead to local frustration in the binding pockets. Such local frustration disrupts the hydrogen bond network around the binding pockets, which causes lowered barrier height for the opening of the enzyme conformations and expedited release of the bottleneck product ADP. Our results directly demonstrated from the atomistic level that the local frustration in the active sites of the enzyme can be utilized to facilitate the key physical steps of the enzymatic cycle, providing numerical evidence to the predictions of the previous theoretical work.

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.

Halogen Bonding Interactions of Polychlorinated Biphenyls and the Potential for Thyroid Disruption

Polychlorinated biphenyl (PCB) flame retardants are persistent pollutants and inhibit neurodevelopment, particularly in the early stages of life. Halogen bonding (XB) to the iodothyronine deiodinases (Dio) that modulate thyroid hormones (THs) is a potential mechanism for endocrine disruption. Cl⋅⋅⋅Se XB interactions of PCBs with SeMe- , a small model of the Dio active site selenocysteine, are compared with previous results on polybrominated diphenylethers (PBDEs) and THs using density functional theory. PCBs generally display weaker XB interactions compared to PBDEs and THs, consistent with the dependence of XB strength on the size of the halogen (I>Br>Cl). PCBs also do not meet a proposed energy threshold for substrates to undergo dehalogenation, suggesting they may behave as competitive inhibitors of Dio in addition to other mechanisms of endocrine disruption. XB interactions in PCBs are position-dependent, with ortho interactions slightly more favorable than meta and para interactions, suggesting that PCBs may have a greater effect on certain classes of Dio. Flexibility of PCBs around the biphenyl C-C bond is limited by ortho substitutions relative to the biphenyl linkage, which may contribute to the ability to inhibit Dio and other TH-related proteins.

Large-Scale Biomolecular Conformational Transitions Explored by a Combined Elastic Network Model and Enhanced Sampling Molecular Dynamics

Biomolecules often undergo large-scale conformational transitions when carrying out their functions. However, it is still challenging for conventional molecular dynamics simulations to provide adequate structural dynamics information to interpret associated mechanisms. Here, we present a combined elastic network model and enhanced sampling-based strategy (iterANM-IaMD) by adopting iterANM to construct initial conformation space and enhanced sampling IaMD to explore the free energy landscape along specific large-scale conformational transitions. We applied this strategy to three functionally and structurally distinct proteins (adenylate kinase, calmodulin, and p38α kinase), which undergo striking conformational change upon ligand binding. The simulation results for both free and ligand-bound proteins show qualitative and quantitative agreement with existing studies, suggesting iterANM-IaMD as an accurate and efficient tool to investigate structural dynamics involved in complicated biological processes. Our work also provides insights into the relationship between the dynamics and functionality of biomolecules.

Substrate-enzyme affinity-based surface modification strategy for endothelial cell-specific binding under shear stress

Establishing an endothelial cell (EC) monolayer on top of the blood contacting surface of grafts is considered to be a promising approach for creating a hemocompatible surface. Here we utilized the high affinity interactions between the EC plasma membrane expressed enzyme called endothelin converting enzyme-1 (ECE-1) and its corresponding substrate big Endothelin-1 (bigET-1) to engineer an EC-specific binding surface. Since enzymatic cleavage of substrates require physical interaction between the enzyme and its corresponding substrate, it was hypothesized that a surface with chemically immobilized synthetic bigET-1 will preferentially attract ECs over other types of cells found in vascular system such as vascular smooth muscle cells (VSMCs). First, the expression of ECE-1 was significantly higher in ECs, and ECs processed synthetic bigET-1 to produce ET-1 in a cell number-dependent manner. Such interaction between ECs and synthetic bigET-1 was also detectible in blood. Next, vinyl-terminated self-assembled monolayers (SAMs) were established, oxidized and activated on a glass substrate as a model to immobilize synthetic bigET-1 via amide bonds. The ECs cultured on the synthetic bigET-1-immobilized surface processed larger amount of synthetic bigET-1 to produce ET-1 compared to VSMCs (102.9±5.13 vs. 9.75±0.74 pg/ml). The number of ECs bound to the synthetic bigET-1-immobilized surface during 1 h of shearing (5dyne/cm2) was approximately 3-fold higher than that of VSMCs (46.25±12.61 vs. 15.25±3.69 cells/100×HPF). EC-specific binding of synthetic bigET-1-immobilized surface over a surface modified with collagen, a common substance for cell adhesion, was also observed. The present study demonstrated that using the substrate-enzyme affinity (SEA) of cell type-specific enzyme and its corresponding substrate can be an effective method to engineer a surface preferentially binds specific type of cells. This novel strategy might open a new route toward rapid endothelialization under dynamic conditions supporting the long-term patency of cardiovascular implants.

Exploring Reactive Conformations of Coenzyme A during Binding and Unbinding to Pyruvate Formate-Lyase

Pyruvate formate-lyase (PFL) is a glycyl radical enzyme that converts pyruvate and coenzyme A (CoA) into formate and acetyl-CoA in two half-reactions. Recently, we showed that the acetylation of the PFL active site in the first half-reaction induces subtle conformational changes, leading to the opening of a potential channel for CoA entry. Entry of CoA into the active site is crucial for the second half-reaction, involving the acetyl transfer to CoA, and the completion of the catalytic cycle. Using steered molecular dynamics (SMD) simulations, performed on acetylated and nonacetylated monomeric PFL model systems, we first of all investigate the possible entry/exit pathways of CoA with respect to the active site through the previously identified channel. We then perform umbrella sampling simulations on multiple snapshots from SMD trajectories as well as unrestrained molecular dynamics simulations starting from the final structures obtained from entry SMD, with a view to identifying possible bound states of CoA in the near vicinity of the active site. Detailed study of the unrestrained dissociation processes reveals the presence of stable and reactive bound states of CoA close to the active site, one of which is in an ideal position for triggering the second half-reaction. Examination of the spatial distributions associated with the reactive bound states allows us to discuss the free energy barriers. Umbrella sampling, performed on snapshots from unrestrained dynamics confirms the above findings. The significance of the results for the catalysis are discussed for both acetylated and nonacetylated systems.

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.

High-speed AFM reveals accelerated binding of agitoxin-2 to a K+ channel by induced fit

Agitoxin-2 (AgTx2) from scorpion venom is a potent blocker of K⁺ channels. The docking model has been elucidated, but it remains unclear whether binding dynamics are described by a two-state model (AgTx2-bound and AgTx2-unbound) or a more complicated mechanism, such as induced fit or conformational selection. Here, we observed the binding dynamics of AgTx2 to the KcsA channel using high-speed atomic force microscopy. From images of repeated binding and dissociation of AgTx2 to the channel, single-molecule kinetic analyses revealed that the affinity of the channel for AgTx2 increased during persistent binding and decreased during persistent dissociation. We propose a four-state model, including high- and low-affinity states of the channel, with relevant rate constants. An induced-fit pathway was dominant and accelerated binding by 400 times. This is the first analytical imaging of scorpion toxin binding in real time, which is applicable to various biological dynamics including channel ligands, DNA-modifier proteins, and antigen-antibody complexes.

Overcoming the Bottleneck of the Enzymatic Cycle by Steric Frustration

The enormous catalytic power of natural enzymes relies on the ability to overcome the bottleneck event in the enzymatic cycle, yet the underlying physical mechanisms are not fully understood. Here, by performing molecular simulations of the whole enzymatic cycle for a model multisubstrate enzyme with a dynamic energy landscape model, we show that multisubstrate enzymes can utilize steric frustration to facilitate the rate-limiting product-release step. During the enzymatic cycles, the bottleneck product is actively squeezed out by the binding of a new substrate at the neighboring site through the population of a substrate-product cobound complex, in which the binding pockets are frustrated due to steric incompatibility. Such steric frustration thereby enables an active mechanism of product release driven by substrate-binding energy, facilitating the enzymatic cycle.

Uncorrelated Effect of Interdomain Contact on Pin1 Isomerase Activity Reveals Positive Catalytic Cooperativity

Pin1 is a two-domain peptidyl-prolyl isomerase (PPIase) associated with neurodegeneration and tumorigenesis. The two domains, a WW and a PPIase domain, are connected by a flexible linker, making Pin1 adopt various conformations ranging from compact to extended, wherein Pin1 exhibits different extents of interdomain contact. Previous studies have shown that weakening interdomain contact increases the isomerase activity of Pin1. Here, we propose an NMR chemical shift correlation-analysis-based method that will be general for two-domain proteins to gauge two-state populations of Pin1, and we report a linker-modified mutant of Pin1 with enhanced interdomain contact and increased isomerase activity, with the latter suggesting an uncorrelated effect of interdomain contact on isomerase activity. Thus, although bindings of different substrates in the WW domain impose opposite effects on interdomain contact, in both cases, it may promote isomerization, implying cooperativity between substrate binding in the WW domain and isomerization in the PPIase domain.

The Inescapable Effects of Ribosomes on In-Cell NMR Spectroscopy and the Implications for Regulation of Biological Activity

The effects of RNA on in-cell NMR spectroscopy and ribosomes on the kinetic activity of several metabolic enzymes are reviewed. Quinary interactions between labelled target proteins and RNA broaden in-cell NMR spectra yielding apparent megadalton molecular weights in-cell. The in-cell spectra can be resolved by using cross relaxation-induced polarization transfer (CRINEPT), heteronuclear multiple quantum coherence (HMQC), transverse relaxation-optimized, NMR spectroscopy (TROSY). The effect is reproduced in vitro by using reconstituted total cellular RNA and purified ribosome preparations. Furthermore, ribosomal binding antibiotics alter protein quinary structure through protein-ribosome and protein-mRNA-ribosome interactions. The quinary interactions of Adenylate kinase, Thymidylate synthase and Dihydrofolate reductase alter kinetic properties of the enzymes. The results demonstrate that ribosomes may specifically contribute to the regulation of biological activity.

Development of Software for the In-Depth Analysis of Protein Dynamics as Determined by MALDI Mass Spectrometry-Based Hydrogen/Deuterium Exchange

Hydrogen/deuterium exchange (HDX) coupled with pepsin digestion is useful for rapidly analyzing the kinetic properties of small amounts of protein. However, the analysis of HDX by matrix-assisted laser desorption/ionization (MALDI) mass spectrometry (MS) is time-consuming due to a lack of dedicated software. Currently available software programs mainly calculate average mass shifts, even though the isotopic distribution width contains information regarding multiple protein conformations. Moreover, HDX reaction samples are typically composed of peptides that contain various numbers of deuterium atoms, which also hinders the rapid and comprehensive analysis of protein dynamics. We report here on the development of a software program "Scipas DX" that can be used to automatically analyze the hydrogen-deuterium isotopic distribution in peaks in HDX spectra and calculate the average number of atoms exchanged, the average deuteration ratio, the abundance ratio for exchanged atoms, and their fitted spectra with a high degree of accuracy within a few minutes. Analysis of the abundance ratio for exchanged atoms of a model protein, adenylate kinase 1, using Scipas DX indicate that the local structure at residues 83-106 and 107-117 are in a slow equilibrium, suggesting that these regions adopt multiple conformations that are involved in the stability and in switching between the active and inactive forms. Furthermore, precise HDX kinetics of the average deuteration ratio both confirmed the known induced conformations of two regions (residues 46-75 and 131-165) that are responsible for ligand binding and verified the novel structural dynamics of residues 107-117 and 166-196 following ligand binding to ligand-binding pockets 1 and 2, respectively. Collectively, these results highlight the usefulness and versatility of Scipas DX in MALDI-MS HDX-based analyses of protein dynamics.

Role of conformational dynamics in the evolution of novel enzyme function

The free energy landscape concept that describes enzymes as an ensemble of differently populated conformational sub-states in dynamic equilibrium is key for evaluating enzyme activity, enantioselectivity, and specificity. Mutations introduced in the enzyme sequence can alter the populations of the pre-existing conformational states, thus strongly modifying the enzyme ability to accommodate alternative substrates, revert its enantiopreferences, and even increase the activity for some residual promiscuous reactions. In this feature article, we present an overview of the current experimental and computational strategies to explore the conformational free energy landscape of enzymes. We provide a series of recent publications that highlight the key role of conformational dynamics for the enzyme evolution towards new functions and substrates, and provide some perspectives on how conformational dynamism should be considered in future computational enzyme design protocols.

Ligand-Induced Variations in Structural and Dynamical Properties Within an Enzyme Superfamily

Enzyme catalysis is a complex process involving several steps along the reaction coordinates, including substrate recognition and binding, chemical transformation, and product release. Evidence continues to emerge linking the functional and evolutionary role of conformational exchange processes in optimal catalytic activity. Ligand binding changes the conformational landscape of enzymes, inducing long-range conformational rearrangements. Using functionally distinct members of the pancreatic ribonuclease superfamily as a model system, we characterized the structural and conformational changes associated with the binding of two mononucleotide ligands. By combining NMR chemical shift titration experiments with the chemical shift projection analysis (CHESPA) and relaxation dispersion experiments, we show that biologically distinct members of the RNase superfamily display discrete chemical shift perturbations upon ligand binding that are not conserved even in structurally related members. Amino acid networks exhibiting coordinated chemical shift displacements upon binding of the two ligands are unique to each of the RNases analyzed. Our results reveal the contribution of conformational rearrangements to the observed chemical shift perturbations. These observations provide important insights into the contribution of the different ligand binding specificities and effects of conformational exchange on the observed perturbations associated with ligand binding for functionally diverse members of the pancreatic RNase superfamily.

A gratuitous β-Lactamase inducer uncovers hidden active site dynamics of the Staphylococcus aureus BlaR1 sensor domain

Increasing evidence shows that active sites of proteins have non-trivial conformational dynamics. These dynamics include active site residues sampling different local conformations that allow for multiple, and possibly novel, inhibitor binding poses. Yet, active site dynamics garner only marginal attention in most inhibitor design efforts and exert little influence on synthesis strategies. This is partly because synthesis requires a level of atomic structural detail that is frequently missing in current characterizations of conformational dynamics. In particular, while the identity of the mobile protein residues may be clear, the specific conformations they sample remain obscure. Here, we show how an appropriate choice of ligand can significantly sharpen our abilities to describe the interconverting binding poses (conformations) of protein active sites. Specifically, we show how 2-(2'-carboxyphenyl)-benzoyl-6-aminopenicillanic acid (CBAP) exposes otherwise hidden dynamics of a protein active site that binds β-lactam antibiotics. When CBAP acylates (binds) the active site serine of the β-lactam sensor domain of BlaR1 (BlaRS), it shifts the time scale of the active site dynamics to the slow exchange regime. Slow exchange enables direct characterization of inter-converting protein and bound ligand conformations using NMR methods. These methods include chemical shift analysis, 2-d exchange spectroscopy, off-resonance ROESY of the bound ligand, and reduced spectral density mapping. The active site architecture of BlaRS is shared by many β-lactamases of therapeutic interest, suggesting CBAP could expose functional motions in other β-lactam binding proteins. More broadly, CBAP highlights the utility of identifying chemical probes common to structurally homologous proteins to better expose functional motions of active sites.

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.

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.