Allostery and conformational free energy changes in human tryptophanyl‐tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of allostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side‐chains and then to compute the ligand‐induced population shift. Finally, we obtain the free‐energy landscape of the protein in equilibrium, characterizing the free‐energy minima accessed by the protein complexes. We have chosen human tryptophanyl‐tRNA synthetase (hTrpRS), a protein responsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Stability and binding properties of a modified thrombin binding aptamer

Aptamer-based drugs represent an attractive approach in pharmacological therapy. The most studied aptamer, thrombin binding aptamer (TBA), folds into a well-defined quadruplex structure and binds to its target with good specificity and affinity. Modified aptamers with improved biophysical properties could constitute a new class of therapeutic aptamers. In this study we show that the modified thrombin binding aptamer (mTBA), ^3′GGT^5′-^5′TGGTGTGGTTGG^3′, which also folds into a quadruplex structure, is more stable than its unmodified counterpart and shows a higher thrombin affinity. The stability of the modified aptamer was investigated using differential scanning calorimetry, and the energetics of mTBA and TBA binding to thrombin was characterized by means of isothermal titration calorimetry (ITC). ITC data revealed that TBA/thrombin and mTBA/thrombin binding stoichiometry is 1:2 for both interactions. Structural models of the two complexes of thrombin with TBA and with mTBA were also obtained and subjected to molecular dynamics simulations in explicit water. Analysis of the models led to an improvement of the understanding of the aptamer-thrombin recognition at a molecular level.     

Allostery in pharmacology: thermodynamics, evolution and design

This review focuses on basic models of allostery, the ambiguous application of the allosteric term in pharmacology illustrated by receptors, the role of thermodynamics in allosteric mechanisms, evolution and design of allostery. The initial step of ligand activation is closure of the agonist-binding cavity. Large entropy increases accompany the agonist-elicited conformational changes of pentameric ligand-gated ion channels due to cavity closure and rearrangement of transmembrane helices. The effects of point mutations on thermodynamic parameters of binding and function can reveal energetic coupling of neighbouring (and distant) amino acid residues in activation. High-order double-mutant cycle analysis and rate-equilibrium linear free-energy relationships can identify the trajectory and conformational spread of activation.Protein assembly and allostery can be deduced from colocalization and physicochemical principles. Molecular evolution has led from homooligomerization of protomers to heterotropic cooperativity and to allosteric regulation. Examples are discussed such as similar paths of protein (dis)assembly and evolution, irreversible evolution, statistical analysis of sequence homology revealing coevolution, different impacts of adaptation and evolution on hemoglobin, and the flagellar motor switch of bacteria. The driving force of dynamic allostery is associated with funnel-like free energy landscapes of protein binding and shifts in conformational fluctuations upon binding. Allostery can be designed based on our increasing knowledge of natural allosteric mechanisms and evolution. The allosteric principle has been applied for various bio/macro/molecular and signal transduction systems as well as in cognitive sciences.     

The molecular basis of thrombin allostery revealed by a 1.8 A structure of the "slow" form

Thrombin participates in its own positive and negative feedback loops, and its allosteric state helps determine the hemostatic balance. Here we present the 1.8 A crystallographic structure of S195A thrombin in two conformational states: active site occupied and active site free. The active site-occupied form shows how thrombin can accommodate substrates, such as protein C. The active site-free form is in a previously unobserved closed conformation of thrombin, which satisfies all the conditions of the so-called "slow" form. A mechanism of allostery is revealed, which relies on the concerted movement of the disulphide bond between Cys168 and 182 and aromatic residues Phe227, Trp215, and Trp60d. These residues constitute an allosteric switch, which is flipped directly through sodium binding, resulting in the fast form with an open active site.     

Exploring Mechanisms of Allosteric Regulation and Communication Switching in the Multiprotein Regulatory Complexes of the Hsp90 Chaperone with Cochaperones and Client Proteins: Atomistic Insights from Integrative Biophysical Modeling and Network Analysis of Conformational Landscapes

Understanding molecular principles underlying Hsp90 chaperone functions and modulation of client activity is fundamental to dissect activation mechanisms of many proteins. In this work, we performed a computational investigation of the Hsp90-Hsp70-Hop-CR client complex to examine allosteric regulatory mechanisms underlying dynamic chaperone interactions and principles of chaperone-dependent client recognition and remodeling. Conformational dynamics analysis using high-resolution coarse-grained simulations and ensemble-based local frustration analysis suggest that the Hsp90 chaperone could recognize and recruit the GR client by invoking reciprocal dynamic exchanges near the intermolecular interfaces with the client. Using mutational scanning of the intermolecular residues in the Hsp90-Hsp70-Hop-GR complex, we identified binding energy hotspots in the regulatory complex. Perturbation-based network analysis and dynamic fluctuations-based modeling of allosteric residue potentials are employed for a detailed analysis of allosteric interaction networks and identification of conformational communication switches. We found that allosteric interactions between the Hsp90, the client-bound Hsp70 and Hop cochaperone can define two allosteric residue clusters that control client recruitment in which the intrinsic Hsp70 allostery is exploited to mediate integration of the Hsp70-bound client into the Hsp90 chaperone system. The results suggest a model of dynamics-driven allostery that enables efficient client recruitment and loading through allosteric couplings between intermolecular interfaces and communication switch centers. This study showed that the Hsp90 interactions with client proteins may operate under dynamic-based allostery in which ensembles of preexisting conformational states and intrinsic allosteric pathways present in the Hsp90 and Hsp70 chaperones can be exploited for recognition and integration of substrate proteins.     

Protein dynamics and allostery: an NMR view

Allostery, the process by which distant sites within a protein system are energetically coupled, is an efficient and ubiquitous mechanism for activity regulation. A purely mechanical view of allostery invoking only structural changes has developed over the decades as the classical view of the phenomenon. However, a fast growing list of examples illustrate the intimate link between internal motions over a wide range of time scales and function in protein-ligand interactions. Proteins respond to perturbations by redistributing their motions and they use fluctuating conformational states for binding and conformational entropy as a carrier of allosteric energy to modulate association with ligands. In several cases allosteric interactions proceed with minimal or no structural changes. We discuss emerging paradigms for the central role of protein dynamics in allostery.     

Structure-Based Statistical Mechanical Model Accounts for the Causality and Energetics of Allosteric Communication

Allostery is one of the pervasive mechanisms through which proteins in living systems carry out enzymatic activity, cell signaling, and metabolism control. Effective modeling of the protein function regulation requires a synthesis of the thermodynamic and structural views of allostery. We present here a structure-based statistical mechanical model of allostery, allowing one to observe causality of communication between regulatory and functional sites, and to estimate per residue free energy changes. Based on the consideration of ligand free and ligand bound systems in the context of a harmonic model, corresponding sets of characteristic normal modes are obtained and used as inputs for an allosteric potential. This potential quantifies the mean work exerted on a residue due to the local motion of its neighbors. Subsequently, in a statistical mechanical framework the entropic contribution to allosteric free energy of a residue is directly calculated from the comparison of conformational ensembles in the ligand free and ligand bound systems. As a result, this method provides a systematic approach for analyzing the energetics of allosteric communication based on a single structure. The feasibility of the approach was tested on a variety of allosteric proteins, heterogeneous in terms of size, topology and degree of oligomerization. The allosteric free energy calculations show the diversity of ways and complexity of scenarios existing in the phenomenology of allosteric causality and communication. The presented model is a step forward in developing the computational techniques aimed at detecting allosteric sites and obtaining the discriminative power between agonistic and antagonistic effectors, which are among the major goals in allosteric drug design.     

An ensemble view of thrombin allostery

Abstract Thrombin is the central protease of the coagulation cascade. Its activity is tightly regulated to ensure rapid blood clotting while preventing uncontrolled thrombosis. Thrombin interacts with multiple substrates and cofactors and is critically involved in both pro- and anticoagulant pathways of the coagulation network. Its allosteric regulation, especially by the monovalent cation Na+, has been the focus of research for more than 30 years. It is believed that thrombin can adopt an anticoagulant ('slow') conformation and, after Na+ binding, a structurally distinct procoagulant ('fast') state. In the past few years, however, the general view of allostery has evolved from one of rigid structural changes towards thermodynamic ensembles of conformational states. With this background, the view of the allosteric regulation of thrombin has also changed. The static view of the two-state model has been dismissed in favor of a more dynamic view of thrombin allostery. Herein, we review recent data that demonstrate that apo-thrombin is zymogen-like and exists as an ensemble of conformations. Furthermore, we describe how ligand binding to thrombin allosterically stabilizes conformations on the continuum from zymogen to protease.     

Protein mechanics: How force regulates molecular function

Background

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

Scope

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

Major conclusions

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

General significance

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

The role of dynamics in allosteric regulation

The biomolecular conformational changes often associated with allostery are, by definition, dynamic processes. Recent publications have disclosed the role of pre-existing equilibria of conformational substates in this process. In addition, the role of dynamics as an entropic carrier of free energy of allostery has been investigated. Recent work thus shows that dynamics is pivotal to allostery, and that it constitutes much more than just the move from the 'T'-state to the 'R'-state. Emerging computational studies have described the actual pathways of allosteric change.     

An electrostatic network and long-range regulation of Src kinases

The regulatory mechanism of Src tyrosine kinases includes conformational activation by a change in the catalytic domain tertiary structure and in domain-domain contacts between the catalytic domain and the SH2/SH3 regulatory domains. The kinase is activated when tyrosine phosphorylation occurs on the activation loop, but without phosphorylation of the C-terminal tail. Activation also occurs by allostery when contacts between the catalytic domain (CD) and the regulatory SH3 and SH2 domains are released as a result of exogenous protein binding. The aim of this work is to examine the proposed role of an electrostatic network in the conformational transition and to elucidate the molecular mechanism for long-range, allosteric conformational activation by using a combination of experimental enzyme kinetics and nonequilibrium molecular dynamics simulations. Salt dependence of the induction phase is observed in kinetic assays and supports the role of an electrostatic network in the transition. In addition, simulations provide evidence that allosteric activation involves a concerted motion coupling highly conserved residues, and spanning several nanometers from the catalytic site to the regulatory domain interface to communicate between the CD and the regulatory domains.     

Emerging Methods and Applications to Decrypt Allostery in Proteins and Nucleic Acids

Many large protein-nucleic acid complexes exhibit allosteric regulation. In these systems, the propagation of the allosteric signaling is strongly coupled to conformational dynamics and catalytic function, challenging state-of-the-art analytical methods. Here, we review established and innovative approaches used to elucidate allosteric mechanisms in these complexes. Specifically, we report network models derived from graph theory and centrality analyses in combination with molecular dynamics (MD) simulations, introducing novel schemes that implement the synergistic use of graph theory with enhanced simulations methods and ab-initio MD. Accelerated MD simulations are used to construct “enhanced network models”, describing the allosteric response over long timescales and capturing the relation between allostery and conformational changes. “Ab-initio network models” combine graph theory with ab-initio MD and quantum mechanics/molecular mechanics (QM/MM) simulations to describe the allosteric regulation of catalysis by following the step-by-step dynamics of biochemical reactions. This approach characterizes how the allosteric regulation changes from reactants to products and how it affects the transition state, revealing a tense-to-relaxed allosteric regulation along the chemical step. Allosteric models and applications are showcased for three paradigmatic examples of allostery in protein-nucleic acid complexes: (i) the nucleosome core particle, (ii) the CRISPR-Cas9 genome editing system and (iii) the spliceosome. These methods and applications create innovative protocols to determine allosteric mechanisms in protein-nucleic acid complexes that show tremendous promise for medicine and bioengineering.     

Allosteric Effect of Adenosine Triphosphate on Peptide Recognition by 3′5′-Cyclic Adenosine Monophosphate Dependent Protein Kinase Catalytic Subunits

The allosteric influence of adenosine triphosphate (ATP) on the binding effectiveness of a series of peptide inhibitors with the catalytic subunit of 3′5′-cyclic adenosine monophosphate dependent protein kinase was investigated, and the dependence of this effect on peptide structure was analyzed. The allosteric effect was calculated as ratio of peptide binding effectiveness with the enzyme-ATP complex and with the free enzyme, quantified by the competitive inhibition of the enzyme in the presence of ATP excess, and by the enzyme-peptide complex denaturation assay, respectively It was found that the principle “better binding—stronger allostery” holds for interactions of the studied peptides with the enzyme, indicating that allostery and peptide binding with the free enzyme are governed by the same specificity pattern. This means that the allosteric regulation does not include new ligand–protein interactions, but changes the intensity (strength) of the interatomic forces that govern the complex formation in the case of each individual ligand. We propose that the allosteric regulation can be explained by the alteration of the intrinsic dynamics of the protein by ligand binding, and that this phenomenon, in turn, modulates the ligand off-rate from its binding site as well as the binding affinity. The positive allostery could therefore be induced by a reduction in the enzyme’s overall intrinsic dynamics.     

Recombinant FimH Adhesin Demonstrates How the Allosteric Catch Bond Mechanism Can Support Fast and Strong Bacterial Attachment in the Absence of Shear

The FimH protein of Escherichia coli is a model two-domain adhesin that is able to mediate an allosteric catch bond mechanism of bacterial cell attachment, where the mannose-binding lectin domain switches from an ‘inactive’ conformation with fast binding to mannose to an ‘active’ conformation with slow detachment from mannose. Because mechanical tensile force favors separation of the domains and, thus, FimH activation, it has been thought that the catch bonds can only be manifested in a fluidic shear-dependent mode of adhesion. Here, we used recombinant FimH variants with a weakened inter-domain interaction and show that a fast and sustained allosteric activation of FimH can also occur under static, non-shear conditions. Moreover, it appears that lectin domain conformational activation happens intrinsically at a constant rate, independently from its ability to interact with the pilin domain or mannose. However, the latter two factors control the rate of FimH deactivation. Thus, the allosteric catch bond mechanism can be a much broader phenomenon involved in both fast and strong cell-pathogen attachments under a broad range of hydrodynamic conditions. This concept that allostery can enable more effective receptor-ligand interactions is fundamentally different from the conventional wisdom that allostery provides a mechanism to turn binding off under specific conditions.     

Ligand-induced stabilization of the aptamer terminal helix in the add adenine riboswitch

Riboswitches are structured mRNA elements that modulate gene expression. They undergo conformational changes triggered by highly specific interactions with sensed metabolites. Among the structural rearrangements engaged by riboswitches, the forming and melting of the aptamer terminal helix, the so-called P1 stem, is essential for genetic control. The structural mechanisms by which this conformational change is modulated upon ligand binding mostly remain to be elucidated. Here, we used pulling molecular dynamics simulations to study the thermodynamics of the P1 stem in the add adenine riboswitch. The P1 ligand-dependent stabilization was quantified in terms of free energy and compared with thermodynamic data. This comparison suggests a model for the aptamer folding in which direct P1-ligand interactions play a minor role on the conformational switch when compared with those related to the ligand-induced aptamer preorganization.     

Investigation of interaction of thrombin-binding aptamer with thrombin and prethrombin-2 by simulation of molecular dynamics

Thrombin is a major component of blood clotting and involved in the formation of a fibrin clot. One of the precursors during thrombin maturation is prethrombin-2, with the presence of Arg363-Ile364 bond between the light and heavy chain of protein, the only distinction from thrombin. Prethrombin-2 is able to interact with less efficiency with a 15-mer thrombin-binding aptamer (TBA). We investigate the interaction of both known conformers of TBA with thrombin and prethrombin-2 by simulation of molecular dynamics. It was shown that TBA could interact with thrombin in both conformations with similar efficiency, although a stable complex of prethrombin-2 with TBA was found only in conformation identical with the aptamer structure, pdb 1HAO. Analysis of molecular dynamics of complexes offered an assumption that the motion of the exosite-1 forming loop Lys428-Ile438 determined the difference in affinity of the complexes of TBA with thrombin and prethrombin-2.     

The mechanism of high affinity pentasaccharide binding to antithrombin,insights from Gaussian accelerated molecular dynamics simulations

Abstract

The activity of antithrombin (AT), a serpin protease inhibitor, is enhanced by heparin and heparin analogs against its target proteases, mainly thrombin, factors Xa and IXa. Considerable amount of information is available on the multistep mechanism of the heparin pentasaccharide binding and conformational activation. However, much of the details were inferred from ‘static’ structures obtained by X-ray diffraction. Moreover, limited information is available for the early steps of binding mechanism other than kinetic studies with various ligands. To gain insights into these processes, we performed enhanced sampling molecular dynamics (MD) simulations using the Gaussian Accelerated Molecular Dynamics (GAMD) method, applied previously in drug binding studies. We were able to observe the binding of the pentasaccharide idraparinux to a ‘non-activated’ AT conformation in two separate trajectories with low root mean square deviation (RMSD) values compared to X-ray structures of the bound state. These trajectories along with further simulations of the AT-pentasaccharide complex provided insights into the mechanisms of multiple conformational transitions, including the expulsion of the hinge region, the extension of helix D and the conformational behavior of the reactive center loop (RCL). We could also confirm the high stability of helix P in non-activated AT conformations, such states might play an important role in heparin binding. ‘Generalized correlation’ matrices revealed possible paths of allosteric signal propagation to the binding sites for the target proteases, factors Xa and IXa. Enhanced MD simulations of ligand binding to AT may assist the design of new anticoagulant drugs.

Communicated by Ramaswamy H. Sarma     

The changing landscape of protein allostery

It is becoming increasingly clear that the fundamental capacity to undergo conformational change in response to ligand binding is intrinsic to proteins. This property confers on proteins the ability to be allosterically modulated in order to shift substrate binding affinities, alter enzymatic activity or regulate protein-protein interaction. How this allosteric modulation occurs--the pathways of communication, the shifting of conformational ensembles and the altered molecular dynamics--has received considerable attention during the past two years. Recent progress has helped outline the molecular origins of allostery in proteins as diverse as Hsp70 molecular chaperones and signal integrating proteins, such as WASP. In addition, allosteric properties have been successfully engineered into proteins for drug design or the development of novel biosensors. Methodological advances have provided exciting prospects for new insights and new biological roles of allosteric systems have been uncovered.     

Role of water-bridged interactions in metal ion coupled protein allostery

Allosteric communication between distant parts of proteins controls many cellular functions, in which metal ions are widely utilized as effectors to trigger the allosteric cascade. Due to the involvement of strong coordination interactions, the energy landscape dictating the metal ion binding is intrinsically rugged. How metal ions achieve fast binding by overcoming the landscape ruggedness and thereby efficiently mediate protein allostery is elusive. By performing molecular dynamics simulations for the Ca²⁺ binding mediated allostery of the calmodulin (CaM) domains, each containing two Ca²⁺ binding helix-loop-helix motifs (EF-hands), we revealed the key role of water-bridged interactions in Ca²⁺ binding and protein allostery. The bridging water molecules between Ca²⁺ and binding residue reduces the ruggedness of ligand exchange landscape by acting as a lubricant, facilitating the Ca²⁺ coupled protein allostery. Calcium-induced rotation of the helices in the EF-hands, with the hydrophobic core serving as the pivot, leads to exposure of hydrophobic sites for target binding. Intriguingly, despite being structurally similar, the response of the two symmetrically arranged EF-hands upon Ca²⁺ binding is asymmetric. Breakage of symmetry is needed for efficient allosteric communication between the EF-hands. The key roles that water molecules play in driving allosteric transitions are likely to be general in other metal ion mediated protein allostery.