Facilitate Ligand Binding Research Articles

Chemical probing profiling of 2'-deoxyguanosine sensing RNA reveals magnesium ion-dependent conformational switching.

Maintaining the homeostasis of cellular metabolites is essential for all living organisms. Noncoding mRNA in procaryotes strategically regulates the expression of downstream genes which are linked to the expression of proteins involved in specific metabolite salvage and biosynthesis, in response to binding to that specific metabolite. Such noncoding mRNAs are known as the Riboswitches. A class of riboswitches that recognizes purine analogs were identified. In this present contribution, we report the influence of magnesium and temperature on the conformation of the aptamer domain (the ligand binding domain) of 2’-deoxyguanosine riboswitch from Mesoplasma florum. Our comprehensive SHAPE (Selective 2′ Hydroxyl Acylation analyzed by Primer Extension) probing experiments have shown that the aptamer domain exhibited a strong dependence on the Mg2+ ion concentration for the efficient ligand binding and the concomitant structural change. Mg2+ ions reduce the conformational flexibility mainly by organizing P2 and P3 helices and thereby facilitate ligand binding by shifting the conformational ensemble of the ligand-free aptamer towards a pre-organized, ligand-binding competent conformation. The subtle conformational switching in the presence and absence of 2’-deoxyguanosine is observed at the (A) regulatory helix (P1) and (B) the interhelical junctions (which serve as the ligand binding pocket) of the aptamer. We have observed the formation of a regulatory helix, the folding of the binding pocket, and the subsequent encapsulation of the ligand is accomplished only in the simultaneous presence of Mg2+ and 2’-deoxyguanosine. We have also identified the existence of two conformations in equilibrium in a ligand-free state. Finally, we have shown that for the 2’-dG sensing riboswitch from Mesoplasma forum, the conformation of the full-length transcripts populates the functional OFF-state regardless of the presence or absence of ligand, consistent with a kinetic riboswitch mechanism.

Mechanism of ligand binding signal transduction in a transcriptional riboswitch.

Riboswitches are noncoding RNA molecules that regulate gene expression in bacteria. The riboswitch's aptamer domain (AD) senses the cognate ligand and induces a structural change in the expression platform (EP) to signal on/off to the gene expression machinery. Despite extensive studies, the question of how the ligand binding signal is transmitted from the AD to EP resulting in the structural transitions in the EP remains unanswered. We studied the ligand binding signal transduction from AD to EP in the preQ1-sensing transcriptional riboswitch. In the absence of preQ1, both AD and EP individually form stem-loops with no structural overlap between the two stem-loops. The segment between the AD and EP stem loops is rich in adenine. In the presence of preQ1, AD forms a pseudoknot to encapsulate preQ1. We performed coarse-grained molecular dynamics simulations of full-length preQ1 riboswitch in the presence of Mg2+ and preQ1 ligand. In the absence of preQ1 and variation in Mg2+ concentration ([Mg2+]), we observed that AD populates only stem-loop and does not form the pseudoknot required to capture the ligand. Only in the presence of preQ1 AD attains the pseudoknot conformation by wrapping the A-rich region around the pre-existing stem-loop to facilitate ligand binding. The pseudoknot formation at AD allows EP to form the terminator helix and sequester the poly-U through strand invasion. The wrapping of the A-rich region around the AD stem-loop forces the stem-loop at EP to rupture and slide by 10 nucleotides towards 3'-end and form the terminator helix. The propensity of pseudoknot and terminator formation through strand invasion in the presence of preQ1 is enhanced with increment in [Mg2+]. We show that the competition between base-pair formation and helix sliding rate is crucial for the ligand binding signal transduction.

Dynamical characterization and multiple unbinding paths of two PreQ1 ligands in one pocket.

Riboswitches naturally regulate gene expression in bacteria by binding to specific small molecules. Class 1 preQ1 riboswitch aptamer is an important model not only for RNA folding but also as a target for designing small molecule antibiotics due to its well-known minimal aptamer domain. Here, we ran a total of 62.4 μs conventional and enhanced-sampling molecular dynamics (MD) simulations to characterize the determinants underlying the binding of the preQ1-II riboswitch aptamer to two preQ1 ligands in one binding pocket. Decomposition of binding free energy suggested that preQ1 ligands at α and β sites interact with four nucleotides (G5, C17, C18, and A30) and two nucleotides (A12 and C31), respectively. Mg2+ ions play a crucial role in both stabilizing the binding pocket and facilitating ligand binding. The flexible preQ1 ligand at the β site leads to the top of the binding pocket loosening and thus pre-organizes the riboswitch for ligand entry. Enhanced sampling simulations further revealed that the preQ1 ligand at the α site unbinds through two orthogonal pathways, which are dependent on whether or not a β site preQ1 ligand is present. One of the two preQ1 ligands has been identified in the binding pocket, which will aid to identify the second preQ1 Ligand. Our work provides new information for designing robust ligands.

Emerging role of exosomes in hematological malignancies.

Hematological malignancies are a heterogeneous group of neoplasms in the blood characterized by dysregulated hematopoiesis and classified as leukemia, lymphoma, and myeloma. The occurrence and progression of hematological malignancies depend on transformed hematopoietic stem cells, which refract to chemotherapy and often cause relapse. In recent years, monoclonal antibody therapies are preferred for hematopoietic cancers, owing to their inherent mechanisms of action and improved outcomes. However, efficient drug delivery methods and the establishment of novel biomarkers are currently being investigated and warranted to improve the outcome of patients with hematological malignancies. For instance, non-viral-mediated, natural carriers have been suggested for latent intracellular drug delivery. In this purview, repurposing small vesicles (e.g., exosomes) is considered a latent approach for myeloma therapy. Exosomes (nano-vesicles) have many advantages in that they are secreted by various animals and plants and become sought after for therapeutic and diagnostic purposes. The size of the cellular membrane of exosomes (30-150nm) facilitates ligand binding and targeted delivery of the loaded molecules. Furthermore, exosomes can be modified to express specific target moiety on their cell membrane and can also be featured with desired biological activity, thereby potentially employed for various convoluted diseases, including hematological malignancies. To advance the current knowledge, this review is focused on the source, composition, function and surface engineering of exosomes pertaining to hematological malignancies.

Allosteric modulation of the chemokine receptor-chemokine CXCR4-CXCL12 complex by tyrosine sulfation

The chemokine receptor CXCR4 and its cognate ligand CXCL12 mediate pathways that lead to cell migration and chemotaxis. Although the structural details of related receptor-ligand complexes have been resolved, the roles of the N-terminal domain of the receptor and post-translational sulfation that are determinants of ligand selectivity and affinity remain unclear. Here, we analyze the structural dynamics induced by receptor sulfation by combining molecular dynamics, docking and network analysis. The sulfotyrosine residues, 7YsN-term, 12YsN-term and 21YsN-term allow the N-terminal domain of the apo-sulfated receptor to adopt an “open” conformation that appears to facilitate ligand binding. The overall topology of the CXCR4-CXCL12 complex is independent of the sulfation state, but an extensive network of protein-protein interactions characterizes the sulfated receptor, in line with its increased ligand affinity. The altered interactions of sulfotyrosine residues, such as 21YsN-term-47RCXCL12 replacing the 21YN-term-13FCXCL12 interaction, propagate via allosteric pathways towards the receptor lumen. In particular, our results suggest that the experimentally-reported receptor-ligand interactions 262D6.58-8RCXCL12 and 277E7.28-12RCXCL12 could be dependent on the sulfation state of the receptor and need to be carefully analyzed. Our work is an important step in understanding chemokine-receptor interactions and how post-translational modifications could modulate receptor-ligand complexes.

Computational Cosolvent Mapping Analysis Leads to Identify Salicylic Acid Analogs as Weak Inhibitors of ST2 and IL33 Binding.

Cytokine signaling initiated by the binding of the cytokine receptors to cytokines plays important roles in immune regulation and diseases. Structurally, cytokine receptors interact with cytokines via an extensive, rugged interface that represents a challenge in inhibitor development. Our computational analysis has previously indicated that butyric acid, mimicking acidic residues, preferentially binds to sites in ST2 (Stimulation-2) that interact with acidic residues of IL33, the endogenous cytokine for ST2. To investigate if a charged group in small molecules facilitates ligand binding to ST2, we developed a biochemical homogeneous time resolved fluorescence assay to determine the inhibition of ST2/IL33 binding by five molecules containing an aromatic ring and a charged group. Three molecules, including niacin, salicylic acid, and benzamidine, exhibit inhibition activities at millimolar concentrations. We further employed the computational cosolvent mapping analysis to identify a shared mode of interaction between niacin, salicylic acid, and ST2. The mode of interaction was further confirmed by four analogous compounds that exhibited similar or improved activities. Our study provided the evidence of inhibition of ST2 and IL33 binding by salicylic acid and analogs. The results suggest that biological activity of salicylic acid may be partly mediated through modulating extracellular cytokine receptors and cytokine interaction.

Using ConSurf to Detect Functionally Important Regions in RNA.

The ConSurf web server (https://consurf.tau.ac.il/) for using evolutionary data to detect functional regions is useful for analyzing proteins. The analysis is based on the premise that functional regions, which may for example facilitate ligand binding and catalysis, often evolve slowly. The analysis requires finding enough effective, i.e., non-redundant, sufficiently remote homologs. Indeed, the ConSurf pipeline, which is based on state-of-the-art protein sequence databases and analysis tools, is highly valuable for protein analysis. ConSurf also allows evolutionary analysis of RNA, but the analysis often fails due to insufficient data, particularly the inability of the current pipeline to detect enough effective RNA homologs. This is because the RNA search tools and databases offered are not as good as those used for protein analysis. Fortunately, ConSurf also allows importing external collections of homologs in the form of a multiple sequence alignment (MSA). Leveraging this, here we describe various protocols for constructing MSAs for successful ConSurf analysis of RNA queries. We report the level of success of these protocols on an exemplary set comprising a dozen RNA molecules of diverse structure and function. © 2021 Wiley Periodicals LLC. Basic Protocol 1: Standard ConSurf evolutionary conservation analysis of an RNA query. Basic Protocol 2: ConSurf evolutionary conservation analysis of an RNA query with external MSA. Support Protocol 1: Construction of an MSA for an RNA query using other online servers. Support Protocol 2: Construction of an MSA for an RNA query using nHMMER locally.

The Oxoglutarate Binding Site and Regulatory Mechanism Are Conserved in Ammonium Transporter Inhibitors GlnKs from Methanococcales.

Protein inhibition is a natural regulatory process to control cellular metabolic fluxes. PII-family signal-transducing effectors are in this matter key regulators of the nitrogen metabolism. Their interaction with their various targets is governed by the cellular nitrogen level and the energy charge. Structural studies on GlnK, a PII-family inhibitor of the ammonium transporters (Amt), showed that the T-loops responsible for channel obstruction are displaced upon the binding of 2-oxoglutarate, magnesium and ATP in a conserved cleft. However, GlnK from Methanocaldococcus jannaschii was shown to bind 2-oxoglutarate on the tip of its T-loop, causing a moderate disruption to GlnK–Amt interaction, raising the question if methanogenic archaea use a singular adaptive strategy. Here we show that membrane fractions of Methanothermococcus thermolithotrophicus released GlnKs only in the presence of Mg-ATP and 2-oxoglutarate. This observation led us to structurally characterize the two GlnK isoforms apo or in complex with ligands. Together, our results show that the 2-oxoglutarate binding interface is conserved in GlnKs from Methanococcales, including Methanocaldococcus jannaschii, emphasizing the importance of a free carboxy-terminal group to facilitate ligand binding and to provoke the shift of the T-loop positions.

Evidence that GPVI is Expressed as a Mixture of Monomers and Dimers, and that the D2 Domain is not Essential for GPVI Activation.

Collagen has been proposed to bind to a unique epitope in dimeric glycoprotein VI (GPVI) and the number of GPVI dimers has been reported to increase upon platelet activation. However, in contrast, the crystal structure of GPVI in complex with collagen-related peptide (CRP) showed binding distinct from the site of dimerization. Further fibrinogen has been reported to bind to monomeric but not dimeric GPVI. In the present study, we have used the advanced fluorescence microscopy techniques of single-molecule microscopy, fluorescence correlation spectroscopy (FCS) and bioluminescence resonance energy transfer (BRET), and mutagenesis studies in a transfected cell line model to show that GPVI is expressed as a mixture of monomers and dimers and that dimerization through the D2 domain is not critical for activation. As many of these techniques cannot be applied to platelets to resolve this issue, due to the high density of GPVI and its anucleate nature, we used Förster resonance energy transfer (FRET) to show that endogenous GPVI is at least partially expressed as a dimer on resting and activated platelet membranes. We propose that GPVI may be expressed as a monomer on the cell surface and it forms dimers in the membrane through diffusion, giving rise to a mixture of monomers and dimers. We speculate that the formation of dimers facilitates ligand binding through avidity.

Ligand Gaussian Accelerated Molecular Dynamics (LiGaMD) for Characterization of Ligand Binding Thermodynamics and Kinetics

Calculations of ligand binding free energies and kinetic rates are important for drug design. However, such tasks have proven challenging in computational chemistry and biophysics. To address this challenge, we have developed a new computational method “LiGaMD”, which selectively boosts the ligand non-bonded interaction potential energy based on the Gaussian accelerated molecular dynamics (GaMD) enhanced sampling technique. Another boost potential could be applied to the remaining potential energy of the entire system in a dual-boost algorithm (LiGaMD_Dual) to facilitate ligand binding.

Evidence that immunization with TP0751, a bipartite Treponema pallidum lipoprotein with an intrinsically disordered region and lipocalin fold, fails to protect in the rabbit model of experimental syphilis.

Deconvolution of syphilis pathogenesis and selection of candidate syphilis vaccinogens requires detailed knowledge of the molecular architecture of the Treponema pallidum outer membrane (OM). The T. pallidum OM contains a low density of integral OM proteins, while the spirochete’s many lipoprotein immunogens are periplasmic. TP0751, a lipoprotein with a lipocalin fold, is reportedly a surface-exposed protease/adhesin and protective antigen. The rapid expansion of calycin/lipocalin structures in the RCSB PDB database prompted a comprehensive reassessment of TP0751. Small angle X-ray scattering analysis of full-length protein revealed a bipartite topology consisting of an N-terminal, intrinsically disordered region (IDR) and the previously characterized C-terminal lipocalin domain. A DALI server query using the lipocalin domain yielded 97 hits, 52 belonging to the calycin superfamily, including 15 bacterial lipocalins, but no Gram-negative surface proteins. Surprisingly, Tpp17 (TP0435) was identified as a structural ortholog of TP0751. In silico docking predicted that TP0751 can bind diverse ligands along the rim of its eight-stranded β-barrel; high affinity binding of one predicted ligand, heme, to the lipocalin domain was demonstrated. qRT-PCR and immunoblotting revealed very low expression of TP0751 compared to other T. pallidum lipoproteins. Immunoblot analysis of immune rabbit serum failed to detect TP0751 antibodies, while only one of five patients with secondary syphilis mounted a discernible TP0751-specific antibody response. In opsonophagocytosis assays, neither TP0751 nor Tpp17 antibodies promoted uptake of T. pallidum by rabbit peritoneal macrophages. Rabbits immunized with intact, full-length TP0751 showed no protection against local or disseminated infection following intradermal challenge with T. pallidum. Our data argue that, like other lipoprotein lipocalins in dual-membrane bacteria, TP0751 is periplasmic and binds small molecules, and we propose that its IDR facilitates ligand binding by and offloading from the lipocalin domain. The inability of TP0751 to elicit opsonic or protective antibodies is consistent with a subsurface location.

Ligand Gaussian Accelerated Molecular Dynamics (LiGaMD): Characterization of Ligand Binding Thermodynamics and Kinetics.

Calculations of ligand binding free energies and kinetic rates are important for drug design. However, such tasks have proven challenging in computational chemistry and biophysics. To address this challenge, we have developed a new computational method, ligand Gaussian accelerated molecular dynamics (LiGaMD), which selectively boosts the ligand nonbonded interaction potential energy based on the Gaussian accelerated molecular dynamics (GaMD) enhanced sampling technique. Another boost potential could be applied to the remaining potential energy of the entire system in a dual-boost algorithm (LiGaMD_Dual) to facilitate ligand binding. LiGaMD has been demonstrated on host-guest and protein-ligand binding model systems. Repetitive guest binding and unbinding in the β-cyclodextrin host were observed in hundreds-of-nanosecond LiGaMD_Dual simulations. The calculated guest binding free energies agreed excellently with experimental data with <1.0 kcal/mol errors. Compared with converged microsecond-time scale conventional molecular dynamics simulations, the sampling errors of LiGaMD_Dual simulations were also <1.0 kcal/mol. Accelerations of ligand kinetic rate constants in LiGaMD simulations were properly estimated using Kramers' rate theory. Furthermore, LiGaMD allowed us to capture repetitive dissociation and binding of the benzamidine inhibitor in trypsin within 1 μs simulations. The calculated ligand binding free energy and kinetic rate constants compared well with the experimental data. In summary, LiGaMD provides a powerful enhanced sampling approach for characterizing ligand binding thermodynamics and kinetics simultaneously, which is expected to facilitate computer-aided drug design.

The Observation of Ligand-Binding-Relevant Open States of Fatty Acid Binding Protein by Molecular Dynamics Simulations and a Markov State Model.

As a member of the fatty acids transporter family, the heart fatty acid binding proteins (HFABPs) are responsible for many important biological activities. The binding mechanism of fatty acid with FABP is critical to the understanding of FABP functions. The uncovering of binding-relevant intermediate states and interactions would greatly increase our knowledge of the binding process. In this work, all-atom molecular dynamics (MD) simulations were performed to characterize the structural properties of nativelike intermediate states. Based on multiple 6 μs MD simulations and Markov state model (MSM) analysis, several “open” intermediate states were observed. The transition rates between these states and the native closed state are in good agreement with the experimental measurements, which indicates that these intermediate states are binding relevant. As a common property in the open states, the partially unfolded α2 helix generates a larger portal and provides the driving force to facilitate ligand binding. On the other side, there are two kinds of open states for the ligand-binding HFABP: one has the partially unfolded α2 helix, and the other has the looser β-barrel with disjointing βD-βE strands. Our results provide atomic-level descriptions of the binding-relevant intermediate states and could improve our understanding of the binding mechanism.

Noncanonical role for the binding protein in substrate uptake by the MetNI methionine ATP Binding Cassette (ABC) transporter

The Escherichia coli methionine ABC transporter MetNI exhibits both high-affinity transport toward l-methionine and broad specificity toward methionine derivatives, including d-methionine. In this work, we characterize the transport of d-methionine derivatives by the MetNI transporter. Unexpectedly, the N229A substrate-binding deficient variant of the cognate binding protein MetQ was found to support high MetNI transport activity toward d-selenomethionine. We determined the crystal structure at 2.95 Å resolution of the ATPγS-bound MetNIQ complex in the outward-facing conformation with the N229A apo MetQ variant. This structure revealed conformational changes in MetQ providing substrate access through the binding protein to the transmembrane translocation pathway. MetQ likely mediates uptake of methionine derivatives through two mechanisms: in the methionine-bound form delivering substrate from the periplasm to the transporter (the canonical mechanism) and in the apo form by facilitating ligand binding when complexed to the transporter (the noncanonical mechanism). This dual role for substrate-binding proteins is proposed to provide a kinetic strategy for ABC transporters to transport both high- and low-affinity substrates present in a physiological concentration range.

A molecular dynamics investigation of CDK8/CycC and ligand binding: conformational flexibility and implication in drug discovery.

Abnormal activity of cyclin-dependent kinase 8 (CDK8) along with its partner protein cyclin C (CycC) is a common feature of many diseases including colorectal cancer. Using molecular dynamics (MD) simulations, this study determined the dynamics of the CDK8-CycC system and we obtained detailed breakdowns of binding energy contributions for four type-I and five type-II CDK8 inhibitors. We revealed system motions and conformational changes that will affect ligand binding, confirmed the essentialness of CycC for inclusion in future computational studies, and provide guidance in development of CDK8 binders. We employed unbiased all-atom MD simulations for 500ns on twelve CDK8-CycC systems, including apoproteins and protein-ligand complexes, then performed principal component analysis (PCA) and measured the RMSF of key regions to identify protein dynamics. Binding pocket volume analysis identified conformational changes that accompany ligand binding. Next, H-bond analysis, residue-wise interaction calculations, and MM/PBSA were performed to characterize protein-ligand interactions and find the binding energy. We discovered that CycC is vital for maintaining a proper conformation of CDK8 to facilitate ligand binding and that the system exhibits motion that should be carefully considered in future computational work. Surprisingly, we found that motion of the activation loop did not affect ligand binding. Type-I and type-II ligand binding is driven by van der Waals interactions, but electrostatic energy and entropic penalties affect type-II binding as well. Binding of both ligand types affects protein flexibility. Based on this we provide suggestions for development of tighter-binding CDK8 inhibitors and offer insight that can aid future computational studies.

Interplay of self-association and conformational flexibility in regulating protein function.

Many functional roles have been attributed to homodimers, the most common mode of protein self-association, notably in the regulation of enzymes, ion channels, transporters and transcription factors. Here we review findings that offer new insights into the different roles conformational flexibility plays in regulating homodimer function. Intertwined homodimers of two-domain proteins and their related family members display significant conformational flexibility, which translates into concerted motion between structural domains. This flexibility enables the corresponding proteins to regulate function across family members by modulating the spatial positions of key recognition surfaces of individual domains, to either maintain subunit interfaces, alter them or break them altogether, leading to a variety of functional consequences. Many proteins may exist as monomers but carry out their biological function as homodimers or higher-order oligomers. We present early evidence that in such systems homodimer formation primes the protein for its functional role. It does so by inducing elevated mobility in protein regions corresponding to the binding epitopes of functionally important ligands. In some systems this process acts as an allosteric response elicited by the self-association reaction itself. Our analysis furthermore suggests that the induced extra mobility likely facilitates ligand binding through the mechanism of conformational selection.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

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.

Kinetic Dissection of the Pre-Existing Conformational Equilibrium in the Trypsin Fold

Structural biology has recently documented the conformational plasticity of the trypsin fold for both the protease and zymogen in terms of a pre-existing equilibrium between closed (E∗) and open (E) forms of the active site region. How such plasticity is manifested in solution and affects ligand recognition by the protease and zymogen is poorly understood in quantitative terms. Here we dissect the E∗-E equilibrium with stopped-flow kinetics in the presence of excess ligand or macromolecule. Using the clotting protease thrombin and its zymogen precursor prethrombin-2 as relevant models we resolve the relative distribution of the E∗ and E forms and the underlying kinetic rates for their interconversion. In the case of thrombin, the E∗ and E forms are distributed in a 1:4 ratio and interconvert on a time scale of 45 ms. In the case of prethrombin-2, the equilibrium is shifted strongly (10:1 ratio) in favor of the closed E∗ form and unfolds over a faster time scale of 4.5 ms. The distribution of E∗ and E forms observed for thrombin and prethrombin-2 indicates that zymogen activation is linked to a significant shift in the pre-existing equilibrium between closed and open conformations that facilitates ligand binding to the active site. These findings broaden our mechanistic understanding of how conformational transitions control ligand recognition by thrombin and its zymogen precursor prethrombin-2 and have direct relevance to other members of the trypsin fold.

Kinetic Dissection of the Pre-existing Conformational Equilibrium in the Trypsin Fold

Structural biology has recently documented the conformational plasticity of the trypsin fold for both the protease and zymogen in terms of a pre-existing equilibrium between closed (E*) and open (E) forms of the active site region. How such plasticity is manifested in solution and affects ligand recognition by the protease and zymogen is poorly understood in quantitative terms. Here we dissect the E*-E equilibrium with stopped-flow kinetics in the presence of excess ligand or macromolecule. Using the clotting protease thrombin and its zymogen precursor prethrombin-2 as relevant models we resolve the relative distribution of the E* and E forms and the underlying kinetic rates for their interconversion. In the case of thrombin, the E* and E forms are distributed in a 1:4 ratio and interconvert on a time scale of 45 ms. In the case of prethrombin-2, the equilibrium is shifted strongly (10:1 ratio) in favor of the closed E* form and unfolds over a faster time scale of 4.5 ms. The distribution of E* and E forms observed for thrombin and prethrombin-2 indicates that zymogen activation is linked to a significant shift in the pre-existing equilibrium between closed and open conformations that facilitates ligand binding to the active site. These findings broaden our mechanistic understanding of how conformational transitions control ligand recognition by thrombin and its zymogen precursor prethrombin-2 and have direct relevance to other members of the trypsin fold.

Extended molecular dynamics of a c-kit promoter quadruplex.

The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.