Machine-driven parameter screen of biochemical reactions.

The development of complex methods in molecular biology is a laborious, costly, iterative and often intuition-bound process where optima are sought in a multidimensional parameter space through step-by-step optimizations. The difficulty of miniaturizing reactions under the microliter volumes usually handled in multiwell plates by robots, plus the cost of the experiments, limit the number of parameters and the dynamic ranges that can be explored. Nevertheless, because of non-linearities of the response of biochemical systems to their reagent concentrations, broad dynamic ranges are necessary. Here we use a high-performance nanoliter handling platform and computer generation of liquid transfer programs to explore in quadruplicates 648 combinations of 4 parameters of a biochemical reaction, the reverse-transcription, which lead us to uncover non-linear responses, parameter interactions and novel mechanistic insights. With the increased availability of computer-driven laboratory platforms for biotechnology, our results demonstrate the feasibility and advantage of methods development based on reproducible, computer-aided exhaustive characterization of biochemical systems.

BAFfling pathologies: Alterations of BAF complexes in cancer.

To activate or repress specific genes, chromatin is constantly modified by chromatin-remodeling complexes. Among these complexes, the SWItch/Sucrose Non-Fermenting (SWI/SNF) complex, also referred to as BRG1-Associated Factor (BAF) complex, moves the nucleosome along chromatin using energy provided by ATP hydrolysis. In mammalian organisms, the SWI/SNF complex is composed of 10-15 subunits, depending on cell type, and a defect in one of these subunits can have dramatic consequences. In this review we will focus on the alterations identified in the SWI/SNF (BAF) complex subunits that lead to cancerous pathologies. While SMARCB1 was the first mutated subunit to be reported in a majority of malignant rhabdoid tumors, the advent of next-generation sequencing allowed the discovery of mutations in various SWI/SNF subunits within a broad spectrum of cancers. In most cases, the mutation leads to a loss of expression or to a truncated subunit unable to perform its function. Even though it is now commonly acknowledged that approximately 20% of all cancers present a mutation in a SWI/SNF subunit, some cancers are associated to a specific alteration of a SWI/SNF subunit, which acts either as tumor suppressor genes or as oncogenes, and therefore constitute diagnostic or prognostic biomarkers. Consistently, therapeutic strategies targeting SWI/SNF subunits or the genes affected downstream have been revealed to treat cancers.

NanoCAGE: A Method for the Analysis of Coding and Noncoding 5'-Capped Transcriptomes.

Transcripts in all eukaryotes are characterized by the 5'-end specific cap structure in mRNAs. Cap Analysis Gene Expression or CAGE makes use of these caps to specifically obtain cDNA fragments from the 5'-end of RNA and sequences those at high throughput for transcript identification and genome-wide mapping of transcription start sites for coding and noncoding genes. Here, we provide an improved version of our nanoCAGE protocol that has been developed for preparing CAGE libraries from as little as 50 ng of total RNA within three standard working days. Key steps in library preparation have been improved over our previously published protocol to obtain libraries having a good 5'-end selection and a more equal size distribution for higher sequencing efficiency on Illumina MiSeq and HiSeq sequencers. We recommend nanoCAGE as the method of choice for transcriptome profiling projects even from limited amounts of RNA, and as the best approach for genome-wide mapping of transcription start sites within promoter regions.

Single-Cell-Based Analysis Highlights a Surge in Cell-to-Cell Molecular Variability Preceding Irreversible Commitment in a Differentiation Process.

In some recent studies, a view emerged that stochastic dynamics governing the switching of cells from one differentiation state to another could be characterized by a peak in gene expression variability at the point of fate commitment. We have tested this hypothesis at the single-cell level by analyzing primary chicken erythroid progenitors through their differentiation process and measuring the expression of selected genes at six sequential time-points after induction of differentiation. In contrast to population-based expression data, single-cell gene expression data revealed a high cell-to-cell variability, which was masked by averaging. We were able to show that the correlation network was a very dynamical entity and that a subgroup of genes tend to follow the predictions from the dynamical network biomarker (DNB) theory. In addition, we also identified a small group of functionally related genes encoding proteins involved in sterol synthesis that could act as the initial drivers of the differentiation. In order to assess quantitatively the cell-to-cell variability in gene expression and its evolution in time, we used Shannon entropy as a measure of the heterogeneity. Entropy values showed a significant increase in the first 8 h of the differentiation process, reaching a peak between 8 and 24 h, before decreasing to significantly lower values. Moreover, we observed that the previous point of maximum entropy precedes two paramount key points: an irreversible commitment to differentiation between 24 and 48 h followed by a significant increase in cell size variability at 48 h. In conclusion, when analyzed at the single cell level, the differentiation process looks very different from its classical population average view. New observables (like entropy) can be computed, the behavior of which is fully compatible with the idea that differentiation is not a "simple" program that all cells execute identically but results from the dynamical behavior of the underlying molecular network.

Targeted reduction of highly abundant transcripts using pseudo-random primers.

Transcriptome studies based on quantitative sequencing can estimate levels of gene expression by measuring target RNA abundance in sequencing libraries. Sequencing costs are proportional to the total number of sequenced reads, and in order to cover rare RNAs, considerable quantities of abundant and identical reads are needed. This major limitation can be addressed by depleting a proportion of the most abundant sequences from the library. However, such depletion strategies involve either extra handling of the input RNA sample or use of a large number of reverse transcription primers, termed not-so-random (NSR) primers, which are costly to synthesize. Taking advantage of the high tolerance of reverse transcriptase to mis-prime, we found that it is possible to use as few as 40 pseudo-random (PS) reverse transcription primers to decrease the rate of undesirable abundant sequences within a library without affecting the overall transcriptome diversity. PS primers are simple to design and can be used to deplete several undesirable RNAs simultaneously, thus creating a flexible tool for enriching transcriptome libraries for rare transcript sequences.

Temperature-induced variation in gene expression burst size in metazoan cells.

BACKGROUND: Gene expression is an inherently stochastic process, owing to its dynamic molecular nature. Protein amount distributions, which can be acquired by cytometry using a reporter gene, can inform about the mechanisms of the underlying microscopic molecular system. RESULTS: By using different clones of chicken erythroid progenitor cells harboring different integration sites of a CMV-driven mCherry protein, we investigated the dynamical behavior of such distributions. We show that, on short term, clone distributions can be quickly regenerated from small population samples with a high accuracy. On longer term, on the contrary, we show variations manifested by correlated fluctuation in the Mean Fluorescence Intensity. In search for a possible cause of this correlation, we demonstrate that in response to small temperature variations cells are able to adjust their gene expression rate: a modest (2 °C) increase in external temperature induces a significant down regulation of mean expression values, with a reverse effect observed when the temperature is decreased. Using a two-state model of gene expression we further demonstrate that temperature acts by modifying the size of transcription bursts, while the burst frequency of the investigated promoter is less systematically affected. CONCLUSIONS: For the first time, we report that transcription burst size is a key parameter for gene expression that metazoan cells from homeotherm animals can modify in response to an external thermal stimulus.

Stochastic fluctuations and distributed control of gene expression impact cellular memory.

Despite the stochastic noise that characterizes all cellular processes the cells are able to maintain and transmit to their daughter cells the stable level of gene expression. In order to better understand this phenomenon, we investigated the temporal dynamics of gene expression variation using a double reporter gene model. We compared cell clones with transgenes coding for highly stable mRNA and fluorescent proteins with clones expressing destabilized mRNA-s and proteins. Both types of clones displayed strong heterogeneity of reporter gene expression levels. However, cells expressing stable gene products produced daughter cells with similar level of reporter proteins, while in cell clones with short mRNA and protein half-lives the epigenetic memory of the gene expression level was completely suppressed. Computer simulations also confirmed the role of mRNA and protein stability in the conservation of constant gene expression levels over several cell generations. These data indicate that the conservation of a stable phenotype in a cellular lineage may largely depend on the slow turnover of mRNA-s and proteins.

Understanding polyspecificity within the substrate-binding cavity of the human multidrug resistance P-glycoprotein.

Human P-glycoprotein (P-gp) controls drugs bioavailability by pumping structurally unrelated drugs out of cells. The X-ray structure of the mouse P-gp ortholog has been solved, with two SSS enantiomers or one RRR enantiomer of the selenohexapeptide inhibitor QZ59, found within the putative drug-binding pocket (Aller SG, Yu J, Ward A, Weng Y, Chittaboina S, Zhuo R, Harrell PM, Trinh YT, Zhang Q, Urbatsch IL et al. (2009). Science 323, 1718-1722). This offered the first opportunity to localize the well-known H and R drug-binding sites with respect to the QZ59 inhibition mechanisms of Hoechst 33342 and daunorubicin transports, characterized here in cellulo. We found that QZ59-SSS competes efficiently with both substrates, with K(I,app) values of 0.15 and 0.3 µM, which are 13 and 2 times lower, respectively, than the corresponding K(m,app) values. In contrast, QZ59-RRR non-competitively inhibited daunorubicin transport with moderate efficacy (K(I,app) = 1.9 µM); it also displayed a mixed-type inhibition of the Hoechst 33342 transport, resulting from a main non-competitive tendency (K(i2,app) = 1.6 µM) and a limited competitive tendency (K(i1,app) = 5 µM). These results suggest a positional overlap of QZ59 and drugs binding sites: full for the SSS enantiomer and partial for the RRR enantiomer. Crystal structure analysis suggests that the H site overlaps both QZ59-SSS locations while the R site overlaps the most embedded location.

Methoxy stilbenes as potent, specific, untransported, and noncytotoxic inhibitors of breast cancer resistance protein.

The ABCG2 multidrug transporter is known to confer cancer cell multidrug resistance by causing the efflux of anticancer drugs; therefore, selective inhibitors have the potential to improve chemotherapeutic treatments. Here, various methoxy derivatives of resveratrol are shown to be potent inhibitors of mitoxantrone efflux by ABCG2: among a series of 11 derivatives, compound 9 (3,5,3',4'-tetramethoxy trans-stilbene) had an IC(50) of 0.16 µM and showed a maximal inhibition of 75%, as measured by flow cytometry. It was not transported, as shown by HPLC fractionation and mass spectrometry titration and the lack of any cross-resistance in cell survival experiments. Compound 9 had a very low intrinsic cytotoxicity and was able to chemosensitize the growth of resistant ABCG2-transfected HEK293 cells at submicromolar concentrations. Drug-efflux inhibition was specific for ABCG2 since very low effects were observed with ABCB1 and ABCC1. The action mechanism of compound 9 was different from that of GF120918, which produced a complete and partly competitive but not ABCG2-specific inhibition, since ABCB1 was even more strongly inhibited. The two inhibitors also displayed different effects on the ABCG2 vanadate-sensitive ATPase activity, suggesting that they either bound to distinct sites or induced different conformational changes. Mitoxantrone efflux was fully inhibited by combining low concentrations of compound 9 with either GF120918 or a transport substrate such as prazosin or nilotinib. We conclude that methoxy derivatives of stilbene are good candidates for investigating future in vivo modulation of ABCG2 drug-efflux activity.

The multidrug resistance half-transporter ABCG2 is purified as a tetramer upon selective extraction from membranes.

ABCG2 is a human membrane ATP-binding cassette half-transporter that hydrolyzes ATP to efflux a large number of chemotherapeutic agents. Several oligomeric states of ABCG2 from homodimers to dodecamers have been reported depending on the overexpression systems and/or the protocols used for purification. Here, we compared the oligomeric state of His(6)-ABCG2 expressed in Sf9 insect cells and in human Flp-In-293/ABCG2 cells after solubilization in mild detergents. His(6)-ABCG2 was purified through a new approach involving its specific recognition onto a functionalized lipid layer containing a Ni-NTA lipid. This approach allowed the purification of His-ABCG2 in presence of all solubilized membrane components that might be involved in the stabilisation of native oligomers and without requiring any additional washing or concentration passages. ABCG2 purified onto the NiNTA lipid surfaces were directly analyzed by electron microscopy and by biochemical assays. Altogether, our data are consistent with a tetrameric organization of ABCG2 when expressed in either heterologous Sf9 insect cells or in human homologous cells.

Potent and fully noncompetitive peptidomimetic inhibitor of multidrug resistance P-glycoprotein.

N(α)-Boc-l-Asp(OBn)-l-Lys(Z)-OtBu (reversin 121, 1), an inhibitor of the P-gp ABC transporter, was used to conceive compounds inhibiting the drug efflux occurring through the Hoechst 33342 and daunorubicin transport sites of P-gp, respectively H and R sites. Replacement of the aspartyl residue by trans-4-hydroxy-l-proline (4(R)Hyp) gave compounds 11 and 15 characterized by half-maximal inhibitory concentrations (IC(50)) of 0.6 and 0.2 µM, which are 2- and 7-fold lower than that of the parent molecule. The difference in IC(50) between 11 and 15 rests on the carbonyl group of the peptidyl bond, reduced in 15. Those compounds are rather specific of P-gp, having no or limited activity on MRP1 and BCRP. 15 displayed no marked cytotoxicity up to 10-fold its IC(50). Importantly, 15 equally inhibited the Hoechst 33342 and daunorubicin effluxes through a typical noncompetitive inhibition mechanism, suggesting its binding to a site different from the H and R drug-transport sites.

ABCG2 transports and transfers heme to albumin through its large extracellular loop.

ABCG2 is an ATP-binding cassette (ABC) transporter preferentially expressed by immature human hematopoietic progenitors. Due to its role in drug resistance, its expression has been correlated with a protection role against protoporhyrin IX (PPIX) accumulation in stem cells under hypoxic conditions. We show here that zinc mesoporphyrin, a validated fluorescent heme analog, is transported by ABCG2. We also show that the ABCG2 large extracellular loop ECL3 constitutes a porphyrin-binding domain, which strongly interacts with heme, hemin, PPIX, ZnPPIX, CoPPIX, and much less efficiently with pheophorbide a, but not with vitamin B12. K(d) values are in the range 0.5-3.5 µm, with heme displaying the highest affinity. Nonporphyrin substrates of ABCG2, such as mitoxantrone, doxo/daunorubicin, and riboflavin, do not bind to ECL3. Single-point mutations H583A and C603A inside ECL3 prevent the binding of hemin but hardly affect that of iron-free PPIX. The extracellular location of ECL3 downstream from the transport sites suggests that, after membrane translocation, hemin is transferred to ECL3, which is strategically positioned to release the bound porphyrin to extracellular partners. We show here that human serum albumin could be one of these possible partners as it removes hemin bound to ECL3 and interacts with ABCG2, with a K(d) of about 3 µm.