Fine structure of conformational ensembles in adenylate kinase

Allostery can convert binding free energies into concerted domain motions in enzymes

Enzymatic mechanisms are typically inferred from structural data. However, understanding enzymes require unravelling the intricate dynamic interplay between dynamics, conformational substates, and multiple protein structures. Here, we use single-molecule nanopore analysis to investigate the catalytic conformational changes of adenylate kinase (AK), an enzyme that catalyzes the interconversion of various adenosine phosphates (ATP, ADP, and AMP). Kinetic analysis validated by hidden Markov models unravels the details of domain motions during catalysis. Our findings reveal that allosteric interactions between ligands and cofactor enable converting binding energies into directional conformational changes of the two catalytic domains of AK. These coordinated motions emerged to control the exact sequence of ligand binding and the affinity for the three different substrates, thereby guiding the reactants along the reaction coordinates. Interestingly, we find that about 10% of enzymes show altered allosteric regulation and ligand affinities, indicating that a subset of enzymes folds in alternative catalytically active forms. Since molecules or proteins might be able to selectively stabilize one of the folds, this observation suggests an evolutionary path for allostery in enzymes. In AK, this complex catalytic framework has likely emerged to prevent futile ATP/ADP hydrolysis and to regulate the enzyme for different energy needs of the cell.

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

Knowledge of the structures formed by proteins and small molecules is key to understand the molecular principles of chemotherapy and for designing new and more effective drugs. During the early stage of a drug discovery program, it is customary to predict ligand-protein complexes in silico , particularly when screening large compound databases. 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 to generate bound-like conformations of very flexible and allosteric proteins bearing multiple binding sites by exploiting only information on the unbound structure and the putative binding sites. The protocol is validated on the paradigm enzyme adenylate kinase, for which we generated a significant fraction of bound-like structures. A fraction of these conformations, employed in ensemble-docking calculations, allowed to find native-like poses of substrates and inhibitors (binding to the active form of the enzyme), as well as catalytically incompetent analogs (binding the inactive form). Our protocol provides a general framework for the generation of bound-like conformations of challenging drug targets 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, in the prediction of the impact of amino acid mutations on structure and dynamics, and in protein engineering.

Role of Repeated Conformational Transitions in Substrate Binding of Adenylate Kinase

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

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

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

Converting the E. coli Isochorismatase Nicotinamidase into γ-Lactamase

Nicotinamidase (Nic) (E.C.3.5.1.19) is a representative protein of the isochorismatase superfamily from Escherichia coli. Despite showing no (+) γ-lactamase activity, its active site constellations (ASCs) are very similar to those of two other known (+) γ-lactamases (Mhpg and RutB), indicating that it could be a latent (+) γ-lactamase. In this study, the primary sequences of the five representative proteins of the isochorismatase superfamily from E. coli were aligned, and a "lid"-like unit of a six-residue loop (112GENPLV117) was established. The Nic protein was converted to a (+) γ-lactamase by eliminating the loop. A conversion mechanism was proposed in which a more compact binding pocket is formed after lid deletion. In addition, the "shrunk" binding pocket stabilized the small substrate and the catalysis intermediate, which triggered catalysis. Moreover, we identified another latent (+) γ-lactamase in the E. coli isochorismatase superfamily and successfully converted it into an active (+) γ-lactamase. In summary, the isochorismatase superfamily is potentially a good candidate for obtaining novel (+) γ-lactamases. IMPORTANCE γ-Lactamases are important enzymatic catalysts in preparing optically pure γ-lactam enantiomers, which are high-value chiral intermediates. Different studies have presumed that the isochorismatase superfamily is a candidate to obtain novel (+) γ-lactamases. By engineering its substrate entrance tunnel, Nic, a representative protein of the isochorismatase superfamily, is converted to a (+) γ-lactamase. Tunnel engineering has proven effective in enhancing enzyme promiscuity. Therefore, the latent or active γ-lactamase activities of the isochorismatase superfamily members indicate their evolutionary path positions.

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.

Moving beyond static snapshots: Protein dynamics and the Protein Data Bank

Proteins are the molecular machines of living systems. Their dynamics are an intrinsic part of their evolutionary selection in carrying out their biological functions. Although the dynamics are more difficult to observe than a static, average structure, we are beginning to observe these dynamics and form sound mechanistic connections between structure, dynamics, and function. This progress is highlighted in case studies from myoglobin and adenylate kinase to the ribosome and molecular motors where these molecules are being probed with a multitude of techniques across many timescales. New approaches to time-resolved crystallography are allowing simple “movies” to be taken of proteins in action, and new methods of mapping the variations in cryo-electron microscopy are emerging to reveal a more complete description of life’s machines. The results of these new methods are aided in their dissemination by continual improvements in curation and distribution by the Protein Data Bank and their partners around the world.

Identifying Conformational Isomers of Organic Molecules in Solution via Unsupervised Clustering

We present a systematic approach for the identification of statistically relevant conformational macrostates of organic molecules from molecular dynamics trajectories. The approach applies to molecules characterized by an arbitrary number of torsional degrees of freedom and enables the transferability of the macrostates definition across different environments. We formulate a dissimilarity measure between molecular configurations that incorporates information on the characteristic energetic cost associated with transitions along all relevant torsional degrees of freedom. Such metric is employed to perform unsupervised clustering of molecular configurations based on the Fast Search and Find of Density Peaks algorithm. We apply this method to investigate the equilibrium conformational ensemble of Sildenafil, a conformationally complex pharmaceutical compound, in different environments including the crystal bulk, the gas phase, and three different solvents (acetonitrile, 1-butanol, and toluene). We demonstrate that while Sildenafil can adopt more than 100 metastable conformational configurations, only 12 are significantly populated across all of the environments investigated. Despite the complexity of the conformational space, we find that the most abundant conformers in solution are the closest to the conformers found in the most common Sildenafil crystal phase.

Crystal structure of adenylate kinase from an extremophilic archaeon Aeropyrum pernix with ATP and AMP

The crystal structure of an adenylate kinase from an extremophilic archaeon Aeropyrum pernix was determined in complex with full ligands, ATP-Mg2+ and AMP, at a resolution of 2.0 Å. The protein forms a trimer as found for other adenylate kinases from archaea. Interestingly, the reacting three atoms, two phosphorus and one oxygen atoms, were located almost in line, supporting the SN2 nucleophilic substitution reaction mechanism. Based on the crystal structure obtained, the reaction coordinate was estimated by the quantum mechanics calculations combined with molecular dynamics. It was found that the reaction undergoes two energy barriers; the steps for breaking the bond between the oxygen and γ-phosphorus atoms of ATP to produce a phosphoryl fragment and creating the bond between the phosphoryl fragment and the oxygen atom of the β-phosphate group of ADP. The reaction coordinate analysis also suggested the role of amino-acid residues for the catalysis of adenylate kinase.

Host-Induced Gene Silencing of an Adenylate Kinase Gene Involved in Fungal Energy Metabolism Improves Plant Resistance to Verticillium dahliae

Verticillium wilt, caused by the ascomycete fungus Verticillium dahliae (Vd), is a devastating disease of numerous plant species. However, the pathogenicity/virulence-related genes in this fungus, which may be potential targets for improving plant resistance, remain poorly elucidated. For the study of these genes in Vd, we used a well-established host-induced gene silencing (HIGS) approach and identified 16 candidate genes, including a putative adenylate kinase gene (VdAK). Transiently VdAK-silenced plants developed milder wilt symptoms than control plants did. VdAK-knockout mutants were more sensitive to abiotic stresses and had reduced germination and virulence on host plants. Transgenic Nicotiana benthamiana and Arabidopsis thaliana plants that overexpressed VdAK dsRNAs had improved Vd resistance than the wild-type. RT-qPCR results showed that VdAK was also crucial for energy metabolism. Importantly, in an analysis of total small RNAs from Vd strains isolated from the transgenic plants, a small interfering RNA (siRNA) targeting VdAK was identified in transgenic N. benthamiana. Our results demonstrate that HIGS is a promising strategy for efficiently screening pathogenicity/virulence-related genes of Vd and that VdAK is a potential target to control this fungus.

Adenylate Kinase: A Ubiquitous Enzyme Correlated with Medical Conditions

Adenylate kinase is a small, usually monomeric, enzyme found in every living thing due to its crucial role in energetic metabolism. This paper outlines the most relevant data about adenylate kinases isoforms, and the connection between dysregulation or mutation of human adenylate kinase and medical conditions. The following datadases were consulted: National Centre for Biotechnology Information, Protein Data Bank, and Mouse Genomic Informatics. The SmartBLAST tool, EMBOSS Needle Program, and Clustal Omega Program were used to analyze the best protein match, and to perform pairwise sequence alignment and multiple sequence alignment. Human adenylate kinase genes are located on different chromosomes, six of them being on the chromosomes 1 and 9. The adenylate kinases' intracellular localization and organ distribution explain their dysregulation in many diseases. The cytosolic isoenzyme 1 and the mitochondrial isoenzyme 2 are the main adenylate kinases that are integrated in the vast network of inflammatory modulators. The cytosolic isoenzyme 5 is correlated with limbic encephalitis and Leu673Pro mutation of the isoenzyme 7 leads to primary male infertility due to impairment of the ciliary function. The impairment of the mitochondrial isoenzymes 2 and 4 is demonstrated in neuroblastoma or glioma. The adenylate kinases are disease modifier that can assess the risk of diseases where oxidative stress plays a crucial role in pathogenesis like metabolic syndrome or neurodegenerative diseases. Because adenylate kinases has ATP as substrate, they are integrated in the global network of energetic process of any organism therefore are valid target for new pharmaceutical compounds.

Molecular Docking Evaluation of (E)-5-arylidene-2-thioxothiazolidin-4-one Derivatives as Selective Bacterial Adenylate Kinase Inhibitors

Multi-drug resistant microorganism infections with emerging problems that require not only a prevention strategy, but also the development of new inhibitory compounds. Six previously synthesized 5-arylidene-2-thioxothiazolidin-4-one derivatives 1a–f, were screened for inhibitory activity on adenylate kinases of different origins by molecular docking. The compounds 1c and 1d were the most efficient inhibitors of bacterial and some archean adenylate kinases. Hydrogen bond interactions were observed with the residues belonging to the ATP binding site. Moreover human adenylate kinases are poor targets, suggesting that this selectivity offers promising prospectives for refining the structure of our compounds.