Regularization for sparsity in statistical analysis and machine learning

Pragmatism is the leading motivation of regularization. We can understand regularization as a modification of the maximum-likelihood estimator so that a reasonable answer could be given in an unstable or ill-posed situation. To mention some typical examples, this happens when fitting parametric or non-parametric models with more parameters than data or when estimating large covariance matrices. Regularization is usually used, in addition, to improve the bias-variance tradeoff of an estimation. Then, the definition of regularization is quite general, and, although the introduction of a penalty is probably the most popular type, it is just one out of multiple forms of regularization. In this dissertation, we focus on the applications of regularization for obtaining sparse or parsimonious representations, where only a subset of the inputs is used. A particular form of regularization, L1-regularization, plays a key role for reaching sparsity. Most of the contributions presented here revolve around L1-regularization, although other forms of regularization are explored (also pursuing sparsity in some sense). In addition to present a compact review of L1-regularization and its applications in statistical and machine learning, we devise methodology for regression, supervised classification and structure induction of graphical models. Within the regression paradigm, we focus on kernel smoothing learning, proposing techniques for kernel design that are suitable for high dimensional settings and sparse regression functions. We also present an application of regularized regression techniques for modeling the response of biological neurons. Supervised classification advances deal, on the one hand, with the application of regularization for obtaining a na¨ıve Bayes classifier and, on the other hand, with a novel algorithm for brain-computer interface design that uses group regularization in an efficient manner. Finally, we present a heuristic for inducing structures of Gaussian Bayesian networks using L1-regularization as a filter. El pragmatismo es la principal motivación de la regularización. Podemos entender la regularización como una modificación del estimador de máxima verosimilitud, de tal manera que se pueda dar una respuesta cuando la configuración del problema es inestable. A modo de ejemplo, podemos mencionar el ajuste de modelos paramétricos o no paramétricos cuando hay más parámetros que casos en el conjunto de datos, o la estimación de grandes matrices de covarianzas. Se suele recurrir a la regularización, además, para mejorar el compromiso sesgo-varianza en una estimación. Por tanto, la definición de regularización es muy general y, aunque la introducción de una función de penalización es probablemente el método más popular, éste es sólo uno de entre varias posibilidades. En esta tesis se ha trabajado en aplicaciones de regularización para obtener representaciones dispersas, donde sólo se usa un subconjunto de las entradas. En particular, la regularización L1 juega un papel clave en la búsqueda de dicha dispersión. La mayor parte de las contribuciones presentadas en la tesis giran alrededor de la regularización L1, aunque también se exploran otras formas de regularización (que igualmente persiguen un modelo disperso). Además de presentar una revisión de la regularización L1 y sus aplicaciones en estadística y aprendizaje de máquina, se ha desarrollado metodología para regresión, clasificación supervisada y aprendizaje de estructura en modelos gráficos. Dentro de la regresión, se ha trabajado principalmente en métodos de regresión local, proponiendo técnicas de diseño del kernel que sean adecuadas a configuraciones de alta dimensionalidad y funciones de regresión dispersas. También se presenta una aplicación de las técnicas de regresión regularizada para modelar la respuesta de neuronas reales. Los avances en clasificación supervisada tratan, por una parte, con el uso de regularización para obtener un clasificador naive Bayes y, por otra parte, con el desarrollo de un algoritmo que usa regularización por grupos de una manera eficiente y que se ha aplicado al diseño de interfaces cerebromáquina. Finalmente, se presenta una heurística para inducir la estructura de redes Bayesianas Gaussianas usando regularización L1 a modo de filtro.

The Optimal combination: Grammatical Swarm, Particle Swarm Optimization and Neural Networks.

Social behaviour is mainly based on swarm colonies, in which each individual shares its knowledge about the environment with other individuals to get optimal solutions. Such co-operative model differs from competitive models in the way that individuals die and are born by combining information of alive ones. This paper presents the particle swarm optimization with differential evolution algorithm in order to train a neural network instead the classic back propagation algorithm. The performance of a neural network for particular problems is critically dependant on the choice of the processing elements, the net architecture and the learning algorithm. This work is focused in the development of methods for the evolutionary design of artificial neural networks. This paper focuses in optimizing the topology and structure of connectivity for these networks.

Comparative study of ion beam-irradiation effects on silica and α-quartz: evidences for excitonic mechanisms

Introduction - SiO 2 •Simple composition and structure; Crystalline and amorphous phases •Adequate for atomistic simulations •Abundant in nature. Relevant for many technologies -Irradiation with swift heavy ions: •They provide EXTREME physical conditions •Very high excitation densities similar to high power lasers •Very high local temperatures •By playing with high energy and heavy mass (SHI) : •One can go from low electronic excitations (collisions regime) to high electronic excitations (electronic regime

Electrodynamic Tethers For Planetary And De-orbiting Missions

Nuevas aplicaciones tecnológicas y científicas mediante amarras electrodinámicas son analizadas para misiones planetarias. i) Primero, se considera un conjunto de amarras cilíndricas en paralelo (veleros electrosolares) para una misión interplanetaria. Los iones provenientes del viento solar son repelidos por el alto potencial de dichas amarras generando empuje sobre el velero. Para conocer el intercambio de momento que provocan los iones sobre las amarras se ha considerado un modelo de potencial estacionario. Se ha analizado la transferencia orbital de la Tierra a Júpiter siguiendo un método de optimización de trayectoria indirecto. ii) Una vez que el velero se encuentra cerca de Júpiter, se ha considerado el despliegue de una amarra para diferentes objetivos científicos. iia) Una amarra podría ser utilizada para diagnóstico de plasmas, al ser una fuente efectiva de ondas, y también como un generador de auroras artificiales. Una amarra conductora que orbite en la magnetosfera jovial es capaz de producir ondas. Se han analizado las diferentes ondas radiadas por un conductor por el que circula una corriente constante que sigue una órbita polar de alta excentricidad y bajo apoápside, como ocurre en la misión Juno de la NASA. iib) Además, se ha estudiado una misión tentativa que sigue una órbita ecuatorial (LJO) por debajo de los intensos cinturones de radiación. Ambas misiones requiren potencia eléctrica para los sistemas de comunicación e instrumentos científicos. Las amarras pueden generar potencia de manera más eficiente que otros sistemas que utlizan paneles solares o sistemas de potencia de radioisótopos (RPS). La impedancia de radiación es necesaria para determinar la corriente que circula por todo el circuito de la amarra. En un modelo de plasma frío, la radiación ocurre principalmente en los modos de Alfven y magnetosónica rápida, mostrando un elevado índice de refracción. Se ha estudiado la impedancia de radiación en amarras con recubrimiento aislante para los dos modos de radiación y cada una de las misiones. A diferencia del caso ionosférico terrestre, la baja densidad y el intenso campo magnético que aparecen en el entorno de Júpiter consiguen que la girofrecuencia de los electrones sea mucho mayor que la frecuencia del plasma; esto hace que el espectro de potencia para cada modo se modifique substancialmente, aumentando la velocidad de Alfven. Se ha estimado también la impedancia de radiación para amarras sin aislante conductor. En la misión LJO, un vehículo espacial bajando lentamente la altitud de su órbita permitiría estudiar la estructura del campo magnético y composición atmosférica para entender la formación, evolución, y estructura de Júpiter. Adicionalmente, si el contactor (cátodo) se apaga, se dice que la amarra flota eléctricamente, permitiendo emisión de haz de electrones que generan auroras. El continuo apagado y encendido produce pulsos de corriente dando lugar a emisiones de señales, que pueden ser utilizadas para diagnóstico del plasma jovial. En Órbita Baja Jovial, los iones que impactan contra una amarra polarizada negativamente producen electrones secundarios, que, viajando helicoidalmente sobre las líneas de campo magnético de Júpiter, son capaces de alcanzar su atmósfera más alta, y, de esta manera, generar auroras. Se han identificado cuáles son las regiones donde la amarra sería más eficiente para producir auroras. iic) Otra aplicación científica sugerida para la misión LJO es la detección de granos cargados que orbitan cerca de Júpiter. Los electrones de alta energía en este ambiente pueden ser modelados por una distribucción no Maxwelliana conocida como distribución kappa. En escenarios con plasmas complejos, donde los campos eléctricos en Júpiter pueden acelerar las cargas hasta velocidades que superen la velocidad térmica, este tipo de distribuciones son muy útiles. En este caso las colas de las distribuciones de electrones siguen una ley de potencias. Se han estudiado las fluctuaciones de granos cargados para funciones de distribución kappa. iii) La tesis concluye con el análisis para deorbitar satélites con amarras electrodinámicas que siguen una Órbita Baja Terrestre (LEO). Una amarra debe presentar una baja probabilidad de corte por pequeño debris y además debe ser suficientemente ligero para que el cociente entre la masa de la amarra y el satélite sea muy pequeño. En este trabajo se estiman las medidas de la longitud, anchura y espesor que debe tener una amarra para minimizar el producto de la probabilidad de corte por el cociente entre las masas de la amarra y el satélite. Se presentan resultados preliminares del diseño de una amarra con forma de cinta para deorbitar satélites relativamente ligeros como Cryosat y pesados como Envisat. Las misiones espaciales a planetas exteriores y en el ámbito terrestre plantean importantes retos científico-tecnológicos que deben ser abordados y solucionados. Por ello, desde el inicio de la era espacial se han diseñando novedosos métodos propulsivos, sistemas de guiado, navegación y control más robustos, y nuevos materiales para mejorar el rendimiento de los vehículos espaciales (SC). En un gran número de misiones interplanetarias y en todas las misiones a planetas exteriores se han empleado sistemas de radioisótopos (RPS) para generar potencia eléctrica en los vehículos espaciales y en los rovers de exploración. Estos sistemas emplean como fuente de energía el escaso y costoso plutonio-238. La NASA, por medio de un informe de la National Academy of Science (5 de Mayo del 2009), expresó una profunda preocupación por la baja cantidad de plutonio almacenado, insuficiente para desarrollar todas las misiones de exploración planetaria planeadas en el futuro [81, 91]. Esta circustancia ha llevado a dicha Agencia tomar la decisión de limitar el uso de estos sistemas RPS en algunas misiones de especial interés científico y una recomendación de alta prioridad para que el Congreso de los EEUU apruebe el reestablecimiento de la producción de plutonio-238, -son necesarios cerca de 5 kg de este material radiactivo al año-, para salvaguardar las misiones que requieran dichos sistemas de potencia a partir del año 2018. Por otro lado, la Agencia estadounidense ha estado considerando el uso de fuentes de energía alternativa; como la fisión nuclear a través del ambicioso proyecto Prometheus, para llevar a cabo una misión de exploración en el sistema jovial (JIMO). Finalmente, dicha misión fue desestimada por su elevado coste. Recientemente se han estado desarrollando sistemas que consigan energía a través de los recursos naturales que nos aporta el Sol, mediante paneles solares -poco eficientes para misiones a planetas alejados de la luz solar-. En este contexto, la misión JUNO del programa Nuevas Fronteras de la NASA, cuyo lanzamiento fue realizado con éxito en Agosto de 2011, va a ser la primera misión equipada con paneles solares que sobrevolará Júpiter en el 2015 siguiendo una órbita polar. Anteriormente se habían empleado los antes mencionados RPS para las misiones Pioneer 10,11, Voyager 1,2, Ulysses, Cassini-Huygens y Galileo (todas sobrevuelos excepto Galileo). Dicha misión seguirá una órbita elíptica de alta excentricidad con un periápside muy cercano a Júpiter, y apoápside lejano, evitando que los intensos cinturones de radiación puedan dañar los instrumentos de navegación y científicos. Un tether o amarra electrodinámica es capaz de operar como sistema propulsivo o generador de potencia, pero también puede ser considerado como solución científicotecnológica en misiones espaciales tanto en LEO (Órbita Baja Terrestre) como en planetas exteriores. Siguiendo una perspectiva histórica, durante las misiones terrestres TSS-1 (1992) y TSS-1R (1996) se emplearon amarras estandard con recubrimiento aislante en toda su longitud, aplicando como terminal anódico pasivo un colector esférico para captar electrones. En una geometría alternativa, propuesta por J. R. Sanmartín et al. (1993) [93], se consideró dejar la amarra sin recubrimiento aislante (“bare tether”), y sin colector anódico esférico, de forma que recogiera electrones a lo largo del segmento que resulta polarizado positivo, como si se tratara de una sonda de Langmuir de gran longitud. A diferencia de la amarra estandard, el “bare tether” es capaz de recoger electrones a lo largo de una superficie grande ya que este segmento es de varios kilómetros de longitud. Como el radio de la amarra es del orden de la longitud de Debye y pequeño comparado con el radio de Larmor de los electrones, permite una recolección eficiente de electrones en el régimen OML (Orbital Motion Limited) de sondas de Langmuir. La corriente dada por la teoría OML varía en función del perímetro y la longitud. En el caso de una cinta delgada, el perímetro depende de la anchura, que debe ser suficientemente grande para evitar cortes producidos por debris y micrometeoritos, y suficientemente pequeño para que la amarra funcione en dicho régimen [95]. En el experimento espacial TSS-1R mencionado anteriormente, se identificó una recolección de corriente más elevada que la que predecía el modelo teórico de Parker- Murphy, debido posiblemente a que se utilizaba un colector esférico de radio bastante mayor que la longitud de Debye [79]. En el caso de una amarra “bare”, que recoge electrones a lo largo de gran parte de su longitud, se puede producir un fenómeno conocido como atrapamiento adiabático de electrones (adiabatic electron trapping) [25, 40, 60, 73, 74, 97]. En el caso terrestre (LEO) se da la condición mesotérmica en la que la amarra se mueve con una velocidad muy superior a la velocidad térmica de los iones del ambiente y muy inferior a la velocidad térmica de los electrones. J. Laframboise y L. Parker [57] mostraron que, para una función de distribución quasi-isotrópica, la densidad de electrones debe entonces ser necesariamente inferior a la densidad ambiente. Por otra parte, debido a su flujo hipersónico y a la alta polarización positiva de la amarra, la densidad de los iones es mayor que la densidad ambiente en una vasta región de la parte “ram” del flujo, violando la condición de cuasi-neutralidad,-en una región de dimensión mayor que la longitud de Debye-. La solución a esta paradoja podría basarse en el atrapamiento adiabático de electrones ambiente en órbitas acotadas entorno al tether. ABSTRACT New technological and scientific applications by electrodynamic tethers for planetary missions are analyzed: i) A set of cylindrical, parallel tethers (electric solar sail or e-sail) is considered for an interplanetary mission; ions from the solar wind are repelled by the high potential of the tether, providing momentum to the e-sail. An approximated model of a stationary potential for a high solar wind flow is considered. With the force provided by a negative biased tether, an indirect method for the optimization trajectory of an Earth-to-Jupiter orbit transfer is analyzed. ii) The deployment of a tether from the e-sail allows several scientific applications in Jupiter. iia) It might be used as a source of radiative waves for plasma diagnostics and artificial aurora generator. A conductive tether orbiting in the Jovian magnetosphere produces waves. Wave radiation by a conductor carrying a steady current in both a polar, highly eccentric, low perijove orbit, as in NASA’s Juno mission, and an equatorial low Jovian orbit (LJO) mission below the intense radiation belts, is considered. Both missions will need electric power generation for scientific instruments and communication systems. Tethers generate power more efficiently than solar panels or radioisotope power systems (RPS). The radiation impedance is required to determine the current in the overall tether circuit. In a cold plasma model, radiation occurs mainly in the Alfven and fast magnetosonic modes, exhibiting a large refraction index. The radiation impedance of insulated tethers is determined for both modes and either mission. Unlike the Earth ionospheric case, the low-density, highly magnetized Jovian plasma makes the electron gyrofrequency much larger than the plasma frequency; this substantially modifies the power spectrum for either mode by increasing the Alfven velocity. An estimation of the radiation impedance of bare tethers is also considered. iib) In LJO, a spacecraft orbiting in a slow downward spiral under the radiation belts would allow determining magnetic field structure and atmospheric composition for understanding the formation, evolution, and structure of Jupiter. Additionally, if the cathodic contactor is switched off, a tether floats electrically, allowing e-beam emission that generate auroras. On/off switching produces bias/current pulses and signal emission, which might be used for Jovian plasma diagnostics. In LJO, the ions impacting against the negative-biased tether do produce secondary electrons, which racing down Jupiter’s magnetic field lines, reach the upper atmosphere. The energetic electrons there generate auroral effects. Regions where the tether efficiently should produce secondary electrons are analyzed. iic) Other scientific application suggested in LJO is the in-situ detection of charged grains. Charged grains naturally orbit near Jupiter. High-energy electrons in the Jovian ambient may be modeled by the kappa distribution function. In complex plasma scenarios, where the Jovian high electric field may accelerate charges up superthermal velocities, the use of non-Maxwellian distributions should be considered. In these cases, the distribution tails fit well to a power-law dependence for electrons. Fluctuations of the charged grains for non-Mawellian distribution function are here studied. iii) The present thesis is concluded with the analysis for de-orbiting satellites at end of mission by electrodynamic tethers. A de-orbit tether system must present very small tether-to-satellite mass ratio and small probability of a tether cut by small debris too. The present work shows how to select tape dimensions so as to minimize the product of those two magnitudes. Preliminary results of tape-tether design are here discussed to minimize that function. Results for de-orbiting Cryosat and Envisat are also presented.

Tether radiation in Juno-type and circular-equatorial Jovian orbits

Wave radiation by a conductor carrying a steady current in both a polar, highly eccentric, low perijove orbit, as in NASA's planned Juno mission, and an equatorial low Jovian orbit (LJO) mission below the intense radiation belts, is considered. Both missions will need electric power generation for scientific instruments and communication systems. Tethers generate power more efficiently than solar panels or radioisotope power systems (RPS). The radiation impedance is required to determine the current in the overall tether circuit. In a cold plasma model, radiation occurs mainly in the Alfven and fast magnetosonic modes, exhibiting a large refraction index. The radiation impedance of insulated tethers is determined for both modes and either mission. Unlike the Earth ionospheric case, the low-density, highly magnetized Jovian plasma makes the electron gyrofrequency much larger than the plasma frequency; this substantially modifies the power spectrum for either mode by increasing the Alfven velocity. Finally, an estimation of the radiation impedance of bare tethers is considered. In LJO, a spacecraft orbiting in a slow downward spiral under the radiation belts would allow determining magnetic field structure and atmospheric composition for understanding the formation, evolution, and structure of Jupiter. Additionally, if the cathodic contactor is switched off, a tether floats electrically, allowing e-beam emission that generate auroras. On/off switching produces bias/current pulses and signal emission, which might be used for Jovian plasma diagnostics.

Chemical and structural analysis of Eucalyptus globulus and E. camaldulensis leaf cuticles: a lipidized cell wall region

The plant cuticle has traditionally been conceived as an independent hydrophobic layer that covers the external epidermal cell wall. Due to its complexity, the existing relationship between cuticle chemical composition and ultra-structure remains unclear to date. This study aimed to examine the link between chemical composition and structure of isolated, adaxial leaf cuticles of Eucalyptus camaldulensis and E. globulus by the gradual extraction and identification of lipid constituents (cutin and soluble lipids), coupled to spectroscopic and microscopic analyses. The soluble compounds and cutin monomers identified could not be assigned to a concrete internal cuticle ultra-structure. After cutin depolymerization, a cellulose network resembling the cell wall was observed, with different structural patterns in the regions ascribed to the cuticle proper and cuticular layer, respectively. Our results suggest that the current cuticle model should be revised, stressing the presence and major role of cell wall polysaccharides. It is concluded that the cuticle may be interpreted as a modified cell wall region which contains additional lipids. The major heterogeneity of the plant cuticle makes it difficult to establish a direct link between cuticle chemistry and structure with the existing methodologies.

Genetic characterization of chestnut (Castanea sativa Mill.) orchards and traditional nut varieties in El Bierzo, a glacial refuge and major cultivation site in Northwestern Spain

This paper presents a detailed genetic study of Castanea sativa in El Bierzo, a major nut production region with interesting features. It is located within a glacial refuge at one extreme of the distribution area (northwest Spain); it has a centenary tradition of chestnut management; and more importantly, it shows an unusual degree of genetic isolation. Seven nuclear microsatellite markers were selected to analyze the genetic variability and structure of 169 local trees grafted for nut production. We analyzed in the same manner 62 local nuts. The selected loci were highly discriminant for the genotypes studied, giving a combined probability of identity of 6.1 × 10−6. An unprecedented density of trees was sampled for this project over the entire region, and nuts were collected representing 18 cultivars marketed by local producers. Several instances of misclassification by local growers were detected. Fixation index estimates and analysis of molecular variance (AMOVA) data are supportive of an unexpectedly high level of genetic differentiation in El Bierzo, larger than that estimated in a previous study with broader geographical scope but based on limited local sampling (Pereira-Lorenzo et al., Tree Genet Genomes 6: 701–715, 2010a). Likewise, we have determined that clonality due to grafting had been previously overestimated. In line with these observations, no significant spatial structure was found using both a model-based Bayesian procedure and Mantel’s tests. Taken together, our results evidence the need for more fine-scale genetic studies if conservation strategies are to be efficiently improved.

Cuticle Structure in Relation to Chemical Composition: Re-assessing the Prevailing Model

The surface of most aerial plant organs is covered with a cuticle that provides protection against multiple stress factors including dehydration. Interest on the nature of this external layer dates back to the beginning of the 19th century and since then, several studies facilitated a better understanding of cuticular chemical composition and structure. The prevailing undertanding of the cuticle as a lipidic, hydrophobic layer which is independent from the epidermal cell wall underneath stems from the concept developed by Brongniart and von Mohl during the first half of the 19th century. Such early investigations on plant cuticles attempted to link chemical composition and structure with the existing technologies, and have not been directly challenged for decades. Beginning with a historical overview about the development of cuticular studies, this review is aimed at critically assessing the information available on cuticle chemical composition and structure, considering studies performed with cuticles and isolated cuticular chemical components. The concept of the cuticle as a lipid layer independent from the cell wall is subsequently challenged, based on the existing literature, and on new findings pointing toward the cell wall nature of this layer, also providing examples of different leaf cuticle structures. Finally, the need for a re-assessment of the chemical and structural nature of the plant cuticle is highlighted, considering its cell wall nature and variability among organs, species, developmental stages, and biotic and abiotic factors during plant growth.

High-precision 40Ar/39Ar geochronology and the advent of North America’s Late Cretaceous terrestrial fauna

A densely sampled, diverse new fauna from the uppermost Cedar Mountain Formation, Utah, indicates that the basic pattern of faunal composition for the Late Cretaceous of North America was already established by the Albian-Cenomanian boundary. Multiple, concordant 40Ar/39Ar determinations from a volcanic ash associated with the fauna have an average age of 98.39 ± 0.07 million years. The fauna of the Cedar Mountain Formation records the first global appearance of hadrosaurid dinosaurs, advanced lizard (e.g., Helodermatidae), and mammal (e.g., Marsupialia) groups, and the first North American appearance of other taxa such as tyrannosaurids, pachycephalosaurs, and snakes. Although the origin of many groups is unclear, combined biostratigraphic and phylogenetic evidence suggests an Old World, specifically Asian, origin for some of the taxa, an hypothesis that is consistent with existing evidence from tectonics and marine invertebrates. Large-bodied herbivores are mainly represented by low-level browsers, ornithopod dinosaurs, whose radiations have been hypothesized to be related to the initial diversification of angiosperm plants. Diversity at the largest body sizes (>106 g) is low, in contrast to both preceding and succeeding faunas; sauropods, which underwent demise in the Northern hemisphere coincident with the radiation of angiosperms, apparently went temporarily unreplaced by other megaherbivores. Morphologic and taxonomic diversity among small, omnivorous mammals, multituberculates, is also low. A later apparent increase in diversity occurred during the Campanian, coincident with the appearance of major fruit types among angiosperms, suggesting the possibility of adaptive response to new resources.

Role of brain allopregnanolone in the plasticity of γ-aminobutyric acid type A receptor in rat brain during pregnancy and after delivery

The relation between changes in brain and plasma concentrations of neurosteroids and the function and structure of γ-aminobutyric acid type A (GABAA) receptors in the brain during pregnancy and after delivery was investigated in rats. In contrast with plasma, where all steroids increased in parallel, the kinetics of changes in the cerebrocortical concentrations of progesterone, allopregnanolone (AP), and allotetrahydrodeoxycorticosterone (THDOC) diverged during pregnancy. Progesterone was already maximally increased between days 10 and 15, whereas AP and allotetrahydrodeoxycorticosterone peaked around day 19. The stimulatory effect of muscimol on 36Cl− uptake by cerebrocortical membrane vesicles was decreased on days 15 and 19 of pregnancy and increased 2 days after delivery. Moreover, the expression in cerebral cortex and hippocampus of the mRNA encoding for γ2L GABAA receptor subunit decreased during pregnancy and had returned to control values 2 days after delivery. Also α1,α2, α3, α4, β1, β2, β3, and γ2S mRNAs were measured and failed to change during pregnancy. Subchronic administration of finasteride, a 5α-reductase inhibitor, to pregnant rats reduced the concentrations of AP more in brain than in plasma as well as prevented the decreases in both the stimulatory effect of muscimol on 36Cl− uptake and the decrease of γ2L mRNA observed during pregnancy. These results indicate that the plasticity of GABAA receptors during pregnancy and after delivery is functionally related to fluctuations in endogenous brain concentrations of AP whose rate of synthesis/metabolism appears to differ in the brain, compared with plasma, in pregnant rats.

Patterned library analysis: A method for the quantitative assessment of hypotheses concerning the determinants of protein structure

Site-directed mutagenesis and combinatorial libraries are powerful tools for providing information about the relationship between protein sequence and structure. Here we report two extensions that expand the utility of combinatorial mutagenesis for the quantitative assessment of hypotheses about the determinants of protein structure. First, we show that resin-splitting technology, which allows the construction of arbitrarily complex libraries of degenerate oligonucleotides, can be used to construct more complex protein libraries for hypothesis testing than can be constructed from oligonucleotides limited to degenerate codons. Second, using eglin c as a model protein, we show that regression analysis of activity scores from library data can be used to assess the relative contributions to the specific activity of the amino acids that were varied in the library. The regression parameters derived from the analysis of a 455-member sample from a library wherein four solvent-exposed sites in an α-helix can contain any of nine different amino acids are highly correlated (P < 0.0001, R2 = 0.97) to the relative helix propensities for those amino acids, as estimated by a variety of biophysical and computational techniques.

Cryptocyanin, a crustacean molting protein: Evolutionary link with arthropod hemocyanins and insect hexamerins

Cryptocyanin, a copper-free hexameric protein in crab (Cancer magister) hemolymph, has been characterized and the amino acid sequence has been deduced from its cDNA. It is markedly similar in sequence, size, and structure to hemocyanin, the copper-containing oxygen-transport protein found in many arthropods. Cryptocyanin does not bind oxygen, however, and lacks three of the six highly conserved copper-binding histidine residues of hemocyanin. Cryptocyanin has no phenoloxidase activity, although a phenoloxidase is present in the hemolymph. The concentration of cryptocyanin in the hemolymph is closely coordinated with the molt cycle and reaches levels higher than hemocyanin during premolt. Cryptocyanin resembles insect hexamerins in the lack of copper, molt cycle patterns of biosynthesis, and potential contributions to the new exoskeleton. Phylogenetic analysis of sequence similarities between cryptocyanin and other members of the hemocyanin gene family shows that cryptocyanin is closely associated with crustacean hemocyanins and suggests that cryptocyanin arose as a result of a hemocyanin gene duplication. The presence of both hemocyanin and cryptocyanin in one animal provides an example of how insect hexamerins might have evolved from hemocyanin. Our results suggest that multiple members of the hemocyanin gene family—hemocyanin, cryptocyanin, phenoloxidase, and hexamerins—may participate in two vital functions of molting animals, oxygen binding and molting. Cryptocyanin may provide important molecular data to further investigate evolutionary relationships among all molting animals.

Efficient recognition of substrates and substrate analogs by the adenine glycosylase MutY requires the C-terminal domain

The Escherichia coli DNA repair enzyme MutY plays an important role in the prevention of DNA mutations by removing misincorporated adenine residues from 7,8-dihydro-8-oxo-2′-deoxyguanosine:2′-deoxyadenosine (OG:A) mispairs. The N-terminal domain of MutY (Stop 225, Met1–Lys225) has a sequence and structure that is characteristic of a superfamily of base excision repair glycosylases; however, MutY and its homologs contain a unique C-terminal domain. Previous studies have shown that the C-terminal domain confers specificity for OG:A substrates over G:A substrates and exhibits homology to the d(OG)TPase MutT, suggesting a role in OG recognition. In order to provide additional information on the importance of the C-terminal domain in damage recognition, we have investigated the kinetic properties of a form lacking this domain (Stop 225) under multiple- and single-turnover conditions. In addition, the interaction of Stop 225 with a series of non-cleavable substrate and product analogs was evaluated using gel retardation assays and footprinting experiments. Under multiple-turnover conditions Stop 225 exhibits biphasic kinetic behavior with both OG:A and G:A substrates, likely due to rate-limiting DNA product release. However, the rate of turnover of Stop 225 was increased 2-fold with OG:A substrates compared to the wild-type enzyme. In contrast, the intrinsic rate for adenine removal by Stop 225 from both G:A and OG:A substrates is significantly reduced (10- to 25-fold) compared to the wild-type. The affinity of Stop 225 for substrate analogs was dramatically reduced, as was the ability to discriminate between substrate analogs paired with OG over G. Interestingly, similar hydroxyl radical and DMS footprinting patterns are observed for Stop 225 and wild-type MutY bound to DNA duplexes containing OG opposite an abasic site mimic or a non-hydrogen bonding A analog, suggesting that similar regions of the DNA are contacted by both enzyme forms. Importantly, Stop 225 has a reduced ability to prevent DNA mutations in vivo. This implies that the reduced adenine glycosylase activity translates to a reduced capacity of Stop 225 to prevent DNA mutations in vivo.

PALI—a database of Phylogeny and ALIgnment of homologous protein structures

PALI (release 1.2) contains three-dimensional (3-D) structure-dependent sequence alignments as well as structure-based phylogenetic trees of homologous protein domains in various families. The data set of homologous protein structures has been derived by consulting the SCOP database (release 1.50) and the data set comprises 604 families of homologous proteins involving 2739 protein domain structures with each family made up of at least two members. Each member in a family has been structurally aligned with every other member in the same family (pairwise alignment) and all the members in the family are also aligned using simultaneous super­position (multiple alignment). The structural alignments are performed largely automatically, with manual interventions especially in the cases of distantly related proteins, using the program STAMP (version 4.2). Every family is also associated with two dendrograms, calculated using PHYLIP (version 3.5), one based on a structural dissimilarity metric defined for every pairwise alignment and the other based on similarity of topologically equivalent residues. These dendrograms enable easy comparison of sequence and structure-based relationships among the members in a family. Structure-based alignments with the details of structural and sequence similarities, superposed coordinate sets and dendrograms can be accessed conveniently using a web interface. The database can be queried for protein pairs with sequence or structural similarities falling within a specified range. Thus PALI forms a useful resource to help in analysing the relationship between sequence and structure variation at a given level of sequence similarity. PALI also contains over 653 ‘orphans’ (single member families). Using the web interface involving PSI_BLAST and PHYLIP it is possible to associate the sequence of a new protein with one of the families in PALI and generate a phylogenetic tree combining the query sequence and proteins of known 3-D structure. The database with the web interfaced search and dendrogram generation tools can be accessed at http://pa uling.mbu.iisc.ernet.in/~pali.

Potent, structurally constrained agonists and competitive antagonists of corticotropin-releasing factor.

Predictive methods, physicochemical measurements, and structure activity relationship studies suggest that corticotropin-releasing factor (CRF; corticoliberin), its family members, and competitive antagonists (resulting from N-terminal deletions) usually assume an alpha-helical conformation when interacting with the CRF receptor(s). To test this hypothesis further, we have scanned the whole sequence of the CRF antagonist [D-Phe12,Nle21,38]r/hCRF-(12-41) (r/hCRF, rat/human CRF; Nle, norleucine) with an i-(i + 3) bridge consisting of the Glu-Xaa-Xaa-Lys scaffold. We have found astressin [cyclo(30-33)[D-Phe12,Nle21,38,Glu30,Lys33]r/ hCRF(12-41)] to be approximately 30 times more potent than [D-Phe12,Nle21,38]r/hCRF-(12-41), our present standard, and 300 times more potent than the corresponding linear analog in an in vitro pituitary cell culture assay. Astressin has low affinity for the CRF binding protein and high affinity (Ki = 2 nM) for the cloned pituitary receptor. Radioiodinated [D-125I-Tyr12]astressin was found to be a reliable ligand for binding assays. In vivo, astressin is significantly more potent than any previously tested antagonist in reducing hypophyseal corticotropin (ACTH) secretion in stressed or adrenalectomized rats. The cyclo(30-33)[Ac-Pro4,D-Phe12,Nle21,38,Glu30,Lys33++ +]r/hCRF-(4-41) agonist and its linear analog are nearly equipotent, while the antagonist astressin and its linear form vary greatly in their potencies. This suggests that the lactam cyclization reinstates a structural constraint in the antagonists that is normally induced by the N terminus of the agonist.