Selected-fit versus induced-fit protein binding: kinetic differences and mutational analysis

Recognition Mechanisms between a Nanobody and Disordered Epitopes of the Human Prion Protein: An Integrative Molecular Dynamics Study

Immunotherapy using antibodies to target the aggregation of flexible proteins holds promise for therapeutic interventions in neurodegenerative diseases caused by protein misfolding. Prions or PrP^Sc, the causal agents of transmissible spongiform encephalopathies (TSE), represent a model target for immunotherapies as TSE are prototypical protein misfolding diseases. The X-ray crystal structure of the wild-type (WT) human prion protein (HuPrP) bound to a camelid antibody fragment, denoted as Nanobody 484 (Nb484), has been previously solved. Nb484 was found to inhibit prion aggregation in vitro through a unique mechanism of structural stabilization of two disordered epitopes, that is, the palindromic motif (residues 113-120) and the β2-α2 loop region (residues 164-185). The study of the structural basis for antibody recognition of flexible proteins requires appropriate sampling techniques for the identification of conformational states occurring in disordered epitopes. To elucidate the Nb484-HuPrP recognition mechanisms, here we applied molecular dynamics (MD) simulations complemented with available NMR and X-ray crystallography data collected on the WT HuPrP to describe the conformational spaces occurring on HuPrP prior to Nb484 binding. We observe the experimentally determined binding competent conformations within the ensembles of pre-existing conformational states in solution before binding. We also described the Nb484 recognition mechanisms in two HuPrP carrying a polymorphism (E219K) and a TSE-causing mutation (V210I). Our hybrid approaches allow the identification of dynamic conformational landscapes existing on HuPrP and highly characterized by molecular disorder to identify physiologically relevant and druggable transitions.

Mechanism of activation and the rewired network: New drug design concepts

Precision oncology benefits from effective early phase drug discovery decisions. Recently, drugging inactive protein conformations has shown impressive successes, raising the cardinal questions of which targets can profit and what are the principles of the active/inactive protein pharmacology. Cancer driver mutations have been established to mimic the protein activation mechanism. We suggest that the decision whether to target an inactive (or active) conformation should largely rest on the protein mechanism of activation. We next discuss the recent identification of double (multiple) same-allele driver mutations and their impact on cell proliferation and suggest that like single driver mutations, double drivers also mimic the mechanism of activation. We further suggest that the structural perturbations of double (multiple) in cis mutations may reveal new surfaces/pockets for drug design. Finally, we underscore the preeminent role of the cellular network which is deregulated in cancer. Our structure-based review and outlook updates the traditional Mechanism of Action, informs decisions, and calls attention to the intrinsic activation mechanism of the target protein and the rewired tumor-specific network, ushering innovative considerations in precision medicine.

Cyclosporin Structure and Permeability: From A to Z and Beyond

Cyclosporins are natural or synthetic undecapeptides with a wide range of actual and potential pharmaceutical applications. Several members of the cyclosporin compound family have remarkably high passive membrane permeabilities that are not well-described by simple structural metrics. Here we review experimental studies of cyclosporin structure and permeability, including cyclosporin-metal complexes. We also discuss models for the conformation-dependent permeability of cyclosporins and similar compounds. Finally, we identify current knowledge gaps in the literature and provide recommendations regarding future avenues of exploration.

What's in a name? From "fluctuation fit" to "conformational selection": rediscovery of a concept

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.

Measuring and Analyzing Binding Kinetics of Coupled Folding and Binding Reactions Under Pseudo-First-Order Conditions

Many intrinsically disordered proteins (IDPs) adopt a well-defined structure upon binding to their interaction partners. Kinetic characterization is a requirement for the investigation of the dynamics and mechanisms of these folding-upon-binding reactions. Here a protocol is described for the investigation of binding kinetics of bimolecular binding and folding reactions of IDPs to their ligand partner under pseudo-first-order conditions using stopped-flow mixing and fluorescence detection.

Protein structural heterogeneity: A hypothesis for the basis of proteolytic recognition by the main protease of SARS-CoV and SARS-CoV-2

The main protease (M^pro) 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 M^pro 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₁ subpocket and segments incorporating methionine 49 and glutamine 189. The structural evidence introduces “structure-based recognition” as a new paradigm in substrate proteolysis by M^pro. 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.

Bringing the Ca2+ sensitivity of myristoylated recoverin into the physiological range

The prototypical Ca²⁺-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 Ca²⁺ 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 [Ca²⁺]. Instead, the simultaneous presence of GRK1 and membrane allows the T to R transition to occur in a physiological range of [Ca²⁺] 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.

Mechanisms of ligand binding

Many processes in chemistry and biology involve interactions of a ligand with its molecular target. Interest in the mechanism governing such interactions has dominated theoretical and experimental analysis for over a century. The interpretation of molecular recognition has evolved from a simple rigid body association of the ligand with its target to appreciation of the key role played by conformational transitions. Two conceptually distinct descriptions have had a profound impact on our understanding of mechanisms of ligand binding. The first description, referred to as induced fit, assumes that conformational changes follow the initial binding step to optimize the complex between the ligand and its target. The second description, referred to as conformational selection, assumes that the free target exists in multiple conformations in equilibrium and that the ligand selects the optimal one for binding. Both descriptions can be merged into more complex reaction schemes that better describe the functional repertoire of macromolecular systems. This review deals with basic mechanisms of ligand binding, with special emphasis on induced fit, conformational selection, and their mathematical foundations to provide rigorous context for the analysis and interpretation of experimental data. We show that conformational selection is a surprisingly versatile mechanism that includes induced fit as a mathematical special case and even captures kinetic properties of more complex reaction schemes. These features make conformational selection a dominant mechanism of molecular recognition in biology, consistent with the rich conformational landscape accessible to biological macromolecules being unraveled by structural biology.

Thermodynamic and Mechanistic Insights into Coupled Binding and Unwinding of Collagen by Matrix Metalloproteinase 1

Local unwinding of the collagen triple helix is a necessary step for initiating the collagen degradation cascade in extracellular matrices. A few matrix metalloproteinases (MMPs) are known to support this key process, but its energetic aspects remain unknown. Here, we captured the thermodynamics of the triple helix unwinding by monitoring interactions between a collagen peptide and MMP-1(E200A) - an active-site mutant of an archetypal vertebrate collagenase - at increasing temperatures, using isothermal titration calorimetry (ITC). Coupled binding and unwinding manifests as a curved relationship between the total enthalpy change and temperature of the reaction, producing increasingly negative heat capacity change (ΔΔC_p ≈ -36.3 kcal/molK²). A specially designed solid-phase binding and cleavage assay (SPBCA) reported strain in the catalytically relevant unwound state, suggesting that this state is distinct from the horizon of sampled conformations of the collagenase-susceptible site. MMP-1 appears to blend selected fit with induced fit mechanisms to catalyse collagen unwinding prior to cleavage of individual collagen chains.

Structure dictates the mechanism of ligand recognition in the histidine and maltose binding proteins

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.

•

MBP operates by induced fit, HisJ by the conformational selection mechanism.

•

Dual structure based models (dSBMs) encode two structures of a protein.

•

MD simulations of dSBMs can identify the mechanism of conformational transitions.

•

Locks, absent in HisJ, hold MBP open with ligand contacts required for closing.

•

Binding mechanisms can be modified by altering such structural locks.

Prospecting for new catechol-O-methyltransferase (COMT) inhibitors as a potential treatment for Parkinson's disease: a study by molecular dynamics and structure-based virtual screening

Parkinson's disease (PD) is a neurodegenerative, chronic, and progressive disease, common in the elderly. The catechol-O-methyltransferase (COMT) is a monomeric enzyme involved in dopamine (DA) degradation, the neurotransmitter in deficit in patients with PD. The reference treatment of PD consists of levodopa (L-dopa) administration, which is the precursor of DA. The inhibition of COMT is an adjuvant treatment in PD since it keeps DA levels constant. The goal of this study was to identify drug candidates capable of inhibiting COMT for the treatment of PD and identify important fragments of these molecules. Initially, we analyzed the flexibility of COMT and defined its main conformations in solution regarding the absence (system I) and presence of the S-adenosyl-L-methionine (SAM) cofactor (system II) through molecular dynamics (MD) simulations. Two regions in these structures were selected for molecular docking, firstly the entire cavity where the cofactor and substrates are bound and secondly the specific biding region of the enzyme substrates. Based on the conformations of the MD, the virtual screening (VS) was performed against FDA Approved and Zinc Natural Products databases aiming at the selection of the best compounds. Subsequently, the absorption, distribution, metabolization, excretion, and toxicity (ADMET) properties, as well as drug-score and drug-likeness indexes of the most promising compounds were analyzed. After a detailed analysis of the compounds selected by structure-based VS, it was possible to highlight the fragments most frequently involved in their stability: 2,3,4,9-tetrahydro-1H-pyrido[3,4-b]indole, 9H-Benz(c)indole(3,2,1-ij)(1,5)naphthyridin-9-one and (10R,13S)-10,13-dimethyl-1,2,6,7,8,9,11,12,14,15,16,17dodecahydrocyclopenta[a]phenanthren-3-one. The identification of these potential fragments is essential for the prospection of more specific inhibitors against COMT using the technique of Fragment-based lead discovery (FBLD). Besides, this study allowed us to identify the potential COMT inhibitors through a complete understanding of molecular-level interactions based on the flexibility of this protein.Communicated by Ramaswamy H. Sarma.

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

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

Structure-Based Analysis of Cryptic-Site Opening

Many proteins in their unbound structures have cryptic sites that are not appropriately sized for drug binding. We consider here 32 proteins from the recently published CryptoSite set with validated cryptic sites, and study whether the sites remain cryptic in all available X-ray structures of the proteins solved without any ligand bound near the sites. It was shown that only few of these proteins have binding pockets that never form without ligand binding. Sites that are cryptic in some structures but spontaneously form in others are also rare. In most proteins the forming of pockets is affected by mutations or ligand binding at locations far from the cryptic site. To further explore these mechanisms, we applied adiabatic biased molecular dynamics simulations to guide the proteins from their ligand-free structures to ligand-bound conformations, and studied the distribution of druggability scores of the pockets located at the cryptic sites.

Features of molecular recognition of intrinsically disordered proteins via coupled folding and binding

The sequence-structure-function paradigm of proteins has been revolutionized by the discovery of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs). In contrast to traditional ordered proteins, IDPs/IDRs are unstructured under physiological conditions. The absence of well-defined three-dimensional structures in the free state of IDPs/IDRs is fundamental to their function. Folding upon binding is an important mode of molecular recognition for IDPs/IDRs. While great efforts have been devoted to investigating the complex structures and binding kinetics and affinities, our knowledge on the binding mechanisms of IDPs/IDRs remains very limited. Here, we review recent advances on the binding mechanisms of IDPs/IDRs. The structures and kinetic parameters of IDPs/IDRs can vary greatly, and the binding mechanisms can be highly dependent on the structural properties of IDPs/IDRs. IDPs/IDRs can employ various combinations of conformational selection and induced fit in a binding process, which can be templated by the target and/or encoded by the IDP/IDR. Further studies should provide deeper insights into the molecular recognition of IDPs/IDRs and enable the rational design of IDP/IDR binding mechanisms in the future.

The kink-turn in the structural biology of RNA

The kink-turn (k-turn) is a widespread structural motif found in functional RNA species. It typically comprises a three-nucleotide bulge followed by tandem trans sugar edge-Hoogsteen G:A base pairs. It introduces a sharp kink into the axis of duplex RNA, juxtaposing the minor grooves. Cross-strand H-bonds form at the interface, accepted by the conserved adenine nucleobases of the G:A basepairs. Alternative acceptors for one of these divides the k-turns into two conformational classes N3 and N1. The base pair that follows the G:A pairs (3b:3n) determines which conformation is adopted by a given k-turn. k-turns often mediate tertiary contacts in folded RNA species and frequently bind proteins. Common k-turn binding proteins include members of the L7Ae family, such as the human 15·5k protein. A recognition helix within these proteins binds in the widened major groove on the outside of the k-turn, that makes specific H-bonds with the conserved guanine nucleobases of the G:A pairs. L7Ae binds with extremely high affinity, and single-molecule data are consistent with folding by conformational selection. The standard, simple k-turn can be elaborated in a variety of ways, that include the complex k-turns and the k-junctions. In free solution in the absence of added metal ions or protein k-turns do not adopt the tightly-kinked conformation. They undergo folding by the binding of proteins, by the formation of tertiary contacts, and some (but not all) will fold on the addition of metal ions. Whether or not folding occurs in the presence of metal ions depends on local sequence, including the 3b:3n position, and the -1b:-1n position (5' to the bulge). In most cases -1b:-1n = C:G, so that the 3b:3n position is critical since it determines both folding properties and conformation. In general, the selection of these sequence matches a given k-turn to its biological requirements. The k-turn structure is now very well understood, to the point at which they can be used as a building block for the formation of RNA nano-objects, including triangles and squares.

Binding of a Telomestatin Derivative Changes the Mechanical Anisotropy of a Human Telomeric G-Quadruplex

Mechanical anisotropy is an essential property for biomolecules to assume structural and functional roles in mechanobiology. However, there is insufficient information on the mechanical anisotropy of ligand-biomolecule complexes. Herein, we investigated the mechanical property of individual human telomeric G-quadruplexes bound to telomestatin, using optical tweezers. Stacking of the ligand to the G-tetrad planes changes the conformation of the G-quadruplex, which resembles a balloon squeezed in certain directions. Such a squeezed balloon effect strengthens the G-tetrad planes, but dislocates and weakens the loops in the G-quadruplex upon ligand binding. These dynamic interactions indicate that the binding between the ligand and G-quadruplex follows the induced-fit model. We anticipate that the altered mechanical anisotropy of the ligand-G-quadruplex complex can add additional level of regulations on the motor enzymes that process DNA or RNA molecules.

An integrated perspective on RNA aptamer ligand-recognition models: clearing muddy waters

Riboswitches are short RNA motifs that sensitively and selectively bind cognate ligands to modulate gene expression. Like protein receptor-ligand pairs, their binding dynamics are traditionally categorized as following one of two paradigmatic mechanisms: conformational selection and induced fit. In conformational selection, ligand binding stabilizes a particular state already present in the receptor's dynamic ensemble. In induced fit, ligand-receptor interactions enable the system to overcome the energetic barrier into a previously inaccessible state. In this article, we question whether a polarized division of RNA binding mechanisms truly meets the conceptual needs of the field. We will review the history behind this classification of RNA-ligand interactions, and the way induced fit in particular has been rehabilitated by single-molecule studies of RNA aptamers. We will highlight several recent results from single-molecule experimental studies of riboswitches that reveal gaps or even contradictions between common definitions of the two terms, and we will conclude by proposing a more robust framework that considers the range of RNA behaviors unveiled in recent years as a reality to be described, rather than an increasingly unwieldy set of exceptions to the traditional models.

Network of Conformational Transitions Revealed by Molecular Dynamics Simulations of the Carbonic Anhydrase II Apo-Enzyme

Human carbonic anhydrase II (HCA II) is an enzyme that catalyzes the reversible hydration of CO2 into bicarbonate (HCO3-) and a proton (H+) as well as other reactions at an extremely high rate. This enzyme plays fundamental roles in human physiology/pathology, such as controlling the pH level in cells and so on. However, the binding mechanism between apo-HCA II and CO2 or other ligands as well as related conformational changes remains poorly understood, and atomic investigation into it could promote our understanding of related internal physiological/pathological mechanisms. In this study, long-time atomic molecular dynamics simulations as well as the clustering and free-energy analysis were performed to reveal the dynamics of apo-HCA II as well as the mechanism upon ligand binding. Our simulations indicate that the crystallographic B-factors considerably underestimate the loop dynamics: multiple conformations can be adopted by loops 1 and 2, especially for loop 1 because loop 1 is one side of the binding pocket, and its left-to-right movement can compress or extend the binding pocket, leading to one inactive (closed) state, three intermediate (semiopen) states, and one active (open) state; CO2 cannot get into the binding pocket of the inactive state but can get into those of intermediate and active states. The coexistence of multiple conformational states proposes a possible conformational selection model for the binding mechanism between apo-HCA II and CO2 or other ligands, revising our previous view of its functional mechanism of conformational change upon ligand binding and offering valuable structural insights into the workings of HCA II.

Intact Transition Epitope Mapping (ITEM)

Intact transition epitope mapping (ITEM) enables rapid and accurate determination of protein antigen-derived epitopes by either epitope extraction or epitope excision. Upon formation of the antigen peptide-containing immune complex in solution, the entire mixture is electrosprayed to translate all constituents as protonated ions into the gas phase. There, ions from antibody-peptide complexes are separated from unbound peptide ions according to their masses, charges, and shapes either by ion mobility drift or by quadrupole ion filtering. Subsequently, immune complexes are dissociated by collision induced fragmentation and the ion signals of the "complex-released peptides," which in effect are the epitope peptides, are recorded in the time-of-flight analyzer of the mass spectrometer. Mixing of an antibody solution with a solution in which antigens or antigen-derived peptides are dissolved is, together with antigen proteolysis, the only required in-solution handling step. Simplicity of sample handling and speed of analysis together with very low sample consumption makes ITEM faster and easier to perform than other experimental epitope mapping methods. Graphical Abstract ᅟ.

A Network of Conformational Transitions in the Apo Form of NDM-1 Enzyme Revealed by MD Simulation and a Markov State Model

New Delhi metallo-β-lactamase-1 (NDM-1) is a novel β-lactamase enzyme that confers enteric bacteria with nearly complete resistance to all β-lactam antibiotics, so it raises a formidable and global threat to human health. However, the binding mechanism between apo-NDM-1 and antibiotics as well as related conformational changes remains poorly understood, which largely hinders the overcoming of its antibiotic resistance. In our study, long-time conventional molecular dynamics simulation and Markov state models were applied to reveal both the dynamical and conformational landscape of apo-NDM-1: the MD simulation demonstrates that loop L3, which is responsible for antibiotic binding, is the most flexible and undergoes dramatic conformational changes; moreover, the Markov state model built from the simulation maps four metastable states including open, semiopen, and closed conformations of loop L3 as well as frequent transitions between the states. Our findings propose a possible conformational selection model for the binding mechanism between apo-NDM-1 and antibiotics, which facilitates the design of novel inhibitors and antibiotics.

Protein-ligand (un)binding kinetics as a new paradigm for drug discovery at the crossroad between experiments and modelling

In the last three decades, protein and nucleic acid structure determination and comprehension of the mechanisms, leading to their physiological and pathological functions, have become a cornerstone of biomedical sciences. A deep understanding of the principles governing the fates of cells and tissue at the molecular level has been gained over the years, offering a solid basis for the rational design of drugs aimed at the pharmacological treatment of numerous diseases. Historically, affinity indicators (i.e. Kd and IC50/EC50) have been assumed to be valid indicators of the in vivo efficacy of a drug. However, recent studies pointed out that the kinetics of the drug-receptor binding process could be as important or even more important than affinity in determining the drug efficacy. This eventually led to a growing interest in the characterisation and prediction of the rate constants of protein-ligand association and dissociation. For instance, a drug with a longer residence time can kinetically select a given receptor over another, even if the affinity for both receptors is comparable, thus increasing its therapeutic index. Therefore, understanding the molecular features underlying binding and unbinding processes is of central interest towards the rational control of drug binding kinetics. In this review, we report the theoretical framework behind protein-ligand association and highlight the latest advances in the experimental and computational approaches exploited to investigate the binding kinetics.

A Network of Conformational Transitions Revealed by Molecular Dynamics Simulations of the Binary Complex of Escherichia coli 6-Hydroxymethyl-7,8-dihydropterin Pyrophosphokinase with MgATP

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the first reaction in the folate biosynthetic pathway. Comparison of its X-ray and nuclear magnetic resonance structures suggests that the enzyme undergoes significant conformational change upon binding to its substrates, especially in three catalytic loops. Experimental research has shown that, in its binary form, even bound by analogues of MgATP, loops 2 and 3 remain rather flexible; this raises questions about the putative large-scale induced-fit conformational change of the HPPK-MgATP binary complex. In this work, long-time all-atomic molecular dynamics simulations were conducted to investigate the loop dynamics in this complex. Our simulations show that, with loop 3 closed, multiple conformations of loop 2, including the open, semiopen, and closed forms, are all accessible to the binary complex. These results provide valuable structural insights into the details of conformational changes upon 6-hydroxymethyl-7,8-dihydropterin (HP) binding and biological activities of HPPK. Conformational network analysis and principal component analysis related to the loops are also discussed.

How to Distinguish Conformational Selection and Induced Fit Based on Chemical Relaxation Rates

Protein binding often involves conformational changes. Important questions are whether a conformational change occurs prior to a binding event (‘conformational selection’) or after a binding event (‘induced fit’), and how conformational transition rates can be obtained from experiments. In this article, we present general results for the chemical relaxation rates of conformational-selection and induced-fit binding processes that hold for all concentrations of proteins and ligands and, thus, go beyond the standard pseudo-first-order approximation of large ligand concentration. These results allow to distinguish conformational-selection from induced-fit processes—also in cases in which such a distinction is not possible under pseudo-first-order conditions—and to extract conformational transition rates of proteins from chemical relaxation data.

The function of proteins is affected by their conformational dynamics, i.e. by transitions between lower-energy ground-state conformations and higher-energy excited-state conformations of the proteins. Advanced NMR and single-molecule experiments indicate that higher-energy conformations in the unbound state of proteins can be similar to ground-state conformations in the bound state, and vice versa. These experiments illustrate that the conformational change of a protein during binding may occur before a binding event, rather than being induced by this binding event. However, determining the temporal order of conformational transitions and binding events typically requires additional information from chemical relaxation experiments that probe the relaxation kinetics of a mixture of proteins and ligands into binding equilibrium. These chemical relaxation experiments are usually performed and analysed at ligand concentrations that are much larger than the protein concentrations. At such high ligand concentrations, the temporal order of conformational transitions and binding events can only be inferred in special cases. In this article, we present general equations that describe the dominant chemical relaxation kinetics for all protein and ligand concentrations. Our general equations allow to clearly infer from relaxation data whether a conformational transition occurs prior to a binding event, or after the binding event.

Toward High-Throughput Predictive Modeling of Protein Binding/Unbinding Kinetics

One of the unaddressed challenges in drug discovery is that drug potency determined in vitro is not a reliable indicator of drug activity in vivo. Accumulated evidence suggests that in vivo activity is more strongly correlated with the binding/unbinding kinetics than the equilibrium thermodynamics of protein-ligand interactions (PLIs). However, existing experimental and computational techniques are insufficient in studying the molecular details of kinetics processes of PLIs on a large scale. Consequently, we not only have limited mechanistic understanding of the kinetic processes but also lack a practical platform for high-throughput screening and optimization of drug leads on the basis of their kinetic properties. For the first time, we address this unmet need by integrating coarse-grained normal mode analysis with multitarget machine learning (MTML). To test our method, HIV-1 protease is used as a model system. We find that computational models based on the residue normal mode directionality displacement of PLIs can not only recapitulate the results from all-atom molecular dynamics simulations but also predict protein-ligand binding/unbinding kinetics accurately. When this is combined with energetic features, the accuracy of combined kon and koff prediction reaches 74.35%. Furthermore, our integrated model provides us with new insights into the molecular determinants of the kinetics of PLIs. We propose that the coherent coupling of conformational dynamics and thermodynamic interactions between the receptor and the ligand may play a critical role in determining the kinetic rate constants of PLIs. In conclusion, we demonstrate that residue normal mode directionality displacement can serve as a kinetic fingerprint to capture long-time-scale conformational dynamics of the binding/unbinding kinetics. When this is coupled with MTML, it is possible to screen and optimize compounds on the basis of their binding/unbinding kinetics in a high-throughput fashion. The further development of such computational tools will bridge one of the critical missing links between in vitro compound screening and in vivo drug activity.

Evolutionary Conserved Positions Define Protein Conformational Diversity

Conformational diversity of the native state plays a central role in modulating protein function. The selection paradigm sustains that different ligands shift the conformational equilibrium through their binding to highest-affinity conformers. Intramolecular vibrational dynamics associated to each conformation should guarantee conformational transitions, which due to its importance, could possibly be associated with evolutionary conserved traits. Normal mode analysis, based on a coarse-grained model of the protein, can provide the required information to explore these features. Herein, we present a novel procedure to identify key positions sustaining the conformational diversity associated to ligand binding. The method is applied to an adequate refined dataset of 188 paired protein structures in their bound and unbound forms. Firstly, normal modes most involved in the conformational change are selected according to their corresponding overlap with structural distortions introduced by ligand binding. The subspace defined by these modes is used to analyze the effect of simulated point mutations on preserving the conformational diversity of the protein. We find a negative correlation between the effects of mutations on these normal mode subspaces associated to ligand-binding and position-specific evolutionary conservations obtained from multiple sequence-structure alignments. Positions whose mutations are found to alter the most these subspaces are defined as key positions, that is, dynamically important residues that mediate the ligand-binding conformational change. These positions are shown to be evolutionary conserved, mostly buried aliphatic residues localized in regular structural regions of the protein like β-sheets and α-helix.

Discrimination between conformational selection and induced fit protein-ligand binding using Integrated Global Fit analysis

Molecular recognition between proteins and small molecule ligands is at the heart of biological function in cellular systems and the basis of modern rational drug development. Therefore, the mechanisms governing protein-ligand interaction have been objects of research for many decades. The last 15 years has seen a revival of a discussion whether conformational selection (CS) or induced fit (IF) is the most relevant binding mechanism. A decreasing observed rate constant, k obs, with increasing ligand concentration was considered to be a hallmark of CS, but according to contemporary knowledge, a positive saturating behavior of k obs can be explained by both CS and IF mechanisms. The only currently recognized kinetic method to differentiate between both binding mechanisms includes the measurement of two separate series of binding kinetics with variation of either protein or ligand under pseudo-first-order conditions. This study avoids the disadvantage of high protein concentrations and provides evidence that a comprehensive Integrated Global Fit analysis of sets of binding kinetics with just varied ligand concentration in combination with equilibrium data and optional displacement kinetics can effectively differentiate between CS and IF binding mechanisms. The limiting situation, when physical binding dominates over the previous (CS) or subsequent (IF) conformational changes, is carefully analyzed. Finally, the relevance of kinetic methods and the elucidation of more complex binding mechanisms are discussed for advanced rational selection and optimization of drug candidates.

What is so special about smell? Olfaction as a model system in neurobiology

Neurobiology studies mechanisms of cell signalling. A key question is how cells recognise specific signals. In this context, olfaction has become an important experimental system over the past 25 years. The olfactory system responds to an array of structurally diverse stimuli. The discovery of the olfactory receptors (ORs), recognising these stimuli, established the olfactory pathway as part of a greater group of signalling mechanisms mediated by G-protein-coupled receptors (GPCRs). GPCRs are the largest protein family in the mammalian genome and involved in numerous fundamental physiological processes. The OR family exhibits two characteristics that make them an excellent model system to understand GPCRs: its size and the structural diversity of its members. Research on the OR binding site investigates what amino acid sequences determine the receptor-binding capacity. This promises a better understanding of how the basic genetic makeup of GPCRs relates to their diversification in ligand-binding capacities.

Sampling of conformational ensemble for virtual screening using molecular dynamics simulations and normal mode analysis

Aim: Molecular dynamics simulations and normal mode analysis are well-established approaches to generate receptor conformational ensembles (RCEs) for ligand docking and virtual screening. Here, we report new fast molecular dynamics-based and normal mode analysis-based protocols combined with conformational pocket classifications to efficiently generate RCEs. Materials & Methods: We assessed our protocols on two well-characterized protein targets showing local active site flexibility, dihydrofolate reductase and large collective movements, CDK2. The performance of the RCEs was validated by distinguishing known ligands of dihydrofolate reductase and CDK2 among a dataset of diverse chemical decoys. Results & discussion: Our results show that different simulation protocols can be efficient for generation of RCEs depending on different kind of protein flexibility.[Formula: see text]

Molecular dynamics simulations of the Escherichia coli HPPK apo-enzyme reveal a network of conformational transitions

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the first reaction in the folate biosynthetic pathway. Comparison of its X-ray and nuclear magnetic resonance structures suggests that the enzyme undergoes significant conformational change upon binding to its substrates, especially in three catalytic loops. Experimental research has shown that even when confined by crystal contacts, loops 2 and 3 remain rather flexible when the enzyme is in its apo form, raising questions about the putative large-scale induced-fit conformational change of HPPK. To investigate the loop dynamics in a crystal-free environment, we performed conventional molecular dynamics simulations of the apo-enzyme at two different temperatures (300 and 350 K). Our simulations show that the crystallographic B-factors considerably underestimate the loop dynamics; multiple conformations of loops 2 and 3, including the open, semi-open, and closed conformations that an enzyme must adopt throughout its catalytic cycle, are all accessible to the apo-enzyme. These results revise our previous view of the functional mechanism of conformational change upon MgATP binding and offer valuable structural insights into the workings of HPPK. In this paper, conformational network analysis and principal component analysis related to the loops are discussed to support the presented conclusions.

Biochemical Basis for Dominant Inheritance, Variable Penetrance, and Maternal Effects in RBP4 Congenital Eye Disease

Gestational vitamin A (retinol) deficiency poses a risk for ocular birth defects and blindness. We identified missense mutations in RBP4, encoding serum retinol binding protein, in three families with eye malformations of differing severity, including bilateral anophthalmia. The mutant phenotypes exhibit dominant inheritance, but incomplete penetrance. Maternal transmission significantly increases the probability of phenotypic expression. RBP normally delivers retinol from hepatic stores to peripheral tissues, including the placenta and fetal eye. The disease mutations greatly reduce retinol binding to RBP, yet paradoxically increase the affinity of RBP for its cell surface receptor, STRA6. By occupying STRA6 nonproductively, the dominant-negative proteins disrupt vitamin A delivery from wild-type proteins within the fetus, but also, in the case of maternal transmission, at the placenta. These findings establish a previously uncharacterized mode of maternal inheritance, distinct from imprinting and oocyte-derived mRNA, and define a group of hereditary disorders plausibly modulated by dietary vitamin A.

Biochemical characterization of smoothened receptor antagonists by binding kinetics against drug-resistant mutant

Hedgehog (Hh) signaling critical for development, differentiation, and cell growth is involved in several cancers, including medulloblastoma and basal cell carcinoma. Although antagonism of the smoothened receptor (SMO), which mediates Hh signaling, is an attractive therapeutic target, a drug-resistant mutation in SMO (SMO-D473H) was identified in a clinical trial of the approved drug vismodegib. TAK-441 potently inhibits SMO-D473H, unlike vismodegib and another SMO antagonist, cyclopamine, whereas the differences in binding modes between these antagonists remain unknown. Here we report the biochemical characterization of TAK-441, vismodegib, and cyclopamine by binding kinetics. The association (kon) and dissociation (koff) rates were determined by kinetic binding studies using [(3)H]TAK-441, and dissociation was confirmed by label-free affinity selection-mass spectrometry (AS-MS). In the [(3)H]TAK-441 competition assay, TAK-441 but not vismodegib and cyclopamine showed time-dependent inhibition. Quantitative kinetic binding analysis revealed that koff of TAK-441 was >10-fold smaller than those of vismodegib and cyclopamine. To further assess the binding mode of antagonists, kinetic binding analysis was performed against SMO-D473H. The D473H mutation affected koff of TAK-441 but not kon. In contrast, only kon was changed by the D473H mutation in the case of vismodegib and cyclopamine. These results suggest that the difference in antagonist efficacy against D473H is associated with the binding mode of antagonists. These findings provide a new insight into the drug action of SMO antagonists and help develop potential therapeutics for drug-resistant mutants.

Protein Flexibility in Drug Discovery: From Theory to Computation

Nowadays it is widely accepted that the mechanisms of biomolecular recognition are strongly coupled to the intrinsic dynamic of proteins. In past years, this evidence has prompted the development of theoretical models of recognition able to describe ligand binding assisted by protein conformational changes. On a different perspective, the need to take into account protein flexibility in structure-based drug discovery has stimulated the development of several and extremely diversified computational methods. Herein, on the basis of a parallel between the major recognition models and the simulation strategies used to account for protein flexibility in ligand binding, we sort out and describe the most innovative and promising implementations for structure-based drug discovery.

Demonstration of a folding after binding mechanism in the recognition between the measles virus NTAIL and X domains

In the past decade, a wealth of experimental data has demonstrated that a large fraction of proteins, while functional, are intrinsically disordered at physiological conditions. Many intrinsically disordered proteins (IDPs) undergo a disorder-to-order transition upon binding to their biological targets, a phenomenon known as induced folding. Induced folding may occur through two extreme mechanisms, namely conformational selection and folding after binding. Although the pre-existence of ordered structures in IDPs is a prerequisite for conformational selection, it does not necessarily commit to this latter mechanism, and kinetic studies are needed to discriminate between the two possible scenarios. So far, relatively few studies have addressed this issue from an experimental perspective. Here, we analyze the interaction kinetics between the intrinsically disordered C-terminal domain of the measles virus nucleoprotein (NTAIL) and the X domain (XD) of the viral phosphoprotein. Data reveal that NTAIL recognizes XD by first forming a weak encounter complex in a disordered conformation, which is subsequently locked-in by a folding step; i.e., binding precedes folding. The implications of our kinetic results, in the context of previously reported equilibrium data, are discussed. These results contribute to enhancing our understanding of the molecular mechanisms by which IDPs recognize their partners and represent a paradigmatic example of the need of kinetic methods to discriminate between reaction mechanisms.

Disruption of integrin-fibronectin complexes by allosteric but not ligand-mimetic inhibitors

Failure of Arg-Gly-Asp (RGD)-based inhibitors to reverse integrin-ligand binding has been reported, but the prevalence of this phenomenon among integrin heterodimers is currently unknown. In the present study we have investigated the interaction of four different RGD-binding integrins (α5β1, αVβ1, αVβ3 and αVβ6) with fibronectin (FN) using surface plasmon resonance. The ability of inhibitors to reverse ligand binding was assessed by their capacity to increase the dissociation rate of pre-formed integrin-FN complexes. For all four receptors we showed that RGD-based inhibitors (such as cilengitide) were completely unable to increase the dissociation rate. Formation of the non-reversible state occurred very rapidly and did not rely on the time-dependent formation of a high-affinity state of the integrin, or the integrin leg regions. In contrast with RGD-based inhibitors, Ca2+ (but not Mg2+) was able to greatly increase the dissociation rate of integrin-FN complexes, with a half-maximal response at ~0.4 mM Ca2+ for αVβ3-FN. The effect of Ca2+ was overcome by co-addition of Mn2+, but not Mg2+. A stimulatory anti-β1 monoclonal antibody (mAb) abrogated the effect of Ca2+ on α5β1-FN complexes; conversely, a function-blocking mAb mimicked the effect of Ca2+. These results imply that Ca2+ acts allosterically, probably through binding to the adjacent metal-ion-dependent adhesion site (ADMIDAS), and that the α1 helix in the β subunit I domain is the key element affected by allosteric modulators. The data suggest an explanation for the limited clinical efficacy of RGD-based integrin antagonists, and we propose that allosteric antagonists could prove to be of greater therapeutic benefit.

Anadenanthera colubrina (Vell.) Brenan produces steroidal substances that are active against Alternaria alternata (Fr.) Keissler and that may bind to oxysterol-binding proteins

BACKGROUND

In previous studies, the extract from Anadenanthera colubrina was active against Alternaria alternata in vitro and reduced the disease caused by this fungus on Murcott tangor fruits to levels that have been obtained using commercial fungicides. Therefore, the goal of the present work was to isolate and identify the active substances in this extract and identify in silico their protein target in the fungus.

RESULTS

The bioguided fractionation of the methanol extract from the fruits of A. colubrina resulted in the isolation of β-sitosterol and β-sitosteryl linoleate, which had minimal inhibitory concentrations (MICs) of 250 and 500 µg mL(-1) , respectively, against A. alternata. Under the same conditions, the MICs for two commercial fungicides were 1250 and 19 µg mL(-1) . In silico studies showed that these steroidal substances bind well to oxysterol-binding proteins from Saccharomyces cerevisiae.

CONCLUSION

β-Sitosterol and β-sitosteryl linoleate, produced by A. colubrina, are active against A. alternata. In silico studies suggest that these substances may act by binding to oxysterol-binding proteins. Therefore, both substances and these proteins have potential use in the development of new steroidal structures and analogues to control the disease caused by A. alternata.

Structure-based modeling of protein: DNA specificity

Protein:DNA interactions are essential to a range of processes that maintain and express the information encoded in the genome. Structural modeling is an approach that aims to understand these interactions at the physicochemical level. It has been proposed that structural modeling can lead to deeper understanding of the mechanisms of protein:DNA interactions, and that progress in this field can not only help to rationalize the observed specificities of DNA-binding proteins but also to allow researchers to engineer novel DNA site specificities. In this review we discuss recent developments in the structural description of protein:DNA interactions and specificity, as well as the challenges facing the field in the future.

Conformational selection in protein binding and function

Protein binding and function often involves conformational changes. Advanced nuclear magnetic resonance (NMR) experiments indicate that these conformational changes can occur in the absence of ligand molecules (or with bound ligands), and that the ligands may "select" protein conformations for binding (or unbinding). In this review, we argue that this conformational selection requires transition times for ligand binding and unbinding that are small compared to the dwell times of proteins in different conformations, which is plausible for small ligand molecules. Such a separation of timescales leads to a decoupling and temporal ordering of binding/unbinding events and conformational changes. We propose that conformational-selection and induced-change processes (such as induced fit) are two sides of the same coin, because the temporal ordering is reversed in binding and unbinding direction. Conformational-selection processes can be characterized by a conformational excitation that occurs prior to a binding or unbinding event, while induced-change processes exhibit a characteristic conformational relaxation that occurs after a binding or unbinding event. We discuss how the ordering of events can be determined from relaxation rates and effective on- and off-rates determined in mixing experiments, and from the conformational exchange rates measured in advanced NMR or single-molecule fluorescence resonance energy transfer experiments. For larger ligand molecules such as peptides, conformational changes and binding events can be intricately coupled and exhibit aspects of conformational-selection and induced-change processes in both binding and unbinding direction.

Kink turn sRNA folding upon L7Ae binding using molecular dynamics simulations

The kink-turn sRNA motif in archaea, whose combination with protein L7Ae initializes the assembly of small ribonucleoprotein particles (sRNPs), plays a key role in ribosome maturation and the translation process. Although many studies have been reported on this motif, the mechanism of sRNA folding coupled with protein binding is still poorly understood. Here, room and high temperature molecular dynamics (MD) simulations were performed on the complex of 25-nt kink-turn sRNA and L7Ae. The average RMSD values between the bound and corresponding apo structures and Kolmogorov-Smirnov P test analysis indicate that sRNA may follow an induced fit mechanism upon binding with L7Ae, both locally and globally. These conclusions are further supported by high-temperature unfolding kinetic analysis. Principal component analysis (PCA) found both closing and opening motions of the kink-turn sRNA. This might play a key role in the sRNP assembly and methylation catalysis. These combined computational methods can be used to study the specific recognition of other sRNAs and proteins.

Molecular biology: Protein binding cannot subdue a lively RNA

Ribosomes, the cell’s protein-synthesis machines, are assembled from their components in a defined order. It emerges that the first assembly step must overcome dynamic structural rearrangements.

Multiple conformational selection and induced fit events take place in allosteric propagation

The fact that we observe a single conformational selection event during binding does not necessarily mean that only a single conformational selection event takes place, even though this is the common assumption. Here we suggest that conformational selection takes place not once in a given binding/allosteric event, but at every step along the allosteric pathway. This view generalizes conformational selection and makes it applicable also to other allosteric events, such as post-translational modifications (PTMs) and photon absorption. Similar to binding, at each step along a propagation pathway, conformational selection is coupled with induced fit which optimizes the interactions. Thus, as in binding, the allosteric effects induced by PTMs and light relate not only to population shift; but to conformational selection as well. Conformational selection and population shift take place conjointly.