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.

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.

Insights into the origin of the high energy-conversion efficiency of F1-ATPase.

Our understanding of the rotary-coupling mechanism of F1-ATPase has been greatly enhanced in the last decade by advances in X-ray crystallography, single-molecular imaging, and theoretical models. Recently, Volkán-Kacsó and Marcus [S. Volkán-Kacsó, R. A. Marcus, Proc. Natl. Acad. Sci. U.S.A. 112, 14230 (2015)] presented an insightful thermodynamic model based on the Marcus reaction theory coupled with an elastic structural deformation term to explain the observed γ-rotation angle dependence of the adenosine triphosphate (ATP)/adenosine diphosphate (ADP) exchange rates of F1-ATPase. Although the model is successful in correlating single-molecule data, it is not in agreement with the available theoretical results. We describe a revision of the model, which leads to consistency with the simulation results and other experimental data on the F1-ATPase rotor compliance. Although the free energy liberated on ATP hydrolysis by F1-ATPase is rapidly dissipated as heat and so cannot contribute directly to the rotation, we show how, nevertheless, F1-ATPase functions near the maximum possible efficiency. This surprising result is a consequence of the differential binding of ATP and its hydrolysis products ADP and Pi along a well-defined pathway.

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.

Development of Noonan syndrome by deregulation of allosteric SOS autoactivation.

Ras family proteins play an essential role in several cellular functions, including growth, differentiation, and survival. The mechanism of action of Ras mutants in Costello syndrome and cancers has been identified, but the contribution of Ras mutants to Noonan syndrome, a genetic disorder that prevents normal development in various parts of the body, is unknown. Son of Sevenless (SOS) is a Ras guanine nucleotide exchange factor. In response to Ras-activating cell signaling, SOS autoinhibition is released and is followed by accelerative allosteric feedback autoactivation. Here, using mutagenesis-based kinetic and pulldown analyses, we show that Noonan syndrome Ras mutants I24N, T50I, V152G, and D153V deregulate the autoactivation of SOS to populate their active form. This previously unknown process has been linked so far only to the development of Noonan syndrome. In contrast, other Noonan syndrome Ras mutants-V14I, T58I, and G60E-populate their active form by deregulation of the previously documented Ras GTPase activities. We propose a novel mechanism responsible for the deregulation of SOS autoactivation, where I24N, T50I, V152G, and D153V Ras mutants evade SOS autoinhibition. Consequently, they are capable of forming a complex with the SOS allosteric site, thus aberrantly promoting SOS autoactivation, resulting in the population of active Ras mutants in cells. The results of this study elucidate the molecular mechanism of the Ras mutant-mediated development of Noonan syndrome.

A simplified charge projection scheme for long-range electrostatics in ab initio QM/MM calculations.

In a previous work [Pan et al., Molecules 23, 2500 (2018)], a charge projection scheme was reported, where outer molecular mechanical (MM) charges [>10 Å from the quantum mechanical (QM) region] were projected onto the electrostatic potential (ESP) grid of the QM region to accurately and efficiently capture long-range electrostatics in ab initio QM/MM calculations. Here, a further simplification to the model is proposed, where the outer MM charges are projected onto inner MM atom positions (instead of ESP grid positions). This enables a representation of the long-range MM electrostatic potential via augmentary charges (AC) on inner MM atoms. Combined with the long-range electrostatic correction function from Cisneros et al. [J. Chem. Phys. 143, 044103 (2015)] to smoothly switch between inner and outer MM regions, this new QM/MM-AC electrostatic model yields accurate and continuous ab initio QM/MM electrostatic energies with a 10 Å cutoff between inner and outer MM regions. This model enables efficient QM/MM cluster calculations with a large number of MM atoms as well as QM/MM calculations with periodic boundary conditions.

Structural Basis for GTP versus ATP Selectivity in the NMP Kinase AK3.

ATP and GTP are exceptionally important molecules in biology with multiple, and often discrete, functions. Therefore, enzymes that bind to either of them must develop robust mechanisms to selectively utilize one or the other. Here, this specific problem is addressed by molecular studies of the human NMP kinase AK3, which uses GTP to phosphorylate AMP. AK3 plays an important role in the citric acid cycle, where it is responsible for GTP/GDP recycling. By combining a structural biology approach with functional experiments, we present a comprehensive structural and mechanistic understanding of the enzyme. We discovered that AK3 functions by recruitment of GTP to the active site, while ATP is rejected and nonproductively bound to the AMP binding site. Consequently, ATP acts as an inhibitor with respect to GTP and AMP. The overall features with specific recognition of the correct substrate and nonproductive binding by the incorrect substrate bear a strong similarity to previous findings for the ATP specific NMP kinase adenylate kinase. Taken together, we are now able to provide the fundamental principles for GTP and ATP selectivity in the large NMP kinase family. As a side-result originating from nonlinearity of chemical shifts in GTP and ATP titrations, we find that protein surfaces offer a general and weak binding affinity for both GTP and ATP. These nonspecific interactions likely act to lower the available intracellular GTP and ATP concentrations and may have driven evolution of the Michaelis constants of NMP kinases accordingly.

Rhizolutin, a Novel 7/10/6-Tricyclic Dilactone, Dissociates Misfolded Protein Aggregates and Reduces Apoptosis/Inflammation Associated with Alzheimer's Disease.

Rhizolutin (1) was discovered as a natural product of ginseng-rhizospheric Streptomyces sp. WON17. Its structure features an unprecedented 7/10/6-tricyclic dilactone carbon skeleton composed of dimethylcyclodecatriene flanked by a 7-membered and a 6-membered lactone ring based on spectroscopic analysis. During an unbiased screening of natural product libraries, this novel compound was found to dissociate amyloid-ß (Aß) plaques and tau tangles, which are key pathological hallmarks of Alzheimer's disease (AD). Rhizolutin treatment of APP/PS1 double transgenic mice with AD significantly dissociated hippocampal plaques. In vitro, rhizolutin substantially decreased Aß-induced apoptosis and inflammation in neuronal and glial cells. Our findings introduce a unique chemical entity that targets Aß and tau concurrently by mimicking misfolded protein clearance mechanisms of immunotherapy, which is prominently investigated in clinical trials.

Trapping the ATP binding state leads to a detailed understanding of the F1-ATPase mechanism.

The rotary motor enzyme FoF1-ATP synthase uses the proton-motive force across a membrane to synthesize ATP from ADP and Pi (H2PO4(-)) under cellular conditions that favor the hydrolysis reaction by a factor of 2 × 10(5). This remarkable ability to drive a reaction away from equilibrium by harnessing an external force differentiates it from an ordinary enzyme, which increases the rate of reaction without shifting the equilibrium. Hydrolysis takes place in the neighborhood of one conformation of the catalytic moiety F1-ATPase, whose structure is known from crystallography. By use of molecular dynamics simulations we trap a second structure, which is rotated by 40° from the catalytic dwell conformation and represents the state associated with ATP binding, in accord with single-molecule experiments. Using the two structures, we show why Pi is not released immediately after ATP hydrolysis, but only after a subsequent 120° rotation, in agreement with experiment. A concerted conformational change of the α3ß3 crown is shown to induce the 40° rotation of the γ-subunit only when the ßE subunit is empty, whereas with Pi bound, ßE serves as a latch to prevent the rotation of γ. The present results provide a rationalization of how F1-ATPase achieves the coupling between the small changes in the active site of ßDP and the 40° rotation of γ.

Identification and Molecular Interaction Studies of Thyroid Hormone Receptor Disruptors among Household Dust Contaminants.

Thyroid hormone disrupting chemicals (THDCs), often found abundantly in the environment, interfere with normal thyroid hormone signaling and induce physiological malfunctions, possibly by affecting thyroid hormone receptors (THRs). Indoor dust ingestion is a significant human exposure route of THDCs, raising serious concerns for human health. Here, we developed a virtual screening protocol based on an ensemble of X-ray crystallographic structures of human THRß1 and the generalized Born solvation model to identify potential THDCs targeting the human THRß1 isoform. The protocol was applied to virtually screen an in-house indoor dust contaminant inventory, yielding 31 dust contaminants as potential THRß1 binders. Five predicted binders and one negative control were tested using isothermal titration calorimetry, of which four, i.e., 2,4,5-trichlorophenoxyacetic acid (2,4,5-T), bisphenol A (3-chloro-2-hydroxypropyl) (2,3-dihydroxypropyl) ether (BADGE-HCl-H2O), 2,2',4,4'-tetrahydroxybenzophenone (BP2), and 2,4-dichlorophenoxyacetic acid (2,4-D), were identified as THRß1 binders with binding affinities ranging between 60 µM and 460 µM. Molecular dynamics (MD) simulations were employed to examine potential binding modes of these binders and provided a rationale for explaining their specific recognition by THRß1. The combination of in vitro binding affinity measurements and MD simulations allowed identification of four new potential THR-targeting THDCs that have been found in household dust. We suggest using the developed structure-based virtual screening protocol to identify and prioritize testing of potential THDCs.

The spontaneous replication error and the mismatch discrimination mechanisms of human DNA polymerase ß.

To provide molecular-level insights into the spontaneous replication error and the mismatch discrimination mechanisms of human DNA polymerase ß (polß), we report four crystal structures of polß complexed with dG•dTTP and dA•dCTP mismatches in the presence of Mg2+ or Mn2+. The Mg(2+)-bound ground-state structures show that the dA•dCTP-Mg2+ complex adopts an 'intermediate' protein conformation while the dG•dTTP-Mg2+ complex adopts an open protein conformation. The Mn(2+)-bound 'pre-chemistry-state' structures show that the dA•dCTP-Mn2+ complex is structurally very similar to the dA•dCTP-Mg2+ complex, whereas the dG•dTTP-Mn2+ complex undergoes a large-scale conformational change to adopt a Watson-Crick-like dG•dTTP base pair and a closed protein conformation. These structural differences, together with our molecular dynamics simulation studies, suggest that polß increases replication fidelity via a two-stage mismatch discrimination mechanism, where one is in the ground state and the other in the closed conformation state. In the closed conformation state, polß appears to allow only a Watson-Crick-like conformation for purine•pyrimidine base pairs, thereby discriminating the mismatched base pairs based on their ability to form the Watson-Crick-like conformation. Overall, the present studies provide new insights into the spontaneous replication error and the replication fidelity mechanisms of polß.

Role of Protein Dynamics in Allosteric Control of the Catalytic Phosphoryl Transfer of Insulin Receptor Kinase.

The catalytic and allosteric mechanisms of insulin receptor kinase (IRK) are investigated by a combination of ab initio and semiempirical quantum mechanical and molecular mechanical (QM/MM) methods and classical molecular dynamics (MD) simulations. The simulations reveal that the catalytic reaction proceeds in two steps, starting with the transfer of a proton from substrate Tyr to the catalytic Asp1132, followed by the phosphoryl transfer from ATP to substrate Tyr. The enhancement of the catalytic rate of IRK upon phosphorylations in the enzyme's activation loop is found to occur mainly via changes to the free energy landscape of the proton transfer step, favoring the proton transfer in the fully phosphorylated enzyme. In contrast, the effects of the phosphorylations on the phosphoryl transfer are smaller. Equilibrium MD simulations show that IRK phosphorylations affect the protein dynamics of the enzyme before the proton transfer to Asp1132 with only a minor effect after the proton transfer. This finding is consistent with the large change in the proton transfer free energy and the smaller change in the free energy barrier of phosphoryl transfer found by QM/MM simulations. Taken together, the present results provide details on how IRK phosphorylation exerts allosteric control of the catalytic activity via modifications of protein dynamics and free energy landscape of catalytic reaction. The results also highlight the importance of protein dynamics in connecting protein allostery and catalysis to control catalytic activity of enzymes.

Encounter and extrusion of an intrahelical lesion by a DNA repair enzyme.

How living systems detect the presence of genotoxic damage embedded in a million-fold excess of undamaged DNA is an unresolved question in biology. Here we have captured and structurally elucidated a base-excision DNA repair enzyme, MutM, at the stage of initial encounter with a damaged nucleobase, 8-oxoguanine (oxoG), nested within a DNA duplex. Three structures of intrahelical oxoG-encounter complexes are compared with sequence-matched structures containing a normal G base in place of an oxoG lesion. Although the protein-DNA interfaces in the matched complexes differ by only two atoms-those that distinguish oxoG from G-their pronounced structural differences indicate that MutM can detect a lesion in DNA even at the earliest stages of encounter. All-atom computer simulations show the pathway by which encounter of the enzyme with the lesion causes extrusion from the DNA duplex, and they elucidate the critical free energy difference between oxoG and G along the extrusion pathway.

Strandwise translocation of a DNA glycosylase on undamaged DNA.

Base excision repair of genotoxic nucleobase lesions in the genome is critically dependent upon the ability of DNA glycosylases to locate rare sites of damage embedded in a vast excess of undamaged DNA, using only thermal energy to fuel the search process. Considerable interest surrounds the question of how DNA glycosylases translocate efficiently along DNA while maintaining their vigilance for target damaged sites. Here, we report the observation of strandwise translocation of 8-oxoguanine DNA glycosylase, MutM, along undamaged DNA. In these complexes, the protein is observed to translocate by one nucleotide on one strand while remaining untranslocated on the complementary strand. We further report that alterations of single base-pairs or a single amino acid substitution (R112A) can induce strandwise translocation. Molecular dynamics simulations confirm that MutM can translocate along DNA in a strandwise fashion. These observations reveal a previously unobserved mode of movement for a DNA-binding protein along the surface of DNA.

Unsuspected pathway of the allosteric transition in hemoglobin.

Large conformational transitions play an essential role in the function of many proteins, but experiments do not provide the atomic details of the path followed in going from one end structure to the other. For the hemoglobin tetramer, the transition path between the unliganded (T) and tetraoxygenated (R) structures is not known, which limits our understanding of the cooperative mechanism in this classic allosteric system, where both tertiary and quaternary changes are involved. The conjugate peak refinement algorithm is used to compute an unbiased minimum energy path at atomic detail between the two end states. Although the results confirm some of the proposals of Perutz [Perutz MF (1970) Stereochemistry of cooperative effects in haemoglobin. Nature 228:726-734], the subunit motions do not follow the textbook description of a simple rotation of one αß-dimer relative to the other. Instead, the path consists of two sequential quaternary rotations, each involving different subdomains and axes. The quaternary transitions are preceded and followed by phases of tertiary structural changes. The results explain the recent photodissociation measurements, which suggest that the quaternary transition has a fast (2 µs) as well as a slow (20 µs) component and provide a testable model for single molecule FRET experiments.

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.

CHARMM-GUI QM/MM Interfacer for a Quantum Mechanical and Molecular Mechanical (QM/MM) Simulation Setup: 1. Semiempirical Methods.

Quantum mechanical (QM) treatments, when combined with molecular mechanical (MM) force fields, can effectively handle enzyme-catalyzed reactions without significantly increasing the computational cost. In this context, we present CHARMM-GUI QM/MM Interfacer, a web-based cyberinfrastructure designed to streamline the preparation of various QM/MM simulation inputs with ligand modification. The development of QM/MM Interfacer has been achieved through integration with existing CHARMM-GUI modules, such as PDB Reader and Manipulator, Solution Builder, and Membrane Builder. In addition, new functionalities have been developed to facilitate the one-stop preparation of QM/MM systems and enable interactive and intuitive ligand modifications and QM atom selections. QM/MM Interfacer offers support for a range of semiempirical QM methods, including AM1(+/d), PM3(+/PDDG), MNDO(+/d, +/PDDG), PM6, RM1, and SCC-DFTB, tailored for both AMBER and CHARMM. A nontrivial setup related to ligand modification, link-atom insertion, and charge distribution is automatized through intuitive user interfaces. To illustrate the robustness of QM/MM Interfacer, we conducted QM/MM simulations of three enzyme-substrate systems: dihydrofolate reductase, insulin receptor kinase, and oligosaccharyltransferase. In addition, we have created three tutorial videos about building these systems, which can be found at https://www.charmm-gui.org/demo/qmi. QM/MM Interfacer is expected to be a valuable and accessible web-based tool that simplifies and accelerates the setup process for hybrid QM/MM simulations.

Enforced presentation of an extrahelical guanine to the lesion recognition pocket of human 8-oxoguanine glycosylase, hOGG1.

A poorly understood aspect of DNA repair proteins is their ability to identify exceedingly rare sites of damage embedded in a large excess of nearly identical undamaged DNA, while catalyzing repair only at the damaged sites. Progress toward understanding this problem has been made by comparing the structures and biochemical behavior of these enzymes when they are presented with either a target lesion or a corresponding undamaged nucleobase. Trapping and analyzing such DNA-protein complexes is particularly difficult in the case of base extrusion DNA repair proteins because of the complexity of the repair reaction, which involves extrusion of the target base from DNA followed by its insertion into the active site where glycosidic bond cleavage is catalyzed. Here we report the structure of a human 8-oxoguanine (oxoG) DNA glycosylase, hOGG1, in which a normal guanine from DNA has been forcibly inserted into the enzyme active site. Although the interactions of the nucleobase with the active site are only subtly different for G versus oxoG, hOGG1 fails to catalyze excision of the normal nucleobase. This study demonstrates that even if hOGG1 mistakenly inserts a normal base into its active site, the enzyme can still reject it on the basis of catalytic incompatibility.

The water R1(ω) NMRD profiles of a hydrated protein from molecular dynamics simulation.

The hydration of a protein, Peroxiredoxin 5, is obtained from a molecular dynamics simulation and compared with the picture of hydration which is obtained by analysing the water proton R1 NMRD profiles using a generally accepted relaxation model [K. Venu, V. P. Denisov and B. Halle, J. Am. Chem. Soc., 1997, 119, 3122]. The discrepancy between the hydration pictures derived from the water R1(ω0)-NMRD profiles and MD is relevant in a discussion of the factors behind the stretched NMRD profile, the distribution of orientational order parameters and residence times of buried water used in the NMRD model.

Protein dynamics: The future is bright and complicated!

Biological life depends on motion, and this manifests itself in proteins that display motion over a formidable range of time scales spanning from femtoseconds vibrations of atoms at enzymatic transition states, all the way to slow domain motions occurring on micro to milliseconds. An outstanding challenge in contemporary biophysics and structural biology is a quantitative understanding of the linkages among protein structure, dynamics, and function. These linkages are becoming increasingly explorable due to conceptual and methodological advances. In this Perspective article, we will point toward future directions of the field of protein dynamics with an emphasis on enzymes. Research questions in the field are becoming increasingly complex such as the mechanistic understanding of high-order interaction networks in allosteric signal propagation through a protein matrix, or the connection between local and collective motions. In analogy to the solution to the "protein folding problem," we argue that the way forward to understanding these and other important questions lies in the successful integration of experiment and computation, while utilizing the present rapid expansion of sequence and structure space. Looking forward, the future is bright, and we are in a period where we are on the doorstep to, at least in part, comprehend the importance of dynamics for biological function.