An updated evolutionary classification of CRISPR–Cas systems

The evolution of CRISPR–cas loci, which encode adaptive immune systems in archaea and bacteria, involves rapid changes, in particular numerous rearrangements of the locus architecture and horizontal transfer of complete loci or individual modules. These dynamics complicate straightforward phylogenetic classification, but here we present an approach combining the analysis of signature protein families and features of the architecture of cas loci that unambiguously partitions most CRISPR–cas loci into distinct classes, types and subtypes. The new classification retains the overall structure of the previous version but is expanded to now encompass two classes, five types and 16 subtypes. The relative stability of the classification suggests that the most prevalent variants of CRISPR–Cas systems are already known. However, the existence of rare, currently unclassifiable variants implies that additional types and subtypes remain to be characterized.

Path-integral and Coarse-graining Strategies for Complex Molecular Phenomena

Molecular simulation provides a powerful tool for connecting molecular-level processes to physical observables. However, the facility to make those connections relies upon the application and development of theoretical methods that permit appropriate descriptions of the systems or processes to be studied. In this thesis, we utilize molecular simulation to study and predict two phenomena with very different theoretical challenges, beginning with (1) lithium-ion transport behavior in polymers and following with (2) equilibrium isotope effects with relevance to position-specific and clumped isotope studies. In the case of ion transport in polymers, there is motivation to use molecular simulation to provide guidance in polymer electrolyte design, but the length and timescales relevant for ion diffusion in polymers preclude the use of direct molecular dynamics simulation to compute ion diffusivities in more than a handful of candidate systems. In the case of equilibrium isotope effects, the thermodynamic driving forces for isotopic fractionation are often fundamentally quantum mechanical in nature, and the high precision of experimental instruments demands correspondingly accurate theoretical approaches. Herein, we describe respectively coarse-graining and path-integral strategies to address outstanding questions in these two subject areas.

PARAMETRIC STUDY OF REAXFF SIMULATION PARAMETERS FOR MOLECULAR DYNAMICS MODELING OF REACTIVE CARBON GASES

Abstract The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although the Reax Force Field (ReaxFF) can be used to simulate the chemical behavior of carbon-based systems, the simulation settings required for accurate predictions have not been fully explored. Using the ReaxFF, molecular dynamics (MD) simulations are used to simulate the chemical behavior of pure carbon and hydrocarbon reactive gases that are involved in the formation of carbon structures such as graphite, buckyballs, amorphous carbon, and carbon nanotubes. It is determined that the maximum simulation time step that can be used in MD simulations with the ReaxFF is dependent on the simulated temperature and selected parameter set, as are the predicted reaction rates. It is also determined that different carbon-based reactive gases react at different rates, and that the predicted equilibrium structures are generally the same for the different ReaxFF parameter sets, except in the case of the predicted formation of large graphitic structures with the Chenoweth parameter set under specific conditions.

MOLECULAR MODELING OF EPON 862-DETDA / CARBON COMPOSITES

The thermoset epoxy resin EPON 862, coupled with the DETDA hardening agent, are utilized as the polymer matrix component in many graphite (carbon fiber) composites. Because it is difficult to experimentally characterize the interfacial region, computational molecular modeling is a necessary tool for understanding the influence of the interfacial molecular structure on bulk-level material properties. The purpose of this research is to investigate the many possible variables that may influence the interfacial structure and the effect they will have on the mechanical behavior of the bulk level composite. Molecular models are established for EPON 862-DETDA polymer in the presence of a graphite surface. Material characteristics such as polymer mass-density, residual stresses, and molecular potential energy are investigated near the polymer/fiber interface. Because the exact degree of crosslinking in these thermoset systems is not known, many different crosslink densities (degrees of curing) are investigated. It is determined that a region exists near the carbon fiber surface in which the polymer mass density is different than that of the bulk mass density. These surface effects extend ~10 Å into the polymer from the center of the outermost graphite layer. Early simulations predict polymer residual stress levels to be higher near the graphite surface. It is also seen that the molecular potential energy in polymer atoms decreases with increasing crosslink density. New models are then established in order to investigate the interface between EPON 862-DETDA polymer and graphene nanoplatelets (GNPs) of various atomic thicknesses. Mechanical properties are extracted from the models using Molecular Dynamics techniques. These properties are then implemented into micromechanics software that utilizes the generalized method of cells to create representations of macro-scale composites. Micromechanics models are created representing GNP doped epoxy with varying number of graphene layers and interfacial polymer crosslink densities. The initial micromechanics results for the GNP doped epoxy are then taken to represent the matrix component and are re-run through the micromechanics software with the addition of a carbon fiber to simulate a GNP doped epoxy/carbon fiber composite. Micromechanics results agree well with experimental data, and indicate GNPs of 1 to 2 atomic layers to be highly favorable. The effect of oxygen bonded to the surface of the GNPs is lastly investigated. Molecular Models are created for systems with varying graphene atomic thickness, along with different amounts of oxygen species attached to them. Models are created for graphene containing hydroxyl groups only, epoxide groups only, and a combination of epoxide and hydroxyl groups. Results show models of oxidized graphene to decrease in both tensile and shear modulus. Attaching only epoxide groups gives the best results for mechanical properties, though pristine graphene is still favored.

Understanding ten-eleven translocation-2 in hematological and nervous systems

I proposed the study of two distinct aspects of Ten-Eleven Translocation 2 (TET2) protein for understanding specific functions in different body systems. ^ In Part I, I characterized the molecular mechanisms of Tet2 in the hematological system. As the second member of Ten-Eleven Translocation protein family, TET2 is frequently mutated in leukemic patients. Previous studies have shown that the TET2 mutations frequently occur in 20% myelodysplastic syndrome/myeloproliferative neoplasm (MDS/MPN), 10% T-cell lymphoma leukemia and 2% B-cell lymphoma leukemia. Genetic mouse models also display distinct phenotypes of various types of hematological malignancies. I performed 5-hydroxymethylcytosine (5hmC) chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq) of hematopoietic stem/progenitor cells to determine whether the deletion of Tet2 can affect the abundance of 5hmC at myeloid, T-cell and B-cell specific gene transcription start sites, which ultimately result in various hematological malignancies. Subsequent Exome sequencing (Exome-Seq) showed that disease-specific genes are mutated in different types of tumors, which suggests that TET2 may protect the genome from being mutated. The direct interaction between TET2 and Mutator S Homolog 6 (MSH6) protein suggests TET2 is involved in DNA mismatch repair. Finally, in vivo mismatch repair studies show that the loss of Tet2 causes a mutator phenotype. Taken together, my data indicate that TET2 binds to MSH6 to protect genome integrity. ^ In Part II, I intended to better understand the role of Tet2 in the nervous system. 5-hydroxymethylcytosine regulates epigenetic modification during neurodevelopment and aging. Thus, Tet2 may play a critical role in regulating adult neurogenesis. To examine the physiological significance of Tet2 in the nervous system, I first showed that the deletion of Tet2 reduces the 5hmC levels in neural stem cells. Mice lacking Tet2 show abnormal hippocampal neurogenesis along with 5hmC alternations at different gene promoters and corresponding gene expression downregulation. Through the luciferase reporter assay, two neural factors Neurogenic differentiation 1 (NeuroD1) and Glial fibrillary acidic protein (Gfap) were down-regulated in Tet2 knockout cells. My results suggest that Tet2 regulates neural stem/progenitor cell proliferation and differentiation in adult brain.^

Environmental scanning electron microscope imaging of vesicle systems

The structural characteristics of liposomes have been widely investigated and there is certainly a strong understanding of their morphological characteristics. Imaging of these systems, using techniques such as freeze-fracturing methods, transmission electron microscopy, and cryo-electron imaging, has allowed us to appreciate their bilayer structures and factors which can influence this. However, there are few methods which all us to study these systems in their natural hydrated state; commonly the liposomes are visualized after drying, staining, and/or fixation of the vesicles. Environmental Scanning Electron Microscopy (ESEM) offers the ability to image a liposome in its hydrated state without the need for prior sample preparation. Within our studies we were the first to use ESEM to study liposomes and niosomes and we have been able to dynamically follow the hydration of lipid films and changes in liposome suspensions as water condenses on to, or evaporates from, the sample in real time. This provides insight into the resistance of liposomes to coalescence during dehydration, thereby providing an alternative assay of liposome formulation and stability.

Molecular Origins of Higher Harmonics in Large-Amplitude Oscillatory Shear Flow: Shear Stress Response

Recent work has focused on deepening our understanding of the molecular origins of the higher harmonics that arise in the shear stress response of polymeric liquids in large-amplitude oscillatory shear flow. For instance, these higher harmonics have been explained by just considering the orientation distribution of rigid dumbbells suspended in a Newtonian solvent. These dumbbells, when in dilute suspension, form the simplest relevant molecular model of polymer viscoelasticity, and this model specifically neglects interactions between the polymer molecules [R.B. Bird et al., J Chem Phys, 140, 074904 (2014)]. In this paper, we explore these interactions by examining the Curtiss-Bird model, a kinetic molecular theory designed specifically to account for the restricted motions that arise when polymer chains are concentrated, thus interacting and specifically, entangled. We begin our comparison using a heretofore ignored explicit analytical solution [Fan and Bird, JNNFM, 15, 341 (1984)]. For concentrated systems, the chain motion transverse to the chain axis is more restricted than along the axis. This anisotropy is described by the link tension coefficient, ε, for which several special cases arise: ε = 0 corresponds to reptation, ε > 1/8 to rod-climbing, 1/2 ≥ ε ≥ 3/4 to reasonable predictions for shear-thinning in steady simple shear flow, and ε = 1 to the dilute solution without hydrodynamic interaction. In this paper, we examine the shapes of the shear stress versus shear rate loops for the special cases ε = (0,1/8, 3/8,1) , and we compare these with those of rigid dumbbell and reptation model predictions.

Identification of novel molecular determinants of tissue mineralization in fish

Tese de doutoramento, Ciências Biomédicas, Departamento de Ciências Biomédicas e Medicina, Universidade do Algarve, 2015

Physicochemical complexity in complex chemical systems

Self-replication and compartmentalization are two central properties thought to be essential for minimal life, and understanding how such processes interact in the emergence of complex reaction networks is crucial to exploring the development of complexity in chemistry and biology. Autocatalysis can emerge from multiple different mechanisms such as formation of an initiator, template self-replication and physical autocatalysis (where micelles formed from the reaction product solubilize the reactants, leading to higher local concentrations and therefore higher rates). Amphiphiles are also used in artificial life studies to create protocell models such as micelles, vesicles and oil-in-water droplets, and can increase reaction rates by encapsulation of reactants. So far, no template self-replicator exists which is capable of compartmentalization, or transferring this molecular scale phenomenon to micro or macro-scale assemblies. Here a system is demonstrated where an amphiphilic imine catalyses its own formation by joining a non-polar alkyl tail group with a polar carboxylic acid head group to form a template, which was shown to form reverse micelles by Dynamic Light Scattering (DLS). The kinetics of this system were investigated by 1H NMR spectroscopy, showing clearly that a template self-replication mechanism operates, though there was no evidence that the reverse micelles participated in physical autocatalysis. Active oil droplets, composed from a mixture of insoluble organic compounds in an aqueous sub-phase, can undergo processes such as division, self-propulsion and chemotaxis, and are studied as models for minimal cells, or protocells. Although in most cases the Marangoni effect is responsible for the forces on the droplet, the behaviour of the droplet depends heavily on the exact composition. Though theoretical models are able to calculate the forces on a droplet, to model a mixture of oils on an aqueous surface where compounds from the oil phase are dissolving and diffusing through the aqueous phase is beyond current computational capability. The behaviour of a droplet in an aqueous phase can only be discovered through experiment, though it is determined by the droplet's composition. By using an evolutionary algorithm and a liquid handling robot to conduct droplet experiments and decide which compositions to test next, entirely autonomously, the composition of the droplet becomes a chemical genome capable of evolution. The selection is carried out according to a fitness function, which ranks the formulation based on how well it conforms to the chosen fitness criteria (e.g. movement or division). Over successive generations, significant increases in fitness are achieved, and this increase is higher with more components (i.e. greater complexity). Other chemical processes such as chemiluminescence and gelation were investigated in active oil droplets, demonstrating the possibility of controlling chemical reactions by selective droplet fusion. Potential future applications for this might include combinatorial chemistry, or additional fitness goals for the genetic algorithm. Combining the self-replication and the droplet protocells research, it was demonstrated that the presence of the amphiphilic replicator lowers the interfacial tension between droplets of a reaction mixture in organic solution and the alkaline aqueous phase, causing them to divide. Periodic sampling by a liquid handling robot revealed that the extent of droplet fission increased as the reaction progressed, producing more individual protocells with increased self-replication. This demonstrates coupling of the molecular scale phenomenon of template self-replication to a macroscale physicochemical effect.

The dehydrogenation reactions between LiH and organic amines for solid state hydrogen storage systems

Hydrogen is considered as an appealing alternative to fossil fuels in the pursuit of sustainable, secure and prosperous growth in the UK and abroad. However there exists a persisting bottleneck in the effective storage of hydrogen for mobile applications in order to facilitate a wide implementation of hydrogen fuel cells in the fossil fuel dependent transportation industry. To address this issue, new means of solid state chemical hydrogen storage are proposed in this thesis. This involves the coupling of LiH with three different organic amines: melamine, urea and dicyandiamide. In principle, thermodynamically favourable hydrogen release from these systems proceeds via the deprotonation of the protic N-H moieties by the hydridic metal hydride. Simultaneously hydrogen kinetics is expected to be enhanced over heavier hydrides by incorporating lithium ions in the proposed binary hydrogen storage systems. Whilst the concept has been successfully demonstrated by the results obtained in this work, it was observed that optimising the ball milling conditions is central in promoting hydrogen desorption in the proposed systems. The theoretical amount of 6.97 wt% by dry mass of hydrogen was released when heating a ball milled mixture of LiH and melamine (6:1 stoichiometry) to 320 °C. It was observed that ball milling introduces a disruption in the intermolecular hydrogen bonding network that exists in pristine melamine. This effect extends to a molecular level electron redistribution observed as a function of shifting IR bands. It was postulated that stable phases form during the first stages of dehydrogenation which contain the triazine skeleton. Dehydrogenation of this system yields a solid product Li2NCN, which has been rehydrogenated back to melamine via hydrolysis under weak acidic conditions. On the other hand, the LiH and urea system (4:1 stoichiometry) desorbed approximately 5.8 wt% of hydrogen, from the theoretical capacity of 8.78 wt% (dry mass), by 270 °C accompanied by undesirable ammonia and trace amount of water release. The thermal dehydrogenation proceeds via the formation of Li(HN(CO)NH2) at 104.5 °C; which then decomposes to LiOCN and unidentified phases containing C-N moieties by 230 °C. The final products are Li2NCN and Li2O (270 °C) with LiCN and Li2CO3 also detected under certain conditions. It was observed that ball milling can effectively supress ammonia formation. Furthermore results obtained from energetic ball milling experiments have indicated that the barrier to full dehydrogenation between LiH and urea is principally kinetic. Finally the dehydrogenation reaction between LiH and dicyandiamide system (4:1 stoichiometry) occurs through two distinct pathways dependent on the ball milling conditions. When ball milled at 450 RPM for 1 h, dehydrogenation proceeds alongside dicyandiamide condensation by 400 °C whilst at a slower milling speed of 400 RPM for 6h, decomposition occurs via a rapid gas desorption (H2 and NH3) at 85 °C accompanied by sample foaming. The reactant dicyandiamide can be generated by hydrolysis using the product Li2NCN.

Applications of metrological techniques for clinical implementation of dosimetry and radiobiology in molecular radiotherapy

Molecular radiotherapy (MRT) is a fast developing and promising treatment for metastasised neuroendocrine tumours. Efficacy of MRT is based on the capability to selectively "deliver" radiation to tumour cells, minimizing administered dose to normal tissues. Outcome of MRT depends on the individual patient characteristics. For that reason, personalized treatment planning is important to improve outcomes of therapy. Dosimetry plays a key role in this setting, as it is the main physical quantity related to radiation effects on cells. Dosimetry in MRT consists in a complex series of procedures ranging from imaging quantification to dose calculation. This doctoral thesis focused on several aspects concerning the clinical implementation of absorbed dose calculations in MRT. Accuracy of SPECT/CT quantification was assessed in order to determine the optimal reconstruction parameters. A model of PVE correction was developed in order to improve the activity quantification in small volume, such us lesions in clinical patterns. Advanced dosimetric methods were compared with the aim of defining the most accurate modality, applicable in clinical routine. Also, for the first time on a large number of clinical cases, the overall uncertainty of tumour dose calculation was assessed. As part of the MRTDosimetry project, protocols for calibration of SPECT/CT systems and implementation of dosimetry were drawn up in order to provide standard guidelines to the clinics offering MRT. To estimate the risk of experiencing radio-toxicity side effects and the chance of inducing damage on neoplastic cells is crucial for patient selection and treatment planning. In this thesis, the NTCP and TCP models were derived based on clinical data as help to clinicians to decide the pharmaceutical dosage in relation to the therapy control and the limitation of damage to healthy tissues. Moreover, a model for tumour response prediction based on Machine Learning analysis was developed.

Peptide-based low molecular weight gelators for the preparation of self-assembled materials

Low molecular weight gelators (LMWGs) based on pseudo-peptides are here studied for the preparation of supramolecular materials. These compounds can self-assemble through non-covalent interactions such as hydrogen bonds and π-π stacking, forming fibres and gels. A wide variety of materials can be prepared starting from these building blocks, which can be tuned and functionalised depending on the application. In this work, derivatives of the three aromatic amino acids L-Phenylalanine, L-Tyrosine and L-DOPA (3,4-dihydroxiphenylalanine) were synthesised and tested as gelators for water or organic solvents. First, the optimal gelating conditions were studied for each compound, varying concentration, solvent and trigger. Then the materials were characterised in terms of mechanical properties and morphology. Water remediation from dye pollution was the first focus of this work. Organogels were studied as absorbent of dyes from contaminated water. Hydrogels functionalised with TiO2 nanoparticles and graphene platelets were proposed as efficient materials for the photo-degradation of dyes. An efficient method for the incorporation of graphene inside hydrogels using the gelator itself as dispersant was proposed. In these materials a high storage modulus coexists with good self-healing and biocompatibility. The incorporation of a mineral phase inside the gel matrix was then investigated, leading to the preparation of composite organic/inorganic materials. In a first study, the growth of calcium carbonate crystals was achieved inside the hydrogel, which preserved its structure after crystal formation. Then the self-assembled fibres made of LMWGs were used for the first time instead of the polymeric ones as reinforcement inside calcium phosphate cements (CPCs) for bone regeneration. Gel-to-crystal transitions occurring with time in a metastable gel were also examined. The formation of organic crystals in gels can be achieved in multicomponent systems, in which a second gelator constitutes the independent gel network. Finally, some compounds unable to gelate were tested as underwater adhesives.

Studying the alignment of 5-cyanobiphenyl on a polydimethylsiloxane surface through Molecular Dynamics simulations

For 40 years, at the University of Bologna, a group of researchers coordinated by professor Claudio Zannoni has been studying liquid crystals by employing computational techniques. They have developed effective models of these interesting, and still far from being completely understood, systems. They were able to reproduce with simulations important features of some liquid crystal molecules, such as transition temperature. Then they focused their attention on the interactions that these molecules have with different kinds of surface, and how these interactions affect the alignment of liquid crystals. The group studied the behaviour of liquid crystals in contact with different kinds of surfaces, from silica, either amorphous and crystalline, to organic self assembled monolayers (SAMs) and even some common polymers, such as polymethylmethacrylate (PMMA) and polystyrene (PS). Anyway, a library of typical surfaces is still far from being complete, and a lot of work must be done to investigate the cases which have not been analyzed yet. A hole that must be filled is represented by polydimethylsiloxane (PDMS), a polymer on which the interest of industry has enormously grown up in the last years, thanks to its peculiar features, allowing it to be employed in many fields of applications. It has been observed experimentally that PDMS causes 4-cyano-4’-pentylbiphenyl (well known as 5CB), one of the most common liquid crystal molecules, to align homeotropically (i.e. perpendicular) with respect to a surface made of this polymer. Even though some hypothesis have been presented to rationalize the effect, a clear explanation of this phenomenon has not been given yet. This dissertation shows the work I did during my internship in the group of professor Zannoni. The challenge that I had to tackle was to investigate, via Molecular Dynamics (MD) simulations, the reasons of 5CB homeotropic alignment on a PDMS surface, as the group had previously done for other surfaces.

Molecular simulations and modeling of pharmaceutically relevant RNA-processing enzymes

The two-metal-ion architecture is a structural feature found in a variety of RNA processing metalloenzymes or ribozymes (RNA-based enzymes), which control the biogenesis and the metabolism of vital RNAs, including non-coding RNAs (ncRNAs). Notably, such ncRNAs are emerging as key players for the regulation of cellular homeostasis, and their altered expression has been often linked to the development of severe human pathologies, from cancer to mental disorders. Accordingly, understanding the biological processing of ncRNAs is foundational for the development of novel therapeutic strategies and tools. Here, we use state-of the-art molecular simulations, complemented with X-ray crystallography and biochemical experiments, to characterize the RNA processing cycle as catalyzed by two two-metal-ion enzymes: the group II intron ribozymes and the RNase H1. We show that multiple and diverse cations are strategically recruited at and timely released from the enzymes’ active site during catalysis. Such a controlled cations’ trafficking leads to the recursive formation and disruption of an extended two-metal ion architecture that is functional for RNA-hydrolysis – from substrate recruitment to product release. Importantly, we found that these cations’ binding sites are conserved among other RNA-processing machineries, including the human spliceosome and CRISPR-Cas systems, suggesting that an evolutionarily-converged catalytic strategy is adopted by these enzymes to process RNA molecules. Thus, our findings corroborate and sensibly extend the current knowledge of two-metal-ion enzymes, and support the design of novel drugs targeting RNA-processing metalloenzymes or ribozymes as well as the rational engineering of novel programmable gene-therapy tools.

Photoactive (supra)molecular devices in double layer membranes

Biological systems are complex and highly organized architectures governed by non-covalent interactions responsible for the regulation of essential tasks in all living organisms. These systems are a constant source of inspiration for supramolecular chemists aiming to design multicomponent molecular assemblies able to perform elaborated tasks, thanks to the role and action of the components that constitute them. Artificial supramolecular systems exploit non-covalent interactions to mimic naturally occurring events. In this context, stimuli-responsive supramolecular systems have attracted attention due to the possibility to control macroscopic effects through modifications at the nanoscale. This thesis is divided in three experimental chapters, characterized by a progressive increase in molecular complexity. Initially, the preparation and studies of liposomes functionalized with a photoactive guest such as azobenzene in the bilayer were tackled, in order to evaluate the effect of such photochrome on the vesicle properties. Subsequently, the synthesis and studies of thread-like molecules comprising an azobenzene functionality was reported. Such molecules were conceived to be intercalated in the bilayer membrane of liposomes with the aim to be used as components for photoresponsive transmembrane molecular pumps. Finally, a [3]rotaxane was developed and studied in solution. This system is composed of two crown ether rings interlocked with an axle containing three recognition sites for the macrocycles, i.e. two pH-switchable ammonium stations and a permanent triazolium station. Such molecule was designed to achieve a change in the ratio between the recognition sites and the crown ethers as a consequence of acid-base inputs. This leads to the formation of rotaxanes containing a number of recognition sites respectively larger, equal or lower than the number of interlocked rings and connected by a network of acid-base reactions.