Folding funnels and binding mechanisms

Angewandte Chemie International Edition, 2011

Sample preparation was optimized for dielectric measurements. Problems with direct current effects originating from high salt concentration were alleviated by using a 3 mM sample of human ubiquitin (wild type) dissolved in deionized ultrapure water (SG Wasseraufbereitung und Regenerierstation GmbH). The salt content of the ubiquitin sample was minimized by excessive dialyzing against ultrapure water, exchanging the water five times every eight hours. A final NaCl concentration of 2.2 mM was estimated from the zero frequency DC. Subsequently the dialyzed sample was lyophilized and redissolved at 3 mM concentration in ultrapure water. The integrity of such samples was checked by recording 2D 15 N, 1 H-heteronuclear single quantum correlation (HSQC) spectra of a 15 N-labeled protein ) and by dynamic light scattering . The experiments were done on a DynaPro Titan instrument (Wyatt Technology) at 308 K. Data analysis was performed with the built-in software DYNAMICS from 10 successive measurements.

Journal of The Royal Society Interface, 2006

The protein structure prediction (PSP) problem is concerned with the prediction of the folded, native, tertiary structure of a protein given its sequence of amino acids. It is a challenging and computationally open problem, as proven by the numerous methodological attempts and the research effort applied to it in the last few years. The potential energy functions used in the literature to evaluate the conformation of a protein are based on the calculations of two different interaction energies: local (bond atoms) and non-local (non-bond atoms). In this paper, we show experimentally that those types of interactions are in conflict, and do so by using the potential energy function Chemistry at HARvard Macromolecular Mechanics. A multi-objective formulation of the PSP problem is introduced and its applicability studied. We use a multi-objective evolutionary algorithm as a search procedure for exploring the conformational space of the PSP problem.

European Biophysics Journal, 2011

We review the role conformational ensembles can play in the analysis of biomolecular dynamics, molecular recognition, and allostery. We introduce currently available methods for generating ensembles of biomolecules and illustrate their application with relevant examples from the literature. We show how, for binding, conformational ensembles provide a way of distinguishing the competing models of induced fit and conformational selection. For allostery we review the classic models and show how conformational ensembles can play a role in unravelling the intricate pathways of communication that enable allostery to occur. Finally, we discuss the limitations of conformational ensembles and highlight some potential applications for the future.

European Journal of Biochemistry, 2004

Two tetracycline repressor (TetR) sequence variants sharing 63% identical amino acids were investigated in terms of their recognition specificity for tetracycline and anhydrotetracycline. Thermodynamic complex stabilities determined by urea-dependent unfolding reveal that tetracycline stabilizes both variants to a similar extent but that anhydrotetracycline discriminates between them significantly. Isofunctional TetR hybrid proteins of these sequence variants were constructed and their denaturation profiles identified residues 57 and 61 as the complex stability determinant. Association kinetics reveal different recognition of these TetR variants by anhydrotetracycline, but the binding constants indicate similar stabilization. The identified residues connect to an internal water network, which suggests that the discrepancy in the observed thermodynamics may be caused by an entropy effect. Exchange of these interacting residues between the two TetR variants appears to influence the flexibility of this water organization, demonstrating the importance of buried, structural water molecules for ligand recognition and protein function. Therefore, this structural module seems to be a key requisite for the plasticity of the multiple ligand binding protein TetR.

Nature Chemical Biology, 2009

Molecular recognition is central to all biological processes. For the past fifty years, Koshland's 'induced fit' hypothesis has been the textbook explanation for molecular recognition events. However, recent experimental evidence supports an alternative mechanism. 'Conformational selection' postulates that all protein conformations pre-exist, and the ligand selects the most favored conformation. Following binding the ensemble undergoes a population shift, redistributing the conformational states. Both conformational selection and induced fit appear to play roles. Following binding by a primary conformational selection event, optimization of sidechain and backbone interactions is likely to proceed by an induced fit mechanism. Conformational selection has been observed for protein-ligand, protein-protein, protein-DNA, protein-RNA and RNA-ligand interactions. These data support a new molecular recognition paradigm for processes as diverse as signaling, catalysis, gene regulation, and protein aggregation in disease, which has the potential to significantly impact our views and strategies in drug design, biomolecular engineering and molecular evolution.

BMC Biology, 2012

Cell proliferation, differentiation, migration and adhesion are essential processes in development, morphing cells into critical anatomical structures. When cells approach each other, they may have a seemingly simple choice between adhesion and repulsion; however, the precise positioning of cells, as in the case of vascular patterning or in controlling axon growth in the assembly of topographic neural maps, is highly complex.

Journal of Biological Chemistry, 2010

The relationship between amyloid deposition and cellular toxicity is still controversial. In addition to fibril-forming oligomers, other soluble Aβ forms (amyloid β-derived diffusible ligands (ADDLs)) were also suggested to form and to present different morphologies and mechanisms of toxicity. One ADDL type, the "globulomer," apparently forms independently of the fibril aggregation pathway. Even though many studies argue that such soluble Aβ oligomers are off fibril formation pathways, they may nonetheless share some structural similarity with protofibrils. NMR data of globulomer intermediates, "preglobulomers," suggested parallel in-register C-terminal β-sheets, with different N-terminal conformations. Based on experimental data, we computationally investigate four classes of Aβ dodecamers: fibril, fibril oligomer, prefibril/preglobulomer cluster, and globulomer models. Our simulations of the solvent protection of double-layered fibril and globulomer models reproduce experimental observations. Using a single layer Aβ fibril oligomer β-sheet model, we found that the C-terminal β-sheet in the fibril oligomer is mostly curved, preventing it from quickly forming a fibril and leading to its breaking into shorter pieces. The simulations also indicate that β-sheets packed orthogonally could be the most stable species for Aβ dodecamers. The major difference between fibril-forming oligomers and ADDL-like oligomers (globulomers) could be the exposure of Met-35 patches. Although the Met-35 patches are necessarily exposed in fibril-forming oligomers to allow their maturation into fibrils, the Met-35 patches in the globulomer are covered by other residues in the orthogonally packed Aβ peptides. Our results call attention to the possible existence of certain "critical intermediates" that can lead to both seeds and other soluble ADDL-like oligomers.

Proteins: Structure, Function, and Genetics, 2001

Traditionally, molecular disorder has been viewed as local or global instability. Molecules or regions displaying disorder have been considered inherently unstructured. The term has been routinely applied to cases for which no atomic coordinates can be derived from crystallized molecules. Yet, even when it appears that the molecules are disordered, prevailing conformations exist, with population times higher than those of all alternate conformations. Disordered molecules are the outcome of rugged energy landscapes away from the native state around the bottom of the funnel. Ruggedness has a biological function, creating a distribution of structured conformers that bind via conformational selection, driving association and multimolecular complex formation, whether chainlinked in folding or unlinked in binding. We classify disordered molecules into two types. The first type possesses a hydrophobic core. Here, even if the native conformation is unstable, it still has a large enough population time, enabling its experimental detection. In the second type, no such hydrophobic core exists. Hence, the native conformations of molecules belonging to this category have shorter population times, hindering their experimental detection. Although there is a continuum of distribution of hydrophobic cores in proteins, an empirical, statistically based hydrophobicity function may be used as a guideline for distinguishing the two disordered molecule types. Furthermore, the two types relate to steps in the protein folding reaction. With respect to protein design, this leads us to propose that engineering-optimized specific electrostatic interactions to avoid electrostatic repulsion would reduce the type I disordered state, driving the molten globule (MG) 3 native (N) state. In contrast, for overcoming the type II disordered state, in addition to specific interactions, a stronger hydrophobic core is also indicated, leading to the denatured 3 MG 3 N state. Proteins 2001;44:418 -427.

Proceedings of the National Academy of Sciences, 2005

Protein-protein binding usually involves structural changes that may extend beyond the rearrangements on a local scale, and cannot be explained by a classical lock-and-key mechanism. Several models have been advanced to explain the flexible binding of proteins such as the induced fit mechanism where the ligand is postulated to induce a conformational change at the interaction site upon binding, or the preexisting equilibrium hypothesis that assumes that protein samples an ensemble of conformations at equilibrium conditions and that the ligand binds selectively to an active conformation. We explored the equilibrium motions of proteins that exhibit relatively large (nonlocal) conformational changes upon protein binding using the Gaussian network model and the anisotropic network model of protein dynamics. For four complexes, LIR-1͞HLA-A2, Actin͞DNase I, CDK2͞cyclin, and CDK6͞ p16 INK4a , the motions calculated for the monomer exhibiting the largest conformational change, in its unbound (free) form, correlate with the experimentally observed structural changes upon binding. This study emphasizes the preexisting equilibrium͞conformational selection as a mechanism for protein-protein interaction and lends support the concept that proteins, in their native conformation, are predisposed to undergo conformational fluctuations that are relevant to, or even required for, their biological functions.

Journal of Biomolecular Structure and Dynamics, 2009

Anisotropic network model is used to generate a set of distinct conformations for cylophilin A (CypA). The native structure is deformed to different extents along each of the lowestfrequency modes (first 7 modes) both in negative and positive directions. Each node of the elastic network represents either a single atom in the high-resolution model or a single residue in the low-resolution model. Realistic conformations with energies close to or lower than the crystal structure and with satisfactory internal geometry are recovered by energy minimization using implicit solvation model. These conformations are then used for ensemble docking to the ligand cyclosporin A for both a further test of accuracy of generated conformers and exploration of different binding modes. Higher number of correctly docked ligands are obtained for conformations with low deformation factors as a result of lower root mean square distances with respect to crystal structure. Yet, surprisingly, the lowest binding energy is obtained for one of the highly deformed conformations as a result of its special contact with arginine side chain oriented towards binding site. Considering the fact that the cyclic ligand's backbone and protein's side chains are held rigid during docking, the conformers generated by high-and low-resolution elastic network models are almost equally successful in providing the correct binding mode. The shape of the binding pocket that incorporates crucial interaction sites for hydrogen bond formation is found to be another important determining factor for the success of the dock. Also, the small backbone variations of a few Ångstroms in magnitude at the loop regions surrounding the binding pocket can cause amino acids' side chains to be displaced by magnitudes of up to 10 Å and therefore have a strong influence on the efficiency of the conformational search during docking.

Biopolymers, 2004

Bovine seminal ribonuclease (BS-RNase) is a unique member of the pancreatic-like ribonuclease superfamily. This enzyme exists as two conformational isomers with distinctive biological properties. The structure of the major isomer is characterized by the swapping of the N-terminal segment (MϫM BS-RNase). In this article, the crystal structures of the ligand-free MϫM BS-RNase and its complex with 2Ј-deoxycitidylyl(3Ј,5Ј)-2Ј-deoxyadenosine derived from isomorphous crystals have been refined. Interestingly, the comparison between this novel ligandfree form and the previously published sulfate-bound structure reveals significant differences. In particular, the ligand-free MϫM BS-RNase is closer to the structure of MϫM BS-RNase productive complexes than to the sulfate-bound form. These results reveal that MϫM BS-RNase presents a remarkable flexibility, despite the structural constraints of the interchain disulfide bridges and the Correspondence to: Lelio Mazzarella,

Journal of Biological Chemistry, 2002

Studies in many rhodopsin-like G-protein-coupled receptors are providing a general scheme of the structural processes underlying receptor activation. Microdomains in several receptors have been identified that appear to function as activation switches. However, evidence is emerging that these receptor proteins exist in multiple conformational states. To study the molecular control of this switching process, we investigated the function of a microdomain involving the conserved helix 7 tyrosine in the serotonin 5HT2C receptor. This tyrosine of the NPXXY motif was substituted for all naturally occurring amino acids. Three distinct constitutively active receptor phenotypes were found: moderate, high, and "locked-on" constitutive activity. In contrast to the activity of the other receptor mutants, the high basal signaling of the locked-on Y7.53N mutant was neither increased by agonists nor decreased by inverse agonists. The Y7.53F mutant was uncoupled. Computational modeling based on the rhodopsin crystal structure suggested that Y7.53 interacts with the conserved aromatic ring at position 7.60 in the recently identified helix 8 domain. This provided a basis for seeking revertant mutations to correct the defective function of the Y7.53F receptor. When the Y7.53F receptor was mutated at position 7.60, the wild-type phenotype was restored. These results suggest that Y7.53 and Y7.60 contribute to a common functional microdomain connecting helices 7 and 8 that influences the switching of the 5HT2C receptor among multiple active and inactive conformations.

PLoS ONE, 2010

The cholera toxin B pentamer (CtxB 5 ), which belongs to the AB 5 toxin family, is used as a model study for protein assembly. The effect of the pH on the reassembly of the toxin was investigated using immunochemical, electrophoretic and spectroscopic methods. Three pH-dependent steps were identified during the toxin reassembly: (i) acquisition of a fully assembly-competent fold by the CtxB monomer, (ii) association of CtxB monomer into oligomers, (iii) acquisition of the native fold by the CtxB pentamer. The results show that CtxB 5 and the related heat labile enterotoxin LTB 5 have distinct mechanisms of assembly despite sharing high sequence identity (84%) and almost identical atomic structures. The difference can be pinpointed to four histidines which are spread along the protein sequence and may act together. Thus, most of the toxin B amino acids appear negligible for the assembly, raising the possibility that assembly is driven by a small network of amino acids instead of involving all of them.

Proteins: Structure, Function, and Bioinformatics, 2012

PLoS Computational Biology, 2014

Many protein-protein interactions (PPIs) are compelling targets for drug discovery, and in a number of cases can be disrupted by small molecules. The main goal of this study is to examine the mechanism of binding site formation in the interface region of proteins that are PPI targets by comparing ligand-free and ligand-bound structures. To avoid any potential bias, we focus on ensembles of ligand-free protein conformations obtained by nuclear magnetic resonance (NMR) techniques and deposited in the Protein Data Bank, rather than on ensembles specifically generated for this study. The measures used for structure comparison are based on detecting binding hot spots, i.e., protein regions that are major contributors to the binding free energy. The main tool of the analysis is computational solvent mapping, which explores the surface of proteins by docking a large number of small ''probe'' molecules. Although we consider conformational ensembles obtained by NMR techniques, the analysis is independent of the method used for generating the structures. Finding the energetically most important regions, mapping can identify binding site residues using ligand-free models based on NMR data. In addition, the method selects conformations that are similar to some peptide-bound or ligand-bound structure in terms of the properties of the binding site. This agrees with the conformational selection model of molecular recognition, which assumes such pre-existing conformations. The analysis also shows the maximum level of similarity between unbound and bound states that is achieved without any influence from a ligand. Further shift toward the bound structure assumes protein-peptide or protein-ligand interactions, either selecting higher energy conformations that are not part of the NMR ensemble, or leading to induced fit. Thus, forming the sites in protein-protein interfaces that bind peptides and can be targeted by small ligands always includes conformational selection, although other recognition mechanisms may also be involved.

Physical Biology, 2013

Protein rigidity and flexibility can be analyzed accurately and efficiently using the program FIRST. Previous studies using FIRST were designed to analyze the rigidity and flexibility of proteins using a single static (snapshot) structure. It is however well known that proteins can undergo spontaneous sub-molecular unfolding and refolding, or conformational dynamics, even under conditions that strongly favour a well-defined native structure. These (local) unfolding events result in a large number of conformers that differ from each other very slightly. In this context, proteins are better represented as a thermodynamic ensemble of 'native-like' structures, and not just as a single static low-energy structure.

mAbs

We analyzed pairs of protein-binding, peptide-binding and hapten-binding antibodies crystallized as complex and in the absence of the antigen with and without conformational differences upon binding in the complementarity-determining region (CDR)-H3 loop. Here, we introduce a molecular dynamicsbased approach to capture a diverse conformational ensemble of the CDR-H3 loop in solution. The results clearly indicate that the inherently flexible CDR-H3 loop indeed needs to be characterized as a conformational ensemble. The conformational changes of the CDR-H3 loop in all antibodies investigated follow the paradigm of conformation selection, because we observe the experimentally determined binding competent conformation without the presence of the antigen within the ensemble of pre-existing conformational states in solution before binding. We also demonstrate for several examples that the conformation observed in the antibody crystal structure without antigen present is actually selected to bind the carboxyterminal tail region of the antigen-binding fragment (Fab). Thus, special care must be taken when characterizing antibody CDR-H3 loops by Fab X-ray structures, and the possibility that pre-existing conformations are present should always be considered.

Journal of Computational Biology, 2004

FlexProt is a novel technique for the alignment of exible proteins. Unlike all previous algorithms designed to solve the problem of structural comparisons allowing hinge-bending motions, FlexProt does not require an a priori knowledge of the location of the hinge(s). Flex-Prot carries out the exible alignment, superimposing the matching rigid subpart pairs, and detects the exible hinge regions simultaneously. A large number of methods are available to handle rigid structural alignment. However, proteins are exible molecules, which may appear in different conformations. Hence, protein structural analysis requires algorithms that can deal with molecular exibility. Here, we present a method addressing speci cally a exible protein alignment task. First, the method ef ciently detects maximal congruent rigid fragments in both molecules. Transforming the task into a graph theoretic problem, our method proceeds to calculate the optimal arrangement of previously detected maximal congruent rigid fragments. The fragment arrangement does not violate the protein sequence order. A clustering procedure is performed on fragment-pairs which have the same 3-D rigid transformation regardless of insertions and deletions (such as loops and turns) which separate them. Although the theoretical worst case complexity of the algorithm is O.n 6 /, in practice FlexProt is highly ef cient. It performs a structural comparison of a pair of proteins 300 amino acids long in about seven seconds on a standard desktop PC (400 MHz Pentium II processor with 256MB internal memory). We have performed extensive experiments with the algorithm. An assortment of these results is presented here. FlexProt can be accessed via WWW at bioinfo3d.cs.tau.ac.il/FlexProt/.

Physical Biology, 2011

A frequently neglected aspect of protein-protein interactions is flexibility. Small-scale fluctuations are present even in globular proteins, and alternative conformations can have a significant influence on the binding process. However, flexibility becomes highly prominent in complexes involving intrinsically disordered proteins. The importance of disordered regions in protein interactions has been recognized only relatively recently. In this survey we examine the basic properties of the complexes of disordered and ordered proteins from three different directions. The comparison of the interface properties shows that although disordered proteins can also adopt well-defined conformations in their bound form, their inherently dynamic nature is cast into their complexes. Furthermore, an overview of prediction methods indicates that disordered proteins as well as their binding regions can be recognized from the amino acid sequence by capturing the basic biophysical properties of these segments. Finally, we propose the generalization of the 'energy landscape model' for the description of complex formation that can help to put the various types of protein associations on a common ground.

BMC genomics, 2016

Structural excursions of a protein at equilibrium are key to biomolecular recognition and function modulation. Protein modeling research is driven by the need to aid wet laboratories in characterizing equilibrium protein dynamics. In principle, structural excursions of a protein can be directly observed via simulation of its dynamics, but the disparate temporal scales involved in such excursions make this approach computationally impractical. On the other hand, an informative representation of the structure space available to a protein at equilibrium can be obtained efficiently via stochastic optimization, but this approach does not directly yield information on equilibrium dynamics. We present here a novel methodology that first builds a multi-dimensional map of the energy landscape that underlies the structure space of a given protein and then queries the computed map for energetically-feasible excursions between structures of interest. An evolutionary algorithm builds such maps w...

Biomedicines, 2017

Targeted human cytolytic fusion proteins (hCFPs) are humanized immunotoxins for selective treatment of different diseases including cancer. They are composed of a ligand specifically binding to target cells genetically linked to a human apoptosis-inducing enzyme. hCFPs target cancer cells via an antibody or derivative (scFv) specifically binding to e.g., tumor associated antigens (TAAs). After internalization and translocation of the enzyme from endocytosed endosomes, the human enzymes introduced into the cytosol are efficiently inducing apoptosis. Under in vivo conditions such enzymes are subject to tight regulation by native inhibitors in order to prevent inappropriate induction of cell death in healthy cells. Tumor cells are known to up-regulate these inhibitors as a survival mechanism resulting in escape of malignant cells from elimination by immune effector cells. Cytosolic inhibitors of Granzyme B and Angiogenin (Serpin P9 and RNH1, respectively), reduce the efficacy of hCFPs with these enzymes as effector domains, requiring detrimentally high doses in order to saturate inhibitor binding and rescue cytolytic activity. Variants of Granzyme B and Angiogenin might feature reduced affinity for their respective inhibitors, while retaining or even enhancing their catalytic activity. A powerful tool to design hCFPs mutants with improved potency is given by in silico methods. These include molecular dynamics (MD) simulations and enhanced sampling methods (ESM). MD and ESM allow predicting the enzyme-protein inhibitor binding stability and the associated conformational changes, provided that structural information is available. Such "high-resolution" detailed description enables the elucidation of interaction domains and the identification of sites where particular point mutations may modify those interactions. This review discusses recent advances in the use of MD and ESM for hCFP development from the viewpoints of scientists involved in both fields.

Biophysical Reviews, 2013

Docking methodology aims to predict the experimental binding modes and affinities of small molecules within the binding site of particular receptor targets and is currently used as a standard computational tool in drug design for lead compound optimisation and in virtual screening studies to find novel biologically active molecules. The basic tools of a docking methodology include a search algorithm and an energy scoring function for generating and evaluating ligand poses. In this review, we present the search algorithms and scoring functions most commonly used in current molecular docking methods that focus on protein-ligand applications. We summarise the main topics and recent computational and methodological advances in protein-ligand docking. Protein flexibility, multiple ligand binding modes and the free-energy landscape profile for binding affinity prediction are important and interconnected challenges to be overcome by further methodological developments in the docking field.

PLoS Computational Biology, 2009

Conformational ensembles are increasingly recognized as a useful representation to describe fundamental relationships between protein structure, dynamics and function. Here we present an ensemble of ubiquitin in solution that is created by sampling conformational space without experimental information using ''Backrub'' motions inspired by alternative conformations observed in sub-Angstrom resolution crystal structures. Backrub-generated structures are then selected to produce an ensemble that optimizes agreement with nuclear magnetic resonance (NMR) Residual Dipolar Couplings (RDCs). Using this ensemble, we probe two proposed relationships between properties of protein ensembles: (i) a link between native-state dynamics and the conformational heterogeneity observed in crystal structures, and (ii) a relation between dynamics of an individual protein and the conformational variability explored by its natural family. We show that the Backrub motional mechanism can simultaneously explore protein native-state dynamics measured by RDCs, encompass the conformational variability present in ubiquitin complex structures and facilitate sampling of conformational and sequence variability matching those occurring in the ubiquitin protein family. Our results thus support an overall relation between protein dynamics and conformational changes enabling sequence changes in evolution. More practically, the presented method can be applied to improve protein design predictions by accounting for intrinsic native-state dynamics.

Scientific reports, 2017

Crucial for immune and anti-inflammatory cellular responses, signal transducer and activator of transcription 6 (STAT6) regulates transcriptional activation in response to interleukin-4 and -13 -induced tyrosine phosphorylation by direct interaction with coactivators. The interaction of STAT6 with nuclear coactivator 1 (NCoA1) is mediated by a short region of the STAT6 transactivation domain that includes the motif LXXLL and interacts with the PAS-B domain of NCoA1. Despite the availability of an X-ray structure of the PAS-B domain/ Leu794-Gly814-STAT6 complex, the mechanistic details of this interaction are still poorly understood. Here, we determine the structure of the NCoA1257-385/STAT6783-814 complex using Nuclear Magnetic Resonance (NMR) and X-ray crystallography. The STAT6783-814 peptide binds with additional N-terminal amino acids to NCoA1257-385, compared to the STAT6794-814 peptide, explaining its higher affinity. Secondary and tertiary structures existing in the free pept...

PloS one, 2011

The growing interest in quantifying the molecular basis of protein kinase activation and allosteric regulation by cancer mutations has fueled computational studies of allosteric signaling in protein kinases. In the present study, we combined computer simulations and the energy landscape analysis of protein kinases to characterize the interplay between oncogenic mutations and locally frustrated sites as important catalysts of allostetric kinase activation. While structurally rigid kinase core constitutes a minimally frustrated hub of the catalytic domain, locally frustrated residue clusters, whose interaction networks are not energetically optimized, are prone to dynamic modulation and could enable allosteric conformational transitions. The results of this study have shown that the energy landscape effect of oncogenic mutations may be allosteric eliciting global changes in the spatial distribution of highly frustrated residues. We have found that mutation-induced allosteric signaling...

The FEBS Journal

Tumor necrosis factor-a (TNF-a) is a pleiotropic inflammatory cytokine that exerts potent cytotoxic effects on solid tumor cells, while not affecting their normal counterparts. It is also known that TNF-a exerts many of its biological functions via interaction with specific receptors. To understand the potential roles of intrinsic disorder in the functioning of this important cytokine, we explored the peculiarities of intrinsic disorder distribution in human TNF-a and its homologs from various species, ranging from zebrafish to chimpanzee. We also studied the peculiarities of intrinsic disorder distribution in human TNF-a receptors, TNFR1 and TNFR2. Analysis revealed that cytoplasmic domains of TNF-a and its receptors are expected to be highly disordered. Furthermore, although the sequence identities of analyzed TNF-a homologs range from 99.57% (between human and chimpanzee proteins) to 22.33% (between frog and fish proteins), their intrinsic disorder profiles are characterized by a remarkable similarity. These observations indicate that the peculiarities of distribution of the intrinsic disorder propensity within the amino acid sequences are evolutionary conserved, and therefore could be of functional importance for this family of proteins. We also show that disordered and flexible regions of human TNF-a and its TNFR1 and TNFR2 receptors are crucial for some of their biological activities.

BMC genomics, 2018

The protein energy landscape underscores the inherent nature of proteins as dynamic molecules interconverting between structures with varying energies. Reconstructing a protein's energy landscape holds the key to characterizing a protein's equilibrium conformational dynamics and its relationship to function. Many pathogenic mutations in protein sequences alter the equilibrium dynamics that regulates molecular interactions and thus protein function. In principle, reconstructing energy landscapes of a protein's healthy and diseased variants is a central step to understanding how mutations impact dynamics, biological mechanisms, and function. Recent computational advances are yielding detailed, sample-based representations of protein energy landscapes. In this paper, we propose and describe two novel methods that leverage computed, sample-based representations of landscapes to reconstruct them and extract from them informative local structures that reveal the underlying org...

Protein & Cell

Arrestins are soluble relatively small 44-46 kDa proteins that specifically bind hundreds of active phosphorylated GPCRs and dozens of non-receptor partners. There are binding partners that demonstrate preference for each of the known arrestin conformations: free, receptor-bound, and microtubule-bound. Recent evidence suggests that conformational flexibility in every functional state is the defining characteristic of arrestins. Flexibility, or plasticity, of proteins is often described as structural disorder, in contrast to the fixed conformational order observed in high-resolution crystal structures. However, protein-protein interactions often involve highly flexible elements that can assume many distinct conformations upon binding to different partners. Existing evidence suggests that arrestins are no exception to this rule: their flexibility is necessary for functional versatility. The data on arrestins and many other multi-functional proteins indicate that in many cases, "order" might be artificially imposed by highly non-physiological crystallization conditions and/or crystal packing forces. In contrast, conformational flexibility (and its extreme case, intrinsic disorder) is a more natural state of proteins, representing true biological order that underlies their physiologically relevant functions.

History and Philosophy of the Life Sciences

Rediscoveries are not rare in biology. A recent example is the re-birth of the "fluctuation fit" concept developed by F. B. Straub and G. Szabolcsi in the sixties of the last century, under various names, the most popular of which is the "conformational selection". This theory offers an alternative to the "induced fit" concept by Koshland for the interpretation of the mechanism of protein—ligand interactions. A central question is whether the ligand induces a conformational change (as described by the induced fit model) or rather selects and stabilizes a complementary conformation from a pre-existing equilibrium of various states of the protein (according to the fluctuation fit/conformational selection model). Straub and Szabolcsi’s role and the factors hindering the spread of the fluctuation fit theory are discussed in the context of the history of the Hungarian biology in the 1950s and 1960s.

Molecular BioSystems, 2010

Intrinsic disorder is believed to contribute to the ability of some proteins to interact with multiple partners which is important for protein functional promiscuity and regulation of the cross-talk between pathways. To better understand the mechanisms of molecular recognition through disordered regions, here, we systematically investigate the coupling between disorder and binding within domain families in a structure interaction network and in terminal and inter-domain linker regions. We showed that the canonical domain-domain interaction model should take into account contributions of N-and C-termini and inter-domain linkers, which may form all or part of the binding interfaces. For the majority of proteins, binding interfaces on domain and terminal regions were predicted to be less disordered than non-interface regions. Analysis of all domain families revealed several exceptions, such as kinases, DNA/RNA binding proteins, certain enzymes, and regulatory proteins, which are candidates for disorder-to-order transitions that can occur upon binding. Domain interfaces that bind single or multiple partners do not exhibit significant difference in disorder content if normalized by the number of interactions. In general, protein families with more diverse interactions exhibit less average disorder over all members of the family. Our results shed light on recent controversies regarding the relationship between disorder and binding of multiple partners at common interfaces. In particular, they support the hypothesis that protein domains with many interacting partners should have a pleiotropic effect on functional pathways and consequently might be more constrained in evolution.

PLoS Computational Biology, 2009

We perform a large-scale study of intrinsically disordered regions in proteins and protein complexes using a non-redundant set of hundreds of different protein complexes. In accordance with the conventional view that folding and binding are coupled, in many of our cases the disorder-to-order transition occurs upon complex formation and can be localized to binding interfaces. Moreover, analysis of disorder in protein complexes depicts a significant fraction of intrinsically disordered regions, with up to one third of all residues being disordered. We find that the disorder in homodimers, especially in symmetrical homodimers, is significantly higher than in heterodimers and offer an explanation for this interesting phenomenon. We argue that the mechanisms of regulation of binding specificity through disordered regions in complexes can be as common as for unbound monomeric proteins. The fascinating diversity of roles of disordered regions in various biological processes and protein oligomeric forms shown in our study may be a subject of future endeavors in this area.

Biophysical Reviews

Autoinhibition is an effective mechanism that guards proteins against spurious activation. Despite its ubiquity, the distinct organizations of the autoinhibited states and their release mechanisms differ. Signaling is most responsive to the cell environment only if a small shift in the equilibrium is required to switch the system from an inactive (occluded) to an active (exposed) state. Ras signaling follows this paradigm. This underscores the challenge in pharmacological intervention to exploit and enhance autoinhibited states. Here, we review autoinhibition and release mechanisms at the membrane focusing on three representative Ras effectors, Raf protein kinase, PI3Kα lipid kinase, and NORE1A (RASSF5) tumor suppressor, and point to the ramifications to drug discovery. We further touch on Ras upstream and downstream signaling, Ras activation, and the Ras superfamily in this light, altogether providing a broad outlook of the principles and complexities of autoinhibition.

Biomolecules, 2020

WASp-interacting protein (WIP), a regulator of actin cytoskeleton assembly and remodeling, is a cellular multi-tasker and a key member of a network of protein–protein interactions, with significant impact on health and disease. Here, we attempt to complement the well-established understanding of WIP function from cell biology studies, summarized in several reviews, with a structural description of WIP interactions, highlighting works that present a molecular view of WIP’s protein–protein interactions. This provides a deeper understanding of the mechanisms by which WIP mediates its biological functions. The fully disordered WIP also serves as an intriguing example of how intrinsically disordered proteins (IDPs) exert their function. WIP consists of consecutive small functional domains and motifs that interact with a host of cellular partners, with a striking preponderance of proline-rich motif capable of interactions with several well-recognized binding partners; indeed, over 30% of ...

Proceedings of the National Academy of Sciences, 2020

There is ongoing debate regarding the mechanism through which cation/proton antiporters (CPAs), like Thermus thermophilus NapA (TtNapA) and Escherichia coli NapA (EcNhaA), alternate between their outward- and inward-facing conformations in the membrane. CPAs comprise two domains, and it is unclear whether the transition is driven by their rocking-bundle or elevator motion with respect to each other. Here we address this question using metadynamics simulations of TtNapA, where we bias conformational sampling along two axes characterizing the two proposed mechanisms: angular and translational motions, respectively. By applying the bias potential for the two axes simultaneously, as well as to the angular, but not the translational, axis alone, we manage to reproduce each of the two known states of TtNapA when starting from the opposite state, in support of the rocking-bundle mechanism as the driver of conformational change. Next, starting from the inward-facing conformation of EcNhaA, ...

Biomolecules, 2020

Robust biological systems are able to adapt to internal and environmental perturbations. This is ensured by a thick crosstalk between metabolism and signal transduction pathways, through which cell cycle progression, cell metabolism and growth are coordinated. Although several reports describe the control of cell signaling on metabolism (mainly through transcriptional regulation and post-translational modifications), much fewer information is available on the role of metabolism in the regulation of signal transduction. Protein-metabolite interactions (PMIs) result in the modification of the protein activity due to a conformational change associated with the binding of a small molecule. An increasing amount of evidences highlight the role of metabolites of the central metabolism in the control of the activity of key signaling proteins in different eukaryotic systems. Here we review the known PMIs between primary metabolites and proteins, through which metabolism affects signal transd...

The Journal of Chemical Physics, 2020

Molecular dynamics simulations are an invaluable tool to characterize the dynamic motions of proteins in atomistic detail. However, the accuracy of models derived from simulations inevitably relies on the quality of the underlying force field. Here, we present an evaluation of current non-polarizable and polarizable force fields (AMBER ff14SB, CHARMM 36m, GROMOS 54A7, and Drude 2013) based on the longstanding biophysical challenge of protein folding. We quantify the thermodynamics and kinetics of the β-hairpin formation using Markov state models of the fast-folding mini-protein CLN025. Furthermore, we study the (partial) folding dynamics of two more complex systems, a villin headpiece variant and a WW domain. Surprisingly, the polarizable force field in our set, Drude 2013, consistently leads to destabilization of the native state, regardless of the secondary structure element present. All non-polarizable force fields, on the other hand, stably characterize the native state ensembles in most cases even when starting from a partially unfolded conformation. Focusing on CLN025, we find that the conformational space captured with AMBER ff14SB and CHARMM 36m is comparable, but the ensembles from CHARMM 36m simulations are clearly shifted toward disordered conformations. While the AMBER ff14SB ensemble overstabilizes the native fold, CHARMM 36m and GROMOS 54A7 ensembles both agree remarkably well with experimental state populations. In addition, GROMOS 54A7 also reproduces experimental folding times most accurately. Our results further indicate an over-stabilization of helical structures with AMBER ff14SB. Nevertheless, the presented investigations strongly imply that reliable (un)folding dynamics of small proteins can be captured in feasible computational time with current additive force fields.

The FASEB Journal, 2019

The emergence of precision medicine has underscored the pressing need for reliable and comprehensive strategies for identification of mutations that drive cancer (eg, Refs. 1-13). The mutations can be at functional sites and directly abolish (or, strengthen) interactions. We dub these "orthosteric" driver mutations. Alternatively, they can be allosteric, in which case their mechanism is by destabilizing the inactive protein state; stabilizing the active state; or both. 14 These so-called "activating" (or, "deactivating," in repressors) mutations are allosteric drivers. 15 Since orthosteric drivers are at protein-protein or protein-ligand interaction sites, searches to detect them can be guided by these sites. 16-19 That however is not the case for allosteric mutations, making detection of rare allosteric driver mutations a challenging and daunting task. Rather than directly impede (or promote) interactions, allosteric drivers destabilize

The Journal of Physical Chemistry B, 2020

Allosteric molecules provide a powerful means to modulate protein function. However, the effect of such ligands on distal orthosteric sites cannot be easily described by classical docking methods. Here, we applied machine learning (ML) approaches to expose the links between local dynamic patterns and different degrees of allosteric inhibition of the ATPase function in the molecular chaperone TRAP1. We focused on 11 novel allosteric modulators with similar affinities to the target but with inhibitory efficacy between the 26.3 and 76%. Using a set of experimentally related local descriptors, ML enabled us to connect the molecular dynamics (MD) accessible to ligand-bound (perturbed) and unbound (unperturbed) systems to the degree of ATPase allosteric inhibition. The ML analysis of the comparative perturbed ensembles revealed a redistribution of dynamic states in the inhibitor-bound versus inhibitor-free systems following allosteric binding. Linear regression models were built to quantify the percentage of experimental variance explained by the predicted inhibitor-bound TRAP1 states. Our strategy provides a comparative MD−ML framework to infer allosteric ligand functionality. Alleviating the time scale issues which prevent the routine use of MD, a combination of MD and ML represents a promising strategy to support in silico mechanistic studies and drug design.

Journal of Computer-Aided Molecular Design, 2022

The discovery of new drugs is a time consuming and expensive process. Methods such as virtual screening, which can filter out ineffective compounds from drug libraries prior to expensive experimental study, have become popular research topics. As the computational drug discovery community has grown, in order to benchmark the various advances in methodology, organizations such as the Drug Design Data Resource have begun hosting blinded grand challenges seeking to identify the best methods for ligand pose-prediction, ligand affinity ranking, and free energy calculations. Such open challenges offer a unique opportunity for researchers to partner with junior students (e.g., high school and undergraduate) to validate basic yet fundamental hypotheses considered to be uninteresting to domain experts. Here, we, a group of high school-aged students and their mentors, present the results of our participation in Grand Challenge 4 where we predicted ligand affinity rankings for the Cathepsin S ...

Structure, 2021

Challenges in modelling in chemistry 1. Degrees of freedom governing chemical processes are electronic, nuclear, atomic, molecular and supra-molecular, and the corresponding particles have different masses and sizes 2. Interactions are governed by quantum mechanics: Dirac or Schrödinger equation, for not too small mass or too low temperature by classical equations of motion 3. At non-zero temperatures, the behavior of particles is governed by statistical mechanics: Fermi-Dirac, Bose-Einstein, or Boltzmann ensembles of configurations are to be considered, not single structures 4. The Coulomb interaction is rather long-ranged (distance-1) and induces many-body effects that make accurate modelling rather expensive 5. (Free) energy differences that drive processes can be very small compared to the total energy of the interacting particles 6. Time scales of processes easily span 15 orders of magnitude (fsec => sec) Challenge: striking an appropriate balance between accuracy and computational cost while maintaining a physically correct mechanism of the process of interest 11 W.F. van Gunsteren/New Orleans 080413/4

Nature Communications

Partner recognition in protein binding is critical for all biological functions, and yet, delineating its mechanism is challenging, especially when recognition happens within microseconds. We present a theoretical and experimental framework based on straight-forward nuclear magnetic resonance relaxation dispersion measurements to investigate protein binding mechanisms on sub-millisecond timescales, which are beyond the reach of standard rapid-mixing experiments. This framework predicts that conformational selection prevails on ubiquitin’s paradigmatic interaction with an SH3 (Src-homology 3) domain. By contrast, the SH3 domain recognizes ubiquitin in a two-state binding process. Subsequent molecular dynamics simulations and Markov state modeling reveal that the ubiquitin conformation selected for binding exhibits a characteristically extended C-terminus. Our framework is robust and expandable for implementation in other binding scenarios with the potential to show that conformationa...

Frontiers in Immunology, 2021

Antibodies have emerged as one of the fastest growing classes of biotherapeutic proteins. To improve the rational design of antibodies, we investigate the conformational diversity of 16 different germline combinations, which are composed of 4 different kappa light chains paired with 4 different heavy chains. In this study, we systematically show that different heavy and light chain pairings strongly influence the paratope, interdomain interaction patterns and the relative VH-VLinterface orientations. We observe changes in conformational diversity and substantial population shifts of the complementarity determining region (CDR) loops, resulting in distinct dominant solution structures and differently favored canonical structures. Additionally, we identify conformational changes in the structural diversity of the CDR-H3 loop upon different heavy and light chain pairings, as well as upon changes in sequence and structure of the neighboring CDR loops, despite having an identical CDR-H3 ...

Frontiers in Immunology, 2020

T-cell receptors are an important part in the adaptive immune system as they are responsible for detecting foreign proteins presented by the major histocompatibility complex (MHC). The affinity is predominantly determined by structure and sequence of the complementarity determining regions (CDRs), of which the CDR3 loops are responsible for peptide recognition. We present a kinetic classification of T-cell receptor CDR3 loops with different loop lengths into canonical and non-canonical solution structures. Using molecular dynamics simulations, we do not only sample available X-ray structures, but we also observe a substantially broader CDR3 loop ensemble with various distinct kinetic minima in solution. Our results strongly imply, that for given CDR3 loop sequences several canonical structures have to be considered to characterize the conformational diversity of these loops. Our suggested dominant solution structures could extend the repertoire of available canonical clusters by including kinetic minimum structures present in solution. Thus, the CDR3 loops need to be characterized as conformational ensembles in solution. Furthermore, the conformational changes of the CDR3 loops follow the paradigm of conformational selection, because the experimentally determined binding competent state is present within this ensemble of pre-existing conformations without the presence of the antigen. We also identify strong correlations between the CDR3 loops and include combined state descriptions. Additionally, we observe a strong dependency of the CDR3 loop conformations on the relative Vα-Vβ interdomain orientations, revealing that certain CDR3 loop states favor specific interface orientations.

Scientific Reports, 2020

We investigated T-cell receptor variable β chains binding to the superantigen staphylococcal enterotoxin C3 (SEC 3) with structure information in different stages of affinity maturation. Metadynamics in combination with molecular dynamics simulations allow to access the micro-to-millisecond timescale and reveal a strong effect of energetically significant mutations on the flexibility of the antigen-binding site. The observed changes in dynamics of the complementarity determining region (CDR) loops, especially the CDR 2, and HV 4 loop on this specific pathway of affinity maturation are reflected in their structural diversity, thermodynamics of conformations and kinetics of structural transitions. In addition, this affinity maturation pathway follows the concept of conformational selection, because even without the presence of the antigen the binding competent state is present in this pre-existing ensemble of conformations. In all stages of this affinity maturation process we observe ...

Frontiers in Genetics, 2021

The multi-level organization of nature is self-evident: proteins do interact among them to give rise to an organized metabolism, while in the same time each protein (a single node of such interaction network) is itself a network of interacting amino-acid residues allowing coordinated motion of the macromolecule and systemic effect as allosteric behavior. Similar pictures can be drawn for structure and function of cells, organs, tissues, and ecological systems. The majority of biologists are used to think that causally relevant events originate from the lower level (the molecular one) in the form of perturbations, that “climb up” the hierarchy reaching the ultimate layer of macroscopic behavior (e.g., causing a specific disease). Such causative model, stemming from the usual genotype-phenotype distinction, is not the only one. As a matter of fact, one can observe top-down, bottom-up, as well as middle-out perturbation/control trajectories. The recent complex network studies allow to ...

Antibiotics target key biological processes that include protein synthesis. Bacteria respond by developing resistance, which increases rapidly due to antibiotics overuse. Mupirocin, a clinically used natural antibiotic, inhibits isoleucyl-tRNA synthetase (IleRS), an enzyme that links isoleucine to its tRNAIlefor protein synthesis. Two IleRSs, mupirocin-sensitive IleRS1 and resistant IleRS2, coexist in bacteria. The latter may also be found in resistantStaphylococcus aureusclinical isolates. Here, we describe the structural basis of mupirocin resistance and unravel a mechanism of hyper-resistance evolved by some IleRS2 proteins. We surprisingly find that an up to 103-fold increase in resistance originates from alteration of the HIGH motif, a signature motif of the class I aminoacyl-tRNA synthetases to which IleRSs belong. The structural analysis demonstrates how an altered HIGH motif could be adopted in IleRS2 but not IleRS1, providing insight into an elegant mechanism for coevolutio...

PLoS Computational Biology, 2011

Molecular recognition is determined by the structure and dynamics of both a protein and its ligand, but it is difficult to directly assess the role of each of these players. In this study, we use Markov State Models (MSMs) built from atomistic simulations to elucidate the mechanism by which the Lysine-, Arginine-, Ornithine-binding (LAO) protein binds to its ligand. We show that our model can predict the bound state, binding free energy, and association rate with reasonable accuracy and then use the model to dissect the binding mechanism. In the past, this binding event has often been assumed to occur via an induced fit mechanism because the protein's binding site is completely closed in the bound state, making it impossible for the ligand to enter the binding site after the protein has adopted the closed conformation. More complex mechanisms have also been hypothesized, but these have remained controversial. Here, we are able to directly observe roles for both the conformational selection and induced fit mechanisms in LAO binding. First, the LAO protein tends to form a partially closed encounter complex via conformational selection (that is, the apo protein can sample this state), though the induced fit mechanism can also play a role here. Then, interactions with the ligand can induce a transition to the bound state. Based on these results, we propose that MSMs built from atomistic simulations may be a powerful way of dissecting ligand-binding mechanisms and may eventually facilitate a deeper understanding of allostery as well as the prediction of new protein-ligand interactions, an important step in drug discovery.

Chemical Science, 2021

Probing the energy landscape and thermodynamics of biomolecules for drug design.

PLoS Computational Biology, 2009

Structural and functional studies of the ABL and EGFR kinase domains have recently suggested a common mechanism of activation by cancer-causing mutations. However, dynamics and mechanistic aspects of kinase activation by cancer mutations that stimulate conformational transitions and thermodynamic stabilization of the constitutively active kinase form remain elusive. We present a large-scale computational investigation of activation mechanisms in the ABL and EGFR kinase domains by a panel of clinically important cancer mutants ABL-T315I, ABL-L387M, EGFR-T790M, and EGFR-L858R. We have also simulated the activating effect of the gatekeeper mutation on conformational dynamics and allosteric interactions in functional states of the ABL-SH2-SH3 regulatory complexes. A comprehensive analysis was conducted using a hierarchy of computational approaches that included homology modeling, molecular dynamics simulations, protein stability analysis, targeted molecular dynamics, and molecular docking. Collectively, the results of this study have revealed thermodynamic and mechanistic catalysts of kinase activation by major cancer-causing mutations in the ABL and EGFR kinase domains. By using multiple crystallographic states of ABL and EGFR, computer simulations have allowed one to map dynamics of conformational fluctuations and transitions in the normal (wild-type) and oncogenic kinase forms. A proposed multi-stage mechanistic model of activation involves a series of cooperative transitions between different conformational states, including assembly of the hydrophobic spine, the formation of the Src-like intermediate structure, and a cooperative breakage and formation of characteristic salt bridges, which signify transition to the active kinase form. We suggest that molecular mechanisms of activation by cancer mutations could mimic the activation process of the normal kinase, yet exploiting conserved structural catalysts to accelerate a conformational transition and the enhanced stabilization of the active kinase form. The results of this study reconcile current experimental data with insights from theoretical approaches, pointing to general mechanistic aspects of activating transitions in protein kinases.

Journal of Biological Chemistry, 2011

Lombricine kinase is a member of the phosphagen kinase family and a homolog of creatine and arginine kinases, enzymes responsible for buffering cellular ATP levels. Structures of lombricine kinase from the marine worm Urechis caupo were determined by x-ray crystallography. One form was crystallized as a nucleotide complex, and the other was substrate-free. The two structures are similar to each other and more similar to the substrate-free forms of homologs than to the substrate-bound forms of the other phosphagen kinases. Active site specificity loop 309-317, which is disordered in substrate-free structures of homologs and is known from the NMR of arginine kinase to be inherently dynamic, is resolved in both lombricine kinase structures, providing an improved basis for understanding the loop dynamics. Phosphagen kinases undergo a segmented closing on substrate binding, but the lombricine kinase ADP complex is in the open form more typical of substrate-free homologs. Through a comparison with prior complexes of intermediate structure, a correlation was revealed between the overall enzyme conformation and the substrate interactions of His 178. Comparative modeling provides a rationale for the more relaxed specificity of these kinases, of which the natural substrates are among the largest of the phosphagen substrates.