Conformation and Stereodynamic of Hindered Aromatic System

The preparation of conformationally hindered molecules and their study by DNMR and computational methods are my thesis’s core. In the first chapter, the conformations and the stereodynamics of symmetrically ortho-disubstituted aryl carbinols and aryl ethers are described. In the second chapter, the structures of axially chiral atropisomers of hindered biphenyl carbinols are studied. In the third chapter, the steric barriers and the -barrier of 1,8-di-aylbiphenylenes are determined. Interesting atropisomers are found in the cases of arylanthrones, arylanthraquinones and arylanthracenes and are reported in the fourth chapter. By the combined use of dynamic NMR, ECD spectroscopy and DFT computations, the conformations and the absolute configurations of 2-Naphthylalkylsulfoxides are studied in the fifth chapter. In the last chapter, a new synthetic route to ,’-arylated secondary or tertiary alcohols by lithiated O-benzyl-carbamates carrying an N-aryl substituent and DFT calculations to determinate the cyclic intermediate are reported. This work was done in the research group of Prof. Jonathan Clayden, at the University of Manchester.

Quantum effects and dynamics in hydrogen-bonded systems: a first-principles approach to spectroscopic experiments

Computer simulations have become an important tool in physics. Especially systems in the solid state have been investigated extensively with the help of modern computational methods. This thesis focuses on the simulation of hydrogen-bonded systems, using quantum chemical methods combined with molecular dynamics (MD) simulations. MD simulations are carried out for investigating the energetics and structure of a system under conditions that include physical parameters such as temperature and pressure. Ab initio quantum chemical methods have proven to be capable of predicting spectroscopic quantities. The combination of these two features still represents a methodological challenge. Furthermore, conventional MD simulations consider the nuclei as classical particles. Not only motional effects, but also the quantum nature of the nuclei are expected to influence the properties of a molecular system. This work aims at a more realistic description of properties that are accessible via NMR experiments. With the help of the path integral formalism the quantum nature of the nuclei has been incorporated and its influence on the NMR parameters explored. The effect on both the NMR chemical shift and the Nuclear Quadrupole Coupling Constants (NQCC) is presented for intra- and intermolecular hydrogen bonds. The second part of this thesis presents the computation of electric field gradients within the Gaussian and Augmented Plane Waves (GAPW) framework, that allows for all-electron calculations in periodic systems. This recent development improves the accuracy of many calculations compared to the pseudopotential approximation, which treats the core electrons as part of an effective potential. In combination with MD simulations of water, the NMR longitudinal relaxation times for 17O and 2H have been obtained. The results show a considerable agreement with the experiment. Finally, an implementation of the calculation of the stress tensor into the quantum chemical program suite CP2K is presented. This enables MD simulations under constant pressure conditions, which is demonstrated with a series of liquid water simulations, that sheds light on the influence of the exchange-correlation functional used on the density of the simulated liquid.

Subject-specific musculoskeletal models of the lower limbs for the prediction of skeletal loads during motion

The determination of skeletal loading conditions in vivo and their relationship to the health of bone tissues, remain an open question. Computational modeling of the musculoskeletal system is the only practicable method providing a valuable approach to muscle and joint loading analyses, although crucial shortcomings limit the translation process of computational methods into the orthopedic and neurological practice. A growing attention focused on subject-specific modeling, particularly when pathological musculoskeletal conditions need to be studied. Nevertheless, subject-specific data cannot be always collected in the research and clinical practice, and there is a lack of efficient methods and frameworks for building models and incorporating them in simulations of motion. The overall aim of the present PhD thesis was to introduce improvements to the state-of-the-art musculoskeletal modeling for the prediction of physiological muscle and joint loads during motion. A threefold goal was articulated as follows: (i) develop state-of-the art subject-specific models and analyze skeletal load predictions; (ii) analyze the sensitivity of model predictions to relevant musculotendon model parameters and kinematic uncertainties; (iii) design an efficient software framework simplifying the effort-intensive phases of subject-specific modeling pre-processing. The first goal underlined the relevance of subject-specific musculoskeletal modeling to determine physiological skeletal loads during gait, corroborating the choice of full subject-specific modeling for the analyses of pathological conditions. The second goal characterized the sensitivity of skeletal load predictions to major musculotendon parameters and kinematic uncertainties, and robust probabilistic methods were applied for methodological and clinical purposes. The last goal created an efficient software framework for subject-specific modeling and simulation, which is practical, user friendly and effort effective. Future research development aims at the implementation of more accurate models describing lower-limb joint mechanics and musculotendon paths, and the assessment of an overall scenario of the crucial model parameters affecting the skeletal load predictions through probabilistic modeling.

Chiroptical properties of bioactive molecules: sensitivity to conformation and solvation

Chiroptical spectroscopies play a fundamental role in pharmaceutical analysis for the stereochemical characterisation of bioactive molecules, due to the close relationship between chirality and optical activity and the increasing evidence of stereoselectivity in the pharmacological and toxicological profiles of chiral drugs. The correlation between chiroptical properties and absolute stereochemistry, however, requires the development of accurate and reliable theoretical models. The present thesis will report the application of theoretical chiroptical spectroscopies in the field of drug analysis, with particular emphasis on the huge influence of conformational flexibility and solvation on chiroptical properties and on the main computational strategies available to describe their effects by means of electronic circular dichroism (ECD) spectroscopy and time-dependent density functional theory (TD-DFT) calculations. The combination of experimental chiroptical spectroscopies with state-of-the-art computational methods proved to be very efficient at predicting the absolute configuration of a wide range of bioactive molecules (fluorinated 2-arylpropionic acids, β-lactam derivatives, difenoconazole, fenoterol, mycoleptones, austdiol). The results obtained for the investigated systems showed that great care must be taken in describing the molecular system in the most accurate fashion, since chiroptical properties are very sensitive to small electronic and conformational perturbations. In the future, the improvement of theoretical models and methods, such as ab initio molecular dynamics, will benefit pharmaceutical analysis in the investigation of non-trivial effects on the chiroptical properties of solvated systems and in the characterisation of the stereochemistry of complex chiral drugs.

Metodi computazionali per l'annotazione di genomi e proteomi

Il progresso tecnologico nel campo della biologia molecolare, pone la comunità scientifica di fronte all’esigenza di dare un’interpretazione all’enormità di sequenze biologiche che a mano a mano vanno a costituire le banche dati, siano esse proteine o acidi nucleici. In questo contesto la bioinformatica gioca un ruolo di primaria importanza. Un nuovo livello di possibilità conoscitive è stato introdotto con le tecnologie di Next Generation Sequencing (NGS), per mezzo delle quali è possibile ottenere interi genomi o trascrittomi in poco tempo e con bassi costi. Tra le applicazioni del NGS più rilevanti ci sono senza dubbio quelle oncologiche che prevedono la caratterizzazione genomica di tessuti tumorali e lo sviluppo di nuovi approcci diagnostici e terapeutici per il trattamento del cancro. Con l’analisi NGS è possibile individuare il set completo di variazioni che esistono nel genoma tumorale come varianti a singolo nucleotide, riarrangiamenti cromosomici, inserzioni e delezioni. Va però sottolineato che le variazioni trovate nei geni vanno in ultima battuta osservate dal punto di vista degli effetti a livello delle proteine in quanto esse sono le responsabili più dirette dei fenotipi alterati riscontrabili nella cellula tumorale. L’expertise bioinformatica va quindi collocata sia a livello dell’analisi del dato prodotto per mezzo di NGS ma anche nelle fasi successive ove è necessario effettuare l’annotazione dei geni contenuti nel genoma sequenziato e delle relative strutture proteiche che da esso sono espresse, o, come nel caso dello studio mutazionale, la valutazione dell’effetto della variazione genomica. È in questo contesto che si colloca il lavoro presentato: da un lato lo sviluppo di metodologie computazionali per l’annotazione di sequenze proteiche e dall’altro la messa a punto di una pipeline di analisi di dati prodotti con tecnologie NGS in applicazioni oncologiche avente come scopo finale quello della individuazione e caratterizzazione delle mutazioni genetiche tumorali a livello proteico.

Beyond simple proton transfer processes: domino reactions between 4-substituted indoles and nitroethene as a new gateway to ergot alkaloids

The multimodal biology activity of ergot alkaloids is known by humankind since middle ages. Synthetically modified ergot alkaloids are used for the treatment of various medical conditions. Despite the great progress in organic syntheses, the total synthesis of ergot alkaloids remains a great challenge due to the complexity of their polycyclic structure with multiple stereogenic centres. This project has developed a new domino reaction between indoles bearing a Michael acceptor at the 4 position and nitroethene, leading to potential ergot alkaloid precursors in highly enantioenriched form. The reaction was optimised and applied to a large variety of substrate with good results. Even if unfortunately all attempts to further modify the obtained polycyclic structure failed, it was found a reaction able to produce the diastereoisomer of the polycyclic product in excellent yields. The compounds synthetized were characterized by NMR and ESIMS analysis confirming the structure and their enantiomeric excess was determined by chiral stationary phase HPLC. The mechanism of the reaction was evaluated by DFT calculations, showing the formation of a key bicoordinated nitronate intermediate, and fully accounting for the results observed with all substrates. The relative and absolute configuration of the adducts were determined by a combination of NMR, ECD and computational methods.

EPR spectroscopy as a tool for the elucidation of non-covalent structuring principles in complex soft matter

This thesis aims at connecting structural and functional changes of complex soft matter systems due to external stimuli with non-covalent molecular interaction profiles. It addresses the problem of elucidating non-covalent forces as structuring principle of mainly polymer-based systems in solution. The structuring principles of a wide variety of complex soft matter types are analyzed. In many cases this is done by exploring conformational changes upon the exertion of external stimuli. The central question throughout this thesis is how a certain non-covalent interaction profile leads to solution condition-dependent structuring of a polymeric system.rnTo answer this question, electron paramagnetic resonance (EPR) spectroscopy is chosen as the main experimental method for the investigation of the structure principles of polymers. With EPR one detects only the local surroundings or environments of molecules that carry an unpaired electron. Non-covalent forces are normally effective on length scales of a few nanometers and below. Thus, EPR is excellently suited for their investigations. It allows for detection of interactions on length scales ranging from approx. 0.1 nm up to 10 nm. However, restriction to only one experimental technique likely leads to only incomplete pictures of complex systems. Therefore, the presented studies are frequently augmented with further experimental and computational methods in order to yield more comprehensive descriptions of the systems chosen for investigation.rnElectrostatic correlation effects in non-covalent interaction profiles as structuring principles in colloid-like ionic clusters and DNA condensation are investigated first. Building on this it is shown how electrostatic structuring principles can be combined with hydrophobic ones, at the example of host-guest interactions in so-called dendronized polymers (denpols).rnSubsequently, the focus is shifted from electrostatics in dendronized polymers to thermoresponsive alkylene oxide-based materials, whose structuring principles are based on hydrogen bonds and counteracting hydrophobic interactions. The collapse mechanism in dependence of hydrophilic-hydrophobic balance and topology of these polymers is elucidated. Complementarily the temperature-dependent phase behavior of elastin-like polypeptides (ELPs) is investigated. ELPs are the first (and so far only) class of compounds that is shown to feature a first-order inverse phase transition on nanoscopic length scales.rnFinally, this thesis addresses complex biological systems, namely intrinsically disordered proteins (IDPs). It is shown that the conformational space of the IDPs Osteopontin (OPN), a cytokine involved in metastasis of several kinds of cancer, and BASP1 (brain acid soluble protein one), a protein associated with neurite outgrowth, is governed by a subtle interplay between electrostatic forces, hydrophobic interaction, system entropy and hydrogen bonds. Such, IDPs can even sample cooperatively folded structures, which have so far only been associated with globular proteins.

Hypersonic phonon propagation in mesoscopic systems by Brillouin spectroscopy

Phononic crystals, capable to block or direct the propagation of elastic/acoustic waves, have attracted increasing interdisciplinary interest across condensed matter physics and materials science. As of today, no generalized full description of elastic wave propagation in phononic structures is available, mainly due to the large number of variables determining the band diagram. Therefore, this thesis aims for a deeper understanding of the fundamental concepts governing wave propagation in mesoscopic structures by investigation of appropriate model systems. The phononic dispersion relation at hypersonic frequencies is directly investigated by the non-destructive technique of high-resolution spontaneous Brillouin light scattering (BLS) combined with computational methods. Due to the vector nature of the elastic wave propagation, we first studied the hypersonic band structure of hybrid superlattices. These 1D phononic crystals composed of alternating layers of hard and soft materials feature large Bragg gaps. BLS spectra are sensitive probes of the moduli, photo-elastic constants and structural parameters of the constituent components. Engineering of the band structure can be realized by introduction of defects. Here, cavity layers are employed to launch additional modes that modify the dispersion of the undisturbed superlattice, with extraordinary implications to the band gap region. Density of states calculations in conjunction with the associated deformation allow for unambiguous identication of surface and cavity modes, as well as their interaction with adjacent defects. Next, the role of local resonances in phononic systems is explored in 3D structures based on colloidal particles. In turbid media BLS records the particle vibration spectrum comprising resonant modes due to the spatial confinement of elastic energy. Here, the frequency and lineshapes of the particle eigenmodes are discussed as function of increased interaction and departure from spherical symmetry. The latter is realized by uniaxial stretching of polystyrene spheres, that can be aligned in an alternating electric field. The resulting spheroidal crystals clearly exhibit anisotropic phononic properties. Establishing reliable predictions of acoustic wave propagation, necessary to advance, e.g., optomechanics and phononic devices is the ultimate aim of this thesis.

Esperienza di creazione di entrate lessicografiche combinatorie: metodi e dati dal progetto CombiNet

The present dissertation aims at simulating the construction of lexicographic layouts for an Italian combinatory dictionary based on real linguistic data, extracted from corpora by using computational methods. This work is based on the assumption that the intuition of the native speaker, or the lexicographer, who manually extracts and classifies all the relevant data, are not adequate to provide sufficient information on the meaning and use of words. Therefore, a study of the real use of language is required and this is particularly true for dictionaries that collect the combinatory behaviour of words, where the task of the lexicographer is to identify typical combinations where a word occurs. This study is conducted in the framework of the CombiNet project aimed at studying Italian Word Combinationsand and at building an online, corpus-based combinatory lexicographic resource for the Italian language. This work is divided into three chapters. Chapter 1 describes the criteria considered for the classification of word combinations according to the work of Ježek (2011). Chapter 1 also contains a brief comparison between the most important Italian combinatory dictionaries and the BBI Dictionary of Word Combinations in order to describe how word combinations are considered in these lexicographic resources. Chapter 2 describes the main computational methods used for the extraction of word combinations from corpora, taking into account the advantages and disadvantages of the two methods. Chapter 3 mainly focuses on the practical word carried out in the framework of the CombiNet project, with reference to the tools and resources used (EXTra, LexIt and "La Repubblica" corpus). Finally, the data extracted and the lexicographic layout of the lemmas to be included in the combinatory dictionary are commented, namely the words "acqua" (water), "braccio" (arm) and "colpo" (blow, shot, stroke).

The binding mode of side chain- and C3-modified epothilones to tubulin

The tubulin-binding mode of C3- and C15-modified analogues of epothilone A (Epo A) was determined by NMR spectroscopy and computational methods and compared with the existing structural models of tubulin-bound natural Epo A. Only minor differences were observed in the conformation of the macrocycle between Epo A and the C3-modified analogues investigated. In particular, 3-deoxy- (compound 2) and 3-deoxy-2,3-didehydro-Epo A (3) were found to adopt similar conformations in the tubulin-binding cleft as Epo A, thus indicating that the 3-OH group is not essential for epothilones to assume their bioactive conformation. None of the available models of the tubulin-epothilone complex is able to fully recapitulate the differences in tubulin-polymerizing activity and microtubule-binding affinity between C20-modified epothilones 6 (C20-propyl), 7 (C20-butyl), and 8 (C20-hydroxypropyl). Based on the results of transferred NOE experiments in the presence of tubulin, the isomeric C15 quinoline-based Epo B analogues 4 and 5 show very similar orientations of the side chain, irrespective of the position of the nitrogen atom in the quinoline ring. The quinoline side chain stacks on the imidazole moiety of beta-His227 with equal efficiency in both cases, thus suggesting that the aromatic side chain moiety in epothilones contributes to tubulin binding through strong van der Waals interactions with the protein rather than hydrogen bonding involving the heteroaromatic nitrogen atom. These conclusions are in line with existing tubulin polymerization and microtubule-binding data for 4, 5, and Epo B.

The Search for Low Energy Conformational Families of Small Peptides: Searching for Active Conformations of Small Peptides in the Absence of a Known Receptor

Breast cancer is the most common cancer among women. Tamoxifen is the preferred drug for estrogen receptor-positive breast cancer treatment, yet many of these cancers are intrinsically resistant to tamoxifen or acquire resistance during treatment. Therefore, scientists are searching for breast cancer drugs that have different molecular targets. Previous work revealed that 8-mer and cyclic 9-mer peptides inhibit breast cancer in mouse and rat model systems, interacting with an unknown receptor, while peptides smaller than eight amino acids did not inhibit breast cancer. We have shown that the use of replica exchange molecular dynamics predicts structure and dynamics of active peptides, leading to the discovery of smaller peptides with full biological activity. These simulations identified smaller peptide analogs with a conserved turn, a β-turn formed in the larger peptides. These analogs inhibit estrogen-dependent cell growth in a mouse uterine growth assay, a test showing reliable correlation with human breast cancer inhibition. We outline the computational methods that were tried and used with the experimental information that led to the successful completion of this research.

Evidence for Long-lived Excited States of CnH22+ Carbodications

The energetics, structures, stabilities and reactivities of[CnH2]2+ ions have been investigated using computational methods and experimental mass spectrometric techniques. Spontaneous decompositions of [CnH2]2+ into [CnH]+ + H+ products, observed for ions with odd-n values, have been explained by invoking the formation of excited triplet states. Even-n [CnH]+ ions possess triplet ground states with low-lying excited states, whereas odd-n ions have triplet states with energies several eV above ground singlet states. Radiationless transitions of vibrationally excited long-lived triplet state ions into singlet state continua are suggested as possible mechanisms for spontaneous deprotonation processes of odd-n [CnH2]2+ ions. Evidence for these long-lived excited states has been obtained in bimolecular single electron transfer reactions.

Mechanical modeling of Intraneural Ganglion Cyst

Intraneural Ganglion Cyst is a 200 year old mystery related to nerve injury which is yet to be solved. Current treatments for the above problem are relatively simple procedures related to removal of cystic contents from the nerve. However, these treatments may result into neuropathic pain and recurrence of the cyst. The articular theory proposed by Spinner et al., (Spinner et al. 2003) takes into consideration the neurological deficit in Common Peroneal Nerve (CPN) branch of the sciatic nerve and affirms that in addition to the above treatments, ligation of articular branch results into foolproof eradication of the deficit. Mechanical Modeling of the Affected Nerve Cross Section will reinforce the articular theory (Spinner et al. 2003). As the cyst propagates, it compresses the neighboring fascicles and the nerve cross section appears like a signet ring. Hence, in order to mechanically model the affected nerve cross section; computational methods capable of modeling excessively large deformations are required. Traditional FEM produces distorted elements while modeling such deformations, resulting into inaccuracies and premature termination of the analysis. The methods described in this Master’s Thesis are effective enough to be able to simulate such deformations. The results obtained from the model adequately resemble the MRI image obtained at the same location and shows an appearance of a signet ring. This Master’s Thesis describes the neurological deficit in brief followed by detail explanation of the advanced computational methods used to simulate this problem. Finally, qualitative results show the resemblance of mechanical model to MRI images of the Nerve Cross Section at the same location validating the capability of these methods to study this neurological deficit.

DISSECTION OF STRESS RESPONSE NETWORKS REGULATING MULTIPLE STRESSES IN RICE

Important food crops like rice are constantly exposed to various stresses that can have devastating effect on their survival and productivity. Being sessile, these highly evolved organisms have developed elaborate molecular machineries to sense a mixture of stress signals and elicit a precise response to minimize the damage. However, recent discoveries revealed that the interplay of these stress regulatory and signaling molecules is highly complex and remains largely unknown. In this work, we conducted large scale analysis of differential gene expression using advanced computational methods to dissect regulation of stress response which is at the heart of all molecular changes leading to the observed phenotypic susceptibility. One of the most important stress conditions in terms of loss of productivity is drought. We performed genomic and proteomic analysis of epigenetic and miRNA mechanisms in regulation of drought responsive genes in rice and found subsets of genes with striking properties. Overexpressed genesets included higher number of epigenetic marks, miRNA targets and transcription factors which regulate drought tolerance. On the other hand, underexpressed genesets were poor in above features but were rich in number of metabolic genes with multiple co-expression partners contributing majorly towards drought resistance. Identification and characterization of the patterns exhibited by differentially expressed genes hold key to uncover the synergistic and antagonistic components of the cross talk between stress response mechanisms. We performed meta-analysis on drought and bacterial stresses in rice and Arabidopsis, and identified hundreds of shared genes. We found high level of conservation of gene expression between these stresses. Weighted co-expression network analysis detected two tight clusters of genes made up of master transcription factors and signaling genes showing strikingly opposite expression status. To comprehensively identify the shared stress responsive genes between multiple abiotic and biotic stresses in rice, we performed meta-analyses of microarray studies from seven different abiotic and six biotic stresses separately and found more than thirteen hundred shared stress responsive genes. Various machine learning techniques utilizing these genes classified the stresses into two major classes' namely abiotic and biotic stresses and multiple classes of individual stresses with high accuracy and identified the top genes showing distinct patterns of expression. Functional enrichment and co-expression network analysis revealed the different roles of plant hormones, transcription factors in conserved and non-conserved genesets in regulation of stress response.