Global Transitions of Proteins Explored by a Multiscale Hybrid Methodology: Application to Adenylate Kinase

Efficient and accurate mapping of transition pathways is a challenging problem in allosteric proteins. We propose here a to our knowledge new methodology called collective molecular dynamics (coMD). coMD takes advantage of the collective modes of motions encoded by the fold, simultaneously evaluating the interactions and energetics via a full-atomic MD simulation protocol. The basic approach is to deform the structure collectively along the modes predicted by the anisotropic network model, upon selecting them via a Monte Carlo/Metropolis algorithm from among the complete pool of all accessible modes. Application to adenylate kinase, an allosteric enzyme composed of three domains, CORE, LID, and NMP, shows that both open-to-closed and closed-to-open transitions are readily sampled by coMD, with large-scale motions of the LID dominating. An energy-barrier crossing occurs during the NMP movements. The energy barrier originates from a switch between the salt bridges K136-D118 at the LID-CORE interface and K57-E170 and D33-R156 at the CORE-NMP and LID-NMP interfaces, respectively. Despite its simplicity and computing efficiency, coMD yields ensembles of transition pathways in close accord with detailed full atomic simulations, lending support to its utility as a multiscale hybrid method for efficiently exploring the allosteric transitions of multidomain or multimeric proteins.     

Global Transitions of Proteins Explored by a Multiscale Hybrid Methodology: Application to Adenylate Kinase

Efficient and accurate mapping of transition pathways is a challenging problem in allosteric proteins. We propose here a to our knowledge new methodology called collective molecular dynamics (coMD). coMD takes advantage of the collective modes of motions encoded by the fold, simultaneously evaluating the interactions and energetics via a full-atomic MD simulation protocol. The basic approach is to deform the structure collectively along the modes predicted by the anisotropic network model, upon selecting them via a Monte Carlo/Metropolis algorithm from among the complete pool of all accessible modes. Application to adenylate kinase, an allosteric enzyme composed of three domains, CORE, LID, and NMP, shows that both open-to-closed and closed-to-open transitions are readily sampled by coMD, with large-scale motions of the LID dominating. An energy-barrier crossing occurs during the NMP movements. The energy barrier originates from a switch between the salt bridges K136-D118 at the LID-CORE interface and K57-E170 and D33-R156 at the CORE-NMP and LID-NMP interfaces, respectively. Despite its simplicity and computing efficiency, coMD yields ensembles of transition pathways in close accord with detailed full atomic simulations, lending support to its utility as a multiscale hybrid method for efficiently exploring the allosteric transitions of multidomain or multimeric proteins.     

Probing molecular mechanisms of the Hsp90 chaperone: biophysical modeling identifies key regulators of functional dynamics

Deciphering functional mechanisms of the Hsp90 chaperone machinery is an important objective in cancer biology aiming to facilitate discovery of targeted anti-cancer therapies. Despite significant advances in understanding structure and function of molecular chaperones, organizing molecular principles that control the relationship between conformational diversity and functional mechanisms of the Hsp90 activity lack a sufficient quantitative characterization. We combined molecular dynamics simulations, principal component analysis, the energy landscape model and structure-functional analysis of Hsp90 regulatory interactions to systematically investigate functional dynamics of the molecular chaperone. This approach has identified a network of conserved regions common to the Hsp90 chaperones that could play a universal role in coordinating functional dynamics, principal collective motions and allosteric signaling of Hsp90. We have found that these functional motifs may be utilized by the molecular chaperone machinery to act collectively as central regulators of Hsp90 dynamics and activity, including the inter-domain communications, control of ATP hydrolysis, and protein client binding. These findings have provided support to a long-standing assertion that allosteric regulation and catalysis may have emerged via common evolutionary routes. The interaction networks regulating functional motions of Hsp90 may be determined by the inherent structural architecture of the molecular chaperone. At the same time, the thermodynamics-based "conformational selection" of functional states is likely to be activated based on the nature of the binding partner. This mechanistic model of Hsp90 dynamics and function is consistent with the notion that allosteric networks orchestrating cooperative protein motions can be formed by evolutionary conserved and sparsely connected residue clusters. Hence, allosteric signaling through a small network of distantly connected residue clusters may be a rather general functional requirement encoded across molecular chaperones. The obtained insights may be useful in guiding discovery of allosteric Hsp90 inhibitors targeting protein interfaces with co-chaperones and protein binding clients.     

Structure-Based Network Analysis of Activation Mechanisms in the ErbB Family of Receptor Tyrosine Kinases: The Regulatory Spine Residues Are Global Mediators of Structural Stability and Allosteric Interactions

The ErbB protein tyrosine kinases are among the most important cell signaling families and mutation-induced modulation of their activity is associated with diverse functions in biological networks and human disease. We have combined molecular dynamics simulations of the ErbB kinases with the protein structure network modeling to characterize the reorganization of the residue interaction networks during conformational equilibrium changes in the normal and oncogenic forms. Structural stability and network analyses have identified local communities integrated around high centrality sites that correspond to the regulatory spine residues. This analysis has provided a quantitative insight to the mechanism of mutation-induced “superacceptor” activity in oncogenic EGFR dimers. We have found that kinase activation may be determined by allosteric interactions between modules of structurally stable residues that synchronize the dynamics in the nucleotide binding site and the αC-helix with the collective motions of the integrating αF-helix and the substrate binding site. The results of this study have pointed to a central role of the conserved His-Arg-Asp (HRD) motif in the catalytic loop and the Asp-Phe-Gly (DFG) motif as key mediators of structural stability and allosteric communications in the ErbB kinases. We have determined that residues that are indispensable for kinase regulation and catalysis often corresponded to the high centrality nodes within the protein structure network and could be distinguished by their unique network signatures. The optimal communication pathways are also controlled by these nodes and may ensure efficient allosteric signaling in the functional kinase state. Structure-based network analysis has quantified subtle effects of ATP binding on conformational dynamics and stability of the EGFR structures. Consistent with the NMR studies, we have found that nucleotide-induced modulation of the residue interaction networks is not limited to the ATP site, and may enhance allosteric cooperativity with the substrate binding region by increasing communication capabilities of mediating residues.     

Coupling between global dynamics and signal transduction pathways: a mechanism of allostery for chaperonin GroEL

Despite significant efforts toward understanding the molecular basis of allosteric communication, the mechanisms by which local energetic and conformational changes cooperatively diffuse from ligand-binding sites to distal regions across the 3-dimensional structure of allosteric proteins remain to be established. Recent experimental and theoretical evidence supports the view that allosteric communication is facilitated by the intrinsic ability of the biomolecules to undergo collective changes in structure, triggered by ligand binding. Two groups of studies recently proved to provide insights into such intrinsic, structure-induced effects: elastic network models that permit us to visualize the cooperative changes in conformation that are most readily accessible near native state conditions, and information-theoretic approaches that elucidate the most efficient pathways of signal transmission favored by the overall architecture. Using a combination of these two approaches, we highlight, by way of application to the bacterial chaperonin complex GroEL-GroES, how the most cooperative modes of motion play a role in mediating the propagation of allosteric signals. A functional coupling between the global dynamics sampled under equilibrium conditions and the signal transduction pathways inherently favored by network topology appears to control allosteric effects.     

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.     

Maturation dynamics of bacteriophage HK97 capsid

Maturation of the bacteriophage HK97 capsid requires a large conformational change of the virus capsid. Experimental studies have identified several intermediates along this maturation pathway. To gain insights into the molecular mechanisms of capsid maturation, we examined the fluctuation dynamics of the procapsid and mature capsid using a residue-level computational approach. The most cooperative motions of the procapsid are found to be consistent with the observed change in configuration that takes place during maturation. A few dominant modes of motion are sufficient to describe the anisotropic expansion that accompanies maturation. Based upon these modes, maturation is proposed to occur via an overall expansion and reconfiguration of the capsid initiated by puckering of the pentamers, followed by flattening and crosslinking of the hexameric subunits, and finally crosslinking of the pentameric subunits. The highly mobile E loops are stabilized by anchoring to highly stable residues belonging to neighboring subunits.     

Structural dynamics of nucleosome core particle: comparison with nucleosomes containing histone variants

The present study provides insights on the dominant mechanisms of motions of the nucleosome core particle and the changes in its functional dynamics in response to histone variants. Comparative analysis of the global dynamics of nucleosomes with native and variant H2A histones, using normal mode analysis revealed that the dynamics of the nucleosome is highly symmetric, and its interaction with the nucleosomal DNA plays a vital role in its regulation. The collective dynamics of nucleosomes are predicted to be dominated by two types of large-scale motions: (1) a global stretching-compression of nucleosome along the dyad axis by which the nucleosome undergoes a breathing motion with a massive distortion of nucleosomal DNA, modulated by histone-DNA interactions; and (2) the flipping (or bending) of both the sides of the nucleosome in an out-of-plane fashion with respect to the dyad axis, originated by the highly dynamic N-termini of H3 and (H2A.Z-H2B) dimer in agreement with the experimentally observed perturbed dynamics of the particular N-terminus under physiological conditions. In general, the nucleosomes with variant histones exhibit higher mobilities and weaker correlations between internal motions compared to the nucleosome containing ordinary histones. The differences are more pronounced at the L1 and L2 loops of the respective monomers H2B and H2A, and at the N-termini of the monomers H3 and H4, all of which closely interact with the wrapping DNA.     

Differential Modulation of Functional Dynamics and Allosteric Interactions in the Hsp90-Cochaperone Complexes with p23 and Aha1: A Computational Study

Allosteric interactions of the molecular chaperone Hsp90 with a large cohort of cochaperones and client proteins allow for molecular communication and event coupling in signal transduction networks. The integration of cochaperones into the Hsp90 system is driven by the regulatory mechanisms that modulate the progression of the ATPase cycle and control the recruitment of the Hsp90 clientele. In this work, we report the results of computational modeling of allosteric regulation in the Hsp90 complexes with the cochaperones p23 and Aha1. By integrating protein docking, biophysical simulations, modeling of allosteric communications, protein structure network analysis and the energy landscape theory we have investigated dynamics and stability of the Hsp90-p23 and Hsp90-Aha1 interactions in direct comparison with the extensive body of structural and functional experiments. The results have revealed that functional dynamics and allosteric interactions of Hsp90 can be selectively modulated by these cochaperones via specific targeting of the regulatory hinge regions that could restrict collective motions and stabilize specific chaperone conformations. The protein structure network parameters have quantified the effects of cochaperones on conformational stability of the Hsp90 complexes and identified dynamically stable communities of residues that can contribute to the strengthening of allosteric interactions. According to our results, p23-mediated changes in the Hsp90 interactions may provide “molecular brakes” that could slow down an efficient transmission of the inter-domain allosteric signals, consistent with the functional role of p23 in partially inhibiting the ATPase cycle. Unlike p23, Aha1-mediated acceleration of the Hsp90-ATPase cycle may be achieved via modulation of the equilibrium motions that facilitate allosteric changes favoring a closed dimerized form of Hsp90. The results of our study have shown that Aha1 and p23 can modulate the Hsp90-ATPase activity and direct the chaperone cycle by exerting the precise control over structural stability, global movements and allosteric communications in Hsp90.     

Distinct allosteric pathways in imidazole glycerol phosphate synthase from yeast and bacteria

Understanding the relationship between protein structures and their function is still an open question that becomes very challenging when allostery plays an important functional role. Allosteric proteins, in fact, exploit different ranges of motions (from sidechain local fluctuations to long-range collective motions) to effectively couple distant binding sites, and of particular interest is whether allosteric proteins of the same families with similar functions and structures also necessarily share the same allosteric mechanisms. Here, we compared the early dynamics initiating the allosteric communication of a prototypical allosteric enzyme from two different organisms, i.e., the imidazole glycerol phosphate synthase (IGPS) enzymes from the thermophilic bacteria and the yeast, working at high and room temperatures, respectively. By combining molecular dynamics simulations and network models derived from graph theory, we found rather distinct early allosteric dynamics in the IGPS from the two organisms, involving significatively different allosteric pathways in terms of both local and collective motions. Given the successful prediction of key allosteric residues in the bacterial IGPS, whose mutation disrupts its allosteric communication, the outcome of this study paves the way for future experimental studies on the yeast IGPS that could foster therapeutic applications by exploiting the control of IGPS enzyme allostery.     

Weak Intra-Ring Allosteric Communications of the Archaeal Chaperonin Thermosome Revealed by Normal Mode Analysis

Chaperonins are molecular machines that use ATP-driven cycles to assist misfolded substrate proteins to reach the native state. During the functional cycle, these machines adopt distinct nucleotide-dependent conformational states, which reflect large-scale allosteric changes in individual subunits. Distinct allosteric kinetics has been described for the two chaperonin classes. Bacterial (group I) chaperonins, such as GroEL, undergo concerted subunit motions within each ring, whereas archaeal and eukaryotic chaperonins (group II) undergo sequential subunit motions. We study these distinct mechanisms through a comparative normal mode analysis of monomer and double-ring structures of the archaeal chaperonin thermosome and GroEL. We find that thermosome monomers of each type exhibit common low-frequency behavior of normal modes. The observed distinct higher-frequency modes are attributed to functional specialization of these subunit types. The thermosome double-ring structure has larger contribution from higher-frequency modes, as it is found in the GroEL case. We find that long-range intersubunit correlation of amino-acid pairs is weaker in the thermosome ring than in GroEL. Overall, our results indicate that distinct allosteric behavior of the two chaperonin classes originates from different wiring of individual subunits as well as of the intersubunit communications.     

Simulating nanoscale functional motions of biomolecules

We are describing efficient dynamics simulation methods for the characterization of functional motion of biomolecules on the nanometer scale. Multivariate statistical methods are widely used to extract and enhance functional collective motions from molecular dynamics (MD) simulations. A dimension reduction in MD is often realized through a principal component analysis (PCA) or a singular value decomposition (SVD) of the trajectory. Normal mode analysis (NMA) is a related collective coordinate space approach, which involves the decomposition of the motion into vibration modes based on an elastic model. Using the myosin motor protein as an example we describe a hybrid technique termed amplified collective motions (ACM) that enhances sampling of conformational space through a combination of normal modes with atomic level MD. Unfortunately, the forced orthogonalization of modes in collective coordinate space leads to complex dependencies that are not necessarily consistent with the symmetry of biological macromolecules and assemblies. In many biological molecules, such as HIV-1 protease, reflective or rotational symmetries are present that are broken using standard orthogonal basis functions. We present a method to compute the plane of reflective symmetry or the axis of rotational symmetry from the trajectory frames. Moreover, we develop an SVD that best approximates the given trajectory while respecting the symmetry. Finally, we describe a local feature analysis (LFA) to construct a topographic representation of functional dynamics in terms of local features. The LFA representations are low-dimensional, and provide a reduced basis set for collective motions, but unlike global collective modes they are sparsely distributed and spatially localized. This yields a more reliable assignment of essential dynamics modes across different MD time windows.     

THz time scale structural rearrangements and binding modes in lysozyme-ligand interactions

Predicting the conformational changes in proteins that are relevant for substrate binding is an ongoing challenge in the aim of elucidating the functional states of proteins. The motions that are induced by protein-ligand interactions are governed by the protein global modes. Our measurements indicate that the detected changes in the global backbone motion of the enzyme upon binding reflect a shift from the large-scale collective dominant mode in the unbound state towards a functional twisting deformation that assists in closing the binding cleft. Correlated motion in lysozyme has been implicated in enzyme function in previous studies, but detailed characterization of the internal fluctuations that enable the protein to explore the ensemble of conformations that ultimately foster large-scale conformational change is yet unknown. For this reason, we use THz spectroscopy to investigate the picosecond time scale binding modes and collective structural rearrangements that take place in hen egg white lysozyme (HEWL) when bound by the inhibitor (NAG) ₃. These protein thermal motions correspond to fluctuations that have a role in both selecting and sampling from the available protein intrinsic conformations that communicate function. Hence, investigation of these fast, collective modes may provide knowledge about the mechanism leading to the preferred binding process in HEWL-(NAG) ₃. Specifically, in this work we find that the picosecond time scale hydrogen-bonding rearrangements taking place in the protein hydration shell with binding modify the packing density within the hydrophobic core on a local level. These localized, intramolecular contact variations within the protein core appear to facilitate the large cooperative movements within the interfacial region separating the α- and β- domain that mediate binding. The THz time-scale fluctuations identified in the protein-ligand system may also reveal a molecular mechanism for substrate recognition.     

Insights into the evolution of allosteric properties. The NADH binding site of hexameric type II citrate synthases

Study of the hexameric and allosterically regulated citrate synthases (type II CS) provides a rare opportunity to gain not only an understanding of a novel allosteric mechanism but also insight into how such properties can evolve from an unregulated structural platform (the dimeric type I CS). To address both of these issues, we have determined the structure of the complex of NADH (a negative allosteric effector) with the F383A variant of type II Escherichia coli CS. This variant was chosen because its kinetics indicate it is primarily in the T or inactive allosteric conformation, the state that strongly binds to NADH. Our structural analyses show that the six NADH binding sites in the hexameric CS complex are located at the interfaces between dimer units such that most of each site is formed by one subunit, but a number of key residues are drawn from the adjacent dimer. This arrangement of interactions serves to explain why NADH allosteric regulation is a feature only of hexameric type II CS. Surprisingly, in both the wild-type enzyme and the NADH complex, the two subunits of each dimer within the hexameric conformation are similar but not identical in structure, and therefore, while the general characteristics of NADH binding interactions are similar in each subunit, the details of these are somewhat different between subunits. Detailed examination of the observed NADH binding sites indicates that both direct charged interactions and the overall cationic nature of the sites are likely responsible for the ability of these sites to discriminate between NADH and NAD(+). A particularly novel characteristic of the complex is the horseshoe conformation assumed by NADH, which is strikingly different from the extended conformation found in its complexes with most proteins. Sequence homology studies suggest that this approach to binding NADH may arise out of the evolutionary need to add an allosteric regulatory function to the base CS structure. Comparisons of the amino acid sequences of known type II CS enzymes, from different Gram-negative bacteria taxonomic groups, show that the NADH-binding residues identified in our structure are strongly conserved, while hexameric CS molecules that are insensitive to NADH have undergone key changes in the sequence of this part of the protein.     

The biological functions of low-frequency vibrations (phonons) 5. A phenomenological theory

Low-frequency internal motions of a biomacromolecule are thought to possess significant biological function from the dynamic point of view. In this paper, a general phenomenological theory is established by which it is clearly verified that low-frequency resonance plays a central role in the energy transmission required during the cooperative interaction between subunits in a protein oligomer. According to the present theory, it is found that the energy transmission between a pair of diagonal subunits in a protein oligomer with a polygon arrangement is the most efficient, so as to in a sense further predict that after a ligand is bound to a subunit by random collision, its diagonal subunit in the same protein oligomer will possess the greatest probability of binding with the next ligand. Furthermore, based on the concept of the 'resonance-controlled trigger' derived from the phenomenological theory, it is feasible to estimate the lower time limit of allosteric transition from one subunit to the other. Such a time limit depends on the dominant low-frequency mode of each subunit, the ratio of the coupling force constant to the corresponding inherent force constant, as well as the geometrical arrangement of subunits in a protein oligomer. So far none of the allosteric transitions observed in proteins has exceeded the time limit as defined here, indicating a logical consistency between our theory and the experiments.     

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.