Multiple Binding Modes Research Articles

Characterizing protein-DNA binding event subtypes in ChIP-exo data.

Regulatory proteins associate with the genome either by directly binding cognate DNA motifs or via protein-protein interactions with other regulators. Each recruitment mechanism may be associated with distinct motifs and may also result in distinct characteristic patterns in high-resolution protein-DNA binding assays. For example, the ChIP-exo protocol precisely characterizes protein-DNA crosslinking patterns by combining chromatin immunoprecipitation (ChIP) with 5' → 3' exonuclease digestion. Since different regulatory complexes will result in different protein-DNA crosslinking signatures, analysis of ChIP-exo tag enrichment patterns should enable detection of multiple protein-DNA binding modes for a given regulatory protein. However, current ChIP-exo analysis methods either treat all binding events as being of a uniform type or rely on motifs to cluster binding events into subtypes. To systematically detect multiple protein-DNA interaction modes in a single ChIP-exo experiment, we introduce the ChIP-exo mixture model (ChExMix). ChExMix probabilistically models the genomic locations and subtype memberships of binding events using both ChIP-exo tag distribution patterns and DNA motifs. We demonstrate that ChExMix achieves accurate detection and classification of binding event subtypes using in silico mixed ChIP-exo data. We further demonstrate the unique analysis abilities of ChExMix using a collection of ChIP-exo experiments that profile the binding of key transcription factors in MCF-7 cells. In these data, ChExMix identifies possible recruitment mechanisms of FoxA1 and ERα, thus demonstrating that ChExMix can effectively stratify ChIP-exo binding events into biologically meaningful subtypes. ChExMix is available from https://github.com/seqcode/chexmix. Supplementary data are available at Bioinformatics online.

Evolutionary repurposing of a sulfatase: A new Michaelis complex leads to efficient transition state charge offset

The recruitment and evolutionary optimization of promiscuous enzymes is key to the rapid adaptation of organisms to changing environments. Our understanding of the precise mechanisms underlying enzyme repurposing is, however, limited: What are the active-site features that enable the molecular recognition of multiple substrates with contrasting catalytic requirements? To gain insights into the molecular determinants of adaptation in promiscuous enzymes, we performed the laboratory evolution of an arylsulfatase to improve its initially weak phenylphosphonate hydrolase activity. The evolutionary trajectory led to a 100,000-fold enhancement of phenylphosphonate hydrolysis, while the native sulfate and promiscuous phosphate mono- and diester hydrolyses were only marginally affected (≤50-fold). Structural, kinetic, and in silico characterizations of the evolutionary intermediates revealed that two key mutations, T50A and M72V, locally reshaped the active site, improving access to the catalytic machinery for the phosphonate. Measured transition state (TS) charge changes along the trajectory suggest the creation of a new Michaelis complex (E•S, enzyme-substrate), with enhanced leaving group stabilization in the TS for the promiscuous phosphonate (βleavinggroup from -1.08 to -0.42). Rather than altering the catalytic machinery, evolutionary repurposing was achieved by fine-tuning the molecular recognition of the phosphonate in the Michaelis complex, and by extension, also in the TS. This molecular scenario constitutes a mechanistic alternative to adaptation solely based on enzyme flexibility and conformational selection. Instead, rapid functional transitions between distinct chemical reactions rely on the high reactivity of permissive active-site architectures that allow multiple substrate binding modes.

Elucidation of the orientation of selected drugs with 2-hydroxylpropyl-β-cyclodextrin using 2D-NMR spectroscopy and molecular modeling

This project aims to study the nature of interaction and orientation of selected drugs such as dexamethorphan HBr (DXM), diphenhydramine HCl (DPH), and lidocaine HCl (LDC) inclusion complexes with hydroxyl-propyl ß-cyclodextrin (HP-ß-CD) using 1HNMR spectroscopy, 2D-NMR ROESY and molecular-modeling techniques. Freeze-drying technique was used to formulate the inclusion complexes between DXM, DPH and LDC with HP-ß-CD (1:1 M ratio) in solid state. Inclusion complex formation was initially characterized by Fourier transform-infrared spectroscopy (FT-IR), differential scanning calorimetry (DSC), X-ray diffraction (XRD) and scanning electron microscopy (SEM) techniques. Further characterization of inclusion complexes to determine the interaction of DXM, DPH and LDC with HP-β-CD was performed using the 1HNMR spectroscopy, 2D-NMR ROESY and molecular modeling techniques. Inclusion complexes of DXM, DPH and LDC with HP-ß-CD were successfully prepared using the freeze-drying technique. Preliminary studies with FT-IR, DSC, XRD and SEM indicated the formation of inclusion complexes of DXM, DPH and LDC with HP-β-CD at 1:1 M ratio. 1HNMR study showed a change in proton chemical shift upon complexation. 2D-NMR ROESY (two-dimensional) spectroscopy gave an insight into the spatial arrangement between the host and guest atoms. 2D-ROESY experiments further predicted the direction of orientation of guest molecules, indicating the probability that amino moieties of DXM, DPH and LDC are inside the hydrophobic HP-ß-CD cavity. Cross-peaks of inclusion complexes demonstrated intermolecular nuclear Overhauser effects (NOE) between the amino protons in DXM, DPH and LDC and H-atoms of HP-ß-CD. Molecular modeling studies further confirmed the NMR data, providing a structural basis of the individual complex formations. Microsecond time-level molecular dynamics and metadynamics simulations indicate much stronger binding of DXM to HP-ß-CD and more dynamic behavior for DPH and LDC. In particular, LDC can exhibit multiple binding modes, and even spent some time (∼1–2%) out of the carrier, proving the dynamic nature of the complex. To conclude, 2D-NMR and molecular dynamic simulations elucidate the formation of inclusion complexes and intermolecular interactions of DXM, DPH and LDC with HP-ß-CD.

Prediction of Accurate Binding Modes Using Combination of Classical and Accelerated Molecular Dynamics and Free-Energy Perturbation Calculations: An Application to Toxicity Studies.

Estimating the correct binding modes of ligands in protein–ligand complexes is crucial not only in the drug discovery process but also for elucidating potential toxicity mechanisms. In the current paper, we propose a computational modeling workflow using the combination of docking, classical molecular dynamics (cMD), accelerated molecular dynamics (aMD) and free-energy perturbation (FEP+ protocol) for identification of possible ligand binding modes. It was applied for investigation of selected perfluorocarboxyl acids (PFCAs) in the PPARγ nuclear receptor. Although both regular and induced fit docking failed to reproduce the experimentally determined binding mode of the ligands when docked into a non-native X-ray structure, cMD and aMD simulations successfully identified the most probable binding conformations. Moreover, multiple binding modes were identified for all of these compounds and the shorter-chain PFCAs continuously moved between a few energetically favorable binding conformations. On the basis of MD predictions of binding conformations, we applied the default and also redesigned FEP+ sampling protocols, which accurately reproduced experimental differences in the binding energies. Thus, the preliminary MD simulations can also provide helpful information about correct setup of the FEP+ calculations. These results show that the PFCA binding modes were accurately predicted and that the FEP+ protocol can be used to estimate free energies of binding of flexible ligands that are not typical druglike compounds. Our in silico workflow revealed the specific ligand–residue interactions within the ligand binding domain and the main characteristics of the PFCAs, and it was concluded that these compounds are week PPARγ partial agonists. This work also suggests a common pipeline for identification of ligand binding modes, ligand–protein dynamics description, and relative free-energy calculations.

Accurate and sensitive quantification of protein-DNA binding affinity

Transcription factors (TFs) control gene expression by binding to genomic DNA in a sequence-specific manner. Mutations in TF binding sites are increasingly found to be associated with human disease, yet we currently lack robust methods to predict these sites. Here, we developed a versatile maximum likelihood framework named No Read Left Behind (NRLB) that infers a biophysical model of protein-DNA recognition across the full affinity range from a library of in vitro selected DNA binding sites. NRLB predicts human Max homodimer binding in near-perfect agreement with existing low-throughput measurements. It can capture the specificity of the p53 tetramer and distinguish multiple binding modes within a single sample. Additionally, we confirm that newly identified low-affinity enhancer binding sites are functional in vivo, and that their contribution to gene expression matches their predicted affinity. Our results establish a powerful paradigm for identifying protein binding sites and interpreting gene regulatory sequences in eukaryotic genomes.

Mutational analysis of the extracellular disulphide bridges of the atypical chemokine receptor ACKR3/CXCR7 uncovers multiple binding and activation modes for its chemokine and endogenous non-chemokine agonists

The atypical chemokine receptor ACKR3/CXCR7 plays crucial roles in numerous physiological processes but also in viral infection and cancer. ACKR3 shows strong propensity for activation and, unlike classical chemokine receptors, can respond to chemokines from both the CXC and CC families as well as to the endogenous peptides BAM22 and adrenomedullin. Moreover, despite belonging to the G protein coupled receptor family, its function appears to be mainly dependent on β-arrestin. ACKR3 has also been shown to continuously cycle between the plasma membrane and the endosomal compartments, suggesting a possible role as a scavenging receptor. So far, the molecular basis accounting for these atypical binding and signalling properties remains elusive. Noteworthy, ACKR3 extracellular domains bear three disulphide bridges. Two of them lie on top of the two main binding subpockets and are conserved among chemokine receptors, and one, specific to ACKR3, forms an intra-N terminus four-residue-loop of so far unknown function. Here, by mutational and functional studies, we examined the impact of the different disulphide bridges for ACKR3 folding, ligand binding and activation. We showed that, in contrast to most classical chemokine receptors, none of the extracellular disulphide bridges was essential for ACKR3 function. However, the disruption of the unique ACKR3 N-terminal loop drastically reduced the binding of CC chemokines whereas it only had a mild impact on CXC chemokine binding. Mutagenesis also uncovered that chemokine and endogenous non-chemokine ligands interact and activate ACKR3 according to distinct binding modes characterized by different transmembrane domain subpocket occupancy and N-terminal loop contribution, with BAM22 mimicking the binding mode of CC chemokine N terminus.

Binding Modes of Ligands Using Enhanced Sampling (BLUES): Rapid Decorrelation of Ligand Binding Modes via Nonequilibrium Candidate Monte Carlo.

Accurately predicting protein-ligand binding affinities and binding modes is a major goal in computational chemistry, but even the prediction of ligand binding modes in proteins poses major challenges. Here, we focus on solving the binding mode prediction problem for rigid fragments. That is, we focus on computing the dominant placement, conformation, and orientations of a relatively rigid, fragment-like ligand in a receptor, and the populations of the multiple binding modes which may be relevant. This problem is important in its own right, but is even more timely given the recent success of alchemical free energy calculations. Alchemical calculations are increasingly used to predict binding free energies of ligands to receptors. However, the accuracy of these calculations is dependent on proper sampling of the relevant ligand binding modes. Unfortunately, ligand binding modes may often be uncertain, hard to predict, and/or slow to interconvert on simulation time scales, so proper sampling with current techniques can require prohibitively long simulations. We need new methods which dramatically improve sampling of ligand binding modes. Here, we develop and apply a nonequilibrium candidate Monte Carlo (NCMC) method to improve sampling of ligand binding modes. In this technique, the ligand is rotated and subsequently allowed to relax in its new position through alchemical perturbation before accepting or rejecting the rotation and relaxation as a nonequilibrium Monte Carlo move. When applied to a T4 lysozyme model binding system, this NCMC method shows over 2 orders of magnitude improvement in binding mode sampling efficiency compared to a brute force molecular dynamics simulation. This is a first step toward applying this methodology to pharmaceutically relevant binding of fragments and, eventually, drug-like molecules. We are making this approach available via our new Binding modes of ligands using enhanced sampling (BLUES) package which is freely available on GitHub.

Probing Multiple Binding Modes of DNA Hybridization: A Comparison between Single-Molecule Observations and Ensemble Measurements.

Interactions between biomolecules are generally analyzed by ensemble measurements, assuming that the interactions occur in a single binding manner. However, such interactions may occur via multiple binding modes. We investigated the kinetics of DNA hybridization as a multiple dynamic model of biomolecular interactions. Two kinetic analyses were performed with a single-molecule observation using total internal reflection fluorescence microscopy (TIRFM) and with ensemble measurements using a quartz-crystal microbalance (QCM) biosensor. We observed the DNA hybridization of 8 and 12 bp DNAs with random sequences and dA12–dT12 and calculated the kinetic parameters, including the dissociation rate constant (koff). Hybridization of 8 bp DNA proceeded mainly via a single binding mode. However, hybridization of 12 bp DNA indicated at least two different binding modes and dA12–dT12 hybridization showed multiple binding modes. For the multiple binding interactions, the kinetic parameters obtained from TIRFM and QCM were different because kinetic parameters obtained from QCM indicate average number of molecules, whereas those from TIRFM indicate average association time. The present study revealed the details of multiple interactions, which can be utilized for better understanding of not only DNA hybridization but also biomolecular interaction mechanisms.

Fluorescent single-stranded DNA-binding protein from Plasmodium falciparum as a biosensor for single-stranded DNA.

Single-stranded DNA (ssDNA) is a product of many cellular processes that involve double-stranded DNA, for example during DNA replication and repair, and is formed transiently in many others. Measurement of ssDNA formation is fundamental for understanding many such processes. The availability of a fluorescent biosensor for the determination of single-stranded DNA provides an important route to achieve this. Single-stranded DNA binding proteins (SSBs) protect ssDNA from degradation, but can be displaced to allow processing of the ssDNA. Their tight binding of ssDNA means that they are very good candidates for the development of a biosensor. Previously, the single stranded DNA binding protein from Escherichia coli, labeled with a fluorophore, (DCC-EcSSB) was developed and used for this purpose. However, the multiple binding modes of this protein meant that interpretation of DCC-EcSSB fluorescence was potentially complex in terms of determining the amount of ssDNA. Here, we present an improved biosensor, developed using the tetrameric SSB from Plasmodium falciparum as a new scaffold for fluorophore attachment. Each subunit of this tetrameric SSB was labeled with a diethylaminocoumarin fluorophore at a single site on its surface, such that there is a very large, 20-fold, fluorescence increase when it binds to ssDNA. This adduct can be used as a biosensor to report ssDNA formation. Because SSB from this organism has a single mode of binding ssDNA, namely 65–70 bases per tetramer, over a wide range of conditions, the fluorescent SSB allows simple quantitation of ssDNA. The binding is fast, possibly diffusion controlled, and tight (dissociation constant for DCC-PfSSB <5 pM). Its suitability for real-time assays of ssDNA formation was demonstrated by measurement of AddAB helicase activity, unwinding double-stranded DNA.

Multiple binding modes of a moderate ice-binding protein from a polar microalga.

Ice-binding proteins (IBPs) produced by cold-tolerant organisms interact with ice and strongly control crystal growth. The molecular basis for the different magnitudes of activity displayed by various IBPs (moderate and hyperactive) has not yet been clarified. Previous studies questioned whether the moderate activity of some IBPs relies on their weaker binding modus to the ice surface, compared to hyperactive IBPs, rather than relying on binding only to selected faces of the ice crystal. We present the structure of one moderate IBP from the sea-ice diatom Fragilariopsis cylindrus (fcIBP) as determined by X-ray crystallography and investigate the protein's binding modes to the growing ice-water interface using molecular dynamics simulations. The structure of fcIBP is the IBP-1 fold, defined by a discontinuous β-solenoid delimitated by three faces (A, B and C-faces) and braced by an α-helix. The fcIBP structure shows capping loops on both N- and C-terminal parts of the solenoid. We show that the protein adsorbs on both the prism and the basal faces of ice crystals, confirming experimental results. The fcIBP binds irreversibly to the prism face using the loop between the B and the C-faces, involving also the B-face in water immobilization despite its irregular structure. The α-helix attaches the protein to the basal face with a partly reversible modus. Our results suggest that fcIBP has a looser attachment to ice and that this weaker binding modus is the basis to explain the moderate activity of fcIBP.

Multiple binding modes of an unconjugated bis(pyridine) ligand stabilize low-valent [Cp*Rh] complexes.

The ligand 2,2'-bipyridine (bpy) can support metal centers in low formal oxidation states by delocalization of electron density into its π-system. We show that, in a model rhodium complex supported by the pentamethylcyclopentadienyl ligand (Cp*), the analogous dimethyldipyridylmethane ligand (Me2dpma) enforces a bpy-like coordination environment but disrupts the inter-ring conjugation responsible for charge delocalization upon metal reduction. As a result, reduction proceeds in discrete one-electron steps (Rh(iii) to Rh(ii) to Rh(i)), contrasting with the 2e- chemistry engendered by bpy. Upon reduction to Rh(i), the Me2dpma ligand rearranges to activate strong π-backbonding via facial coordination of one pyridine motif. Structural and spectroscopic studies confirm stabilization of the Rh(i) center in this compound, revealing a mode of metal-ligand cooperation that represents a useful counterpoint to charge delocalization in conjugated poly(pyridyl) ligands.

Computational Design of Antiangiogenic Peptibody by Fusing Human IgG1 Fc Fragment and HRH Peptide: Structural Modeling, Energetic Analysis, and Dynamics Simulation of Its Binding Potency to VEGF Receptor.

Peptibodies represent a new class of biological therapeutics with combination of peptide activity and antibody-like properties. Previously, we discovered a novel peptide HRH that exhibited a dose-dependent angiogenesis-suppressing effect by targeting vascular endothelial growth factor receptors (VEGFRs). Here, we computationally designed an antiangiogenic peptibody, termed as PbHRH, by fusing the HRH peptide to human IgG1 Fc fragment using the first approved peptibody drug Romiplostim as template. The biologically active peptide of Romiplostim is similar with HRH peptide; both of them have close sequence lengths and can fold into a α-helical conformation in free state. Molecular dynamics simulations revealed that the HRH functional domain is highly flexible, which is functionally independent of Fc fragment in the designed PbHRH peptibody. Subsequently, the intermolecular interactions between VEGFR-1 domain 2 (D2) and PbHRH were predicted, clustered and refined into three representatives. Conformational analysis and energetic evaluation unraveled that the PbHRH can adopt multiple binding modes to block the native VEGF-A binding site of VEGFR-1 D2 with its HRH functional domain, although the binding effectiveness of HRH segments in peptibody context seems to be moderately decreased relative to that of free HRH peptide. Overall, it is suggested that integrating HRH peptide into PbHRH peptibody does not promote the direct intermolecular interaction between VEGFR-1 D2 and HRH. Instead, the peptibody may indirectly help to improve the pharmacokinetic profile and bioavailability of HRH.

X-ray Absorption Spectroscopy Combined with Time-Dependent Density Functional Theory Elucidates Differential Substitution Pathways of Au(I) and Au(III) with Zinc Fingers.

A combination of two elements' (Au, Zn) X-ray absorption spectroscopy (XAS) and time-dependent density functional theory (TD-DFT) allowed the elucidation of differential substitution pathways of Au(I) and Au(III) compounds reacting with biologically relevant zinc fingers (ZnFs). Gold L3-edge XAS probed the interaction of gold and the C-terminal Cys2HisCys finger of the HIV-1 nucleocapsid protein NCp7, and the Cys2His2 human transcription factor Sp1. The use of model compounds helped assign oxidation states and the identity of the gold-bound ligands. The computational studies accurately reproduced the experimental XAS spectra and allowed the proposition of structural models for the interaction products at early time points. The direct electrophilic attack on the ZnF by the highly thiophilic Au(I) resulted in a linear P-Au-Cys coordination sphere after zinc ejection whereas for the Sp1, loss of PEt3 results in linear Cys-Au-Cys or Cys-Au-His arrangements. Reactions with Au(III) compounds, on the other hand, showed multiple binding modes. Prompt reaction between [AuCl(dien)]2+ and [Au(dien)(DMAP)]3+ with Sp1 showed a partially reduced Au center and a final linear His-Au-His coordination. Differently, in the presence of NCp7, [AuCl(dien)]2+ readily reduces to Au(I) and changes from square-planar to linear geometry with Cys-Au-His coordination, while [Au(dien)(DMAP)]3+ initially maintains its Au(III) oxidation state and square-planar geometry and the same first coordination sphere. The latter is the first observation of a "noncovalent" interaction of a Au(III) complex with a zinc finger and confirms early hypotheses that stabilization of Au(III) occurs with N-donor ligands. Modification of the zinc coordination sphere, suggesting full or partial zinc ejection, is observed in all cases, and for [Au(dien)(DMAP)]3+ this represents a novel mechanism for nucleocapsid inactivation. The combination of XAS and TD-DFT presents the first direct experimental observation that not only compound reactivity, but also ZnF core specificity, can be modulated on the basis of the coordination sphere of Au(III) compounds.

Sampling and energy evaluation challenges in ligand binding protein design.

The steroid hormone 17α‐hydroxylprogesterone (17‐OHP) is a biomarker for congenital adrenal hyperplasia and hence there is considerable interest in development of sensors for this compound. We used computational protein design to generate protein models with binding sites for 17‐OHP containing an extended, nonpolar, shape‐complementary binding pocket for the four‐ring core of the compound, and hydrogen bonding residues at the base of the pocket to interact with carbonyl and hydroxyl groups at the more polar end of the ligand. Eight of 16 designed proteins experimentally tested bind 17‐OHP with micromolar affinity. A co‐crystal structure of one of the designs revealed that 17‐OHP is rotated 180° around a pseudo‐two‐fold axis in the compound and displays multiple binding modes within the pocket, while still interacting with all of the designed residues in the engineered site. Subsequent rounds of mutagenesis and binding selection improved the ligand affinity to nanomolar range, while appearing to constrain the ligand to a single bound conformation that maintains the same “flipped” orientation relative to the original design. We trace the discrepancy in the design calculations to two sources: first, a failure to model subtle backbone changes which alter the distribution of sidechain rotameric states and second, an underestimation of the energetic cost of desolvating the carbonyl and hydroxyl groups of the ligand. The difference between design model and crystal structure thus arises from both sampling limitations and energy function inaccuracies that are exacerbated by the near two‐fold symmetry of the molecule.

Fragment optimization for GPCRs by molecular dynamics free energy calculations: Probing druggable subpockets of the A2A adenosine receptor binding site

Fragment-based lead discovery is becoming an increasingly popular strategy for drug discovery. Fragment screening identifies weakly binding compounds that require optimization to become high-affinity leads. As design of leads from fragments is challenging, reliable computational methods to guide optimization would be invaluable. We evaluated using molecular dynamics simulations and the free energy perturbation method (MD/FEP) in fragment optimization for the A2A adenosine receptor, a pharmaceutically relevant G protein-coupled receptor. Optimization of fragments exploring two binding site subpockets was probed by calculating relative binding affinities for 23 adenine derivatives, resulting in strong agreement with experimental data (R2 = 0.78). The predictive power of MD/FEP was significantly better than that of an empirical scoring function. We also demonstrated the potential of the MD/FEP to assess multiple binding modes and to tailor the thermodynamic profile of ligands during optimization. Finally, MD/FEP was applied prospectively to optimize three nonpurine fragments, and predictions for 12 compounds were evaluated experimentally. The direction of the change in binding affinity was correctly predicted in a majority of the cases, and agreement with experiment could be improved with rigorous parameter derivation. The results suggest that MD/FEP will become a powerful tool in structure-driven optimization of fragments to lead candidates.

Endocrine-disrupting activity of per- and polyfluoroalkyl substances: Exploring combined approaches of ligand and structure based modeling

Exposure to perfluorinated and polyfluoroalkyl substances (PFCs/PFASs), endocrine disrupting halogenated pollutants, has been linked to various diseases including thyroid toxicity in human populations across the globe. PFASs can compete with thyroxine (T4) for binding to the human thyroid hormone transport protein transthyretin (TTR) which may lead to reduce thyroid hormone levels leading to endocrine disrupting adverse effects. Environmental fate and endocrine-disrupting activity of PFASs has initiated several research projects, but the amount of experimental data available for these pollutants is limited. In this study, experimental data for T4-TTR competing potency of 24 PFASs obtained in a radioligand-binding assay were modeled using classification- and regression-based quantitative structure-activity relationship (QSAR) tools with simple molecular descriptors obtained from chemical structure of these compounds in order to identify the responsible structural features and fragments of the studied PFASs for endocrine disruption activity. Additionally, docking studies were performed employing the crystal structure complex of TTR with bound 2′, 6′-difluorobiphenyl-4-carboxylic acid (PDB: 2F7I) in order to constitute the receptor model for human TTR. The results corroborate evidence for these binding interactions and indicate multiple high-affinity modes of binding. The developed in silico models therefore advance our understanding of important structural attributes of these chemicals and may provide important information for the design of future synthesis of PFASs as well as may serve as an efficient query tool for virtual screening of large PFAS databases to check their endocrine toxicity profile.

Identification of Transthyretin Fibril Formation Inhibitors Using Structure-Based Virtual Screening.

Transthyretin (TTR) is the primary carrier for thyroxine (T4 ) in cerebrospinal fluid and a secondary carrier in blood. TTR is a stable homotetramer, but certain factors, genetic or environmental, could promote its degradation to form amyloid fibrils. A docking study using crystal structures of wild-type TTR was planned; our aim was to design new ligands that are able to inhibit TTR fibril formation. The computational protocol was thought to overcome the multiple binding modes of the ligands induced by the peculiarity of the TTR binding site and by the pseudosymmetry of the site pockets, which generally weaken such structure-based studies. Two docking steps, one that is very fast and a subsequent step that is more accurate, were used to screen the Aldrich Market Select database. Five compounds were selected, and their activity toward inhibiting TTR fibril formation was assessed. Three compounds were observed to be actives, two of which have the same potency as the positive control, and the other was found to be a promising lead compound. These results validate a computational protocol that is able to archive information on the key interactions between database compounds and TTR, which is valuable for supporting further studies.

Binding Mechanism between CrkII and cAbl Kinase

The interaction between non-receptor tyrosine kinase cAbl and signaling adaptor protein CrkII is implicated in diverse cellular processes, such as cell migration, proliferation, and bacterial infection. Because of the modular architecture of their structures, the interaction between cAbl and CrkII can undergo multiple alternative binding modes. This series of interactions starts with the binding of the N-terminal Src homology 3 (nSH3) domain of CrkII to proline-rich motifs (PRMs) in the intrinsically disordered region of cAbl kinase. Despite their importance, the detailed binding mechanism and functional significance of these interactions remains elusive. In this presentation, we report the detailed structure, thermodynamics, and kinetics of the interactions between the nSH3 domain of CrkII and PRMs of cAbl kinase. Our dynamics and structural studies, using NMR spectroscopy and crystallography, indicate that these modular interactions act as a binding catalyst that accelerates the association between CrkII and cAbl kinase. We also present our recent finding of a potent inhibitor of the cAbl-CrkII interaction.

Copper binding and redox chemistry of the Aβ16 peptide and its variants: insights into determinants of copper-dependent reactivity.

The metal-binding sites of Aβ peptides are dictated primarily by the coordination preferences of the metal ion. Consequently, Cu(i) is typically bound with two His ligands in a linear mode while Cu(ii) forms a pseudo-square planar stereochemistry with the N-terminal amine nitrogen acting as an anchoring ligand. Several distinct combinations of other groups can act as co-ligands for Cu(ii). A population of multiple binding modes is possible with the equilibrium position shifting sensitively with solution pH and the nature of the residues in the N-terminal region. This work examined the Cu(ii) chemistry of the Aβ16 peptide and several variants that targeted these binding modes. The results are consistent with: (i) at pH < 7.8, the square planar site in CuII-Aβ16 consists primarily of a bidentate ligand provided by the carboxylate sidechain of Asp1 and the N-terminal amine supported by the imidazole sidechains of two His residues (designated here as component IA); it is in equilibrium with a less stable component IB in which the carboxylate ligand is substituted by the Asp1-Ala2 carbonyl oxygen. (ii) Both IA and IB convert to a common component II (apparent transition pKa ∼7.8 for IA and ∼6.5 for IB, respectively) featuring a tridentate ligand consisting of the N-terminal amine, the Asp1-Ala2 amide and the Ala2-Pro3 carbonyl; this stereochemistry is stabilized by two five-membered chelate rings. (iii) Component IA is stabilized for variant Aβ16-D1H, components I (both IA and IB) are imposed on Aβ16-A2P while the less stable IB is enforced on Aβ16-D1A (which is converted to component II at pH ∼6.5); (iv) components IA and IB share two His ligands with Cu(i) and are more reactive in redox catalysis than component II that features a highly covalent and less reactive amide N- ligand. The redox activity of IA is further enhanced for peptides with a His1 N-terminus that may act as a ligand for either Cu(i) or Cu(ii) with lower re-organization energy required for redox-shuttling. This study provided insights into the determinants that regulate the reactivity of Cu-Aβ complexes.

Two-Dimensional Crowding Uncovers a Hidden Conformation of α-Synuclein.

The intrinsically disordered protein (IDP), α-synuclein (αS), is well-known for phospholipid membrane binding-coupled folding into tunable helical conformers. Here, using single-molecule experiments in conjunction with ensemble assays and a theoretical model, we present a unique case demonstrating that the interaction-folding landscape of αS can be tuned by two-dimensional (2D) crowding through simultaneous binding of a second protein on the bilayer surface. Unexpectedly, the experimental data show a clear deviation from a simple competitive inhibition model, but are consistent with a bimodal inhibition mechanism wherein membrane binding of a second protein (a membrane interacting chaperone, Hsp27, in this case) differentially inhibits two distinct modules of αS-membrane interaction. As a consequence, αS molecules are forced to access a hidden conformational state on the phospholipid bilayer in which only the higher-affinity module remains membrane-bound. Our results demonstrate that macromolecular crowding in two dimensions can play a significant role in shaping the conformational landscape of membrane-binding IDPs with multiple binding modes.