Inference Models in DNA Computing

La computación molecular es una disciplina que se ocupa del diseño e implementación de dispositivos para el procesamiento de información sobre un sustrato biológico, como el ácido desoxirribonucleico (ADN), el ácido ribonucleico (ARN) o las proteínas. Desde que Watson y Crick descubrieron en los años cincuenta la estructura molecular del ADN en forma de doble hélice, se desencadenaron otros descubrimientos, como las enzimas de restricción o la reacción en cadena de la polimerasa (PCR), contribuyendo de manera determinante a la irrupción de la tecnología del ADN recombinante. Gracias a esta tecnología y al descenso vertiginoso de los precios de secuenciación y síntesis del ADN, la computación biomolecular pudo abandonar su concepción puramente teórica. El trabajo presentado por Adleman (1994) logró resolver un problema de computación NP-completo (El Problema del Camino de Hamilton dirigido) utilizando únicamente moléculas de ADN. La gran capacidad de procesamiento en paralelo ofrecida por las técnicas del ADN recombinante permitió a Adleman ser capaz de resolver dicho problema en tiempo polinómico, aunque a costa de un consumo exponencial de moléculas de ADN. Utilizando algoritmos de fuerza bruta similares al utilizado por Adleman se logró resolver otros problemas NP-completos, como por ejemplo el de Satisfacibilidad de Fórmulas Lógicas / SAT (Lipton, 1995). Pronto se comprendió que la computación biomolecular no podía competir en velocidad ni precisión con los ordenadores de silicio, por lo que su enfoque y objetivos se centraron en la resolución de problemas con aplicación biomédica (Simmel, 2007), dejando de lado la resolución de problemas clásicos de computación. Desde entonces se han propuesto diversos modelos de dispositivos biomoleculares que, de forma autónoma (sin necesidad de un bio-ingeniero realizando operaciones de laboratorio), son capaces de procesar como entrada un sustrato biológico y proporcionar una salida también en formato biológico: procesadores que aprovechan la extensión de la polimerasa (Hagiya et al., 1997), autómatas que funcionan con enzimas de restricción (Benenson et al., 2001) o con deoxiribozimas (Stojanovic et al., 2002), o circuitos de hibridación competitiva (Yurke et al., 2000). Esta tesis presenta un conjunto de modelos de dispositivos de ácidos nucleicos capaces de implementar diversas operaciones de computación lógica aprovechando técnicas de computación biomolecular (hibridación competitiva del ADN y reacciones enzimáticas) con aplicaciones en diagnóstico genético. El primer conjunto de modelos, presentados en el Capítulo 5 y publicados en Sainz de Murieta and Rodríguez-Patón (2012b), Rodríguez-Patón et al. (2010a) y Sainz de Murieta and Rodríguez-Patón (2010), define un tipo de biosensor que usa hebras simples de ADN para codificar reglas sencillas, como por ejemplo "SI hebra-ADN-1 Y hebra-ADN-2 presentes, ENTONCES enfermedad-B". Estas reglas interactúan con señales de entrada (ADN o ARN de cualquier tipo) para producir una señal de salida (también en forma de ácido nucleico). Dicha señal de salida representa un diagnóstico, que puede medirse mediante partículas fluorescentes técnicas FRET) o incluso ser un tratamiento administrado en respuesta a un conjunto de síntomas. El modelo presentado en el Capítulo 5, publicado en Rodríguez-Patón et al. (2011), es capaz de ejecutar cadenas de resolución sobre fórmulas lógicas en forma normal conjuntiva. Cada cláusula de una fórmula se codifica en una molécula de ADN. Cada proposición p se codifica asignándole una hebra simple de ADN, y la correspondiente hebra complementaria a la proposición ¬p. Las cláusulas se codifican incluyendo distintas proposiciones en la misma hebra de ADN. El modelo permite ejecutar programas lógicos de cláusulas Horn aplicando múltiples iteraciones de resolución en cascada, con el fin de implementar la función de un nanodispositivo autónomo programable. Esta técnica también puede emplearse para resolver SAP sin ayuda externa. El modelo presentado en el Capítulo 6 se ha publicado en publicado en Sainz de Murieta and Rodríguez-Patón (2012c), y el modelo presentado en el Capítulo 7 se ha publicado en (Sainz de Murieta and Rodríguez-Patón, 2013c). Aunque explotan métodos de computación biomolecular diferentes (hibridación competitiva de ADN en el Capítulo 6 frente a reacciones enzimáticas en el 7), ambos modelos son capaces de realizar inferencia Bayesiana. Funcionan tomando hebras simples de ADN como entrada, representando la presencia o la ausencia de un indicador molecular concreto (una evidencia). La probabilidad a priori de una enfermedad, así como la probabilidad condicionada de una señal (o síntoma) dada la enfermedad representan la base de conocimiento, y se codifican combinando distintas moléculas de ADN y sus concentraciones relativas. Cuando las moléculas de entrada interaccionan con las de la base de conocimiento, se liberan dos clases de hebras de ADN, cuya proporción relativa representa la aplicación del teorema de Bayes: la probabilidad condicionada de la enfermedad dada la señal (o síntoma). Todos estos dispositivos pueden verse como elementos básicos que, combinados modularmente, permiten la implementación de sistemas in vitro a partir de sensores de ADN, capaces de percibir y procesar señales biológicas. Este tipo de autómatas tienen en la actualidad una gran potencial, además de una gran repercusión científica. Un perfecto ejemplo fue la publicación de (Xie et al., 2011) en Science, presentando un autómata biomolecular de diagnóstico capaz de activar selectivamente el proceso de apoptosis en células cancerígenas sin afectar a células sanas.

Single-pair fluorescence resonance energy transfer on freely diffusing molecules: Observation of Förster distance dependence and subpopulations

Photon bursts from single diffusing donor-acceptor labeled macromolecules were used to measure intramolecular distances and identify subpopulations of freely diffusing macromolecules in a heterogeneous ensemble. By using DNA as a rigid spacer, a series of constructs with varying intramolecular donor-acceptor spacings were used to measure the mean and distribution width of fluorescence resonance energy transfer (FRET) efficiencies as a function of distance. The mean single-pair FRET efficiencies qualitatively follow the distance dependence predicted by Förster theory. Possible contributions to the widths of the FRET efficiency distributions are discussed, and potential applications in the study of biopolymer conformational dynamics are suggested. The ability to measure intramolecular (and intermolecular) distances for single molecules implies the ability to distinguish and monitor subpopulations of molecules in a mixture with different distances or conformational states. This is demonstrated by monitoring substrate and product subpopulations before and after a restriction endonuclease cleavage reaction. Distance measurements at single-molecule resolution also should facilitate the study of complex reactions such as biopolymer folding. To this end, the denaturation of a DNA hairpin was examined by using single-pair FRET.

Preassembly of interleukin 2 (IL-2) receptor subunits on resting Kit 225 K6 T cells and their modulation by IL-2, IL-7, and IL-15: A fluorescence resonance energy transfer study

Assembly and mutual proximities of α, β, and γc subunits of the interleukin 2 receptors (IL-2R) in plasma membranes of Kit 225 K6 T lymphoma cells were investigated by fluorescence resonance energy transfer (FRET) using fluorescein isothiocyanate- and Cy3-conjugated monoclonal antibodies (mAbs) that were directed against the IL-2Rα, IL-2Rβ, and γc subunits of IL-2R. The cell-surface distribution of subunits was analyzed at the nanometer scale (2–10 nm) by FRET on a cell-by-cell basis. The cells were probed in resting phase and after coculture with saturating concentrations of IL-2, IL-7, and IL-15. FRET data from donor- and acceptor-labeled IL-2Rβ-α, γ-α, and γ-β pairs demonstrated close proximity of all subunits to each other in the plasma membrane of resting T cells. These mutual proximities do not appear to represent mAb-induced microaggregation, because FRET measurements with Fab fragments of the mAbs gave similar results. The relative proximities were meaningfully modulated by binding of IL-2, IL-7, and IL-15. Based on FRET analysis the topology of the three subunits at the surface of resting cells can be best described by a “triangular model” in the absence of added interleukins. IL-2 strengthens the bridges between the subunits, making the triangle more compact. IL-7 and IL-15 act in the opposite direction by opening the triangle possibly because they associate their private specific α receptors with the β and/or γc subunits of the IL-2R complex. These data suggest that IL-2R subunits are already colocalized in resting T cells and do not require cytokine-induced redistribution. This colocalization is significantly modulated by binding of relevant interleukins in a cytokine-specific manner.

Single-molecule protein folding: Diffusion fluorescence resonance energy transfer studies of the denaturation of chymotrypsin inhibitor 2

We report single-molecule folding studies of a small, single-domain protein, chymotrypsin inhibitor 2 (CI2). CI2 is an excellent model system for protein folding studies and has been extensively studied, both experimentally (at the ensemble level) and theoretically. Conformationally assisted ligation methodology was used to synthesize the proteins and site-specifically label them with donor and acceptor dyes. Folded and denatured subpopulations were observed by fluorescence resonance energy transfer (FRET) measurements on freely diffusing single protein molecules. Properties of these subpopulations were directly monitored as a function of guanidinium chloride concentration. It is shown that new information about different aspects of the protein folding reaction can be extracted from such subpopulation properties. Shifts in the mean transfer efficiencies are discussed, FRET efficiency distributions are translated into potentials, and denaturation curves are directly plotted from the areas of the FRET peaks. Changes in stability caused by mutation also are measured by comparing pseudo wild-type CI2 with a destabilized mutant (K17G). Current limitations and future possibilities and prospects for single-pair FRET protein folding investigations are discussed.

Role of metallothionein in nitric oxide signaling as revealed by a green fluorescent fusion protein

Although the function of metallothionein (MT), a 6- to 7-kDa cysteine-rich metal binding protein, remains unclear, it has been suggested from in vitro studies that MT is an important component of intracellular redox signaling, including being a target for nitric oxide (NO). To directly study the interaction between MT and NO in live cells, we generated a fusion protein consisting of MT sandwiched between two mutant green fluorescent proteins (GFPs). In vitro studies with this chimera (FRET-MT) demonstrate that fluorescent resonance energy transfer (FRET) can be used to follow conformational changes indicative of metal release from MT. Imaging experiments with live endothelial cells show that agents that increase cytoplasmic Ca2+ act via endogenously generated NO to rapidly and persistently release metal from MT. A role for this interaction in intact tissue is supported by the finding that the myogenic reflex of mesenteric arteries is absent in MT knockout mice (MT−/−) unless endogenous NO synthesis is blocked. These results are the first application of intramolecular green fluorescent protein (GFP)-based FRET in a native protein and demonstrate the utility of FRET-MT as an intracellular surrogate indicator of NO production. In addition, an important role of metal thiolate clusters of MT in NO signaling in vascular tissue is revealed.

A quantitative method for evaluating the stabilities of nucleic acids

We report a general method for screening, in solution, the impact of deviations from canonical Watson-Crick composition on the thermodynamic stability of nucleic acid duplexes. We demonstrate how fluorescence resonance energy transfer (FRET) can be used to detect directly free energy differences between an initially formed “reference” duplex (usually a Watson-Crick duplex) and a related “test” duplex containing a lesion/alteration of interest (e.g., a mismatch, a modified, a deleted, or a bulged base, etc.). In one application, one titrates into a solution containing a fluorescently labeled, FRET-active, reference duplex, an unlabeled, single-stranded nucleic acid (test strand), which may or may not compete successfully to form a new duplex. When a new duplex forms by strand displacement, it will not exhibit FRET. The resultant titration curve (normalized fluorescence intensity vs. logarithm of test strand concentration) yields a value for the difference in stability (free energy) between the newly formed, test strand-containing duplex and the initial reference duplex. The use of competitive equilibria in this assay allows the measurement of equilibrium association constants that far exceed the magnitudes accessible by conventional titrimetric techniques. Additionally, because of the sensitivity of fluorescence, the method requires several orders of magnitude less material than most other solution methods. We discuss the advantages of this method for detecting and characterizing any modification that alters duplex stability, including, but not limited to, mutagenic lesions. We underscore the wide range of accessible free energy values that can be defined by this method, the applicability of the method in probing for a myriad of nucleic acid variations, such as single nucleotide polymorphisms, and the potential of the method for high throughput screening.

Spatiotemporal dynamics of guanosine 3′,5′-cyclic monophosphate revealed by a genetically encoded, fluorescent indicator

To investigate the dynamics of guanosine 3′,5′-cyclic monophosphate (cGMP) in single living cells, we constructed genetically encoded, fluorescent cGMP indicators by bracketing cGMP-dependent protein kinase (cGPK), minus residues 1–77, between cyan and yellow mutants of green fluorescent protein. cGMP decreased fluorescence resonance energy transfer (FRET) and increased the ratio of cyan to yellow emissions by up to 1.5-fold with apparent dissociation constants of ≈2 μM and >100:1 selectivity for cGMP over cAMP. To eliminate constitutive kinase activity, Thr516 of cGPK was mutated to Ala. Emission ratio imaging of the indicators transfected into rat fetal lung fibroblast (RFL)-6 showed cGMP transients resulting from activation of soluble and particulate guanylyl cyclase, respectively, by nitric oxide (NO) and C-type natriuretic peptide (CNP). Whereas all naive cells tested responded to CNP, only 68% responded to NO. Both sets of signals showed large and variable (0.5–4 min) latencies. The phosphodiesterase (PDE) inhibitor 3-isobutyl-1-methylxanthine (IBMX) did not elevate cGMP on its own but consistently amplified responses to NO or CNP, suggesting that basal activity of guanylate cyclase is very low and emphasizing the importance of PDEs in cGMP recycling. A fraction of RFL cells showed slowly propagating tides of cGMP spreading across the cell in response to delocalized application of NO. Biolistically transfected Purkinje neurons showed cGMP responses to parallel fiber activity and NO donors, confirming that single-cell increases in cGMP occur under conditions appropriate to cause synaptic plasticity.

Non-Arrhenius kinetics for the loop closure of a DNA hairpin

Intramolecular chain diffusion is an elementary process in the conformational fluctuations of the DNA hairpin-loop. We have studied the temperature and viscosity dependence of a model DNA hairpin-loop by FRET (fluorescence resonance energy transfer) fluctuation spectroscopy (FRETfs). Apparent thermodynamic parameters were obtained by analyzing the correlation amplitude through a two-state model and are consistent with steady-state fluorescence measurements. The kinetics of closing the loop show non-Arrhenius behavior, in agreement with theoretical prediction and other experimental measurements on peptide folding. The fluctuation rates show a fractional power dependence (β = 0.83) on the solution viscosity. A much slower intrachain diffusion coefficient in comparison to that of polypeptides was derived based on the first passage time theory of SSS [Szabo, A., Schulten, K. & Schulten, Z. (1980) J. Chem. Phys. 72, 4350–4357], suggesting that intrachain interactions, especially stacking interaction in the loop, might increase the roughness of the free energy surface of the DNA hairpin-loop.

The synaptic activity of HsDmc1, a human recombination protein specific to meiosis

Human Dmc1 protein, a meiosis-specific homolog of Escherichia coli RecA protein, has previously been shown to promote DNA homologous pairing and strand-exchange reactions that are qualitatively similar to those of RecA protein and Rad51. Human and yeast Rad51 proteins each form a nucleoprotein filament that is very similar to the filament formed by RecA protein. However, recent studies failed to find a similar filament made by Dmc1 but showed instead that this protein forms octameric rings and stacks of rings. These observations stimulated further efforts to elucidate the mechanism by which Dmc1 promotes the recognition of homology. Dmc1, purified to a state in which nuclease and helicase activities were undetectable, promoted homologous pairing and strand exchange as measured by fluorescence resonance energy transfer (FRET). Observations on the intermediates and products, which can be distinguished by FRET assays, provided direct evidence of a three-stranded synaptic intermediate. The effects of helix stability and mismatched base pairs on the recognition of homology revealed further that human Dmc1, like human Rad51, requires the preferential breathing of A⋅T base pairs for recognition of homology. We conclude that Dmc1, like human Rad51 and E. coli RecA protein, promotes homologous pairing and strand exchange by a “synaptic pathway” involving a three-stranded nucleoprotein intermediate, rather than by a “helicase pathway” involving the separation and reannealing of DNA strands.

Probing the interaction between two single molecules: fluorescence resonance energy transfer between a single donor and a single acceptor.

We extend the sensitivity of fluorescence resonance energy transfer (FRET) to the single molecule level by measuring energy transfer between a single donor fluorophore and a single acceptor fluorophore. Near-field scanning optical microscopy (NSOM) is used to obtain simultaneous dual color images and emission spectra from donor and acceptor fluorophores linked by a short DNA molecule. Photodestruction dynamics of the donor or acceptor are used to determine the presence and efficiency of energy transfer. The classical equations used to measure energy transfer on ensembles of fluorophores are modified for single-molecule measurements. In contrast to ensemble measurements, dynamic events on a molecular scale are observable in single pair FRET measurements because they are not canceled out by random averaging. Monitoring conformational changes, such as rotations and distance changes on a nanometer scale, within single biological macromolecules, may be possible with single pair FRET.

Why so few students do maths and science. CEPS Policy Brief No. 313, 10 January 2014

Policy-makers often fret about the low number of university graduates in the fields of science, technology, engineering and mathematics (STEM). Proposed solutions often focus on providing better information for students and parents about the employability or average wages of different fields to emphasise that STEM professions pay. This paper argues that, from a personal point of view, students are actually making rational decisions, if all benefits and costs are factored into the equation. The authors conclude, therefore, that public policy needs to change the incentives to induce students to enter these fields and not just provide information about them.

Enseignement et recherche

Le vote du 9 février sur l'immigration a eu des répercussions sur les programmes de formation et de recherche académique, Eramsus+ et Horizon 2020. - L'initiative « Protection contre la sexualisation à l'école maternelle et primaire » a récolté le nombre de signatures requises. - La question de l'équivalence des diplômes de la formation professionnelle supérieure a fait débat. - La loi sur la formation continue a abouti. - Le contre-projet indirect à l'initative sur les bourses d'étude propose une réforme de la loi sur les aides à la formation. - La participation de la Suisse à la « Source européenne de spallation ESS » à passer le cap du premier Conseil.

The brook /

For voice and piano.

Œuvres complètes de Voltaire.

["Vie de Voltaire, par Condorcet": v. 1, p. 1-30.]