Methods for the analysis of multi-scale cell dynamics from fluorescence microscopy images

La embriogénesis es el proceso mediante el cual una célula se convierte en un ser un vivo. A lo largo de diferentes etapas de desarrollo, la población de células va proliferando a la vez que el embrión va tomando forma y se configura. Esto es posible gracias a la acción de varios procesos genéticos, bioquímicos y mecánicos que interaccionan y se regulan entre ellos formando un sistema complejo que se organiza a diferentes escalas espaciales y temporales. Este proceso ocurre de manera robusta y reproducible, pero también con cierta variabilidad que permite la diversidad de individuos de una misma especie. La aparición de la microscopía de fluorescencia, posible gracias a proteínas fluorescentes que pueden ser adheridas a las cadenas de expresión de las células, y los avances en la física óptica de los microscopios han permitido observar este proceso de embriogénesis in-vivo y generar secuencias de imágenes tridimensionales de alta resolución espacio-temporal. Estas imágenes permiten el estudio de los procesos de desarrollo embrionario con técnicas de análisis de imagen y de datos, reconstruyendo dichos procesos para crear la representación de un embrión digital. Una de las más actuales problemáticas en este campo es entender los procesos mecánicos, de manera aislada y en interacción con otros factores como la expresión genética, para que el embrión se desarrolle. Debido a la complejidad de estos procesos, estos problemas se afrontan mediante diferentes técnicas y escalas específicas donde, a través de experimentos, pueden hacerse y confrontarse hipótesis, obteniendo conclusiones sobre el funcionamiento de los mecanismos estudiados. Esta tesis doctoral se ha enfocado sobre esta problemática intentando mejorar las metodologías del estado del arte y con un objetivo específico: estudiar patrones de deformación que emergen del movimiento organizado de las células durante diferentes estados del desarrollo del embrión, de manera global o en tejidos concretos. Estudios se han centrado en la mecánica en relación con procesos de señalización o interacciones a nivel celular o de tejido. En este trabajo, se propone un esquema para generalizar el estudio del movimiento y las interacciones mecánicas que se desprenden del mismo a diferentes escalas espaciales y temporales. Esto permitiría no sólo estudios locales, si no estudios sistemáticos de las escalas de interacción mecánica dentro de un embrión. Por tanto, el esquema propuesto obvia las causas de generación de movimiento (fuerzas) y se centra en la cuantificación de la cinemática (deformación y esfuerzos) a partir de imágenes de forma no invasiva. Hoy en día las dificultades experimentales y metodológicas y la complejidad de los sistemas biológicos impiden una descripción mecánica completa de manera sistemática. Sin embargo, patrones de deformación muestran el resultado de diferentes factores mecánicos en interacción con otros elementos dando lugar a una organización mecánica, necesaria para el desarrollo, que puede ser cuantificado a partir de la metodología propuesta en esta tesis. La metodología asume un medio continuo descrito de forma Lagrangiana (en función de las trayectorias de puntos materiales que se mueven en el sistema en lugar de puntos espaciales) de la dinámica del movimiento, estimado a partir de las imágenes mediante métodos de seguimiento de células o de técnicas de registro de imagen. Gracias a este esquema es posible describir la deformación instantánea y acumulada respecto a un estado inicial para cualquier dominio del embrión. La aplicación de esta metodología a imágenes 3D + t del pez zebra sirvió para desvelar estructuras mecánicas que tienden a estabilizarse a lo largo del tiempo en dicho embrión, y que se organizan a una escala semejante al del mapa de diferenciación celular y con indicios de correlación con patrones de expresión genética. También se aplicó la metodología al estudio del tejido amnioserosa de la Drosophila (mosca de la fruta) durante el cierre dorsal, obteniendo indicios de un acoplamiento entre escalas subcelulares, celulares y supracelulares, que genera patrones complejos en respuesta a la fuerza generada por los esqueletos de acto-myosina. En definitiva, esta tesis doctoral propone una estrategia novedosa de análisis de la dinámica celular multi-escala que permite cuantificar patrones de manera inmediata y que además ofrece una representación que reconstruye la evolución de los procesos como los ven las células, en lugar de como son observados desde el microscopio. Esta metodología por tanto permite nuevas formas de análisis y comparación de embriones y tejidos durante la embriogénesis a partir de imágenes in-vivo. ABSTRACT The embryogenesis is the process from which a single cell turns into a living organism. Through several stages of development, the cell population proliferates at the same time the embryo shapes and the organs develop gaining their functionality. This is possible through genetic, biochemical and mechanical factors that are involved in a complex interaction of processes organized in different levels and in different spatio-temporal scales. The embryogenesis, through this complexity, develops in a robust and reproducible way, but allowing variability that makes possible the diversity of living specimens. The advances in physics of microscopes and the appearance of fluorescent proteins that can be attached to expression chains, reporting about structural and functional elements of the cell, have enabled for the in-vivo observation of embryogenesis. The imaging process results in sequences of high spatio-temporal resolution 3D+time data of the embryogenesis as a digital representation of the embryos that can be further analyzed, provided new image processing and data analysis techniques are developed. One of the most relevant and challenging lines of research in the field is the quantification of the mechanical factors and processes involved in the shaping process of the embryo and their interactions with other embryogenesis factors such as genetics. Due to the complexity of the processes, studies have focused on specific problems and scales controlled in the experiments, posing and testing hypothesis to gain new biological insight. However, methodologies are often difficult to be exported to study other biological phenomena or specimens. This PhD Thesis is framed within this paradigm of research and tries to propose a systematic methodology to quantify the emergent deformation patterns from the motion estimated in in-vivo images of embryogenesis. Thanks to this strategy it would be possible to quantify not only local mechanisms, but to discover and characterize the scales of mechanical organization within the embryo. The framework focuses on the quantification of the motion kinematics (deformation and strains), neglecting the causes of the motion (forces), from images in a non-invasive way. Experimental and methodological challenges hamper the quantification of exerted forces and the mechanical properties of tissues. However, a descriptive framework of deformation patterns provides valuable insight about the organization and scales of the mechanical interactions, along the embryo development. Such a characterization would help to improve mechanical models and progressively understand the complexity of embryogenesis. This framework relies on a Lagrangian representation of the cell dynamics system based on the trajectories of points moving along the deformation. This approach of analysis enables the reconstruction of the mechanical patterning as experienced by the cells and tissues. Thus, we can build temporal profiles of deformation along stages of development, comprising both the instantaneous events and the cumulative deformation history. The application of this framework to 3D + time data of zebrafish embryogenesis allowed us to discover mechanical profiles that stabilized through time forming structures that organize in a scale comparable to the map of cell differentiation (fate map), and also suggesting correlation with genetic patterns. The framework was also applied to the analysis of the amnioserosa tissue in the drosophila’s dorsal closure, revealing that the oscillatory contraction triggered by the acto-myosin network organized complexly coupling different scales: local force generation foci, cellular morphology control mechanisms and tissue geometrical constraints. In summary, this PhD Thesis proposes a theoretical framework for the analysis of multi-scale cell dynamics that enables to quantify automatically mechanical patterns and also offers a new representation of the embryo dynamics as experienced by cells instead of how the microscope captures instantaneously the processes. Therefore, this framework enables for new strategies of quantitative analysis and comparison between embryos and tissues during embryogenesis from in-vivo images.

Force-mediated kinetics of single P-selectin/ligand complexes observed by atomic force microscopy

Leukocytes roll along the endothelium of postcapillary venules in response to inflammatory signals. Rolling under the hydrodynamic drag forces of blood flow is mediated by the interaction between selectins and their ligands across the leukocyte and endothelial cell surfaces. Here we present force-spectroscopy experiments on single complexes of P-selectin and P-selectin glycoprotein ligand-1 by atomic force microscopy to determine the intrinsic molecular properties of this dynamic adhesion process. By modeling intermolecular and intramolecular forces as well as the adhesion probability in atomic force microscopy experiments we gain information on rupture forces, elasticity, and kinetics of the P-selectin/P-selectin glycoprotein ligand-1 interaction. The complexes are able to withstand forces up to 165 pN and show a chain-like elasticity with a molecular spring constant of 5.3 pN nm−1 and a persistence length of 0.35 nm. The dissociation constant (off-rate) varies over three orders of magnitude from 0.02 s−1 under zero force up to 15 s−1 under external applied forces. Rupture force and lifetime of the complexes are not constant, but directly depend on the applied force per unit time, which is a product of the intrinsic molecular elasticity and the external pulling velocity. The high strength of binding combined with force-dependent rate constants and high molecular elasticity are tailored to support physiological leukocyte rolling.

RecA binding to a single double-stranded DNA molecule: A possible role of DNA conformational fluctuations

Most genetic regulatory mechanisms involve protein–DNA interactions. In these processes, the classical Watson–Crick DNA structure sometimes is distorted severely, which in turn enables the precise recognition of the specific sites by the protein. Despite its key importance, very little is known about such deformation processes. To address this general question, we have studied a model system, namely, RecA binding to double-stranded DNA. Results from micromanipulation experiments indicate that RecA binds strongly to stretched DNA; based on this observation, we propose that spontaneous thermal stretching fluctuations may play a role in the binding of RecA to DNA. This has fundamental implications for the protein–DNA binding mechanism, which must therefore rely in part on a combination of flexibility and thermal fluctuations of the DNA structure. We also show that this mechanism is sequence sensitive. Theoretical simulations support this interpretation of our experimental results, and it is argued that this is of broad relevance to DNA–protein interactions.

Fluorescence and photobleaching dynamics of single light-harvesting complexes

Single light-harvesting complexes LH-2 from Rhodopseudomonas acidophila were immobilized on various charged surfaces under physiological conditions. Polarized light experiments showed that the complexes were situated on the surface as nearly upright cylinders. Their fluorescence lifetimes and photobleaching properties were obtained by using a confocal fluorescence microscope with picosecond time resolution. Initially all molecules fluoresced with a lifetime of 1 ± 0.2 ns, similar to the bulk value. The photobleaching of one bacteriochlorophyll molecule from the 18-member assembly caused the fluorescence to switch off completely, because of trapping of the mobile excitations by energy transfer. This process was linear in light intensity. On continued irradiation the fluorescence often reappeared, but all molecules did not show the same behavior. Some LH-2 complexes displayed a variation of their quantum yields that was attributed to photoinduced confinement of the excited states and thereby a diminution of the superradiance. Others showed much shorter lifetimes caused by excitation energy traps that are only ≈3% efficient. On repeated excitation some molecules entered a noisy state where the fluorescence switched on and off with a correlation time of ≈0.1 s. About 490 molecules were examined.

Controlled unzipping of a bacterial surface layer with atomic force microscopy

We have combined high-resolution atomic force microscopy (AFM) imaging and force spectroscopy to gain insight into the interaction forces between the individual protomers of the hexagonally packed intermediate (HPI) layer of Deinococcus radiodurans. After imaging the HPI layer, the AFM stylus was attached to individual protomers by enforced stylus-sample contact to allow force spectroscopy experiments. Imaging of the HPI layer after recording force-extension curves allowed adhesion forces to be correlated with structural alterations. By using this approach, individual protomers of the HPI layer were found to be removed at pulling forces of ≈300 pN. Furthermore, it was possible to sequentially unzip entire bacterial pores formed by six HPI protomers. The combination of high-resolution AFM imaging of individual proteins with the determination of their intramolecular forces is a method of studying the mechanical stability of supramolecular structures at the level of single molecules.

Three-dimensional structure of human electron transfer flavoprotein to 2.1-Å resolution

Mammalian electron transfer flavoproteins (ETF) are heterodimers containing a single equivalent of flavin adenine dinucleotide (FAD). They function as electron shuttles between primary flavoprotein dehydrogenases involved in mitochondrial fatty acid and amino acid catabolism and the membrane-bound electron transfer flavoprotein ubiquinone oxidoreductase. The structure of human ETF solved to 2.1-Å resolution reveals that the ETF molecule is comprised of three distinct domains: two domains are contributed by the α subunit and the third domain is made up entirely by the β subunit. The N-terminal portion of the α subunit and the majority of the β subunit have identical polypeptide folds, in the absence of any sequence homology. FAD lies in a cleft between the two subunits, with most of the FAD molecule residing in the C-terminal portion of the α subunit. Alignment of all the known sequences for the ETF α subunits together with the putative FixB gene product shows that the residues directly involved in FAD binding are conserved. A hydrogen bond is formed between the N5 of the FAD isoalloxazine ring and the hydroxyl side chain of αT266, suggesting why the pathogenic mutation, αT266M, affects ETF activity in patients with glutaric acidemia type II. Hydrogen bonds between the 4′-hydroxyl of the ribityl chain of FAD and N1 of the isoalloxazine ring, and between αH286 and the C2-carbonyl oxygen of the isoalloxazine ring, may play a role in the stabilization of the anionic semiquinone. With the known structure of medium chain acyl-CoA dehydrogenase, we hypothesize a possible structure for docking the two proteins.

Determination of the lifetime and force dependence of interactions of single bonds between surface-attached CD2 and CD48 adhesion molecules

We studied single molecular interactions between surface-attached rat CD2, a T-lymphocyte adhesion receptor, and CD48, a CD2 ligand found on antigen-presenting cells. Spherical particles were coated with decreasing densities of CD48–CD4 chimeric molecules then driven along CD2-derivatized glass surfaces under a low hydrodynamic shear rate. Particles exhibited multiple arrests of varying duration. By analyzing the dependence of arrest frequency and duration on the surface density of CD48 sites, it was concluded that (i) arrests were generated by single molecular bonds and (ii) the initial bond dissociation rate was about 7.8 s−1. The force exerted on bonds was increased from about 11 to 22 pN; the detachment rate exhibited a twofold increase. These results agree with and extend studies on the CD2–CD48 interaction by surface plasmon resonance technology, which yielded an affinity constant of ≈104 M−1 and a dissociation rate of ≥6 s−1. It is concluded that the flow chamber technology can be an useful complement to atomic force microscopy for studying interactions between isolated biomolecules, with a resolution of about 20 ms and sensitivity of a few piconewtons. Further, this technology might be extended to actual cells.

Real-time visualization of intracellular hydrodynamics in single living cells

Intracellular water concentrations in single living cells were visualized by nonlinear coherent anti-Stokes Raman scattering (CARS) microscopy. In combination with isotopic exchange measurements, CARS microscopy allowed the real-time observation of transient intracellular hydrodynamics at a high spatial resolution. Studies of the hydrodynamics in the microorganism Dictyostelium discoideum indicated the presence of a microscopic region near the plasma membrane where the mobility of water molecules is severely restricted. Modeling the transient hydrodynamics eventuated in the determination of cell-specific cytosolic diffusion and plasma membrane permeability constants. Our experiments demonstrate that CARS microscopy offers an invaluable tool for probing single-cell water dynamics.

Imaging ROMK1 inwardly rectifying ATP-sensitive K+ channel protein using atomic force microscopy.

The inwardly rectifying K+ channel ROMK1 has been implicated as being significant in K+ secretion in the distal nephron. ROMK1 has been shown by immunocytochemistry to be expressed in relevant nephron segments. The development of the atomic force microscope has made possible the production of high resolution images of small particles, including a variety of biological macromolecules. Recently, a fusion protein of glutathione S-transferase (GST) and ROMK1 (ROMK1-GST) has been used to produce a polyclonal antibody for immunolocalization of ROMK1. We have used atomic force microscopy to examine ROMK1-GST and the native ROMK1 polypeptide cleaved from GST. Imaging was conducted with the proteins in physiological solutions attached to mica. ROMK1-GST appears in images as a particle composed of two units of similar size. Analyses of images indicate that the two units have volumes of approximately 118 nm3, which is close to the theoretical volume of a globular protein of approximately 65 kDa (the molecular mass of ROMK1-GST). Native GST exists as a dimer, and the images obtained here are consistent with the ROMK1-GST fusion protein's existence as a heterodimer. In experiments on ROMK1 in aqueous solution, single molecules appear to aggregate, but contact to the mica was maintained. Addition of ATP to the solution produced a change in height of the aggregates. This change (which was reversible) suggests that ATP induces a structural change in the ROMK1 protein. The data show that atomic force microscopy is a useful tool for examination of purified protein molecules under near-physiological conditions, and furthermore, that structural alterations in the proteins may be continuously investigated.

Direct observation of iron-induced conformational changes of mitochondrial DNA by high-resolution field-emission in-lens scanning electron microscopy.

When respiring rat liver mitochondria are incubated in the presence of Fe(III) gluconate, their DNA (mtDNA) relaxes from the supercoiled to the open circular form dependent on the iron dose. Anaerobiosis or antioxidants fail to completely inhibit the unwinding. High-resolution field-emission in-lens scanning electron microscopy imaging, in concert with backscattered electron detection, pinpoints nanometer-range iron colloids bound to mtDNA isolated from iron-exposed mitochondria. High-resolution field-emission in-lens scanning electron microscopy with backscattered electron detection imaging permits simultaneous detailed visual analysis of DNA topology, iron dose-dependent mtDNA unwinding, and assessment of iron colloid formation on mtDNA strands.

Detection and localization of individual antibody-antigen recognition events by atomic force microscopy.

A methodology has been developed for the study of molecular recognition at the level of single events and for the localization of sites on biosurfaces, in combining force microscopy with molecular recognition by specific ligands. For this goal, a sensor was designed by covalently linking an antibody (anti-human serum albumin, polyclonal) via a flexible spacer to the tip of a force microscope. This sensor permitted detection of single antibody-antigen recognition events by force signals of unique shape with an unbinding force of 244 +/- 22 pN. Analysis revealed that observed unbinding forces originate from the dissociation of individual Fab fragments from a human serum albumin molecule. The two Fab fragments of the antibody were found to bind independently and with equal probability. The flexible linkage provided the antibody with a 6-nm dynamical reach for binding, rendering binding probability high, 0.5 for encounter times of 60 ms. This permitted fast and reliable detection of antigenic sites during lateral scans with a positional accuracy of 1.5 nm. It is indicated that this methodology has promise for characterizing rate constants and kinetics of molecular recognition complexes and for molecular mapping of biosurfaces such as membranes.

Mapping nucleosome position at single base-pair resolution by using site-directed hydroxyl radicals.

A base-pair resolution method for determining nucleosome position in vitro has been developed to com- plement existing, less accurate methods. Cysteaminyl EDTA was tethered to a recombinant histone octamer via a mutant histone H4 with serine 47 replaced by cysteine. When assembled into nucleosome core particles, the DNA could be cut site specifically by hydroxyl radical-catalyzed chain scission by using the Fenton reaction. Strand cleavage occurs mainly at a single nucleotide close to the dyad axis of the core particle, and assignment of this location via the symmetry of the nucleosome allows base-pair resolution mapping of the histone octamer position on the DNA. The positions of the histone octamer and H3H4 tetramer were mapped on a 146-bp Lytechinus variegatus 5S rRNA sequence and a twofold-symmetric derivative. The weakness of translational determinants of nucleosome positioning relative to the overall affinity of the histone proteins for this DNA is clearly demonstrated. The predominant location of both histone octamer and H3H4 tetramer assembled on the 5S rDNA is off center. Shifting the nucleosome core particle position along DNA within a conserved rotational phase could be induced under physiologically relevant conditions. Since nucleosome shifting has important consequences for chromatin structure and gene regulation, an approach to the thermodynamic characterization of this movement is proposed. This mapping method is potentially adaptable for determining nucleosome position in chromatin in vivo.

Expression of a single-chain HLA class I molecule in a human cell line: presentation of exogenous peptide and processed antigen to cytotoxic T lymphocytes.

We have synthesized a recombinant gene encoding a single-chain HLA-A2/beta 2-microglobulin (beta 2m) molecule by linking beta 2m through its carboxyl terminus via a short peptide spacer to HLA-A2 (A*0201). This gene has been expressed in the beta 2m-deficient colorectal tumor cell line DLD-1. Transfection of this cell with the single-chain construct was associated with conformationally correct cell surface expression of a class I molecule of appropriate molecular mass. The single-chain HLA class I molecule presented either exogenously added peptide or (after interferon-gamma treatment) endogenously processed antigen to an influenza A matrix-specific, HLA-A2-restricted cytotoxic T-lymphocyte line. The need for interferon gamma for the processing and presentation of endogenous antigen suggests that DLD-1 has an antigen-processing defect that can be up-regulated, a feature that may be found in other carcinomas. Our data indicate that single-chain HLA class I constructs can form functional class I molecules capable of presenting endogenously processed antigens. Such molecules should be of use for functional studies, as well as providing potential anticancer immunotherapeutic agents or vaccines.

A small bispecific antibody construct expressed as a functional single-chain molecule with high tumor cell cytotoxicity.

Construction of a bispecific single-chain antibody derivative is described that consists of two different single-chain Fv fragments joined through a Gly-Ser linker. One specificity of the two Fv fragments is directed against the CD3 antigen of human T cells and the other is directed against the epithelial 17-1A antigen; the latter had been found in a clinical trial to be a suitable target for antibody therapy of minimal residual colorectal cancer. The construct could be expressed in CHO cells as a fully functional protein, while its periplasmic expression in Escherichia coli resulted in a nonfunctional protein only. The antigen-binding properties of the bispecific single-chain antibody are indistinguishable from those of the corresponding univalent single-chain Fv fragments. By redirecting human peripheral T lymphocytes against 17-1A-positive tumor cells, the bispecific antibody proved to be highly cytotoxic at nanomolar concentrations as demonstrated by 51Cr release assay on various cell lines. The described bispecific construct has a molecular mass of 60 kDa and can be easily purified by its C-terminal histidine tail on a Ni-NTA chromatography column. As bispecific antibodies have already been shown to be effective in vivo in experimental tumor systems as well as in phase-one clinical trials, the small CD3/17-1A-bispecific antibody may be more efficacious than intact antibodies against minimal residual cancer cells.