Local Earthquake Tomography by trans-dimensional Monte Carlo Sampling

In this study a new, fully non-linear, approach to Local Earthquake Tomography is presented. Local Earthquakes Tomography (LET) is a non-linear inversion problem that allows the joint determination of earthquakes parameters and velocity structure from arrival times of waves generated by local sources. Since the early developments of seismic tomography several inversion methods have been developed to solve this problem in a linearized way. In the framework of Monte Carlo sampling, we developed a new code based on the Reversible Jump Markov Chain Monte Carlo sampling method (Rj-McMc). It is a trans-dimensional approach in which the number of unknowns, and thus the model parameterization, is treated as one of the unknowns. I show that our new code allows overcoming major limitations of linearized tomography, opening a new perspective in seismic imaging. Synthetic tests demonstrate that our algorithm is able to produce a robust and reliable tomography without the need to make subjective a-priori assumptions about starting models and parameterization. Moreover it provides a more accurate estimate of uncertainties about the model parameters. Therefore, it is very suitable for investigating the velocity structure in regions that lack of accurate a-priori information. Synthetic tests also reveal that the lack of any regularization constraints allows extracting more information from the observed data and that the velocity structure can be detected also in regions where the density of rays is low and standard linearized codes fails. I also present high-resolution Vp and Vp/Vs models in two widespread investigated regions: the Parkfield segment of the San Andreas Fault (California, USA) and the area around the Alto Tiberina fault (Umbria-Marche, Italy). In both the cases, the models obtained with our code show a substantial improvement in the data fit, if compared with the models obtained from the same data set with the linearized inversion codes.

Molecular events associated with trisomy of chromosome 7 in mouse skin chemical carcinogenesis

The primary objective of this study has been to investigate the effects at the molecular level of trisomy of mouse chromosome 7 in chemically induced skin tumors. It was previously proposed that the initiation event in the mouse skin carcinogenesis model is a heterozygous mutation of the Ha-ras-1 gene, mapped to chromosome 7. Previous studies in this laboratory identified trisomy 7 as one of the primary nonrandom cytogenetic abnormalities found in the majority of severely dysplastic papillomas and squamous cell carcinomas induced in SENCAR mice by an initiation-promotion protocol. Therefore, the first hypothesis tested was that trisomy 7 occurs by specific duplication of the chromosome carrying a mutated Ha-ras-1 allele. Results of a quantitative analysis of normal/mutated allelic ratios of the Ha-ras-1 gene confirmed this hypothesis, showing that most of the tumors exhibited overrepresentation of the mutated allele in the form of 1/2, 0/3, and 0/2 (normal/mutated) ratios. In addition, histopathological analysis of the tumors showed an apparent association between the degree of malignancy and the dosage of the mutated Ha-ras-1 allele. To determine the mechanism for loss of the normal Ha-ras-1 allele, found in 30% of the tumors, a comparison of constitutional and tumor genotypes was performed at different informative loci of chromosome 7. By combining Southern blot and polymerase chain reaction fragment length polymorphism analyses of DNAs extracted from squamous cell carcinomas, complete loss of heterozygosity was detected in 15 of 20 tumors at the Hbb locus, and in 5 of 5 tumors at the int-2 locus, both distal to Ha-ras-1. In addition, polymerase chain reaction analysis of DNA extracted from papillomas indicated that loss of heterozygosity occurs in late-stage lesions exhibiting a high degree of dysplasia and areas of microinvasion, suggesting that this event may be associated to the acquisition of the malignant phenotype. Allelic dosage analysis of tumors that had become homozygous at Hbb but retained heterozygosis at Ha-ras-1, indicated that loss of heterozygosity on mouse chromosome 7 occurs by a mitotic recombination mechanism. Overall, these findings suggest the presence of a putative tumor suppressor locus on the 7F1-ter region of mouse chromosome 7. Thus, loss of function by homozygosis at this putative suppressor locus may complement activation of the Ha-ras-1 gene during tumor progression, and might be associated with the malignant conversion stage of mouse skin carcinogenesis. ^

Polytomies and bayesian phylogenetic inference

Bayesian phylogenetic analyses are now very popular in systematics and molecular evolution because they allow the use of much more realistic models than currently possible with maximum likelihood methods. There are, however, a growing number of examples in which large Bayesian posterior clade probabilities are associated with very short edge lengths and low values for non-Bayesian measures of support such as nonparametric bootstrapping. For the four-taxon case when the true tree is the star phylogeny, Bayesian analyses become increasingly unpredictable in their preference for one of the three possible resolved tree topologies as data set size increases. This leads to the prediction that hard (or near-hard) polytomies in nature will cause unpredictable behavior in Bayesian analyses, with arbitrary resolutions of the polytomy receiving very high posterior probabilities in some cases. We present a simple solution to this problem involving a reversible-jump Markov chain Monte Carlo (MCMC) algorithm that allows exploration of all of tree space, including unresolved tree topologies with one or more polytomies. The reversible-jump MCMC approach allows prior distributions to place some weight on less-resolved tree topologies, which eliminates misleadingly high posteriors associated with arbitrary resolutions of hard polytomies. Fortunately, assigning some prior probability to polytomous tree topologies does not appear to come with a significant cost in terms of the ability to assess the level of support for edges that do exist in the true tree. Methods are discussed for applying arbitrary prior distributions to tree topologies of varying resolution, and an empirical example showing evidence of polytomies is analyzed and discussed.

Polytomies and bayesian phylogenetic inference

Bayesian phylogenetic analyses are now very popular in systematics and molecular evolution because they allow the use of much more realistic models than currently possible with maximum likelihood methods. There are, however, a growing number of examples in which large Bayesian posterior clade probabilities are associated with very short edge lengths and low values for non-Bayesian measures of support such as nonparametric bootstrapping. For the four-taxon case when the true tree is the star phylogeny, Bayesian analyses become increasingly unpredictable in their preference for one of the three possible resolved tree topologies as data set size increases. This leads to the prediction that hard (or near-hard) polytomies in nature will cause unpredictable behavior in Bayesian analyses, with arbitrary resolutions of the polytomy receiving very high posterior probabilities in some cases. We present a simple solution to this problem involving a reversible-jump Markov chain Monte Carlo (MCMC) algorithm that allows exploration of all of tree space, including unresolved tree topologies with one or more polytomies. The reversible-jump MCMC approach allows prior distributions to place some weight on less-resolved tree topologies, which eliminates misleadingly high posteriors associated with arbitrary resolutions of hard polytomies. Fortunately, assigning some prior probability to polytomous tree topologies does not appear to come with a significant cost in terms of the ability to assess the level of support for edges that do exist in the true tree. Methods are discussed for applying arbitrary prior distributions to tree topologies of varying resolution, and an empirical example showing evidence of polytomies is analyzed and discussed.

Phosphorylated α(1→4)Glucans as Substrate for Potato Starch-Branching Enzyme I1

The possible involvement of potato (Solanum tuberosum L.) starch-branching enzyme I (PSBE-I) in the in vivo synthesis of phosphorylated amylopectin was investigated in in vitro experiments with isolated PSBE-I using 33P-labeled phosphorylated and 3H end-labeled nonphosphorylated α(1→4)glucans as the substrates. From these radiolabeled substrates PSBE-I was shown to catalyze the formation of dual-labeled (3H/33P) phosphorylated branched polysaccharides with an average degree of polymerization of 80 to 85. The relatively high molecular mass indicated that the product was the result of multiple chain-transfer reactions. The presence of α(1→6) branch points was documented by isoamylase treatment and anion-exchange chromatography. Although the initial steps of the in vivo mechanism responsible for phosphorylation of potato starch remains elusive, the present study demonstrates that the enzyme machinery available in potato has the ability to incorporate phosphorylated α(1→4)glucans into neutral polysaccharides in an interchain catalytic reaction. Potato mini tubers synthesized phosphorylated starch from exogenously supplied 33PO43− and [U-14C]Glc at rates 4 times higher than those previously obtained using tubers from fully grown potato plants. This system was more reproducible compared with soil-grown tubers and was therefore used for preparation of 33P-labeled phosphorylated α(1→4)glucan chains.

Methacryloxyethyl phosphate-grafted expanded polytetrafluoroethylene membranes for biomedical applications

Expanded polytetrafluoroethylene (ePTFE) membranes were modified by graft copolymerization with methacryloxyethyl phosphate (MOEP) in methanol and 2-butanone (methyl ethyl ketone (MEK)) at ambient temperature using gamma irradiation. The effect of dose rate (0.46 and 4.6 kGyh(-1)), monomer concentration (1-40 %) and solvent were studied and the modified membranes were characterized by weight increase, X-ray photoelectron spectroscopy (XPS), Fourier transform infrared spectroscopy (FTIR) and scanning electron microscopy (SEM). XPS was used to determine the % degree of surface coverage using the C-F (ePTFE membrane) and the C-C (MOEP graft copolymer) peaks. Grafting yield, as well as surface coverage, were found to increase with increasing monomer concentration and were significantly higher for samples grafted in MEK than in methanol solution. SEM images showed distinctly different surface morphologies for the membranes grafted in methanol (smooth) and MEK (globular), hence indicating phase separation of the homopolymer in MEK. We propose that in our system, the non-solvent properties of MEK for the homopolymer play a more important role than solvent chain transfer reactions in determining grafting outcomes. (c) 2005 Society of Chemical Industry.

The ring opening polymerization of ring strained cyclic ethers

The kinetics and mechanisms of the ring-opening polymerization of oxetane were studied using cationic and coordinated anionic catalysts. The cationic initiators used were BF30Et2!/ethanol, BF30Et2!/ethanediol and BF30Et2/propantriol. Kinetic determinations with the BF30Et2/diol system indicated that a 1: 1 BF3:0H ratio gave the maximum rate of polymerization and this ratio was employed to detenmne the overall rates of polymerization. An overall second-order dependence was obtained when the system involved ethanediol or propantriol as co-catalyst and a 3/2-order dependence with ethanol, in each case the monomer gave a first-order relationship. This suggested that two mechanisms accounted for the cationic polymerization. These mechanisms were investigated and further evidence for these was obtained from the study of the complex formation of BF30Et2 and the co-catalysts by 1H NMR. Molecular weight studies (using size-exclusion chromatography) indicated that the hydroxyl ion acted as a chain transfer reagent when the [OH] > [BF3]. A linear relationship was observed when the number average molecular weight was plotted against [oxetane] at constant [BF3:0H], and similarly a linear dependency was observed on the BF3:0H 1:1 adduct at constant oxetane concentration. Copolymerization of oxetane and THF was carried out using BF30Et2/ethanol system. The reactivity ratios were calculated as rOXT = 1.2 ± 0.30 and rTHF = 0.14 ± 0.03. These copolymers were random copolymers with no evidence of oligomer formation. The coordinated anionic catalyst, porphinato-aluminium chloride [(TPP)AICl], was used to produce a living polymerization of oxetane. An overall third-order kinetics was obtained, with a second-order with respect to the [(TPP)AICl] and a first-order with respect to the [oxetane] and a mechanism was postulated using these results. The stereochemistry of [(TPP)AlCl] catalyst was investigated using cyclohexene and cyclopentene oxide monomers, using extensive 1H NMR, 2-D COSY and decoupling NMR techniques it was concluded that [(TPP)AlCl] gave rise to stereoregular polymers.

Promoter effects in the polymerisation of a mixed hydrocarbon feed with silica-supported BF3

The activity of a silica-supported BF3–methanol solid acid catalyst in the cationic polymerisation of an industrial aromatic C9 feedstock has been investigated. Reuse has been achieved under continuous conditions. Titration of the catalyst acid sites with triethylphosphine oxide (TEPO) in conjunction with 31P MAS NMR shows the catalyst to have two types of acid sites. Further analysis with 2,6 di-tert-butyl-4-methylpyridine (DBMP) has revealed the majority of these acid sites to be Brønsted in nature. The role of α-methylstyrene in promoting resin polymerisation via chain transfer is proposed.

Bayesian nonparametric latent variable models

L’un des problèmes importants en apprentissage automatique est de déterminer la complexité du modèle à apprendre. Une trop grande complexité mène au surapprentissage, ce qui correspond à trouver des structures qui n’existent pas réellement dans les données, tandis qu’une trop faible complexité mène au sous-apprentissage, c’est-à-dire que l’expressivité du modèle est insuffisante pour capturer l’ensemble des structures présentes dans les données. Pour certains modèles probabilistes, la complexité du modèle se traduit par l’introduction d’une ou plusieurs variables cachées dont le rôle est d’expliquer le processus génératif des données. Il existe diverses approches permettant d’identifier le nombre approprié de variables cachées d’un modèle. Cette thèse s’intéresse aux méthodes Bayésiennes nonparamétriques permettant de déterminer le nombre de variables cachées à utiliser ainsi que leur dimensionnalité. La popularisation des statistiques Bayésiennes nonparamétriques au sein de la communauté de l’apprentissage automatique est assez récente. Leur principal attrait vient du fait qu’elles offrent des modèles hautement flexibles et dont la complexité s’ajuste proportionnellement à la quantité de données disponibles. Au cours des dernières années, la recherche sur les méthodes d’apprentissage Bayésiennes nonparamétriques a porté sur trois aspects principaux : la construction de nouveaux modèles, le développement d’algorithmes d’inférence et les applications. Cette thèse présente nos contributions à ces trois sujets de recherches dans le contexte d’apprentissage de modèles à variables cachées. Dans un premier temps, nous introduisons le Pitman-Yor process mixture of Gaussians, un modèle permettant l’apprentissage de mélanges infinis de Gaussiennes. Nous présentons aussi un algorithme d’inférence permettant de découvrir les composantes cachées du modèle que nous évaluons sur deux applications concrètes de robotique. Nos résultats démontrent que l’approche proposée surpasse en performance et en flexibilité les approches classiques d’apprentissage. Dans un deuxième temps, nous proposons l’extended cascading Indian buffet process, un modèle servant de distribution de probabilité a priori sur l’espace des graphes dirigés acycliques. Dans le contexte de réseaux Bayésien, ce prior permet d’identifier à la fois la présence de variables cachées et la structure du réseau parmi celles-ci. Un algorithme d’inférence Monte Carlo par chaîne de Markov est utilisé pour l’évaluation sur des problèmes d’identification de structures et d’estimation de densités. Dans un dernier temps, nous proposons le Indian chefs process, un modèle plus général que l’extended cascading Indian buffet process servant à l’apprentissage de graphes et d’ordres. L’avantage du nouveau modèle est qu’il admet les connections entres les variables observables et qu’il prend en compte l’ordre des variables. Nous présentons un algorithme d’inférence Monte Carlo par chaîne de Markov avec saut réversible permettant l’apprentissage conjoint de graphes et d’ordres. L’évaluation est faite sur des problèmes d’estimations de densité et de test d’indépendance. Ce modèle est le premier modèle Bayésien nonparamétrique permettant d’apprendre des réseaux Bayésiens disposant d’une structure complètement arbitraire.

p-cresol as a reversible acylium ion scavenger in solid-phase peptide synthesis

In this work we have defined the nature of the p-cresol and p-thiocresol adducts generated from acylium ions during HF cleavage, following contemporary Boc/benzyl solid-phase peptide synthesis. Contrary to the results in previous reports, we found that both p-cresol and p-thiocresol predominantly form. aryl esters under typical cleavage conditions. Initially we investigated a number of small peptides containing either a single glutamate residue or a C-terminal long-chain amino acid which allowed us to unambiguously characterize the scavenged side products. Whereas, the p-cresol esters are stable at 0 degrees C they rearrange irreversibly at higher temperatures (5-20 degrees C) to form aryl ketones. By contrast, p-thiocresol esters do not undergo a Fries rearrangement but readily undergo further additions of p-thiocresol to form ketenebisthioacetals and trithio ortho esters, even at low temperatures. Importantly, we found by LC/MS and FT-ICR MS analysis that peptides containing p-cresol esters at glutamyl side chains are susceptible to amidation and fragmentation reactions at these sites during standard mild base workup procedures. The significance of these side reactions was further demonstrated in the synthesis of neutrophil immobilization factor, a 26-residue peptide, containing four glutamic acid residues. The side reactions were largely avoided by mild hydrogen peroxide-catalyzed hydrolysis which converted the p-cresol adducts to the free carboxylic acids in near quantitative yield. The choice of p-cresol as a reversible acylium ion scavenger when coupled with the simple workup conditions described is broadly applicable to Boc/benzyl peptide synthesis and will significantly enhance the quality of peptides produced.

Marker concentration patterns of labelled leaf and stem particles in the rumen of cattle grazing bermuda grass (Cynodon dactylon) analysed by reference to a raft model

Large (>1600 mum), ingestively masticated particles of bermuda grass (Cynodon dactylon L. Pers.) leaf and stem labelled with Yb-169 and Ce-144 respectively were inserted into the rumen digesta raft of heifers grazing bermuda grass. The concentration of markers in digesta sampled from the raft and ventral rumen were monitored at regular intervals over approximately 144 h. The data from the two sampling sites were simultaneously fitted to two pool (raft and ventral rumen-reticulum) models with either reversible or sequential flow between the two pools. The sequential flow model fitted the data equally as well as the reversible flow model but the reversible flow model was used because of its greater application. The reversible flow model, hereafter called the raft model, had the following features: a relatively slow age-dependent transfer rate from the raft (means for a gamma 2 distributed rate parameter for leaf 0.0740 v. stem 0.0478 h(-1)), a very slow first order reversible flow from the ventral rumen to the raft (mean for leaf and stem 0.010 h(-1)) and a very rapid first order exit from the ventral rumen (mean of leaf and stem 0.44 h(-1)). The raft was calculated to occupy approximately 0.82 total rumen DM of the raft and ventral rumen pools. Fitting a sequential two pool model or a single exponential model individually to values from each of the two sampling sites yielded similar parameter values for both sites and faster rate parameters for leaf as compared with stem, in agreement with the raft model. These results were interpreted as indicating that the raft forms a large relatively inert pool within the rumen. Particles generated within the raft have difficulty escaping but once into the ventral rumen pool they escape quickly with a low probability of return to the raft. It was concluded that the raft model gave a good interpretation of the data and emphasized escape from and movement within the raft as important components of the residence time of leaf and stem particles within the rumen digesta of cattle.

Spectroscopic Studies and Characterization of a Novel Electron-Transfer Chain

A novel two-component enzyme system from Escherichia coli involving a flavorubredoxin (FlRd) and its reductase was studied in terms of spectroscopic, redox, and biochemical properties of its constituents. FlRd contains one FMN and one rubredoxin (Rd) center per monomer. To assess the role of the Rd domain, FlRd and a truncated form lacking the Rd domain (FlRd¢Rd), were characterized. FlRd contains 2.9 ( 0.5 iron atoms/subunit, whereas FlRd¢Rd contains 2.1 ( 0.6 iron atoms/subunit. While for FlRd one iron atom corresponds to the Rd center, the other two irons, also present in FlRd¢Rd, are most probably due to a di-iron site. Redox titrations of FlRd using EPR and visible spectroscopies allowed us to determine that the Rd site has a reduction potential of -140 ( 15 mV, whereas the FMN undergoes reduction via a red-semiquinone, at -140 ( 15 mV (Flox/Flsq) and -180 ( 15 mV (Flsq/Flred), at pH 7.6. The Rd site has the lowest potential ever reported for a Rd center, which may be correlated with specific amino acid substitutions close to both cysteine clusters. The gene adjacent to that encoding FlRd was found to code for an FAD-containing protein, (flavo)rubredoxin reductase (FlRd-reductase), which is capable of mediating electron transfer from NADH to DesulfoVibrio gigas Rd as well as to E. coli FlRd. Furthermore, electron donation was found to proceed through the Rd domain of FlRd as the Rd-truncated protein does not react with FlRd-reductase. In vitro, this pathway links NADH oxidation with dioxygen reduction. The possible function of this chain is discussed considering the presence of FlRd homologues in all known genomes of anaerobes and facultative aerobes.

The role of halouracils in radiotherapy studied by electron transfer in atom-molecule collisions experiments

Dissertação para obtenção do Grau de Doutor em Engenharia Física

Analysis of the activation mechanism of Pseudomonas stutzeri cytochrome c peroxidase through an electron transfer chain

J Biol Inorg Chem (2011) 16:881–888 DOI 10.1007/s00775-011-0785-8

CD8beta endows CD8 with efficient coreceptor function by coupling T cell receptor/CD3 to raft-associated CD8/p56(lck) complexes.

The extraordinary sensitivity of CD8+ T cells to recognize antigen impinges to a large extent on the coreceptor CD8. While several studies have shown that the CD8beta chain endows CD8 with efficient coreceptor function, the molecular basis for this is enigmatic. Here we report that cell-associated CD8alphabeta, but not CD8alphaalpha or soluble CD8alphabeta, substantially increases the avidity of T cell receptor (TCR)-ligand binding. To elucidate how the cytoplasmic and transmembrane portions of CD8beta endow CD8 with efficient coreceptor function, we examined T1.4 T cell hybridomas transfected with various CD8beta constructs. T1.4 hybridomas recognize a photoreactive Plasmodium berghei circumsporozoite (PbCS) peptide derivative (PbCS (4-azidobezoic acid [ABA])) in the context of H-2K(d), and permit assessment of TCR-ligand binding by TCR photoaffinity labeling. We find that the cytoplasmic portion of CD8beta, mainly due to its palmitoylation, mediates partitioning of CD8 in lipid rafts, where it efficiently associates with p56(lck). In addition, the cytoplasmic portion of CD8beta mediates constitutive association of CD8 with TCR/CD3. The resulting TCR-CD8 adducts exhibit high affinity for major histocompatibility complex (MHC)-peptide. Importantly, because CD8alphabeta partitions in rafts, its interaction with TCR/CD3 promotes raft association of TCR/CD3. Engagement of these TCR/CD3-CD8/lck adducts by multimeric MHC-peptide induces activation of p56(lck) in rafts, which in turn phosphorylates CD3 and initiates T cell activation.