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.

Monooxygenase Diversity and Oxygen Activation within Polyketide Antibiotic Biosynthesis

Molecular oxygen (O2) is a key component in cellular respiration and aerobic life. Through the redox potential of O2, the amount of free energy available to organisms that utilize it is greatly increased. Yet, due to the nature of the O2 electron configuration, it is non-reactive to most organic molecules in the ground state. For O2 to react with most organic compounds it must be activated. By activating O2, oxygenases can catalyze reactions involving oxygen incorporation into organic compounds. The oxygen activation mechanisms employed by many oxygenases to have been studied, and they often include transition metals and selected organic compounds. Despite the diversity of mechanisms for O2 activation explored in this thesis, all of the monooxygenases studied in the experimental part activate O2 through a transient carbanion intermediate. One of these enzymes is the small cofactorless monooxygenase SnoaB. Cofactorless monooxygenases are unusual oxygenases that require neither transition metals nor cofactors to activate oxygen. Based on our biochemical characterization and the crystal structure of this enzyme, the mechanism most likely employed by SnoaB relies on a carbanion intermediate to activate oxygen, which is consistent with the proposed substrate-assisted mechanism for this family of enzymes. From the studies conducted on the two-component system AlnT and AlnH, both the functions of the NADH-dependent flavin reductase, AlnH, and the reduced flavin dependent monooxygenase, AlnT, were confirmed. The unusual regiochemistry proposed for AlnT was also confirmed on the basis of the structure of a reaction product. The mechanism of AlnT, as with other flavin-dependent monooxygenases, is likely to involve a caged radical pair consisting of a superoxide anion and a neutral flavin radical formed from an initial carbanion intermediate. In the studies concerning the engineering of the S-adenosyl-L-methionine (SAM) dependent 4-O-methylase DnrK and the homologous atypical 10-hydroxylase RdmB, our data suggest that an initial decarboxylation of the substrate is catalyzed by both of these enzymes, which results in the generation of a carbanion intermediate. This intermediate is not essential for the 4-O-methylation reaction, but it is important for the 10-hydroxylation reaction, since it enables substrate-assisted activation of molecular oxygen involving a single electron transfer to O2 from a carbanion intermediate. The only role for SAM in the hydroxylation reaction is likely to be stabilization of the carbanion through the positive charge of the cofactor. Based on the DnrK variant crystal structure and the characterizations of several DnrK variants, the insertion of a single amino acid in DnrK (S297) is sufficient for gaining a hydroxylation function, which is likely caused by carbanion stabilization through active site solvent restriction. Despite large differences in the three-dimensional structures of the oxygenases and the potential for multiple oxygen activation mechanisms, all the enzymes in my studies rely on carbanion intermediates to activate oxygen from either flavins or their substrates. This thesis provides interesting examples of divergent evolution and the prevalence of carbanion intermediates within polyketide biosynthesis. This mechanism appears to be recurrent in aromatic polyketide biosynthesis and may reflect the acidic nature of these compounds, propensity towards hydrogen bonding and their ability to delocalize π-electrons.

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.

Engineering an improved adenovirus packaging cell line

Recombinant human adenovirus (Ad) vectors are being extensively explored for their use in gene therapy and recombinant vaccines. Ad vectors are attractive for many reasons, including the fact that (1) they are relatively safe, based on their use as live oral vaccines, (2) they can accept large transgene inserts, (3) they can infect dividing and postmitotic cells, and (4) they can be produced to high titers. However, there are also a number of major problems associated with Ad vectors, including transient foreign gene expression due to host cellular immune responses, problems with humoral immunity, and the creation of replication competent adenoviruses (RCA). Most Ad vectors contain deletions in the E1 region that allow for insertion of a transgene. However, the E1 gene products are required for replication and thus must be supplied in trans by a helper ceillille that will allow for the growth and packaging of the defective virus. For this purpose the 293 cell line (Graham et al., 1977) is used most often; however, homologous recombination between the vector and the cell line often results in the generation of RCA. The presence of RCA in batches of adenoviral vectors for clinical use is a safety risk because tlley . may result in the mobilization and spread of the replication-defective vector viruses, and in significant tissue damage and pathogenicity. The present research focused on the alteration of the 293 cell line such that RCA formation can be eliminated. The strategy to modify the 293 cells involved the removal of the first 380 bp of the adenovirus genome through the process of homologous recombination. The first step towards this goal involved identifying and cloning the left-end cellular-viral jUl1ction from 293 cells to assemble sequences required for homologous recombination. Polymerase chain reaction (PCR) was performed to clone the junction, and the clone was verified through sequencing. The plasn1id PAM2 was then constructed, which served as the targeting cassette used to modify the 293 cells. The cassette consisted of (1) the cellular-viral junction as the left-end region of homology, (2) the neo gene to use for positive selection upon tranfection into 293 cells, (3) the adenoviral genome from bp 380 to bp 3438 as the right-end region of homology, and (4) the HSV-tk gene to use for negative selection. The plasmid PAM2 was linearized to produce a double strand break outside the region of homology, and transfected into 293 cells using the calcium-phosphate technique. Cells were first selected for their resistance to the drug G418, and subsequently for their resistance to the drug Gancyclovir (GANC). From 17 transfections, 100 pools of G418f and GANCf cells were picked using cloning lings and expanded for screening. Genomic DNA was isolated from the pools and screened for the presence of the 380 bps using PCR. Ten of the most promising pools were diluted to single cells and expanded in order to isolate homogeneous cell lines. From these, an additional 100 G41Sf and GANef foci were screened. These preliminary screening results appear promising for the detection of the desired cell line. Future work would include further cloning and purification of the promising cell lines that have potentially undergone homologous recombination, in order to isolate a homogeneous cell line of interest.

Isolation of DNA sequences with homology to a degenerate FaRP oligonucleotide from the crayfish Procambarus clarkii

The neuropeptide Th1RFamide with the sequence Phe-Met-Arg-Phe-amide was originally isolated in the clam Macrocallista nimbosa (price and Greenberg, 1977). Since its discovery, a large family ofFl\1RFamide-related peptides termed FaRPs have been found to be present in all major animal phyla with functions ranging from modulation of neuronal activity to alteration of muscular contractions. However, little is known about the genetics encoding these peptides, especially in invertebrates. As FaRP-encoding genes have yet to be investigated in the invertebrate Malacostracean subphylum, the isolation and characterization ofFaRP-encoding DNA and mRNA was pursued in this project. The immediate aims of this thesis were: (1) to amplify mRNA sequences of Procambarus clarkii using a degenerate oligonucleotide primer deduced from the common amino acid sequence ofisolated Procambarus FaRPS, (2) to determine if these amplification products encode FaRP gene sequences, and (3) to create a selective cDNA library of sequences recognized by the degenerate oligonucleotide primer. The polymerase chain reaction - rapid amplification of cDNA ends (PCR-RACE) is a procedure in which a single gene-specific primer is used in conjunction with a generalized 3' or 5' primer to amplify copies ofthe region between a single point in the transcript and the 3' or 5' end of cDNA of interest (Frohman et aI., 1988). PCRRACE reactions were optimized with respect to primers used, buffer composition, cycle number, nature ofgenetic substrate to be amplified, annealing, extension and denaturation temperatures and times, and use of reamplification procedures. Amplification products were cloned into plasmid vectors and recombinant products were isolated, as were the recombinant plaques formed in the selective cDNA library. Labeled amplification products were hybridized to recombinant bacteriophage to determine ligated amplification product presence. When sequenced, the five isolated PCR-RACE amplification products were determined not to possess FaRP-encoding sequences. The 200bp, 450bp, and 1500bp sequences showed homology to the Caenorhabditis elegans cosmid K09A11, which encodes for cytochrome P450; transfer-RNA; transposase; and tRNA-Tyr, while the 500bp and 750bp sequences showed homology with the complete genome of the Vaccinia virus. Under the employed amplification conditions the degenerate oligonucleotide primer was observed to bind to and to amplify sequences with either 9 or 10bp of 17bp identity. The selective cDNA library was obselVed to be of extremely low titre. When library titre was increased, white. plaques were isolated. Amplification analysis of eight isolated Agt11 sequences from these plaques indicated an absence of an insertion sequence. The degenerate 17 base oligonucleotide primer synthesized from the common amino acid sequence ofisolated Procambarus FaRPs was thus determined to be non-specific in its binding under the conditions required for its use, and to be insufficient for the isolation and identification ofFaRP-encoding sequences. A more specific primer oflonger sequence, lower degeneracy, and higher melting temperature (TJ is recommended for further investigation into the FaRP-encoding genes of Procambarlls clarkii.

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 nucleophillic substitution of various chloroanthraquinones of their possible rearrangement via aryne intermediates

The work in this thesis deals mainly with nucleophilic substitution of chloroanthraquinones as a route to various starting materials which might rearrange, via aryne intermediates to afford fused-ring heterocy1ic carboxylic acids. 1-Amino-5-chloroanthraquinone was successfully prepared by reacting 1,5-dichloroanthraquinone with sodium aZide in ref1uxing dimethylsulfoxide (DMSO). It could also be prepared from the same starting material by reaction with ammonia (gas) in DMSO in the presence of potassium fluoride. Treatment of l-amino-5-chloroanthraquinone with potassium amide in liquid ammonia or with potassium t-butoxide in t-butylbenzene returned mainly starting material, although in the latter case some 1-amino-5-hydroxyanthraquinone was also isolated. 1-Hydroxy-5-chloroanthraquinone was ultimately prepared by diazotization of the amino-analog. It was recovered almost quantitatively after treatmenu'with potassium amide in liquid ammonia. The reaction with potassium t-butoxide in t-buty1benzene was anomalous and gave 1-hydroxyanthraquinone as the only iso1able product. Acridines were successfully prepared by the action of 70% sulfuric acid on 1,5-bis(p-toluidino)-anthraquinone and 1-p-toluidino-5- ch10roanthraquinone, and in the latter case, cleavage to give an acridinecarboxylate was attempted. Substituted anthraquinones reacted with sodium azide in sulfuric acid to give azepindiones by -NH insertion. Methods for separating and identifying isomeric mixtures of these compounds were examined. Attempted decarbonylation of selected azepindiones to give acridones gave mainly what were thought to be amino-benzophenone derivatives. Chloroanthraquinones were found to react with hexamethylphosphoramide (HMPA) to give mixtures of the dimethylamino- and methylaminoderivatives. Under the same conditions halogeno-nitrobenzenes and nitrophenols were substituted to give the appropriate dimethyl aminobenzenes, except in two cases. 3-Chloronitrobenzene reacted anomalously to give a small amount of 3,3'-dichloroazobenzene and a trace of 4-dimethylamino-nitrobenzene. Pentachlorophenol reacted to give a pentachlorophenylphosphorodiamidate in good yield.

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