How the Ligand-Induced Reorganization of Protein Internal Energies Is Coupled to Conformational Events

Conformational Dynamics, Energetics, and the Divergent Evolution of Allosteric Regulation: The Case of the Yeast MAPK Family

Allosteric mechanisms provide finely-tuned control over signalling proteins. Proteins of the same family may share high sequence identity and structural similarity but show distinct traits of allosteric control and evolutionary divergent regulation. Revealing the determinants of such properties may be important to understand the molecular bases of different regulatory pathways. Herein, we investigate whether and how evolutionarily-divergent traits of allosteric regulation in homologous proteins can be decoded in terms of internal dynamics and interaction networks that support functionally oriented conformations. In this framework, we start from the comparative analysis of the dynamics and energetics of the yeast MAP Kinases (MAPKs) Fus3 and Kss1 in their native basins. Importantly, distinctive dynamic and energetic stabilization features emerge, which can be related to the two proteins' differential ability to be phosphorylated and engage with the allosteric activator Ste5. We then expanded our study to other evolutionarily-related MAPKs. We show that the dynamical and energetical traits defining the distinct regulatory profiles of Fus3 and Kss1 can be traced along their evolutionary tree. Overall, our approach is able to reconnect (latent) allostery with the principal elements of protein structural stabilization and dynamics, showing how allosteric regulation was encrypted in MAPKs structure well before Ste5 appearance.

Computing allostery: from the understanding of biomolecular regulation and the discovery of cryptic sites to molecular design

The concept of allostery has become a central tenet in the study of biological systems. In parallel, the discovery of allosteric drugs is generating new opportunities to selectively modulate difficult targets involved in pathologic mechanisms. Molecular simulations can provide atomistically detailed insight into the processes involved in allosteric regulation and signaling, and at the same time, they have the potential to unveil regulatory hotspots or cryptic sites that are not immediately evident from the analysis of static structures. In this context, computational approaches should be able to connect the study of allosteric regulation at different scales to the possibility of leveraging this knowledge to expand the chemical space of new, active drugs. Here, we will discuss recent advances in the study of allosteric regulation via computational methods and connect the mechanistic knowledge generated to the possibility of designing new small molecules that can tweak the functions of their receptors.

Integrating Protein Interaction Surface Prediction with a Fragment-Based Drug Design: Automatic Design of New Leads with Fragments on Energy Surfaces

Protein-protein interactions (PPIs) have emerged in the past years as significant pharmacological targets in the development of new therapeutics due to their key roles in determining pathological pathways. Herein, we present fragments on energy surfaces, a simple and general design strategy that integrates the analysis of the dynamic and energetic signatures of proteins to unveil the substructures involved in PPIs, with docking, selection, and combination of drug-like fragments to generate new PPI inhibitor candidates. Specifically, structural representatives of the target protein are used as inputs for the blind physics-based prediction of potential protein interaction surfaces using the matrix of low coupling energy decomposition method. The predicted interaction surfaces are subdivided into overlapping windows that are used as templates to direct the docking and combination of fragments representative of moieties typically found in active drugs. This protocol is then applied and validated using structurally diverse, important PPI targets as test systems. We demonstrate that our approach facilitates the exploration of the molecular diversity space of potential ligands, with no requirement of prior information on the location and properties of interaction surfaces or on the structures of potential lead compounds. Importantly, the hit molecules that emerge from our ab initio design share high chemical similarity with experimentally tested active PPI inhibitors. We propose that the protocol we describe here represents a valuable means of generating initial leads against difficult targets for further development and refinement.

Large-Scale Ligand Perturbations of the Protein Conformational Landscape Reveal State-Specific Interaction Hotspots

Protein flexibility is important for ligand binding but often ignored in drug design. Considering proteins as ensembles rather than static snapshots creates opportunities to target dynamic proteins that lack FDA-approved drugs, such as the human chaperone, heat shock protein 90 (Hsp90). Hsp90α accommodates ligands with a dynamic lid domain, yet no comprehensive analysis relating lid conformations to ligand properties is available. To date, ∼300 ligand-bound Hsp90α crystal structures are deposited in the Protein Data Bank, which enables us to consider ligand binding as a perturbation of the protein conformational landscape. By estimating binding site volumes, we classified structures into distinct major and minor lid conformations. Supported by retrospective docking, each conformation creates unique hotspots that bind chemically distinguishable ligands. Clustering revealed insightful exceptions and the impact of crystal packing. Overall, Hsp90α’s plasticity provides a cautionary tale of overinterpreting individual crystal structures and motivates an ensemble-based view of drug design.

Ligand binding and global adaptation of the GlnPQ substrate binding domain 2 revealed by molecular dynamics simulations

Substrate-binding domains (SBD) are important structural elements of substrate transporters mediating the transport of essential molecules across the cell membrane. The SBD2 domain of the glutamine (GLN) transporter from bacteria consists of two domains D1 and D2 that bind GLN in the space between the domains in a closed conformation. In the absence of ligand, SBD2 adopts an open conformation with increased domain distance. In molecular dynamics (MD) simulations in the absence of ligands, no closing of the open conformation was observed on the MD time scale. Addition of GLN resulted in several reversible binding and unbinding events of GLN at the binding site on the D1 domain but did not induce domain closing indicating that binding and global domain closing do not occur simultaneously. The SBD2 structure remained in a closed state when starting from the GLN-bound closed crystal structure and opened quickly to reach the open state upon removal of the GLN ligand. Free energy simulations to induce opening to closing indicated a barrier for closing in the absence and presence of ligand and a significant penalty for closing without GLN in the binding pocket. Simulations of a Leu480Ala mutation also indicate that an interaction of a C-terminal D1-tail_471-484 with a D2-helix_418-427 (not contacting the substrate-binding region) plays a decisive role for controlling the barrier of conformational switching in the SBD2 protein. The results allow us to derive a model of the molecular mechanism of substrate binding to SBD2 and associated conformational changes.

Simple Model of Protein Energetics To Identify Ab Initio Folding Transitions from All-Atom MD Simulations of Proteins

A fundamental requirement to predict the native conformation, address questions of sequence design and optimization, and gain insights into the folding mechanisms of proteins lies in the definition of an unbiased reaction coordinate that reports on the folding state without the need to compare it to reference values, which might be unavailable for new (designed) sequences. Here, we introduce such a reaction coordinate, which does not depend on previous structural knowledge of the native state but relies solely on the energy partition within the protein: the spectral gap of the pair nonbonded energy matrix (ENergy Gap, ENG). This quantity can be simply calculated along unbiased MD trajectories. We show that upon folding the gap increases significantly, while its fluctuations are reduced to a minimum. This is consistently observed for a diverse set of systems and trajectories. Our approach allows one to promptly identify residues that belong to the folding core as well as residues involved in non-native contacts that need to be disrupted to guide polypeptides to the folded state. The energy gap and fluctuations criteria are then used to develop an automatic detection system which allows us to extract and analyze folding transitions from a generic MD trajectory. We speculate that our method can be used to detect conformational ensembles in dynamic and intrinsically disordered proteins, revealing potential preorganization for binding.

Chemical Perturbation of Oncogenic Protein Folding: from the Prediction of Locally Unstable Structures to the Design of Disruptors of Hsp90-Client Interactions

Protein folding quality control in cells requires the activity of a class of proteins known as molecular chaperones. Heat shock protein-90 (Hsp90), a multidomain ATP driven molecular machine, is a prime representative of this family of proteins. Interactions between Hsp90, its co-chaperones, and client proteins have been shown to be important in facilitating the correct folding and activation of clients. Hsp90 levels and functions are elevated in tumor cells. Here, we computationally predict the regions on the native structures of clients c-Abl, c-Src, Cdk4, B-Raf and Glucocorticoid Receptor, that have the highest probability of undergoing local unfolding, despite being ordered in their native structures. Such regions represent potential ideal interaction points with the Hsp90-system. We synthesize mimics spanning these regions and confirm their interaction with partners of the Hsp90 complex (Hsp90, Cdc37 and Aha1) by Nuclear Magnetic Resonance (NMR). Designed mimics selectively disrupt the association of their respective clients with the Hsp90 machinery, leaving unrelated clients unperturbed and causing apoptosis in cancer cells. Overall, selective targeting of Hsp90 protein-protein interactions is achieved without causing indiscriminate degradation of all clients, setting the stage for the development of therapeutics based on specific chaperone:client perturbation.

Impact of Mutations on NPAC Structural Dynamics: Mechanistic Insights from MD Simulations

NPAC is a cytokine-like nuclear factor involved in chromatin modification and regulation of gene expression. In humans, the C-terminal domain of NPAC has the conserved structure of the β-hydroxyacid dehydrogenases (β-HAD) protein superfamily, which forms a stable tetrameric core scaffold for demethylase enzymes and organizes multiple sites for chromatin interactions. In spite of the close structural resemblance to other β-HAD family members, the human NPAC dehydrogenase domain lacks a highly conserved catalytic lysine, substituted by a methionine. The reintroduction of the catalytic lysine by M437 K mutation results in a significant decrease of stability of the tetramer. Here, we have computationally investigated the molecular determinants of the functional differences between methionine and lysine-containing NPAC proteins. We find that the single mutation can determine strong consequences in terms of dynamics, stability, and ultimately ability to assemble in supramolecular complexes: the higher stability and lower flexibility of the methionine variant structurally preorganizes the monomer for tetramerization, whereas lysine increases flexibility and favors conformations that, while catalytically active, are not optimal for tetrameric assembly. We combine structure-dynamics analysis to an evolutionary study of NPAC sequences, showing that the methionine mutation occurs in a specifically flexible region of the lysine-containing protein, flanked by two domains that concentrate most of the stabilizing interactions. In our model, such separation of stability nuclei and flexible regions appears to favor the functional innovability of the protein.

Advances in the Computational Identification of Allosteric Sites and Pathways in Proteins

With the increasing difficulty to develop new drugs and the emergence of resistance to traditional orthosteric-site inhibitors, the search for alternatives is finally approaching the focus on allosteric sites. Allosteric sites offer opportunities to regulate many pharmacologically targeted pathways by inhibition or activation. In addition, allosteric sites tend to be less conserved than the functional site, which may facilitate the design of specific effectors in the protein families for which specific orthosteric inhibitors have proved difficult to design. Furthermore, recent evidence suggests that all proteins might be susceptible of allosteric regulation, increasing the space of druggable targets. Computational identification of allosteric sites has therefore become an active field of research. The problem can be approached from two sides: (1) the identification of allosteric-communication pathways between the functional site and potential allosteric sites and (2) the functional-site-independent identification of allosteric sites. While the first approach tends to be more laborious and thus restricted to a single protein, the second tends to be more amenable to larger-scale analysis, thus providing tools for the two drug discovery scenarios: the analysis of known targets and the screening for new potential targets. Here, I show some basic concepts and methods useful to the identification of allosteric sites and pathways, in line with these two approaches. I describe them in some detail to build a clear framework, at the risk of losing the interest of experts. Examples of recent studies involving these methods are also illustrated, focusing on the techniques rather than on their findings on allosterism.