Conformational Selection Mechanism Research Articles

Structural determinants of ligand binding in the ternary complex of human ileal bile acid binding protein with glycocholate and glycochenodeoxycholate obtained from solution NMR.

The coordinates of the 10 lowest energy structures of the human I-BABP : GCDA : GCA complex as well as the distance restraints used to calculate the final ensemble have been deposited in the Brookhaven Protein Data Bank with accession number 2MM3.

New open conformation of SMYD3 implicates conformational selection and allostery.

SMYD3 plays a key role in cancer cell viability, adhesion, migration and invasion. SMYD3 promotes formation of inducible regulatory T cells and is involved in reducing autoimmunity. However, the nearly “closed” substrate-binding site and poor in vitro H3K4 methyltransferase activity have obscured further understanding of this oncogenically related protein. Here we reveal that SMYD3 can adopt an “open” conformation using molecular dynamics simulation and small-angle X-ray scattering. This ligand-binding-capable open state is related to the crystal structure-like closed state by a striking clamshell-like inter-lobe dynamics. The two states are characterized by many distinct structural and dynamical differences and the conformational transition pathway is mediated by a reversible twisting motion of the C-terminal domain (CTD). The spontaneous transition from the closed to open states suggests two possible, mutually non-exclusive models for SMYD3 functional regulation and the conformational selection mechanism and allostery may regulate the catalytic or ligand binding competence of SMYD3. This study provides an immediate clue to the puzzling role of SMYD3 in epigenetic gene regulation.

Conformational states of the full-length glucagon receptor.

Class B G protein-coupled receptors are composed of an extracellular domain (ECD) and a seven-transmembrane (7TM) domain, and their signalling is regulated by peptide hormones. Using a hybrid structural biology approach together with the ECD and 7TM domain crystal structures of the glucagon receptor (GCGR), we examine the relationship between full-length receptor conformation and peptide ligand binding. Molecular dynamics (MD) and disulfide crosslinking studies suggest that apo-GCGR can adopt both an open and closed conformation associated with extensive contacts between the ECD and 7TM domain. The electron microscopy (EM) map of the full-length GCGR shows how a monoclonal antibody stabilizes the ECD and 7TM domain in an elongated conformation. Hydrogen/deuterium exchange (HDX) studies and MD simulations indicate that an open conformation is also stabilized by peptide ligand binding. The combined studies reveal the open/closed states of GCGR and suggest that glucagon binds to GCGR by a conformational selection mechanism.

Analysis of the Free Energy Landscapes for the Opening-Closing Dynamics of the Maltose Transporter ATPase MalK2 Using Enhanced-Sampling Molecular Dynamics Simulation.

Protein dynamics are considered significant for many physiological processes, such as metabolism, biomolecular recognition, and the regulation of several vital cellular processes. Due to their flexibility, proteins may stay in different substates with or without the existence of the cognate substrates. To describe these phenomena, two models have been proposed: the "induced fit" and the "conformational selection" mechanisms. In this study, we used MalK2, the subunits that mainly include the nucleotide-binding domains (NBDs) of the maltose transporter from Escherichia coli, as a target to understand the NBD dimerization mechanism. Accelerated and conventional molecular dynamics have been performed. The results revealed that Mg-ATP binding to MalK2 led to a significant change in the free energy profile and thus stabilized the closed conformation. On the contrary, when Mg-ATP was removed, the open conformation would be favored. The fact that ligand binding induces a drastic free energy change leads to a significant inference: MalK2 dimerization would occur through the induced-fit mechanism rather than the conformational selection mechanism. This study sheds new light on the NBD dimerization mechanism and would be of wide applicability to other ABC transporters.

Lys63-linked ubiquitin chain adopts multiple conformational states for specific target recognition.

A polyubiquitin comprises multiple covalently linked ubiquitins and recognizes myriad targets. Free or bound to ligands, polyubiquitins are found in different arrangements of ubiquitin subunits. To understand the structural basis for polyubiquitin quaternary plasticity and to explore the target recognition mechanism, we characterize the conformational space of Lys63-linked diubiquitin (K63-Ub2). Refining against inter-subunit paramagnetic NMR data, we show that free K63-Ub2 exists as a dynamic ensemble comprising multiple closed and open quaternary states. The quaternary dynamics enables K63-Ub2 to be specifically recognized in a variety of signaling pathways. When binding to a target protein, one of the preexisting quaternary states is selected and stabilized. A point mutation that shifts the equilibrium between the different states modulates the binding affinities towards K63-Ub2 ligands. This conformational selection mechanism at the quaternary level may be used by polyubiquitins of different lengths and linkages for target recognition.

New Insights into the Molecular Mechanism of E-Cadherin-Mediated Cell Adhesion by Free Energy Calculations.

Three-dimensional domain swapping is an important mode of protein association leading to the formation of stable dimers. Monomers associating via this mechanism mutually exchange a domain to form a homodimer. Classical cadherins, an increasingly important target for anticancer therapy, use domain swapping to mediate cell adhesion. However, despite its importance, the molecular mechanism of domain swapping is still debated. Here, we study the conformational changes that lead to activation and dimerization via domain swapping of E-cadherin. Using state-of-the-art enhanced sampling atomistic simulations, we reconstruct its conformational free energy landscape, obtaining the free energy profile connecting the inactive and active form. Our simulations predict that the E-cadherin monomer populates the open and closed forms almost equally, which is in agreement with the proposed "selected fit" mechanism in which monomers in an active conformational state bind to form a homodimer, analogous to the conformational selection mechanism often observed in ligand-target binding. Moreover, we find that the open state population is increased in the presence of calcium ions at the extracellular boundary, suggesting their possible role as allosteric activators of the conformational change.

The PRE-Derived NMR Model of the 38.8-kDa Tri-Domain IsdH Protein from Staphylococcus aureus Suggests That It Adaptively Recognizes Human Hemoglobin

Staphylococcus aureus is a medically important bacterial pathogen that, during infections, acquires iron from human hemoglobin (Hb). It uses two closely related iron-regulated surface determinant (Isd) proteins to capture and extract the oxidized form of heme (hemin) from Hb, IsdH and IsdB. Both receptors rapidly extract hemin using a conserved tri-domain unit consisting of two NEAT (near iron transporter) domains connected by a helical linker domain. To gain insight into the mechanism of extraction, we used NMR to investigate the structure and dynamics of the 38.8-kDa tri-domain IsdH protein (IsdHN2N3, A326–D660 with a Y642A mutation that prevents hemin binding). The structure was modeled using long-range paramagnetic relaxation enhancement (PRE) distance restraints, dihedral angle, small-angle X-ray scattering, residual dipolar coupling and inter-domain NOE nuclear Overhauser effect data. The receptor adopts an extended conformation wherein the linker and N3 domains pack against each other via a hydrophobic interface. In contrast, the N2 domain contacts the linker domain via a hydrophilic interface and, based on NMR relaxation data, undergoes inter-domain motions enabling it to reorient with respect to the body of the protein. Ensemble calculations were used to estimate the range of N2 domain positions compatible with the PRE data. A comparison of the Hb-free and Hb-bound forms reveals that Hb binding alters the positioning of the N2 domain. We propose that binding occurs through a combination of conformational selection and induced-fit mechanisms that may promote hemin release from Hb by altering the position of its F helix.

Investigation of arc repressor DNA-binding specificity by comparative molecular dynamics simulations

Transcription factors regulate gene expression through binding to specific DNA sequences. How transcription factors achieve high binding specificity is still not well understood. In this paper, we investigated the role of protein flexibility in protein–DNA-binding specificity by comparative molecular dynamics (MD) simulations. Protein flexibility has been considered as a key factor in molecular recognition, which is intrinsically a dynamic process involving fine structural fitting between binding components. In this study, we performed comparative MD simulations on wild-type and F10V mutant P22 Arc repressor in both free and complex conformations. The F10V mutant has lower DNA-binding specificity though both the bound and unbound main-chain structures between the wild-type and F10V mutant Arc are highly similar. We found that the DNA-binding motif of wild-type Arc is structurally more flexible than the F10V mutant in the unbound state, especially for the six DNA base-contacting residues in each dimer. We demonstrated that the flexible side chains of wild-type Arc lead to a higher DNA-binding specificity through forming more hydrogen bonds with DNA bases upon binding. Our simulations also showed a possible conformational selection mechanism for Arc-DNA binding. These results indicate the important roles of protein flexibility and dynamic properties in protein–DNA-binding specificity.

Kinase Structural Dynamics Enables Tight and Selective Binding of Inhibitors

Protein kinases are obvious drug targets against cancer due to their central role in cellular regulation. With oncologic diseases being the second leading cause of death in the US kinases rapidly gain attention and are likely to become the number one drug target. Using NMR and fast kinetics, we establish a novel model that solves a longstanding question of high selectivity of clinically relevant drug Gleevec that effectively inhibits Abl tyrosine kinase while closely related Src family of kinases is affected much less. Our study of an entirely different family of Ser/Thr Aurora kinases and its specific inhibitors suggests that an energy landscape that provides tight affinity via an induced-fit and binding plasticity via a conformational selection mechanism is likely to be general for many inhibitors.

Spin labeling and Double Electron-Electron Resonance (DEER) to Deconstruct Conformational Ensembles of HIV Protease.

An understanding of macromolecular conformational equilibrium in biological systems is oftentimes essential to understand function, dysfunction, and disease. For the past few years, our lab has been utilizing site-directed spin labeling (SDSL), coupled with electron paramagnetic resonance (EPR) spectroscopy, to characterize the conformational ensemble and ligand-induced conformational shifts of HIV-1 protease (HIV-1PR). The biomedical importance of characterizing the fractional occupancy of states within the conformational ensemble critically impacts our hypothesis of a conformational selection mechanism of drug-resistance evolution in HIV-1PR. The purpose of the following chapter is to give a timeline perspective of our SDSL EPR approach to characterizing conformational sampling of HIV-1PR. We provide detailed instructions for the procedure utilized in analyzing distance profiles for HIV-1PR obtained from pulsed electron-electron double resonance (PELDOR). Specifically, we employ a version of PELDOR known as double electron-electron resonance (DEER). Data are processed with the software package "DeerAnalysis" (http://www.epr.ethz.ch/software), which implements Tikhonov regularization (TKR), to generate a distance profile from electron spin-echo amplitude modulations. We assign meaning to resultant distance profiles based upon a conformational sampling model, which is described herein. The TKR distance profiles are reconstructed with a linear combination of Gaussian functions, which is then statistically analyzed. In general, DEER has proven powerful for observing structural ensembles in proteins and, more recently, nucleic acids. Our goal is to present our advances in order to aid readers in similar applications.

Free Energy Landscapes for a Kainate Receptor Ligand-Binding Domain

The kainate family of ionotropic glutamate receptors (GluK) mediate neurotransmission in pre- and post-synaptic neurons. Previously computed free energy landscapes for iGluR ligand-binding domains (LBD) revealed subtype-dependent conformational differences between NMDA and AMPA iGluR families. Here we report free energy landscapes for the GluK2 ligand-binding domain (LBD) in apo and glutamate-bound states. The GluK2 free energy landscapes reveal both apo and glutamate-bound LBDs preferentially access the closed-cleft conformation and suggest a conformational selection mechanism of binding. These results are consistent with the free energy landscapes for NMDA receptors but not AMPA receptors. Additionally, the open-cleft conformation is more accessible to the GluK2 apo state than the glutamate-bound state, consistent with previous results for both NMDA and AMPA receptors. Finally, large-scale structural dynamics are characterized by principle component analysis.

BAX-Induced Apoptosis Can Be Initiated through a Conformational Selection Mechanism

BAX protein plays a key role in the mitochondria-mediated apoptosis. However, it remains unclear by what mechanism BAX is triggered to initiate apoptosis. Here, we reveal the mechanism using electron spin resonance (ESR) techniques. An inactive BAX monomer was found to exhibit conformational heterogeneity and exist at equilibrium in two conformations, one of which has never been reported. We show that upon apoptotic stimulus by BH3-only peptides, BAX can be induced to convert into either a ligand-bound monomer or an oligomer through a conformational selection mechanism. The kinetics of reaction is studied by means of time-resolved ESR, allowing a direct insitu observation for the transformation of BAX from the native to the bound states. Invitro mitochondrial assays provide further discrimination between the proposed BAX states, thereby revealing a population-shift allosteric mechanism in the process. BAX's apoptotic function is shown to critically depend on excursions between different structural conformations.

Microsecond folding and domain motions of a spider silk protein structural switch.

Web spiders rapidly assemble protein monomers, so-called spidroins, into extraordinarily tough silk fibers. The process involves the pH-triggered self-association of the spidroin N-terminal domain (NTD), which contains a structural switch connecting spidroins to supermolecules. Single-molecule spectroscopy can detect conformational heterogeneity that is hidden to conventional methods, but motions of the NTD are beyond the resolution limit. Here, we engineered probes for 1 nm conformational changes based on the phenomenon of fluorescence quenching by photoinduced electron transfer into the isolated NTD of a spidroin from the nursery web spider Euprosthenops australis. Correlation analysis of single-molecule fluorescence fluctuations uncovered site-dependent nanosecond-to-microsecond movement of secondary and tertiary structure. Kinetic amplitudes were most pronounced for helices that are part of the association interface and where structural studies show large displacements between monomeric and dimeric conformations. A single tryptophan at the center of the five-helix bundle toggled conformations in ∼100 μs and in a pH-dependent manner. Equilibrium denaturation and temperature-jump relaxation experiments revealed cooperative and ultrafast folding in only 60 μs. We deduced a free-energy surface that exhibits native-state ruggedness with apparently similar barrier heights to folding and native motions. Observed equilibrium dynamics within the domain suggest a conformational selection mechanism in the rapid association of spidroins through their NTDs during silk synthesis by web spiders.

Ligand‐Detected Relaxation Dispersion NMR Spectroscopy: Dynamics of preQ1–RNA Binding

AbstractAn NMR‐based approach to characterizing the binding kinetics of ligand molecules to biomolecules, like RNA or proteins, by ligand‐detected Carr‐Purcell‐Meiboom‐Gill (CPMG) relaxation dispersion experiments is described. A 15N‐modified preQ1 ligand is used to acquire relaxation dispersion experiments in the presence of low amounts of the Fsu class I preQ1 aptamer RNA, and increasing ligand concentrations to probe the RNA small molecule interaction. Our experimental data strongly support the conformational selection mechanism postulated. The approach gives direct access to two parameters of a ligand–receptor interaction: the off rate and the population of the small molecule–receptor complex. A detailed description of the kinetics underlying the ligand binding process is of crucial importance to fully understanding a riboswitch’s function and to evaluate potential new antibiotics candidates targeting the noncoding RNA species. Ligand‐detected NMR relaxation dispersion experiments represent a valuable diagnostic tool for the characterization of binding mechanisms.

Ligand-detected relaxation dispersion NMR spectroscopy: dynamics of preQ1-RNA binding.

An NMR-based approach to characterizing the binding kinetics of ligand molecules to biomolecules, like RNA or proteins, by ligand-detected Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments is described. A (15)N-modified preQ1 ligand is used to acquire relaxation dispersion experiments in the presence of low amounts of the Fsu class I preQ1 aptamer RNA, and increasing ligand concentrations to probe the RNA small molecule interaction. Our experimental data strongly support the conformational selection mechanism postulated. The approach gives direct access to two parameters of a ligand-receptor interaction: the off rate and the population of the small molecule-receptor complex. A detailed description of the kinetics underlying the ligand binding process is of crucial importance to fully understanding a riboswitch's function and to evaluate potential new antibiotics candidates targeting the noncoding RNA species. Ligand-detected NMR relaxation dispersion experiments represent a valuable diagnostic tool for the characterization of binding mechanisms.

Dual effects of familial Alzheimer's disease mutations (D7H, D7N, and H6R) on amyloid β peptide: correlation dynamics and zinc binding.

Although the N-terminal region of Amyloid β (Aβ) peptides plays dual roles as metal-coordinating sites and conformational modulator, few studies have been performed to explore the effects of mutations at this region on the overall conformational ensemble of Aβ and the binding propensity of metal ions. In this work, we focus on how three familial Alzheimer's disease mutations (D7H, D7N, and H6R) alter the structural characteristics and thermodynamic stabilities of Aβ42 using molecular dynamics simulations. We observe that each mutation displays increased β-sheet structures in both N and C termini. In particular, both the N terminus and central hydrophobic region of D7H can form stable β-hairpin structures with its C terminus. The conserved turn structure at Val²⁴-Lys²⁸ in all peptides and Zn²⁺-bound Aβ42 is confirmed as the common structural motif to nucleate folding of Aβ. Each mutant can significantly increase the solvation free energy and thus enhance the aggregation of Aβ monomers. The correlation dynamics between Aβ(1-16) and Aβ(17-42) fragments are elucidated by linking the domain motions with the corresponding structured conformations. We characterize the different populations of correlated domain motions for each mutant from a more macroscopic perspective, and unexpectedly find that Zn²⁺-bound Aβ42 ensemble shares the same populations as Aβ42, indicating that the binding of Zn²⁺ to Aβ follows the conformational selection mechanism, and thus is independent of domain motions, even though the structures of Aβ have been modified at a residue level.

Specific recognition mechanism between RNA and the KH3 domain of Nova-2 protein.

The KH3 domain of Nova-2 protein can precisely recognize the sequence-specific target RNA of human glycine receptor α2. However, the recognition mechanism between the protein and its target RNA is still hotly debated. In this study, molecular dynamic simulations in explicit solvent were utilized to understand the recognition mechanism. The structural analysis and the Kolmogorov-Smirnov P-test statistics reveal that the KH3 domain might obey a conformational selection mechanism upon RNA binding. However, the induced fit mechanism could not be completely ruled out. Unfolding kinetics indicates that the folding of RNA and KH3 happens first and then the binding between RNA and KH3 follows. Principle component analysis shows that the invariant Gly-Lys-Gly-Gly loop moves toward to the RNA molecule but the C-terminal domain moves away from the RNA molecule upon binding. These specific dominant motions were hypothesized to stabilize the complex structure. The hydrophobic and hydrogen bonding interactions were found to be the driving forces for the specific recognition, in contrast to the dominant electrostatic interactions for nonspecific recognition.

Computational study of Gleevec and G6G reveals molecular determinants of kinase inhibitor selectivity.

Gleevec is a potent inhibitor of Abl tyrosine kinase but not of the highly homologous c-Src kinase. Because the ligand binds to an inactive form of the protein in which an Asp-Phe-Gly structural motif along the activation loop adopts a so-called DFG-out conformation, it was suggested that binding specificity was controlled by a “conformational selection” mechanism. In this context, the binding affinity displayed by the kinase inhibitor G6G poses an intriguing challenge. Although it possesses a chemical core very similar to that of Gleevec, G6G is a potent inhibitor of both Abl and c-Src kinases. Both inhibitors bind to the DFG-out conformation of the kinases, which seems to be in contradiction with the conformational selection mechanism. To address this issue and display the hidden thermodynamic contributions affecting the binding selectivity, molecular dynamics free energy simulations with explicit solvent molecules were carried out. Relative to Gleevec, G6G forms highly favorable van der Waals dispersive interactions upon binding to the kinases via its triazine functional group, which is considerably larger than the corresponding pyridine moiety in Gleevec. Upon binding of G6G to c-Src, these interactions offset the unfavorable free energy cost of the DFG-out conformation. When binding to Abl, however, G6G experiences an unfavorable free energy penalty due to steric clashes with the phosphate-binding loop, yielding an overall binding affinity that is similar to that of Gleevec. Such steric clashes are absent when G6G binds to c-Src, due to the extended conformation of the phosphate-binding loop.

Using Markov State Models to Develop a Mechanistic Understanding of Protein Kinase A Regulatory Subunit RIα Activation in Response to cAMP Binding

Protein kinase A (PKA) holoenzyme consists of two catalytic (C) subunits and a regulatory (R) subunit dimer (R(2)C(2)). The kinase is activated by the binding of cAMPs to the two cyclic nucleotide binding domains (CBDs), A and B, on each R-subunit. Despite extensive study, details of the allosteric mechanisms underlying the cooperativity of holoenzyme activation remain unclear. Several Markov state models of PKA-RIα were developed to test competing theories of activation for the R(2)C(2) complex. We found that CBD-B plays an essential role in R-C interaction and promotes the release of the first C-subunit prior to the binding to CBD-A. This favors a conformational selection mechanism for release of the first C-subunit of PKA. However, the release of the second C-subunit requires all four cAMP sites to be occupied. These analyses elucidate R-C heterodimer interactions in the cooperative activation of PKA and cAMP binding and represent a new mechanistic model of R(2)C(2) PKA-RIα activation.

Energetic dissection of Gleevec's selectivity toward human tyrosine kinases.

Protein kinases are obvious drug targets against cancer due to their central role in cellular regulation. Since the discovery of Gleevec, a potent and specific inhibitor of Abl kinase, as a highly successful cancer therapeutic, the ability of this drug to distinguish between Abl and other tyrosine kinases like Src has been intensely investigated, but without much success. Using NMR and fast kinetics, we establish a novel model that solves this longstanding question of two tyrosine kinases adopting almost identical structures when bound to Gleevec, yet having vastly different affinities. In contrast to all other proposed models we show that the origin of Abl’s high affinity lies predominantly in a conformational change after binding. An energy landscape that provides tight affinity via an induced-fit and binding plasticity via conformational selection mechanism is likely to be general for many inhibitors.