Deconstructing Best-in-Class Neoglycoclusters as a Tool for Dissecting Key Multivalent Processes in Glycosidase Inhibition.

Multivalency represents an appealing option to modulate selectivity in enzyme inhibition and transform moderate glycosidase inhibitors into highly potent ones. The rational design of multivalent inhibitors is however challenging because global affinity enhancement relies on several interconnected local mechanistic events, whose relative impact is unknown. So far, the largest multivalent effects ever reported for a non-polymeric glycosidase inhibitor have been obtained with cyclopeptoid-based inhibitors of Jack bean α-mannosidase (JBα-man). Here, we report a structure-activity relationship (SAR) study based on the top-down deconstruction of best-in-class multivalent inhibitors. This approach provides a valuable tool to understand the complex interdependent mechanisms underpinning the inhibitory multivalent effect. Combining SAR experiments, binding stoichiometry assessments, thermodynamic modelling and atomistic simulations allowed us to establish the significant contribution of statistical rebinding mechanisms and the importance of several key parameters, including inhitope accessibility, topological restrictions, and electrostatic interactions. Our findings indicate that strong chelate-binding, resulting from the formation of a cross-linked complex between a multivalent inhibitor and two dimeric JBα-man molecules, is not a sufficient condition to reach high levels of affinity enhancements. The deconstruction approach thus offers unique opportunities to better understand multivalent binding and provides important guidelines for the design of potent and selective multiheaded inhibitors.

Titration of Apparent In-Cellula Affinities of Protein-Protein Interactions.

A genetic assay permits simultaneous quantification of two interacting proteins and their bound fraction at the single-cell level using flow cytometry. Apparent in-cellula affinities of protein-protein interactions can be extracted from the acquired data through a titration-like analysis. The applicability of this approach is demonstrated on a diverse set of interactions with proteins from different families and organisms and with in-vitro dissociation constants ranging from picomolar to micromolar.

A Quantitative Tri-fluorescent Yeast Two-hybrid System: From Flow Cytometry to In cellula Affinities.

We present a technological advancement for the estimation of the affinities of Protein-Protein Interactions (PPIs) in living cells. A novel set of vectors is introduced that enables a quantitative yeast two-hybrid system based on fluorescent fusion proteins. The vectors allow simultaneous quantification of the reaction partners (Bait and Prey) and the reporter at the single-cell level by flow cytometry. We validate the applicability of this system on a small but diverse set of PPIs (eleven protein families from six organisms) with different affinities; the dissociation constants range from 117 pm to 17 µm After only two hours of reaction, expression of the reporter can be detected even for the weakest PPI. Through a simple gating analysis, it is possible to select only cells with identical expression levels of the reaction partners. As a result of this standardization of expression levels, the mean reporter levels directly reflect the affinities of the studied PPIs. With a set of PPIs with known affinities, it is straightforward to construct an affinity ladder that permits rapid classification of PPIs with thus far unknown affinities. Conventional software can be used for this analysis. To permit automated analysis, we provide a graphical user interface for the Python-based FlowCytometryTools package.

Structure-activity relationship and molecular modelling studies of quinazolinedione derivatives MMV665916 as potential antimalarial agent.

A series of new quinazolinedione derivatives have been readily synthesized and evaluated for their in vitro antiplasmodial growth inhibition activity. Most of the compounds inhibited P. falciparum FcB1 strain in the low to medium micromolar concentration. The 2-ethoxy 8ag', 2-trifluoromethoxy 8ai' and 4-fluoro-2-methoxy 8ak' showed the best inhibitory activity with EC50 values around 5 µM and were non-toxic to the primary human fibroblast cell line AB943. However, these compounds were less potent than the original hit MMV665916, which showed remarkable growth inhibition with EC50 value of 0.4 µM and presented the highest selectivity index (SI > 250). In addition, a novel approach for determining the docking poses of these quinazolinedione derivatives with their potential protein target, the P. falciparum farnesyltransferase PfFT, was investigated.

Physical and functional interaction between SET1/COMPASS complex component CFP-1 and a Sin3S HDAC complex in C. elegans.

The CFP1 CXXC zinc finger protein targets the SET1/COMPASS complex to non-methylated CpG rich promoters to implement tri-methylation of histone H3 Lys4 (H3K4me3). Although H3K4me3 is widely associated with gene expression, the effects of CFP1 loss vary, suggesting additional chromatin factors contribute to context dependent effects. Using a proteomics approach, we identified CFP1 associated proteins and an unexpected direct link between Caenorhabditis elegans CFP-1 and an Rpd3/Sin3 small (SIN3S) histone deacetylase complex. Supporting a functional connection, we find that mutants of COMPASS and SIN3 complex components genetically interact and have similar phenotypic defects including misregulation of common genes. CFP-1 directly binds SIN-3 through a region including the conserved PAH1 domain and recruits SIN-3 and the HDA-1/HDAC subunit to H3K4me3 enriched promoters. Our results reveal a novel role for CFP-1 in mediating interaction between SET1/COMPASS and a Sin3S HDAC complex at promoters.

Assigning function to natural allelic variation via dynamic modeling of gene network induction.

More and more natural DNA variants are being linked to physiological traits. Yet, understanding what differences they make on molecular regulations remains challenging. Important properties of gene regulatory networks can be captured by computational models. If model parameters can be "personalized" according to the genotype, their variation may then reveal how DNA variants operate in the network. Here, we combined experiments and computations to visualize natural alleles of the yeast GAL3 gene in a space of model parameters describing the galactose response network. Alleles altering the activation of Gal3p by galactose were discriminated from those affecting its activity (production/degradation or efficiency of the activated protein). The approach allowed us to correctly predict that a non-synonymous SNP would change the binding affinity of Gal3p with the Gal80p transcriptional repressor. Our results illustrate how personalizing gene regulatory models can be used for the mechanistic interpretation of genetic variants.

Calcineurin A versus NS5A-TP2/HD domain containing 2: a case study of site-directed low-frequency random mutagenesis for dissecting target specificity of peptide aptamers.

We previously identified a peptide aptamer (named R5G42) via functional selection for its capacity to slow cell proliferation. A yeast two-hybrid screen of human cDNA libraries, using R5G42 as "bait," allowed the identification of two binding proteins with very different functions: calcineurin A (CnA) (PP2B/PPP3CA), a protein phosphatase well characterized for its role in the immune response, and NS5A-TP2/HD domain containing 2, a much less studied protein induced subsequent to hepatitis C virus non-structural protein 5A expression in HepG2 hepatocellular carcinoma cells, with no known activity. Our objective in the present study was to dissect the dual target specificity of R5G42 in order to have tools with which to better characterize the actions of the peptide aptamers toward their individual targets. This was achieved through the selection of random mutants of the variable loop, derived from R5G42, evaluating their specificity toward CnA and NS5A-TP2 and analyzing their sequence. An interdisciplinary approach involving biomolecular computer simulations with integration of the sequence data and yeast two-hybrid binding phenotypes of these mutants yielded two structurally distinct conformers affording the potential molecular basis of the binding diversity of R5G42. Evaluation of the biological impact of CnA- versus NS5A-TP2-specific peptide aptamers indicated that although both contributed to the anti-proliferative effect of R5G42, CnA-binding was essential to stimulate the nuclear translocation of nuclear factor of activated T cells, indicative of the activation of endogenous CnA. By dissecting the target specificity of R5G42, we have generated novel tools with which to study each target individually. Apta-C8 is capable of directly activating CnA independent of binding to NS5A-TP2 and will be an important tool in studying the role of CnA activation in the regulation of different signaling pathways, whereas Apta-E1 will allow dissection of the function of NS5A-TP2, serving as an example of the usefulness of peptide aptamer technology for investigating signaling pathways.

Thermodynamics and kinetics of large-time-step molecular dynamics.

Molecular dynamics (MD) simulations provide essential information about the thermodynamics and kinetics of proteins. Technological advances in both hardware and algorithms have seen this method accessing timescales that used to be unreachable only few years ago. The quest to simulate slow, biologically relevant macromolecular conformational changes, is still open. Here, we present an approximate approach to increase the speed of MD simulations by a factor of ∼4.5. This is achieved by using a large integration time step of 7 fs, in combination with frozen covalent bonds and look-up tables for nonbonded interactions of the solvent. Extensive atomistic MD simulations for a flexible peptide in water show that the approach reproduces the peptide's equilibrium conformational changes, preserving the essential properties of both thermodynamics and kinetics. Comparison of this approximate method with state-of-the-art implicit solvation simulations indicates that the former provides a better description of the underlying free-energy surface. Finally, simulations of a 33-residue peptide show that these fast MD settings are readily applicable to investigate biologically relevant systems.

A single-chain and fast-responding light-inducible Cre recombinase as a novel optogenetic switch.

Optogenetics enables genome manipulations with high spatiotemporal resolution, opening exciting possibilities for fundamental and applied biological research. Here, we report the development of LiCre, a novel light-inducible Cre recombinase. LiCre is made of a single flavin-containing protein comprising the AsLOV2 photoreceptor domain of Avena sativa fused to a Cre variant carrying destabilizing mutations in its N-terminal and C-terminal domains. LiCre can be activated within minutes of illumination with blue light without the need of additional chemicals. When compared to existing photoactivatable Cre recombinases based on two split units, LiCre displayed faster and stronger activation by light as well as a lower residual activity in the dark. LiCre was efficient both in yeast, where it allowed us to control the production of ß-carotene with light, and human cells. Given its simplicity and performances, LiCre is particularly suited for fundamental and biomedical research, as well as for controlling industrial bioprocesses.

In a biologist's toolkit, the Cre protein holds a special place. Naturally found in certain viruses, this enzyme recognises and modifies specific genetic sequences, creating changes that switch on or off whatever gene is close by. Genetically engineering cells or organisms so that they carry Cre and its target sequences allows scientists to control the activation of a given gene, often in a single tissue or organ. However, this relies on the ability to activate the Cre protein 'on demand' once it is in the cells of interest. One way to do so is to split the enzyme into two pieces, which can then reassemble when exposed to blue light. Yet, this involves the challenging step of introducing both parts separately into a tissue. Instead, Duplus-Bottin et al. engineered LiCre, a new system where a large section of the Cre protein is fused to a light sensor used by oats to detect their environment. LiCre is off in the dark, but it starts to recognize and modify Cre target sequences when exposed to blue light. Duplus-Bottin et al. then assessed how LiCre compares to the two-part Cre system in baker's yeast and human kidney cells. This showed that the new protein is less 'incorrectly' active in the dark, and can switch on faster under blue light. The improved approach could give scientists a better tool to study the role of certain genes at precise locations and time points, but also help them to harness genetic sequences for industry or during gene therapy.

A differential fluctuation theorem.

We derive a nonequilibrium thermodynamics identity (the "differential fluctuation theorem") that connects forward and reverse joint probabilities of nonequilibrium work and of arbitrary generalized coordinates corresponding to states of interest. This identity allows us to estimate the free energy difference between domains of these states. Our results follow from a general symmetry relation between averages over nonequilibrium forward and backward path functions derived by Crooks [Crooks, G. E. Phys. Rev. E 2000, 61, 2361-2366]. We show how several existing nonequilibrium thermodynamic identities can be obtained directly from the differential fluctuation theorem. We devise an approach for measuring conformational free energy differences, and we demonstrate its applicability to the analysis of molecular dynamics simulations by estimating the free energy difference between two conformers of the alanine dipeptide model system. We anticipate that these developments can be applied to the analysis of laboratory experiments.

Eicosapentaenoic acid modulates the synergistic action of CREB1 and ID/E2A family members in the rat pup brain and mouse embryonic stem cells.

The aim of this study was to investigate the molecular mechanism by which eicosapentaenoic acid (EPA) may exert neuroprotective effects through an "EPA-cyclic AMP response element-binding protein (CREB)" signaling pathway. The current study reveals that EPA modulates the exquisite interplay of interaction of CREB1 with the inhibitor of DNA binding (ID) and E2A family members, thereby delivering mechanistic insights into specific neural differentiation program. In this scenario, our work provides evidence for the capability of CREB1 to sequester ID:E2A family members in brain tissues and neural differentiating mouse embryonic stem cells (mESCs) through formation of a [CREB1]2:ID2:E47 tetrameric complex.In essence, the molecular function of CREB1 is to dynamically regulate the location-specific assembly or disassembly of basic-helix-loop-helix (bHLH):HLH protein complexes to mediate the activation of neural/glial target genes. Together, these findings support the one-to-many binding mechanism of CREB1 and indicate that EPA treatment potentiates the integration of CREB dependent signaling with HLH/bHLH transcriptional network, adding specificity to the CREB1-mediated gene regulation during neural/glial differentiation. Our current research on the EPA-CREB axis could reveal new molecular targets for treating neurogenerative disease.

Bis(mesitoyl)phosphinic acid: photo-triggered release of metaphosphorous acid in solution.

Bis(mesitoyl)phosphinic acid and its sodium salt display a unique photo-induced reactivity: both derivatives stepwise release two mesitoyl radicals and, remarkably, metaphosphorous acid (previously postulated as transient species in the gas phase), providing a new phosphorus-based reagent.

Electron Transfer through DNA in the Course of Radical-Induced Strand Cleavage.

No benefit from base stacking is observed for rates of electron transfer in DNA. This conclusion was drawn from experiments with a new DNA assay in which a radical cationic site, generated by strand cleavage, can be reduced by the guanine bases in the same DNA (the electron transfer is indicated by arrows in the diagram). The distance dependence of this electron transfer step is determined by the chemical yield of the reduction product.

Automated high-throughput quantification of mitotic spindle positioning from DIC movies of Caenorhabditis embryos.

The mitotic spindle is a microtubule-based structure that elongates to accurately segregate chromosomes during anaphase. Its position within the cell also dictates the future cell cleavage plan, thereby determining daughter cell orientation within a tissue or cell fate adoption for polarized cells. Therefore, the mitotic spindle ensures at the same time proper cell division and developmental precision. Consequently, spindle dynamics is the matter of intensive research. Among the different cellular models that have been explored, the one-cell stage C. elegans embryo has been an essential and powerful system to dissect the molecular and biophysical basis of spindle elongation and positioning. Indeed, in this large and transparent cell, spindle poles (or centrosomes) can be easily detected from simple DIC microscopy by human eyes. To perform quantitative and high-throughput analysis of spindle motion, we developed a computer program ACT for Automated-Centrosome-Tracking from DIC movies of C. elegans embryos. We therefore offer an alternative to the image acquisition and processing of transgenic lines expressing fluorescent spindle markers. Consequently, experiments on large sets of cells can be performed with a simple setup using inexpensive microscopes. Moreover, analysis of any mutant or wild-type backgrounds is accessible because laborious rounds of crosses with transgenic lines become unnecessary. Last, our program allows spindle detection in other nematode species, offering the same quality of DIC images but for which techniques of transgenesis are not accessible. Thus, our program also opens the way towards a quantitative evolutionary approach of spindle dynamics. Overall, our computer program is a unique macro for the image- and movie-processing platform ImageJ. It is user-friendly and freely available under an open-source licence. ACT allows batch-wise analysis of large sets of mitosis events. Within 2 minutes, a single movie is processed and the accuracy of the automated tracking matches the precision of the human eye.

Synthesis, spectroscopic and DNA alkylating properties of malondialdehyde (MDA) bis-imine fluorescent adducts.

The synthesis of a series of malondialdehyde (MDA) fluorescent adducts that mimic the well-known pentamethine cyanine dyes is reported. This new subclass of bis-imino dyes shares some common spectroscopic properties with their polymethine analogues although absorbing at significantly shorter wavelengths. A small library of trimethine and pentamethine cyanine dye bis-imino analogues have been synthesised and characterised that cover a spectral range from blue to orange. Of particular interest is their capacity to act as mono- and bis-alkylating agents of nucleosides in general and of cytidine (and 2'-deoxycytidine) in particular.

Reversible synthesis and characterization of dynamic imino analogues of trimethine and pentamethine cyanine dyes.

A new family of unsymmetrical imine-based trimethine and pentamethine cyanine dye analogues is reported that can form under reversible and thermodynamically controlled conditions from non- or weakly emissive amine and aldehyde building blocks. These dynamic fluorophores show spectroscopic properties comparable to those of their parent cyanine dyes and are responsive to external effectors.

The HSP90 binding mode of a radicicol-like E-oxime determined by docking, binding free energy estimations, and NMR 15N chemical shifts.

We determine the binding mode of a macrocyclic radicicol-like oxime to yeast HSP90 by combining computer simulations and experimental measurements. We sample the macrocyclic scaffold of the unbound ligand by parallel tempering simulations and dock the most populated conformations to yeast HSP90. Docking poses are then evaluated by the use of binding free energy estimations with the linear interaction energy method. Comparison of QM/MM-calculated NMR chemical shifts with experimental shift data for a selective subset of backbone (15)N provides an additional evaluation criteria. As a final test we check the binding modes against available structure-activity-relationships. We find that the most likely binding mode of the oxime to yeast HSP90 is very similar to the known structure of the radicicol-HSP90 complex.

Optimal estimates of free energies from multistate nonequilibrium work data.

We derive the optimal estimates of the free energies of an arbitrary number of thermodynamic states from nonequilibrium work measurements; the work data are collected from forward and reverse switching processes and obey a fluctuation theorem. The maximum likelihood formulation properly reweights all pathways contributing to a free energy difference and is directly applicable to simulations and experiments. We demonstrate dramatic gains in efficiency by combining the analysis with parallel tempering simulations for alchemical mutations of model amino acids.

Role of conformational heterogeneity in domain swapping and adapter function of the Cks proteins.

Cks proteins are adapter molecules that coordinate the assembly of multiprotein complexes. They share the ability to domain swap by exchanging a beta-strand, beta4. Here we use NMR spectroscopy and molecular dynamics simulations to investigate the dynamic properties of human Cks1 and its response on assembly with components of the SCF(Skp2) ubiquitin ligation machinery. In the NMR experiment with the free form of Cks1, a subset of residues displayed elevated R2 values and the cross-peaks of neighboring residues were missing from the spectrum, indicating a substantial conformational exchange contribution on the microsecond to millisecond time scale. Strikingly the region of greatest conformational variability was the beta4-strand that domain swaps to form the dimer. Binding of the ligand common to all Cks proteins, Cdk2, suppressed the conformational heterogeneity. This response was specific to Cdk2 binding; in contrast, binding of Skp2, a ligand unique to human Cks1, did not alter the dynamic behavior. Short time (<5 ns) molecular dynamics simulations indicate that residues of Cks1 that form the binding site for phosphorylated ligands are considerably more flexible in the free form of Cks1 than they are in the Cdk2-Cks1 complex. A cooperative interaction between Cdk2 and Cks1 is suggested, which reduces the configurational entropy of Cks1 and therefore facilitates phosphoprotein binding. Indications of an unusual dynamic behavior of strand beta4 in the free form of Cks1 were obtained from longer time scale (50 ns) dynamics simulations. A spontaneous reversible unzipping of hydrogen bonds between beta4 and beta2 was observed, suggesting an early intermediate structure for unfolding and/or domain swapping. We propose that the dynamic properties of the beta-sheet and its modification upon ligand binding underlie the domain swapping ability and the adapter function of Cks proteins.