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

Role of Repeated Conformational Transitions in Substrate Binding 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.

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

The human estrogen receptor α (hER α ) is involved in the regulation of growth, development, and tissue homeostasis. Agonists that bind to the receptor’s ligand-binding domain (LBD) lead to recruitment of coactivators and the enhancement of gene expression. In contrast, antagonists bind to the LBD and block the binding of coactivators thus decreasing gene expressions. In this work, we carry out simulations using the AWSEM (Associative memory, Water mediated, Structure and Energy Model)-Suite force field along with the 3SPN.2C force field for DNA to predict the structure of hER α and study its dynamics when binding to DNA and coactivators. Using simulations of antagonist-bound hER α and agonist-bound hER α by themselves and also along with bound DNA and coactivators, principal component analyses and free energy landscape analyses capture the pathway of domain–domain communication for agonist-bound hER α . This communication is mediated through the hinge domains that are ordinarily intrinsically disordered. These disordered segments manipulate the hinge domains much like the strings of a marionette as they twist in different ways when antagonists or agonists are bound to the ligand-binding domain. 

The marionette mechanism of domain–domain communication in the antagonist, agonist, and coactivator responses of the estrogen receptor

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. Mutation of multiple conserved residues leads to increased substrate inhibition, suggesting a positive selection during evolution. 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. Significance Statement How conformational dynamics affect the catalytic activity of enzymes remains a topic of active debate. We focus here on the domain closure dynamics of adenylate kinase (AK) and how they are affected by substrate inhibition. By screening an extensive mutant library, we show that this feature of the enzyme is well conserved in evolution. Importantly, domain closure is required in order to bring AK’s substrates close together for their chemical reaction; single-molecule FRET studies directly measure the populations of the open and closed states. We find that overpopulation of the closed state can be detrimental to activity. The results allow us to develop a kinetic model that properly accounts for AK kinetics by combining conformational dynamics and biochemical steps.

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

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

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

Uptake kinetics and delivery mechanisms of nanoparticles (NPs) in crop plants need to be urgently understood for the application of nanotechnology in agriculture as delivery systems for eco-friendly nanoagrochemicals. Here, we investigated the uptake kinetics, translocation pathway, and key internalization process of graphene in wheat (Triticum aestivum L.) by applying three specific hydroponic cultivation methods (submerging, hanging, and split-root). Quantification results on the uptake of carbon-14 radiolabeled graphene in each tissue indicated that graphene could enter the root of wheat and further translocate to the shoot with a low delivery rate (<2%). Transmission electron microscopy (TEM) images showed that internalized graphene was transported to adjacent cells through the plasmodesmata, clearly indicating the symplastic pathway of graphene translocation. The key site for the introduction of graphene into root cells for translocation through the symplastic pathway is evidenced to be the apex of growing root hair, where the newly constructed primary cell wall is much thinner. The confirmation of uptake kinetics and delivery mechanisms is useful for the development of nanotechnology in sustainable agriculture, especially for graphene serving as the delivery vector for pesticides, genes, and sensors.

Root Hair Apex is the Key Site for Symplastic Delivery of Graphene into Plants.

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.

https://pubs.acs.org/doi/pdf/10.1021/acs.jpcb.2c03832

Jiajun Lu

Papers

Role of Repeated Conformational Transitions in Substrate Binding of Adenylate Kinase

Root Hair Apex is the Key Site for Symplastic Delivery of Graphene into Plants.

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