Flow machines are very important to industry, being widely used on various processes. Performance improvements are relevant factors and can be achieved by using optimization methods, such as topology optimization. Thus, this work aims to develop a method to design radial flow machine rotors operating on laminar regime, by implementing a topology optimization formulation based on density model. The design of a rotor involves firstly modelling the fluid flow by using the Navier-Stokes equations on a rotating reference frame and using the Finite Element Method for solving the differential equations. To determine the material distribution on the domain, a porous flow model based on the Darcy equation is employed by using an inverse permeability that interpolates between fluid and solid. In the optimization phase, it is defined a multi-objective function that aims to minimize the viscous energy dissipation, vorticity and power. The optimization problem is implemented using the FEniCS environment and the libraries dolfin-adjoint and pyIpopt. The optimized topologies are verified with the ANSYS software. The resulting topologies are post-processed and a CAD model is created. The rotors are manufactured by using a 3D printer, the complete prototype is built by coupling an electric brushless motor and an experimental characterization is performed by measuring fluid flow and pressure head given by the pumps. Experimental and computational results are compared and the improvement is verified. / Máquinas de fluxo são muito importantes para a indústria, sendo utilizadas em diversos processos. Assim, melhorias de desempenho são fatores relevantes e podem ser alcançadas com a utilização de métodos de otimização, como a otimização topológica. Este trabalho visa desenvolver uma metodologia para projetar rotores de máquinas de fluxo radiais que operam em escoamento laminar implementando-se a formulação de otimização topológica baseada no modelo de densidades. O projeto de rotores envolve, primeiramente, a modelagem do escoamento utilizando-se as equações de Navier-Stokes em um referencial rotativo e a utilização do Método de Elementos Finitos para a resolução das equações diferenciais. A distribuição de material no domínio é feita empregando-se um modelo de escoamento em meio poroso baseado nas equações de Darcy, utilizando-se a permeabilidade inversa que interpola o elemento entre sólido e fluido. Na fase de otimização é definida uma função multi-objetivo, que visa minimizar dissipação de energia viscosa, a vorticidade e a potência. O problema de otimização é implementado utilizando-se o ambiente FEniCS para a resolução do sistema de elementos finitos e as bibliotecas dolfin-adjoint e pyIpopt para o algorithmo de otimização. As topologias otimizadas são verificadas com o software ANSYS. As topologias resultantes são pós-processadas para a criação de um modelo CAD dos rotores. Os rotores são construídos utilizando-se a impressão 3D, o protótipo completo é montado acoplando-se um motor elétrico sem escovas e a caracterização experimental é feita medindo-se a vazão e o ganho de pressão dados pelas bombas. Por fim, os resultados experimentais e computacionais são comparados e uma melhoria de desempenho é observada.

As cisteíno-proteases da família da papaína desempenham funções essenciais em processos biológicos, entre eles o desenvolvimento e crescimento do organismo, vias de sinalização celular e apoptose, invasão de parasitas em células hospedeiras. Assim, trata-se de uma classe de proteínas de grande interesse para as indústrias farmacêuticas, sendo utilizada como alvo para o tratamento de doenças como o câncer e metástases, osteoporose. Disfunções relacionadas ao sistema imune, doenças parasitárias como malária, leishmaniose, doença do sono e doença de Chagas. Esta última é uma enfermidade considerada negligenciada pelas grandes indústrias farmacêuticas, sem nenhum tratamento eficaz e seguro disponível, que gera um problema econômico de mais de sete bilhões de dólares anuais devido à perda de mão de obra e gastos com tratamento para amenizar os efeitos da doença. A cisteíno protease cruzaína de Trypanosoma cruzi, causador da Doença de Chagas, desponta como um alvo validado na busca de novos fármacos contra essa enfermidade. Essa enzima apresenta um par de aspartatos que interagem entre si, para os quais foi predito um pKa de 7, sendo possível a forma desprotonada desse par em condições biológicas. Neste caso, pode levar à exposição de uma nova cavidade por meio do movimento da alça entre os resíduos 57-62, segundo as simulações de dinâmica molecular desse trabalho, que se trata de uma possível candidata a ponto de seletividade de inibidores de cisteínoproteases de parasitos em relação às suas ortólogas em Homo sapiens que não possuem o par de aspartatos. Em pH ácido, foi mostrado por meio de análise de componentes principais de simulações de dinâmica molecular que as cisteíno protease apresentam uma restrição gradual na amostragem conformacional do sítio ativo quando complexadas com as formas não covalente e covalente de inibidores derivados de dipeptidil nitrilas. Isso sugere que esse sistema segue o modelo de seleção conformacional para flexibilidade de proteína. Notou-se também que o perfil de restrição de ligantes que inibem na faixa de nmol.L-1 difere daqueles a µmol.L-1 , o que possibilitou a construção de uma árvore de decisão para identificar os complexos que apresentam afinidade a nmol.L-1 . / The papain-like cysteine proteases are essentials for biological process, performing important roles on the parasite development, growth and also in the parasite invasion process on the host cell, in cellular signaling pathways and apoptosis, among others. Thus, the pharmaceutical industry widely uses this class of protein as target for the development of new drugs, against cancer and metastasis, osteoporosis and immune system disorders, resulting in many approved drugs. Additionally, these enzymes are validated target against parasitic diseases as leishmaniose, malaria and African and American trypanosomiasis. The last one, also known as Chagas\' disease, is neglected disease for which, further a century form this discovery, there is no effective and safe chemotherapy and is responsible for an economic loss of around seven billion dollars in the world per year due to the health care and lost productivity from infected people. Faced with this situation, the Cruzain, a cysteine protease from the Trypanosoma cruzi, the causative agent of Chagas\' disease, is emerging as interesting and validated target to the search for new drugs against this sickness. This enzyme has a pair of interactiong Asp for which was predicted a pKa of 7, by computational methods. By this way, this pair under neutral to alkaline pH adopts the deprotonated form, which exposed a new cavity through a movement of loop of residues 57-62, as we showed here by molecular dynamics simulations. This cavity emerges as a possible selectivity point of the cruzain inhibitors, once Homo sapiens enzynes does not present the aspartic acid - aspartate pair. In condition of acidic pH, principal component analysis of molecular dynamics simulations showed a gradual reduction of the conformational space covered by the active site of cruzain, cathepsin K and cathepsin L in it free form and complexed with dipeptidyl nitrilelike molecules in it noncovalent e covalent forms. This suggests these systems follows the conformational selection model of protein flexibility. Furthermore, we observed the ligands that inhibits the protein at nmol.L-1 induces the protein flexibility in a similar way, while the µmol.L-1 ones leads to another pattern. That made possible the construction of a decision tree which is able to identify nmol.L-1 from µmol.L-1 complexes.

As alquilfosfocolinas (APC) são uma classe de fármacos derivados de fosfolipídios endógenos que apresentam potencial antitumoral. Diferentemente de outros fármacos antitumorais que agem no DNA da célula, as APCs têm como primeiro local de ação a membrana plasmática e proteínas de sinalização, como a PKC. O objetivo desse trabalho é elucidar, através de metodologias computacionais, os possíveis mecanismos de ação das APCs que provocam hemólise, inibição da PKC e interação com membranas celulares. Inicialmente, a toxicidade de um conjunto de 34 APCs foi estudada pelos métodos quimiométricos de Análise de Agrupamentos Hierárquicos (HCA) e Componentes Principais (PCA). As moléculas foram simuladas com dinâmica molecular (DM) e propriedades físico-químicas e estruturais foram calculadas para os confôrmeros de menor energia. Após aplicação de HCA e PCA, as APCs foram divididas em 3 grupos, de acordo com suas características estruturais. Os resultados sugerem que a presença de grupos catiônicos volumosos, ou anéis como adamantila e ciclohexila, aumentam a hemólise de compostos de cadeia alquílica longa. Anéis macrocíclicos como ciclopentadecila parecem ser importantes para o potencial hemolítico de compostos com cadeia alquílica curta. Com relação a compostos sem anéis e de cadeia linear, grupos catiônicos menos volumosos parecem favorecer a hemólise. Na próxima etapa do estudo, 7 derivados de APC, com diferentes grupos catiônicos, foram selecionados e ancorados no domínio C2 da PKCα. O intuito foi mapear resíduos de aminoácidos importantes para a interação dos ligantes com a enzima, e comparar com o modo de ligação do ativador endógeno fosfatidilserina (PS). Mais uma vez, HCA e PCA foram aplicados para extrair informação relevante do mapeamento. Os resultados mostraram que as cadeias laterais de Pro188, Asn189, Arg216, Trp247, Asp249 e Thr250 não permitem a aproximação adequada do ligante, o que impede que a porção fosforila se coordene com um dos átomos de cálcio. A porção catiônica da PS, em contrapartida, consegue estabelecer ligação-hidrogênio com Asn189 de forma a posicionar os oxigênios da fosforila para interagir, ao mesmo tempo, com o átomo de cálcio. Com menos pontos de coordenação, a afinidade de ligação do cálcio pela PKCα diminui e a ativação da enzima fica comprometida, interrompendo toda a cascata de sinalização que depende dela. A parte final desse trabalho se dedicou ao estudo da interação da miltefosina com diferentes bicamadas lipídicas sob o ponto de vista termodinâmico. Oito bicamadas de diferentes fosfolipídios foram simuladas por DM e a interação energética da miltefosina foi calculada por Umbrella Sampling. Os resultados mostraram que a miltefosina apresenta maior partição em bicamadas contendo colesterol, sendo a miscibilidade nesses sistemas cerca de 76 vezes maior que os valores encontrados para bicamadas sem colesterol. Além disso, verificou-se que a internalização da miltefosina é mais fácil em regiões contendo lipídeos poli-insaturados, provavelmente devido ao empacotamento mais frouxo da bicamada. Os dados sugerem que a miltefosina age principalmente em rafts lipídicos e que células contendo mais lipídicos poli-insaturados podem incorporar maior quantidade do fármaco. / Alquilfosfocolines (APCs) comprise a class of drugs with antitumor activity derived from endogenous phospholipids. Differently from other drugs whose primary site of action is the DNA, APCs act firstly in the plasma membrane and signaling proteins, such as PKC. The main objective of this work is to elucidate, via computational approaches, the possible mechanisms of actions that cause hemolysis, PKC inhibition and interaction with cellular membranes. Initially, a set of 34 APCs was studied by means of Hierarchical Cluster Analysis (HCA) and Principal Component Analysis (PCA). The molecules were simulated by means of molecular dynamics simulations (MD) and molecular and structural properties were calculated for the lowest-energy conformer. After HCA and PCA methodologies, the set was divided into 3 groups according to their structural features. The findings suggest that the presence of bulky cationic moieties, or the adamantyl and cyclohexil rings, increase the hemolytic potential of compounds with long alkyl chains. Macrocyclic rings, such as cyclopentadecyl, seem to be important to elevate the hemolysis of compounds with short alkyl chains. Regarding linear carbon chain derivatives with no ring substitution, less bulky cationic head groups seem to favor hemolysis. In the next step of this work, 7 APC derivatives were selected and docked in the C2 domain of PKCα. The aim now was to map the residues relevant for ligands interaction compared to the binding mode of the endogenous activator, phosphatidylserine (PS). HCA and PCA were again applied in order to extract relevant information from the mapping. The results showed that the lateral chains of Pro188, Asn189, Arg216, Trp247, Asp249 and Thr250 do not allow the proper approximation of the ligands, impeding the phosphoryl moiety from coordinating with one of the calcium atoms. On the other hand, the cationic moiety of PS forms hydrogen-bonding with Asn189 in order to position the oxygens to interact, at the same time, with a calcium atom. With less coordination sites, calcium binding affinity diminishes and the enzyme activation is compromised, interrupting the signaling cascade. The final part of this work was dedicated to the study of miltefosine interaction with different lipid bilayers from the thermodynamics standpoint. Eight bilayers were simulated with MD and the energetic interaction was calculated via Umbrella Sampling simulations. The findings showed that miltefosine has higher partition in bilayers containing cholesterol, with miscibility of about 76 times higher than the values referring to bilayers without cholesterol. Moreover, it was observed that the internalization of miltefosine is facilitated in regions containing polyunsaturated lipids, probably due to the looser packing. The data suggest that miltefosine acts primarily in lipid rafts, and that cells containing more polyunsaturated lipids in their membranes can incorporate higher quantities of this drug.

Proteins and DNA are essential to life as we know it and understanding their function is understanding their structure and dynamics. The importance of the latter is being appreciated more in recent years and has led to the development of novel interdisciplinary techniques and approaches to studying protein function. Three techniques to study protein structure and dynamics have been used and combined in different ways in the context of this thesis and have led to a better understanding of the three systems described herein. X-ray crystallography is the oldest and still arguably most popular technique to study macromolecular structures. Nuclear magnetic resonance (NMR) spectroscopy is a not much younger technique that is a powerful tool not only to probe molecular structure but also dynamics. The last technique described herein are molecular dynamics (MD) simulations, which are only just growing out of their infancy. MD simulations are computer simulations of macromolecules based on structures solved by X-ray crystallography or NMR spectroscopy, that can give mechanistic insight into dynamic processes of macromolecules whose amplitudes can be estimated by the former two techniques. MD simulations of the model protein GB3 (B3 immunoglobulin-binding domain of streptococcal protein G) were conducted to identify origins of discrepancies between order parameters derived from different sets of MD simulations and NMR relaxation experiments.The results highlight the importance of time scales as well as sampling when comparing MD simulations to NMR experiments. Discrepancies are seen for unstructured regions like loops and termini and often correspond to nanosecond time scale transitions between conformational substates that are either over- or undersampled in simulation. Sampling biases can be somewhat remedied by running longer (microsecond time scale) simulations. However, some discrepancies persist over even very long trajectories. We show that these discrepancies can be due to the choice of the starting structure and more specifically even differences in protonation procedures. A test for convergence on the nanosecond time scale is shown to be able to correct for many of the observed discrepancies. Next, MD simulations were used to predict in vitro thermostability of members of the bacterial Ribonuclease HI (RNase H) family of endonucleases. Thermodynamic stability is a central requirement for protein function and a goal of protein engineering is improvement of stability, particularly for applications in biotechnology. The temperature dependence of the generalized order parameter, S, for four RNase H homologs, from psychrotrophic, mesophilic and thermophilic organisms, is highly correlated with experimentally determined melting temperatures and with calculated free energies of folding at the midpoint temperature of the simulations. This study provides an approach for in silico mutational screens to improve thermostability of biologically and industrially relevant enzymes. Lastly, we used a combination of X-ray crystallography, NMR spectroscopy and MD simulations to study specificity of the interaction between Drosophila Hox proteins and their DNA target sites. Hox proteins are transcription factors specifying segment identity during embryogenesis of bilaterian animals. The DNA binding homeodomains have been shown to confer specificity to the different Hox paralogs, while being very similar in sequence and structure. Our results underline earlier findings about the importance of the N-terminal arm and linker region of Hox homeodomains, the cofactor Exd, as well as DNA shape, for specificity. A comparison of predicted DNA shapes based on sequence alone with the shapes observed for different DNA target sequences in four crystal structures when in complex with the Drosophila Hox protein AbdB and the cofactor Exd, shows that a combined ”induced fit”/”conformational selection” mechanism is the most likely mechanism by which Hox homeodomains recognize DNA shape and achieve specificity. The minor groove widths for all sequences is close to identical for all ternary complexes found in the different crystal structures, whereas predicted shapes vary between the different DNA sequences. The sequences that have shown higher affinity to AbdB in vitro have a predicted DNA shape that matches the observed DNA shape in the ternary complexes more closely than the sequences that show low in vitro affinity to AbdB. This strongly suggests that the AbdB-Exd complex selects DNA sequences with a higher propensity to adopt the final shape in their unbound form, leading to higher affinity. An additional AbdB monomer binding site with a strongly preformed binding competent shape is observed for one of the oligomers in the reverse complement strand of one of the canonical (weak) Hox-Exd complex binding site. The shape preference seems strong enough for AbdB monomer binding to compete with AbdB-Exd dimer binding to that same oligomer, suggested by the presence of both binding modes in the same crystal. The monomer binding site is essentially able to compete with the dimer binding site, even though binding with the cofactor is not possible, because its shape is very close to the ideal shape. A comparison of different crystal structures solved herein and in the literature as well as a set of molecular dynamics simulations was performed and led to insights about the importance of residues in the Hox N-terminal arm for the preference of certain Hox paralogs to certain DNA shapes. Taken together all these insights contribute to our understanding of Hox specificity in particular as well as protein-DNA interactions in general.

Nowadays, increasing pressure from legislation and customer demands in the automotive industry are forcing manufacturers to produce greener vehicles with lower emissions and fuel consumption.As a result, electrified and hybrid vehicles are a growing popular alternative to traditional internal combustion engines (ICE). The noise from an electric vehicle comes mainly from contact between tyres and road, wind resistance and driveline. The noise emitted from the driveline is for the mostpart related to the gearbox. When developing a driveline, it is a factor of importance to estimate the noise radiating from the gearbox to achieve an acceptable design.Gears are used extensively in the driveline of electric vehicles. As the gears are in mesh, a main intrusive concern is known as gear whine noise. Gear whine noise is an undesired vibroacoustic phenomenon and is likely to originate through the gear contacts and be transferred through themechanical components to the housing where the vibrations are converted into airborne and structure-borne noise. The gear whine noise originates primarily from the excitation coming from transmission error (TE). Transmission error is defined as the difference between the ideal smoothtransfer of motion of a gear and what is in practice due to lack of smoothness.The main objective of this study is to simulate the vibrations generated by the gear whine noise in an electric powertrain line developed by AVL Vicura. The electric transmission used in this study provides only a fixed overall gear ratio, i.e. 9.59, under all operation conditions. It is assumed thatthe system is excited only by the transmission error and the mesh stiffness of the gear contacts. In order to perform NVH analysis under different operating conditions, a multibody dynamics model according to the AVL Excite program has been developed. The dynamic simulations are thencompared with previous experimental measurements provided by AVL Vicura.Two validation criteria have been used to analyse the dynamic behaviour of the AVL Excite model: signal processing using the FFT method and comparison with the experimental measurements.The results from the AVL Excite model show that the FFT criterion is quite successful and all excitation frequencies are properly observed in FFT plots. Nevertheless, when it comes to the second criterion, as long as not all dynamic parameters of the system such as damping or stiffnesscoefficients are provided with certainty in the model, it is too difficult to investigate the accuracy of the AVL Excite model. Another investigation is a numerical design study to analyses how the damping coefficients influence the response. After reducing the damping parameters, the results show that the housing and bearings have the highest influence on the response. If more acceptable results are desired,future studies must be concentrated on these to obtain more acceptable damping values. / För närvarande tvingar ökat tryck från lagstiftning och kundkrav inom bilindustrin tillverkarna attproducera grönare fordon med lägre utsläpp och bränsleförbrukning. Som ett resultat ärelektrifierade och hybridfordon ett växande populärt alternativ till traditionellaförbränningsmotorer (ICE). Bullret från ett elfordon kommer främst från kontakten mellan däckoch väg, vindmotstånd och drivlinan. Bullret från drivlinan är i huvudsak relaterat till växellådan.Vid utveckling av en drivlina är det av betydelse att uppskatta bullret från växellådan för att uppnåen acceptabel design.Utväxlingar används i stor utsträckning i elfordons drivlina. Eftersom kugghjulen är i kontaktuppstår ett huvudproblem som är känt som ett vinande ljud från kugghjulskontakten.Kugghjulsljud är ett oönskat vibro-akustiskt fenomen och uppstår sannolikt på grund avkugghjulkontakterna och överförs via de mekaniska komponenterna till växellådshuset därvibrationerna omvandlas till luftburet och strukturburet ljud. Kugghjulsljudet härstammarhuvudsakligen från exciteringen som kommer från transmissionsfel (TE) i kugghjulskontakten.Överföringsfelet definieras som skillnaden mellan den ideala smidiga rörelseöverföringen hoskugghjulen och rörelsen som sker i verkligheten på grund av ojämnheter.Huvudsyftet med denna studie är att simulera vibrationerna som genereras avkugghjulskontakterna i en elektrisk drivlina utvecklad av AVL Vicura. Den elektriska drivlinan somanvänds i denna studie har endast ett fast utväxlingsförhållande, dvs 9,59, för alladriftsförhållanden. Det antas att systemet är exciterat endast av överföringsfelet och kugghjulensstyvhet i kuggkontakterna. För att kunna utföra NVH-analys under olika driftsförhållanden har enstelkroppsdynamikmodell utvecklats med hjälp av programmet AVL Excite. De dynamiskasimuleringarna jämförs sedan med tidigare experimentella mätningar som tillhandahålls av AVLVicura.Två valideringskriterier har använts för att analysera det dynamiska beteendet hos AVL Excitemodellen:signalbehandling med FFT-metoden och jämförelse med experimentella mätningar.Resultaten från AVL Excite-modellen visar att FFT-kriteriet är ganska framgångsrikt och allaexcitationsfrekvenser observeras korrekt i FFT-diagrammen. Men när det gäller det andra kriteriet,så länge som inte alla dynamiska parametrar i systemet, såsom dämpnings- ellerstyvhetskoefficienter, är tillförlitliga i modellen, är det för svårt att undersöka exaktheten hos AVLExcite-modellen.En annan undersökning som utförts är en numerisk designstudie för att analysera hurdämpningskoefficienterna påverkar responsen. Efter minskning av dämpningsparametrarna visarresultaten att växellådshus och lager har störst inflytande på resultatet. Om mer acceptabla resultatär önskvärda måste framtida studier koncentreras på dessa parametrar för att uppnå mer acceptabladämpningsvärden.

Este trabalho introduz um modelo fisiológico para o reflexo pupilar em função das condições de iluminação (Pupil Light Reflex - PLR), e um modelo baseado em imagem para deformação dos padrões da íris. O modelo para PLR expressa o diâmetro da pupila ao longo do tempo e em função da iluminação ambiental, sendo descrito por uma equação diferencial com atraso, adaptando naturalmente o tamanho da pupila a mudanças bruscas de iluminação. Como os parâmetros do nosso modelo são derivados a partir de modelos baseados em experimentos científicos, ele simula corretamente o comportamento da pupila humana para um indivíduo médio. O modelo é então estendido para dar suporte a diferenças individuais e a hippus, além de utilizar modelos para latência e velocidade de dilatação e contração. Outra contribuição deste trabalho é um modelo para deformação realista dos padrões da íris em função da contração e dilatação da pupila. Após capturar várias imagens de íris de diversos voluntários durante diferentes estágios de dilatação, as trajetórias das estruturas das íris foram mapeadas e foi identificado um comportamento médio para as mesmas. Demonstramos a eficácia e qualidade dos resultados obtidos, comparando-os com fotografias e vídeos capturados de íris reais. Os modelos aqui apresentados produzem efeitos foto-realistas e podem ser utilizados para produzir animações preditivas da pupila e da íris em tempo real, na presença de variações na iluminação. Combinados, os dois modelos permitem elevar a qualidade de animações faciais, mais especificamente, animações da íris humana. / This thesis introduces a physiologically-based model for pupil light reflex (PLR) and an image-based model for iridal pattern deformation. The PLR model expresses the pupil diameter as a function of the environment lighting, naturally adapting the pupil diameter even to abrupt changes in light conditions. Since the parameters of the PLR model were derived from measured data, it correctly simulates the actual behavior of the human pupil. The model is extended to include latency, constriction and dilation velocities, individual differences and some constrained random noise to model hippus. The predictability and quality of the simulations were validated through comparisons of modeled results against measured data derived from experiments also described in this work. Another contribution is a model for realist deformation of the iris pattern as a function of pupil dilation and constriction. The salient features of the iris are tracked in photographs, taken from several volunteers during an induced pupil-dilation process, and an average behavior of the iridal features is defined. The effectiveness and quality of the results are demonstrated by comparing the renderings produced by the models with photographs and videos captured from real irises. The resulting models produce high-fidelity appearance effects and can be used to produce real-time predictive animations of the pupil and iris under variable lighting conditions. Combined, the proposed models can bring facial animation to new photorealistic standards.

Yan, Changqing. / Thesis (M.Phil.)--Chinese University of Hong Kong, 2009. / Includes bibliographical references (leaves 92-97). / Abstracts in English and Chinese. / ABSTRACT (ENGLISH) --- p.iii / ABSTRACT (CHINESE) --- p.iv / ACKNOWLEDGMENTS --- p.v / LIST OF ABBREVIATIONS --- p.vi / TABLE OF CONTENTS --- p.vii / LIST OF FIGURES --- p.ix / LIST OF TABLES --- p.x / Chapter CHAPTER ONE: --- BACKGROUND INFORMATION --- p.1 / Chapter 1.1 --- Introduction --- p.1 / Chapter 1.2 --- Deoxyribonucleic Acid --- p.2 / Chapter 1.3 --- DNA Studies --- p.9 / Chapter 1.4 --- Cisplatin Studies --- p.11 / Chapter 1.5 --- Scope of the Thesis --- p.13 / Chapter CHAPTER TWO: --- METHODOLOY AND COMPUTATION --- p.16 / Chapter 2.1 --- Introduction --- p.16 / Chapter 2.2 --- Molecular Dynamics Simulation --- p.16 / Chapter 2.3 --- Quantum Mechanics Calculation --- p.23 / Chapter 2.4 --- Verification of Methodology --- p.25 / Chapter 2.4.1 --- Backbone Torsion Angles --- p.25 / Chapter 2.4.2 --- N7-N7 Distance --- p.30 / Chapter 2.4.3 --- Location of HOMO --- p.33 / Chapter 2.5 --- Summary --- p.35 / Chapter CHAPTER THREE: --- UNDERSTANDING OF THE CISPLATIN-DNA CROSSLINKS --- p.36 / Chapter 3.1 --- Introduction --- p.36 / Chapter 3.2 --- MO Analysis --- p.37 / Chapter 3.3 --- Potential Binding Products with the Ligand --- p.37 / Chapter 3.3.1 --- "1,2-d(GpG) Intrastrand Crosslink" --- p.43 / Chapter 3.3.2 --- "l,2-d(ApG) Intrastrand Crosslink" --- p.43 / Chapter 3.3.3 --- "l,3-d(GpXpG) Intrastrand Crosslink" --- p.44 / Chapter 3.3.4 --- d(GpC)d(GpC) Interstrand Crosslink --- p.44 / Chapter 3.3.5 --- d(GpXpC)d(GpXpC) Interstrand Crosslink --- p.44 / Chapter 3.3.6 --- Summary --- p.45 / Chapter 3.4 --- Potential Binding Products Analysis --- p.47 / Chapter 3.4.1 --- Site Identification Convention --- p.47 / Chapter 3.4.2 --- Potential Binding Products Analysis --- p.48 / Chapter 3.4.3 --- Applications --- p.53 / Chapter 3.5 --- Cisplatin-DNA Crosslink Products Analysis --- p.56 / Chapter 3.5.1 --- "1,2-d(GpG) and l,2-d(ApG) Intrastrand Crosslinks" --- p.61 / Chapter 3.5.2 --- "l,3-d(GpXpG) Intrastrand and d(GpXpC)d(GpXpC) Interstrand Crosslinks" --- p.62 / Chapter 3.5.3 --- d(GpC)d(GpC) Interstrand Crosslinks --- p.63 / Chapter 3.5.4 --- Platination at Terminal Positions --- p.65 / Chapter 3.6 --- Summary --- p.65 / Chapter CAHPTER FOUR: --- CONCLUDING REMARKS --- p.67 / APPENDIX I: BACKBONE TORSION ANGLES AND SUGAR RING CONFORMATIONS OF THE OPTIMIZED GEOMETRIES --- p.69 / APPENDIX II: BACKBONE TORSION ANGLES OF THE EXPERIMENTAL SEQUENCES FROM NUCLEIC ACID DATABASE (NDB) --- p.77 / REFERENCES --- p.92

The conventional Verlet table neighbor list algorithm is improved to reduce the number of unnecessary inter-atomic distance calculations in molecular simulations involving large amount of atoms. Both of the serial and parallelized performance of molecular dynamics simulation are evaluated using the new algorithm and compared with those using the conventional Verlet table and cell-linked list algorithm. Results show that the new algorithm significantly improved the performance of molecular dynamics simulation compared with conventional neighbor list maintaining and utilizing algorithms in serial programs as well as parallelized programs. / Singapore-MIT Alliance (SMA)

Microcantilever sensor has gained much popularity because of its high sensitivity and selectivity. It consists of a micro-sized cantilever that is usually coated on one side with chemical/biological probe agents to generate strong attraction to target molecules. The interactions between the probe and target molecules induce surface stress that bends the microcantilever. This current work applied the molecular dynamics simulation to study the microcantilever system. Lennard-Jones potentials were used to model the target-target and target-probe interactions and bond bending potentials to model the solid cantilever beam. In addition, this work studied the effect of probe locations on the microcantilever deflection. The simulation results suggest that both target-target and target-probe interactions as well as the probe locations affect the arrangement of the bonds; in term of the bonding number, the area containing the bonded molecules, and the distances between them. All these factors influence the microcantilever deflection.

Carbohydrates are abundant in nature and have functions ranging from energy storage to acting as structural components. Analysis of carbohydrate structures is important and can be used for, for instance, clinical diagnosis of diseases as well as in bacterial studies. The complexity of glycans makes it difficult to determine their structures. NMR spectroscopy is an advanced method that can be used to examine carbohydrates at the atomic level, but full assignments of the signals require much work. Reliable automation of this process would be of great help. Herein studies of Escherichia coli O-antigen polysaccharides are presented, both a structure determination by NMR and also research on glycosyltransferases which assemble the polysaccharides. The computer program CASPER has been improved to assist in carbohydrate studies and in the long run make it possible to automatically determine structures based only on NMR data. Detailed computer studies of glycans can shed light on their interactions with proteins and help find inhibitors to prevent unwanted binding. The WaaG glycosyltransferase is important for the formation of E. coli lipopolysaccharides. Molecular docking analyses of structures confirmed to bind this enzyme have provided information on how inhibitors could be composed. Noroviruses cause gastroenteritis, such as the winter vomiting disease, after binding human histo-blood group antigens. In one of the projects, fragment-based docking, followed by molecular dynamics simulations and binding free energy calculations, was used to find competitive binders to the P domain of the capsid of the norovirus VA387. These novel structures have high affinity and are a very good starting point for developing drugs against noroviruses. The protein targets in these two projects are carbohydrate binding, but the techniques are general and can be applied to other research projects. / At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 4: Submitted. Paper 5: Manuscript. Paper 6. Manuscript.