Quantitative refinement of reaction-based biomodels

In the field of molecular biology, scientists adopted for decades a reductionist perspective in their inquiries, being predominantly concerned with the intricate mechanistic details of subcellular regulatory systems. However, integrative thinking was still applied at a smaller scale in molecular biology to understand the underlying processes of cellular behaviour for at least half a century. It was not until the genomic revolution at the end of the previous century that we required model building to account for systemic properties of cellular activity. Our system-level understanding of cellular function is to this day hindered by drastic limitations in our capability of predicting cellular behaviour to reflect system dynamics and system structures. To this end, systems biology aims for a system-level understanding of functional intraand inter-cellular activity. Modern biology brings about a high volume of data, whose comprehension we cannot even aim for in the absence of computational support. Computational modelling, hence, bridges modern biology to computer science, enabling a number of assets, which prove to be invaluable in the analysis of complex biological systems, such as: a rigorous characterization of the system structure, simulation techniques, perturbations analysis, etc. Computational biomodels augmented in size considerably in the past years, major contributions being made towards the simulation and analysis of large-scale models, starting with signalling pathways and culminating with whole-cell models, tissue-level models, organ models and full-scale patient models. The simulation and analysis of models of such complexity very often requires, in fact, the integration of various sub-models, entwined at different levels of resolution and whose organization spans over several levels of hierarchy. This thesis revolves around the concept of quantitative model refinement in relation to the process of model building in computational systems biology. The thesis proposes a sound computational framework for the stepwise augmentation of a biomodel. One starts with an abstract, high-level representation of a biological phenomenon, which is materialised into an initial model that is validated against a set of existing data. Consequently, the model is refined to include more details regarding its species and/or reactions. The framework is employed in the development of two models, one for the heat shock response in eukaryotes and the second for the ErbB signalling pathway. The thesis spans over several formalisms used in computational systems biology, inherently quantitative: reaction-network models, rule-based models and Petri net models, as well as a recent formalism intrinsically qualitative: reaction systems. The choice of modelling formalism is, however, determined by the nature of the question the modeler aims to answer. Quantitative model refinement turns out to be not only essential in the model development cycle, but also beneficial for the compilation of large-scale models, whose development requires the integration of several sub-models across various levels of resolution and underlying formal representations.

The use of polymerase chain reaction for early diagnosis of tuberculosis in Mycobacterium tuberculosis culture

Early diagnosis plays a vital role in controlling tuberculosis. The conventional methodology is slow, with results taking several weeks, in addition to having low sensitivity, especially in clinical paucibacillary samples. The objective of this study was to evaluate the use of polymerase chain reaction (PCR) on solid medium culture for a rapid diagnosis of tuberculosis, mainly in cases of negative sputum smears. Forty sputum samples were collected from inpatients with tuberculosis treated for less than 2 days. Bacilloscopy, PCR for sputum, culture on Löwestein-Jensen (LJ) solid medium, and daily PCR from culture were performed on each sample. DNA extracted from the BCG vaccine, which contains attenuated bacillus Calmette-Guérin, was used as the positive control. Smear microscopy showed 68.6% sensitivity, 80% specificity, 96% positive predictive value, and 26.7% negative predictive value, with culture on LJ medium as the gold standard. Culture at day 28 showed 74.3% sensitivity and 100% specificity. PCR of DNA extracted from sputum amplified a 1027-bp fragment of the 16s RNA gene, showing 22.9% sensitivity and 60% specificity. PCR performed with DNA extracted from daily culture showed that, from the 17th to the 40th day, the sensitivity (85.7%) and specificity (60%) were constant. We conclude that a 17-day culture is a good choice for rapid diagnosis and to interfere with the transmission chain of tuberculosis.

Identification of fungemia agents using the polymerase chain reaction and restriction fragment length polymorphism analysis

Prompt and specific identification of fungemia agents is important in order to define clinical treatment. However, in most cases conventional culture identification can be considered to be time-consuming and not without errors. The aim of the present study was to identify the following fungemia agents: Candida albicans, Candida parapsilosis, Candida tropicalis, Candida glabrata, Cryptococcus neoformans, Cryptococcus gattii, and Histoplasma capsulatum using the polymerase chain reaction and restriction fragment length polymorphism analysis (PCR/RFLP). More specifically: a) to evaluate 3 different amplification regions, b) to investigate 3 different restriction enzymes, and c) to use the best PCR/RFLP procedure to indentify 60 fungemia agents from a culture collection. All 3 pairs of primers (ITS1/ITS4, NL4/ITS5 and Primer1/Primer2) were able to amplify DNA from the reference strains. However, the size of these PCR products did not permit the identification of all the species studied. Three restriction enzymes were used to digest the PCR products: HaeIII, Ddel and Bfal. Among the combinations of pairs of primers and restriction enzymes, only one (primer pair NL4/ITS5 and restriction enzyme Ddel) produced a specific RFLP pattern for each microorganism studied. Sixty cultures of fungemia agents (selected from the culture collection of Fundação de Medicina Tropical do Amazonas - FMTAM) were correctly identified by PCR/RFLP using the prime pair NL4/ITS5 and Ddel. We conclude that the method proved to be both simple and reproducible, and may offer potential advantages over phenotyping methods.

Use of the polymerase chain reaction to detect Mycobacterium leprae in urine

Leprosy is an infectious disease caused by Mycobacterium leprae. The polymerase chain reaction (PCR) has been applied to detect M. leprae in different clinical samples and urine seems to be attractive for this purpose. PCR was used to improve the sensitivity for diagnosing leprosy by amplifying a 151-bp PCR fragment of the M. leprae pra gene (PCR-Pra) in urine samples. Seventy-three leprosy patients (39 males and 34 females, 14 to 78 years old) were selected for leprosy diagnosis at a reference laboratory in Maringá, PR, Brazil. Of these, 36 were under anti-leprosy multidrug therapy with dapsone and rifampicin for tuberculoid (TT) and dapsone, rifampicin and clofazimine for borderline (BB) and lepromatous (LL) forms. The control group contained 50 healthy individuals without any clinical history of leprosy. DNA isolated from leprosy patients’ urine samples was successfully amplified by PCR-Pra in 46.6% (34/73) of the cases. The positivity of PCR-Pra for patients with the TT form was 75% for both patients under treatment and non-treated patients (P = 0.1306). In patients with the LL form, PCR-Pra positivity was 52 and 30% for patients under treatment and non-treated patients, respectively (P = 0.2386). PCR-Pra showed a statistically significant difference in detecting M. leprae between the TT and LL forms of leprosy in patients under treatment (P = 0.0033). Although the current study showed that the proposed PCR-Pra has some limitations in the detection of M. leprae, this method has the potential to be a useful tool for leprosy diagnosis mainly in TT leprosy where the AFB slit-skin smear is always negative.

Effects of sciatic nerve transection on ultrastructure, NADPH-diaphorase reaction and serotonin-, tyrosine hydroxylase-, c-Fos-, glucose transporter 1- and 3-like immunoreactivities in frog dorsal root ganglion

Frogs have been used as an alternative model to study pain mechanisms. Since we did not find any reports on the effects of sciatic nerve transection (SNT) on the ultrastructure and pattern of metabolic substances in frog dorsal root ganglion (DRG) cells, in the present study, 18 adult male frogs (Rana catesbeiana) were divided into three experimental groups: naive (frogs not subjected to surgical manipulation), sham (frogs in which all surgical procedures to expose the sciatic nerve were used except transection of the nerve), and SNT (frogs in which the sciatic nerve was exposed and transected). After 3 days, the bilateral DRG of the sciatic nerve was collected and used for transmission electron microscopy. Immunohistochemistry was used to detect reactivity for glucose transporter (Glut) types 1 and 3, tyrosine hydroxylase, serotonin and c-Fos, as well as nicotinamide adenine dinucleotide phosphate diaphorase (NADPH-diaphorase). SNT induced more mitochondria with vacuolation in neurons, satellite glial cells (SGCs) with more cytoplasmic extensions emerging from cell bodies, as well as more ribosomes, rough endoplasmic reticulum, intermediate filaments and mitochondria. c-Fos immunoreactivity was found in neuronal nuclei. More neurons and SGCs surrounded by tyrosine hydroxylase-like immunoreactivity were found. No change occurred in serotonin- and Glut1- and Glut3-like immunoreactivity. NADPH-diaphorase occurred in more neurons and SGCs. No sign of SGC proliferation was observed. Since the changes of frog DRG in response to nerve injury are similar to those of mammals, frogs should be a valid experimental model for the study of the effects of SNT, a condition that still has many unanswered questions.

Relative performance of the two hands in simple and choice reaction time tasks

There is evidence that the left hemisphere is more competent for motor control than the right hemisphere. This study investigated whether this hemispheric asymmetry is expressed in the latency/duration of sequential responses performed by the left and/or right hands. Thirty-two right-handed young adults (16 males, 16 females; 18-25 years old) were tested in a simple or choice reaction time task. They responded to a left and/or right visual target by moving their left and/or right middle fingers between two keys on each side of the midline. Right hand reaction time did not differ from left hand reaction time. Submovement times were longer for the right hand than the left hand when the response was bilateral. Pause times were shorter for the right hand than the left hand, both when the responses were unilateral or bilateral. Reaction time results indicate that the putatively more efficient response preparation by the left hemisphere motor mechanisms is not expressed behaviorally. Submovement time and pause time results indicate that the putatively more efficient response execution by the left hemisphere motor mechanisms is expressed behaviorally. In the case of the submovements, the less efficient motor control of the left hand would be compensated by a more intense attention to this hand.

Bio and chemocatalysis for stereo and regioselective one-pot reaction applications

The development of cost efficient, selective and sustainable chemical processes for production of chiral building blocks is of great importance in synthetic and industrial organic chemistry. One way to reach these objectives is to carry out several reactions steps in one vessel at one time. Furthermore, when this kind of one-pot multi step reactions are catalyzed by heterogeneous chemo- and bio-catalysts, which can be separated from the reaction products by filtration, practical access to chiral small molecules for further utilization can be obtained. The initial reactions studied in this thesis are the two step dynamic kinetic resolution of rac-2-hydroxy-1-indanone and the regioselective hydrogenation of 1,2-indanedione. These reactions are then combined in a new heterogeneously catalyzed one-pot reaction sequence enabling simple recovery of the catalysts by filtration, facilitating simple reaction product isolation. Conclusively, the readily available 1,2-indanedione is by the presented one-pot sequence, utilizing heterogeneous enzyme and transition metal based catalysts, transferred with high regio- and stereoselectivity to a useful chiral vicinal hydroxyl ketone structure. Additional and complementary investigation of homogeneous half-sandwich ruthenium complexes for catalyzing the epimerization of chiral secondary alcohols of five natural products containing additional non-functionalized stereocenters was conducted. In principle, this kind of epimerization reactions of single stereocenters could be utilized for converting inexpensive starting materials, containing other stereogenic centers, into diastereomeric mixtures from which more valuable compounds can be isolated by traditional isolation techniques.

Use of the polymerase chain reaction for detection of Fusarium graminearum in bulgur wheat

The detection of mycotoxigenic fungi in foodstuff is important because their presence may indicate the possible associated mycotoxin contamination. Fusarium graminearum is a wheat pathogen and a producer of micotoxins. The polymerase chain reaction (PCR) has been employed for the specific identification of F. graminearum. However, this methodology has not been commonly used for detection of F. graminearum in food. Thus, the objective of the present study was to develop a molecular methodology to detect F. graminearum in commercial samples of bulgur wheat. Two methods were tested. In the first method, a sample of this cereal was contaminated with F. graminearum mycelia. The genomic DNA was extracted from this mixture and used in a F. graminearum specific PCR reaction. The F. graminearum species was detected only in samples that were heavily contaminated. In the second method, samples of bulgur wheat were inoculated on a solid medium, and isolates having F. graminearum culture characteristics were obtained. The DNA extracted from these isolates was tested in F. graminearum specific PCR reactions. An isolate obtained had its trichothecene genotype identified by PCR. The established methodology could be used in surveys of food contamination with F. graminearum.

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.

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.

Towards the enantioselective synthesis of lycoricidine alkaloids /

Two efficient, regio- and stereo controlled synthetic approaches to the synthesis of racemic analogs of pancratistatin have been accomplished and they serve as the model systems for the total synthesis of optically active 7-deoxy-pancratistatin. In the Diels-Alder approach, an efficient [4+2] cycloaddition of 3,4-methylenedioxyco- nitrostyrene with Danishefsky's diene to selectively form an exo-nitro adduct has been developed as the key step in the construction of the C-ring of the target molecule. In the Michael addition approach, the key step was a conjugate addition of an organic zinc-cuprate to the 3,4-methylenedioxy-(B-nitrostyrene, followed by a diastereocontroUed closure to form the cyclohexane C-ring of the target molecule via an intramolecular nitro-aldol cyclization on a neutral alumina surface. A chair-like transition state for such a cyclization has been established and such a chelation controlled transition state can be useful in the prediction of diastereoselectivity in other related 6-exo-trig nitroaldol reactions. Cyclization of the above products fi^om both approaches by using a Bischler-Napieralski type reaction afforded two lycoricidine derivatives 38 and 50 in good yields. The initial results from the above modeling studies as well as the analysis of the synthetic strategy were directed to a chiral pool approach to the total synthesis of optically active 7-deoxy-pancratistatin. Selective monsilylation and iodination of Ltartaric acid provided a chiral precursor for the proposed key Michael transformation. The outlook for the total synthesis of 7-deoxy-pancratistatin by this approach is very promising.A concise synthesis of novel designed, optically pure, Cz-symmetrical disulfonylamide chiral ligands starting from L-tartaric acid has also been achieved. This sequence employs the metallation of indole followed by Sfj2 replacement of a dimesylate as the key step. The activity for this Cz-symmetric chiral disulfonamide ligand in the catalytic enantioselective reaction has been confirmed by nucleophilic addition to benzaldehyde in the disulfonamide-Ti (0-i-Pr)4-diethylzinc system with a 48% yield and a 33% e.e. value. Such a ligand tethered with a suitable metal complex should be also applicable towards the total synthesis of 7-deoxy-pancratistatin.

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.

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.