Protein functional landscapes, dynamics, allostery: a tortuous path towards a universal theoretical framework

Illuminating Protein Allostery by Chemically Accurate Contact Response Analysis (ChACRA)

Decoding allostery at the atomic level is essential for understanding the relationship between a protein's sequence, structure, and dynamics. Recently, we have shown that decomposing temperature responses of inter-residue contacts can reveal allosteric couplings and provide useful insight into the functional dynamics of proteins. The details of this Chemically Accurate Contact Response Analysis (ChACRA) are presented here along with its application to two well-known allosteric proteins. The first protein, IGPS, is a model of ensemble allostery that lacks clear structural differences between the active and inactive states. We show that the application of ChACRA reveals the experimentally identified allosteric coupling between effector and active sites of IGPS. The second protein, ATCase, is a classic example of allostery with distinct active and inactive structural states. Using ChACRA, we directly identify the most significant residue level interactions underlying the enzyme's cooperative behavior. Both test cases demonstrate the utility of ChACRA's unsupervised machine learning approach for dissecting allostery at the residue level.

Reconnaissance of Allostery via the Restoration of Native p53 DNA-Binding Domain Dynamics in Y220C Mutant p53 Tumor Suppressor Protein

Allosteric regulation of protein dynamics infers a long-range deliberate propagation of information via micro- and macroscale interactions. The Y220C structural mutant is one of the most frequent cancerous p53 mutants. The mutation is distally located from the DNA-binding site of the p53 DNA-binding domain yet causes changes in DNA recognition. This system presents a unique opportunity to examine the allosteric control of mutated proteins under a drug design paradigm. We focus on the key case study of p53 Y220C mutation restoration by a series of new compounds suggested to have Y220C reactivation properties in comparison to our previous findings on the restorative potential of PK11000, a compound studied extensively for reactivation in vitro and in vivo. Previously, we implemented all-atom molecular dynamics (MD) simulations and our lab's techniques of MD-Sectors and MD-Markov state models on the wild type, the Y220C mutant, and Y220C with PK11000 to characterize the effector's restorative properties in terms of conformational dynamics and hydrogen bonding. In this study, we turn to probing the effects made by docking the battery of a new but less well-tested set of aminobenzothiazole derivative compounds reported by Baud et al., which show promise of Y220C rescue. We find that while complete and precise reconstitution of p53 WT molecular dynamics may not be observed as was the case with PK11000, dispersed local reconstitution of loop dynamics provides evidence of rescuing effects by aminobenzothiazole derivative N,2-dihydroxy-3,5-diiodo-4-(1H-pyrrol-1-yl)benzamide, Effector 22, like what we observed for PK11000. Generalizable insights into the mutation and allosteric reactivation of p53 by various effectors by reconstitution of WT dynamics observed in statistical conformational ensemble analysis and network inference are discussed, considering the development of allosteric drug design rooted in first principles.

Intermediate-aided allostery mechanism for α-glucosidase by Xanthene-11v as an inhibitor using residue interaction network analysis

Exploring allosteric inhibition and the discovery of new inhibitor binding sites are important studies in protein regulation mechanisms and drug discovery. Structural and network-based analyses of trajectories resulting from molecular dynamics (MD) simulations have been developed to discover protein dynamics, landscape, functions, and allosteric regions. Here, an experimentally suggested non-competitive inhibitor, xanthene-11v, was considered to explore its allosteric inhibition mechanism in α-glucosidase MAL12. Comparative structural and network analyses were applied to eight 250 ns independent MD simulations, four of which were performed in the free state and four of which were performed in ligand-bound forms. Projected two-dimensional free energy landscapes (FEL) were constructed from the probabilistic distribution of conformations along the first two principal components. The post-simulation analyses of the coordinates, side-chain torsion angles, non-covalent interaction networks, network communities, and their centralities were performed on α-glucosidase conformations and the intermediate sub-states. Important communities of residues have been found that connect the allosteric site to the active site. Some of these residues like Thr307, Arg312, TYR344, ILE345, Phe357, Asp406, Val407, Asp408, and Leu436 are the key messengers in the transition pathway between allosteric and active sites. Evaluating the probability distribution of distances between gate residues including Val407 in one community and Phe158, and Pro65 in another community depicted the closure of this gate due to the inhibitor binding. Six macro states of protein were deduced from the topology of FEL and analysis of conformational preference of free and ligand-bound systems to these macro states shows a combination of lock-and-key, conformational selection, and induced fit mechanisms are effective in ligand binding. All these results reveal structural states, allosteric mechanisms, and key players in the inhibition pathway of α-glucosidase by xanthene-11v.

Locating dynamic contributions to allostery via determining rates of vibrational energy transfer

Determining rates of energy transfer across non-covalent contacts for different states of a protein can provide information about dynamic and associated entropy changes during transitions between states. We investigate the relationship between rates of energy transfer across polar and nonpolar contacts and contact dynamics for the β₂-adrenergic receptor, a rhodopsin-like G-protein coupled receptor, in an antagonist-bound inactive state and agonist-bound active state. From structures sampled during molecular dynamics (MD) simulations, we find the active state to have, on average, a lower packing density, corresponding to generally more flexibility and greater entropy than the inactive state. Energy exchange networks (EENs) are computed for the inactive and active states from the results of the MD simulations. From the EENs, changes in the rates of energy transfer across polar and nonpolar contacts are found for contacts that remain largely intact during activation. Change in dynamics of the contact, and entropy associated with the dynamics, can be estimated from the change in rates of energy transfer across the contacts. Measurement of change in the rates of energy transfer before and after the transition between states thereby provides information about dynamic contributions to activation and allostery.

Variable Regions of p53 Isoforms Allosterically Hard Code DNA Interaction

Allosteric regulation of protein activity pervades biology as the "second secret of life." We have been examining the allosteric regulation and mutant reactivation of the tumor suppressor protein p53. We have found that generalizing the definition of allosteric effector to include entire proteins and expanding the meaning of binding site to include the interface of a transcription factor with its DNA to be useful in understanding the modulation of protein activity. Here, we cast the variable regions of p53 isoforms as allosteric regulators of p53 interactions with its consensus DNA. We implemented molecular dynamics simulations and our lab's new techniques of molecular dynamics (MD) sectors and MD-Markov state models to investigate the effects of nine naturally occurring splice variant isoforms of p53. We find that all of the isoforms differ from wild type in their dynamic properties and how they interact with the DNA. We consider the implications of these findings on allostery and cancer treatment.

Molecular Dynamics Simulations Establish the Molecular Basis for the Broad Allostery Hotspot Distributions in the Tetracycline Repressor

It is imperative to identify the network of residues essential to the allosteric coupling for the purpose of rationally engineering allostery in proteins. Deep mutational scanning analysis has emerged as a function-centric approach for identifying such allostery hotspots in a comprehensive and unbiased fashion, leading to observations that challenge our understanding of allostery at the molecular level. Specifically, a recent deep mutational scanning study of the tetracycline repressor (TetR) revealed an unexpectedly broad distribution of allostery hotspots throughout the protein structure. Using extensive molecular dynamics simulations (up to 50 μs) and free energy computations, we establish the molecular and energetic basis for the strong anticooperativity between the ligand and DNA binding sites. The computed free energy landscapes in different ligation states illustrate that allostery in TetR is well described by a conformational selection model, in which the apo state samples a broad set of conformations, and specific ones are selectively stabilized by either ligand or DNA binding. By examining a range of structural and dynamic properties of residues at both local and global scales, we observe that various analyses capture different subsets of experimentally identified hotspots, suggesting that these residues modulate allostery in distinct ways. These results motivate the development of a thermodynamic model that qualitatively explains the broad distribution of hotspot residues and their distinct features in molecular dynamics simulations. The multifaceted strategy that we establish here for hotspot evaluations and our insights into their mechanistic contributions are useful for modulating protein allostery in mechanistic and engineering studies.

The effect of dimerization and ligand binding on the dynamics of Kaposi's sarcoma-associated herpesvirus protease

The Kaposi's sarcoma-associated herpesvirus protease is essential for virus maturation. This protease functions under allosteric regulation that establishes its enzymatic activity upon dimerization. It exists in equilibrium between an inactive monomeric state and an active, weakly associating, dimeric state that is stabilized upon ligand binding. The dynamics of the protease dimer and its monomer were studied using the Gaussian network model and the anisotropic network model , and its role in mediating the allosteric regulation is demonstrated. We show that the dimer is composed of five dynamical domains. The central domain is formed upon dimerization and composed of helix five of each monomer, in addition to proximal and distal domains of each monomer. Dimerization reduces the mobility of the central domains and increases the mobility of the distal domains, in particular the binding site within them. The three slowest ANM modes of the dimer assist the protease in ligand binding, motion of the conserved Arg142 and Arg143 toward the oxyanion, and reducing the activation barrier for the tetrahedral transition state by stretching the bond that is cleaved by the protease. In addition, we show that ligand binding reduces the motion of helices α1 and α5 at the interface and explain how ligand binding can stabilize the dimer.

Quantifying Energetic and Entropic Pathways in Molecular Systems

When examining dynamics occurring at nonzero temperatures, both energy and entropy must be taken into account to describe activated barrier crossing events. Furthermore, good reaction coordinates need to be constructed to describe different metastable states and the transition mechanisms between them. Here we use a physics-based machine learning method called state predictive information bottleneck (SPIB) to find nonlinear reaction coordinates for three systems of varying complexity. SPIB is able to correctly predict an entropic bottleneck for an analytical flat-energy double-well system and identify the entropy- and energy-dominated pathways for an analytical four-well system. Finally, for a simulation of benzoic acid permeation through a lipid bilayer, SPIB is able to discover the entropic and energetic barriers to the permeation process. Given these results, we thus establish that SPIB is a reasonable and robust method for finding the important entropy, energy, and enthalpy barriers in physical systems, which can then be used to enhance the understanding and sampling of different activated mechanisms.

Identifying the leading dynamics of ubiquitin: A comparison between the tICA and the LE4PD slow fluctuations in amino acids' position

Molecular Dynamics (MD) simulations of proteins implicitly contain the information connecting the atomistic molecular structure and proteins' biologically relevant motion, where large-scale fluctuations are deemed to guide folding and function. In the complex multiscale processes described by MD trajectories, it is difficult to identify, separate, and study those large-scale fluctuations. This problem can be formulated as the need to identify a small number of collective variables that guide the slow kinetic processes. The most promising method among the ones used to study the slow leading processes in proteins' dynamics is the time-structure based on time-lagged independent component analysis (tICA), which identifies the dominant components in a noisy signal. Recently, we developed an anisotropic Langevin approach for the dynamics of proteins, called the anisotropic Langevin Equation for Protein Dynamics or LE4PD-XYZ. This approach partitions the protein's MD dynamics into mostly uncorrelated, wavelength-dependent, diffusive modes. It associates with each mode a free-energy map, where one measures the spatial extension and the time evolution of the mode-dependent, slow dynamical fluctuations. Here, we compare the tICA modes' predictions with the collective LE4PD-XYZ modes. We observe that the two methods consistently identify the nature and extension of the slowest fluctuation processes. The tICA separates the leading processes in a smaller number of slow modes than the LE4PD does. The LE4PD provides time-dependent information at short times and a formal connection to the physics of the kinetic processes that are missing in the pure statistical analysis of tICA.

Physics of biomolecular recognition and conformational dynamics

Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.

Coarse-Grained Modeling of Multiple Pathways in Conformational Transitions of Multi-Domain Proteins

Large-scale conformational transitions in multi-domain proteins are often essential for their functions. To investigate the transitions, it is necessary to explore multiple potential pathways, which involve different intermediate structures. Here, we present a multi-basin (MB) coarse-grained (CG) structure-based Go̅ model for describing transitions in proteins with more than two moving domains. This model is an extension of our dual-basin Go̅ model in which system-dependent parameters are determined systematically using the multistate Bennett acceptance ratio method. In the MB Go̅ model for multi-domain proteins, we assume that intermediate structures may have partial inter-domain native contacts. This approach allows us to search multiple transition pathways that involve distinct intermediate structures using the CG molecular dynamics (MD) simulations. We apply this scheme to an enzyme, adenylate kinase (AdK), which has three major domains and can move along two different pathways. Using the optimized mixing parameters for each pathway, AdK shows frequent transitions between the Open, Closed, and the intermediate basins and samples a wide variety of conformations within each basin. The explored multiple transition pathways could be compared with experimental data and examined in more detail by atomistic MD simulations.

Frustration and Direct-Coupling Analyses to Predict Formation and Function of Adeno-Associated Virus

Adeno-associated virus (AAV) is a promising gene therapy vector because of its efficient gene delivery and relatively mild immunogenicity. To improve delivery target specificity, researchers use combinatorial and rational library design strategies to generate novel AAV capsid variants. These approaches frequently propose high proportions of nonforming or noninfective capsid protein sequences that reduce the effective depth of synthesized vector DNA libraries, thereby raising the discovery cost of novel vectors. We evaluated two computational techniques for their ability to estimate the impact of residue mutations on AAV capsid protein-protein interactions and thus predict changes in vector fitness, reasoning that these approaches might inform the design of functionally enriched AAV libraries and accelerate therapeutic candidate identification. The Frustratometer computes an energy function derived from the energy landscape theory of protein folding. Direct-coupling analysis (DCA) is a statistical framework that captures residue coevolution within proteins. We applied the Frustratometer to select candidate protein residues predicted to favor assembled or disassembled capsid states, then predicted mutation effects at these sites using the Frustratometer and DCA. Capsid mutants were experimentally assessed for changes in virus formation, stability, and transduction ability. The Frustratometer-based metric showed a counterintuitive correlation with viral stability, whereas a DCA-derived metric was highly correlated with virus transduction ability in the small population of residues studied. Our results suggest that coevolutionary models may be able to elucidate complex capsid residue-residue interaction networks essential for viral function, but further study is needed to understand the relationship between protein energy simulations and viral capsid metastability.

Nano-Surveillance: Tracking Individual Molecules in a Sea of Chromatin

Chromatin is the epigenomic platform for diverse nuclear processes such as DNA repair, replication, transcription, telomere, and centromere function. In cancer cells, mutations in key processes result in DNA amplification, chromosome translocations, and chromothripsis, severely distorting the natural chromatin state. In normal and diseased states, dozens of chromatin effectors alter the physical integrity and dynamics of chromatin at the level of both single nucleosomes and arrays of nucleosomes folded into 3-dimensional shapes. Integrating these length scales, from the 10 nm sized nucleosome to mitotic chromosomes, whilst jostling within the crowded environment of the cell, cannot yet be achieved by a single technology. In this review, we discuss tools that have proven powerful in the investigation of nucleosome and chromatin fiber dynamics. We also provide a deeper focus into atomic force microscopy (AFM) applications that can bridge diverse length and time scales. Using time course AFM, we observe that chromatin condensation by H1.5 is dynamic, whereas using nano-indentation force spectroscopy we observe that both histone variants and nucleosome binding partners alter material properties of individual nucleosomes. Finally, we demonstrate how high-speed AFM can visualize plasmid DNA dynamics, intermittent nucleosome-nucleosome contacts, and changes in nucleosome phasing along a contiguous chromatin fiber. Altogether, the development of innovative technologies holds the promise of revealing the secret lives of nucleosomes, potentially bridging the gaps in our understanding of how chromatin works within living cells and tissues.

Allosteric Regulation at the Crossroads of New Technologies: Multiscale Modeling, Networks, and Machine Learning

Allosteric regulation is a common mechanism employed by complex biomolecular systems for regulation of activity and adaptability in the cellular environment, serving as an effective molecular tool for cellular communication. As an intrinsic but elusive property, allostery is a ubiquitous phenomenon where binding or disturbing of a distal site in a protein can functionally control its activity and is considered as the "second secret of life." The fundamental biological importance and complexity of these processes require a multi-faceted platform of synergistically integrated approaches for prediction and characterization of allosteric functional states, atomistic reconstruction of allosteric regulatory mechanisms and discovery of allosteric modulators. The unifying theme and overarching goal of allosteric regulation studies in recent years have been integration between emerging experiment and computational approaches and technologies to advance quantitative characterization of allosteric mechanisms in proteins. Despite significant advances, the quantitative characterization and reliable prediction of functional allosteric states, interactions, and mechanisms continue to present highly challenging problems in the field. In this review, we discuss simulation-based multiscale approaches, experiment-informed Markovian models, and network modeling of allostery and information-theoretical approaches that can describe the thermodynamics and hierarchy allosteric states and the molecular basis of allosteric mechanisms. The wealth of structural and functional information along with diversity and complexity of allosteric mechanisms in therapeutically important protein families have provided a well-suited platform for development of data-driven research strategies. Data-centric integration of chemistry, biology and computer science using artificial intelligence technologies has gained a significant momentum and at the forefront of many cross-disciplinary efforts. We discuss new developments in the machine learning field and the emergence of deep learning and deep reinforcement learning applications in modeling of molecular mechanisms and allosteric proteins. The experiment-guided integrated approaches empowered by recent advances in multiscale modeling, network science, and machine learning can lead to more reliable prediction of allosteric regulatory mechanisms and discovery of allosteric modulators for therapeutically important protein targets.

Exploring Conformational Space with Thermal Fluctuations Obtained by Normal-Mode Analysis

Proteins in their native states can be represented as ensembles of conformers in dynamical equilibrium. Thermal fluctuations are responsible for transitions between these conformers. Normal-modes analysis (NMA) using elastic network models (ENMs) provides an efficient procedure to explore global dynamics of proteins commonly associated with conformational transitions. In the present work, we present an iterative approach to explore protein conformational spaces by introducing structural distortions according to their equilibrium dynamics at room temperature. The approach can be used either to perform unbiased explorations of conformational space or to explore guided pathways connecting two different conformations, e.g., apo and holo forms. In order to test its performance, four proteins with different magnitudes of structural distortions upon ligand binding have been tested. In all cases, the conformational selection model has been confirmed and the conformational space between apo and holo forms has been encompassed. Different strategies have been tested that impact on the efficiency either to achieve a desired conformational change or to achieve a balanced exploration of the protein conformational multiplicity.

The Objectivity of Organizational Functions

We critique the organizational account of biological functions by showing how its basis in the closure of constraints fails to be objective. While the account treats constraints as objective features of physical systems, the number and relationship of potential constraints are subject to potentially arbitrary redescription by investigators. For example, we show that self-maintaining systems such as candle flames can realize closure on a more thorough analysis of the case, contradicting the claim that these "simple" systems lack functional organization. This also raises problems for Moreno and Mossio's associated theory of biological autonomy, which asserts that living beings are distinguished by their possession of a closed system of constraints that channel and regulate their metabolic processes.

Protein unfolding by SDS: the microscopic mechanisms and the properties of the SDS-protein assembly

The effects of detergent sodium dodecyl sulfate (SDS) on protein structure and dynamics are fundamental to the most common laboratory technique used to separate proteins and determine their molecular weights: polyacrylamide gel electrophoresis. However, the mechanism by which SDS induces protein unfolding and the microstructure of protein-SDS complexes remain largely unknown. Here, we report a detailed account of SDS-induced unfolding of two proteins-I27 domain of titin and β-amylase-obtained through all-atom molecular dynamics simulations. Both proteins were found to spontaneously unfold in the presence of SDS at boiling water temperature on the time scale of several microseconds. The protein unfolding was found to occur via two distinct mechanisms in which specific interactions of individual SDS molecules disrupt the protein's secondary structure. In the final state of the unfolding process, the proteins are found to wrap around SDS micelles in a fluid necklace-and-beads configuration, where the number and location of bound micelles changes dynamically. The global conformation of the protein was found to correlate with the number of SDS micelles bound to it, whereas the number of SDS molecules directly bound to the protein was found to define the relaxation time scale of the unfolded protein. Our microscopic characterization of SDS-protein interactions sets the stage for future refinement of SDS-enabled protein characterization methods, including protein fingerprinting and sequencing using a solid-state nanopore.

Mosaic Cooperativity in Slow Polypeptide Topological Isomerization Revealed by Residue-Specific NMR Thermodynamic Analysis

Slow polypeptide conformational changes on time scales of >1 s are generally assumed to be highly cooperative two-state transitions, reflecting the high energy barrier. However, few experimental characterizations have tested the validity of this assumption. We performed residue-specific NMR thermodynamic analysis of the 27-residue lantibiotic peptide, nukacin ISK-1, to characterize the isomerization between two topological states on the second time scale. Unexpectedly, the thermal transition behaviors were distinct among peptide regions, indicating that the topological isomerization process is a mosaic of different degrees of cooperativity. The conformational change path between the two NMR structures was deduced by a targeted molecular dynamics simulation. The unique side-chain threading motions through the monosulfide rings are the structural basis of the high energy barrier, and the nonlocal interactions in the hydrophobic core are the structural basis of the cooperativity. Taken together, we provide an energetic description of the topological isomerization of nukacin ISK-1.

Integrated Computational Approaches and Tools forAllosteric Drug Discovery

Understanding molecular mechanisms underlying the complexity of allosteric regulationin proteins has attracted considerable attention in drug discovery due to the benefits and versatilityof allosteric modulators in providing desirable selectivity against protein targets while minimizingtoxicity and other side effects. The proliferation of novel computational approaches for predictingligand-protein interactions and binding using dynamic and network-centric perspectives has ledto new insights into allosteric mechanisms and facilitated computer-based discovery of allostericdrugs. Although no absolute method of experimental and in silico allosteric drug/site discoveryexists, current methods are still being improved. As such, the critical analysis and integration ofestablished approaches into robust, reproducible, and customizable computational pipelines withexperimental feedback could make allosteric drug discovery more efficient and reliable. In this article,we review computational approaches for allosteric drug discovery and discuss how these tools can beutilized to develop consensus workflows for in silico identification of allosteric sites and modulatorswith some applications to pathogen resistance and precision medicine. The emerging realization thatallosteric modulators can exploit distinct regulatory mechanisms and can provide access to targetedmodulation of protein activities could open opportunities for probing biological processes and insilico design of drug combinations with improved therapeutic indices and a broad range of activities.

Protein tolerance to random circular permutation correlates with thermostability and local energetics of residue-residue contacts

Adenylate kinase (AK) orthologs with a range of thermostabilities were subjected to random circular permutation, and deep mutational scanning was used to evaluate where new protein termini were nondisruptive to activity. The fraction of circularly permuted variants that retained function in each library correlated with AK thermostability. In addition, analysis of the positional tolerance to new termini, which increase local conformational flexibility, showed that bonds were either functionally sensitive to cleavage across all homologs, differentially sensitive, or uniformly tolerant. The mobile AMP-binding domain, which displays the highest calculated contact energies, presented the greatest tolerance to new termini across all AKs. In contrast, retention of function in the lid and core domains was more dependent upon AK melting temperature. These results show that family permutation profiling identifies primary structure that has been selected by evolution for dynamics that are critical to activity within an enzyme family. These findings also illustrate how deep mutational scanning can be applied to protein homologs in parallel to differentiate how topology, stability, and local energetics govern mutational tolerance.

Large-Scale Conformational Changes and Protein Function: Breaking the in silico Barrier

Large-scale conformational changes are essential to link protein structures with their function at the cell and organism scale, but have been elusive both experimentally and computationally. Over the past few years developments in cryo-electron microscopy and crystallography techniques have started to reveal multiple snapshots of increasingly large and flexible systems, deemed impossible only short time ago. As structural information accumulates, theoretical methods become central to understand how different conformers interconvert to mediate biological function. Here we briefly survey current in silico methods to tackle large conformational changes, reviewing recent examples of cross-validation of experiments and computational predictions, which show how the integration of different scale simulations with biological information is already starting to break the barriers between the in silico, in vitro, and in vivo worlds, shedding new light onto complex biological problems inaccessible so far.

Network analysis of dynamically important residues in protein structures mediating ligand-binding conformational changes

According to the generalized conformational selection model, ligand binding involves the co-existence of at least two conformers with different ligand-affinities in a dynamical equilibrium. Conformational transitions between them should be guaranteed by intramolecular vibrational dynamics associated to each conformation. These motions are, therefore, related to the biological function of a protein. Positions whose mutations are found to alter these vibrations the most can be defined as key positions, that is, dynamically important residues that mediate the ligand-binding conformational change. In a previous study, we have shown that these positions are evolutionarily conserved. They correspond to buried aliphatic residues mostly localized in regular structured regions of the protein like β-sheets and α-helices. In the present paper, we perform a network analysis of these key positions for a large dataset of paired protein structures in the ligand-free and ligand-bound form. We observe that networks of interactions between these key positions present larger and more integrated networks with faster transmission of the information. Besides, networks of residues result that are robust to conformational changes. Our results reveal that the conformational diversity of proteins seems to be guaranteed by a network of strongly interconnected key positions rather than individual residues.

Atomistic Modeling of the ABL Kinase Regulation by Allosteric Modulators Using Structural Perturbation Analysis and Community-Based Network Reconstruction of Allosteric Communications

In this work, we have examined the molecular mechanisms of allosteric regulation of the ABL tyrosine kinase at the atomic level. Atomistic modeling of the ABL complexes with a panel of allosteric modulators has been performed using a combination of molecular dynamics simulations, structural residue perturbation scanning, and a novel community analysis of the residue interaction networks. Our results have indicated that allosteric inhibitors and activators may exert a differential control on allosteric signaling between the kinase binding sites and functional regions. While the inhibitor binding can strengthen the closed ABL state and induce allosteric communications directed from the allosteric pocket to the ATP binding site, the DPH activator may induce a more dynamic open form and activate allosteric couplings between the ATP and substrate binding sites. By leveraging a network-centric theoretical framework, we have introduced a novel community analysis method and global topological parameters that have unveiled the hierarchical modularity and the intercommunity bridging sites in the residue interaction network. We have found that allosteric functional hotspots responsible for the kinase regulation may serve the intermodular bridges in the global interaction network. The central conclusion from this analysis is that the regulatory switch centers play a fundamental role in the modular network organization of ABL as the unique intercommunity bridges that connect the SH2 and SH3 domains with the catalytic core into a functional kinase assembly. The hierarchy of network organization in the ABL regulatory complexes may allow for the synergistic action of dense intercommunity links required for the robust signal transfer in the catalytic core and sparse network bridges acting as the regulatory control points that orchestrate allosteric transitions between the inhibited and active kinase forms.

Advances in NMR Methods to Identify Allosteric Sites and Allosteric Ligands

NMR allows assessment of protein structure in solution. Unlike conventional X-ray crystallography that provides snapshots of protein conformations, all conformational states are simultaneously accessible to analysis by NMR. This is a significant advantage for discovery and characterization of allosteric effects. These effects are observed when binding at one site of the protein affects another distinct site through conformational transitions. Allosteric regulation of proteins has been observed in multiple physiological processes in health and disease, providing an opportunity for the development of allosteric inhibitors. These compounds do not directly interact with the orthosteric site of the protein but influence its structure and function. In this book chapter, we provide an overview on how NMR methods are utilized to identify allosteric sites and to discover novel inhibitors, highlighting examples from the field. We also describe how NMR has contributed to understanding of allosteric mechanisms and propose that it is likely to play an important role in clarification and further development of key concepts of allostery.

Why the Energy Landscape of Barnase Is Hierarchical

We have used NMR and computational methods to characterize the dynamics of the ribonuclease barnase over a wide range of timescales in free and inhibitor-bound states. Using temperature- and denaturant-dependent measurements of chemical shift, we show that barnase undergoes frequent and highly populated hinge bending. Using relaxation dispersion, we characterize a slower and less populated motion with a rate of 750 ± 200 s⁻¹, involving residues around the lip of the active site, which occurs in both free and bound states and therefore suggests conformational selection. Normal mode calculations characterize correlated hinge bending motions on a very rapid timescale. These three measurements are combined with previous measurements and molecular dynamics calculations on barnase to characterize its dynamic landscape on timescales from picoseconds to milliseconds and length scales from 0.1 to 2.5 nm. We show that barnase has two different large-scale fluctuations: one on a timescale of 10⁻⁹−10⁻⁶ s that has no free energy barrier and is a hinge bending that is determined by the architecture of the protein; and one on a timescale of milliseconds (i.e., 750 s⁻¹) that has a significant free energy barrier and starts from a partially hinge-bent conformation. These two motions can be described as hierarchical, in that the more highly populated faster motion provides a platform for the slower (less probable) motion. The implications are discussed. The use of temperature and denaturant is suggested as a simple and general way to characterize motions on the intermediate ns-μs timescale.

Dynamics based clustering of globin family members

A methodology to cluster proteins based on their dynamics’ similarity is presented. For each pair of proteins from a dataset, the structures are superimposed, and the Anisotropic Network Model modes of motions are calculated. The twelve slowest modes from each protein are matched using a local mode alignment algorithm based on the local sequence alignment algorithm of Smith–Waterman. The dynamical similarity distance matrix is calculated based on the top scoring matches of each pair and the proteins are clustered using a hierarchical clustering algorithm. The utility of this method is exemplified on a dataset of protein chains from the globin family and a dataset of tetrameric hemoglobins. The results demonstrate the effect of the quaternary structure of globin members on their intrinsic dynamics and show good ability to distinguish between different states of hemoglobin, revealing the dynamical relations between them.

Dynamical comparison between myoglobin and hemoglobin

Myoglobin and hemoglobin are globular hemeproteins, when the former is a monomer and the latter a heterotetramer. Despite the structural similarity of myoglobin to α and β subunits of hemoglobin, there is a functional difference between the two proteins, owing to the quaternary structure of hemoglobin. The effect of the quaternary structure of hemoglobin on the intrinsic dynamics of its subunits is explored by dynamical comparison of the two proteins. Anisotropic Network Model modes of motion were calculated for hemoglobin and myoglobin. Dynamical comparison between the proteins was performed using global and local Anisotropic Network Model mode alignment algorithms based on the algorithms of Smith-Waterman and Needleman-Wunsch for sequence comparison. The results indicate that the quaternary structure of hemoglobin substantially alters the intrinsic dynamics of its subunits, an effect that may contribute to the functional difference between the two proteins. Local dynamics similarity between the proteins is still observed at the major exit route of the ligand.

Direct examination of the relevance for folding, binding and electron transfer of a conserved protein folding intermediate

Near the minimum free energy basin of proteins where the native ensemble resides, partly unfolded conformations of slightly higher energy can be significantly populated under native conditions. It has been speculated that they play roles in molecular recognition and catalysis, but they might represent contemporary features of the evolutionary process without functional relevance. Obtaining conclusive evidence on these alternatives is difficult because it requires comparing the performance of a given protein when populating and when not populating one such intermediate, in otherwise identical conditions. Wild type apoflavodoxin populates under native conditions a partly unfolded conformation (10% of molecules) whose unstructured region includes the binding sites for the FMN cofactor and for redox partner proteins. We recently engineered a thermostable variant where the intermediate is no longer detectable. Using the wild type and variant, we assess the relevance of the intermediate comparing folding kinetics, cofactor binding kinetics, cofactor affinity, X-ray structure, intrinsic dynamics, redox potential of the apoflavodoxin-cofactor complex (Fld), its affinity for partner protein FNR, and electron transfer rate within the Fld/FNR physiological complex. Our data strongly suggest the intermediate state, conserved in long-chain apoflavodoxins, is not required for the correct assembly of flavodoxin nor does it contribute to shape its electron transfer properties. This analysis can be applied to evaluate other native basin intermediates.

Resolving the NFκB Heterodimer Binding Paradox: Strain and Frustration Guide the Binding of Dimeric Transcription Factors

Many eukaryotic transcription factors function after forming oligomers. The choice of protein partners is a nonrandom event that has distinct functional consequences for gene regulation. In the present work we examine three dimers of transcription factors in the NFκB family: p50p50, p50p65, and p65p65. The NFκB dimers bind to a myriad of genomic sites and switch the targeted genes on or off with precision. The p65p50 heterodimer of NFκB is the strongest DNA binder, and its unbinding is controlled kinetically by molecular stripping from the DNA induced by IκB. In contrast, the homodimeric forms of NFκB, p50p50 and p65p65, bind DNA with significantly less affinity, which places the DNA residence of the homodimers under thermodynamic rather than kinetic control. It seems paradoxical that the heterodimer should bind more strongly than either of the symmetric homodimers since DNA is a nearly symmetric target. Using a variety of energy landscape analysis tools, here we uncover the features in the molecular architecture of NFκB dimers that are responsible for these drastically different binding free energies. We show that frustration in the heterodimer interface gives the heterodimer greater conformational plasticity, allowing the heterodimer to better accommodate the DNA. We also show how the elastic energy and mechanical strain in NFκB dimers can be found by extracting the principal components of the fluctuations in Cartesian coordinates as well as fluctuations in the space of physical contacts, which are sampled via simulations with a predictive energy landscape Hamiltonian. These energetic contributions determine the specific detailed mechanisms of binding and stripping for both homo- and heterodimers.

Insights on the conformational dynamics of human frataxin through modifications of loop-1

Human frataxin (FXN) is a highly conserved mitochondrial protein involved in iron homeostasis and activation of the iron-sulfur cluster assembly. FXN deficiency causes the neurodegenerative disease Friedreich's Ataxia. Here, we investigated the effect of alterations in loop-1, a stretch presumably essential for FXN function, on the conformational stability and dynamics of the native state. We generated four loop-1 variants, carrying substitutions, insertions and deletions. All of them were stable and well-folded proteins. Fast local motions (ps-ns) and slower long-range conformational dynamics (μs-ms) were altered in some mutants as judged by NMR. Particularly, loop-1 modifications impact on the dynamics of a distant region that includes residues from the β-sheet, helix α1 and the C-terminal. Remarkably, all the mutants retain the ability to activate cysteine desulfurase, even when two of them exhibit a strong decrease in iron binding, revealing a differential sensitivity of these functional features to loop-1 perturbation. Consequently, we found that even for a small and relatively rigid protein, engineering a loop segment enables to alter conformational dynamics through a long-range effect, preserving the native-state structure and important aspects of function.

Binding mechanism and dynamic conformational change of C subunit of PKA with different pathways

The catalytic subunit of PKA (PKAc) exhibits three major conformational states (open, intermediate, and closed) during the biocatalysis process. Both ATP and substrate/inhibitor can effectively induce the conformational changes of PKAc from open to closed states. Aiming to explore the mechanism of this allosteric regulation, we developed a coarse-grained model and analyzed the dynamics of conformational changes of PKAc during binding by performing molecular dynamics simulations for apo PKAc, binary PKAc (PKAc with ATP, PKAc with PKI), and ternary PKAc (PKAc with ATP and PKI). Our results suggest a mixed binding mechanism of induced fit and conformational selection, with the induced fit dominant. The ligands can drive the movements of Gly-rich loop as well as some regions distal to the active site in PKAc and stabilize them at complex state. In addition, there are two parallel pathways (pathway with PKAc-ATP as an intermediate and pathway PKAc-PKI as an intermediate) during the transition from open to closed states. By molecular dynamics simulations and rate constant analyses, we find that the pathway through PKAc-ATP intermediate is the main binding route from open to closed state because of the fact that the bound PKI will hamper ATP from successful binding and significantly increase the barrier for the second binding subprocess. These findings will provide fundamental insights of the mechanisms of PKAc conformational change upon binding.

The Exact Nuclear Overhauser Enhancement: Recent Advances

Although often depicted as rigid structures, proteins are highly dynamic systems, whose motions are essential to their functions. Despite this, it is difficult to investigate protein dynamics due to the rapid timescale at which they sample their conformational space, leading most NMR-determined structures to represent only an averaged snapshot of the dynamic picture. While NMR relaxation measurements can help to determine local dynamics, it is difficult to detect translational or concerted motion, and only recently have significant advances been made to make it possible to acquire a more holistic representation of the dynamics and structural landscapes of proteins. Here, we briefly revisit our most recent progress in the theory and use of exact nuclear Overhauser enhancements (eNOEs) for the calculation of structural ensembles that describe their conformational space. New developments are primarily targeted at increasing the number and improving the quality of extracted eNOE distance restraints, such that the multi-state structure calculation can be applied to proteins of higher molecular weights. We then review the implications of the exact NOE to the protein dynamics and function of cyclophilin A and the WW domain of Pin1, and finally discuss our current research and future directions.

Molecular Dynamics-Markov State Model of Protein Ligand Binding and Allostery in CRIB-PDZ: Conformational Selection and Induced Fit

Conformational selection and induced fit are well-known contributors to ligand binding and allosteric effects in proteins. Molecular dynamics (MD) simulations now enable the theoretical study of protein-ligand binding in terms of ensembles of interconverting microstates and the population shifts characteristic of "dynamical allostery." Here we investigate protein-ligand binding and allostery based on a Markov state model (MSM) with states and rates obtained from all-atom MD simulations. As an exemplary case, we consider the single domain protein par-6 PDZ with and without ligand and allosteric effector. This is one of the smallest proteins in which allostery has been experimentally observed. In spite of the increased complexity intrinsic to a statistical ensemble perspective, we find that conformational selection and induced fit mechanisms can be readily identified in the analysis. In the nonallosteric pathway, MD-MSM shows that PDZ binds ligand via conformational selection. However, the allosteric pathway requires an activation step that involves a conformational change induced by the allosteric effector Cdc42. Once in the allosterically activated state, we find that ligand binding can proceed by conformational selection. Our MD-MSM model predicts that allostery in this and possibly other systems involves both induced fit and conformational selection, not just one or the other.

Dynamical differences of hemoglobin and the ionotropic glutamate receptor in different states revealed by a new dynamics alignment method

A new algorithm for comparison of protein dynamics is presented. Compared protein structures are superposed and their modes of motions are calculated using the anisotropic network model. The obtained modes are aligned using the dynamic programming algorithm of Needleman and Wunsch, commonly used for sequence alignment. Dynamical comparison of hemoglobin in the T and R2 states reveals that the dynamics of the allosteric effector 2,3-bisphosphoglycerate binding site is different in the two states. These differences can contribute to the selectivity of the effector to the T state. Similar comparison of the ionotropic glutamate receptor in the kainate+(R,R)-2b and ZK bound states reveals that the kainate+(R,R)-2b bound states slow modes describe upward motions of ligand binding domain and the transmembrane domain regions. Such motions may lead to the opening of the receptor. The upper lobes of the LBDs of the ZK bound state have a smaller interface with the amino terminal domains above them and have a better ability to move together. The present study exemplifies the use of dynamics comparison as a tool to study protein function. Proteins 2017; 85:1507-1517. © 2014 Wiley Periodicals, Inc.

How Intrinsic Dynamics Mediates the Allosteric Mechanism in the ABC Transporter Nucleotide Binding Domain Dimer

A protein's architecture facilitates specific motions-intrinsic dynamic modes-that are employed to effect function. Here we used molecular dynamics (MD) simulations to investigate the dynamics of the MJ0796 ABC transporter nucleotide-binding domain (NBD). ABC transporter NBDs form a rotationally symmetric dimer whereby two equivalent active sites are formed at their interface; in complex with a dimer of transmembrane domains they hydrolyze ATP to energize translocation of substrates across cellular membranes. Our data suggest the ABC NBD's ensemble of functional states can be understood predominately in terms of conformational changes between its major subdomains, occurring along two orthogonal dynamic modes. The data show that ligands and oligomeric interactions modulate the equilibrium conformation of the NBD with respect to these motions, suggesting that allostery is achieved by affecting the energetic profile along these two modes. The observed dynamics and allostery integrate consonantly and logically within a mechanistic framework for the ABC NBD dimer, which is supported by a large body of experimental and theoretical data, providing a higher resolution view of the enzyme's dynamic cycle. Our study shows how valuable mechanistic inferences can be derived from accessible short-time scale MD simulations of an enzyme's substructures.

Comparing allosteric transitions in the domains of calmodulin through coarse-grained simulations

Calmodulin (CaM) is a ubiquitous Ca(2+)-binding protein consisting of two structurally similar domains with distinct stabilities, binding affinities, and flexibilities. We present coarse grained simulations that suggest that the mechanism for the domain's allosteric transitions between the open and closed conformations depends on subtle differences in the folded state topology of the two domains. Throughout a wide temperature range, the simulated transition mechanism of the N-terminal domain (nCaM) follows a two-state transition mechanism while domain opening in the C-terminal domain (cCaM) involves unfolding and refolding of the tertiary structure. The appearance of the unfolded intermediate occurs at a higher temperature in nCaM than it does in cCaM consistent with nCaM's higher thermal stability. Under approximate physiological conditions, the simulated unfolded state population of cCaM accounts for 10% of the population with nearly all of the sampled transitions (approximately 95%) unfolding and refolding during the conformational change. Transient unfolding significantly slows the domain opening and closing rates of cCaM, which can potentially influence its Ca(2+)-binding mechanism.

Tunable allosteric library of caspase-3 identifies coupling between conserved water molecules and conformational selection

The native ensemble of caspases is described globally by a complex energy landscape where the binding of substrate selects for the active conformation, whereas targeting an allosteric site in the dimer interface selects an inactive conformation that contains disordered active-site loops. Mutations and posttranslational modifications stabilize high-energy inactive conformations, with mostly formed, but distorted, active sites. To examine the interconversion of active and inactive states in the ensemble, we used detection of related solvent positions to analyze 4,995 waters in 15 high-resolution (<2.0 Å) structures of wild-type caspase-3, resulting in 450 clusters with the most highly conserved set containing 145 water molecules. The data show that regions of the protein that contact the conserved waters also correspond to sites of posttranslational modifications, suggesting that the conserved waters are an integral part of allosteric mechanisms. To test this hypothesis, we created a library of 19 caspase-3 variants through saturation mutagenesis in a single position of the allosteric site of the dimer interface, and we show that the enzyme activity varies by more than four orders of magnitude. Altogether, our database consists of 37 high-resolution structures of caspase-3 variants, and we demonstrate that the decrease in activity correlates with a loss of conserved water molecules. The data show that the activity of caspase-3 can be fine-tuned through globally desolvating the active conformation within the native ensemble, providing a mechanism for cells to repartition the ensemble and thus fine-tune activity through conformational selection.

Prediction of allosteric sites and mediating interactions through bond-to-bond propensities

Allostery is a fundamental mechanism of biological regulation, in which binding of a molecule at a distant location affects the active site of a protein. Allosteric sites provide targets to fine-tune protein activity, yet we lack computational methodologies to predict them. Here we present an efficient graph-theoretical framework to reveal allosteric interactions (atoms and communication pathways strongly coupled to the active site) without a priori information of their location. Using an atomistic graph with energy-weighted covalent and weak bonds, we define a bond-to-bond propensity quantifying the non-local effect of instantaneous bond fluctuations propagating through the protein. Significant interactions are then identified using quantile regression. We exemplify our method with three biologically important proteins: caspase-1, CheY, and h-Ras, correctly predicting key allosteric interactions, whose significance is additionally confirmed against a reference set of 100 proteins. The almost-linear scaling of our method renders it suitable for high-throughput searches for candidate allosteric sites.

Coupling Protein Dynamics with Proton Transport in Human Carbonic Anhydrase II

The role of protein dynamics in enzyme catalysis is one of the most highly debated topics in enzymology. The main controversy centers around what may be defined as functionally significant conformational fluctuations and how, if at all, these fluctuations couple to enzyme catalyzed events. To shed light on this debate, the conformational dynamics along the transition path surmounting the highest free energy barrier have been herein investigated for the rate limiting proton transport event in human carbonic anhydrase (HCA) II. Special attention has been placed on whether the motion of an excess proton is correlated with fluctuations in the surrounding protein and solvent matrix, which may be rare on the picosecond and subpicosecond time scales of molecular motions. It is found that several active site residues, which do not directly participate in the proton transport event, have a significant impact on the dynamics of the excess proton. These secondary participants are shown to strongly influence the active site environment, resulting in the creation of water clusters that are conducive to fast, moderately slow, or slow proton transport events. The identification and characterization of these secondary participants illuminates the role of protein dynamics in the catalytic efficiency of HCA II.

Mapping transiently formed and sparsely populated conformations on a complex energy landscape

Determining the structures, kinetics, thermodynamics and mechanisms that underlie conformational exchange processes in proteins remains extremely difficult. Only in favourable cases is it possible to provide atomic-level descriptions of sparsely populated and transiently formed alternative conformations. Here we benchmark the ability of enhanced-sampling molecular dynamics simulations to determine the free energy landscape of the L99A cavity mutant of T4 lysozyme. We find that the simulations capture key properties previously measured by NMR relaxation dispersion methods including the structure of a minor conformation, the kinetics and thermodynamics of conformational exchange, and the effect of mutations. We discover a new tunnel that involves the transient exposure towards the solvent of an internal cavity, and show it to be relevant for ligand escape. Together, our results provide a comprehensive view of the structural landscape of a protein, and point forward to studies of conformational exchange in systems that are less characterized experimentally.

DOI:http://dx.doi.org/10.7554/eLife.17505.001

Proteins are the workhorses of cells, where they perform a wide range of roles. To do so, they adopt specific three-dimensional structures that enable them to interact with other molecules as necessary. Often a protein needs to be able to shift between different states with distinct structures as it goes about its job. To fully understand how a protein works, it is important to be able to characterize these different structures and how the protein changes between them.

Many of the experimental techniques used to study protein structure rely on isolating the individual structural forms of a protein. Since many structures only exist briefly, this can be very difficult. To complement experimental results, computer simulations allow researchers to model how atoms behave within a molecule. However, a number of factors limit how well these models represent what happens experimentally, such as the accuracy of the physical description used for the modeling.

Wang et al. set out to test and benchmark how well computer simulations could model changes in structure for a protein called T4 lysozyme, which has been studied extensively using experimental techniques. T4 lysozyme exists in two different states that have distinct structures. By comparing existing detailed experimental measurements with the results of their simulations, Wang et al. found that the simulations could capture key aspects of how T4 lysozyme changes its shape.

The simulations described the structure of the protein in both states and accurately determined the relative proportion of molecules that are found in each state. They could also determine how long it takes for a molecule to change its shape from one state to the other. The findings allowed Wang et al. to describe in fine detail – down to the level of individual atoms – how the protein changes its shape and how mutations in the protein affect its ability to do so. A key question for future studies is whether these insights can be extended to other proteins that are less well characterized experimentally than T4 lysozyme.

DOI:http://dx.doi.org/10.7554/eLife.17505.002

Direct Investigation of Slow Correlated Dynamics in Proteins via Dipolar Interactions

The synchronization of native state motions as they transition between microstates influences catalysis kinetics, mediates allosteric interactions, and reduces the conformational entropy of proteins. However, it has proven difficult to describe native microstates because they are usually minimally frustrated and may interconvert on the micro- to millisecond time scale. Direct observation of concerted equilibrium fluctuations would therefore be an important tool for describing protein native states. Here we propose a strategy that relates NMR cross-correlated relaxation (CCR) rates between dipolar interactions to residual dipolar couplings (RDCs) of individual consecutive H(N)-N and H(α)-C(α) bonds, which act as a proxy for the peptide planes and the side chains, respectively. Using Xplor-NIH ensemble structure calculations restrained with the RDC and CCR data, we observe collective motions on time scales slower than nanoseconds in the backbone for GB3. To directly access the correlations from CCR, we develop a structure-free data analysis. The resulting dynamic correlation map is consistent with the ensemble-restrained simulations and reveals a complex network. In general, we find that the bond motions are on average slightly correlated and that the local environment dominates many observations. Despite this, some patterns are typical over entire secondary structure elements. In the β-sheet, nearly all bonds are weakly correlated, and there is an approximately binary alternation in correlation intensity corresponding to the solvent exposure/shielding alternation of the side chains. For α-helices, there is also a weak correlation in the H(N)-N bonds. The degree of correlation involving H(α)-C(α) bonds is directly affected by side-chain fluctuations, whereas loops show complex and nonuniform behavior.