Intrinsic dynamics of enzymes in the unbound state and relation to allosteric regulation

Catch-and-hold activation of muscle acetylcholine receptors having transmitter binding site mutations

Agonists turn on receptors because their target sites have a higher affinity in the active versus resting conformation of the protein. We used single-channel electrophysiology to measure the lower-affinity (LA) and higher-affinity (HA) equilibrium dissociation constants for acetylcholine in adult-type muscle mouse nicotinic receptors (AChRs) having mutations of agonist binding site amino acids. For a series of agonists and for all mutations of αY93, αG147, αW149, αY190, αY198, εW55, and δW57, the change in LA binding energy was approximately half that in HA binding energy. The results were analyzed as a linear free energy relationship between LA and HA agonist binding, the slope of which (κ) gives the fraction of the overall binding chemical potential where the LA complex is established. The linear correlation between LA and HA binding energies suggests that the overall binding process is by an integrated mechanism (catch-and-hold). For the agonist and the above mutations, κ ∼ 0.5, but side-chain substitutions of two residues had a slope that was significantly higher (0.90; αG153) or lower (0.25; εP121). The results suggest that backbone rearrangements in loop B, loop C, and the non-α surface participate in both LA binding and the LA ↔ HA affinity switch. It appears that all of the intermediate steps in AChR activation comprise a single, energetically coupled process.

Antigen exposure leads to rigidification of germline antibody combining site

Immune complexes involving diverse antigens and corresponding antibodies were analyzed for mapping conformational transitions of an antibody before antigen binding, upon antigen binding and after antigen release. Molecular dynamics simulations of the two comprehensive datasets consisting of the antigen-free and antigen-bound structures of the germline antibodies 36-65 and BBE6.12H3 provided mechanistic model of antigen encounter by primary antibodies. While native germline antibodies exhibit substantial mobility in the antigen-combining sites, their antigen-bound states exhibit relatively rigid conformations, even in the absence of the antigen suggesting preservation of the structural state after antigen release. It is proposed that acquired rigidity by a germline antibody upon antigen binding may be the first step in affinity maturation in favor of that antigen.

ATPase subdomain IA is a mediator of interdomain allostery in Hsp70 molecular chaperones

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

Transient electrostatic interactions dominate the conformational equilibrium sampled by multidomain splicing factor U2AF65: a combined NMR and SAXS study

Multidomain proteins containing intrinsically disordered linkers exhibit large-scale dynamic modes that play key roles in a multitude of molecular recognition and signaling processes. Here, we determine the conformational space sampled by the multidomain splicing factor U2AF65 using complementary nuclear magnetic resonance spectroscopy and small-angle scattering data. Available degrees of conformational freedom are initially stochastically sampled and experimental data then used to delineate the potential energy landscape in terms of statistical probability. The spatial distribution of U2AF65 conformations is found to be highly anisotropic, comprising significantly populated interdomain contacts that appear to be electrostatic in origin. This hypothesis is supported by the reduction of signature PREs reporting on expected interfaces with increasing salt concentration. The described spatial distribution reveals the complete spectrum of the unbound forms of U2AF65 that coexist with the small percentage of a preformed RNA-bound domain arrangement required for polypyrimidine-tract recognition by conformational selection. More generally, the proposed approach to describing conformational equilibria of multidomain proteins can be further combined with other experimental data that are sensitive to domain dynamics.

Statistical model selection for Markov models of biomolecular dynamics

Markov state models provide a powerful framework for the analysis of biomolecular conformation dynamics in terms of their metastable states and transition rates. These models provide both a quantitative and comprehensible description of the long-time scale dynamics of large molecular dynamics with a Master equation and have been successfully used to study protein folding, protein conformational change, and protein-ligand binding. However, to achieve satisfactory performance, existing methodologies often require expert intervention when defining the model's discrete state space. While standard model selection methodologies focus on the minimization of systematic bias and disregard statistical error, we show that by consideration of the states' conditional distribution over conformations, both sources of error can be balanced evenhandedly. Application of techniques that consider both systematic bias and statistical error on two 100 μs molecular dynamics trajectories of the Fip35 WW domain shows agreement with existing techniques based on self-consistency of the model's relaxation time scales with more suitable results in regimes in which those time scale-based techniques encourage overfitting. By removing the need for expert tuning, these methods should reduce modeling bias and lower the barriers to entry in Markov state model construction.

A theoretical view of protein dynamics

Proteins are fascinating supramolecular structures, which are able to recognize ligands transforming binding information into chemical signals. They can transfer information across the cell, can catalyse complex chemical reactions, and are able to transform energy into work with much more efficiency than any human engine. The unique abilities of proteins are tightly coupled with their dynamic properties, which are coded in a complex way in the sequence and carefully refined by evolution. Despite its importance, our experimental knowledge of protein dynamics is still rather limited, and mostly derived from theoretical calculations. I will review here, in a systematic way, the current state-of-the-art theoretical approaches to the study of protein dynamics, emphasizing the most recent advances, examples of use and the expected lines of development in the near future.

An RNA aptamer possessing a novel monovalent cation-mediated fold inhibits lysozyme catalysis by inhibiting the binding of long natural substrates

RNA aptamers are being developed as inhibitors of macromolecular and cellular function, diagnostic tools, and potential therapeutics. Our understanding of the physical nature of this emerging class of nucleic acid-protein complexes is limited; few atomic resolution structures have been reported for aptamers bound to their protein target. Guided by chemical mapping, we systematically minimized an RNA aptamer (Lys1) selected against hen egg white lysozyme. The resultant 59-nucleotide compact aptamer (Lys1.2minE) retains nanomolar binding affinity and the ability to inhibit lysozyme's catalytic activity. Our 2.0-Å crystal structure of the aptamer-protein complex reveals a helical stem stabilizing two loops to form a protein binding platform that binds lysozyme distal to the catalytic cleft. This structure along with complementary solution analyses illuminate a novel protein-nucleic acid interface; (1) only 410 Å(2) of solvent accessible surface are buried by aptamer binding; (2) an unusually small fraction (∼18%) of the RNA-protein interaction is electrostatic, consistent with the limited protein phosphate backbone contacts observed in the structure; (3) a single Na(+) stabilizes the loops that constitute the protein-binding platform, and consistent with this observation, Lys1.2minE-lysozyme complex formation takes up rather than displaces cations at low ionic strength; (4) Lys1.2minE inhibits catalysis of large cell wall substrates but not catalysis of small model substrates; and (5) the helical stem of Lys1.2minE can be shortened to four base pairs (Lys1.2minF) without compromising binding affinity, yielding a 45-nucleotide aptamer whose structure may be an adaptable protein binding platform.

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)3. 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)3. 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.

Output-input ratio in thermally fluctuating biomolecular machines

Biological molecular machines are proteins that operate under isothermal conditions and hence are referred to as free energy transducers. They can be formally considered as enzymes that simultaneously catalyze two chemical reactions: the free energy-donating (input) reaction and the free energy-accepting (output) one. Most if not all biologically active proteins display a slow stochastic dynamics of transitions between a variety of conformational substates composing their native state. This makes the description of the enzymatic reaction kinetics in terms of conventional rate constants insufficient. In the steady state, upon taking advantage of the assumption that each reaction proceeds through a single pair (the gate) of transition conformational substates of the enzyme-substrates complex, the degree of coupling between the output and the input reaction fluxes has been expressed in terms of the mean first-passage times on a conformational transition network between the distinguished substates. The theory is confronted with the results of random-walk simulations on the five-dimensional hypercube. The formal proof is given that, for single input and output gates, the output-input degree of coupling cannot exceed unity. As some experiments suggest such exceeding, looking for the conditions for increasing the degree of coupling value over unity challenges the theory. Performed simulations of random walks on several model networks involving more extended gates indicate that the case of the degree of coupling value higher than 1 is realized in a natural way on critical branching trees extended by long-range shortcuts. Such networks are scale-free and display the property of the small world. For short-range shortcuts, the networks are scale-free and fractal, representing a reasonable model for biomolecular machines displaying tight coupling, i.e., the degree of coupling equal exactly to unity. A hypothesis is stated that the protein conformational transition networks, as just as higher-level biological networks, the protein interaction network, and the metabolic network, have evolved in the process of self-organized criticality.

Effect of intracellular loop 3 on intrinsic dynamics of human β2-adrenergic receptor

Background

To understand the effect of the long intracellular loop 3 (ICL3) on the intrinsic dynamics of human β₂-adrenergic receptor, molecular dynamics (MD) simulations were performed on two different models, both of which were based on the inactive crystal structure in complex with carazolol (after removal of carazolol and T4-lysozyme). In the so-called loop model, the ICL3 region that is missing in available crystal structures was modeled as an unstructured loop of 32-residues length, whereas in the clipped model, the two open ends were covalently bonded to each other. The latter model without ICL3 was taken as a reference, which has also been commonly used in recent computational studies. Each model was embedded into POPC bilayer membrane with explicit water and subjected to a 1 μs molecular dynamics (MD) simulation at 310 K.

Results

After around 600 ns, the loop model started a transition to a “very inactive” conformation, which is characterized by a further movement of the intracellular half of transmembrane helix 6 (TM6) towards the receptor core, and a close packing of ICL3 underneath the membrane completely blocking the G-protein’s binding site. Concurrently, the binding site at the extracellular part of the receptor expanded slightly with the Ser207-Asp113 distance increasing to 18 Å from 11 Å, which was further elaborated by docking studies.

Conclusions

The essential dynamics analysis indicated a strong coupling between the extracellular and intracellular parts of the intact receptor, implicating a functional relevance for allosteric regulation. In contrast, no such transition to the “very inactive” state, nor any structural correlation, was observed in the clipped model without ICL3. Furthermore, elastic network analysis using different conformers for the loop model indicated a consistent picture on the specific ICL3 conformational change being driven by global modes.

Specialized Dynamical Properties of Promiscuous Residues Revealed by Simulated Conformational Ensembles

The ability to interact with different partners is one of the most important features in proteins. Proteins that bind a large number of partners (hubs) have been often associated with intrinsic disorder. However, many examples exist of hubs with an ordered structure, and evidence of a general mechanism promoting promiscuity in ordered proteins is still elusive. An intriguing hypothesis is that promiscuous binding sites have specific dynamical properties, distinct from the rest of the interface and pre-existing in the protein isolated state. Here, we present the first comprehensive study of the intrinsic dynamics of promiscuous residues in a large protein data set. Different computational methods, from coarse-grained elastic models to geometry-based sampling methods and to full-atom Molecular Dynamics simulations, were used to generate conformational ensembles for the isolated proteins. The flexibility and dynamic correlations of interface residues with a different degree of binding promiscuity were calculated and compared considering side chain and backbone motions, the latter both on a local and on a global scale. The study revealed that (a) promiscuous residues tend to be more flexible than nonpromiscuous ones, (b) this additional flexibility has a higher degree of organization, and (c) evolutionary conservation and binding promiscuity have opposite effects on intrinsic dynamics. Findings on simulated ensembles were also validated on ensembles of experimental structures extracted from the Protein Data Bank (PDB). Additionally, the low occurrence of single nucleotide polymorphisms observed for promiscuous residues indicated a tendency to preserve binding diversity at these positions. A case study on two ubiquitin-like proteins exemplifies how binding promiscuity in evolutionary related proteins can be modulated by the fine-tuning of the interface dynamics. The interplay between promiscuity and flexibility highlighted here can inspire new directions in protein–protein interaction prediction and design methods.

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.

Large-scale analysis of the dynamics of enzymes

Protein enzymes enable the cell to execute chemical reactions in short time by accelerating the rate of the reactions in a selective manner. The motions or dynamics of the enzymes are essential for their function. Comparison of the dynamics of a set of 1247 nonhomologous enzymes was performed. For each enzyme, the slowest modes of motion are calculated using the Gaussian network model (GNM) and they are globally aligned. Alignment is done using the dynamic programming algorithm of Needleman and Wunsch, commonly used for sequence alignment. Only 96 pairs of proteins were identified to have three similar GNM slow modes with 63 of them having a similar structure. The most frequent slowest mode of motion describes a two domains anticorrelated motion that characterizes at least 23% of the enzymes. Therefore, dynamics uniqueness cannot be accounted for by the slowest mode itself but rather by the combination of several slow modes. Different quaternary structure packing can restrain the motion of enzyme subunits differently and may serve as another mechanism that increases the dynamics uniqueness.

Allostery and instability in the functional plasticity of synaptotagmin I

Synaptotagmin I (Syt I) is the calcium ion sensor for regulated release of neurotransmitter. How Syt I mediates this cellular event has been a question of extensive study for decades and yet, a clear understanding of the protein’s diverse functionality has remained elusive. Using tools of thermodynamics, we have identified two intrinsic properties that may account for Syt I’s functional plasticity: marginal stability and negative coupling. These two intrinsic properties have the potential to provide great conformational flexibility and suggest that Syt I’s functional plasticity stems in part from subtle rearrangements in the protein’s conformational ensemble. This model for Syt I function is discussed within the context of the nervous system’s overall plasticity.

A mechanistic understanding of allosteric immune escape pathways in the HIV-1 envelope glycoprotein

The HIV-1 envelope (Env) spike, which consists of a compact, heterodimeric trimer of the glycoproteins gp120 and gp41, is the target of neutralizing antibodies. However, the high mutation rate of HIV-1 and plasticity of Env facilitates viral evasion from neutralizing antibodies through various mechanisms. Mutations that are distant from the antibody binding site can lead to escape, probably by changing the conformation or dynamics of Env; however, these changes are difficult to identify and define mechanistically. Here we describe a network analysis-based approach to identify potential allosteric immune evasion mechanisms using three known HIV-1 Env gp120 protein structures from two different clades, B and C. First, correlation and principal component analyses of molecular dynamics (MD) simulations identified a high degree of long-distance coupled motions that exist between functionally distant regions within the intrinsic dynamics of the gp120 core, supporting the presence of long-distance communication in the protein. Then, by integrating MD simulations with network theory, we identified the optimal and suboptimal communication pathways and modules within the gp120 core. The results unveil both strain-dependent and -independent characteristics of the communication pathways in gp120. We show that within the context of three structurally homologous gp120 cores, the optimal pathway for communication is sequence sensitive, i.e. a suboptimal pathway in one strain becomes the optimal pathway in another strain. Yet the identification of conserved elements within these communication pathways, termed inter-modular hotspots, could present a new opportunity for immunogen design, as this could be an additional mechanism that HIV-1 uses to shield vulnerable antibody targets in Env that induce neutralizing antibody breadth.

A fast approximate method of identifying paths of allosteric communication in proteins

Fluctuations of the distance between a pair of residues i and j may be correlated with the fluctuations of the distance between another pair k and l. In this case, information may be transmitted among these four residues. Allosteric activity is postulated to proceed through such correlated paths. In this short communication a fast method for calculating correlations among all possible pairs ij and kl leading to a pathway of correlated residues of a protein is proposed. The method is based on the alpha carbon centered Gaussian Network Model. The model is applied to Glutamine Amidotransferase and pathways of allosteric activity are identified and compared with literature.

Epidemiology and molecular identification of Anisakis pegreffii (Nematoda: Anisakidae) in the horse mackerel Trachurus trachurus from northern Morocco

Anisakis infection parameters were studied in horse mackerel (Trachurus trachurus) taken from two areas of northern Morocco (Tetouan and Tangier), which showed a mean prevalence of 54.9%. Identification of the ITS1–5.8S–ITS2 fragment by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) showed A. pegreffii to be the dominant species; no A. simplex s.s. were detected. The presence of A. pegreffii in horse mackerel was not influenced by the sex (P= 0.46) or catch area (Atlantic versus Mediterranean, P= 0.52) of the fish, but was significantly related to their length, weight, liver weight and gonad weight, and to the season of their capture (P< 0.05). A low prevalence (8.6%) and mean intensity (1.0) was detected in the muscle, probably related to the reduced ability of A. pegreffii to penetrate muscle. The risk of the presence of A. pegreffii in the muscle was fivefold higher in fish caught during the summer than during any other season. Susceptible members of the human population can minimize the risk of infection by avoiding the consumption of larger horse mackerel specimens during the summer.

Structural and dynamic control of T-cell receptor specificity, cross-reactivity, and binding mechanism

Over the past two decades, structural biology has shown how T-cell receptors engage peptide/major histocompatibility complex (MHC) complexes and provided insight into the mechanisms underlying antigen specificity and cross-reactivity. Here we review and contextualize our contributions, which have emphasized the influence of structural changes and molecular flexibility. A repeated observation is the presence of conformational melding, in which the T-cell receptor (TCR), peptide, and in some cases, MHC protein cooperatively adjust in order for recognition to proceed. The structural changes reflect the intrinsic dynamics of the unligated proteins. Characterization of the dynamics of unligated TCR shows how binding loop motion can influence TCR cross-reactivity as well as specificity towards peptide and MHC. Examination of peptide dynamics indicates not only peptide-specific variation but also a peptide dependence to MHC flexibility. This latter point emphasizes that the TCR engages a composite peptide/MHC surface and that physically the receptor makes little distinction between the peptide and MHC. Much additional evidence for this can be found within the database of available structures, including our observations of a peptide dependence to the TCR binding mode and structural compensations for altered interatomic interactions, in which lost TCR-peptide interactions are replaced with TCR-MHC interactions. The lack of a hard-coded physical distinction between peptide and MHC has implications not only for specificity and cross-reactivity but also the mechanisms underlying MHC restriction as well as attempts to modulate and control TCR recognition.

Modeling cellular compartmentation in one-carbon metabolism

Folate-mediated one-carbon metabolism (FOCM) is associated with risk for numerous pathological states including birth defects, cancers, and chronic diseases. Although the enzymes that constitute the biological pathways have been well described and their interdependency through the shared use of folate cofactors appreciated, the biological mechanisms underlying disease etiologies remain elusive. The FOCM network is highly sensitive to nutritional status of several B-vitamins and numerous penetrant gene variants that alter network outputs, but current computational approaches do not fully capture the dynamics and stochastic noise of the system. Combining the stochastic approach with a rule-based representation will help model the intrinsic noise displayed by FOCM, address the limited flexibility of standard simulation methods for coarse-graining the FOCM-associated biochemical processes, and manage the combinatorial complexity emerging from reactions within FOCM that would otherwise be intractable.

Thermodynamics of allostery paves a way to allosteric drugs

We overview here our recent work on the thermodynamic view of allosteric regulation and communication. Starting from the geometry-based prediction of regulatory binding sites in a static structure, we move on to exploring a connection between ligand binding and the intrinsic dynamics of the protein molecule. We describe here two recently introduced measures, binding leverage and leverage coupling, which allow one to analyze the molecular basis of allosteric regulation. We discuss the advantages of these measures and show that they work universally in proteins of different sizes, oligomeric states, and functions. We also point the problems that have to be solved before completing an atomic level description of allostery, and briefly discuss ideas for computational design of allosteric drugs. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

Engineering allosteric control to an unregulated enzyme by transfer of a regulatory domain

Allosteric regulation of protein function is a critical component of metabolic control. Its importance is underpinned by the diversity of mechanisms and its presence in all three domains of life. The first enzyme of the aromatic amino acid biosynthesis, 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase, shows remarkable variation in allosteric response and machinery, and both contemporary regulated and unregulated orthologs have been described. To examine the molecular events by which allostery can evolve, we have generated a chimeric protein by joining the catalytic domain of an unregulated 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase with the regulatory domain of a regulated enzyme. We demonstrate that this simple gene fusion event on its own is sufficient to confer functional allostery to the unregulated enzyme. The fusion protein shares structural similarities with its regulated parent protein and undergoes an analogous major conformational change in response to the binding of allosteric effector tyrosine to the regulatory domain. These findings help delineate a remarkably facile mechanism for the evolution of modular allostery by domain recruitment.

Protein dynamics governed by interfaces of high polarity and low packing density

The folding pathway, three-dimensional structure and intrinsic dynamics of proteins are governed by their amino acid sequences. Internal protein surfaces with physicochemical properties appropriate to modulate conformational fluctuations could play important roles in folding and dynamics. We show here that proteins contain buried interfaces of high polarity and low packing density, coined as LIPs: Light Interfaces of high Polarity, whose physicochemical properties make them unstable. The structures of well-characterized equilibrium and kinetic folding intermediates indicate that the LIPs of the corresponding native proteins fold late and are involved in local unfolding events. Importantly, LIPs can be identified using very fast and uncomplicated computational analysis of protein three-dimensional structures, which provides an easy way to delineate the protein segments involved in dynamics. Since LIPs can be retained while the sequences of the interacting segments diverge significantly, proteins could in principle evolve new functional features reusing pre-existing encoded dynamics. Large-scale identification of LIPS may contribute to understanding evolutionary constraints of proteins and the way protein intrinsic dynamics are encoded.

Elastic network normal modes provide a basis for protein structure refinement

It is well recognized that thermal motions of atoms in the protein native state, the fluctuations about the minimum of the global free energy, are well reproduced by the simple elastic network models (ENMs) such as the anisotropic network model (ANM). Elastic network models represent protein dynamics as vibrations of a network of nodes (usually represented by positions of the heavy atoms or by the C(α) atoms only for coarse-grained representations) in which the spatially close nodes are connected by harmonic springs. These models provide a reliable representation of the fluctuational dynamics of proteins and RNA, and explain various conformational changes in protein structures including those important for ligand binding. In the present paper, we study the problem of protein structure refinement by analyzing thermal motions of proteins in non-native states. We represent the conformational space close to the native state by a set of decoys generated by the I-TASSER protein structure prediction server utilizing template-free modeling. The protein substates are selected by hierarchical structure clustering. The main finding is that thermal motions for some substates, overlap significantly with the deformations necessary to reach the native state. Additionally, more mobile residues yield higher overlaps with the required deformations than do the less mobile ones. These findings suggest that structural refinement of poorly resolved protein models can be significantly enhanced by reduction of the conformational space to the motions imposed by the dominant normal modes.

Dynamics of endoglucanase catalytic domains: implications towards thermostability

Thermostable endoglucanases play a crucial role in the production of biofuels to breakdown plant cellulose. Analyzing their structure-dynamics relationship can inform about the origins of their thermostability. Although tertiary structures of many endoglucanase proteins are available, the relationship between thermostability, structure, and dynamics is not explored fully. We have generated elastic network models for thermostable and mesostable endoglucanases with the (αβ)₈ fold in substrate bound and unbound states. The comparative analyses shed light on the relation between protein dynamics, thermostability, and substrate binding. We observed specific differences in the dynamic behavior of catalytic residues in slow modes: while both the nucleophile and the acid/base donor residues show positively correlated motions in the thermophile, their dynamics is uncoupled in the mesophile. Our proof-of-concept comparison study suggests that global dynamics can be harnessed to further our understanding of thermostability.

Conditional disorder in chaperone action

Protein disorder remains an intrinsically fuzzy concept. Its role in protein function is difficult to conceptualize and its experimental study is challenging. Although a wide variety of roles for protein disorder have been proposed, establishing that disorder is functionally important, particularly in vivo, is not a trivial task. Several molecular chaperones have now been identified as conditionally disordered proteins; fully folded and chaperone-inactive under non-stress conditions, they adopt a partially disordered conformation upon exposure to distinct stress conditions. This disorder appears to be vital for their ability to bind multiple aggregation-sensitive client proteins and to protect cells against the stressors. The study of these conditionally disordered chaperones should prove useful in understanding the functional role for protein disorder in molecular recognition.

Possible origin of life between mica sheets: does life imitate mica?

The mica hypothesis for the origin of life proposes that life originated between the sheets of muscovite mica. This paper elaborates on two ways that life resembles what might have originated between mica sheets. First, enzymes: The configurations and dynamics of enzymes, with their substrates, cofactors, and sometimes transition metal ions, often resemble mica sheets, with their open-and-shut motions, acting on small molecules between them, sometimes assisted by transition metal ions. Second, organisms: Mica world had the potential to be a community or ecosystem of prebiotic organisms in a way unlike other models for the origin of life.

Tunable, mixed-resolution modeling using library-based Monte Carlo and graphics processing units

Building on our recently introduced library-based Monte Carlo (LBMC) approach, we describe a flexible protocol for mixed coarse-grained (CG)/all-atom (AA) simulation of proteins and ligands. In the present implementation of LBMC, protein side chain configurations are pre-calculated and stored in libraries, while bonded interactions along the backbone are treated explicitly. Because the AA side chain coordinates are maintained at minimal run-time cost, arbitrary sites and interaction terms can be turned on to create mixed-resolution models. For example, an AA region of interest such as a binding site can be coupled to a CG model for the rest of the protein. We have additionally developed a hybrid implementation of the generalized Born/surface area (GBSA) implicit solvent model suitable for mixed-resolution models, which in turn was ported to a graphics processing unit (GPU) for faster calculation. The new software was applied to study two systems: (i) the behavior of spin labels on the B1 domain of protein G (GB1) and (ii) docking of randomly initialized estradiol configurations to the ligand binding domain of the estrogen receptor (ERα). The performance of the GPU version of the code was also benchmarked in a number of additional systems.

Light on the structural communication in Ras GTPases

The graph theory was combined with fluctuation dynamics to investigate the structural communication in four small G proteins, Arf1, H-Ras, RhoA, and Sec4. The topology of small GTPases is such that it requires the presence of the nucleotide to acquire a persistent structural network. The majority of communication paths involves the nucleotide and does not exist in the unbound forms. The latter are almost devoid of high-frequency paths. Thus, small Ras GTPases acquire the ability to transfer signals in the presence of nucleotide, suggesting that it modifies the intrinsic dynamics of the protein through the establishment of regions of hyperlinked nodes with high occurrence of correlated motions. The analysis of communication paths in the inactive (S(GDP)) and active (S(GTP)) states of the four G proteins strengthened the separation of the Ras-like domain into two dynamically distinct lobes, i.e. lobes 1 and 2, representing, respectively, the N-terminal and C-terminal halves of the domain. In the framework of this separation, interfunctional states and interfamily differences could be inferred. The structure network undergoes a reshaping depending on the bound nucleotide. Nucleotide-dependent divergences in structural communication reach the maximum in Arf1 and the minimum in RhoA. In Arf1, the nucleotide-dependent paths essentially express a communication between the G box 4 (G4) and distal portions of lobe 1. In the S(GDP) state, the G4 communicates with the N-term, while, in the S(GTP) state, the G4 communicates with the switch II. Clear differences could be also found between Arf1 and the other three G proteins. In Arf1, the nucleotide tends to communicate with distal portions of lobe 1, whereas in H-Ras, RhoA, and Sec4 it tends to communicate with a cluster of aromatic/hydrophobic amino acids in lobe 2. These differences may be linked, at least in part, to the divergent membrane anchoring modes that would involve the N-term for the Arf family and the C-term for the Rab/Ras/Rho families.

Coupled motions during dynamics reveal a tunnel toward the active site regulated by the N-terminal α-helix in an acylaminoacyl peptidase

Acylaminoacyl peptidase (AAP) subfamily belongs to the prolyl oligopeptidase (POP) family of serine-proteases. There is a great interest in the definition of molecular mechanisms related to the activity and substrate recognition of these complex multi-domain enzymes. The active site relies at the interface between the C-terminal catalytic domain and the β-propeller domain, whose N-terminal region acts as a bridge to the hydrolase domain. In AAP, the N-terminal extension is characterized by a structurally conserved α1-helix, which is known to affect thermal stability and thermal dependence of the catalytic activity. In the present contribution, results from hundreds nanosecond all-atom molecular dynamics simulations, along with analyses of the networks of cross-correlated motions of a member of the AAP subfamily are discussed. The MD investigation identifies a tunnel that from the surrounding of the N-terminal α1-helix bring to the catalytic site. This cavity seems to be regulated by conformational changes of the α1-helix itself during the dynamics. The evidence here provided can be a useful guide for a better understanding of the mechanistic aspects related to AAP activity, but also for drug design purposes.

Biomolecular dynamics: order-disorder transitions and energy landscapes

While the energy landscape theory of protein folding is now a widely accepted view for understanding how relatively weak molecular interactions lead to rapid and cooperative protein folding, such a framework must be extended to describe the large-scale functional motions observed in molecular machines. In this review, we discuss (1) the development of the energy landscape theory of biomolecular folding, (2) recent advances toward establishing a consistent understanding of folding and function and (3) emerging themes in the functional motions of enzymes, biomolecular motors and other biomolecular machines. Recent theoretical, computational and experimental lines of investigation have provided a very dynamic picture of biomolecular motion. In contrast to earlier ideas, where molecular machines were thought to function similarly to macroscopic machines, with rigid components that move along a few degrees of freedom in a deterministic fashion, biomolecular complexes are only marginally stable. Since the stabilizing contribution of each atomic interaction is on the order of the thermal fluctuations in solution, the rigid body description of molecular function must be revisited. An emerging theme is that functional motions encompass order-disorder transitions and structural flexibility provides significant contributions to the free energy. In this review, we describe the biological importance of order-disorder transitions and discuss the statistical-mechanical foundation of theoretical approaches that can characterize such transitions.

A rigid network of long-range contacts increases thermostability in a mutant endoglucanase

Thermodynamic stability of a protein at elevated temperatures is a key factor for thermostable enzymes to catalyze their specific reactions. Yet our understanding of biological determinants of thermostability is far from complete. Many different atomistic factors have been suggested as possible means for such proteins to preserve their activity at high temperatures. Among these factors are specific local interatomic interactions or enrichment of specific amino acid types. The case of glycosyl hydrolase family endoglucanase of Trichoderma reesei defies current hypotheses for thermostability because a single mutation far from the active site (A35 V) converts this mesostable protein into a thermostable protein without significant change in the protein structure. This substantial change in enzymatic activity cannot be explained on the basis of local intramolecular interactions alone. Here we present a more global view of the induced thermostability and show that the A35 V mutation affects the underlying structural rigidity of the whole protein via a number of long-range, non-local interactions. Our analysis of this structure reveals a precisely tuned, rigid network of atomic interactions. This cooperative, allosteric effect promotes the transformation of this mesostable protein into a thermostable one.

Coarse-grained simulations of transitions in the E2-to-E1 conformations for Ca ATPase (SERCA) show entropy-enthalpy compensation

SERCA is a membrane transport protein that has been extensively studied. There are a large number of highly resolved X-ray structures and several hundred mutations that have been characterized functionally. Despite this, the molecular details of the catalytic cycle, a cycle that includes large conformational changes, is not fully understood. In this computational study, we provide molecular dynamics descriptions of conformational changes during the E2→E1 transitions. The motivating point for these calculations was a series of insertion mutants in the A-M3 linker region that led to significant shifts in measured rates between the E2 and E1 states, as shown by experimental characterization. Using coarse-grained dynamic importance sampling within the context of a population shift framework, we sample on the intermediates along the transition pathway to address the mechanism for the conformational changes and the effects of the insertion mutations on the kinetics of the transition. The calculations define an approximation for the relative changes in entropy and enthalpy along the transition. These are found to be important for understanding the experimentally observed differences in rates. In particular, the interactions between cytoplasmic domains, water interactions, and the shifts in protein degrees of freedom with the insertion mutations show mutual compensation for the E2→E1 transitions in wild-type and mutant systems.

The role of domain: domain interactions versus domain: water interactions in the coarse-grained simulations of the E1P to E2P transitions in Ca-ATPase (SERCA)

SERCA is an important model system for understanding the molecular details of conformational change in membrane transport systems. This reflects the large number of solved X-ray structures and the equally large database of mutations that have been assayed. In this computational study, we provide a molecular dynamics description of the conformational changes during the E1P → E2P transitions. This set of states further changes with insertion mutants in the A-M3 linker region. These mutants were experimentally shown to lead to significant shifts in rates between the E1P → E2P states. Using the population shift framework and dynamic importance sampling method along with coarse-grained representations of the protein, lipid, and water, we suggest why these changes are found. The calculations sample on intermediates and suggest that changes in interactions, individual helix interactions, and water behavior are key elements in the molecular compositions that underlie shifts in kinetics. In particular, as the insertion length grows, it attracts more water and disrupts domain interactions, creating changes as well at the sites of key helix interactions between the A-Domain and the P-Domain. This provides a conceptual picture that aids understanding of the experimental results.

Simulations of temperature and salt concentration effects on bZIP, a basic region leucine zipper

Basic region leucine zipper (bZIP) transcription factors are dimeric proteins that recognize DNA. The monomers consist of a leucine zipper subdomain responsible for dimerization and a highly basic DNA recognition subdomain. Twelve explicit solvent molecular dynamics (MD) trajectories were run on the GCN4 bZIP transcriptional factor in the absence of DNA at three temperatures and two ion concentrations (0 mM with Cl(-) ions to neutralize the bZIP and 200 mM with additional Na(+) and Cl(-) ions) to probe the conformational ensemble that the basic region samples. In most trajectories, the basic region exhibits an alligator-jaw-like opening and closing (only one monomer moves), versus scissor-like motion, by a mainly rigid body, hinge motion centered on three "fork" residues that span the basic region to the coiled coil. In this motion, the α-helical character of the basic region monomers is mostly maintained. A broad range of distances is accessed, consistent with the absence of particular interactions for the basic region monomers. In two of the trajectories, the basic region monomers "collapse" to form a stable state. The coiled coil, leucine zipper subdomain is very stable for all of the trajectories. Ion solvation of the charged residue side chains is transient, on the scale of a few picoseconds. There is no evidence for persistent specific ion salt bridges to charged residues. For 0 mM, only certain basic region positively charged residues are substantially Cl(-) ion salt bridged. For 200 mM, in addition, some basic region negatively (positively) charged residues are salt bridged to Na(+) (Cl(-)) ions. The different ion solvation patterns at the two ion concentrations are not greatly temperature sensitive, and the conformational sampling found in the MD is remarkably unperturbed by ion concentration and/or temperature.

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.

Insights into the fold organization of TIM barrel from interaction energy based structure networks

There are many well-known examples of proteins with low sequence similarity, adopting the same structural fold. This aspect of sequence-structure relationship has been extensively studied both experimentally and theoretically, however with limited success. Most of the studies consider remote homology or “sequence conservation” as the basis for their understanding. Recently “interaction energy” based network formalism (Protein Energy Networks (PENs)) was developed to understand the determinants of protein structures. In this paper we have used these PENs to investigate the common non-covalent interactions and their collective features which stabilize the TIM barrel fold. We have also developed a method of aligning PENs in order to understand the spatial conservation of interactions in the fold. We have identified key common interactions responsible for the conservation of the TIM fold, despite high sequence dissimilarity. For instance, the central beta barrel of the TIM fold is stabilized by long-range high energy electrostatic interactions and low-energy contiguous vdW interactions in certain families. The other interfaces like the helix-sheet or the helix-helix seem to be devoid of any high energy conserved interactions. Conserved interactions in the loop regions around the catalytic site of the TIM fold have also been identified, pointing out their significance in both structural and functional evolution. Based on these investigations, we have developed a novel network based phylogenetic analysis for remote homologues, which can perform better than sequence based phylogeny. Such an analysis is more meaningful from both structural and functional evolutionary perspective. We believe that the information obtained through the “interaction conservation” viewpoint and the subsequently developed method of structure network alignment, can shed new light in the fields of fold organization and de novo computational protein design.

Proteins are polymers of amino-acids that fold into unique three-dimensional structures to perform cellular functions. This structure formation has been shown to depend on the amino-acid sequences. But examples of proteins with diverse sequences retaining a similar structural fold are quite substantial that we can no longer consider such phenomenon as exceptions. Therefore, this non-canonical relationship has been studied extensively mostly by studying the remote sequence similarities between proteins. Here we have attempted to address the above-mentioned problem by analyzing the similarities in the spatial interactions among amino-acids. Since the protein structure is a resultant of different interactions, we have considered the proteins as networks of interacting amino-acids to derive the common interactions within a popular structural fold called the TIM barrel fold. We were able to find common interactions among different families of the TIM fold and generalize the patterns of interactions by which the fold is being maintained despite sequence diversity. The results substantiate our hypothesis that interaction conservation might by a driving factor in fold formation and this new outlook can be used extensively in engineering proteins with better biophysical characteristics.

Novel inhibitor discovery through virtual screening against multiple protein conformations generated via ligand-directed modeling: a maternal embryonic leucine zipper kinase example

Kinase targets have been demonstrated to undergo major conformational reorganization upon ligand binding. Such protein conformational plasticity remains a significant challenge in structure-based virtual screening methodology and may be approximated by screening against an ensemble of diverse protein conformations. Maternal embryonic leucine zipper kinase (MELK), a member of serine-threonine kinase family, has been recently found to be involved in the tumerogenic state of glioblastoma, breast, ovarian, and colon cancers. We therefore modeled several conformers of MELK utilizing the available chemogenomic and crystallographic data of homologous kinases. We carried out docking pose prediction and virtual screening enrichment studies with these conformers. The performances of the ensembles were evaluated by their ability to reproduce known inhibitor bioactive conformations and to efficiently recover known active compounds early in the virtual screen when seeded with decoy sets. A few of the individual MELK conformers performed satisfactorily in reproducing the native protein-ligand pharmacophoric interactions up to 50% of the cases. By selecting an ensemble of a few representative conformational states, most of the known inhibitor binding poses could be rationalized. For example, a four conformer ensemble is able to recover 95% of the studied actives, especially with imperfect scoring function(s). The virtual screening enrichment varied considerably among different MELK conformers. Enrichment appears to improve by selection of a proper protein conformation. For example, several holo and unliganded active conformations are better to accommodate diverse chemotypes than ATP-bound conformer. These results prove that using an ensemble of diverse conformations could give a better performance. Applying this approach, we were able to screen a commercially available library of half a million compounds against three conformers to discover three novel inhibitors of MELK, one from each template. Among the three compounds validated via experimental enzyme inhibition assays, one is relatively potent (15; K(d) = 0.37 μM), one moderately active (12; K(d) = 3.2 μM), and one weak but very selective (9; K(d) = 18 μM). These novel hits may be utilized to assist in the development of small molecule therapeutic agents useful in diseases caused by deregulated MELK, and perhaps more importantly, the approach demonstrates the advantages of choosing an appropriate ensemble of a few conformers in pursuing compound potency, selectivity, and novel chemotypes over using single target conformation for structure-based drug design in general.

Structural insights into the autoactivation mechanism of p21-activated protein kinase

p21-activated kinases (PAKs) play an important role in diverse cellular processes. Full activation of PAKs requires autophosphorylation of a critical threonine/serine located in the activation loop of the kinase domain. Here we report crystal structures of the phosphorylated and unphosphorylated PAK1 kinase domain. The phosphorylated PAK1 kinase domain has a conformation typical of all active protein kinases. Interestingly, the structure of the unphosphorylated PAK1 kinase domain reveals an unusual dimeric arrangement expected in an authentic enzyme-substrate complex, in which the activation loop of the putative "substrate" is projected into the active site of the "enzyme." The enzyme is bound to AMP-PNP and has an active conformation, whereas the substrate is empty and adopts an inactive conformation. Thus, the structure of the asymmetric homodimer mimics a trans-autophosphorylation complex, and suggests that unphosphorylated PAK1 could dynamically adopt both the active and inactive conformations in solution.

Dynamics, flexibility and ligand-induced conformational changes in biological macromolecules: a computational approach

Biomolecules possess important dynamical properties that enable them to adapt and alternate their conformation as a response to environmental stimuli. Recent advancements in computational resources and methodology allow a higher capability to mimic in vitro conditions and open up the possibility of studying large systems over longer timescales. Here, we describe commonly used computational approaches for studying the dynamic properties of proteins. We review a selected set of simulation studies on ligand-induced changes in the chaperonin GroEL-GroES, a molecular folding machine, maltose-binding protein, a prototypical member of the periplasmic binding proteins, and the bacterial ribosomal A-site, focusing on aminoglycoside antibiotic recognition. We also discuss a recent quantitative reconstruction of the binding process of benzamidine and trypsin. These studies contribute to the understanding and further development of the medicinal regulation of large biomolecular systems.