Search

Myo-inositol-1-phosphate synthase (MIPS) catalyzes the NAD+-dependent isomerization of glucose-6-phosphate (G6P) into inositol-1-phosphate (IMP), controlling the rate-limiting step of the inositol pathway. Previous structural studies focused on the detailed molecular mechanism, neglecting large-scale conformational changes that drive the function of this 240 kDa homotetrameric complex. In this study, we identified the active, endogenous MIPS in cell extracts from the thermophilic fungus Thermochaetoides thermophila. By resolving the native structure at 2.48 Å (FSC = 0.143), we revealed a fully populated active site. Utilizing 3D variability analysis, we uncovered conformational states of MIPS, enabling us to directly visualize an order-to-disorder transition at its catalytic center. An acyclic intermediate of G6P occupied the active site in two out of the three conformational states, indicating a catalytic mechanism where electrostatic stabilization of high-energy intermediates plays a crucial role. Examination of all isomerases with known structures revealed similar fluctuations in secondary structure within their active sites. Based on these findings, we established a conformational selection model that governs substrate binding and eventually inositol availability. In particular, the ground state of MIPS demonstrates structural configurations regardless of substrate binding, a pattern observed across various isomerases. These findings contribute to the understanding of MIPS structure-based function, serving as a template for future studies targeting regulation and potential therapeutic applications. [ABSTRACT FROM AUTHOR]

Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom twodimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C+, U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin. [ABSTRACT FROM AUTHOR]

Protein homodimers have been classified as three-state or two-state dimers depending on whether a folded monomer forms before association, but the details of the folding–binding mechanisms are poorly understood. Kinetic transition networks of conformational states have provided insight into the folding mechanisms of monomeric proteins, but extending such a network to two protein chains is challenging as all the relative positions and orientations of the chains need to be included, greatly increasing the number of degrees of freedom. Here, we present a simplification of the problem by grouping all states of the two chains into two layers: a dissociated and an associated layer. We combined our two-layer approach with the Wako–Saito–Muñoz–Eaton method and used Transition Path Theory to investigate the dimer formation kinetics of eight homodimers. The analysis reveals a remarkable diversity of dimer formation mechanisms. Induced folding, conformational selection, and rigid docking are often simultaneously at work, and their contribution depends on the protein concentration. Pre-folded structural elements are always present at the moment of association, and asymmetric binding mechanisms are common. Our two-layer network approach can be combined with various methods that generate discrete states, yielding new insights into the kinetics and pathways of flexible binding processes. [ABSTRACT FROM AUTHOR]

The β2 adrenergic receptor (β2AR) is an archetypal G protein coupled receptor (GPCR). One structural signature of GPCR activation is a large-scale movement (ca. 6 to 14 Å) of transmembrane helix 6 (TM6) to a conformation which binds and activates a cognate G protein. The β2AR exhibits a low level of agonist-independent G protein activation. The structural origin of this basal activity and its suppression by inverse agonists is unknown but could involve a unique receptor conformation that promotes G protein activation. Alternatively, a conformational selection model proposes that a minor population of the canonical active receptor conformation exists in equilibrium with inactive forms, thus giving rise to basal activity of the ligand-free receptor. Previous spin-labeling and fluorescence resonance energy transfer experiments designed to monitor the positional distribution of TM6 did not detect the presence of the active conformation of ligand-free β2AR. Here we employ spin-labeling and pressure-resolved double electron-electron resonance spectroscopy to reveal the presence of a minor population of unliganded receptor, with the signature outward TM6 displacement, in equilibrium with inactive conformations. Binding of inverse agonists suppresses this population. These results provide direct structural evidence in favor of a conformational selection model for basal activity in β2AR and provide a mechanism for inverse agonism. In addition, they emphasize 1) the importance of minor populations in GPCR catalytic function; 2) the use of spin-labeling and variable-pressure electron paramagnetic resonance to reveal them in a membrane protein; and 3) the quantitative evaluation of their thermodynamic properties relative to the inactive forms, including free energy, partial molar volume, and compressibility.

Macugen is a therapeutic RNA aptamer against vascular endothelial growth factor (VEGF)‐165, the VEGF isoform primarily responsible for angiogenesis. It has been reported that Macugen inhibits angiogenesis by specifically binding to the heparin binding domain (HBD) of VEGF165. The mechanism of the molecular recognition between HBD and Macugen is investigated here using all‐atom molecular dynamics simulations. We find that Macugen recognizes HBD by an induced‐fit mechanism with major conformational changes in Macugen and almost no changes in the structure of HBD, whereas HBD recognizes Macugen by a conformational selection mechanism. [ABSTRACT FROM AUTHOR]

Functional oligopeptides derived from natural proteins are often used as targeting groups or therapeutic drugs but always suffer from a dramatic bioaffinity decrease due to the loss of conformational restriction in proteins, which increases their entropy state. It is worth noting that the linear Cys-Arg-Gly-Asp-Ser (CRGDS) oligopeptide can spontaneously restore its natural conformation by anchoring one terminus and forming hydrogen bonds on another terminus, on protein-like zwitterionic fifth-generation polyamide-amine (G5 PAMAM) nanocarriers, thus exhibiting high bioaffinity. However, the different physicochemical properties of nanocarriers have an impact on the bioaffinity recovery of linear CRGDS. In this work, it is found that the bioaffinity decreases with the reduction of generation of PAMAM and zwitterionic G3 PAMAM is a turning point for bioaffinity recovery via a series of in vitro experiments. This indicates that the spontaneous restoration of the native conformation of the RGD segment requires proper surface group density and structural rigidity of the zwitterionic nanocarrier for reducing the entropic energy of CRGDS peptides before binding. This approach paves a way for functional reproduction of protein molecules using low-cost synthetic polymers for extensive biomedical applications. [ABSTRACT FROM AUTHOR]

"Biased" G protein-coupled receptor (GPCR) agonists preferentially activate pathways mediated by G proteins or β-arrestins. Here, we use double electron-electron resonance spectroscopy to probe the changes that ligands induce in the conformational distribution of the angiotensin II type I receptor. Monitoring distances between 10 pairs of nitroxide labels distributed across the intracellular regions enabled mapping of four underlying sets of conformations. Ligands from different functional classes have distinct, characteristic effects on the conformational heterogeneity of the receptor. Compared to angiotensin II, the endogenous agonist, agonists with enhanced Gq coupling more strongly stabilize an "open" conformation with an accessible transducer-binding site. β-arrestin-biased agonists deficient in Gq coupling do not stabilize this open conformation but instead favor two more occluded conformations. These data suggest a structural mechanism for biased ligand action at the angiotensin receptor that can be exploited to rationally design GPCR-targeting drugs with greater specificity of action.

Binding kinetics has become a popular topic in pharmacology due to its potential contribution to the selectivity and duration of drug action. Yet, the overall kinetic aspects of complex binding mechanisms are still merely described in terms of elaborate algebraic equations. Interestingly, it has been recommended some 10 years ago to examine such mechanisms in terms of binding fluxes instead of the conventional rate constants. Alike the velocity of product formation in enzymology, those fluxes refer to the velocity by which one target species converts into another one. Novel binding flux‐based approaches are utilized to get a better visual insight into the "competition" between two drugs/ligands for a single target as well as between induced fit‐ and conformational selection pathways for a single ligand within a thermodynamic cycle. The present data were obtained by differential equation‐based simulations. Early on, the ligand‐binding steps "race" to equilibrium (i.e., when their forward and reverse fluxes are equal) at their individual pace. The overall/global equilibrium is only reached later on. For the competition association assays, this parting might produce a transient "overshoot" of one of the bound target species. A similar overshoot may also show up within a thermodynamic cycle and, at first glance, suggest that the induced fit pathway dominates. Yet, present findings show that under certain circumstances, it could rather be the other way round. Novel binding flux‐based approaches offer visually attractive insights into crucial aspects of "complex" binding mechanisms under non‐equilibrium conditions. [ABSTRACT FROM AUTHOR]

[Display omitted] • Cryptic sites potentially expand the space of druggable proteins. • Cryptic sites formed solely by side chain motion don't provide nanomolar ligand binding affinity. • However, cryptic sites formed by loop and domain motion can be valuable drug targets. • The difference in affinity is due to the much shorter time scale of side chain motion. • The different behavior is explained in terms of inhibition models and the kinetics of binding. Cryptic sites can expand the space of druggable proteins, but the potential usefulness of such sites needs to be investigated before any major effort. Given that the binding pockets are not formed, the druggability of such sites is not well understood. The analysis of proteins and their ligands shows that cryptic sites that are formed primarily by the motion of side chains moving out of the pocket to enable ligand binding generally do not bind drug-sized molecules with sufficient potency. By contrast, sites that are formed by loop or hinge motion are potentially valuable drug targets. Arguments are provided to explain the underlying causes in terms of classical enzyme inhibition theory and the kinetics of side chain motion and ligand binding. [ABSTRACT FROM AUTHOR]

Protein homodimers have been classified as three-state or two-state dimers depending on whether a folded monomer forms before association, but the details of the folding–binding mechanisms are poorly understood. Kinetic transition networks of conformational states have provided insight into the folding mechanisms of monomeric proteins, but extending such a network to two protein chains is challenging as all the relative positions and orientations of the chains need to be included, greatly increasing the number of degrees of freedom. Here, we present a simplification of the problem by grouping all states of the two chains into two layers: a dissociated and an associated layer. We combined our two-layer approach with the Wako–Saito–Muñoz–Eaton method and used Transition Path Theory to investigate the dimer formation kinetics of eight homodimers. The analysis reveals a remarkable diversity of dimer formation mechanisms. Induced folding, conformational selection, and rigid docking are often simultaneously at work, and their contribution depends on the protein concentration. Pre-folded structural elements are always present at the moment of association, and asymmetric binding mechanisms are common. Our two-layer network approach can be combined with various methods that generate discrete states, yielding new insights into the kinetics and pathways of flexible binding processes.

How proteins switch between various ligand-free and ligand-bound structures has been a key biophysical question ever since the postulation of the Monod-Wyman-Changeux and Koshland-Nemethy-Filmer models over six decades ago. The ability of NMR spectroscopy to provide structural and kinetic information on biomolecular conformational exchange places it in a unique position as an analytical tool to interrogate the mechanisms of biological processes such as protein folding and biomolecular complex formation. In addition, recent methodological developments in the areas of saturation transfer and relaxation dispersion have expanded the scope of NMR for probing the mechanics of transitions in systems where one or more states constituting the exchange process are sparsely populated and 'invisible' in NMR spectra. In this review, we highlight some of the strategies available from NMR spectroscopy for examining the nature of multi-site conformational exchange, using five case studies that have employed NMR, either in isolation, or in conjunction with other biophysical tools.

The neuropeptide vasopressin (VP) and its three G protein-coupled receptors (V1aR, V1bR and V2R) are of high interest in a wide array of drug discovery programs. V1aR is of particular importance due to its cardiovascular functions and diverse roles in the central nervous system. The structure–activity relationships underpinning ligand-receptor interactions remain however largely unclear, hindering rational drug design. This is not least due to the high structural flexibility of VP in its free as well as receptor-bound states. In this work, we developed a novel approach to reveal features of conformational selectivity upon VP-V1aR complex formation. We employed virtual screening strategies to probe VP’s conformational space for transiently adopted structures that favor binding to V1aR. To this end, we dissected the VP conformational space into three sub-ensembles, each containing distinct structural sets for VP’s three-residue C-terminal tail. We validated the computational results with experimental nuclear magnetic resonance (NMR) data and docked each sub-ensemble to V1aR. We observed that the conformation of VP’s three-residue tail significantly modulated the complex dissociation constants. Solvent-exposed and proline trans-configured VP tail conformations bound to the receptor with three-fold enhanced affinities compared to compacted or cis-configured conformations. The solvent-exposed and more flexible structures facilitated unique interaction patterns between VP and V1aR transmembrane helices 3, 4, and 6 which led to high binding energies. The presented “virtual conformational space screening” approach, integrated with NMR spectroscopy, thus enabled identification and characterization of a conformational selection-type complex formation mechanism that confers novel perspectives on targeting the VP-V1aR interactions at the level of the encounter complex – an aspect that opens novel research avenues for understanding the functionality of the evolutionary selected conformational properties of VP, as well as guidance for ligand design strategies to provide more potent and selective VP analogues.

Glutamate is the major excitatory neurotransmitter in the mammalian central nervous system. After its release from presynaptic nerve terminals, glutamate is quickly removed from the synaptic cleft by excitatory amino acid transporters (EAATs) 1–5, a subfamily of glutamate transporters. The five proteins utilize a complex transport stoichiometry that couples glutamate transport to the symport of three Na+ ions and one H+ in exchange with one K+ to accumulate glutamate against up to 106-fold concentration gradients. They are also anion-selective channels that open and close during transitions along the glutamate transport cycle. EAATs belong to a larger family of secondary-active transporters, the SLC1 family, which also includes purely Na+- or H+-coupled prokaryotic transporters and Na+-dependent neutral amino acid exchangers. In recent years, molecular cloning, heterologous expression, cellular electrophysiology, fluorescence spectroscopy, structural approaches, and molecular simulations have uncovered the molecular mechanisms of coupled transport, substrate selectivity, and anion conduction in EAAT glutamate transporters. Here we review recent findings on EAAT transport mechanisms, with special emphasis on the highly conserved hairpin 2 gate, which has emerged as the central processing unit in many of these functions. [ABSTRACT FROM AUTHOR]

Two mechanisms, induced fit (IF) and conformational selection (CS), have been proposed to explain ligand recognition coupled conformational changes. The histidine binding protein (HisJ) adopts the CS mechanism, in which a pre-equilibrium is established between the open and the closed states with the ligand binding to the closed state. Despite being structurally similar to HisJ, the maltose binding protein (MBP) adopts the IF mechanism, in which the ligand binds the open state and induces a transition to the closed state. To understand the molecular determinants of this difference, we performed molecular dynamics (MD) simulations of coarse-grained dual structure based models. We find that intra-protein contacts unique to the closed state are sufficient to promote the conformational transition in HisJ, indicating a CS-like mechanism. In contrast, additional ligand-mimicking contacts are required to “induce” the conformational transition in MBP suggesting an IF-like mechanism. In agreement with experiments, destabilizing modifications to two structural features, the spine helix (SH) and the balancing interface (BI), present in MBP but absent in HisJ, reduce the need for ligand-mimicking contacts indicating that SH and BI act as structural restraints that keep MBP in the open state. We introduce an SH like element into HisJ and observe that this can impede the conformational transition increasing the importance of ligand-mimicking contacts. Similarly, simultaneous mutations to BI and SH in MBP reduce the barrier to conformational transitions significantly and promote a CS-like mechanism. Together, our results show that structural restraints present in the protein structure can determine the mechanism of conformational transitions and even simple models that correctly capture such structural features can predict their positions. MD simulations of such models can thus be used, in conjunction with mutational experiments, to regulate protein ligand interactions, and modulate ligand binding affinities.

The present study aims at numerically describing to what extent substrate - enzyme complexes in solution may change over time as a natural process of conformational changes for a liganded enzyme in comparison to those movements which occur independently from substrate interaction, i.e. without a ligand. To this end, we selected structurally known pairs of liganded / unliganded CYP450 3A4 enzymes with different geometries hinting at induced fit events. We carried out molecular dynamics simulations (MD) comparing the trajectories in a “cross-over” protocol: (i) we added the ligand to the unliganded crystal form which should adopt geometries similar to the known geometry of the liganded crystal structure during MD, and – conversely – (ii) we removed the bound ligand form the known liganded complex to test if a geometry similar to the known unliganded (apo-) form can be adopted during MD. To compare continues changes we measured root means square deviations and frequencies. Results for case (i) hint at larger conformational changes required for accepting the substrate during its approach to final position – in contrast to case (ii) when mobility is fairly reduced by ligand binding (strain energy). In conclusion, a larger conformational sampling prior to ligand binding and the freezing-in (rigidity) of conformations for bound ligands can be interpreted as two conditions linked to induced-fit.

Binding affinity is a fundamental parameter in drug design, describing the strength of the interaction between a molecule and its target protein. Accurately predicting binding affinity is crucial for the rapid development of novel therapeutics, the prioritization of promising candidates, and the optimization of their properties through rational design strategies. Binding affinity is determined by the mechanism of recognition between proteins and ligands. Various models, including the lock and key, induced fit, and conformational selection, have been proposed to explain this recognition process. However, current computational strategies to predict binding affinity, which are based on these models, have yet to produce satisfactory results. This article explores the connection between binding affinity and these protein-ligand interaction models, highlighting that they offer an incomplete picture of the mechanism governing binding affinity. Specifically, current models primarily center on the binding of the ligand and do not address its dissociation. In this context, the concept of ligand trapping is introduced, which models the mechanisms of dissociation. When combined with the current models, this concept can provide a unified theoretical framework that may allow for the accurate determination of the ligands' binding affinity.

Proteases are key enzymes in viral replication, and interfering with these targets is the basis for therapeutic interventions. We previously introduced a hypothesis about conformational selection in the protease of dengue virus and related flaviviruses, based on conformational plasticity noted in X-ray structures. The present work presents the first functional evidence for alternate recognition by the dengue protease, in a mechanism based primarily on conformational selection rather than induced-fit. Recognition of distinct substrates and inhibitors in proteolytic and binding assays varies to a different extent, depending on factors reported to influence the protease structure. The pH, salinity, buffer type, and temperature cause a change in binding, proteolysis, or inhibition behavior. Using representative inhibitors with distinct structural scaffolds, we identify two contrasting binding profiles to dengue protease. Noticeable effects are observed in the binding assay upon inclusion of a non-ionic detergent in comparison to the proteolytic assay. The findings highlight the impact of the selection of testing conditions on the observed ligand affinity or inhibitory potency. From a broader scope, the dengue protease presents an example, where the induced-fit paradigm appears insufficient to explain binding events with the biological target. Furthermore, this protein reveals the complexity of comparing or combining biochemical assay data obtained under different conditions. This can be particularly critical for artificial intelligence (AI) approaches in drug discovery that rely on large datasets of compounds activity, compiled from different sources using non-identical testing procedures. In such cases, mismatched results will compromise the model quality and its predictive power., Competing Interests: Declaration of competing interest The authors declare no competing interests., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

A new mathematical formula was derived for near equilibrium relaxation processes of enzyme reactions including the conformational selection (CS) modes. CS is one of the most accepted molecular recognition mechanisms, in which protein conformers (CS conformers) are in an equilibrium with varying degrees of ligand binding affinity so that ligands select a particular conformer among them to bind. Using computer simulation techniques, our previous study (Egawa and Callender in Math Biosci 313: 61–70, 2019) predicted that the rate constant for the near equilibrium relaxation processes (kNER) and the concentration-sum of substrate and product (CLt) of enzyme reactions uniquely related to the presence of CS steps in a manner that 1/kNER versus CLt plot transformed from linear to quadratic as the elementary rate constants of inter-conversions among the CS conformers were becoming smaller relative to the rate constants at other steps of the enzymatic reaction system. Thus our previous work could provide a potential tool to detect the presence of CS steps in an enzyme reaction simply by assays using only trace amount of enzyme samples, although logical basis to have the quadratic deformation in the 1/kNER versus CLt plot in response to the presence of CS steps could not be clarified. Employing mathematical approaches that were alternative to those in our previous study, this study succeeded in deriving a theoretical equation that fully explained why and how the CS modes caused the quadratic characters in the 1/kNER versus CLt plot. [ABSTRACT FROM AUTHOR]

Induced fit- (IF) and conformational selection (CS) binding mechanisms have long been regarded to be mutually exclusive. Yet, they are now increasingly considered to produce the final ligand-target complex alongside within a thermodynamic cycle. This viewpoint benefited from the introduction of binding fluxes as a tool for analyzing the overall behavior of such cycle. This study aims to provide more vivid and applicable insights into this emerging field. In this respect, combining differential equation-based simulations and hitherto little explored alternative modes of calculation provide concordant information about the intricate workings of such cycle. In line with previous reports, we observe that the relative contribution of IF increases with the ligand concentration at equilibrium. Yet the baseline contribution may vary from one case to another and simulations as well as calculations show that this parameter is essentially regulated by the dissociation rate of both pathways. Closer attention should be paid to how the contributions of IF and CS compare at physiologically relevant drug/ligand concentrations. To this end, a simple equation discloses how changing a limited set of “microscopic” rate constants can extend the concentration range at which CS contributes most effectively. Finally, it could also be beneficial to extend the utilization of flux- based approaches to more physiologically relevant time scales and alternative binding models. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Abstract Induced fit‐ (IF) and conformational selection (CS) binding mechanisms have long been regarded to be mutually exclusive. Yet, they are now increasingly considered to produce the final ligand‐target complex alongside within a thermodynamic cycle. This viewpoint benefited from the introduction of binding fluxes as a tool for analyzing the overall behavior of such cycle. This study aims to provide more vivid and applicable insights into this emerging field. In this respect, combining differential equation‐ based simulations and hitherto little explored alternative modes of calculation provide concordant information about the intricate workings of such cycle. In line with previous reports, we observe that the relative contribution of IF increases with the ligand concentration at equilibrium. Yet the baseline contribution may vary from one case to another and simulations as well as calculations show that this parameter is essentially regulated by the dissociation rate of both pathways. Closer attention should be paid to how the contributions of IF and CS compare at physiologically relevant drug/ligand concentrations. To this end, a simple equation discloses how changing a limited set of “microscopic” rate constants can extend the concentration range at which CS contributes most effectively. Finally, it could also be beneficial to extend the utilization of flux‐ based approaches to more physiologically relevant time scales and alternative binding models.

Biomolecular recognition often involves conformational changes as a prerequisite for binding (i.e., conformational selection) or concurrently with binding (i.e., induced-fit). Recent advances in structural and kinetic approaches have enabled the detailed characterization of protein motions at atomic resolution. However, to fully understand the role of the conformational dynamics in molecular recognition, studies on the binding transition state are needed. Here, we investigate the binding transition state between nonstructural protein 1 (NS1) of the pandemic 1918 influenza A virus and the human p85β subunit of PI3K. 1918 NS1 binds to p85β via conformational selection. We present the free-energy mapping of the transition and bound states of the 1918 NS1:p85β interaction using linear free energy relationship and ϕ-value analyses. We find that the binding transition state of 1918 NS1 and p85β is structurally similar to the bound state with well-defined binding orientation and hydrophobic interactions. Our finding provides a detailed view of how protein motion contributes to the development of intermolecular interactions along the binding reaction coordinate.

This report describes the implementation of a laboratory exercise for an advanced biochemistry or enzyme kinetics class at the undergraduate or graduate level, designed to improve understanding of protein conformational changes associated with the binding of a ligand. Students measure the fluorescence changes induced by the conformational transition of a glycoprotein (the Na,K-ATPase) upon addition of different ligands (Pi and BeF3−) and analyse the results in order to determine the mechanism of the process. The results show that Pi and BeF3− present opposite effects on the observed rate constants (kobs) with ligand concentration: kobs decreases with [Pi] and increases with [BeF3−]. This observation, together with the frequently used assumption that binding occurs under rapid equilibrium, led to propose different models for ligand-induced conformational transitions: a conformational selection for Pi and an induced fit for BeF3−. In this paper, we show that if the rapid-equilibrium approximation for ligand binding is not assumed, a conformational selection mechanism can account for the effects of both ligands. This active-learning exercise serves as the basis for discussing the consequences of not being extremely cautious when invoking approximations about not-very-well-known systems and the importance of a correct understanding of models assigned to chemical processes. [ABSTRACT FROM AUTHOR]

The proper understanding of biomolecular recognition mechanisms that take place in a drug target is of paramount importance to improve the efficiency of drug discovery and development. The intrinsic dynamic character of proteins has a strong influence on biomolecular recognition mechanisms and models such as conformational selection have been widely used to account for this dynamic association process. However, conformational changes occurring in the receptor prior and upon association with other molecules are diverse and not obvious to predict when only a few structures of the receptor are available. In view of the prominent role of protein flexibility in ligand binding and its implications for drug discovery, it is of great interest to identify receptor conformations that play a major role in biomolecular recognition before starting rational drug design efforts. In this review, we discuss a number of recent advances in computer-aided drug discovery techniques that have been proposed to incorporate receptor flexibility into structure-based drug design. The allowance for receptor flexibility provided by computational techniques such as molecular dynamics simulations or enhanced sampling techniques helps to improve the accuracy of methods used to estimate binding affinities and, thus, such methods can contribute to the discovery of novel drug leads.

Sharks and other cartilaginous fish are the phylogenetically oldest living organisms that have antibodies as part of their adaptive immune system. As part of their humoral adaptive immune response, they produce an immunoglobulin, the so-called immunoglobulin new antigen receptor (IgNAR), a heavy-chain only antibody. The variable domain of an IgNAR, also known as VNAR, binds the antigen as an independent soluble domain. In this study, we structurally and dynamically characterized the affinity maturation mechanism of the germline and somatically matured (PBLA8) VNAR to better understand their function and their applicability as therapeutics. We observed a substantial rigidification upon affinity maturation, which is accompanied by a higher number of contacts, thereby contributing to the decrease in flexibility. Considering the static x-ray structures, the observed rigidification is not obvious, as especially the mutated residues undergo conformational changes during the simulation, resulting in an even stronger network of stabilizing interactions. Additionally, the simulations of the VNAR in complex with the hen egg-white lysozyme show that the VNAR antibodies evidently follow the concept of conformational selection, as the binding-competent state already preexisted even without the presence of the antigen. To have a more detailed description of antibody–antigen recognition, we also present here the binding/unbinding mechanism between the hen egg-white lysozyme and both the germline and matured VNARs. Upon maturation, we observed a substantial increase in the resulting dissociation-free energy barrier. Furthermore, we were able to kinetically and thermodynamically describe the binding process and did not only identify a two-step binding mechanism, but we also found a strong population shift upon affinity maturation toward the native binding pose.

The main protease (Mpro) of SARS-CoV and SARS-CoV-2 is a key enzyme in viral replication and a promising target for the development of antiviral therapeutics. The understanding of this protein is based on a number of observations derived from earlier x-ray structures, which mostly consider substrates or ligands as the main reason behind modulation of the active site. This lead to the concept of substrate-induced subsite cooperativity as an initial attempt to explain the dual binding specificity of this enzyme in recognizing the cleavage sequences at its N- and C-termini, which are important processing steps in obtaining the mature protease. The presented hypothesis proposes that structural heterogeneity is a property of the enzyme, independent of the presence of a substrate or ligand. Indeed, the analysis of Mpro structures of SARS-CoV and SARS-CoV-2 reveals a conformational diversity for the catalytically competent state in ligand-free structures. Variation in the binding site appears to result from flexibility at residues lining the S 1 subpocket and segments incorporating methionine 49 and glutamine 189. The structural evidence introduces "structure-based recognition" as a new paradigm in substrate proteolysis by Mpro. In this concept, the binding space in subpockets of the enzyme varies in a non-cooperative manner, causing distinct conformations, which recognize and process different cleavage sites, as the N- and C-termini. Insights into the recognition basis of the protease provide explanation to the ordered processing of cleavage sites. The hypothesis expands the conformational space of the enzyme and consequently opportunities for drug development and repurposing efforts. Image 1 • Structures of SARS-CoV and SARS-CoV-2 Mpro show conformational heterogeneity. • Variations in the active site of Mpro are ligand- and substrate-independent. • Proteolytic recognition in Mpro is based on structural heterogeneity. • Substrate processing takes place by "structure-based recognition". • The finding expands opportunities for structure-based drug design and repurposing. [ABSTRACT FROM AUTHOR]

The prototypical Ca2+-sensor protein recoverin (Rec) is thought to regulate the activity of rhodopsin kinase (GRK1) in photoreceptors by switching from a relaxed (R) disc membrane-bound conformation in the dark to a more compact, cytosol-diffusing tense (T) conformation upon cell illumination. However, the apparent affinity for Ca2+ of its physiologically relevant form (myristoylated recoverin) is almost two orders of magnitude too low to support this mechanism in vivo. In this work, we compared the individual and synergistic roles of the myristic moiety, the GRK1 target and the disc membrane in modulating the calcium sensitivity of Rec. We show that the sole presence of the target or the disc membrane alone are not sufficient to achieve a physiological response to changes in intracellular [Ca2+]. Instead, the simultaneous presence of GRK1 and membrane allows the T to R transition to occur in a physiological range of [Ca2+] with high cooperativity via a conformational selection mechanism that drives the structural transitions of Rec in the presence of multiple ligands. Our conclusions may apply to other sensory transduction systems involving protein complexes and biological membranes.

Histone deacetylases (HDACs), together with histone acetyltransferases (HATs), regulate gene expression by modulating the acetylation level of chromatin. HDAC3 is implicated in many important cellular processes, particularly in cancer cell proliferation and metastasis, making inhibition of HDAC3 a promising epigenetic treatment for certain cancers. HDAC3 is activated upon complex formation with both inositol tetraphosphate (IP4) and the deacetylase-activating domain (DAD) of multi-protein nuclear receptor corepressor complexes. In previous studies, we have shown that binding of DAD and IP4 to HDAC3 significantly restricts its conformational space towards its stable ternary complex conformation, and suggest this to be the active conformation. Here, we report a single mutation of HDAC3 that is capable of mimicking the stabilizing effects of DAD and IP4, without the presence of either. This mutation, however, results in a total loss of deacetylase activity, prompting a closer evaluation of our understanding of the activation of HDAC3.

Activation of T cells triggers the expression of regulatory molecules like the programmed cell death 1 (PD1) protein. The association of PD1 with the natural ligands PDL1 and PDL2 induces an inhibitory signal that prevents T cells from proliferating and exerting effector functions. However, little is known about how the binding of the ligands induce the PD1 inhibitory signal over T cells effector functions. Here, we explore the dynamics of PD1 free, and in complex with different PDL1 variants as well as the therapeutic antibodies nivolumab and pembrolizumab in order to assess the conformational changes in PD1 related to the signaling process. Our simulations suggest a pre‐conformational selection mechanism for the binding of the different PDL1 variants, while an induced‐fit model fits better for the molecular recognition process of the therapeutic antibodies. A deep analysis of the changes on PD1 movement upon the binding to different ligands revealed that as larger is the difference in the conformation adopted by loop C′D with respect to the complex with PDL1 is higher the ligand ability to reduce the PD1 inhibitory signaling. This behavior suggests that targeting specific conformations of this loop can be useful for designing therapies able to recover T cells effector functions. [ABSTRACT FROM AUTHOR]

The β2 adrenergic receptor (β2AR) is an archetypal G protein coupled receptor (GPCR). One structural signature of GPCR activation is a large-scale movement (ca. 6 to 14 Å) of transmembrane helix 6 (TM6) to a conformation which binds and activates a cognate G protein. The β2AR exhibits a low level of agonist-independent G protein activation. The structural origin of this basal activity and its suppression by inverse agonists is unknown but could involve a unique receptor conformation that promotes G protein activation. Alternatively, a conformational selection model proposes that a minor population of the canonical active receptor conformation exists in equilibrium with inactive forms, thus giving rise to basal activity of the ligand-free receptor. Previous spin-labeling and fluorescence resonance energy transfer experiments designed to monitor the positional distribution of TM6 did not detect the presence of the active conformation of ligand-free β2AR. Here we employ spin-labeling and pressure-resolved double electron-electron resonance spectroscopy to reveal the presence of a minor population of unliganded receptor, with the signature outward TM6 displacement, in equilibrium with inactive conformations. Binding of inverse agonists suppresses this population. These results provide direct structural evidence in favor of a conformational selection model for basal activity in β2AR and provide a mechanism for inverse agonism. In addition, they emphasize 1) the importance of minor populations in GPCR catalytic function; 2) the use of spin-labeling and variable-pressure electron paramagnetic resonance to reveal them in a membrane protein; and 3) the quantitative evaluation of their thermodynamic properties relative to the inactive forms, including free energy, partial molar volume, and compressibility. [ABSTRACT FROM AUTHOR]

A decade ago, many high-affinity drugs were thought to bind to their target via an induced-fit pathway instead of conformational selection. Yet, both pathways make up part of a thermodynamic cycle, and, owing to binding flux-based approaches, it is now rather considered that they act in parallel and also that their relative contribution to the final ligand–target complex depends on the ligand concentration. Those approaches are of increasing interest, but published data still merely refer to the peculiar situation of equilibrium binding. This article draws attention to the benefit of extending those approaches to address more physiological nonequilibrium binding conditions and in vivo situations. For the presented example, they help to apprehend transient experimental manifestations of a 'conventional' thermodynamic cycle. Induced-fit and conformational selection mechanisms for ligand binding were long considered to be mutually exclusive, and most of the slowly dissociating drugs were proposed to act in an induced-fit-like manner. Yet, both pathways form a thermodynamic cycle. Inspection with the aid of binding flux-based approaches contributed to the current view that they act in unison and even that the dominance of one or the other depends on the ligand concentration. Those are now increasingly adopted and even extended. Yet, published data still refer to equilibrium binding, which is hardly met in vivo. This calls for extending those approaches to nonequilibrium situations. Current cycle-based simulations illustrate the utility of such extensions. They highlight their benefit for apprehending transient and occasionally even unexpected manifestations of complex binding mechanisms in vitro and in vivo. [ABSTRACT FROM AUTHOR]

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.

Mutations in RAS oncogenes occur in ~ 30% of human cancers, with KRAS being the most frequently altered isoform. RAS proteins comprise a conserved GTPase domain and a C-terminal lipid-modified tail that is unique to each isoform. The GTPase domain is a 'switch' that regulates multiple signaling cascades that drive cell growth and proliferation when activated by binding GTP, and the signal is terminated by GTP hydrolysis. Oncogenic RAS mutations disrupt the GTPase cycle, leading to accumulation of the activated GTP-bound state and promoting proliferation. RAS is a key target in oncology, however it lacks classic druggable pockets and has been extremely challenging to target. RAS signaling has thus been targeted indirectly, by harnessing key downstream effectors as well as upstream regulators, or disrupting the proper membrane localization required for signaling, by inhibiting either lipid modification or 'carrier' proteins. As a small (20 kDa) protein with multiple conformers in dynamic equilibrium, RAS is an excellent candidate for NMR-driven characterization and screening for direct inhibitors. Several molecules have been discovered that bind RAS and stabilize shallow pockets through conformational selection, and recent compounds have achieved substantial improvements in affinity. NMR-derived insight into targeting the RAS-membrane interface has revealed a new strategy to enhance the potency of small molecules, while another approach has been development of peptidyl inhibitors that bind through large interfaces rather than deep pockets. Remarkable progress has been made with mutation-specific covalent inhibitors that target the thiol of a G12C mutant, and these are now in clinical trials. Here we review the history of RAS inhibitor development and highlight the utility of NMR and integrated biophysical approaches in RAS drug discovery. [ABSTRACT FROM AUTHOR]

Chromatin post‐translational modifications are thought to be important for epigenetic effects on gene expression. Methylation of histone N‐terminal tail lysine residues constitutes one of many such modifications, executed by families of histone lysine methyltransferase (HKMTase). One such protein is ASHH2 from the flowering plant Arabidopsis thaliana, equipped with the interaction domain, CW, and the HKMTase domain, SET. The CW domain of ASHH2 is a selective binder of monomethylation at lysine 4 on histone H3 (H3K4me1) and likely helps the enzyme dock correctly onto chromatin sites. The study of CW and related interaction domains has so far been emphasizing lock–key models, missing important aspects of histone‐tail CW interactions. We here present an analysis of the ASHH2 CW‐H3K4me1 complex using NMR and molecular dynamics, as well as mutation and affinity studies of flexible coils. β‐augmentation and rearrangement of coils coincide with changes in the flexibility of the complex, in particular the η1, η3 and C‐terminal coils, but also in the β1 and β2 strands and the C‐terminal part of the ligand. Furthermore, we show that mutating residues with outlier dynamic behaviour affect the complex binding affinity despite these not being in direct contact with the ligand. Overall, the binding process is consistent with conformational selection. We propose that this binding mechanism presents an advantage when searching for the correct post‐translational modification state among the highly modified and flexible histone tails, and also that the binding shifts the catalytic SET domain towards the nucleosome. Databases: Structural data are available in the PDB database under the accession code 6QXZ. Resonance assignments for CW42 in its apo‐ and holo‐forms are available in the BMRB database under the accession code 27251. [ABSTRACT FROM AUTHOR]

In this work, we investigated how activity and oligomeric state are related in a purified GH1 β‐glucosidase from Spodoptera frugiperda (Sfβgly). Gel filtration chromatography coupled to a multiple angle light scattering detector allowed separation of the homodimer and monomer states and determination of the dimer dissociation constant (KD), which was in the micromolar range. Enzyme kinetic parameters showed that the dimer is on average 2.5‐fold more active. Later, we evaluated the kinetics of homodimerization, scanning the changes in the Sfβgly intrinsic fluorescence over time when the dimer dissociates into the monomer after a large dilution. We described how the rate constant of monomerization (koff) is affected by temperature, revealing the enthalpic and entropic contributions to the process. We also evaluated how the rate constant (kobs) by which equilibrium is reached after dimer dilution behaves when varying the initial Sfβgly concentration. These data indicated that Sfβgly dimerizes through the conformational selection mechanism, in which the monomer undergoes a conformational exchange and then binds to a similar monomer, forming a more active homodimer. Finally, we noted that conformational selection reports and experiments usually rely on a ligand whose concentration is in excess, but for homodimerization, this approach does not hold. Hence, since our approach overcomes this limitation, this study not only is a new contribution to the comprehension of GH1 β‐glucosidases, but it can also help to elucidate protein interaction pathways. [ABSTRACT FROM AUTHOR]

Aminoglycoside acetyltransferases (AACs) catalyze the transfer of an acetyl group between acetyl‐CoA and an aminoglycoside, producing CoA and an acetylated aminoglycoside. AAC(6′)‐Ii enzymes target the amino group linked to the 6′ C atom in an aminoglycoside. Several structures of the AAC(6′)‐Ii from Enterococcus faecium [Ef‐AAC(6′)‐Ii] have been reported to date. However, the detailed mechanism of its enzymatic function remains elusive. In this study, the crystal structure of Ef‐AAC(6′)‐Ii was determined in a novel substrate‐free form. Based on structural analysis, it is proposed that Ef‐AAC(6′)‐Ii sequentially undergoes conformational selection and induced fit for substrate binding. These results therefore provide a novel viewpoint on the mechanism of action of Ef‐AAC(6′)‐Ii. [ABSTRACT FROM AUTHOR]

The formation of specific protein‐protein interactions is often a key to a protein's function. During complex formation, each protein component will undergo a change in the conformational state, for some these changes are relatively small and reside primarily at the sidechain level; however, others may display notable backbone adjustments. One of the classic problems in the protein‐docking field is to be able to a priori predict the extent of such conformational changes. In this work, we investigated three protocols to find the most suitable input structure conformations for cross‐docking, including a robust sampling approach in normal mode space. Counterintuitively, knowledge of the theoretically best combination of normal modes for unbound‐bound transitions does not always lead to the best results. We used a novel spatial partitioning library, Aether Engine (see Supplementary Materials), to efficiently search the conformational states of 56 receptor/ligand pairs, including a recent CAPRI target, in a systematic manner and selected diverse conformations as input to our automated docking server, SwarmDock, a server that allows moderate conformational adjustments during the docking process. In essence, here we present a dynamic cross‐docking protocol, which when benchmarked against the simpler approach of just docking the unbound components shows a 10% uplift in the quality of the top docking pose. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The first x-ray structures of flaviviral proteases defined two conformational states, open and closed, depending on the relative position of NS2B with respect to NS3, a feature that affects the shape of the binding site. The degree of flexibility in the active site was limited to changes in the fold of NS2B rather than NS3 and an induced-fit mechanism was regarded as the main factor for ligand binding. A minor degree of conformational plasticity in NS3 is observed in the two protein chains in the asymmetric unit for the structure of Zika protease with a dipeptide boronate, synthesized in our group. We hypothesize that the NS3 fold has a crucial influence on the shape of the binding site and that a reevaluation of the induced-fit interpretation is warranted. A comparison of flaviviral protease structures identifies conformational dynamics of NS3 and their unexpected role in controlling the depth of the, otherwise shallow, active site. The structural changes of NS3 are mediated by conserved residues and reveal a subpocket, which we denote as subpocket B, extending beyond the catalytic aspartate 75 towards the allosteric binding site, providing a unique connection between the orthosteric and allosteric sites in the protease. The structural evidence supports a molecular recognition based primarily on conformational selection and population shift rather than induced-fit. Besides the implications on protease studies and drug development, this hypothesis provides an interpretation for the alternate binding modes with respect to the catalytic serine, which are observed for recently developed beta-lactam inhibitors incorporating benzyloxyphenylglycine. Image 1 • Conformation dynamics of a conserved domain in the flaviviral NS3 protease allow the creation of a deep-seated cavity. • This added cavity is denoted as subpocket B and connects the allosteric site and the S 1 subpocket in the active site. • The induced-fit model and the open/closed classification based solely on NS2B conformation appear to be of limited utility. • We propose conformational selection as the primary component of ligand recognition of the flaviviral protease. • The hypothesis provides an explanation for the variable positioning of inhibitors relative to the catalytic serine. [ABSTRACT FROM AUTHOR]

Detoxication, or 'drug‐metabolizing', enzymes and drug transporters exhibit remarkable substrate promiscuity and catalytic promiscuity. In contrast to substrate‐specific enzymes that participate in defined metabolic pathways, individual detoxication enzymes must cope with substrates of vast structural diversity, including previously unencountered environmental toxins. Presumably, evolution selects for a balance of 'adequate' kcat/KM values for a wide range of substrates, rather than optimizing kcat/KM for any individual substrate. However, the structural, energetic, and metabolic properties that achieve this balance, and hence optimize detoxication, are not well understood. Two features of detoxication enzymes that are frequently cited as contributions to promiscuity include the exploitation of highly reactive versatile cofactors, or cosubstrates, and a high degree of flexibility within the protein structure. This review examines these intuitive mechanisms in detail and clarifies the contributions of the classic ligand binding models 'induced fit' (IF) and 'conformational selection' (CS) to substrate promiscuity. The available literature data for drug metabolizing enzymes and transporters suggest that IF is exploited by these promiscuous detoxication enzymes, as it is with substrate‐specific enzymes, but the detoxication enzymes uniquely exploit 'IFs' to retain a wide range of substrates at their active sites. In contrast, whereas CS provides no catalytic advantage to substrate‐specific enzymes, promiscuous enzymes may uniquely exploit it to recruit a wide range of substrates. The combination of CS and IF, for recruitment and retention of substrates, can potentially optimize the promiscuity of drug metabolizing enzymes and drug transporters. [ABSTRACT FROM AUTHOR]

Intrinsically disordered proteins (IDPs) are widespread in eukaryotes and participate in a variety of important cellular processes. Numerous studies using state-of-the-art experimental and theoretical methods have advanced our understanding of IDPs and revealed that disordered regions engage in a large repertoire of intra- and intermolecular interactions through their conformational dynamics, thereby regulating many intracellular functions in concert with folded domains. The mechanisms by which IDPs interact with their partners are diverse, depending on their conformational propensities, and include induced fit, conformational selection, and their mixtures. In addition, IDPs are implicated in many diseases, and progress has been made in designing inhibitors of IDP-mediated interactions. Here we review these recent advances with a focus on the dynamics and interactions of IDPs. [ABSTRACT FROM AUTHOR]

The genome in the human cell is extraordinarily compacted in the nucleus. As a result, much of the DNA is inaccessible and functionally inert. Notwithstanding the highly efficient packaging, mechanisms have evolved to render DNA sites accessible that then enable a multitude of factors to carry out ongoing and vital functions. The compaction is derived from DNA complexation within nucleosomes, which can further consolidate into a higher-order chromatin structure. The nucleosome and nucleosomal DNA are not static in nature, but are dynamic, undergoing structural and functional changes as the cell responds to stresses and/or metabolic or environmental cues. We are only beginning to understand the forces and the complexes that engage the nucleosome to unearth the tightly bound and inaccessible DNA sequences and provide an opening to more accessible target sites. In many cases, current findings support a major role for the action of ATP-dependent chromatin remodeling complexes (CRCs) in providing an avenue to factor accessibility that leads to the activation of transcription. The estrogen receptor α (ERα) does not bind to the estrogen response element (ERE) in the canonical nucleosome. However, evidence will be presented that HMGB1 restructures the nucleosome in an ATP-independent manner and also facilitates access and strong binding of ERα to ERE. The features that appear important in the mechanism of action for HMGB1 will be highlighted, in addition to the characteristic features of the restructured nucleosome. These findings, together with previous evidence, suggest a collaborative role for HMGB1 in the step-wise transcription of estrogen-responsive genes. In addition, alternate mechanistic pathways will be discussed, with consideration that “HMGB1 restructuring” of the nucleosome may generally be viewed as a perturbation of the equilibrium of an ensemble of nearly isoenergetic nucleosome states in an energy landscape that is driven by conformational selection by HMGB1.

We present an approach to assess antibody CDR-H3 loops according to their dynamic properties using molecular dynamics simulations. We selected six antibodies in three pairs differing substantially in their individual promiscuity respectively specificity. For two pairs of antibodies crystal structures are available in different states of maturation and used as starting structures for the analyses. For a third pair we chose two antibody CDR sequences obtained from a synthetic library and predicted the respective structures. For all three pairs of antibodies we performed metadynamics simulations to overcome the limitations in conformational sampling imposed by high energy barriers. Additionally, we used classic molecular dynamics simulations to describe nano- to microsecond flexibility and to estimate up to millisecond kinetics of captured conformational transitions. The methodology represents the antibodies as conformational ensembles and allows comprehensive analysis of structural diversity, thermodynamics of conformations and kinetics of structural transitions. Referring to the concept of conformational selection we investigated the link between promiscuity and flexibility of the antibodies' binding interfaces. The obtained detailed characterization of the binding interface clearly indicates a link between structural flexibility and binding promiscuity for this set of antibodies.

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