Dr Brendan Howlin

Senior Lecturer in Computational Chemistry
+44 (0)1483 686248
23 AZ 02
9 to 5:30

Academic and research departments

Department of Chemistry.


Areas of specialism

Molecular Modelling

University roles and responsibilities

  • Chemistry Athena Swan Lead
  • Member of IT Governance Committee
  • Member of FEPS Equality and Diversity Committee

    My qualifications

    BSc Biological Chemistry
    University of Essex
    PhD Chemistry
    University of Essex


    Research interests

    Research projects

    Research collaborations


    Postgraduate research supervision

    Completed postgraduate research projects I have supervised

    Courses I teach on



    My publications


    Howlin BJ, Crawford AO, Cavalli G, Hamerton I (2016)Investigation of structure property relationships in liquid processible, solvent free, thermally stable bismaleimide-triazine (BT) resins, In: Reactive and Functional Polymers Elsevier
    Three cyanate ester monomer or oligomer species: 2,2-bis(4-cyanatophenyl)propane 1, 1-1-bis(4-cyanatophenyl)ethane (2), and the oligomeric phenolic cyanate (PrimasetTM PT30) (3), are blended in various ratios with bis(4-maleimidophenyl)methane, (4), to form binary and ternary mixtures (11 in total) and cured, in the absence of catalysts (3 K/min to 150 °C + 1 hour; 3 K/min to 200 °C + 3 hours), followed by a post cure (3 K/min to 260 °C + 1 hour). The use of liquid monomer, (2), offers the possibility of liquid processing in blends containing minority compositions of bismaleimide. Glycidylmethacrylate is explored as a reactive diluent (2.5-10 wt %) to linked interpenetrating network polymer structures comprising cyanate ester and bismaleimide components with glass transition temperatures of 267-275 ºC, depending on composition; the onset of thermo-oxidative degradation ranges from 386-397 ºC. When a binary blend of (2) and (3) (with the former in the minority) is co-cured with (4), an excellent balance of properties is achieved with liquid processing, a Tg greater than 400 C and onset of degradation of 425 ºC in static air. Kinetic analysis of DSC data using Ozawa and Kissinger methods yield activation energies of between 107- 112 kJ/mole for a binary blend of (1)90-(4)10, which is in good agreement with literature. Molecular dynamics simulation of the same blend in cured form gave a simulated glass transition temperature of 250 C that is in very close agreement with empirical DMTA data.
    Hamerton I, Heald CR, Howlin BJ (1996)Molecular modelling of the physical and mechanical properties of two polycyanurate network polymers, In: JOURNAL OF MATERIALS CHEMISTRY6(3)pp. 311-314
    Hamerton I, Howlin BJ, Jones JR, Liu SY, Barton JM (1996)Effect of complexation with copper(II) on cured neat resin properties of a commercial epoxy resin using modified imidazole curing agents, In: JOURNAL OF MATERIALS CHEMISTRY6(3)pp. 305-310
    This study uses Molecular Operating Environment (MOE) to generate models to calculate the char yield of polybenzoxazines. A series of benzoxazine monomers were constructed to which a variety of parameters relating to the structure (e.g. water accessible surface, negative van der Waals surface area and hydrophobic volume, etc.) were obtained and a quantitative structure property relationships (QSPR) model was generated. The model was used to generate data for a new benzoxazine monomer and a comparison was made of predictions based on the QSPR models with the experimental data. This study shows the quality of predictive models and confirms how useful computational screening is prior to synthesis. In order to do that, the QSPR models were tested over a series of internal and external validation tests to explore their internal and external predictivity, prior to experimental validations which were performed later and reported in Chapter 7. The internal and external validations found out that the discrepancy in the general model (GM) which was initially thought to be a drawback to the model’s performance was actually not, as it does not compromise the model’s prediction accuracy, both internally and externally. The validation process also found that one of the structure-specific models, Ph-M (aniline-based benzoxazines) is externally predictive whilst another structure-specific model, the Ace-M (acetylenic-based polybenzoxazines) is not internally and externally predictive due to the too small training set that affects its predictivity performance. An acetylenic-based polybenzoxazine, poly(BA-apa) and a benzylamine-based polybenzoxazine, poly(BO-ba) have been successfully synthesised in this work. Both materials have been characterised using Fourier Transform – Infra Red Spectroscopy (FT-IR), Nuclear Magnetic Resonance (NMR) spectroscopy (both 1H and 13C) and Liquid Chromatography-Mass Spectrometry (LC-MS) to confirm their structures. These materials were analysed using Differential Scanning Calorimetry (DSC) to study their polymerisation behaviour and were later cured and taken further to Thermogravimetric Analysis (TGA) in order to investigate their thermal properties and the amount of char yield formed upon heating at 800 oC under an inert (nitrogen) atmosphere – which then will be used for experimental validation of the QSPR models. The study of DSC thermograms showed that both polymers exhibit a distinct polymerisation behaviour e.g. BA-apa went through two polymerisation reactions simultaneously (the oxazine ring opening polymerisation and the acetylene addition reaction) whilst BO-ba only polymerised via the ring opening reaction from the oxazine rings. It was also found that BA-apa has a lower polymerisation activation energy, consistent to its lower polymerisation temperature in comparison to the BO-ba. TGA analysis revealed that poly(BO-ba) formed an average of 44.35 % char yield and poly(BA-apa) on the other hand formed approximately 10 % higher char which is 56.28 %. The analysis also discovered that poly(BA-apa) synthesised in this work formed 15 % less char yield than previously reported in the literature (56.28 % vs. 71 %1) due to the shorter curing schedule. The final QSPR validation which is the experimental validation found that the char yield of poly(BO-ba) was predicted very well within 5-7 % error by both GM model and Ph-M. Ace-M which was reported earlier as not internally and externally predictive, has made a nearly accurate prediction towards the char yield of poly(BA-apa), close to the literature value of 71 %. The GM model has also made a close prediction to the Ace-M model, but these predictions deviated 15-17 % from the experimental poly(BA-apa) char yield measured in this work.
    Binks Fiona C, Cavalli Gabriel, Henningsen Michael, Howlin Brendan, Hamerton Ian (2018)Examining the effects of storage on the initiation behaviour of ionic liquids towards the cure of epoxy resins, In: Reactive and Functional Polymers133pp. 9-20 Elsevier
    Four structurally related ionic liquids (1-ethyl-3-methylimidazolium acetate, 1-ethyl-3-methylimidazolium diethyl phosphate, 1-ethyl-3-methylimidazolium dicyanamide, and 1-ethyl-3-methylimidazolium thiocyanate) are examined for their storage characteristics and its effect on their ability to initiate the cure of epoxy resins. At ambient temperature, epoxy formulations containing 1-ethyl-3-methylimidazolium acetate and 1-ethyl-3-methylimidazolium thiocyanate display marked colour changes to yield dark red samples with greatly increased viscosity after one day; after six days both samples have undergone vitrification. The epoxy formulation containing 1-ethyl-3-methylimidazolium acetate continued to polymerise even at sub-zero temperatures. Storage in dark bottles retarded the reaction during the 30-minute period that the sample is removed from the freezer prior to an aliquot being taken, but once the autocatalytic low temperature reaction has started, the dark glass no longer provides effective protection. Samples of 1-ethyl-3-methylimidazolium dicyanamide/epoxy were also stored and sampled in the same manner, but no differences were exhibited between the samples in clear and dark brown glass bottles. Infrared and nuclear magnetic resonance studies confirmed that the hygroscopic ionic liquids pick up water readily (coordinating to the H atom at the 2-position on the imidazolium ring), but once dried the initiating ability is lost.
    Hamerton I, Howlin BJ, Smith ER (2013)Using POSS reagents to reduce hydrophobic character in polypropylene nanocomposites, In: Journal of Materials Chemistry A1(41)pp. 12971-12980 Royal Society of Chemistry
    Hamerton I, Howlin BJ (1999)Computer modelling and chemical kinetics of the cure of aerospace resin systems, In: AIRCRAFT ENGINEERING AND AEROSPACE TECHNOLOGY71(5)pp. 470-478 MCB UNIV PRESS LTD
    Hamerton Ian, Howlin Brendan, Tilbrook DA (2013)Examining thermal stability and structure property relationships in coatings based on linear aromatic poly(methoxy-thiocyanurate)s, In: Polymer Degradation and Stability98(11)pp. 2201-2207
    A series of copoly(methoxy-thiocyanurate)s is prepared in good yield and purity, and fully characterised. Many of the resulting polymers, formed at room temperature using phase transfer catalysis, can be cast into films with good resilience and thermal stability (some examples suffer practically no mass loss when held isothermally at 190 °C and only display appreciable losses when held continuously at 225 °C). Char yields of 61-64% are achieved in nitrogen depending on backbone structure. Some problems were encountered with solubility, particularly with copolymers, which limited molecular weights analysis, but values of Mn = 7000-10,000 g mol-1 were obtained for the polycyanurate and polythiocyanurate homopolymers. DSC reveals polymerisation exotherms with maxima at 197-207 °C (ΔHp = 39-48 kJ/mol), which are believed to be due to isomerisation of the (activation energies span 172-205 kJ/mol), since X-ray powder diffraction measurements reveal no evidence of crystalline structure in the resulting product.
    Crawford AO, Hamerton I, Cavalli G, Howlin BJ (2012)Quantifying the Effect of Polymer Blending through Molecular Modelling of Cyanurate Polymers., In: PLoS One7(9)pp. e44487-?
    Modification of polymer properties by blending is a common practice in the polymer industry. We report here a study of blends of cyanurate polymers by molecular modelling that shows that the final experimentally determined properties can be predicted from first principles modelling to a good degree of accuracy. There is always a compromise between simulation length, accuracy and speed of prediction. A comparison of simulation times shows that 125ps of molecular dynamics simulation at each temperature provides the optimum compromise for models of this size with current technology. This study opens up the possibility of computer aided design of polymer blends with desired physical and mechanical properties.
    Hamerton I, Howlin BJ (2017)Modern Developments Using Molecular Simulation to Predict the Physical and Mechanical Properties of Polybenzoxazines, In: Ishida H, Froimowicz P (eds.), Advanced and Emerging Polybenzoxazine Science and Technology(8)pp. 111-129 Elsevier
    The use of molecular simulation techniques in connection with structural polymers is now becoming more widely accepted and is employed as a legitimate way to not only visualize complex structures, but also to replicate empirically determined parameters. In this chapter, the authors present the state-of-the-art uses of molecular simulation (molecular mechanics and molecular dynamics) based on some 24 years of practical experience in the field.
    Hamerton I, Howlin BJ (2017)Predictive Methodology and Properties of Polybenzoxazines, In: Ishida H, Froimowic P (eds.), Advanced and Emerging Polybenzoxazine Science and Technology(9)pp. 131-145 Elsevier
    In order to understand the structural basis of the observed properties of benzoxazine materials, the technique of quantitative structure-property relationships can be applied to sets of monomers or polymers and the understanding gained can then be used to design new monomers or polymers with the desired properties. In the chapter, this technique is applied to some datasets to predict char yield and glass transition temperature, and to illustrate that, despite the limited nature of the datasets currently available, useful predictions can be made.
    Hay JN, Hamerton I, Howlin BJ, Howgate GJ, O'Gara PM (2001)Electron-beam cure of phenylethynyl-terminated oligomers., In: Repecka L, Saremi FF, (eds.), 2001: A MATERIALS AND PROCESSES ODYSSEY, BOOKS 1 AND 246pp. 2140-2146
    Bisbenzoxazines form highly crosslinked, network structures with relatively high Tg (180oC+) and thermal oxidative stability but inherent brittleness (G1C = 168 J/m2) - limiting their applications. We have used synthetic polymer chemistry and molecular simulation to probe the structure/property relationships in model and commercial systems with the aim of understanding the factors governing toughness in these materials. Our previous experience, gained over 20 years, leads us to ensure that the molecular model represents as accurately as possible the real network, e.g. incorporating molecular defects derived from empirical data. We report the preparation and characterisation of selected bisbenzoxazine monomers chosen to exemplify both brittle and tough behaviour to provide a wide spectrum of mechanical characteristics and produce a wider understanding of the factors affecting brittleness. This enables us to design and execute the synthesis of new monomers with the potential for enhanced fracture toughness based on rational molecular design.
    Crawford AO, Howlin BJ, Cavalli G, Hamerton I (2012)Examining the thermo-mechanical properties of novel cyanate ester blends through empirical measurement and simulation, In: Reactive and Functional Polymers72(9)pp. 596-605
    Three cyanate ester monomer or oligomer species: 2,2-bis(4-cyanatophenyl) propane 1, 1-1-bis(4-dicyanatophenyl)ethane 2, and the oligomeric phenolic cyanate (Primaset™ PT30) 3, are blended in various ratios to form binary mixtures, formulated with copper(II) acetylacetonate (200 ppm) in dodecylphenol (1% w/v active copper suspension) and cured (2 K/min to 150 °C + 1 h; 2 K/min to 200 °C + 3 h) followed by a post cure (2 K/min to 260 °C + 1 h). Thermal analysis using DSC reveals good agreement with literature data for the homopolymers: typical polymerisation enthalpies of ca. 97-98 kJ/mol. cyanate are obtained for 1 and 2, with slightly lower values (ca. 80-90 kJ/mol.) obtained for Primaset™ PT30. DMTA data show the possibility of using binary blends of the polymers to yield novel materials with similar thermal and mechanical properties to Primaset™ PT30, while improving the processability of the more highly aromatic oligomer. Two of the homopolymers (1 and 2) and a binary (1:1) blend of the same were simulated. Molecular dynamics experiments reveal good agreement with empirical data generated using DSC, DMTA and TGA. © 2012 Elsevier Ltd. All rights reserved.
    Mortimer S, Hamerton I, Hay JN, Howlin BJ, Jepson P (2001)The development of controllable complex curing agents for epoxy resins. I. Preparation, characterization, and storage behavior of transition metal-diamine complexes, In: JOURNAL OF APPLIED POLYMER SCIENCE80(9)pp. 1489-1503 JOHN WILEY & SONS INC
    Molecular simulation is becoming an important tool for both understanding polymeric structures and predicting their physical and mechanical properties. In this study, temperature ramped molecular dynamics simulations are used to predict two physical properties (i.e., glass transition temperature and thermal degradation temperature) of a previously synthesised and published telechelic benzoxazine. Plots of simulated density versus temperature show decreases in density within the same temperature range as experimental values for the thermal degradation. The predicted value for the thermal degradation temperature for the cured polybenzoxazine based on the telechelic polyetherketone (PEK) monomer was ca. 400°C, in line with the experimental thermal degradation temperature range of 450°C to 500°C. Mechanical Properties of both the unmodified PEK and the telechelic benzoxazines are simulated and compared to experimental values (where available). The introduction of the benoxazine moieties are predicted to increase the elastic moduli in line with the increase of crosslinking in the system.
    Hamerton I, Barton JM, Chaplin A, Howlin BJ, Shaw SJ (2001)The development of novel functionalised aryl cyanate esters. Part 2. Mechanical properties of the polymers and composites, In: POLYMER42(6)pp. 2307-2319 ELSEVIER SCI LTD
    Hamerton I, Hay JN, Herman H, Howlin BJ, Jepson P, Gillies DG (2002)The development of controllable complex curing agents for epoxy resins. II. Examining the dissociation and thermal behavior of transition metal-diamine complex-epoxy blends, In: JOURNAL OF APPLIED POLYMER SCIENCE84(13)pp. 2411-2424 JOHN WILEY & SONS INC
    Hall SA, Hamerton I, Howlin BJ, Mitchell AL (2008)Validating software and force fields for predicting the mechanical and physical properties of poly(bisbenzoxazine)s, In: MOLECULAR SIMULATION34(10-15)PII 90pp. 1259-1266 TAYLOR & FRANCIS LTD
    Binks Fiona C., Cavalli Gabriel, Henningsen Michael, Howlin Brendan, Hamerton Ian (2019)On the use of benzaldehyde to improve the storage stability of one-pot, epoxy ionic liquid formulations, In: EUROPEAN POLYMER JOURNAL112pp. 126-136 PERGAMON-ELSEVIER SCIENCE LTD
    A series of adducts were prepared based on the reaction of 1-ethyl-3-methylimidazolium acetate and benzal- dehyde in various stoichiometries (from equimolar reaction to benzaldehyde in 10-fold excess) and the resulting adducts were characterized using nuclear magnetic resonance spectroscopy (¹H, ¹³C, DEPT, and HQSC experi- ments). Differential scanning calorimetry was used to examine the initiating behaviour of the adducts towards mono- and di-functional epoxy resins and the data were used to determine kinetic parameters for the poly- merization. The lower temperature peak, due to carbene formation, is sensitive to adduct concentration; the residual ionic liquid in the adduct mixture contributes towards the initiation of the curing reaction. When a monofunctional epoxy and the 1:1 adduct was subjected to a 2-week period of storage at room temperature and sub-zero temperatures in the freezer, the profiles of the thermograms for the frozen samples do not change considerably over the storage period and the formulation retains a light yellow colour (rather than the viscous, dark red appearance of the formulation stored at room temperature).
    This study uses the Molecular Operating Environment software (MOE) to generate models to calculate the char yield of polybenzoxazines (PBz). A series of benzoxazine (Bz) monomers were constructed to which a variety of parameters relating to the structure (e.g., water accessible surface, negative van der Waals surface area and hydrophobic volume, etc.) were obtained and a quantitative structure property relationships (QSPR) model was generated. The model was used to generate data for new Bz monomers with desired properties and a comparison was made of predictions based on the QSPR model with the experimental data. This study shows the quality of predictive models and confirms how useful computational screening is prior to synthesis.
    Hamerton I, Hay JN, Howlin BJ, Jones JR, Lu SY, Webb GA, Bader MG (1997)A tritium-hydrogen exchange procedure for determining the concentration of exchangeable sites on PAN based carbon fibres, In: HIGH PERFORMANCE POLYMERS9(3)pp. 281-290 IOP PUBLISHING LTD
    Allington RD, Attwood D, Hay JN, Howlin BJ, Hamerton I (2004)Developing improved models of oxidatively treated carbon fibre surfaces, using molecular simulation, In: COMPOSITES PART A-APPLIED SCIENCE AND MANUFACTURING35(10)pp. 1161-1173 ELSEVIER SCI LTD
    Mhlanga P, Wan Hassan WA, Hamerton I, Howlin BJ (2013)Using combined computational techniques to predict the glass transition temperatures of aromatic polybenzoxazines., In: PLoS One8(1)pp. e53367-? PUBLIC LIBRARY SCIENCE
    The Molecular Operating Environment software (MOE) is used to construct a series of benzoxazine monomers for which a variety of parameters relating to the structures (e.g. water accessible surface area, negative van der Waals surface area, hydrophobic volume and the sum of atomic polarizabilities, etc.) are obtained and quantitative structure property relationships (QSPR) models are formulated. Three QSPR models (formulated using up to 5 descriptors) are first used to make predictions for the initiator data set (n = 9) and compared to published thermal data; in all of the QSPR models there is a high level of agreement between the actual data and the predicted data (within 0.63-1.86 K of the entire dataset). The water accessible surface area is found to be the most important descriptor in the prediction of T(g). Molecular modelling simulations of the benzoxazine polymer (minus initiator) carried out at the same time using the Materials Studio software suite provide an independent prediction of T(g). Predicted T(g) values from molecular modelling fall in the middle of the range of the experimentally determined T(g) values, indicating that the structure of the network is influenced by the nature of the initiator used. Hence both techniques can provide predictions of glass transition temperatures and provide complementary data for polymer design.
    Allington RD, Attwood D, Hamerton I, Hay JN, Howlin BJ (2001)New force-field parameters for use in molecular simulations of s-triazine and cyanurate-containing systems. 1 - derivation and molecular structure synopsis, In: COMPUTATIONAL AND THEORETICAL POLYMER SCIENCE11(6)pp. 467-473 ELSEVIER SCI LTD
    Hamerton I, Howlin BJ, Mitchell AL, McNamara LT, Takeda S (2012)Systematic examination of thermal, mechanical and dielectrical properties of aromatic polybenzoxazines, In: Reactive and Functional Polymers72(10)pp. 736-744
    A series of bis-benzoxazines, prepared in high yield and purity using two synthetic procedures, is reported. Differential scanning calorimetry reveals similar temperatures for the onset of polymerisation (162-180 °C); the higher values representing monomers containing polar bridges or rigid backbones. Dynamic viscoelasticity data reveal glass transition temperatures for the polybenzoxazines ranging from 187 °C to 235 °C; a fluorinated polybenzoxazine consistently yields the highest T of the polymers studied. The latter is interesting since it is superior to many commercial benzoxazines with a relatively high T (235 °C), flexural modulus (5.0 GPa) and flexural strength (146.7 MPa), but coupled with a breaking strain (3.06%) that is uncharacteristically high for polybenzoxazines. The incorporation of fluorine results in a low dielectric loss properties (D = 3.71-4.12 at 10 MHz, D = 0.0109 - 0.0980 at 10 MHz), which are comparable with commercial polybenzoxazines, FR4 and aerospace epoxy resins and superior to commercial bismaleimides.
    Purse Marcus, Edmund Grace, Hall Stephen, Howlin Brendan, Hamerton Ian, Till Stephen (2019)Reactive Molecular Dynamics Study of the Thermal Decomposition of Phenolic Resins, In: Journal of Composites Science3(2) MDPI
    The thermal decomposition of polyphenolic resins was studied by reactive molecular dynamics (RMD) simulation at elevated temperatures. Atomistic models of the polyphenolic resins to be used in the RMD were constructed using an automatic method which calls routines from the software package Materials Studio. In order to validate the models, simulated densities and heat capacities were compared with experimental values. The most suitable combination of force field and thermostat for this system was the Forcite force field with the Nosé–Hoover thermostat, which gave values of heat capacity closest to those of the experimental values. Simulated densities approached a final density of 1.05–1.08 g/cm3 which compared favorably with the experimental values of 1.16–1.21 g/cm3 for phenol-formaldehyde resins. The RMD calculations were run using LAMMPS software at temperatures of 1250 K and 3000 K using the ReaxFF force field and employing an in-house routine for removal of products of condensation. The species produced during RMD correlated with those found experimentally for polyphenolic systems and rearrangements to form cyclopropane moieties were observed. At the end of the RMD simulations a glassy carbon char resulted.
    Tilbrook DA, Howlin BJ, Hamerton I (2001)Assessment of two methods for application in the prediction of the infrared spectra of polymers, In: COMPUTATIONAL AND THEORETICAL POLYMER SCIENCE11(4)pp. 287-302 ELSEVIER SCI LTD
    Edmund GH, Lewis DF, Howlin BJ (2013)Modelling species selectivity in rat and human cytochrome P450 2D enzymes., In: PLoS One8(5)pp. e63335-?
    Updated models of the Rat Cytochrome P450 2D enzymes are produced based on the recent x-ray structures of the Human P450 2D6 enzyme both with and without a ligand bound. The differences in species selectivity between the epimers quinine and quinidine are rationalised using these models and the results are discussed with regard to previous studies. A close approach to the heme is not observed in this study. The x-ray structure of the enzyme with a ligand bound is shown to be a better model for explaining the observed experimental binding of quinine and quinidine. Hence models with larger closed binding sites are recommended for comparative docking studies. This is consistent with molecular recognition in Cytochrome P450 enzymes being the result of a number of non-specific interactions in a large binding site.
    Howlin Brendan J, Hamerton Ian, Hall Stephen A, Mitchell Amy L, McNamara Lisa (2011)Using Molecular Simulation to Predict the Physical and Mechanical Properties of Polybenzoxazines, In: Ishida H, Agag T (eds.), Handbook of Benzoxazine Resinsn/app. 127-142 Elsevier
    The aim of this chapter is to introduce the reader to the practical applications of modern molecular simulation techniques with literature examples drawn specifically from the field of polybenzoxazine research. The increases in computational power ensure that it is possible to apply molecular mechanics and molecular dynamics techniques to the visualisation and simulation of comparatively large model structures comprising in some cases more than six thousand atoms (constructed from a repeat unit containing ca. 250 atoms). This, in turn, offers the potential to replicate a variety of physical and mechanical characteristics with a high degree of accuracy and precision. However, the apparent ease with which modelling may be carried out using modern software is beguiling; the need to validate simulations with real, empirical data is essential to ensure that the researcher obtains meaningful results.
    Phagocyte superoxide production by a multicomponent NADPH oxidase is important in host defense against microbial invasion. However inappropriate NADPH oxidase activation causes inflammation. Endothelial cells express NADPH oxidase and endothelial oxidative stress due to prolonged NADPH oxidase activation predisposes many diseases. Discovering the mechanism of NADPH oxidase activation is essential for developing novel treatment of these diseases. The p47(phox) is a key regulatory subunit of NADPH oxidase; however, due to the lack of full protein structural information, the mechanistic insight of p47(phox) phosphorylation in NADPH oxidase activation remains incomplete. Based on crystal structures of three functional domains, we generated a computational structural model of the full p47(phox) protein. Using a combination of in silico phosphorylation, molecular dynamics simulation and protein/protein docking, we discovered that the C-terminal tail of p47(phox) is critical for stabilizing its autoinhibited structure. Ser-379 phosphorylation disrupts H-bonds that link the C-terminal tail to the autoinhibitory region (AIR) and the tandem Src homology 3 (SH3) domains, allowing the AIR to undergo phosphorylation to expose the SH3 pocket for p22(phox) binding. These findings were confirmed by site-directed mutagenesis and gene transfection of p47(phox-/-) coronary microvascular cells. Compared with wild-type p47(phox) cDNA transfected cells, the single mutation of S379A completely blocked p47(phox) membrane translocation, binding to p22(phox) and endothelial O2 (⨪) production in response to acute stimulation of PKC. p47(phox) C-terminal tail plays a key role in stabilizing intramolecular interactions at rest. Ser-379 phosphorylation is a molecular switch which initiates p47(phox) conformational changes and NADPH oxidase-dependent superoxide production by cells.
    Hamerton I, Howlin BJ, Pratt CR (2010)Simulation of the free energy of mixing for blend components in a new family of flexible polycyanurates, In: Polymer51(24)pp. 5857-5868 Elsevier
    The miscibility of the dicyanate of bisphenol A, with three dicyanate monomers, with aryl/alkylene ether backbones is studied at different compositions of a binary blend. Solubility parameters are calculated for dicyanate monomers and selected oligomers using the methods of Small and Fedors to predict compatibility. The results are evaluated in terms of the accuracy of the model in reproducing observed data. Gibbs free energy of mixing (ΔGmix) values for selected blends are calculated using the BLENDS module of Cerius2. Empirical data (HPLC and MS) are used to inform the construction of selected models to represent different stages of polymer conversion. DMTA analysis is performed to examine the thermo-mechanical properties of the resulting blends and compared with the simulated blend data.
    Mooring Lyndsey, Thompson Scott, Hall Stephen A., Pani Silvia, Zioupos Peter, Swan Martin, Stone Corinne, Howlin Brendan, Hamerton Ian (2018)‘Phoenix polymers’: fire induced nanohardness in fibril-forming aromatic cyanate esters, In: RSC Advances8(63)pp. 36264-36271 Royal Society of Chemistry
    For the first time we present nanoindentation analysis of charred, cured aromatic cyanate esters, which exhibit outstanding mechanical properties when analysed under applied loads of 0.1–300 mN. Following charring (900 C for 10 minutes to achieve graphitised structures), the samples display a remarkable combination of a modulus of elasticity of around 25 GPa and nanohardness of 300 kgf mm2, making them some 30–40% stiffer than bone and practically as hard as tooth enamel. At the same time we find that under the same conditions the chars are highly resilient, displaying complete elastic recovery with very little plastic deformation. When cured in the presence of copper(II) acetylacetonate (200 ppm) in dodecylphenol (1% w/v active copper suspension) to form a polycyanurate, compound (2) forms a dense, consolidated structure compared with compound (1) under the same conditions. At high magnification, the presence of a nanoscale, fibrillar structure is observed, accounting for the high resilience.
    Binks Fiona C., Cavalli Gabriel, Henningsen Michael, Howlin Brendan J., Hamerton Ian (2018)Examining the nature of network formation during epoxy polymerisation initiated with ionic liquids, In: Polymer150pp. 318-325 Elsevier
    A commercial diglycidyl ether of bisphenol A monomer (Baxxores™ ER 2200, eew 182 g/mole, DGEBA) is thermally polymerized in the presence of an ionic liquid, 1-ethyl-3-methylimidazolium acetate at a variety of loadings (5–45 wt %). The loss modulus data for cured samples, containing 5 wt % initiator, display at least two thermal transitions and the highest storage modulus occurs in the sample that has been cured for the shortest time at the lowest temperature. Samples that are exposed to higher temperatures (140, 150 °C) yield more heterogenous networks, whereas following exposure to a much shorter/lower temperature cure schedule (80 °C) exhibits a considerably higher damping ability than the other samples, coupled with a lower glass transition temperature. Differential scanning calorimetry reveals that the latter sample achieves a conversion of 95%, while crosslink densities for the DGEBA samples containing 5 wt % and 15 wt % are respectively 9.5 × 10-3 mol. dm-3 and 1.2 × 10-3 mol. dm-3 (when cured to 80 °C) and 2.0 × 10-2 mol. dm-3 and 2.4 × 10-3 mol. dm-3 (when cured to 140 °C).
    Webb SJ, Hanser T, Howlin B, Krause P, Vessey JD (2014)Feature combination networks for the interpretation of statistical machine learning models: application to Ames mutagenicity, In: JOURNAL OF CHEMINFORMATICS6ARTN 8pp. ?-? BIOMED CENTRAL LTD
    Hamerton I, Hay JN, Howlin BJ, Jones JR, Lu SY, Webb GA, Bader MG (1997)Molecular modelling of interactions at the composite interfaces between electrolytically surface-treated carbon fibre and epoxy resin, In: JOURNAL OF MATERIALS CHEMISTRY7(1)pp. 169-174 ROYAL SOC CHEMISTRY
    Binks Fiona C., Cavalli Gabriel, Henningsen Michael, Howlin Brendan J., Hamerton Ian (2019)Examining the Influence of Anion Nucleophilicity on the Polymerisation Initiation Mechanism of Phenyl Glycidyl Ether, In: Polymers11(4) MDPI
    The reaction of phenyl glycidyl ether (PGE) with 1-ethyl-3-methylimidazolium acetate and 1-ethyl-3-methylimidazolium thiocyanate to initiate the polyetherification reaction was examined using thermal and spectral analysis techniques. The influence of the nucleophilicity of the anions on the deprotonation of the 1-ethyl-3-methylimidazolium cation determined the reaction pathway. The thermal degradation of the ionic liquid liberated the acetate ion and led, subsequently, to the deprotonation of the acidic proton in the imidazole ring. Thus, polymerisation of PGE occurred via a carbene intermediate. The more nucleophilic thiocyanate anion was not sufficiently basic to deprotonate the 1-ethyl-3-methylimidazolium cation, and thus proceeded through direct reaction with the PGE, unless the temperature was elevated and a competing carbene mechanism ensued.
    Hassan WAW, Liu J, Howlin BJ, Ishida H, Hamerton I (2016)Examining the influence of bisphenol A on the polymerisation and network properties of an aromatic benzoxazine, In: Polymer88pp. 52-62 ELSEVIER SCI LTD
    A series of reactive blends, comprising a commercial benzoxazine monomer, 2,2-bis(3,4-dihydro-3-phenyl-2H-1,3-benzoxazine)propane, and bisphenol A is prepared and characterized. Thermal analysis and dynamic rheology reveal how the introduction of up to 15 wt % bisphenol A lead to a significant increase in reactivity (the exothermic peak maximum of thermal polymerization is reduced from 245 °C to 215 °C), with a small penalty in glass transition temperature (reduction of 15 K), but similar thermal stability (onset of degradation = 283 °C, char yield = 26%). With higher concentrations of bisphenol A (e.g. 25 wt %), a significantly more reactive blend is produced (exothermic peak maximum = 192 °C), but with a significantly lower thermal stability (onset of degradation = 265 °C, char yield = 22%) and glass transition temperature (128 °C). Attempts to produce a cured plaque containing 35 wt % bisphenol A were unsuccessful, due to brittleness. Molecular modelling is used to replicate successfully the glass transition temperatures (measured using thermal analysis) of a range of copolymers.
    Lockley WJS, Sfyrakis K, Howlin BJ, Jones JR, Wilkinson DJ (2007)Some IT and data processing applications for H, H,C,C and C-labelling, In: Journal of Labelled Compounds and Radiopharmaceuticals50(5-6)pp. 532-534 Wiley
    Thompson Scott, Stone Corinne, Howlin Brendan, Hamerton Ian (2018)Exploring Structure-Property Relationships in Aromatic Polybenzoxazines through Molecular Simulation, In: Polymers MDPI
    A series of commercial difunctional benzoxazine monomers are characterised carefully using thermal and thermo-mechanical techniques before constructing representative polymer networks using molecular simulation techniques. Good agreement is obtained between replicate analyses and for the kinetic parameters obtained from differential scanning calorimetry data (and determined using the methods of Kissinger and Ozawa). Activation energies range from 85-108 kJ/mol (Kissinger) and 89-110 kJ/mol (Ozawa) for the uncatalysed thermal polymerisation reactions, which achieve conversions of 85-97 %. Glass transition temperatures determined from differential scanning calorimetry and dynamic mechanical thermal analysis are coomparable, ranging from BA-a (151 °C, crosslink density 3.6 x 10-3 mol cm-3) containing the bisphenol A moiety to BP-a, based on a phenolphthalein bridge (239-256 °C, crosslink density 5.5-18.4 x 10-3 mol cm-3, depending on formulation). Molecular dynamics simulations of the polybenzoxazines generally agree well with empirical data, indicating that representative networks have been modelled.
    Cunningham ID, Brownhill A, Hamerton I, Howlin BJ (1997)The ene reaction between maleimides and allyl-substituted aromatics, In: TETRAHEDRON53(39)pp. 13473-13494 PERGAMON-ELSEVIER SCIENCE LTD
    Danil de Namor Angela F., Alharthi Salman, Howlin Brendan, Al Hakawati Nawal A selective and easily recyclable dimer based on a calix[4]pyrrole derivative for the removal of mercury(ii) from water, In: RSC Advances10(6)pp. 3060-3071 Royal Society of Chemistry
    A recyclable mercury (II) selective dimer based on a calix[4]pyrrole derivative has been synthesised and characterised by Mass and FT- IR Spectrometry, Scanning Electron Microscopy (SEM) and Energy Dispersive X-ray (EDX). Information regarding the ability of the dimer to interact with metal cations was obtained from FTIR and SEM-EDX analyses. A striking feature of micrographs of the loaded dimer is the change of morphology with the cation. Based on these results, optimal conditions for removing cations from water were assessed under different experimental conditions. Results obtained demonstrate that the removal process is fast. Capacity values and selectivity factors show that the dimer is selective for Hg(II) in one and in multiple component metal solutions relative to other cations. Single-ion transfer Gibbs energies from water to a solvent containing common functionalities to those of the dimer were used to assess the counter-ion effect on the removal process. Agreement is found between these data and energy calculations derived from Molecular Simulation studies. Studies on polluted water in the presence of normal water components in addition to toxic metal cations are reported. Further experimental work on wastewater from the mining industry is in progress.
    Aspin IP, Hamerton I, Howlin BJ, Jones JR, Parker MJ (1998)The computer modelling and the chemical kinetics of the rate of cure of epoxy resins, In: JOCCA-SURFACE COATINGS INTERNATIONAL81(2)pp. 68-+ OIL COLOUR CHEM ASSN
    Hamerton I, Heap K, Howlin BJ, McNamara LT (2013)Effects of thermal history on the polymerisation mechanism and network development in aromatic polybenzoxazines, In: Reactive and Functional Polymers
    The effect of heating rate (2, 8 and 15 K min-1) during the initial stages of cure of 2,2-bis(3,4-dihydro-3- phenyl-2H-1,3-benzoxazine)propane is examined. The rate of heating has a marked effect on the observed modulus, measured by DMTA, with the higher heating rate giving rise to an increase in storage modulus of ca. 1000 MPa, although this is not accompanied by an increase in glass transition temperature. The thermal stability of the resulting polybenzoxazines also differs with the slower heating rate giving rise to less thermally stable structures. Data obtained from Raman spectroscopy (when combined with principal components analysis) suggest subtle changes in the mechanism during the early stages of reaction associated C–N–C and C–O moieties, some of which persist following a higher temperature postcure step leading to a crosslinked network with higher aliphatic character.
    Hamerton I, Hay JN, Howlin BJ, Jones JR, Lu SY, Webb GA, Bader MG, Brown AM, Watts JF (1997)ToF SIMS and XPS studies of carbon fiber surface during electrolytic oxidation in O-17/O-18 enriched aqueous electrolytes, In: CHEMISTRY OF MATERIALS9(9)pp. 1972-1977 AMER CHEMICAL SOC
    Malkhasian Aramice Y.S., Howlin Brendan J. (2019)Automated drug design of kinase inhibitors to treat Chronic Myeloid Leukemia, In: Journal of Molecular Graphics and Modelling91pp. 52-60 Elsevier
    Medicinal chemistry has in the past been dominated by learned insights from experienced organic chemists. However, with the advent of computer based methods, computer aided drug design has become prominent. We have compared here the ability of locally sourced expert medicinal chemists to purely automated methods and found that the automated method produces a better potential candidate drug than the expert input. The example chosen is based on inhibitors to Abl-kinase and the successful anti-leukaemic drug imatinib. The proposed molecule is a simple modification of nilotinib and has a docking energy of 4.2 kJ/mol better than the best intuitive molecule.
    Chaaban Jinane, Hakawati Nawal Al, Howlin Brendan, Bance-Soualhi Rachida, Danil de Namor Angela An Asymmetric N-Rim Partially Substituted Calix[4]pyrrole: Its Affinity for Ag(I) and its Destruction by Hg(II), In: Arabian Journal of Chemistry Elsevier
    A new partially substituted calix[4] pyrrole derivative obtained by the introduction of three thioamide functionalities in the N-rim has been synthesised and fully characterised by 1H,13C, HSQC, ROESY NMR and mass spectroscopy. Computer modelling suggested an alternate conformation which was confirmed through ROESY 1H NMR. The receptor interacts only with the silver cation as shown by 1H NMR. The strength of interaction is quantitatively assessed by titration calorimetry. N-rim modification eliminates the possibility of interaction with anions. Unlike calix[4] pyrrole derivatives obtained by the introduction of functionalities through the meso-position, addition of Hg(NO3)2 leads to the degeneration of the receptor as demonstrated by 1H NMR, FTIR and XPS analyses. This is for the first time reported. Molecular simulation studies show significant strain in the mercury bound ligand in bonds, angles, torsions leading to the destruction of the receptor. Given the negative environmental impact produced by the availability of silver ions in aquatic organisms, the fundamental studies indicate that this receptor offers potential applications for monitoring silver (ion selective electrode) or indeed as a decontaminating material for removing silver ions from water.
    Molecular Operating Environment (MOE) software has great potential when combined with the Quantitative Structure-Property Relationship (QSPR) approach, and was proven to be useful to make good prediction models for series of polybenzoxazines [1–3]. However, the effect of heterogeneities in the crosslinked network to the prediction accuracy is yet to be tested. It was found that polybenzoxazines with polymerisable functional group (e.g. acetylene-based benzoxazines) form up to 40% higher char yield compared to their analogue polybenzoxazines due to the contribution of the polymerisable functional group (e.g. ethynyl triple bond) in the cross-linked network. In order to investigate the effect of the inconsistent cross-linking network, a data set consisting of thirty-three benzoxazines containing various structures of benzoxazines was subdivided into two smaller data sets based on their functional group, either benzoxazines with polymerisable functional group (acetylene-based benzoxazines set (Ace-M)) or non-polymerisable functional group (aniline-based benzoxazines (Ani-M)). Char yield predictions for the polybenzoxazines for these data sets (Ace-M and Ani-M) were compared with the larger thirty-three polybenzoxazines data set (GM) to investigate the effect of the inconsistency in crosslink network on the quality of prediction afforded by the model. Prediction performed by Ace-M and Ani-M were found to be more accurate when compared with the GM with total prediction error of 3.15% from both models compared to the GM (4.81%). Ace-M and Ani-M are each better at predicting the char yields of similar polybenzoxazines (i.e. one model is specific for a polymerisable functional group; the other for non-polymerisable functional group), but GM is more practical as it has greater ‘general’ utility and is applicable to numerous structures. The error shown by GM is considerably small and therefore it is still a good option for prediction and should not be underestimated.
    Hamerton I, Howlin BJ, Yeung S-YC (2013)Designing thermoplastic oligomers with programmed degradation mechanisms using a combined empirical and simulation approach, In: Polymer Degradation and Stability98(4)pp. 829-838
    Molecular simulation is used to probe the structure property relationships displayed by polysulphone (PS) and polyethersulphone (PES) and reproduces closely the temperatures at which thermal degradation occurs (and the glass transition temperatures). Both data sets agree well with those obtained empirically using TGA. The examination of the thermal and thermo-oxidative stability of thermoplastic oligomers (M = 5454-33,866 g mol , PDI 1.33-1.82) based on PS, PES, polyetherimide (PEI) and poly(amide-imide) (PAI), is reported. TGA reveals the least thermally stable polymer is PES (T = ∼250 °C), while PAI (T = ∼350 °C) is the highest: the materials usually display two-step decomposition patterns: scission of bridging group and degradation of backbone structure. A possible mechanism for the degradation of a PAI is proposed on the basis of the empirical and simulation data. This work provides a general method for the prediction of the thermal stability of oligomeric modifiers (and high molecular weight polymers). © 2013 Elsevier Ltd. All rights reserved.
    de Namor AFD, Al Nuaim M, Villanueva Salas JA, Bryant S, Howlin Brendan (2017)A calix[4]arene derivative and its selective interaction with drugs (clofibric acid, diclofenac and aspirin), In: European Journal of Pharmaceutical Sciences100pp. 1-8 Elvesier
    The synthesis and characterisation of a partially substituted calix[4]arene, namely, 5,11,17,23- tetra-tert-butyl,25,27-bis[aminoethoxy] 26,28-dihydroxycalix[4]arene are reported. Its interaction with commonly used pharmaceuticals (clofibric acid, diclofenac and aspirin) was investigated by spectroscopic (1 H NMR and UV), electrochemical (conductance measurements) and thermal (titration calorimetry) techniques. It is concluded on the basis of the experimental work and molecular simulation studies that the receptor interacts selectively with these drugs. Preliminary studies on the selective extraction of these pharmaceuticals from water by the calix receptor are reported and the potential for a carrier mediated sensor based on this ligand for ‘on site’ monitoring of pharmaceuticals is discussed.
    Hamerton Ian, Howlin Brendan, Tilbrook DA (2014)Examining the preparation and characterization of coatings based on linear aromatic terpoly(methoxy-cyanurate-thiocyanurate)s, In: Polymer International63(1)pp. 60-71
    Two series of terpoly(methoxy-cyanurate-thiocyanurate)s based on thiodiphenol and dithiodiphenyl sulphide, and dihydroxydiphenylether and dithiodiphenyl ether, are prepared in good yield and purity, and fully characterised. Most of the resulting polymers, formed at room temperature using phase transfer catalysis, can be cast into films with good resilience and thermal stability. Two series of terpoly(methoxy-cyanurate-thiocyanurate)s based on thiodiphenol and dithiodiphenyl sulfide and on dihydroxydiphenyl ether and dithiodiphenyl ether, were prepared in good yield and purity and fully characterized. Most of the resulting polymers, formed at room temperature using phase transfer catalysis, can be cast into films with good resilience and thermal stability (some examples suffer practically no mass loss when held isothermally at 190 °C and only display appreciable losses when held continuously at 225 °C). Char yields of 53%-61% are achieved in nitrogen depending on backbone structure. Some problems were encountered with solubility, particularly with copolymers, which limited molecular weight analysis, but values of Mn = 8000-13000 g mol-1 were obtained for the polymers based on thiodiphenol and dithiodiphenyl sulfide, and Mn = 5000-13000 g mol-1 for the polymers based on dihydroxydiphenyl ether and dithiodiphenyl ether. DSC reveals polymerization exotherms with maxima at 184-207 °C (ΔHp = 43-59 kJ mol-1), which are believed to be due to isomerization of the cyanurate to the isocyanurate (activation energies span 159-195 kJ mol-1). Molecular simulation shows that diphenylether and diphenylsulfide display very similar conformational energy surfaces and would therefore be expected to adopt similar conformations, but the diphenylsulfide offers less resistance to deformations that increase the proximity of the two phenyl rings and results in more resilient films. © 2013 Society of Chemical Industry.
    Barton JM, Chaplin A, Hamerton I, Howlin BJ (1999)A new synthetic route for the preparation of alkenyl functionalized aryl cyanate ester monomers, In: POLYMER40(19)pp. 5421-5427 ELSEVIER SCI LTD
    Hamerton I, Hay JN, Howlin BJ, Jones JR, Lu SY, Webb GA, Bader MG (1997)Probing the cure of C-13 labelled bisphenol A dicyanate ester in carbon fibre reinforced composites using solid state C-13 NMR, SEM and FTIR, In: POLYMER BULLETIN38(4)pp. 433-438 SPRINGER VERLAG
    McNamara LT, Howlin BJ, Smith PA, Cross P, Ward S, Hamerton I (2013)Examining the Initiation of the Polymerization Mechanism and Network Development in Aromatic Polybenzoxazines., In: Macromolecules46(13)pp. 5117-5132
    Three bis-benzoxazine monomers based on the aniline derivatives of bisphenol A (BA-a), bisphenol F (BF-a), and 3,3'-thiodiphenol (BT-a) are examined using a variety of spectroscopic, chromatographic, and thermomechanical techniques. The effect on the polymerization of the monomers is compared using two common compounds, 3,3'-thiodiphenol (TDP) and 3,3'-thiodipropionic acid (TDA), at a variety of loadings. It is found that the diacid has a greater effect on reducing the onset of polymerization and increasing cross-link density and T g for a given benzoxazine. However, the addition of >5 wt % of the diacid had a detrimental effect on the cross-link density, T g, and thermal stability of the polymer. The kinetics of the polymerization of BA-a were found to be well described using an autocatalytic model for which values of n = 1.64 and m = 2.31 were obtained for the early and later stages of reaction (activation energy = 81 kJ/mol). Following recrystallization the same monomer yielded values n = 1.89, m = 0.89, and E a = 94 kJ/mol (confirming the influence of higher oligomers on reactivity). The choice of additive (in particular the magnitude of its pK a) appears to influence the nature of the network formation from a linear toward a more clusterlike growth mechanism.
    Baggott A, Bass JR, Hall SA, Hamerton I, Howlin BJ, Mooring L, Sparks D (2014)At the limits of simulation: A new method to predict thermal degradation behavior in cyanate esters and nanocomposites using molecular dynamics simulation, In: Macromolecular Theory and Simulations23(6)pp. 369-372
    Despite their inability to model bond breaking molecular dynamics simulations are shown to predict thermal degradation temperatures of polycyanurate (cyanate ester) homopolymers and nanocomposites in very close agreement with experimental data. Simulated polymer density, used to predict T also shows a reduction within the same temperature range as experimental values for the thermal degradation. © 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.
    The preparation of a series of oligomeric engineering thermoplastics (PS, PES, PEI and PAI) is reported. H and C NMR and FT-IR/ATR spectroscopic techniques are combined to determine the chemical structure of the synthetic polymers, which are produced in good yield and purity. GPC measurements show the weight average molecular weight (M ) of the synthesised thermoplastics fall in the range 5454-33,866 g mol and display polydispersity indices in the range 1.33-1.82. Glass transition temperatures (T) values measured by DSC, occur between 107 and 257 °C and fall in the pattern PS < PES < PEI < PAI. Molecular simulation is used to probe the structure property relationships displayed by PS and PES and reproduces the elastic properties of PS well within the range of the literature; while the values of PES are less well reproduced. The simulated T values of both oligomers agree well with those obtained empirically using DSC. © 2014 Elsevier Ltd. All rights reserved.
    Binks F, Cavalli Gabriel, Henningsen M, Howlin Brendan, Hamerton Ian (2018)Investigating the mechanism through which ionic liquids initiate the polymerisation of epoxy resins, In: Polymer139pp. 163-176 Elsevier
    The mechanism of reaction between 1-ethyl-3-methylimidazolium acetate and the difunctional diglycidyl ether of bisphenol A (DGEBA) is explored using thermal and spectroscopic methods. Investigation of the 1,3-dialkylimidazolium based ionic liquids comprising the common cation (1-ethyl-3-methylimidazolium) and different anions (acetate, diethyl phosphate, dicyanamide or thiocyanate) via thermogravimetric analysis revealed 1-ethyl-3-methylimidazolium acetate to be the least thermally stable, both in air and nitrogen, and 1-ethyl-3-methylimidazolium dicyanamide to be the most thermally stable. Dynamic differential scanning calorimetry reveals the formulations comprising DGEBA and ionic liquid where it was revealed that the lowest and highest temperature for the onset of reaction were observed for formulations with 1-ethyl-3-methylimidazolium acetate and 1-ethyl-3-methylimidazolium dicyanamide respectively. 1-Ethyl-3-methylimidazolium acetate was shown, via nuclear magnetic resonance (NMR) spectroscopy and residual gas analysis, to degrade at 150 °C to yield dealkylated products including methyl acetate and ethyl acetate as well as 1-methylimidazole and 1- ethylimidazole. The dealkylated imidazole ring is proposed as a route for initiation of the epoxy ring. Adduct formation between 1-ethyl-3-methylimidazoloium acetate and benzaldehyde at room temperature was observed leading to the proposal of the generation of a carbene species as a route for initiation of the epoxy ring in formulations with the acetate anion. NMR analysis of formulations comprising 1-ethyl-3-methylimidazolium thiocyanate and epoxy are believed, at room temperature, to initiate via reaction of the thiocyanate anion with the epoxy ring.At elevated temperatures, it is proposed that a second, competing reaction, involving deprotonation of the imidazolium ring, also becomes active. The three proposed reaction pathways, namely the carbene route, the imidazole route and the counter-ion route, are all proposed to occur when an ionic liquid is used to initiate an epoxy resin.
    The p22phox is a key component of the cytochrome b558 of the NADPH oxidase (Nox), which by generating reactive oxygen species (ROS) is involved in the pathogenesis of cardiovascular disease. A p22phox polymorphism (C242T) has been found to reduce oxidative stress in the cardiovascular system and is negatively associated with the incidence of coronary heart disease (CHD). However, the mechanism involved remains unknown. In this study we used computer molecular modelling and bioinformatics to investigate the potential effect of C242T polymorphism on the 3-D protein structure of the p22phox. Based on the published sequence data of p22phox and the principle of regulated prediction algorithms, we found that p22phox consists of two domains: an N-terminal transmembrane domain (124 a.a.) and a C-terminal cytoplasmic domain (71 a.a.). In its most stable form, it has three e26 Heart September 2010 Vol 96 No 17 BSCR Spring 2010 Meeting Abstracts Downloaded from heart.bmj.com on January 31, 2012 - Published by group.bmj.com transmembrane helices leading to an extracellular N-terminus and an extracellular loop between helices two and three. The C242T polymorphism causes a change of His72 to Tyr72. This change results in significant morphological changes of the extracellular loop of the p22phox, which is in the putative interactive region of the p22phox with the catalytic subunit (Nox2). This may interfere with their association, and subsequently result in a reduced cytochrome b function and a reduced ROS production by NADPH oxidase. These results give us insight into the molecular mechanism of this polymorphism in reducing vascular oxidative stress and may explain how this polymorphism is linked with reduced incidence of CHD.
    Howlin BJ, Meijles DN, Li JM (2011)Computer molecular modelling of the p22phox protein structural changes linked to C242T polymorphism, In: British Journal of Clinical Pharmacology71(6)pp. 971-994
    The p22phox is a key component of the cytochrome b558 of the NADPH oxidase (Nox), which by generating reactive oxygen species (ROS) is involved in the pathogenesis of cardiovascular disease. A p22phox polymorphism (C242T) has been found to reduce oxidative stress in the cardiovascular system and is negatively associated with the incidence of coronary heart disease (CHD). However, the mechanism involved remains unknown. In this study we used computer molecular modelling and bioinformatics to investigate the potential effect of C242T polymorphism on the 3-D protein structure of the p22phox. Based on the published sequence data of p22phox and the principle of regulated prediction algorithms, we found that p22phox consists of two domains: an N-terminal transmembrane domain (124 a.a.) and a C-terminal cytoplasmic domain (71 a.a.). In its most stable form, it has three transmembrane helices leading to an extracellular N-terminus and an extracellular loop between helices two and three. The C242T polymorphism causes a change of His72 to Tyr72. This change results in signifi cant morphological changes of the extracellular loop of the p22phox which is in the putative interactive region of the p22phox with the catalytic subunit (Nox2). This may interfere with their association, and subsequently result in a reduced cytochrome b function and a reduced ROS production by NADPH oxidase. These results give us insight into the molecular mechanism of this polymorphism in reducing vascular oxidative stress and may explain how this polymorphism is linked with reduced incidence of CHD.
    Hamerton I, Hay JN, Howlin BJ, Jones J R , Lu S Y , Webb G A , Bader M G (1996)Surface analysis of carbon fibres using a tritium-hydrogen exchange procedure, In: 7th European Conference on Composite Materials 'Realising their Commercial Potential'2(1)pp. 435-439
    Danil de Namor Angela, Abou Hamdan W, Webb O, Bance- Soualhi R, Howlin Brendan, Al Hakawati N (2018)Calix(4)arene urea derivatives: the pathway from fundamental studies to the selective removal of fluorides and phosphates from water, In: Journal of Hazardous Materials Elsevier
    Wet processes of phosphoric acids produce untreated wastewater containing large amounts of fluoride leading to serious environmental problems. This paper reports fundamental studies on two lower rim functionalised calix[4]arene based receptors namely 5, 11, 17, 23 tetra-tert-butyl, 25, 27 bis [diethylphenylurea]ethoxy, 26, 28 dihydroxycalix[4]arene, 1 and 25, 27 bis[diethylphenylurea]ethoxy, 26, 28 dihydroxycalix[4]arene, 2 and their ionic interactions. It is shown that these receptors only interact with fluoride and phosphate in acetonitrile. Their receptive properties are higher for phosphate (2:1 anion:receptor complex) relative to fluoride (1:1 complex).However thermodynamics shows that these receptors are more selective for fluoride relative to phosphate in the formation of the 1:1 complex. The pathway from fundamental studies to the use of these receptors for removing these anions from water has been tested. The receptive properties of 1 for phosphate are held when the extraction involves aqueous solutions containing individual ions. However in mixtures containing both anions, the kinetics of the process and the selectivity of 1 for fluoride predominate and as a result, fluoride is better extracted than phosphate. The counter-ion effect on the removal process is assessed from Molecular Simulation studies. The removal of fluoride from phosphate is discussed taking into account existing technologies.
    Hamerton I, Hay JN, Howlin BJ, Jones JR, Lu SY (1998)Molecular modelling of interactions at the composite interface between surface-treated carbon fibre and polymer matrices: the influence of surface functional groups, In: JOURNAL OF MATERIALS CHEMISTRY8(6)pp. 1333-1337 ROYAL SOC CHEMISTRY
    Deazle AS, Hamerton I, Heald CR, Howlin BJ (1996)Molecular modelling of high performance polymers, In: POLYMER INTERNATIONAL41(2)pp. 151-157 JOHN WILEY & SONS LTD
    Howlin BJ (2012)Predicting Glass Transition Temperatures of Polyarylethersulphones Using QSPR Methods, In: Howlin BJ, Hamerton I, Kamyszek G (eds.), PLoS One7(6)e38424
    The technique of Quantitative Structure Property Relationships has been applied to the glass transition temperatures of polyarylethersulphones. A general equation is reported that calculates the glass transition temperatures with acceptable accuracy (correlation coefficients of between 90–67%, indicating an error of 10–30% with regard to experimentally determined values) for a series of 42 reported polyarylethersulphones. This method is quite simple in assumption and relies on a relatively small number of parameters associated with the structural unit of the polymer: the number of rotatable bonds, the dipole moment, the heat of formation, the HOMO eigenvalue, the molar mass and molar volume. For smaller subsets of the main group (based on families of derivatives containing different substituents) the model can be simplified further to an equation that uses the volume of the substituents as the principal variable.
    Meijles DN, Fan LM, Ghazaly MM, Howlin Brendan, Krönke M, Brooks G, Li J-M (2016)p22phox C242T SNP Inhibits Inflammatory Oxidative Damage to Endothelial Cells and Vessels, In: Circulationpp. 2391-2403 American Heart Association
    Background—T NADPH oxidase, by generating reactive oxygen species, is involved in the pathophysiology of many cardiovascular diseases and represents a therapeutic target for the development of novel drugs. A single-nucleotide polymorphism (SNP) C242T of the p22phox subunit of NADPH oxidase has been reported to be negatively associated with coronary heart disease (CHD) and may predict disease prevalence. However, the underlying mechanisms remain unknown. Methods and Results—Using computer molecular modelling we discovered that C242T SNP causes significant structural changes in the extracellular loop of p22phox and reduces its interaction stability with Nox2 subunit. Gene transfection of human pulmonary microvascular endothelial cells showed that C242T p22phox reduced significantly Nox2 expression but had no significant effect on basal endothelial O2 .- production or the expression of Nox1 and Nox4. When cells were stimulated with TNFĮ (or high glucose), C242T p22phox inhibited significantly TNFĮ- induced Nox2 maturation, O2 .- production, MAPK and NFțB activation and inflammation (all p<0.05). These C242T effects were further confirmed using p22phox shRNA engineered HeLa cells and Nox2-/- coronary microvascular endothelial cells. Clinical significance was investigated using saphenous vein segments from non CHD subjects after phlebectomies. TT (C242T) allele was common (prevalence of ~22%) and compared to CC, veins bearing TT allele had significantly lower levels of Nox2 expression and O2 .- generation in response to high glucose challenge. Conclusions—C242T SNP causes p22phox structural changes that inhibit endothelial Nox2 activation and oxidative response to TNFĮ or high glucose stimulation. C242T SNP may represent a natural protective mechanism against inflammatory cardiovascular diseases.
    Bone Matthew A., Macquart Terence, Hamerton Ian, Howlin Brendan J. (2020)A Novel Approach to Atomistic Molecular Dynamics Simulation of Phenolic Resins Using Symthons, In: Polymers12(4) MDPI AG
    Materials science is beginning to adopt computational simulation to eliminate laboratory trial and error campaigns—much like the pharmaceutical industry of 40 years ago. To further computational materials discovery, new methodology must be developed that enables rapid and accurate testing on accessible computational hardware. To this end, the authors utilise a novel methodology concept of intermediate molecules as a starting point, for which they propose the term ‘symthon’a rather than conventional monomers. The use of symthons eliminates the initial monomer bonding phase, reducing the number of iterations required in the simulation, thereby reducing the runtime. A novel approach to molecular dynamics, with an NVT (Canonical) ensemble and variable unit cell geometry, was used to generate structures with differing physical and thermal properties. Additional script methods were designed and tested, which enabled a high degree of cure in all sampled structures. This simulation has been trialled on large-scale atomistic models of phenolic resins, based on a range of stoichiometric ratios of formaldehyde and phenol. Density and glass transition temperature values were produced, and found to be in good agreement with empirical data and other simulated values in the literature. The runtime of the simulation was a key consideration in script design; cured models can be produced in under 24 h on modest hardware. The use of symthons has been shown as a viable methodology to reduce simulation runtime whilst generating accurate models.
    Mas Claret Eduard, Al Yahyaei Balqees, Chu Shuyu, Elliott Ruan M., Imperato Manuel, Lopez Arnaud, Meira Lisiane B., Howlin Brendan J., Whelligan Daniel K. (2020)An aza-nucleoside, fragment-like inhibitor of the DNA repair enzyme alkyladenine glycosylase (AAG), In: Bioorganic & Medicinal Chemistry28(11)115507 Elsevier
    The DNA repair enzyme AAG has been shown in mice to promote tissue necrosis in response to ischaemic reperfusion or treatment with alkylating agents. A chemical probe inhibitor is required for investigations of the biological mechanism causing this phenomenon and as a lead for drugs that are potentially protective against tissue damage from organ failure and transplantation, and alkylative chemotherapy. Herein, we describe the rationale behind the choice of arylmethylpyrrolidines as appropriate aza-nucleoside mimics for an inhibitor followed by their synthesis and the first use of a microplate-based assay for quantification of their inhibition of AAG. We finally report the discovery of an imidazol-4-ylmethylpyrrolidine as a fragment-sized, weak inhibitor of AAG.
    Hamerton I, Chaplin A, Howlin BJ (1994)Thermal properties of functionalised aryl cyanate ester oligomers, In: Proceedings of the American Chemical Society Division of PMSE71pp. 811-812
    Hamerton I, Howlin BJ, Klewpatinond P, Takeda S (2002)Conformational studies of polycyanurates: a study of internal stress versus molecular structure, In: POLYMER43(17)PII S0032-pp. 4599-4604 ELSEVIER SCI LTD
    The construction of molecular models of crosslinked polymers is an area of some difficulty and considerable interest. We report here a new method of constructing these models and validate the method by modelling three epoxy systems based on the epoxy monomers bisphenol F diglycidyl ether (BFDGE) and triglycidyl-p-amino phenol (TGAP) with the curing agent diamino diphenyl sulphone (DDS). The main emphasis of the work concerns the improvement of the techniques for the molecular simulation of these epoxies and specific attention is paid towards model construction techniques, including automated model building and prediction of glass transition temperatures (T(g)). Typical models comprise some 4200-4600 atoms (ca. 120-130 monomers). In a parallel empirical study, these systems have been cast, cured and analysed by dynamic mechanical thermal analysis (DMTA) to measure T(g). Results for the three epoxy systems yield good agreement with experimental T(g) ranges of 200-220°C, 270-285°C and 285-290°C with corresponding simulated ranges of 210-230°C, 250-300°C, and 250-300°C respectively.
    Lisgarten DR, Palmer RA, Lobley CMC, Naylor CE, Chowdhry BZ, Al‑Kurdi ZI, Badwan AA, Howlin Brendan, Gibbons NCJ, Saldanha JW, Lisgarten JN, Basak AK (2017)Ultra‑high resolution X‑ray structures of two forms of human recombinant insulin at 100 K, In: Chemistry Central Journal11(73) Springer
    The crystal structure of a commercially available form of human recombinant (HR) insulin, Insugen (I), used in the treatment of diabetes has been determined to 0.92 Å resolution using low temperature, 100 K, synchrotron X-ray data collected at 16,000 keV (λ = 0.77 Å). Refinement carried out with anisotropic displacement parameters, removal of main-chain stereochemical restraints, inclusion of H atoms in calculated positions, and 220 water molecules, converged to a final value of R = 0.1112 and Rfree = 0.1466. The structure includes what is thought to be an ordered propanol molecule (POL) only in chain D(4) and a solvated acetate molecule (ACT) coordinated to the Zn atom only in chain B(2). Possible origins and consequences of the propanol and acetate molecules are discussed. Three types of amino acid representation in the electron density are examined in detail: (i) sharp with very clearly resolved features; (ii) well resolved but clearly divided into two conformations which are well behaved in the refinement, both having high quality geometry; (iii) poor density and difficult or impossible to model. An example of type (ii) is observed for the intra-chain disulphide bridge in chain C(3) between Sγ6–Sγ11 which has two clear conformations with relative refined occupancies of 0.8 and 0.2, respectively. In contrast the corresponding S–S bridge in chain A(1) shows one clearly defined conformation. A molecular dynamics study has provided a rational explanation of this difference between chains A and C. More generally, differences in the electron density features between corresponding residues in chains A and C and chains B and D is a common observation in the Insugen (I) structure and these effects are discussed in detail. The crystal structure, also at 0.92 Å and 100 K, of a second commercially available form of human recombinant insulin, Intergen (II), deposited in the Protein Data Bank as 3W7Y which remains otherwise unpublished is compared here with the Insugen (I) structure. In the Intergen (II) structure there is no solvated propanol or acetate molecule. The electron density of Intergen (II), however, does also exhibit the three types of amino acid representations as in Insugen (I). These effects do not necessarily correspond between chains A and C or chains B and D in Intergen (II), or between corresponding residues in Insugen (I). The results of this comparison are reported.
    Thompson Scott, Howlin Brendan J., Stone Corinne A., Hamerton Ian (2018)Exploring the thermal degradation mechanisms of some polybenzoxazines under ballistic heating conditions in helium and air, In: Polymer Degradation and Stability156pp. 180-192 Elsevier
    The degradation behaviour of five polybenzoxazines (PBZs) is studied using pyrolysis-GC/MS. Upon heating to 800 °C in helium the PBZs generate a variety of similar pyrolysis products including aniline (the major product in all cases), substituted phenols, acridine, and 9-vinylcarbazole. During the initial stages of heating (200–300 °C) aniline is the dominant pyrolysis product; from 350 °C onwards substituted phenols are released, particularly 2-methylphenol and 2,6-dimethyl phenol. The same major species are produced on heating in air, but in addition isocyanatobenzene is observed which results from the oxidation of Mannich bridges, along with a number of sulphurous species from the monomer containing a thioether bridge. This suggests that sulphur is more likely to be retained in the char in a helium atmosphere, but takes part in oxidative reactions to form pyrolysis fragments in air. During the ramped temperature cycles in both air and helium atmospheres the release of aniline was observed to rise, fall and then to rise again. This may be due a combination of the very high heating rate, poor thermal conduction of the polymer and the availability of the Mannich bridges to undergo breakdown.
    Howlin BJ, Malkhasian AYS (2015)Docking and DFT Studies on Ligand Binding to Quercetin 2,3-Dioxygenase, In: Journal of Biomolecular Structure and Dynamics Taylor & Francis
    Simple molecular docking calculations on quercetin, kojic acid and diethylcarbamatodithoic acid using the software package MOE are shown to be close to the geometries reported in the X-ray crystal structures of the protein co-crystallised with the respective ligands. Furthermore DFT optimization of the docked conformations is shown to reproduce the essential features of previous studies on quercetin, showing that docking can be used to provide good starting structures for mechanistic study. The flavone ligand, lacking the hydroxyl group of the quercetin is shown by docking to be unable to approach closely the copper atom, indicating the necessity of the presence of the hydroxyl group and providing a prediction of the likely binding environment of this ligand.
    Crawford AO, Cavalli Gabriel, Howlin Brendan, Hammerton I (2016)Improving the Hydrolytic Stability of Aryl Cyanate Esters by Examining the Effects of Extreme Environments on Polycyanurate Copolymers, In: REACTIVE & FUNCTIONAL POLYMERS109(Dec)pp. 104-111 ELSEVIER SCIENCE BV
    Three cyanate ester monomer or oligomer species: 2,2-bis(4-cyanatophenyl)propane 1, 1-1-bis(4-dicyanatophenyl)ethane 2, and the oligomeric phenolic cyanate 3, are blended in various ratios to form binary mixtures (18 in total), formulated with copper(II) acetylacetonate (200 ppm) in dodecylphenol (1 % w/v active copper suspension) and cured (3 Kmin-1 to 150 °C + 1 hour; 3 Kmin-1 to 200 °C + 3 hours) followed by a post cure (3 Kmin-1 to 260 °C + 1 hour). Cured copolymers were exposed to environments of elevated relative humidity (75 % RH) and parallel immersion testing in H2O, H2SO4 (10 %) and NaOH (10 %) at 25 °C for a period of up to 2 years and accelerated ageing in boiling water (14 days). Periodic measurements are made of moisture gain along with infrared spectra and compared with cured homopolymers. Changes in mass are recorded periodically throughout exposure, prior to destructive thermo-mechanical analyses. Dynamic mechanical thermal analysis data comparing neat and exposed blends demonstrate the detrimental effect of moisture ingress whilst data from thermogravimetric analysis demonstrate no change in degradation onset between neat and exposed materials. An optimised blend of 1:1 of monomer units 1 and 2 was found to absorb less moisture than blends of different stoichiometry or between other respective monomeric units, consequently limiting the deleterious effect of moisture ingress.
    Bi- and tri-aromatic compounds of the formula (I) wherein R1 to RIO and X are as defined, are Nox2 inhibitors that are useful as medicaments for the treatment of a disease or condition selected from: cardiovascular diseases, respiratory diseases, inflammatory diseases, cancers, ageing and age related disorders, kidney diseases, neurodegenerative diseases, diabetes and conditions associated with diabetes. The compounds, their preparation and pharmaceutical compositions comprising them are disclosed.
    The thesis deals with an investigation on calix[4] based receptors, namely calix[4]pyrrole and calix[4]arene derivatives and their use in different environmental and medical applications. The work conducted throughout this study addresses three main topics. i) The first topic discusses the ability of meso-octamethylcalix[4]pyrrole, CP to remove arsenite and arsenate anions from natural waters. Sites and nature of interaction of the receptor with the investigated anions using 1H NMR and molecular simulation studies are identified. Optimal environmental conditions led to maximal uptake of the investigated anions by CP, uptake capacity of the receptor and the kinetics of the process for the investigated anions are determined. CP was used to remove arsenic species from real water samples and its efficiency is demonstrated. ii) The second topic is related to the synthesis and structural characterisation of two thioacetamide calix[4]pyrrole derivatives namely, meso-tetramethyl-tetrakis[(diethylthiocarbamoyl)phenoxy] calix[4]-pyrrole, CPIII, and 21,22,23-tri(N,N-diethyl-thioacetamide)octamethylcalix[4]pyrrole, CPII. The selectivity of these receptors towards mercury (II) cation relative to other cations was investigated through 1H NMR, conductance and thermodynamic complexation measurements in acetonitrile. Based on the outcome of fundamental studies, a CPIII based mercury (II) ion selective electrode was constructed and its optimal working conditions were determined. The role of thermodynamics in the selection of the electrode is discussed in relation to selectivity coefficients of the primary ion relative to other cations. Another work on Hg(II)-ion selective electrode using partially substituted calix[4]arene derivative with thioacetamide moiety, namely, 5,11,17,23-tetra-tert-butyl[25,27bis(diethylthiocarbamoyl)oxy]calix[4]arene, CAII, is discussed. iii) The third topic discusses the bioactivity of CAII against human adenocarcinoma Caco-2 cell line (colorectal cancer). Cytotoxicity and cell death mechanism using flow cytometry, 1H NMR and Raman Confocal spectroscopic techniques are reported. Conclusions and suggestions for further investigations are included.
    Dietary obesity, insulin resistance (IR) and type-2 diabetes mellitus (T2DM) are the major risk factors for the development of cardiovascular diseases (CVDs). Recent studies have shown that reactive oxygen species (ROS)-generating enzyme, NADPH oxidase (Nox), causes endothelial oxidative stress and is a prominent feature underlying vascular abnormalities of patients with IR and T2DM. The Nox2 has been found to be constitutively expressed in endothelial cells (ECs) and can be activated in diabetic condition. It is known that short chain fatty acids (SCFAs) i.e. acetate (NaA), propionate (NaP) and butyrate (NaB) are byproducts of bacterial fermentation of resistant starch in the colon and exert beneficial effects in improving insulin sensitivity, colonic cancer and reducing CVDs development. However, the effects of SCFAs in reducing endothelial oxidative stress and dysfunction are still unknown. In this project, I have investigated high glucose (0-30 mM for 24 h)-induced ROS production, Nox2 activation levels, endothelial function (in terms of cell cycle regulation and capillary formation on matrigels), G protein-coupled receptor 43 levels (GPCR43, a SCFAs receptor) and the potential therapeutic effect of SCFAs in diabetes using cultured human pulmonary microvascular ECs (HPMEC). My results have shown that high glucose (30 mM) induced endothelial Nox2 activation and ROS production, which causes ECs cycle deregulation, cell apoptosis and capillaries broken down on matrigels. High glucose increases the levels of cyclin D, E, A and B, p53 and p16 protein expression which lead to abnormal ECs cell cycle regulation, apoptosis and cell death. Adding NaA (5 mM, GPCR43 agonist) into the culture medium can increase the GPCR43 expression on ECs under high glucose conditions, inhibits Nox2-induced ROS production and activation as well as improves ECs function. In short, high glucose causes endothelial oxidative stress and dysfunction via Nox2-induced ROS production. NaA activated GPCR43 inhibits Nox2 activation which reversed high glucose-induced damage on the ECs. SCFAs such as NaA may have the therapeutic potential for treating oxidative stress-related vascular complications in patients with IR, obesity and diabetes.
    Structure Activity Relationship (SAR) modelling capitalises on techniques developed within the computer science community, particularly in the fields of machine learning and data mining. These machine learning approaches are often developed for the optimisation of model accuracy which can come at the expense of the interpretation of the prediction. Highly predictive models should be the goal of any modeller, however, the intended users of the model and all factors relating to usage of the model should be considered. One such aspect is the clarity, understanding and explanation for the prediction. In some cases black box models which do not provide an interpretation can be disregarded regardless of their predictive accuracy. In this thesis the problem of model interpretation has been tackled in the context of models to predict toxicity of drug like molecules. Firstly a novel algorithm has been developed for the interpretation of binary classification models where the endpoint meets defined criteria: activity is caused by the presence of a feature and inactivity by the lack of an activating feature or the deactivation of all such activating features. This algorithm has been shown to provide a meaningful interpretation of the model’s cause(s) of both active and inactive predictions for two toxicological endpoints: mutagenicity and skin irritation. The algorithm shows benefits over other interpretation algorithms in its ability to not only identify the causes of activity mapped to fragments and physicochemical descriptors but also in its ability to account for combinatorial effects of the descriptors. The interpretation is presented to the user in the form of the impact of features and can be visualised as a concise summary or in a hierarchical network detailing the full elucidation of the models behaviour for a particular query compound. The interpretation output has been capitalised on and incorporated into a knowledge mining strategy. The knowledge mining is able to extract the learned structure activity relationship trends from a model such as a Random Forest, decision tree, k Nearest Neighbour or support vector machine. These trends can be presented to the user focused around the feature responsible for the assessment such as ACTIVATING or DEACTIVATING. Supporting examples are provided along with an estimation of the models predictive performance for a given SAR trend. Both the interpretation and knowledge mining has been applied to models built for the prediction of Ames mutagenicity and skin irritation. The performance of the developed models is strong and comparable to both academic and commercial predictors for these two toxicological activities.