Discovery of intramolecular signal transduction network based on a new protein dynamics model of energy dissipation

A novel approach to reveal intramolecular signal transduction network is proposed in this work. To this end, a new algorithm of network construction is developed, which is based on a new protein dynamics model of energy dissipation. A key feature of this approach is that direction information is specified after inferring protein residue-residue interaction network involved in the process of signal transduction. This enables fundamental analysis of the regulation hierarchy and identification of regulation hubs of the signaling network. A well-studied allosteric enzyme, E. coli aspartokinase III, is used as a model system to demonstrate the new method. Comparison with experimental results shows that the new approach is able to predict all the sites that have been experimentally proved to desensitize allosteric regulation of the enzyme. In addition, the signal transduction network shows a clear preference for specific structural regions, secondary structural types and residue conservation. Occurrence of super-hubs in the network indicates that allosteric regulation tends to gather residues with high connection ability to collectively facilitate the signaling process. Furthermore, a new parameter of propagation coefficient is defined to determine the propagation capability of residues within a signal transduction network. In conclusion, the new approach is useful for fundamental understanding of the process of intramolecular signal transduction and thus has significant impact on rational design of novel allosteric proteins.     

Dynamic energy correlation analysis of E. coli aspartokinase III and alteration of allosteric regulation by manipulating energy transduction pathways

Conformational change associated with allosteric regulation in a protein is ultimately driven by energy transformation. However, little is known about the latter process. In this work, we combined steered molecular dynamics simulations and sequence conservation analysis to investigate the conformational changes and energy transformation in the allosteric enzyme aspartokinase III (AK III) from Escherichia coli. Correlation analysis of energy change at residue level indicated significant transformation between electrostatic energy and dihedral angle energy during the allosteric regulation. Key amino acid residues located in the corresponding energy transduction pathways were identified by dynamic energy correlation analysis. To verify their functions, residues with a high energy correlation in the pathways were altered and their effects on allosteric regulation of AKIII were determined. This study sheds new insights into energy transformation during allosteric regulation of AK III and proposes a strategy to identify key residues that are involved in intramolecular energy transduction and thus in driving the allosteric process.     

Integrating molecular dynamics and co-evolutionary analysis for reliable target prediction and deregulation of the allosteric inhibition of aspartokinase for amino acid production

Deregulation of allosteric inhibition of enzymes is a challenge for strain engineering and has been achieved so far primarily by random mutation and trial-and-error. In this work, we used aspartokinase, an important allosteric enzyme for industrial amino acids production, to demonstrate a predictive approach that combines protein dynamics and evolution for a rational reengineering of enzyme allostery. Molecular dynamic simulation of aspartokinase III (AK3) from Escherichia coli and statistical coupling analysis of protein sequences of the aspartokinase family allowed to identify a cluster of residues which are correlated during protein motion and coupled during the evolution. This cluster of residues forms an interconnected network mediating the allosteric regulation, including most of the previously reported positions mutated in feedback insensitive AK3 mutants. Beyond these mutation positions, we have successfully constructed another twelve targeted mutations of AK3 desensitized toward lysine inhibition. Six threonine-insensitive mutants of aspartokinase I-homoserine dehydrogenase I (AK1-HD1) were also created based on the predictions. The proposed approach can be widely applied for the deregulation of other allosteric enzymes.     

The N-terminal end of the catalytic domain of SRC kinase Hck is a conformational switch implicated in long-range allosteric regulation

Signal transduction in cell growth and proliferation involves regulation of kinases through long-range allostery between remote protein regions. Molecular dynamics free energy calculations are used to clarify the coupling between the catalytic domain of Src kinase Hck and its N-terminal end connecting to the regulatory SH2 and SH3 modules. The N-terminal end is stable in the orientation required for the regulatory modules to remain properly bound only in the inactive catalytic domain. In the active catalytic domain, the N-terminal end prefers a different conformation consistent with dissociation of the regulatory modules. The free energy surface shows that the N-terminal end acts as a reversible two-state conformational switch coupling the catalytic domain to the regulatory modules. Structural analogy with insulin receptor kinase and c-Src suggests that such reversible conformational switching in a critical hinge region could be a common mechanism in long-range allosteric regulation of protein kinase activity.     

Change in allosteric network affects binding affinities of PDZ domains: analysis through perturbation response scanning

The allosteric mechanism plays a key role in cellular functions of several PDZ domain proteins (PDZs) and is directly linked to pharmaceutical applications; however, it is a challenge to elaborate the nature and extent of these allosteric interactions. One solution to this problem is to explore the dynamics of PDZs, which may provide insights about how intramolecular communication occurs within a single domain. Here, we develop an advancement of perturbation response scanning (PRS) that couples elastic network models with linear response theory (LRT) to predict key residues in allosteric transitions of the two most studied PDZs (PSD-95 PDZ3 domain and hPTP1E PDZ2 domain). With PRS, we first identify the residues that give the highest mean square fluctuation response upon perturbing the binding sites. Strikingly, we observe that the residues with the highest mean square fluctuation response agree with experimentally determined residues involved in allosteric transitions. Second, we construct the allosteric pathways by linking the residues giving the same directional response upon perturbation of the binding sites. The predicted intramolecular communication pathways reveal that PSD-95 and hPTP1E have different pathways through the dynamic coupling of different residue pairs. Moreover, our analysis provides a molecular understanding of experimentally observed hidden allostery of PSD-95. We show that removing the distal third alpha helix from the binding site alters the allosteric pathway and decreases the binding affinity. Overall, these results indicate that (i) dynamics plays a key role in allosteric regulations of PDZs, (ii) the local changes in the residue interactions can lead to significant changes in the dynamics of allosteric regulations, and (iii) this might be the mechanism that each PDZ uses to tailor their binding specificities regulation.     

Signaling pathways of PDZ2 domain: a molecular dynamics interaction correlation analysis

PDZ domains are found in many signaling proteins. One of their functions is to provide scaffolds for forming membrane-associated protein complexes by binding to the carboxyl termini of their partners. PDZ domains are thought also to play a signal transduction role by propagating the information that binding has occurred to remote sites. In this study, a molecular dynamics (MD) simulation-based approach, referred to as an interaction correlation analysis, is applied to the PDZ2 domain to identify the possible signal transduction pathways. A residue correlation matrix is constructed from the interaction energy correlations between all residue pairs obtained from the MD simulations. Two continuous interaction pathways, starting at the ligand binding pocket, are identified by a hierarchical clustering analysis of the residue correlation matrix. One pathway is mainly localized at the N-terminal side of helix alpha1 and the adjacent C-terminus of loop beta1-beta2. The other pathway is perpendicular to the central beta-sheet and extends toward the side of PDZ2 domain opposite to the ligand binding pocket. The results complement previous studies based on multiple sequence analysis, NMR, and MD simulations. Importantly, they reveal the energetic origin of the long-range coupling. The PDZ2 results, as well as the earlier rhodopsin analysis, show that the interaction correlation analysis is a robust approach for determining pathways of intramolecular signal transduction.     

Interesting structural and dynamical behaviors exhibited by the AF-6 PDZ domain upon Bcr peptide binding

PDZ (postsynaptic density-95, disks large, zonula occludens-1) domains are small, protein-protein interaction modules that have multiple binding surfaces for the docking of diverse molecules. These domains can propagate signals from ligand-binding site to distal regions of the structure through allosteric communication. Recent works have revealed that picosecond to nanosecond time scale dynamics play a potential role in propagating long-range signals within a protein. Comparison of AF-6 PDZ domain structures in free and complex forms shows a conformation rearrangement of distal surface 2, which is far from the peptide binding groove. The relaxation dispersion experiments detected that the free AF-6 PDZ domain was sampling multiple conformations; millisecond dynamics mapped a network for allostery signal transmission throughout the AF-6 PDZ domain in the weak saturation state, and intramolecular motions were observed in distal surface 1 when the protein was saturated. These results provide evidence that the allosteric process in the AF-6 PDZ domain is not two-state; instead, the millisecond dynamic network provides a mechanism for the transmission of allosteric signals throughout a protein. Interestingly, the two distal surfaces of the AF-6 PDZ domain respond differently to peptide binding; distal surface 1 changes in millisecond dynamics, whereas distal surface 2 undergoes structural rearrangement. The significance of the different response patterns in the signaling pathway and its relevance to the function of the AF-6 PDZ domain should be studied further.     

Energy transduction optical sensor in skeletal myosin

The skeletal myosin cross-bridge in dynamic association with actin is the unitary energy transducer in muscle, converting free energy from ATP hydrolysis into contractile force. Myosin's conserved ATP-sensitive tryptophan (AST) is an energy transduction optical sensor signaling transduction-related transient conformation change by modulating its fluorescence intensity amplitude and relaxation rate. Recently introduced techniques have provided the means of observing the time-resolved intensity decay from this single residue in the native protein to elucidate the mechanism of its ATP sensitivity. AST signal characteristics could be derived from local protein structure by a scenario involving interactions with excited-state tryptophan. This investigation suggests the very different possibility that hypochromism induced in the tryptophan absorption band, a ground-state effect, is a significant structural effector of optical transduction sensing. This possibility makes feasible the interpretation of the transient AST optical signal in terms of dynamical protein structure, thereby raising the empirical signal to the level of a structural determinant. Using the crystallographically based geometry from several myosin structures, the maximum calculated AST hypochromism is <10% to be compared with the value of approximately 30% observed here experimentally. Rationalizing the discrepancy invites further investigation of S1 dynamical structure local to the AST during transduction.     

Nucleotide sequence of lysC gene encoding the lysine-sensitive aspartokinase III of Escherichia coli K12. Evolutionary pathway leading to three isofunctional enzymes

The lysC gene encoding the lysine-sensitive aspartokinase III of Escherichia coli K12 has been cloned and its nucleotide sequence determined. Analysis of the deduced protein sequence (449 amino acid residues) reveals that the entire sequence of aspartokinase III is homologous to the N-terminal part of the two iso- and bifunctional aspartokinase-homoserine dehydrogenases I and II of E. coli. An evolutionary pathway leading to the three molecular species present in the same organism is proposed, and the possible involvement of a highly conserved region in subunit interactions is discussed.     

Ultrafast dynamics-driven biomolecular recognition where fast activities dictate slow events

In general, biological macromolecules require significant dynamical freedom to carry out their different functions, including signal transduction, metabolism, catalysis and gene regulation. Effectors (ligands, DNA and external milieu, etc) are considered to function in a purely dynamical manner by selectively stabilizing a specific dynamical state, thereby regulating biological function. In particular, proteins in presence of these effectors can exist in several dynamical states with distinct binding or enzymatic activity. Here, we have reviewed the efficacy of ultrafast fluorescence spectroscopy to monitor the dynamical flexibility of various proteins in presence of different effectors leading to their biological activity. Recent studies demonstrate the potency of a combined approach involving picosecond-resolved Förster resonance energy transfer, polarisation-gated fluorescence and time-dependent stokes shift for the exploration of ultrafast dynamics in biomolecular recognition of various protein molecules. The allosteric protein–protein recognition following differential protein–DNA interaction is shown to be a consequence of some ultrafast segmental motions at the C-terminal of Gal repressor protein dimer with DNA operator sequences O_E and O_I. Differential ultrafast dynamics at the C-terminal of λ-repressor protein with two different operator DNA sequences for the protein–protein interaction with different strengths is also reviewed. We have also systemically briefed the study on the role of ultrafast dynamics of water molecules on the functionality of enzyme proteins α-chymotrypsin and deoxyribonuclease I. The studies on the essential ultrafast dynamics at the active site of the enzyme α-chymotrypsin by using an anthraniloyl fluorescent extrinsic probe covalently attached to the serine-195 residue for the enzymatic activity at homeothermic condition has also been reviewed. Finally, we have highlighted the evidence that a photoinduced dynamical event dictates the molecular recognition of a photochromic ligand, dihydroindolizine with the serine protease α-chymotrypsin and with a liposome (L-α-phosphatidylcholine).     

A systems biology perspective on protein structural dynamics and signal transduction

The functional dynamics of signal transduction through protein interaction networks are determined both by network topology and by the signal processing properties of component proteins. In order to understand the emergent properties of signal transduction networks in terms of information processing, storage and decision making, we not only need to map the so-called 'interactome' but, perhaps more importantly, we also have to understand how the structural dynamics of constituent proteins shape non-linear responses through cooperativity and allostery. Several in silico methods have been developed to identify networks of cooperative residues in proteins and help infer their mode of action. Applying this type of analysis to important classes of modular signal transduction domains should, in principle, allow the function of these proteins to be abstracted in terms of their information processing characteristics, permitting better comprehension of the systemic properties of biological networks.     

Dynamical probing of allosteric control in nuclear receptors

The dynamical behavior of the nuclear receptor LXR/RXR heterodimer was investigated with molecular dynamics simulations. The simulations reveal correlated motion between residues across the dimer interface that depends significantly on occupation of the ligand binding sites of the monomers. These results are broadly consistent with the observed experimental behavior of the dimers, where structural perturbation is thought to be a key element in signal transduction. Our results provide dynamical support for this model of allosteric control.     

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.     

ATPase Subdomain IA Is a Mediator of Interdomain Allostery in Hsp70 Molecular Chaperones

The versatile functions of the heat shock protein 70 (Hsp70) family of molecular chaperones rely on allosteric interactions between their nucleotide-binding and substrate-binding domains, NBD and SBD. Understanding the mechanism of interdomain allostery is essential to rational design of Hsp70 modulators. Yet, despite significant progress in recent years, how the two Hsp70 domains regulate each other''s activity remains elusive. Covariance data from experiments and computations emerged in recent years as valuable sources of information towards gaining insights into the molecular events that mediate allostery. In the present study, conservation and covariance properties derived from both sequence and structural dynamics data are integrated with results from Perturbation Response Scanning and in vivo functional assays, so as to establish the dynamical basis of interdomain signal transduction in Hsp70s. Our study highlights the critical roles of SBD residues D481 and T417 in mediating the coupled motions of the two domains, as well as that of G506 in enabling the movements of the α-helical lid with respect to the β-sandwich. It also draws attention to the distinctive role of the NBD subdomains: Subdomain IA acts as a key mediator of signal transduction between the ATP- and substrate-binding sites, this function being achieved by a cascade of interactions predominantly involving conserved residues such as V139, D148, R167 and K155. Subdomain IIA, on the other hand, is distinguished by strong coevolutionary signals (with the SBD) exhibited by a series of residues (D211, E217, L219, T383) implicated in DnaJ recognition. The occurrence of coevolving residues at the DnaJ recognition region parallels the behavior recently observed at the nucleotide-exchange-factor recognition region of subdomain IIB. These findings suggest that Hsp70 tends to adapt to co-chaperone recognition and activity via coevolving residues, whereas interdomain allostery, critical to chaperoning, is robustly enabled by conserved interactions.     

Molecular dynamics simulations of aptamer-binding reveal generalized allostery in thrombin

Thrombin is an attractive target for antithrombotic therapy due to its central role in thrombosis and hemostasis as well as its role in inducing tumor growth, metastasis, and tumor invasion. The thrombin-binding DNA aptamer (TBA), is under investigation for anticoagulant drugs. Although aptamer binding experiments have been revealed various effects on thrombin’s enzymatic activities, the detailed picture of the thrombin’s allostery from TBA binding is still unclear. To investigate thrombin’s response to the aptamer-binding at the molecular level, we compare the mechanical properties and free energy landscapes of the free and aptamer-bound thrombin using microsecond-scale all-atom GPU-based molecular dynamics simulations. Our calculations on residue fluctuations and coupling illustrate the allosteric effects of aptamer-binding at the atomic level, highlighting the exosite II, 60s, γ and the sodium loops, and the alpha helix region in the light chains involved in the allosteric changes. This level of details clarifies the mechanisms of previous experimentally demonstrated phenomena, and provides a prediction of the reduced autolysis rate after aptamer-binding. The shifts in thrombin’s ensemble of conformations and free energy surfaces after aptamer-binding demonstrate that the presence of bound-aptamer restricts the conformational freedom of thrombin suggesting that conformational selection, i.e. generalized allostery, is the dominant mechanism of thrombin-aptamer binding. The profound perturbation on thrombin’s mechanical and thermodynamic properties due to the aptamer-binding, which was revealed comprehensively as a generalized allostery in this work, may be exploited in further drug discovery and development.     

Ligand-dependent dynamics and intramolecular signaling in a PDZ domain

Allosteric communication is a fundamental process that proteins use to propagate signals from one site to functionally important distal sites. Although allostery is usually associated with multimeric proteins and enzymes, “long-range” communication may be a fundamental property of proteins. In some cases, communication occurs with minimal structural change. PDZ (post-synaptic density-95/discs large/zonula occludens-1) domains are small, protein-protein binding modules that can use multiple surfaces for docking diverse molecules. Furthermore, these domains have long-range energetic couplings that link the ligand-binding site to distal regions of the structure. Here, we show that allosteric behavior in a representative member of the PDZ domain family may be directly detected using side-chain methyl dynamics measurements. The changes in side-chain dynamics parameters in the second PDZ domain from the human tyrosine phosphatase 1E (hPTP1E) were determined upon binding a peptide target. Long-range dynamic effects were detected that correspond to previously observed pair-wise energetic couplings. These results provide one of the first experimental examples for the potential role of ps-ns timescale dynamics in propagating long-range signals within a protein, and reinforce the idea that dynamic fluctuations in proteins contribute to allosteric signal transduction.     

Analysis of two-component signal transduction by mathematical modeling using the KdpD/KdpE system of Escherichia coli

A mathematical model for the KdpD/KdpE two-component system is presented and its dynamical behavior is analyzed. KdpD and KdpE regulate expression of the kdpFABC operon encoding the high affinity K+ uptake system KdpFABC of Escherichia coli. The model is validated in a two step procedure: (i) the elements of the signal transduction part are reconstructed in vitro. Experiments with the purified sensor kinase and response regulator in presence or absence of DNA fragments comprising the response regulator binding-site are performed. (ii) The mRNA and molecule number of KdpFABC are determined in vivo at various extracellular K+ concentrations. Based on the identified parameters for the in vitro system it is shown, that different time hierarchies appear which are used for model reduction. Then the model is transformed in such a way that a singular perturbation problem is formulated. The analysis of the in vivo system shows that the model can be separated into two parts (submodels which are called functional units) that are connected only in a unidirectional way. Hereby one submodel represents signal transduction while the second submodel describes the gene expression.