Angiotensin Analogs with Divergent Bias Stabilize Distinct Receptor Conformations.

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

Disorder-to-order active site capping regulates the rate-limiting step of the inositol pathway.

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

Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch.

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

Structural basis of protein substrate processing by human mitochondrial high-temperature requirement A2 protease.

The human high-temperature requirement A2 (HtrA2) protein is a trimeric protease that cleaves misfolded proteins to protect cells from stresses caused by toxic, proteinaceous aggregates, and the aberrant function of HtrA2 is closely related to the onset of neurodegenerative disorders. Our methyl-transverse relaxation optimized spectroscopy (TROSY)­based NMR studies using small-peptide ligands have previously revealed a stepwise activation mechanism involving multiple distinct conformational states. However, very little is known about how HtrA2 binds to protein substrates and if the distinct conformational states observed in previous peptide studies might be involved in the processing of protein clients. Herein, we use solution-based NMR spectroscopy to investigate the interaction between the N-terminal Src homology 3 domain from downstream of receptor kinase (drk) with an added C-terminal HtrA2-binding motif (drkN SH3-PDZbm) that exhibits marginal folding stability and serves as a mimic of a physiological protein substrate. We show that drkN SH3-PDZbm binds to HtrA2 via a two-pronged interaction, involving both its C-terminal PDZ-domain binding motif and a central hydrophobic region, with binding occurring preferentially via an unfolded ensemble of substrate molecules. Multivalent interactions between several clients and a single HtrA2 trimer significantly stimulate the catalytic activity of HtrA2, suggesting that binding avidity plays an important role in regulating substrate processing. Our results provide a thermodynamic, kinetic, and structural description of the interaction of HtrA2 with protein substrates and highlight the importance of a trimeric architecture for function as a stress-protective protease that mitigates aggregation.

Transient conformations in the unliganded FK506 binding domain of FKBP51 correspond to two distinct inhibitor-bound states.

Members of the FK506-binding protein (FKBP) family regulate a range of important physiological processes. Unfortunately, current therapeutics such as FK506 and rapamycin exhibit only modest selectivity among these functionally distinct proteins. Recent progress in developing selective inhibitors has been reported for FKBP51 and FKBP52, which act as mutual antagonists in the regulation of steroid hormone signaling. Two structurally similar inhibitors yield distinct protein conformations at the binding site. Localized conformational transition in the binding site of the unliganded FK1 domain of FKBP51 is suppressed by a K58T mutation that also suppresses the binding of these inhibitors. Here, it is shown that the changes in amide hydrogen exchange kinetics arising from this K58T substitution are largely localized to this structural region. Accurate determination of the hydroxide-catalyzed exchange rate constants in both the wildtype and K58T variant proteins impose strong constraints upon the pattern of amide exchange reactivities within either a single or a pair of transient conformations that could give rise to the differences between these two sets of measured rate constants. Poisson-Boltzmann continuum dielectric calculations provide moderately accurate predictions of the structure-dependent hydrogen exchange reactivity for solvent-exposed protein backbone amides. Applying such calculations to the local protein conformations observed in the two inhibitor-bound FKBP51 domains demonstrated that the experimentally determined exchange rate constants for the wildtype domain are robustly predicted by a population-weighted sum of the experimental hydrogen exchange reactivity of the K58T variant and the predicted exchange reactivities in model conformations derived from the two inhibitor-bound protein structures.

Binding Affinity Determination in Drug Design: Insights from Lock and Key, Induced Fit, Conformational Selection, and Inhibitor Trapping Models.

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

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.

Inhibitory and stimulatory micropeptides preferentially bind to different conformations of the cardiac calcium pump.

The ATP-dependent ion pump sarco/endoplasmic reticulum Ca2+-ATPase (SERCA) sequesters Ca2+ in the endoplasmic reticulum to establish a reservoir for cell signaling. Because of its central importance in physiology, the activity of this transporter is tightly controlled via direct interactions with tissue-specific regulatory micropeptides that tune SERCA function to match changing physiological conditions. In the heart, the micropeptide phospholamban (PLB) inhibits SERCA, while dwarf open reading frame (DWORF) stimulates SERCA. These competing interactions determine cardiac performance by modulating the amplitude of Ca2+ signals that drive the contraction/relaxation cycle. We hypothesized that the functions of these peptides may relate to their reciprocal preferences for SERCA binding; SERCA binds PLB more avidly at low cytoplasmic [Ca2+] but binds DWORF better when [Ca2+] is high. In the present study, we demonstrated this opposing Ca2+ sensitivity is due to preferential binding of DWORF and PLB to different intermediate states that SERCA samples during the Ca2+ transport cycle. We show PLB binds best to the SERCA E1-ATP state, which prevails at low [Ca2+]. In contrast, DWORF binds most avidly to E1P and E2P states that are more populated when Ca2+ is elevated. Moreover, FRET microscopy revealed dynamic shifts in SERCA-micropeptide binding equilibria during cellular Ca2+ elevations. A computational model showed that DWORF exaggerates changes in PLB-SERCA binding during the cardiac cycle. These results suggest a mechanistic basis for inhibitory versus stimulatory micropeptide function, as well as a new role for DWORF as a modulator of dynamic oscillations of PLB-SERCA regulatory interactions.

Molecular mechanism of binding between a therapeutic RNA aptamer and its protein target VEGF: A molecular dynamics study.

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

Viewing rare conformations of the ß2 adrenergic receptor with pressure-resolved DEER spectroscopy.

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

Molecular Basis of Coupled Transport and Anion Conduction in Excitatory Amino Acid Transporters.

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

Activation loop dynamics are controlled by conformation-selective inhibitors of ERK2.

Conformational selection by small molecules expands inhibitory possibilities for protein kinases. Nuclear magnetic resonance (NMR) measurements of the mitogen-activated protein (MAP) kinase ERK2 have shown that activation by dual phosphorylation induces global motions involving exchange between two states, L and R. We show that ERK inhibitors Vertex-11e and SCH772984 exploit the small energetic difference between L and R to shift the equilibrium in opposing directions. An X-ray structure of active 2P-ERK2 complexed with AMP-PNP reveals a shift in the Gly-rich loop along with domain closure to position the nucleotide in a more catalytically productive conformation relative to inactive 0P-ERK2:ATP. X-ray structures of 2P-ERK2 complexed with Vertex-11e or GDC-0994 recapitulate this closure, which is blocked in a complex with a SCH772984 analog. Thus, the L→R shift in 2P-ERK2 is associated with movements needed to form a competent active site. Solution measurements by hydrogen-exchange mass spectrometry (HX-MS) reveal distinct binding interactions for Vertex-11e, GDC-0994, and AMP-PNP with active vs. inactive ERK2, where the extent of HX protection correlates with R state formation. Furthermore, Vertex-11e and SCH772984 show opposite effects on HX near the activation loop. Consequently, these inhibitors differentially affect MAP kinase phosphatase activity toward 2P-ERK2. We conclude that global motions in ERK2 reflect conformational changes at the active site that promote productive nucleotide binding and couple with changes at the activation loop to allow control of dephosphorylation by conformationally selective inhibitors.

UvrD helicase activation by MutL involves rotation of its 2B subdomain.

Escherichia coli UvrD is a superfamily 1 helicase/translocase that functions in DNA repair, replication, and recombination. Although a UvrD monomer can translocate along single-stranded DNA, self-assembly or interaction with an accessory protein is needed to activate its helicase activity in vitro. Our previous studies have shown that an Escherichia coli MutL dimer can activate the UvrD monomer helicase in vitro, but the mechanism is not known. The UvrD 2B subdomain is regulatory and can exist in extreme rotational conformational states. By using single-molecule FRET approaches, we show that the 2B subdomain of a UvrD monomer bound to DNA exists in equilibrium between open and closed states, but predominantly in an open conformation. However, upon MutL binding to a UvrD monomer-DNA complex, a rotational conformational state is favored that is intermediate between the open and closed states. Parallel kinetic studies of MutL activation of the UvrD helicase and of MutL-dependent changes in the UvrD 2B subdomain show that the transition from an open to an intermediate 2B subdomain state is on the pathway to helicase activation. We further show that MutL is unable to activate the helicase activity of a chimeric UvrD containing the 2B subdomain of the structurally similar Rep helicase. Hence, MutL activation of the monomeric UvrD helicase is regulated specifically by its 2B subdomain.

Structural and Dynamic Determinants of Molecular Recognition in Bile Acid-Binding Proteins.

Disorders in bile acid transport and metabolism have been related to a number of metabolic disease states, atherosclerosis, type-II diabetes, and cancer. Bile acid-binding proteins (BABPs), a subfamily of intracellular lipid-binding proteins (iLBPs), have a key role in the cellular trafficking and metabolic targeting of bile salts. Within the family of iLBPs, BABPs exhibit unique binding properties including positive binding cooperativity and site-selectivity, which in different tissues and organisms appears to be tailored to the local bile salt pool. Structural and biophysical studies of the past two decades have shed light on the mechanism of bile salt binding at the atomic level, providing us with a mechanistic picture of ligand entry and release, and the communication between the binding sites. In this review, we discuss the emerging view of bile salt recognition in intestinal- and liver-BABPs, with examples from both mammalian and non-mammalian species. The structural and dynamic determinants of the BABP-bile-salt interaction reviewed herein set the basis for the design and development of drug candidates targeting the transcellular traffic of bile salts in enterocytes and hepatocytes.

Exploring the conformational dynamics of PD1 in complex with different ligands: What we can learn for designing novel PD1 signaling blockers?

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

Dynamically encoded reactivity of Ras enzymes: opening new frontiers for drug discovery.

Decoding molecular flexibility in order to understand and predict biological processes-applying the principles of dynamic-structure-activity relationships (DSAR)-becomes a necessity when attempting to design selective and specific inhibitors of a protein that has overlapping interaction surfaces with its upstream and downstream partners along its signaling cascade. Ras proteins are molecular switches that meet this definition perfectly. The close-lying P-loop and the highly flexible switch I and switch II regions are the site of nucleotide-, assisting-, and effector-protein binding. Oncogenic mutations that also appear in this region do not cause easily characterized overall structural changes, due partly to the inherent conformational heterogeneity and pliability of these segments. In this review, we present an overview of the results obtained using approaches targeting Ras dynamics, such as nuclear magnetic resonance (NMR) measurements and experiment-based modeling calculations (mostly molecular dynamics (MD) simulations). These methodologies were successfully used to decipher the mutant- and isoform-specific nature of certain transient states, far-lying allosteric sites, and the internal interaction networks, as well as the interconnectivity of the catalytic and membrane-binding regions. This opens new therapeutic potential: the discovered interaction hotspots present hitherto not targeted, selective sites for drug design efforts in diverse locations of the protein matrix.

Analytical and functional aspects of protein-ligand interactions: Beyond induced fit and conformational selection.

Ligand-dependent changes in protein conformation are foundational to biology. Historical mechanistic models for substrate-specific proteins are induced fit (IF) and conformational selection (CS), which invoke a change in protein conformation after ligand binds or before ligand binds, respectively. These mechanisms have important, but rarely discussed, functional relevance because IF vs. CS can differentially affect a protein's substrate specificity or promiscuity, and its regulatory properties. The modern view of proteins as conformational ensembles in both ligand free and bound states, together with the realization that most proteins exhibit some substrate promiscuity, demands a deeper interpretation of the historical models and provides an opportunity to improve mechanistic analyses. Here we describe alternative analytical strategies for distinguishing the historical models, including the more complex expanded versions of IF and CS. Functional implications of the different models are described. We provide an alternative perspective based on protein ensembles interacting with ligand ensembles that clarifies how a single protein can 'apparently' exploit different mechanisms for different ligands. Mechanistic information about protein ensembles can be optimized when they are probed with multiple ligands.

Unraveling binding mechanism and kinetics of macrocyclic Gαq protein inhibitors.

G proteins represent intracellular switches that transduce signals relayed from G protein-coupled receptors. The structurally related macrocyclic depsipeptides FR900359 (FR) and YM-254890 (YM) are potent, selective inhibitors of the Gαq protein family. We recently discovered that radiolabeled FR and YM display strongly divergent residence times, which translates into significantly longer antiasthmatic effects of FR. The present study is aimed at investigating the molecular basis for this observed disparity. Based on docking studies, we mutated amino acid residues of the Gαq protein predicted to interact with FR or YM, and recombinantly expressed the mutated Gαq proteins in cells in which the native Gαq proteins had been knocked out by CRISPR-Cas9. Both radioligands showed similar association kinetics, and their binding followed a conformational selection mechanism, which was rationalized by molecular dynamics simulation studies. Several mutations of amino acid residues near the putative binding site of the "lipophilic anchors" of FR, especially those predicted to interact with the isopropyl group present in FR but not in YM, led to dramatically accelerated dissociation kinetics. Our data indicate that the long residence time of FR depends on lipophilic interactions within its binding site. The observed structure-kinetic relationships point to a complex binding mechanism of FR, which likely involves snap-lock- or dowel-like conformational changes of either ligand or protein, or both. These experimental data will be useful for the design of compounds with a desired residence time, a parameter that has now been recognized to be of utmost importance in drug development.

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.

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.