Binding-competent Conformation Research Articles

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.

Ligand-induced protein transition state stabilization switches the binding pathway from conformational selection to induced fit

Protein-ligand complex formation is fundamental to biological function. A central question is whether proteins spontaneously adopt binding-competent conformations to which ligands bind conformational selection (CS) or whether ligands induce the binding-competent conformation induced fit (IF). Here, we resolve the CS and IF binding pathways by characterizing protein conformational dynamics over a wide range of ligand concentrations using NMR relaxation dispersion. We determined the relative flux through the two pathways using a four-state binding model that includes both CS and IF. Experiments conducted without ligand show that galectin-3 exchanges between the ground-state conformation and a high-energy conformation similar to the ligand-bound conformation, demonstrating that CS is a plausible pathway. Near-identical crystal structures of the apo and ligand-bound states suggest that the high-energy conformation in solution corresponds to the apo crystal structure. Stepwise additions of the ligand lactose induce progressive changes in the relaxation dispersions that we fit collectively to the four-state model, yielding all microscopic rate constants and binding affinities. The ligand affinity is higher for the bound-like conformation than for the ground state, as expected for CS. Nonetheless, the IF pathway contributes greater than 70% of the total flux even at low ligand concentrations. The higher flux through the IF pathway is explained by considerably higher rates of exchange between the two protein conformations in the ligand-associated state. Thus, the ligand acts to decrease the activation barrier between protein conformations in a manner reciprocal to enzymatic transition-state stabilization of reactions involving ligand transformation.

Structural Characterization of Nanobodies during Germline Maturation.

Camelid heavy-chain antibody variable domains (VHH), nanobodies, are the smallest-known functional antibody fragments with high therapeutic potential. In this study, we investigate a VHH binding to hen egg-white lysozyme (HEL). We structurally and dynamically characterized the conformational diversity of four VHH variants to elucidate the antigen-binding process. For two of these antibodies, not only are the dissociation constants known, but also the experimentally determined crystal structures of the VHH in complex with HEL are available. We performed well-tempered metadynamics simulations in combination with molecular dynamics simulations to capture a broad conformational space and to reconstruct the thermodynamics and kinetics of conformational transitions in the antigen-binding site, the paratope. By kinetically characterizing the loop movements of the paratope, we found that, with an increase in affinity, the state populations shift towards the binding competent conformation. The contacts contributing to antigen binding, and those who contribute to the overall stability, show a clear trend towards less variable but more intense contacts. Additionally, these investigated nanobodies clearly follow the conformational selection paradigm, as the binding competent conformation pre-exists within the structural ensembles without the presence of the antigen.

Multimeric structure enables the acceleration of KaiB-KaiC complex formation induced by ADP/ATP exchange inhibition.

Circadian clocks tick a rhythm with a nearly 24-hour period in a variety of organisms. In the clock proteins of cyanobacteria, KaiA, KaiB, and KaiC, known as a minimum circadian clock, the slow KaiB-KaiC complex formation is essential in determining the clock period. This complex formation, occurring when the C1 domain of KaiC hexamer binds ADP molecules produced by the ATPase activity of C1, is considered to be promoted by accumulating ADP molecules in C1 through inhibiting the ADP/ATP exchange (ADP release) rather than activating the ATP hydrolysis (ADP production). Significantly, this ADP/ATP exchange inhibition accelerates the complex formation together with its promotion, implying a potential role in the period robustness under environmental perturbations. However, the molecular mechanism of this simultaneous promotion and acceleration remains elusive because inhibition of a backward process generally slows down the whole process. In this article, to investigate the mechanism, we build several reaction models of the complex formation with the pre-binding process concerning the ATPase activity. In these models, six KaiB monomers cooperatively and rapidly bind to C1 when C1 binds ADP molecules more than a given threshold while stabilizing the binding-competent conformation of C1. Through comparison among the models proposed here, we then extract three requirements for the simultaneous promotion and acceleration: the stabilization of the binding-competent C1 by KaiB binding, slow ADP/ATP exchange in the binding-competent C1, and relatively fast ADP/ATP exchange occurring in the binding-incompetent C1 in the presence of KaiB. The last two requirements oblige KaiC to form a multimer. Moreover, as a natural consequence, the present models can also explain why the binding of KaiB to C1 reduces the ATPase activity of C1.

Paratope states in solution improve structure prediction and docking

Structure-based antibody design and accurate predictions of antibody-antigen interactions remain major challenges in computational biology. By using molecular dynamics simulations, we show that a single static X-ray structure is not sufficient to identify determinants of antibody-antigen recognition. Here, we investigate antibodies that undergo substantial conformational changes upon antigen binding and have been classified as difficult cases in an extensive benchmark for antibody-antigen docking. We present thermodynamics and transition kinetics of these conformational rearrangements and show that paratope states can be used to improve antibody-antigen docking. By using the unbound antibody X-ray structure as starting structure for molecular dynamics simulations, we retain a binding competent conformation substantially different to the unbound antibody X-ray structure. We also observe that the kinetically dominant antibody paratope conformations are chosen by the bound antigen conformation with the highest probability. Thus, we show that paratope states in solution can improve antibody-antigen docking and structure prediction.

Designer proteins that competitively inhibit Gαq by targeting its effector site

During signal transduction, the G protein, Gαq, binds and activates phospholipase C-β isozymes. Several diseases have been shown to manifest upon constitutively activating mutation of Gαq, such as uveal melanoma. Therefore, methods are needed to directly inhibit Gαq. Previously, we demonstrated that a peptide derived from a helix-turn-helix (HTH) region of PLC-β3 (residues 852–878) binds Gαq with low micromolar affinity and inhibits Gαq by competing with full-length PLC-β isozymes for binding. Since the HTH peptide is unstructured in the absence of Gαq, we hypothesized that embedding the HTH in a folded protein might stabilize the binding-competent conformation and further improve the potency of inhibition. Using the molecular modeling software Rosetta, we searched the Protein Data Bank for proteins with similar HTH structures near their surface. The candidate proteins were computationally docked against Gαq, and their surfaces were redesigned to stabilize this interaction. We then used yeast surface display to affinity mature the designs. The most potent design bound Gαq/i with high affinity in vitro (KD = 18 nM) and inhibited activation of PLC-β isozymes in HEK293 cells. We anticipate that our genetically encoded inhibitor will help interrogate the role of Gαq in healthy and disease model systems. Our work demonstrates that grafting interaction motifs into folded proteins is a powerful approach for generating inhibitors of protein–protein interactions.

Computationally designed peptide macrocycle inhibitors of New Delhi metallo-β-lactamase 1

The rise of antibiotic resistance calls for new therapeutics targeting resistance factors such as the New Delhi metallo-β-lactamase 1 (NDM-1), a bacterial enzyme that degrades β-lactam antibiotics. We present structure-guided computational methods for designing peptide macrocycles built from mixtures of l- and d-amino acids that are able to bind to and inhibit targets of therapeutic interest. Our methods explicitly consider the propensity of a peptide to favor a binding-competent conformation, which we found to predict rank order of experimentally observed IC50 values across seven designed NDM-1- inhibiting peptides. We were able to determine X-ray crystal structures of three of the designed inhibitors in complex with NDM-1, and in all three the conformation of the peptide is very close to the computationally designed model. In two of the three structures, the binding mode with NDM-1 is also very similar to the design model, while in the third, we observed an alternative binding mode likely arising from internal symmetry in the shape of the design combined with flexibility of the target. Although challenges remain in robustly predicting target backbone changes, binding mode, and the effects of mutations on binding affinity, our methods for designing ordered, binding-competent macrocycles should have broad applicability to a wide range of therapeutic targets.

Challenge Accepted – Paratope States in Solution Improve Structure Prediction and Docking

The rise of antibodies as biotherapeutic proteins has attracted the interest in characterizing and understanding the antibody binding interface. Structure-based antibody design and accurate predictions of antibody-antigen interactions remain major challenges in computational biology. By using molecular dynamics simulations in combination with enhanced sampling techniques, we show that a single static X-ray structure is not sufficient to identify determinants of antibody recognition and to characterize functionally relevant conformational rearrangements. Thus, we suggest that antibodies should be described as ensembles of paratope states in solution, which are defined by a combination of correlated loop movements, which favor distinct interface orientations and elbow angles. In this study, we focus our attention on characterizing paratope states in solution of antibodies which undergo substantial conformational changes upon antigen binding. Because of their conformational changes, they have been classified as “difficult” cases for antibody-antigen docking in the Critical Assessment of PRediction of Interactions (CAPRI) blind docking challenge. Here, we present transition kinetics and thermodynamics of these conformational rearrangements of the paratope and show that kinetically relevant states can be used to improve antibody-antigen docking. We use the Rosetta SnugDock docking algorithm to show that by even using the unbound antibody X-ray structure as starting structures for molecular dynamics simulations, a binding competent conformation substantially different to the unbound antibody X-ray structure, can be retained. We also observe that the kinetically dominant antibody conformations are chosen by the bound antigen conformation with the highest probability. Thus, we show that kinetically relevant antibody conformations obtained by a combination of molecular dynamics simulations and Markov-state models can improve antibody-antigen docking and structure prediction.

Zinc-mediated conformational preselection mechanism in the allosteric control of DNA binding to the zinc transcriptional regulator (ZitR)

The zinc transcriptional regulator (ZitR) functions as a metalloregulator that fine tunes transcriptional regulation through zinc-dependent DNA binding. However, the molecular mechanism of zinc-driven allosteric control of the DNA binding to ZitR remains elusive. Here, we performed enhanced sampling accelerated molecular dynamics simulations to figure out the mechanism, revealing the role of protein dynamics in the zinc-induced allosteric control of DNA binding to ZitR. The results suggest that zinc-free ZitR samples distinct conformational states, only a handful of which are compatible with DNA binding. Remarkably, zinc binding reduces the conformational plasticity of the DNA-binding domain of ZitR, promoting the population shift in the ZitR conformational ensemble towards the DNA binding-competent conformation. Further co-binding of DNA to the zinc–ZitR complex stabilizes this competent conformation. These findings suggest that ZitR–DNA interactions are allosterically regulated in a zinc-mediated conformational preselection manner, highlighting the importance of conformational dynamics in the regulation of transcription factor family.

Sodium-induced population shift drives activation of thrombin

The equilibrium between active E and inactive E* forms of thrombin is assumed to be governed by the allosteric binding of a Na+ ion. Here we use molecular dynamics simulations and Markov state models to sample transitions between active and inactive states. With these calculations we are able to compare thermodynamic and kinetic properties depending on the presence of Na+. For the first time, we directly observe sodium-induced conformational changes in long-timescale computer simulations. Thereby, we are able to explain the resulting change in activity. We observe a stabilization of the active form in presence of Na+ and a shift towards the inactive form in Na+-free simulations. We identify key structural features to quantify and monitor this conformational shift. These include the accessibility of the S1 pocket and the reorientation of W215, of R221a and of the Na+ loop. The structural characteristics exhibit dynamics at various timescales: Conformational changes in the Na+ binding loop constitute the slowest observed movement. Depending on its orientation, it induces conformational shifts in the nearby substrate binding site. Only after this shift, residue W215 is able to move freely, allowing thrombin to adopt a binding-competent conformation.

Binding between G Quadruplexes at the Homodimer Interface of the Corn RNA Aptamer Strongly Activates Thioflavin T Fluorescence

Thioflavin T (ThT) is widely used for the detection of amyloids. Many unrelated DNAs and RNAs that contain G-quadruplex motifs also bind ThT and strongly activate its fluorescence. To elucidate the structural basis of ThT binding to G quadruplexes and its fluorescence turn-on, we determined its co-crystal structure with the homodimeric RNA Corn, which contains two G quadruplexes. We found that two ThT molecules bind in the dimer interface, constrained by a G quartet from each protomer into a maximally fluorescent planar conformation. The unliganded Corn homodimer crystal structure reveals a collapsed fluorophore-binding site. In solution, Corn must fluctuate between this and an open, binding-competent conformation. A co-crystal structure with another benzothiazole derivate, thiazole orange (TO), also shows binding at the Corn homodimer interface. As the bound ThT and TO make no interactions with the RNA backbone, their Corn co-crystal structures likely explain their fluorescence activationupon sequence-independent DNA and RNA G-quadruplex binding.

Mechanistic Insights into Microsecond Time-Scale Motion of Solid Proteins Using Complementary 15N and 1H Relaxation Dispersion Techniques.

NMR relaxation dispersion methods provide a holistic way to observe microsecond time-scale protein backbone motion both in solution and in the solid state. Different nuclei (1H and 15N) and different relaxation dispersion techniques (Bloch-McConnell and near-rotary-resonance) give complementary information about the amplitudes and time scales of the conformational dynamics and provide comprehensive insights into the mechanistic details of the structural rearrangements. In this paper, we exemplify the benefits of the combination of various solution- and solid-state relaxation dispersion methods on a microcrystalline protein (α-spectrin SH3 domain), for which we are able to identify and model the functionally relevant conformational rearrangements around the ligand recognition loop occurring on multiple microsecond time scales. The observed loop motions suggest that the SH3 domain exists in a binding-competent conformation in dynamic equilibrium with a sterically impaired ground-state conformation both in solution and in crystalline form. This inherent plasticity between the interconverting macrostates is compatible with a conformational-preselection model and provides new insights into the recognition mechanisms of SH3 domains.

Structure of tRNA-Modifying Enzyme TiaS and Motions of Its Substrate Binding Zinc Ribbon

The accurate modification of the tRNAIle anticodon wobble cytosine 34 is critical for AUA decoding in protein synthesis. Archaeal tRNAIle2 cytosine 34 is modified with agmatine in the presence of ATP by TiaS (tRNAIle2 agmatidine synthetase). However, no structure of apo-form full-length TiaS is available currently. Here, the crystal structures of apo TiaS and a complex of TiaS–agmatine–AMPPCP–Mg are presented, with properly folded zinc ribbon and Cys4-zinc coordination identified. Compared with tRNAIle2-bound form, the architecture of apo TiaS shows a totally different conformation of zinc ribbon. Molecular dynamics simulations of the docking complex between free-state TiaS and tRNAIle2 suggest that zinc ribbon domain is capable of performing large-scale motions to sample substrate binding-competent conformation. Principle component analysis and normal mode analysis show consistent results about the relative directionality of functionally correlated zinc ribbon motions. Apo TiaS and TiaS–agmatine–AMPPCP–Mg/TiaS–AMPCPP–Mg complex structures capture two snapshots of the flexible ATP-Mg binding p2loop step-by-step stabilization. Research from this study provides new insight into TiaS functional mechanism and the dynamic feature of zinc ribbons.

Structures of chaperone-substrate complexes docked onto the export gate in a type III secretion system

The flagellum and the injectisome enable bacterial locomotion and pathogenesis, respectively. These nanomachines assemble and function using a type III secretion system (T3SS). Exported proteins are delivered to the export apparatus by dedicated cytoplasmic chaperones for their transport through the membrane. The structural and mechanistic basis of this process is poorly understood. Here we report the structures of two ternary complexes among flagellar chaperones (FliT and FliS), protein substrates (the filament-capping FliD and flagellin FliC), and the export gate platform protein FlhA. The substrates do not interact directly with FlhA; however, they are required to induce a binding-competent conformation to the chaperone that exposes the recognition motif featuring a highly conserved sequence recognized by FlhA. The structural data reveal the recognition signal in a class of T3SS proteins and provide new insight into the assembly of key protein complexes at the export gate.

Proline Restricts Loop I Conformation of the High Affinity WW Domain from Human Nedd4-1 to a Ligand Binding-Competent Type I β-Turn.

Sequence alignment of the four WW domains from human Nedd4-1 (neuronal precursor cell expressed developmentally down-regulated gene 4-1) reveals that the highest sequence diversity exists in loop I. Three residues in this type I β-turn interact with the PPxY motif of the human epithelial Na+ channel (hENaC) subunits, indicating that peptide affinity is defined by the loop I sequence. The third WW domain (WW3*) has the highest ligand affinity and unlike the other three hNedd4-1 WW domains or other WW domains studied contains the highly statistically preferred proline at the ( i + 1) position found in β-turns. In this report, molecular dynamics simulations and experimental data were combined to characterize loop I stability and dynamics. Exchange of the proline to the equivalent residue in WW4 (Thr) results in the presence of a predominantly open seven residue Ω loop rather than the type I β-turn conformation for the wild-type apo-WW3*. In the presence of the ligand, the structure of the mutated loop I is locked into a type I β-turn. Thus, proline in loop I ensures a stable peptide binding-competent β-turn conformation, indicating that amino acid sequence modulates local flexibility to tune binding preferences and stability of dynamic interaction motifs.

Global Disordering in Stereo-Specific Protein Association

Protein-protein recognition is of fundamental importance for a myriad of biological processes and is ultimately a prerequisite for life as we know it. There exist several established mechanisms that promote formation of stereo-specific protein complexes. Many of these mechanism involve conformational changes of one or both proteins in dimeric assemblies as observed in “conformational selection” and “coupled folding and binding” scenarios. In “coupled folding and binding” events, at least one of the proteins undergoes a global ordering event. By using an integrated computational and experimental approach we have discovered that also global disordering can be a productive route for formation of a stereo-specific protein complex. This mechanism was observed for the chaperone binding domain of the Yersinia effector protein YopH upon binding to its specific chaperone SycH. These two proteins are crucial for type III secretion system mediated infectivity by Yersinia and several other gram negative pathogens. NMR relaxation dispersion experiments demonstrated that the otherwise well folded YopH protein dynamically samples an expanded high-energy state that corresponds to the SycH binding competent conformation. A structure of the protein complex determined from a hybrid SAXS and computational approach revealed that YopH wraps around SycH in a horse shoe like conformation. The binding model was validated by site specific YopH mutations that promoted the disordering event and at the same time displayed improved binding affinity towards SycH. Taken together the data illustrates a tight coupling between a proteins unfolding and functional free energy landscapes and add valuable mechanistic insight into protein-protein recognition.

Substrate Recognition and Activity Regulation of the Escherichia coli mRNA Endonuclease MazF

Escherichia coli MazF (EcMazF) is the archetype of a large family of ribonucleases involved in bacterial stress response. The crystal structure of EcMazF in complex with a 7-nucleotide substrate mimic explains the relaxed substrate specificity of the E. coli enzyme relative to its Bacillus subtilis counterpart and provides a framework for rationalizing specificity in this enzyme family. In contrast to a conserved mode of substrate recognition and a conserved active site, regulation of enzymatic activity by the antitoxin EcMazE diverges from its B. subtilis homolog. Central in this regulation is an EcMazE-induced double conformational change as follows: a rearrangement of a crucial active site loop and a relative rotation of the two monomers in the EcMazF dimer. Both are induced by the C-terminal residues Asp-78-Trp-82 of EcMazE, which are also responsible for strong negative cooperativity in EcMazE-EcMazF binding. This situation shows unexpected parallels to the regulation of the F-plasmid CcdB activity by CcdA and further supports a common ancestor despite the different activities of the MazF and CcdB toxins. In addition, we pinpoint the origin of the lack of activity of the E24A point mutant of EcMazF in its inability to support the substrate binding-competent conformation of EcMazF.

Dissecting the Contribution of Kinase Conformational Reorganization Energies to Inhibitor Selectivity

Kinases are an important target for many cancer therapies. However, they can be difficult to selectively target without hitting one of the other ∼500 kinases in the human genome. Recent experimental and computational evidence suggests that the cost of confining the kinase to the binding-competent conformation could play a large role in determining inhibitor selectivity, specifically in the case of Src and Abl kinase and the ligand imatinib (also known as the paradigm shifting targeted cancer drug ‘Gleevec’). To test this hypothesis on a larger scale, and to dissect the contribution of this reorganization energy to kinase inhibitor affinity and selectivity, we have developed an approach that combines massively distributed molecular simulations with automated experimental biophysical measurements. Atomistic molecular simulations are run on the Folding@home worldwide distributed computing platform to collect ∼ms of simulation data. Markov state models are then constructed of this data to determine the accessible conformations and relative free energies. Alchemical binding free energy calculations are used to assess the ligand binding free energies of selective kinase inhibitors to individual kinetically metastable conformations, and automated biophysical experiments using a novel fluorescence-based assay are used to directly measure binding affinities of kinase inhibitors to a variety of kinases. By analyzing this data together, we believe this work provides a more global picture of the role that kinase conformational reorganization energies play in selectivity, pushing us toward a better understanding of the requirements for selective inhibitor design.

Optical control of protein-protein interactions via blue light-induced domain swapping.

The design of new optogenetic tools for controlling protein function would be facilitated by the development of protein scaffolds that undergo large, well-defined structural changes upon exposure to light. Domain swapping, a process in which a structural element of a monomeric protein is replaced by the same element of another copy of the same protein, leads to a well-defined change in protein structure. We observe domain swapping in a variant of the blue light photoreceptor photoactive yellow protein in which a surface loop is replaced by a well-characterized protein–protein interaction motif, the E-helix. In the domain-swapped dimer, the E-helix sequence specifically binds a partner K-helix sequence, whereas in the monomeric form of the protein, the E-helix sequence is unable to fold into a binding-competent conformation and no interaction with the K-helix is seen. Blue light irradiation decreases the extent of domain swapping (from Kd = 10 μM to Kd = 300 μM) and dramatically enhances the rate, from weeks to <1 min. Blue light-induced domain swapping thus provides a novel mechanism for controlling of protein–protein interactions in which light alters both the stability and the kinetic accessibility of binding-competent states.

Co-Crystallization with Conformation-Specific Designed Ankyrin Repeat Proteins Explains the Conformational Flexibility of BCL-W

BCL-W is a member of the BCL-2 family of anti-apoptotic proteins. A key event in the regulation of apoptosis is the heterodimerization between anti-apoptotic and pro-apoptotic family members, which involves a conserved surface-exposed groove on the anti-apoptotic proteins. Crystal structures of the ligand binding-competent conformation exist for all anti-apoptotic family members, with the exception of BCL-W, due to the flexibility of the BCL-W groove region. Existing structures had suggested major deviations of the BCL-W groove region from the otherwise structurally highly related remaining anti-apoptotic family members. To capture its ligand binding-competent conformation by counteracting the conformational flexibility of the BCL-W groove, we had selected high-affinity groove-binding designed ankyrin repeat proteins (DARPins) using ribosome display. We now determined two high-resolution crystal structures of human BCL-W in complex with different DARPins at resolutions 1.5 and 1.85Å, in which the structure of BCL-W is virtually identical, and BCL-W adopts a conformation extremely similar to the ligand-free conformation of its closest relative BCL-XL in both structures. However, distinct differences to all previous BCL-W structures are evident, notably in the ligand-binding region. We provide the first structural explanation for the conformational flexibility of the BCL-W groove region in comparison to other BCL-2 family members. Due to the importance of the anti-apoptotic BCL-2 family as drug targets, the presented crystal structure of ligand binding-competent BCL-W may serve as a valuable basis for structure-based drug design in the future and provides a missing piece for the structural characterization of this protein family.