Essential role of conformational selection in ligand binding

RNA's Dynamic Conformational Selection and Entropic Allosteric Mechanism in Controlling Cascade Protein Binding Events

In the TAR RNA of immunodeficiency viruses, an allosteric communication exists between a distant loop and a bulge. The bulge interacts with the TAT protein vital for transactivating viral RNA, while the loop interacts with cyclin-T1, contingent on TAT binding. Through extensive atomistic and free energy simulations, we investigate TAR-TAT binding in nonpathogenic bovine immunodeficiency virus (BIV) and pathogenic human immunodeficiency virus (HIV). Thermodynamic analysis reveals enthalpically driven binding in BIV and entropically favored binding in HIV. The broader global basin in HIV is attributed to binding-induced loop fluctuation, corroborated by nuclear magnetic resonance (NMR), indicating classical entropic allostery onset. While this loop fluctuation affects the TAT binding affinity, it generates a binding-competent conformation that aids subsequent effector (cyclin-T1) binding. This study underscores how two structurally similar apo-RNA scaffolds adopt distinct conformational selection mechanisms to drive enthalpic and entropic allostery, influencing protein affinity in the signaling cascade.

Elucidating Dynamics of Adenylate Kinase from Enzyme Opening to Ligand Release

This study explores ligand-driven conformational changes in adenylate kinase (AK), which is known for its open-to-close conformational transitions upon ligand binding and release. By utilizing string free energy simulations, we determine the free energy profiles for both enzyme opening and ligand release and compare them with profiles from the apoenzyme. Results reveal a three-step ligand release process, which initiates with the opening of the adenosine triphosphate-binding subdomain (ATP lid), followed by ligand release and concomitant opening of the adenosine monophosphate-binding subdomain (AMP lid). The ligands then transition to nonspecific positions before complete dissociation. In these processes, the first step is energetically driven by ATP lid opening, whereas the second step is driven by ATP release. In contrast, the AMP lid opening and its ligand release make minor contributions to the total free energy for enzyme opening. Regarding the ligand binding mechanism, our results suggest that AMP lid closure occurs via an induced-fit mechanism triggered by AMP binding, whereas ATP lid closure follows conformational selection. This difference in the closure mechanisms provides an explanation with implications for the debate on ligand-driven conformational changes of AK. Additionally, we determine an X-ray structure of an AK variant that exhibits significant rearrangements in the stacking of catalytic arginines, explaining its reduced catalytic activity. In the context of apoenzyme opening, the sequence of events is different. Here, the AMP lid opens first while the ATP lid remains closed, and the free energy associated with ATP lid opening varies with orientation, aligning with the reported AK opening and closing rate heterogeneity. Finally, this study, in conjunction with our previous research, provides a comprehensive view of the intricate interplay between various structural elements, ligands, and catalytic residues that collectively contribute to the robust catalytic power of the enzyme.

A substrate-induced gating mechanism is conserved among Gram-positive IgA1 metalloproteases

The mucosal adaptive immune response is dependent on the production of IgA antibodies and particularly IgA1, yet opportunistic bacteria have evolved mechanisms to specifically block this response by producing IgA1 proteases (IgA1Ps). Our lab was the first to describe the structures of a metal-dependent IgA1P (metallo-IgA1P) produced from Gram-positive Streptococcus pneumoniae both in the absence and presence of its IgA1 substrate through cryo-EM single particle reconstructions. This prior study revealed an active-site gating mechanism reliant on substrate-induced conformational changes to the enzyme that begged the question of whether such a mechanism is conserved among the wider Gram-positive metallo-IgA1P subfamily of virulence factors. Here, we used cryo-EM to characterize the metallo-IgA1P of a more distantly related family member from Gemella haemolysans, an emerging opportunistic pathogen implicated in meningitis, endocarditis, and more recently bacteremia in the elderly. While the substrate-free structures of these two metallo-IgA1Ps exhibit differences in the relative starting positions of the domain responsible for gating substrate, the enzymes have similar domain orientations when bound to IgA1. Together with biochemical studies that indicate these metallo-IgA1Ps have similar binding affinities and activities, these data indicate that metallo-IgA1P binding requires the specific IgA1 substrate to open the enzymes for access to their active site and thus, largely conform to an "induced fit" model.

Getting Deeper into the Molecular Events of Heme Binding Mechanisms: A Comparative Multi-level Computational Study of HasAsm and HasAyp Hemophores

Many biological systems obtain their activity by the inclusion of metalloporphyrins into one or several binding pockets. However, decoding the molecular mechanism under which these compounds bind to their receptors is something that has not been widely explored and is a field with open questions. In the present work, we apply computational techniques to unravel and compare the mechanisms of two heme-binding systems, concretely the HasA hemophores from Gram negative bacteria Serratiamarcescens (HasAsm) and Yersinia pestis (HasAyp). Despite the high sequence identity between both systems, the comparison between the X-ray structures of their apo and holo forms suggests different heme-binding mechanisms. HasAyp has extremely similar structures for heme-free and heme-bound forms, while HasAsm presents a very large displacement of a loop that ultimately leads to an additional coordination to the metal with respect to HasAyp. We combined Gaussian accelerated molecular dynamics simulations (GaMDs) in explicit solvent and protein–ligand docking optimized for metalloligands. GaMDs were first carried out on heme-free forms of both hemophores. Then, protein–ligand dockings of the heme were performed on cluster representatives of these simulations and the best poses were then subjected to a new series of GaMDs. A series of analyses reveal the following: (1) HasAyp has a conformational landscape extremely similar between heme-bound and unbound states with no to limited impact on the binding of the cofactor, (2) HasAsm presents as a slightly broader conformational landscape in its apo state but can only visit conformations similar to the X-ray of the holo form when the heme has been bound. Such behavior results from a complex cascade of changes in interactions that spread from the heme-binding pocket to the flexible loop previously mentioned. This study sheds light on the diversity of molecular mechanisms of heme-binding and discusses the weight between the pre-organization of the receptor as well as the induced motions resulting in association.

Heme-containing enzymes and proteins are important for many biological and biotechnological processes. However, very little is known about heme-binding mechanisms. To shed light on this, we report a multi-level approach combining Gaussian accelerated molecular dynamics and protein−ligand dockings optimized for metallic moieties. The protocol unveils the difference in heme recruitment between HasAsm and HasAyp hemophores and shows its possible applicability to other heme-binding proteins.

Large-Scale Ligand Perturbations of the Protein Conformational Landscape Reveal State-Specific Interaction Hotspots

Protein flexibility is important for ligand binding but often ignored in drug design. Considering proteins as ensembles rather than static snapshots creates opportunities to target dynamic proteins that lack FDA-approved drugs, such as the human chaperone, heat shock protein 90 (Hsp90). Hsp90α accommodates ligands with a dynamic lid domain, yet no comprehensive analysis relating lid conformations to ligand properties is available. To date, ∼300 ligand-bound Hsp90α crystal structures are deposited in the Protein Data Bank, which enables us to consider ligand binding as a perturbation of the protein conformational landscape. By estimating binding site volumes, we classified structures into distinct major and minor lid conformations. Supported by retrospective docking, each conformation creates unique hotspots that bind chemically distinguishable ligands. Clustering revealed insightful exceptions and the impact of crystal packing. Overall, Hsp90α’s plasticity provides a cautionary tale of overinterpreting individual crystal structures and motivates an ensemble-based view of drug design.

Conformational Selection Governs Carrier Domain Positioning in Staphylococcus aureus Pyruvate Carboxylase

Biotin-dependent enzymes employ a carrier domain to efficiently transport reaction intermediates between distant active sites. The translocation of this carrier domain is critical to the interpretation of kinetic and structural studies, but there have been few direct attempts to investigate the dynamic interplay between ligand binding and carrier domain positioning in biotin-dependent enzymes. Pyruvate carboxylase (PC) catalyzes the MgATP-dependent carboxylation of pyruvate where the biotinylated carrier domain must translocate ∼70 Å from the biotin carboxylase domain to the carboxyltransferase domain. Many prior studies have assumed that carrier domain movement is governed by ligand-induced conformational changes, but the mechanism underlying this movement has not been confirmed. Here, we have developed a system to directly observe PC carrier domain positioning in both the presence and absence of ligands, independent of catalytic turnover. We have incorporated a cross-linking trap that reports on the interdomain conformation of the carrier domain when it is positioned in proximity to a neighboring carboxyltransferase domain. Cross-linking was monitored by gel electrophoresis, inactivation kinetics, and intrinsic tryptophan fluorescence. We demonstrate that the carrier domain positioning equilibrium is sensitive to substrate analogues and the allosteric activator acetyl-CoA. Notably, saturating concentrations of biotin carboxylase ligands do not prevent carrier domain trapping proximal to the neighboring carboxyltransferase domain, demonstrating that carrier domain positioning is governed by conformational selection. This model of carrier domain translocation in PC can be applied to other multi-domain enzymes that employ large-scale domain motions to transfer intermediates during catalysis.

A litmus test for classifying recognition mechanisms of transiently binding proteins

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

The authors provide a litmus test for the recognition mechanism of transiently binding proteins based on nuclear magnetic resonance and find a conformational selection binding mechanism through concentration-dependent kinetics of ubiquitin and SH3.

Kinetic solvent viscosity effects uncover an internal isomerization of the enzyme-substrate complex in Pseudomonas aeruginosa PAO1 NADH:Quinone oxidoreductase

NAD(P)H:quinone oxidoreductases (NQOs) play an essential protective role as antioxidants in the detoxification of quinones in both Prokaryotes and Eukaryotes. NQO from Pseudomonas aeruginosa PAO1 uses FMN to catalyze the two-electron reduction of various quinones with NADH. In this study, steady-state kinetics, kinetic solvent viscosity effects, and rapid reaction kinetics were used to determine which kinetic steps control the overall turnover of the enzyme with benzoquinone or juglone. The rate constant for flavin reduction (k_red) at pH 6.0 was 12.9 ± 0.3 s^-1, and the K_d for NADH was at least an order of magnitude lower than 90 μM. With benzoquinone, the k_cat value was 11.7 ± 0.3 s^-1, consistent with flavin reduction being almost entirely rate-limiting for overall turnover. With juglone, a k_cat value of 10.0 ± 0.5 s^-1 was recorded. The normalized plot of the relative solvent viscosity effects on the k_cat values established that hydride transfer from NADH to the FMN and quinol product release, with a calculated rate constant (k_P-rel) of 52 s^-1, are partially rate-limiting for the overall turnover of NQO. Kinetic solvent viscosity effects with glucose or sucrose revealed a hyperbolic dependence on the k_cat and k_cat/K_m values with benzoquinone or juglone, respectively, consistent with the presence of a solvent-sensitive internal isomerization of the enzyme-substrate complex (ES). The data demonstrate opposing effects of benzoquinone and juglone on the equilibrium of the NQO ES isomerization with glucose or sucrose. Thus, our study demonstrates how quinol substrate properties alter the equilibrium of NQO ES isomerization.

Insights into Membrane Curvature Sensing and Membrane Remodeling by Intrinsically Disordered Proteins and Protein Regions

Cellular membranes are highly dynamic in shape. They can rapidly and precisely regulate their shape to perform various cellular functions. The protein’s ability to sense membrane curvature is essential in various biological events such as cell signaling and membrane trafficking. As they are bound, these curvature-sensing proteins may also change the local membrane shape by one or more curvature driving mechanisms. Established curvature-sensing/driving mechanisms rely on proteins with specific structural features such as amphipathic helices and intrinsically curved shapes. However, the recent discovery and characterization of many proteins have shattered the protein structure–function paradigm, believing that the protein functions require a unique structural feature. Typically, such structure-independent functions are carried either entirely by intrinsically disordered proteins or hybrid proteins containing disordered regions and structured domains. It is becoming more apparent that disordered proteins and regions can be potent sensors/inducers of membrane curvatures. In this article, we outline the basic features of disordered proteins and regions, the motifs in such proteins that encode the function, membrane remodeling by disordered proteins and regions, and assays that may be employed to investigate curvature sensing and generation by ordered/disordered proteins.

Antigen binding by conformational selection in near-germline antibodies

Conformational flexibility in antibody-combining sites has been hypothesized to facilitate polyspecificity toward multiple unique epitopes and enable the limited germline repertoire to match an overwhelming diversity of potential antigens; however, elucidating the mechanisms of antigen recognition by flexible antibodies has been understandably challenging. Here, multiple liganded and unliganded crystal structures of the near-germline anticarbohydrate antibodies S25–2 and S25–39 are reported, which reveal an unprecedented diversity of complementarity-determining region H3 conformations in apparent equilibrium. These structures demonstrate that at least some germline or near-germline antibodies are flexible entities sensitive to their chemical environments, with conformational selection available as an evolved mechanism that preserves the inherited ability to recognize common pathogens while remaining adaptable to new threats.

Roles of cytochrome P450 enzymes in pharmacology and toxicology: Past, present, and future

The development of the cytochrome P450 (P450) field has been remarkable in the areas of pharmacology and toxicology, particularly in drug development. Today it is possible to use the knowledge base and relatively straightforward assays to make intelligent predictions about drug disposition prior to human dosing. Much is known about the structures, regulation, chemistry of catalysis, and the substrate and inhibitor specificity of human P450s. Many aspects of drug-drug interactions and side effects can be understood in terms of P450s. This knowledge has also been useful in pharmacy practice, as well as in the pharmaceutical industry and medical practice. However, there are still basic and practical questions to address regarding P450s and their roles in pharmacology and toxicology. Another aspect is the discovery of drugs that inhibit P450 to treat diseases.

The active site region plays a critical role in Na+ binding to thrombin

The catalytic activity of thrombin and other enzymes of the blood coagulation and complement cascades is enhanced significantly by binding of Na+ to a site >15 Å away from the catalytic residue S195, buried within the 180 and 220 loops that also contribute to the primary specificity of the enzyme. Rapid kinetics support a binding mechanism of conformational selection where the Na+-binding site is in equilibrium between open (N) and closed (N∗) forms and the cation binds selectively to the N form. Allosteric transduction of this binding step produces enhanced catalytic activity. Molecular details on how Na+ gains access to this site and communicates allosterically with the active site remain poorly defined. In this study, we show that the rate of the N∗→N transition is strongly correlated with the analogous E∗→E transition that governs the interaction of synthetic and physiologic substrates with the active site. This correlation supports the active site as the likely point of entry for Na+ to its binding site. Mutagenesis and structural data rule out an alternative path through the pore defined by the 180 and 220 loops. We suggest that the active site communicates allosterically with the Na+ site through a network of H-bonded water molecules that embeds the primary specificity pocket. Perturbation of the mobility of S195 and its H-bonding capabilities alters interaction with this network and influences the kinetics of Na+ binding and allosteric transduction. These findings have general mechanistic relevance for Na+-activated proteases and allosteric enzymes.

Exploring the Structures and Functions of Macromolecular SLX4-Nuclease Complexes in Genome Stability

All organisms depend on the ability of cells to accurately duplicate and segregate DNA into progeny. However, DNA is frequently damaged by factors in the environment and from within cells. One of the most dangerous lesions is a DNA double-strand break. Unrepaired breaks are a major driving force for genome instability. Cells contain sophisticated DNA repair networks to counteract the harmful effects of genotoxic agents, thus safeguarding genome integrity. Homologous recombination is a high-fidelity, template-dependent DNA repair pathway essential for the accurate repair of DNA nicks, gaps and double-strand breaks. Accurate homologous recombination depends on the ability of cells to remove branched DNA structures that form during repair, which is achieved through the opposing actions of helicases and structure-selective endonucleases. This review focuses on a structure-selective endonuclease called SLX1-SLX4 and the macromolecular endonuclease complexes that assemble on the SLX4 scaffold. First, we discuss recent developments that illuminate the structure and biochemical properties of this somewhat atypical structure-selective endonuclease. We then summarize the multifaceted roles that are fulfilled by human SLX1-SLX4 and its associated endonucleases in homologous recombination and genome stability. Finally, we discuss recent work on SLX4-binding proteins that may represent integral components of these macromolecular nuclease complexes, emphasizing the structure and function of a protein called SLX4IP.

Induced fit versus conformational selection: From rate constants to fluxes… and back to rate constants

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.

Nonspecific Binding-Fundamental Concepts and Consequences for Biosensing Applications

Nature achieves differentiation of specific and nonspecific binding in molecular interactions through precise control of biomolecules in space and time. Artificial systems such as biosensors that rely on distinguishing specific molecular binding events in a sea of nonspecific interactions have struggled to overcome this issue. Despite the numerous technological advancements in biosensor technologies, nonspecific binding has remained a critical bottleneck due to the lack of a fundamental understanding of the phenomenon. To date, the identity, cause, and influence of nonspecific binding remain topics of debate within the scientific community. In this review, we discuss the evolution of the concept of nonspecific binding over the past five decades based upon the thermodynamic, intermolecular, and structural perspectives to provide classification frameworks for biomolecular interactions. Further, we introduce various theoretical models that predict the expected behavior of biosensors in physiologically relevant environments to calculate the theoretical detection limit and to optimize sensor performance. We conclude by discussing existing practical approaches to tackle the nonspecific binding challenge in vitro for biosensing platforms and how we can both address and harness nonspecific interactions for in vivo systems.

Distal Regions Regulate Dihydrofolate Reductase-Ligand Interactions

Protein motions play a fundamental role in enzyme catalysis and ligand binding. The relationship between protein motion and function has been extensively investigated in the model enzyme dihydrofolate reductase (DHFR). DHFR is an essential enzyme that catalyzes the reduction of dihydrofolate to tetrahydrofolate. Numerous experimental and computational methods have been used to probe the motions of DHFR through the catalytic cycle and to investigate the effect of distal mutations on DHFR motions and ligand binding. These experimental investigations have pushed forward the study of protein motions and their role in protein-ligand interactions. The introduction of mutations distal to the active site has been shown to have profound effects on ligand binding, hydride transfer rates and catalytic efficacy and these changes are captured by enzyme kinetics measurements. Distal mutations have been shown to exert their effects through a network of correlated amino acids and these effects have been investigated by NMR, protein dynamics, and analysis of coupled amino acids. The experimental methods and the findings that are reviewed here have broad implications for our understanding of enzyme mechanisms, ligand binding and for the future design and discovery of enzyme inhibitors.

Peptide-based inhibitors of protein-protein interactions: biophysical, structural and cellular consequences of introducing a constraint

Protein-protein interactions (PPIs) are implicated in the majority of cellular processes by enabling and regulating the function of individual proteins. Thus, PPIs represent high-value, but challenging targets for therapeutic intervention. The development of constrained peptides represents an emerging strategy to generate peptide-based PPI inhibitors, typically mediated by α-helices. The approach can confer significant benefits including enhanced affinity, stability and cellular penetration and is ingrained in the premise that pre-organization simultaneously pays the entropic cost of binding, prevents a peptide from adopting a protease compliant β-strand conformation and shields the hydrophilic amides from the hydrophobic membrane. This conceptual blueprint for the empirical design of peptide-based PPI inhibitors is an exciting and potentially lucrative way to effect successful PPI inhibitor drug-discovery. However, a plethora of more subtle effects may arise from the introduction of a constraint that include changes to binding dynamics, the mode of recognition and molecular properties. In this review, we summarise the influence of inserting constraints on biophysical, conformational, structural and cellular behaviour across a range of constraining chemistries and targets, to highlight the tremendous success that has been achieved with constrained peptides alongside emerging design opportunities and challenges.

Host assisted molecular recognition by human serum albumin: Study of molecular recognition controlled protein/drug mimic binding in a microfluidic channel

Human serum albumin (HSA) plays a pivotal role in drug release from its delivery vehicles such as cyclodextrins (CDs) by binding to the drugs. Here molecular recognition and binding of a drug mimic (CD1) to HSA have been explored in a microfluidic channel when CD1 is encapsulated in β-cyclodextrin (βCD) and heptakis(2,3,6-tri-O-methyl)-β-cyclodextrin (TRIMEB), respectively, to investigate whether change of the host vehicle modulate the rate of drug binding to the serum protein. Molecular recognition of βCD encapsulated CD1 by HSA occurs by the conformational selection fit mechanism leading to rapid binding of CD1 to HSA (k₁ ~ 700 s^-11) when the βCD/CD1 complex interacts with HSA. In contrary, HSA recognizes CD1 encapsulated in TRIMEB by an induced fit mechanism leading to a significantly slower binding rate (k₁ ~ 20.8 s^-1) of the drug mimic to the protein. Thus molecular recognition controls the rate of HSA binding by CD1 which in turn modulates the rate of delivery of the drug mimic from its macrocyclic hosts. The remarkable change in the molecular recognition pathway of CD1 by HSA, upon change of the host from βCD to TRIMEB, originates from significantly different conformational flexibility of the host/drug mimic complexes.

Ligand binding free-energy calculations with funnel metadynamics

The accurate resolution of the binding mechanism of a ligand to its molecular target is fundamental to develop a successful drug design campaign. Free-energy calculations, which provide the energy value of the ligand-protein binding complex, are essential for resolving the binding mode of the ligand. The accuracy of free-energy calculation methods is counteracted by their poor user-friendliness, which hampers their broad application. Here we present the Funnel-Metadynamics Advanced Protocol (FMAP), which is a flexible and user-friendly graphical user interface (GUI)-based protocol to perform funnel metadynamics, a binding free-energy method that employs a funnel-shape restraint potential to reveal the ligand binding mode and accurately calculate the absolute ligand-protein binding free energy. FMAP guides the user through all phases of the free-energy calculation process, from preparation of the input files, to production simulation, to analysis of the results. FMAP delivers the ligand binding mode and the absolute protein-ligand binding free energy as outputs. Alternative binding modes and the role of waters are also elucidated, providing a detailed description of the ligand binding mechanism. The entire protocol on the paradigmatic system benzamidine-trypsin, composed of ~105 k atoms, took ~2.8 d using the Cray XC50 piz Daint cluster at the Swiss National Supercomputing Centre.

Energetics and kinetics of substrate analog-coupled staphylococcal nuclease folding revealed by a statistical mechanical approach

Protein conformational changes associated with ligand binding, especially those involving intrinsically disordered proteins, are mediated by tightly coupled intra- and intermolecular events. Such reactions are often discussed in terms of two limiting kinetic mechanisms, conformational selection (CS), where folding precedes binding, and induced fit (IF), where binding precedes folding. It has been shown that coupled folding/binding reactions can proceed along both CS and IF pathways with the flux ratio depending on conditions such as ligand concentration. However, the structural and energetic basis of such complex reactions remains poorly understood. Therefore, we used experimental, theoretical, and computational approaches to explore structural and energetic aspects of the coupled-folding/binding reaction of staphylococcal nuclease in the presence of the substrate analog adenosine-3',5'-diphosphate. Optically monitored equilibrium and kinetic data, combined with a statistical mechanical model, gave deeper insight into the relative importance of specific and Coulombic protein-ligand interactions in governing the reaction mechanism. We also investigated structural aspects of the reaction at the residue level using NMR and all-atom replica-permutation molecular dynamics simulations. Both approaches yielded clear evidence for accumulation of a transient protein-ligand encounter complex early in the reaction under IF-dominant conditions. Quantitative analysis of the equilibrium/kinetic folding revealed that the ligand-dependent CS-to-IF shift resulted from stabilization of the compact transition state primarily by weakly ligand-dependent Coulombic interactions with smaller contributions from specific binding energies. At a more macroscopic level, the CS-to-IF shift was represented as a displacement of the reaction "route" on the free energy surface, which was consistent with a flux analysis.

Ligand-bound glutamine binding protein assumes multiple metastable binding sites with different binding affinities

Protein dynamics plays key roles in ligand binding. However, the microscopic description of conformational dynamics-coupled ligand binding remains a challenge. In this study, we integrate molecular dynamics simulations, Markov state model (MSM) analysis and experimental methods to characterize the conformational dynamics of ligand-bound glutamine binding protein (GlnBP). We show that ligand-bound GlnBP has high conformational flexibility and additional metastable binding sites, presenting a more complex energy landscape than the scenario in the absence of ligand. The diverse conformations of GlnBP demonstrate different binding affinities and entail complex transition kinetics, implicating a concerted ligand binding mechanism. Single molecule fluorescence resonance energy transfer measurements and mutagenesis experiments are performed to validate our MSM-derived structure ensemble as well as the binding mechanism. Collectively, our study provides deeper insights into the protein dynamics-coupled ligand binding, revealing an intricate regulatory network underlying the apparent binding affinity.

Zhang, Wu, Feng et al. show that ligand-bound glutamine binding protein assumes multiple metastable binding sites, presenting a more dynamic energy landscape than its ligand-free form. This study provides insights into the ligand-binding mechanisms coupled with protein dynamics that underly the apparent binding affinity.

Conformational selection is present in ligand binding to cytochrome P450 19A1 lipoprotein nanodiscs

Cytochromes P450 (CYPs) display remarkable plasticity in their ability to bind substrates and catalyze a broad array of chemical reactions. Herein we evaluate binding of androstenedione, testosterone, and 7-hydroxyflavone to CYP19A1, also known as aromatase, in phospholipid nanodiscs by stopped-flow UV-vis spectroscopy. Exponential fitting of the kinetic traces supports the possibility of a multi-step binding mechanism. Subsequent global fitting of the data to the solutions of the coupled differential equations describing the fundamental mechanisms of induced fit and conformational selection, consistently support presence of the latter. To our knowledge, this is the first discrimination of conformational selection from induced fit for a mono-disperse CYP in a native-like membrane environment. In addition, 7-hydroxyflavone binds to CYP19A1 nanodiscs with comparable affinity to the substrates and induces an unusual spectral response likely attributable to hydrogen bonding to, rather than displacement of the heme-coordinated water molecule.

19F NMR reveals the conformational properties of free thrombin and its zymogen precursor prethrombin-2

The conformational properties of trypsin-like proteases and their zymogen forms remain controversial because of a lack of sufficient information on their free forms. Specifically, it is unclear whether the free protease is zymogen-like and shifts to its mature form upon a ligand-induced fit or exists in multiple conformations in equilibrium from which the ligand selects the optimal fit via conformational selection. Here we report the results of ¹⁹F NMR measurements that reveal the conformational properties of a protease and its zymogen precursor in the free form. Using the trypsin-like, clotting protease thrombin as a relevant model system, we show that its conformation is quite different from that of its direct zymogen precursor prethrombin-2 and more similar to that of its fully active Na⁺-bound form. The results cast doubts on recent hypotheses that free thrombin is zymogen-like and transitions to protease-like forms upon ligand binding. Rather, they validate the scenario emerged from previous findings of X-ray crystallography and rapid kinetics supporting a pre-existing equilibrium between open (E) and closed (E*) forms of the active site. In this scenario, prethrombin-2 is more dynamic and exists predominantly in the E* form, whereas thrombin is more rigid and exists predominantly in the E form. Ligand binding to thrombin takes place exclusively in the E form without significant changes in the overall conformation. In summary, these results disclose the structural architecture of the free forms of thrombin and prethrombin-2, consistent with an E*-E equilibrium and providing no evidence that free thrombin is zymogen-like.

Mechanisms of promiscuity among drug metabolizing enzymes and drug transporters

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.

Antibody Structure and Function: The Basis for Engineering Therapeutics

Antibodies and antibody-derived macromolecules have established themselves as the mainstay in protein-based therapeutic molecules (biologics). Our knowledge of the structure–function relationships of antibodies provides a platform for protein engineering that has been exploited to generate a wide range of biologics for a host of therapeutic indications. In this review, our basic understanding of the antibody structure is described along with how that knowledge has leveraged the engineering of antibody and antibody-related therapeutics having the appropriate antigen affinity, effector function, and biophysical properties. The platforms examined include the development of antibodies, antibody fragments, bispecific antibody, and antibody fusion products, whose efficacy and manufacturability can be improved via humanization, affinity modulation, and stability enhancement. We also review the design and selection of binding arms, and avidity modulation. Different strategies of preparing bispecific and multispecific molecules for an array of therapeutic applications are included.

Membrane interactions of intrinsically disordered proteins: The example of alpha-synuclein

Peripheral membrane proteins associate reversibly with biological membranes that, compared to protein binding partners, are structurally labile and devoid of specific binding pockets. Membranes in different subcellular compartments vary primarily in their chemical composition and physical properties, and recognition of these features is therefore critical for allowing such proteins to engage their proper membrane targets. Intrinsically disordered proteins (IDPs) are well-suited to accomplish this task using highly specific and low- to moderate-affinity interactions governed by recognition principles that are both similar to and different from those that mediate the membrane interactions of rigid proteins. IDPs have also evolved multiple mechanisms to regulate membrane (and other) interactions and achieve their impressive functional diversity. Moreover, IDP-membrane interactions may have a kinetic advantage in fast processes requiring rapid control of such interactions, such as synaptic transmission or signaling. Herein we review the biophysics, regulation and functional implications of IDP-membrane interactions and include a brief overview of some of the methods that can be used to study such interactions. At each step, we use the example of alpha-synuclein, a protein involved in the pathogenesis of Parkinson's disease and one of the best characterized membrane-binding IDP, to illustrate some of the principles discussed.

Residues W215, E217 and E192 control the allosteric E*-E equilibrium of thrombin

A pre-existing, allosteric equilibrium between closed (E*) and open (E) conformations of the active site influences the level of activity in the trypsin fold and defines ligand binding according to the mechanism of conformational selection. Using the clotting protease thrombin as a model system, we investigate the molecular determinants of the E*-E equilibrium through rapid kinetics and X-ray structural biology. The equilibrium is controlled by three residues positioned around the active site. W215 on the 215-217 segment defining the west wall of the active site controls the rate of transition from E to E* through hydrophobic interaction with F227. E192 on the opposite 190-193 segment defining the east wall of the active site controls the rate of transition from E* to E through electrostatic repulsion of E217. The side chain of E217 acts as a lever that moves the entire 215-217 segment in the E*-E equilibrium. Removal of this side chain converts binding to the active site to a simple lock-and-key mechanism and freezes the conformation in a state intermediate between E* and E. These findings reveal a simple framework to understand the molecular basis of a key allosteric property of the trypsin fold.

It takes two to tango: The dance of the permease

The lactose permease (LacY) of Escherichia coli is the prototype of the major facilitator superfamily, one of the largest families of membrane transport proteins. Structurally, two pseudo-symmetrical six-helix bundles surround a large internal aqueous cavity. Single binding sites for galactoside and H+ are positioned at the approximate center of LacY halfway through the membrane at the apex of the internal cavity. These features enable LacY to function by an alternating-access mechanism that can catalyze galactoside/H+ symport in either direction across the cytoplasmic membrane. The H+-binding site is fully protonated under physiological conditions, and subsequent sugar binding causes transition of the ternary complex to an occluded intermediate that can open to either side of the membrane. We review the structural and functional evidence that has provided new insight into the mechanism by which LacY achieves active transport against a concentration gradient.

Human cytochrome P450 enzymes bind drugs and other substrates mainly through conformational-selection modes

Cytochrome P450 (P450) enzymes are major catalysts involved in the oxidations of most drugs, steroids, carcinogens, fat-soluble vitamins, and natural products. The binding of substrates to some of the 57 human P450s and other mammalian P450s is more complex than a two-state system and has been proposed to involve mechanisms such as multiple ligand occupancy, induced-fit, and conformational-selection. Here, we used kinetic analysis of binding with multiple concentrations of substrates and computational modeling of these data to discern possible binding modes of several human P450s. We observed that P450 2D6 binds its ligand rolapitant in a mechanism involving conformational-selection. P450 4A11 bound the substrate lauric acid via conformational-selection, as did P450 2C8 with palmitic acid. Binding of the steroid progesterone to P450 21A2 was also best described by a conformational-selection model. Hexyl isonicotinate binding to P450 2E1 could be described by either a conformational-selection or an induced-fit model. Simulation of the binding of the ligands midazolam, bromocriptine, testosterone, and ketoconazole to P450 3A4 was consistent with an induced-fit or a conformational-selection model, but the concentration dependence of binding rates for varying both P450 3A4 and midazolam concentrations revealed discordance in the parameters, indicative of conformational-selection. Binding of the P450s 2C8, 2D6, 3A4, 4A11, and 21A2 was best described by conformational-selection, and P450 2E1 appeared to fit either mode. These findings highlight the complexity of human P450-substrate interactions and that conformational-selection is a dominant feature of many of these interactions.

General mathematical formula for near equilibrium relaxation kinetics of basic enzyme reactions and its applications to find conformational selection steps

A general mathematical formula of basic enzyme reactions was derived with nearly no dependence on conditions nor assumptions on relaxation kinetic processes near equilibrium in a simple single-substrate-single-product enzyme reaction. The new formula gives precise relationships between the rate constants of the elementary reaction steps and the apparent relaxation rate constant, rather than the initial velocity that is generally used to determine enzymatic parameters according to the Michaelis-Menten theory. The present formula is shown to be complementary to the Michaelis-Menten formulae in a sense that the initial velocity and the relaxation rate constant data together could determine the enzyme-substrate dissociation constant K_ES, which has been usually conditionally approximated by the Michaelis constant K_M within the framework of the Michaelis-Menten formulae. We also describe relaxation kinetics of enzyme reactions that include the conformational selection processes, in which only one enzymatic conformer among a conformational ensemble can bind with either the substrate or product. The present mathematical approaches, together with numerical computation analyses, suggested that the presence of conformational selection steps in enzymatic reactions can be experimentally detected simply by enzymatic assays with catalytic amounts of enzyme.

Structural Explorations of NCp7-Nucleic Acid Complexes Give Keys to Decipher the Binding Process

A comprehensive view of all the structural aspects related to NCp7 is essential to understand how this protein, crucial in many steps of the HIV-1 cycle, binds and anneals nucleic acids (NAs), mainly thanks to two zinc fingers, ZF1 and ZF2. Here, we inspected the structural properties of the available experimental models of NCp7 bound to either DNA or RNA molecules, or free of ligand. Our analyses included the characterization of the relative positioning of ZF1 and ZF2, accessibility measurements and the exhaustive, quantitative mapping of the contacts between amino acids and nucleotides by a recent tessellation method, VLDM. This approach unveiled the intimate connection between NA binding process and the conformations explored by the free protein. It also provided new insights into the functional specializations of ZF1 and ZF2. The larger accessibility of ZF2 in free NCp7 and the consistency of the ZF2/NA interface in different models and conditions give ZF2 the lead of the binding process. ZF1 contributes to stabilize the complexes through various organizations of the ZF1/NA interface. This work outcome is a global binding scheme of NCp7 to DNA and RNA, and an example of how protein-NA complexes are stabilized.

Guest-mediated chirality transfer in the host–guest complexes of an atropisomeric perylene bisimide cyclophane host

Chirality transfer upon preferential binding of homochiral guests to one stereoisomer of a conformationally equilibrated atropisomeric cyclophane is reported.

Y-box proteins combine versatile cold shock domains and arginine-rich motifs (ARMs) for pleiotropic functions in RNA biology

Y-box proteins are single-strand DNA- and RNA-binding proteins distinguished by a conserved cold shock domain (CSD) and a variable C-terminal domain organized into alternating short modules rich in basic or acidic amino acids. A huge literature depicts Y-box proteins as highly abundant, staggeringly versatile proteins that interact with all mRNAs and function in most forms of mRNA-specific regulation. The mechanisms by which Y-box proteins recognize mRNAs are unclear, because their CSDs bind a jumble of diverse elements, and the basic modules in the C-terminal domain are considered to bind nonspecifically to phosphates in the RNA backbone. A survey of vertebrate Y-box proteins clarifies the confusing names for Y-box proteins, their domains, and RNA-binding motifs, and identifies several novel conserved sequences: first, the CSD is flanked by linkers that extend its binding surface or regulate co-operative binding of the CSD and N-terminal and C-terminal domains to proteins and RNA. Second, the basic modules in the C-terminal domain are bona fide arginine-rich motifs (ARMs), because arginine is the predominant amino acid and comprises 99% of basic residues. Third, conserved differences in AA (amino acid) sequences between isoforms probably affect RNA-binding specificity. C-terminal ARMs connect with many studies, demonstrating that ARMs avidly bind sites containing specific RNA structures. ARMs crystallize insights into the under-appreciated contributions of the C-terminal domain to site-specific binding by Y-box proteins and difficulties in identifying site-specific binding by the C-terminal domain. Validated structural biology techniques are available to elucidate the mechanisms by which YBXprot (Y-box element-binding protein) CSDs and ARMs identify targets.

Mathematical and Structural Characterization of Strong Nonadditive Structure-Activity Relationship Caused by Protein Conformational Changes

In medicinal chemistry, accurate prediction of additivity-based structure-activity relationship (SAR) analysis rests on three assumptions: (1) a consistent binding pose of the central scaffold, (2) no interaction between the R group substituents, and (3) a relatively rigid binding pocket in which the R group substituents act independently. Previously, examples of nonadditive SAR have been documented in systems that deviate from the first two assumptions. Local protein structural change upon ligand binding, through induced fit or conformational selection, although a well-known phenomenon that invalidates the third assumption, has not been linked to nonadditive SAR conclusively. Here, for the first time, we present clear structural evidence that the formation of a hydrophobic pocket upon ligand binding in PDE2 catalytic site reduces the size of another distinct subpocket and contributes to strong nonadditive SAR between two otherwise distant R groups.

Binding of l-Argininamide to a DNA Aptamer: A Volumetric Study

We use a combination of volumetric and spectroscopic techniques to characterize the binding of l-argininamide to its aptamer, the 24-base DNA hairpin 5'-d(GATCGAAACGTAGCGCCTTCGATC)-3'. The binding causes increases in volume, Δ V, and adiabatic compressibility, Δ K_S, of 12 ± 7 cm³ mol^-1 bar and (73 ± 8) × 10^-4 cm³ mol^-1 bar^-1, respectively. These volumetric results combined with structural data reveal that the binding is accompanied by release of 73 ± 27 waters from the hydration shells of the interacting molecules to the bulk. We use the estimated change in hydration to estimate the hydration, Δ S_hyd, and configurational, Δ S_conf, contributions to the binding entropy. The large and unfavorable change in configurational entropy, Δ S_conf, is nearly compensated by a favorable change in the hydration contribution, Δ S_hyd.

Magnesium Activates Microsecond Dynamics to Regulate Integrin-Collagen Recognition

Integrin receptors bind collagen via metal-mediated interactions that are modulated by magnesium (Mg²⁺) levels in the extracellular matrix. Nuclear magnetic resonance-based relaxation experiments, isothermal titration calorimetry, and adhesion assays reveal that Mg²⁺ functions as both a structural anchor and dynamic switch of the α₁β₁ integrin I domain (α₁I). Specifically, Mg²⁺ binding activates micro- to millisecond timescale motions of residues distal to the binding site, particularly those surrounding the salt bridge at helix 7 and near the metal ion-dependent adhesion site. Mutagenesis of these residues impacts α₁I functional activity, thereby suggesting that Mg-bound α₁I dynamics are important for collagen binding and consequent allosteric rearrangement of the low-affinity closed to high-affinity open conformation. We propose a multistep recognition mechanism for α₁I-Mg-collagen interactions involving both conformational selection and induced-fit processes. Our findings unravel the multifaceted role of Mg²⁺ in integrin-collagen recognition and assist in elucidating the molecular mechanisms by which metals regulate protein-protein interactions.

Molecular mechanisms of AcrB-mediated multidrug export

The over-expression of multidrug efflux pumps belonging to the Resistance-Nodulation-Division (RND) superfamily is one of the main causes of multidrug-resistance (MDR) in Gram-negative pathogenic bacteria. AcrB is the most thoroughly studied RND transporter and has functioned as a model for our understanding of efflux-mediated MDR. This multidrug-exporter can recognize and transport a wide range of structurally unrelated compounds (including antibiotics, dyes, bile salts and detergents), while it shows a strict inhibitor specificity. Here we discuss our current knowledge of AcrB, and include recent advances, regarding its structure, mechanism of drug transport, substrate recognition, different intramolecular entry pathways and the drug export driven by remote conformational coupling.

Self-organization, entropy and allostery

Allostery is a fundamental regulatory mechanism in biology. Although generally accepted that it is a dynamics-driven process, the exact molecular mechanism of allosteric signal transmission is hotly debated. We argue that allostery is as a part of a bigger picture that also includes fractal-like properties of protein interior, hierarchical protein folding and entropy-driven molecular recognition. Although so far all these phenomena were studied separately, they stem from the same common root: self-organization of polypeptide chains and, thus, has to be studied collectively. This merge will allow the cross-referencing of a broad spectrum of multi-disciplinary data facilitating progress in all these fields.

Crystal structures of a bacterial dipeptidyl peptidase IV reveal a novel substrate recognition mechanism distinct from that of mammalian orthologues

Dipeptidyl peptidase IV (DPP IV, DPP4, or DAP IV) preferentially cleaves substrate peptides with Pro or Ala at the P1 position. The substrate recognition mechanism has been fully elucidated for mammalian DPP IV by crystal structure analyses but not for bacterial orthologues. Here, we report the crystal structures of a bacterial DPP IV (PmDAP IV) in its free form and in complexes with two kinds of dipeptides as well as with a non-peptidyl inhibitor at 1.90 to 2.47 Å resolution. Acyl-enzyme intermediates were observed for the dipeptide complexes of PmDAP IV, whereas tetrahedral intermediates were reported for the oligopeptide complexes of mammalian DPP IVs. This variation reflects the different structural environments of the active site Arg residues, which are involved in the recognition of a substrate carbonyl group, of mammalian and bacterial enzymes. A phylogenetic analysis revealed that PmDAP IV is a closer relative of dipeptidyl peptidases 8 and 9 (DPP8 and DPP9, DPP IV-family enzymes) than DPP IV. These results provide new insights into the substrate recognition mechanism of bacterial DAP IVs and may assist in the development of selective inhibitors for DAP IVs from pathogenic asaccharolytic bacteria, which utilise proteins or peptides as an energy source.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

Importance of Rigidity in Designing Small Molecule Drugs To Tackle Protein-Protein Interactions (PPIs) through Stabilization of Desired Conformers

Tackling PPIs, particularly by stabilizing clinically favored conformations of target proteins, with orally available, bona fide small molecules remains a significant but immensely worthwhile challenge for the pharmaceutical industry. Success may be more likely through the application of nature's learnings to build intrinsic rigidity into the design of clinical candidates.

The Too Many Faces of PD-L1: A Comprehensive Conformational Analysis Study

In the current study, we focused on the immune-checkpoints PD-1 pathway and in particular on the ligand PD-L1. We studied the conformational dynamics of PD-L1 through principal component analysis of existing crystal structures combined with classical and accelerated molecular dynamics simulations. We identified the maximum structural displacements that take place in all PD-L1 crystal structures and in the molecular dynamics trajectories. We found that these displacements are attributed to specific flexible regions in the protein. We also investigated the conformational preference for small molecule binding and highlighted a methionine residue at the binding site, which plays a key role in drug binding. The binding mechanism of PD-L1 to other binding partners is also discussed in detail from a computational perspective. We hope that the data presented here support the ongoing efforts to discover effective therapies targeting the PD-1 immune-checkpoint pathway.

Osmolyte-Like Stabilizing Effects of Low GdnHCl Concentrations on d-Glucose/d-Galactose-Binding Protein

The ability of d-glucose/d-galactose-binding protein (GGBP) to reversibly interact with its ligands, glucose and galactose, makes this protein an attractive candidate for sensing elements of glucose biosensors. This potential is largely responsible for attracting researchers to study the conformational properties of this protein. Previously, we showed that an increase in the fluorescence intensity of the fluorescent dye 6-bromoacetyl-2-dimetylaminonaphtalene (BADAN) is linked to the holo-form of the GGBP/H152C mutant in solutions containing sub-denaturing concentrations of guanidine hydrochloride (GdnHCl). It was hypothesized that low GdnHCl concentrations might lead to compaction of the protein, thereby facilitating ligand binding. In this work, we utilize BADAN fluorescence spectroscopy, intrinsic protein UV fluorescence spectroscopy, and isothermal titration calorimetry (ITC) to show that the sub-denaturing GdnHCl concentrations possess osmolyte-like stabilizing effects on the structural dynamics, conformational stability, and functional activity of GGBP/H152C and the wild type of this protein (wtGGBP). Our data are consistent with the model where low GdnHCl concentrations promote a shift in the dynamic distribution of the protein molecules toward a conformational ensemble enriched in molecules with a tighter structure and a more closed conformation. This promotes the increase in the configurational complementarity between the protein and glucose molecules that leads to the increase in glucose affinity in both GGBP/H152C and wtGGBP.