Discarge compilation from The Global River Discharge (RivDIS) Project

The Global River Discharge (RivDIS) data set contains monthly discharge measurements for 1018 stations located throughout the world. The period of record varies widely from station to station, with a mean of 21.5 years. These data were digitized from published UNESCO archives by Charles Voromarty, Balaze Fekete, and B.A. Tucker of the Complex Systems Research Center (CSRC) at the University of New Hampshire. River discharge is typically measured through the use of a rating curve that relates local water level height to discharge. This rating curve is used to estimate discharge from the observed water level. The rating curves are periodically rechecked and recalibrated through on-site measurement of discharge and river stage.

Segmentación de la glotis en imágenes laríngeas usando snakes.

El presente trabajo describe una nueva metodología para la detección automática del espacio glotal de imágenes laríngeas tomadas a partir de 15 vídeos grabados por el servicio ORL del hospital Gregorio Marañón de Madrid con luz estroboscópica. El sistema desarrollado está basado en el modelo de contornos activos (snake). El algoritmo combina en el pre-procesado, algunas técnicas tradicionales (umbralización y filtro de mediana) con técnicas más sofisticadas tales como filtrado anisotrópico. De esta forma, se obtiene una imagen apropiada para el uso de las snakes. El valor escogido para el umbral es del 85% del pico máximo del histograma de la imagen; sobre este valor la información de los píxeles no es relevante. El filtro anisotrópico permite distinguir dos niveles de intensidad, uno es el fondo y el otro es la glotis. La inicialización se basa en obtener el módulo del campo GVF; de esta manera se asegura un proceso automático para la selección del contorno inicial. El rendimiento del algoritmo se valida usando los coeficientes de Pratt y se compara contra una segmentación realizada manualmente y otro método automático basado en la transformada de watershed. SUMMARY: The present work describes a new methodology for the automatic detection of the glottal space from laryngeal images taken from 15 videos recorded by the ENT service of the Gregorio Marañon Hospital in Madrid with videostroboscopic equipment. The system is based on active contour models (snakes). The algorithm combines for the pre-processing, some traditional techniques (thresholding and median filter) with more sophisticated techniques such as anisotropic filtering. In this way, we obtain an appropriate image for the use of snake. The value selected for the threshold is 85% of the maximum peak of the image histogram; over this point the information of the pixels is not relevant. The anisotropic filter permits to distinguish two intensity levels, one is the background and the other one is the glottis. The initialization is based on the obtained magnitude by GVF field; in this manner an automatic process for the initial contour selection will be assured. The performance of the algorithm is tested using the Pratt coefficient and compared against a manual segmentation and another automatic method based on the watershed transformation.

The gastrulation defective gene of Drosophila melanogaster is a member of the serine protease superfamily

The establishment of dorsal–ventral polarity in the oocyte involves two sets of genes. One set belongs to the gurken-torpedo signaling pathway and affects the development of the egg chorion as well as the polarity of the embryo. The second set of genes affects only the dorsal–ventral polarity of the embryo but not the eggshell. gastrulation defective is one of the earliest acting of this second set of maternally required genes. We have cloned and characterized the gastrulation defective gene and determined that it encodes a protein structurally related to the serine protease superfamily, which also includes the Snake, Easter, and Nudel proteins. These data provide additional support for the involvement of a protease cascade in generating an asymmetric signal (i.e., asymmetric Spätzle activity) during establishment of dorsal–ventral polarity in the Drosophila embryo.

Structural origins of the exonuclease resistance of a zwitterionic RNA

Nuclease resistance and RNA affinity are key criteria in the search for optimal antisense nucleic acid modifications, but the origins of the various levels of resistance to nuclease degradation conferred by chemical modification of DNA and RNA are currently not understood. The 2′-O-aminopropyl (AP)-RNA modification displays the highest nuclease resistance among all phosphodiester-based analogues and its RNA binding affinity surpasses that of phosphorothioate DNA by 1°C per modified residue. We found that oligodeoxynucleotides containing AP-RNA residues at their 3′ ends competitively inhibit the degradation of single-stranded DNA by the Escherichia coli Klenow fragment (KF) 3′-5′ exonuclease and snake venom phosphodiesterase. To shed light on the origins of nuclease resistance brought about by the AP modification, we determined the crystal structure of an A-form DNA duplex with AP-RNA modifications at 1.6-Å resolution. In addition, the crystal structures of complexes between short DNA fragments carrying AP-RNA modifications and wild-type KF were determined at resolutions between 2.2 and 3.0 Å and compared with the structure of the complex between oligo(dT) and the D355A/E357A KF mutant. The structural models suggest that interference of the positively charged 2′-O-substituent with the metal ion binding site B of the exonuclease allows AP-RNA to effectively slow down degradation.

Unblocking of chain-terminated primer by HIV-1 reverse transcriptase through a nucleotide-dependent mechanism

HIV-1 replication is inhibited by the incorporation of chain-terminating nucleotides at the 3′ end of the growing DNA chain. Here we show a nucleotide-dependent reaction catalyzed by HIV-1 reverse transcriptase that can efficiently remove the chain-terminating residue, yielding an extendible primer terminus. Radioactively labeled 3′-terminal residue from the primer can be transferred into a product that is resistant to calf intestinal alkaline phosphatase and sensitive to cleavage by snake venom phosphodiesterase. The products formed from different nucleotide substrates have unique electrophoretic migrations and have been identified as dinucleoside tri- or tetraphosphates. The reaction is inhibited by dNTPs that are complementary to the next position on the template (Ki ≈ 5 μM), suggesting competition between dinucleoside polyphosphate synthesis and DNA polymerization. Dinucleoside polyphosphate synthesis was inhibited by an HIV-1 specific non-nucleoside inhibitor and was absent in mutant HIV-1 reverse transcriptase deficient in polymerase activity, indicating that this activity requires a functional polymerase active site. We suggest that dinucleoside polyphosphate synthesis occurs by transfer of the 3′ nucleotide from the primer to the pyrophosphate moiety in the nucleoside di- or triphosphate substrate through a mechanism analogous to pyrophosphorolysis. Unlike pyrophosphorolysis, however, the reaction is nucleotide-dependent, is resistant to pyrophosphatase, and produces dinucleoside polyphosphates. Because it occurs at physiological concentrations of ribonucleoside triphosphates, this reaction may determine the in vivo activity of many nucleoside antiretroviral drugs.

Unusual horizontal transfer of a long interspersed nuclear element between distant vertebrate classes

We have shown previously by Southern blot analysis that Bov-B long interspersed nuclear elements (LINEs) are present in different Viperidae snake species. To address the question as to whether Bov-B LINEs really have been transmitted horizontally between vertebrate classes, the analysis has been extended to a larger number of vertebrate, invertebrate, and plant species. In this paper, the evolutionary origin of Bov-B LINEs is shown unequivocally to be in Squamata. The previously proposed horizontal transfer of Bov-B LINEs in vertebrates has been confirmed by their discontinuous phylogenetic distribution in Squamata (Serpentes and two lizard infra-orders) as well as in Ruminantia, by the high level of nucleotide identity, and by their phylogenetic relationships. The horizontal transfer of Bov-B LINEs from Squamata to the ancestor of Ruminantia is evident from the genetic distances and discontinuous phylogenetic distribution. The ancestor of Colubroidea snakes is a possible donor of Bov-B LINEs to Ruminantia. The timing of horizontal transfer has been estimated from the distribution of Bov-B LINEs in Ruminantia and the fossil data of Ruminantia to be 40–50 My ago. The phylogenetic relationships of Bov-B LINEs from the various Squamata species agrees with that of the species phylogeny, suggesting that Bov-B LINEs have been maintained stably by vertical transmission since the origin of Squamata in the Mesozoic era.

Activation of a protease cascade involved in patterning the Drosophila embryo

Dorsoventral patterning of the Drosophila embryo is initiated by a ventralizing signal. Production of this signal requires the serine proteases Gastrulation Defective (GD), Snake, and Easter, which genetic studies suggest act sequentially in a cascade that is activated locally in response to a ventral cue provided by the pipe gene. Here, we demonstrate biochemically that GD activates Snake, which in turn activates Easter. We also provide evidence that GD zymogen cleavage is important for triggering this cascade but is not spatially localized by pipe. Our results suggest that a broadly, rather than locally, activated protease cascade produces the ventralizing signal, so a distinct downstream step in this cascade must be spatially regulated to restrict signaling to the ventral side of the embryo.

A β-hairpin structure in a 13-mer peptide that binds α-bungarotoxin with high affinity and neutralizes its toxicity

Snake-venom α-bungarotoxin is a member of the α-neurotoxin family that binds with very high affinity to the nicotinic acetylcholine receptor (AChR) at the neuromuscular junction. The structure of the complex between α-bungarotoxin and a 13-mer peptide (WRYYESSLEPYPD) that binds the toxin with high affinity, thus inhibiting its interactions with AChR with an IC50 of 2 nM, has been solved by 1H-NMR spectroscopy. The bound peptide folds into a β-hairpin structure created by two antiparallel β-strands, which combine with the already existing triple-stranded β-sheet of the toxin to form a five-stranded intermolecular, antiparallel β-sheet. Peptide residues Y3P, E5P, and L8P have the highest intermolecular contact area, indicating their importance in the binding of α-bungarotoxin; W1P, R2P, and Y4P also contribute significantly to the binding. A large number of characteristic hydrogen bonds and electrostatic and hydrophobic interactions are observed in the complex. The high-affinity peptide exhibits inhibitory potency that is better than any known peptide derived from AChR, and is equal to that of the whole α-subunit of AChR. The high degree of sequence similarity between the peptide and various types of AChRs implies that the binding mode found within the complex might possibly mimic the receptor binding to the toxin. The design of the high-affinity peptide was based on our previous findings: (i) the detection of a lead peptide (MRYYESSLKSYPD) that binds α-bungarotoxin, using a phage-display peptide library, (ii) the information about the three-dimensional structure of α-bungarotoxin/lead-peptide complex, and (iii) the amino acid sequence analysis of different AChRs.

Crystal structure of an anticoagulant protein in complex with the Gla domain of factor X

The γ-carboxyglutamic acid (Gla) domain of blood coagulation factors is responsible for Ca2+-dependent phospholipid membrane binding. Factor X-binding protein (X-bp), an anticoagulant protein from snake venom, specifically binds to the Gla domain of factor X. The crystal structure of X-bp in complex with the Gla domain peptide of factor X at 2.3-Å resolution showed that the anticoagulation is based on the fact that two patches of the Gla domain essential for membrane binding are buried in the complex formation. The Gla domain thus is expected to be a new target of anticoagulant drugs, and X-bp provides a basis for designing them. This structure also provides a membrane-bound model of factor X.

Two subsites in the binding domain of the acetylcholine receptor: an aromatic subsite and a proline subsite.

The ligand binding site of the nicotinic acetylcholine receptor (AcChoR) is localized in the alpha-subunit within a domain containing the tandem Cys-192 and -193. By analyzing the binding-site region of AcChoR from animal species that are resistant to alpha-neurotoxins, we have previously shown that four residues in this region, at positions 187, 189, 194, and 197, differ between animals sensitive (e.g., mouse) and resistant (e.g., mongoose and snake) to alpha-bungarotoxin (alpha-BTX). In the present study, we performed site-directed mutagenesis on a fragment of the mongoose AcChoR alpha-subunit (residues 122-205) and exchanged residues 187, 189, 194, and 197, either alone or in combination, with those present in the mouse alpha-subunit sequence. Only the mongoose fragment in which all four residues were mutated to the mouse ones exhibited alpha-BTX binding similar to that of the mouse fragment. The mongoose double mutation in which Leu-194 and His-197 were replaced with proline residues, which are present at these positions in the mouse AcChoR and in all other toxin binders, bound alpha-BTX to approximately 60% of the level of binding exhibited by the mouse fragment. In addition, replacement of either Pro-194 or -197 in the mouse fragment with serine and histidine, respectively, markedly decreased alpha-BTX binding. All other mutations resulted in no or just a small increase in alpha-BTX binding. These results have led us to propose two subsites in the binding domain for alpha-BTX: the proline subsite, which includes Pro-194 and -197 and is critical for alpha-BTX binding, and the aromatic subsite, which includes amino acid residues 187 and 189 and determines the extent of alpha-BTX binding.

Novel microbicidal agents derived from Naja atra cardiotoxin 1 (CTX-1)

During my PhD course, I focused my research on antimicrobial peptides (AMPs), in particular on the aspects of their computational design and development. This work led to the development of a new family of AMPs that I designed, starting from the amino acid sequence of a snake venom toxin, the cardiotoxin 1 (CTX-1) of Naja atra. Naja atra atra cardiotoxin 1, produced by Chinese cobra snakes belonging to Elapidae family, is included in the three-finger toxin family and exerts high cytotoxicity and antimicrobial activity too. This toxin family is characterized by specific folding of three beta-sheet loops (“fingers”) extending from the central core and by four conserved disulfide bridges. Using as template the first loop of this toxin, different sequences of 20 amino acids linear cationic peptides have been designed in order to avoid toxic effects but to maintain and strengthen the antimicrobial activity. As a result, the sequence NCP-0 (Naja Cardiotoxin Peptide-0) was designed as ancestor and subsequently other 4 variant sequences of NCP0 were developed. These variant sequences have shown microbicidal activity towards a panel of reference strains of Gram-positive and Gram-negative bacteria, fungi and an enveloped virus. In particular, the sequence designed as NCP-3 (Naja Cardiotoxin Peptide-3) and its variants NCP-3a and NCP-3b have shown the best antimicrobial activity together with low cytotoxicity against eukaryotic cells and low hemolytic activity. Bactericidal activity has been demonstrated by minimum bactericidal concentration (MBC) assay at values below 10 μg/ml for Pseudomonas aeruginosa ATCC 27853, Acinetobacter baumannii ( clinical isolates), Moraxella catharralis ATCC 25238, MRSA ATCC 43400, while towards Staphylococcus aureus ATCC 25923, Enterococcus hirae ATCC 10541 and Streptococcus agalactiae ATCC 13813 the bactericidal activity was demonstrated even below 1.6 μg/ml concentration. This potent antimicrobial activity was confirmed even for unicellular fungi Candida albicans, Candida glabrata and Malassezia pachydermatis (MBC 32.26-6.4 μg/ml), and also against the fast-growing mycobacteria Mycobacterium smegmatis DSMZ 43756 and Mycobacterium fortuitum DSMZ 46621 (MBC 100 μg/ml). Moreover, NCP-3 has shown a virucidal activity on the enveloped virus Bovine Herpesvirus 1 (BoHV1) belonging to herpesviridae family. The bactericidal activity is maintained in a high salt concentration (125 and 250 mM NaCl) medium and PB +20% Mueller Hinton Medium for E. coli, MRSA and Pseudomonas aeruginosa reference strains. Considering these in vitro obtained data, we propose NCP-3 and its variants NCP-3a and NCP-3b as promising antimicrobial candidates. For this reason, the whole novel AMPs family has been protected by a national patent (n°102015000015951).

Video of the labelling of projected Placido rings through the scanning annulus

The Elliptical Scanning Algorithm is an effective method to individually detect and label the projected rings. It consecutively defines an elliptical annulus of one pixel wide which grows pixel by pixel and sweeps the image, from centre to periphery, until it detects and labels each whole ring. In a way, it works like a snake-annealing algorithm. Active contour models (snakes) are energy-minimising curves that deform to fit image features. Elliptical Scanning Algorithm changes its geometry in order to label reflected rings.

Corneal topography reinterpretation through separate analysis of the projected rings

Póster presentado en SPIE Photonics Europe, Brussels, 16-19 April 2012.

Snakes based tracking and texture analysis of microscopic images

Comunicación presentada en el VII Symposium Nacional de Reconocimiento de Formas y Análisis de Imágenes, SNRFAI, Barcelona, abril 1997.