Implementation of Quantitative Reaction-based Model Refinement

Building a computational model for complex biological systems is an iterative process. It starts from an abstraction of the process and then incorporates more details regarding the specific biochemical reactions which results in the change of the model fit. Meanwhile, the model’s numerical properties such as its numerical fit and validation should be preserved. However, refitting the model after each refinement iteration is computationally expensive resource-wise. There is an alternative approach which ensures the model fit preservation without the need to refit the model after each refinement iteration. And this approach is known as quantitative model refinement. The aim of this thesis is to develop and implement a tool called ModelRef which does the quantitative model refinement automatically. It is both implemented as a stand-alone Java application and as one of Anduril framework components. ModelRef performs data refinement of a model and generates the results in two different well known formats (SBML and CPS formats). The development of this tool successfully reduces the time and resource needed and the errors generated as well by traditional reiteration of the whole model to perform the fitting procedure.

Computational models for and from biology : simple gene assembly and reaction systems

The advancement of science and technology makes it clear that no single perspective is any longer sufficient to describe the true nature of any phenomenon. That is why the interdisciplinary research is gaining more attention overtime. An excellent example of this type of research is natural computing which stands on the borderline between biology and computer science. The contribution of research done in natural computing is twofold: on one hand, it sheds light into how nature works and how it processes information and, on the other hand, it provides some guidelines on how to design bio-inspired technologies. The first direction in this thesis focuses on a nature-inspired process called gene assembly in ciliates. The second one studies reaction systems, as a modeling framework with its rationale built upon the biochemical interactions happening within a cell. The process of gene assembly in ciliates has attracted a lot of attention as a research topic in the past 15 years. Two main modelling frameworks have been initially proposed in the end of 1990s to capture ciliates’ gene assembly process, namely the intermolecular model and the intramolecular model. They were followed by other model proposals such as templatebased assembly and DNA rearrangement pathways recombination models. In this thesis we are interested in a variation of the intramolecular model called simple gene assembly model, which focuses on the simplest possible folds in the assembly process. We propose a new framework called directed overlap-inclusion (DOI) graphs to overcome the limitations that previously introduced models faced in capturing all the combinatorial details of the simple gene assembly process. We investigate a number of combinatorial properties of these graphs, including a necessary property in terms of forbidden induced subgraphs. We also introduce DOI graph-based rewriting rules that capture all the operations of the simple gene assembly model and prove that they are equivalent to the string-based formalization of the model. Reaction systems (RS) is another nature-inspired modeling framework that is studied in this thesis. Reaction systems’ rationale is based upon two main regulation mechanisms, facilitation and inhibition, which control the interactions between biochemical reactions. Reaction systems is a complementary modeling framework to traditional quantitative frameworks, focusing on explicit cause-effect relationships between reactions. The explicit formulation of facilitation and inhibition mechanisms behind reactions, as well as the focus on interactions between reactions (rather than dynamics of concentrations) makes their applicability potentially wide and useful beyond biological case studies. In this thesis, we construct a reaction system model corresponding to the heat shock response mechanism based on a novel concept of dominance graph that captures the competition on resources in the ODE model. We also introduce for RS various concepts inspired by biology, e.g., mass conservation, steady state, periodicity, etc., to do model checking of the reaction systems based models. We prove that the complexity of the decision problems related to these properties varies from P to NP- and coNP-complete to PSPACE-complete. We further focus on the mass conservation relation in an RS and introduce the conservation dependency graph to capture the relation between the species and also propose an algorithm to list the conserved sets of a given reaction system.

Prospects of Bitcoin diffusion

Currency is something people deal with every day in their lives. The contemporary society is very much revolving around currencies. Even though technological development has been rapid, the principle of currency has stayed relatively unchanged for a long time. Bitcoin is a digital currency that introduced an alternative to other digital currencies, and to the traditional physical currencies. Bitcoin is peer-to-peer, open source, and it erases the need of a third party in transactions. Bitcoin has since inception gained certain fame, but it has not established itself as a common currency in the world. The purpose of this study was to analyse what kind of potential does Bitcoin have to become a widely accepted currency in day-to-day transactions. The main research question was divided into three sub questions: • What kind of a process is the diffusion of new innovations? • What kinds of factors speak for the wider adoption of Bitcoin? • What kinds of factors speak against the wider adoption of Bitcoin? The purpose of the study was approached by having diffusion of innovations as the theoretical framework. The four elements in diffusion of innovations are, innovation, communication, time, and social system. The theoretical framework is applied to Bitcoin, and the research questions answered by analysing Bitcoin’s potential diffusion prospects. The body of research data consisted of media texts and statistics. In this study, content analysis was the research method. The main findings of the study are that Bitcoin has clear strengths, but it faces a large amount of uncertainty. Bitcoin’s strong areas are the transactions. They are fast, easy, and cheap. From the innovation diffusion perspective Bitcoin is still relatively unknown, and the general public’s attitudes towards it are sceptical. The research findings purport that Bitcoin has potential demand especially when the financial system of a region is dysfunctional, or when there is a financial crisis. Bitcoin is not very trusted, and the majority of people do not see a reason to start using Bitcoin in the future. A large number of people associate it with illegal activities. In general people are largely unaware of what Bitcoin is or what are the strengths and weaknesses. Bitcoin is an innovative alternative currency. However, unless people see a major need for Bitcoin due to a financial crisis, or dysfunctionality in the financial system, Bitcoin will not become much more widespread as it is today. Bitcoin’s underlying technology can be harnessed to multiple uses. Developments in that field in the future are something that future researchers could look into.

Two-dimensional drift-diffusion simulation of organic solar cells

The aim of this master's thesis is to develop a two-dimensional drift-di usion model, which describes charge transport in organic solar cells. The main bene t of a two-dimensional model compared to a one-dimensional one is the inclusion of the nanoscale morphology of the active layer of a bulk heterojunction solar cell. The developed model was used to study recombination dynamics at the donor-acceptor interface. In some cases, it was possible to determine e ective parameters, which reproduce the results of the two-dimensional model in the one-dimensional case. A summary of the theory of charge transport in semiconductors was presented and discussed in the context of organic materials. Additionally, the normalization and discretization procedures required to nd a numerical solution to the charge transport problem were outlined. The charge transport problem was solved by implementing an iterative scheme called successive over-relaxation. The obtained solution is given as position-dependent electric potential, free charge carrier concentrations and current densities in the active layer. An interfacial layer, separating the pure phases, was introduced in order to describe charge dynamics occurring at the interface between the donor and acceptor. For simplicity, an e ective generation of free charge carriers in the interfacial layer was implemented. The pure phases simply act as transport layers for the photogenerated charges. Langevin recombination was assumed in the two-dimensional model and an analysis of the apparent recombination rate in the one-dimensional case is presented. The recombination rate in a two-dimensional model is seen to e ectively look like reduced Langevin recombination at open circuit. Replicating the J-U curves obtained in the two-dimensional model is, however, not possible by introducing a constant reduction factor in the Langevin recombination rate. The impact of an acceptor domain in the pure donor phase was investigated. Two cases were considered, one where the acceptor domain is isolated and another where it is connected to the bulk of the acceptor. A comparison to the case where no isolated domains exist was done in order to quantify the observed reduction in the photocurrent. The results show that all charges generated at the isolated domain are lost to recombination, but the domain does not have a major impact on charge transport. Trap-assisted recombination at interfacial trap states was investigated, as well as the surface dipole caused by the trapped charges. A theoretical expression for the ideality factor n_id as a function of generation was derived and shown to agree with simulation data. When the theoretical expression was fitted to simulation data, no interface dipole was observed.

A reaction engineering approach to homogeneously catalyzed glycerol hydrochlorination

The production of biodiesel through transesterification has created a surplus of glycerol on the international market. In few years, glycerol has become an inexpensive and abundant raw material, subject to numerous plausible valorisation strategies. Glycerol hydrochlorination stands out as an economically attractive alternative to the production of biobased epichlorohydrin, an important raw material for the manufacturing of epoxy resins and plasticizers. Glycerol hydrochlorination using gaseous hydrogen chloride (HCl) was studied from a reaction engineering viewpoint. Firstly, a more general and rigorous kinetic model was derived based on a consistent reaction mechanism proposed in the literature. The model was validated with experimental data reported in the literature as well as with new data of our own. Semi-batch experiments were conducted in which the influence of the stirring speed, HCl partial pressure, catalyst concentration and temperature were thoroughly analysed and discussed. Acetic acid was used as a homogeneous catalyst for the experiments. For the first time, it was demonstrated that the liquid-phase volume undergoes a significant increase due to the accumulation of HCl in the liquid phase. Novel and relevant features concerning hydrochlorination kinetics, HCl solubility and mass transfer were investigated. An extended reaction mechanism was proposed and a new kinetic model was derived. The model was tested with the experimental data by means of regression analysis, in which kinetic and mass transfer parameters were successfully estimated. A dimensionless number, called Catalyst Modulus, was proposed as a tool for corroborating the kinetic model. Reactive flash distillation experiments were conducted to check the commonly accepted hypothesis that removal of water should enhance the glycerol hydrochlorination kinetics. The performance of the reactive flash distillation experiments were compared to the semi-batch data previously obtained. An unforeseen effect was observed once the water was let to be stripped out from the liquid phase, exposing a strong correlation between the HCl liquid uptake and the presence of water in the system. Water has revealed to play an important role also in the HCl dissociation: as water was removed, the dissociation of HCl was diminished, which had a retarding effect on the reaction kinetics. In order to obtain a further insight on the influence of water on the hydrochlorination reaction, extra semi-batch experiments were conducted in which initial amounts of water and the desired product were added. This study revealed the possibility to use the desired product as an ideal “solvent” for the glycerol hydrochlorination process. A co-current bubble column was used to investigate the glycerol hydrochlorination process under continuous operation. The influence of liquid flow rate, gas flow rate, temperature and catalyst concentration on the glycerol conversion and product distribution was studied. The fluid dynamics of the system showed a remarkable behaviour, which was carefully investigated and described. Highspeed camera images and residence time distribution experiments were conducted to collect relevant information about the flow conditions inside the tube. A model based on the axial dispersion concept was proposed and confronted with the experimental data. The kinetic and solubility parameters estimated from the semi-batch experiments were successfully used in the description of mass transfer and fluid dynamics of the bubble column reactor. In light of the results brought by the present work, the glycerol hydrochlorination reaction mechanism has been finally clarified. It has been demonstrated that the reactive distillation technology may cause drawbacks to the glycerol hydrochlorination reaction rate under certain conditions. Furthermore, continuous reactor technology showed a high selectivity towards monochlorohydrins, whilst semibatch technology was demonstrated to be more efficient towards the production of dichlorohydrins. Based on the novel and revealing discoveries brought by the present work, many insightful suggestions are made towards the improvement of the production of αγ-dichlorohydrin on an industrial scale.

Mass transfer and reaction kinetics in single droplet systems

As reactive extraction grown more and more popular in a variety of technological applications, optimizing its performance becomes more and more important. The process of complex formation is affected by a great number of both physical and chemical properties of all the components involved, and sometimes their interference with one another makes improving the effectiveness of such processes very difficult. In this Master’s Theses, the processes of complex formation between the aqueous phase - represented by copper sulfate water solution, and organic phase – represented by Acorga M5640 solvent extractor, were studied in order to establish the effect these components have on reactive extraction performance and to determine which step is bottlenecking the process the most.

Investigation of single tryptophan proteins encapsulated in TEOS-derived sol-gel matrices by fluorescence spectroscopy /

In the work reported here, optically clear, ultrathin TEOS derived sol-gel slides which were suitable for studies of tryptophan (Trp) fluorescence from entrapped proteins were prepared by the sol-gel technique and characterized. The monitoring of intrinsic protein fluorescence provided information about the structure and environment of the entrapped protein, and about the kinetics of the interaction between the entrapped protein and extemal reagents. Initial studies concentrated on the single Trp protein monellin which was entrapped into the sol-gel matrices. Two types of sol-gel slides, termed "wet aged", in which the gels were aged in buffer and "dry-aged", in which the gels were aged in air , were studied in order to compare the effect of the sol-gel matrix on the structure of the protein at different aging stages. Fluorescence results suggested that the mobility of solvent inside the slides was substantially reduced. The interaction of the entrapped protein with both neutral and charged species was examined and indicated response times on the order of minutes. In the case of the neutral species the kinetics were diffusion limited in solution, but were best described by a sum of first order rate constants when the reactions occurred in the glass matrix. For charged species, interactions between the analytes and the negatively charged glass matrix caused the reaction kinetics to become complex, with the overall reaction rate depending on both the type of aging and the charge on the analyte. The stability and conformational flexibility of the entrapped monellin were also studied. These studies indicated that the encapsulation of monellin into dry-aged monoliths caused the thermal unfolding transition to broaden and shift upward by 14°C, and causedthe long-term stability to improve by 12-fold (compared to solution). Chemical stability studies also showed a broader transition for the unfolding of the protein in dry-aged monoliths, and suggested that the protein was present in a distribution of environments. Results indicated that the entrapped proteins had a smaller range of conformational motions compared to proteins in solution, and that entrapped proteins were not able to unfold completely. The restriction of conformational motion, along with the increased structural order of the internal environment of the gels, likely resulted in the improvements in themial and long-term stability that were observed. A second protein which was also studied in this work is the metal binding protein rat oncomodulin. Initially, the unfolding behavior of this protein in aqueous solution was examined. Several single tryptophan mutants of the metal-binding protein rat oncomodulin (OM) were examined; F102W, Y57W, Y65W and the engineered protein CDOM33 which had all 12 residues of the CD loop replaced with a higher affinity binding loop. Both the thermal and the chemical stability were improved upon binding of metal ions with the order apo < Ca^^ < Tb^"^. During thermal denaturation, the transition midpoints (Tun) of Y65W appeared to be the lowest, followed by Y57W and F102W. The placement of the Trp residue in the F-helix in F102W apparently made the protein slightly more thermostable, although the fluorescence response was readily affected by chemical denaturants, which probably acted through the disruption of hydrogen bonds at the Cterminal end of the F-helix. Under both thermal and chemical denaturation, the engineered protein showed the highest stability. This indicated that increasing the number of metal ligating oxygens in the binding site, either by using a metal ion with a higher coordinatenumber (i.e. Tb^*) which binds more carboxylate ligands, or by providing more ligating groups, as in the CDOM33 replacement, produces notable improvements in protein stability. Y57W and CE)OM33 OM were chosen for further studies when encapsulated into sol-gel derived matrices. The kinetics of interaction of terbium with the entrapped proteins, the ability of the entrapped protein to binding terbium, as well as thermal stability of these two entrapped protein were compared with different levels of Ca^"*^ present in the matrix and in solution. Results suggested that for both of the proteins, the response time and the ability to bind terbium could be adjusted by adding excess calcium to the matrix before gelation. However, the less stable protein Y57W only retained at most 45% of its binding ability in solution while the more stable protein CDOM33 was able to retain 100% binding ability. Themially induced denaturation also suggested that CDOM33 showed similar stability to the protein in solution while Y57W was destabilized. All these results suggested that "hard" proteins (i.e. very stable) can easily survive the sol-gel encapsulation process, but "soft" proteins with lower thermodynamic stability may not be able to withstand the sol-gel process. However, it is possible to control many parameters in order to successfully entrap biological molecules into the sol-gel matrices with maxunum retention of activity.

Diffusion Monte Carlo study of electronic properties for H and Be atoms /

We examined three different algorithms used in diffusion Monte Carlo (DMC) to study their precisions and accuracies in predicting properties of isolated atoms, which are H atom ground state, Be atom ground state and H atom first excited state. All three algorithms — basic DMC, minimal stochastic reconfiguration DMC, and pure DMC, each with future-walking, are successfully impletmented in ground state energy and simple moments calculations with satisfactory results. Pure diffusion Monte Carlo with future-walking algorithm is proven to be the simplest approach with the least variance. Polarizabilities for Be atom ground state and H atom first excited state are not satisfactorily estimated in the infinitesimal differentiation approach. Likewise, an approach using the finite field approximation with an unperturbed wavefunction for the latter system also fails. However, accurate estimations for the a-polarizabilities are obtained by using wavefunctions that come from the time-independent perturbation theory. This suggests the flaw in our approach to polarizability estimation for these difficult cases rests with our having assumed the trial function is unaffected by infinitesimal perturbations in the Hamiltonian.

The application of the structural correlation method to determine the reaction coordinate for a putative ring closure reaction

The cr ystal structure of the compound 2-benzoylethylidene-3-(2,4- dibromophenyl)-2,3-dihydro-5-phenyl-l,3,4-thiadiazole* C23H16Br2NZOS (BRMEO) has been determined by using three dimensiona l x-ray diffraction data. The crys tal form is monoclinic, space group P21/c, a = 17.492(4), o -.t' 0 R 0 b =: 16.979(1), c = 14.962(1) A, "X. =o= 90 ',= 106.46(1) , z = 8, graphite-monochromatized Mo~ rad iation, Jl= 0.710J3~, D = 1.62g/cc and o D = 1.65g/cc. The data were col lected on ~ Nonius CAD-4 c diffractometer. The following atoms were made anisotropic: Br, S, N, 0, C7, and C14-C16 for each i ndependent molecu le ; the rest were left isotropic. For 3112 independent refl ec tions with F > 6G\F), R == 0.057. The compound has two independent molecules within the asymmetric unit. Two different conformers were observed which pack well together. /l The S---O interaction distances of 2.493(6) and 2 . 478(7) A were observed for molecules A and B respectively. These values are consistent with earlier findings for 2-benzoylmethylene-3-(2,4-dibromophenyl)- ~~ 2,3-dihydro-5-phenyl-l,3,4-thiadiazole C22H14Br2N20S (BRPHO) and 2-benzoylpropylidene-3-(2,4-dibromophenyl)-2,3-dihydroiii ,'r 5-phenyl-l,3,4-thiadiazole C24H18Br2N20S (BRPETO ) where S---O distances are l ess than the van der Waals (3.251\) but greater than those expected for () a single bond (1.50A). From the results and the literature it appears obvious that the energy/reaction coordinate pathway has a minimum between the end structures (the mono- and bicyclic compounds). * See reference (21) for nomenclature.

The diffusion of Co⠶溩 Zrâ â Tiâ â alloy /|nS. Kumar. -- 260 St Catharines [Ont.] : Dept. of Physics, Brock University,

The diffusion of Co60 in the body centered cubic beta phase of a ZrSOTi SO alloy has been studied at 900°, 1200°, and 1440°C. The results confirm earlier unpublished data obtained by Kidson17 • The temperature dependence of the diffusion coefficient is unusual and suggests that at least two and possibly three mechanisms may be operative Annealing of the specimen in the high B.C.C. region prior to the deposition of the tracer results in a large reduction in the diffusion coefficient. The possible significance of this effect is discussed in terms of rapid transport along dislocation network.

Exact property estimation from diffusion Monte Carlo with minimal stochastic reconfiguration /

Our objective is to develop a diffusion Monte Carlo (DMC) algorithm to estimate the exact expectation values, ($o|^|^o), of multiplicative operators, such as polarizabilities and high-order hyperpolarizabilities, for isolated atoms and molecules. The existing forward-walking pure diffusion Monte Carlo (FW-PDMC) algorithm which attempts this has a serious bias. On the other hand, the DMC algorithm with minimal stochastic reconfiguration provides unbiased estimates of the energies, but the expectation values ($o|^|^) are contaminated by ^, an user specified, approximate wave function, when A does not commute with the Hamiltonian. We modified the latter algorithm to obtain the exact expectation values for these operators, while at the same time eliminating the bias. To compare the efficiency of FW-PDMC and the modified DMC algorithms we calculated simple properties of the H atom, such as various functions of coordinates and polarizabilities. Using three non-exact wave functions, one of moderate quality and the others very crude, in each case the results are within statistical error of the exact values.

Modifications to the Suzuki reaction and mechanistic insights on the NBS mediated cleavage of benzylidene acetals /

One of the most challenging tasks for a synthetic organic chemist today, is the development of chemo, regio, and stereoselective methodologies toward the total synthesis of macromolecules. r . The objective of my thesis was to develop methodologies towards this end. The first part of my project was to develop highly functionalized chirons from D-glucose, a cheap, chiral starting material, to be utilized in this capacity. The second part of the project dealt with modifying the carbon-carbon bond forming Suzuki reaction, which is utilized quite often as a means of combining molecular sub units in total synthesis applications. As previously stated the first area of the project was to develop high value chirons from D-glucose, but the mechanism of their formation was also investigated. The free radical initiated oxidative fragmentation of benzylidene acetals was investigated through the use of several test-case substrates in order to unravel the possible mechanistic pathways. This was performed by reacting the different acetals with N-bromosuccinimide and benzoyl peroxide in chlorobenzene at 70^C in all cases. Of the three mechanistic pathways discussed in the literature, it was determined, from the various reaction products obtained, that the fragmentation of the initial benzylic radical does not occur spontaneously but rather, oxidation proceeds to give the benzyl bromide, which then fragments via a polar pathway. It was also discovered that the regioselectivity of the fragmentation step could be altered through incorporation of an allylic system into the benzylidene acetal. This allows for the acquisition of a new set of densely functionalized. chiral, valuable synthetic intermediates in only a few steps and in high yields from a-Dglucose. The second part of the project was the utilization of the phosphonium salt room temperature ionic liquid tetradecyltrihexylphosphonium chloride (THPC) as an efficient reusable medium for the palladium catalyzed Suzuki cross-coupling reaction of aryl halides, including aryl chlorides, under mild conditions. The cross-coupling reactions were found to proceed in THPC containing small amounts of water and toluene using potassium phosphate and 1% Pd2(dba)3. Variously substituted iodobenzenes, including electron rich derivatives, reacted efficiently in THPC with a variety of arylboronic acids and afforded complete conversion within 1 hour at 50 ^C. The corresponding aryl bromides also reacted under these conditions with the addition of a catalytic amount of triphenylphosphine that allowed for complete conversion and high isolated yields. The reactions involving aryl chlorides were considerably slower, although the addition of triphenylphosphine and heating at 70 ^C allowed high conversion of electron deficient derivatives. Addition of water and hexane to the reaction products results in a triphasic system in which the top hexane phase contained the biaryl products, the palladium catalyst remained fully dissolved in the central THPC layer, while the inorganic salts were extracted into the lower aqueous phase. The catalyst was then recycled by removing the top and bottom layers and adding the reagents to the ionic liquid which was heated again at 50 ^C; resulting in complete turnover of iodobenzene. Repetition of this procedure gave the biphenyl product in 82-97% yield (repeated five times) for both the initial and recycled reaction sequences.

A G. C./M. S. study of the reaction and decomposition products of pentafluorophenyl-grignard and lithium reagents

Decomposition and side reactions of, and the synthetic use of, pentafluorophenylmagnesium bromide and pentafluorophenyllithium have been investigated using G,C9/M.S, techniques• Their reactions with reagents such as CgF^X (X - H, F, CI, Br, 1), C6F4X2 (X - H, CI)f C6F3C13, C6H6. (CgX5)3P (X = H, F), (C6X5)3P=0 (X = H, F), (CgX5)Si (CH3)3 (X = H, F) and (CH0K SiCl , n = 1,2, in ether or ether/n-hexane were studied• In addition to the principal reaction of synthetic use, namely the replacement of a halogen by a pentafluorophenyl group, two types of side reactions were observed* These were (i) intermolecular loss of LiF via a nucleophilic substitution, and (ii) intramolecular loss of LiF, followed by the addition of either inorganic salts such as lithium or magnesium halides, or organometal compounds such as organolithium or organo-Grigaard* G.C«/M.S. techniques were routinely employed to study complicated reaction mixtures. Although mass spectrometry alone has disadvantages for the identification of isomers, deduction of the most probable pathway often helps overcome this problem.

Biochemical and histological investigations of viral localisation in the hypersensitive reaction of Phaselous vulgaris L. var Pinto to tobacco mosaic virus infection

STOBBS, Lorne,W ABSTRACT Biochemical and Histological Investigations of viral localisation in the hypersensitive reaction of Phaseolus vulgaris L. var Pinto to tobacco mosaic virus infection. The infection of Phaseolus vulgaris L. var Pinto with tobacco mosaic virus (TMV) results in the production of distinct necrotic lesions confining the virus to restricted areas of the leaf surface. Biochemical and histological changes in the leaf tissue as a result of infection have been described. Trace accumulations of fluorescent metabolites, detected prior to lesion expression represent metabolites produced, by the cell in response to virus infection. These substances, are considered to undergo oxidation and in diffusing into adjacent cells, react with cellular constituents causing the death of these cells. Such cellular necrosis in advance of infection effectively limits virus spread. Chromatographic studies on extracts from TMV infected Pinto bean leaf tissue suggests that a number of extra-fluorescent metabolites produced on lesion'expression represent end products of phenolic oxidation r,eactionsoccurring earlier in these cells. Inhibition of phenolic oxidation by ascorbate infiltration or elevated temperature treatment resulted in the absence of extra-fluorescent metabolites and the continued movement of virus in the absence of necrosis. Further studies with i ascorbate infiltration indicated that irreversible necrotic events were determined as early as 12 tci 18 hrs after viral inoculation. Histochemical tests indicated that callose formation was initiated at this time, and occurred in response to necrotisation. Inhibition of necrosis by either ascorbate infiltration or elevated temp8rature treatment resulted in the absence of callose deposition. Scanning electron'micrographs of infected tissue revealed severe epidermal and palisade cell damage. Histochemical tests indicated extensive callose formation in cells bordering the lesion, and suggested the role of callose iTh the blockage of intercellular connections limiting virus movement. The significance of these cellular changes is discussed. ii