Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide:can the α3 and β2m domains be neglected?

Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histo-compatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVY-DGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide–MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent-accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the α1–α2 domain, the α3 and β2m domains cannot be neglected. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1803–1813, 2004

Novos protocolos teóricos utilizados nos estudos das ligações hidrogênio-hidrogênio, na estabilidade do tetraedrano, seus derivados e das reações de diels-alder e na quantificação da afinidade de fármacos

This thesis was performed in four chapters, at the theoretical level, focused mainly on electronic density. In the first chapter, we have applied an undergraduate minicourse of Diels-Alder reaction in Federal University of Rio Grande do Norte. By using computational chemistry tools students could build the knowledge by themselves and they could associate important aspects of physical-chemistry with Organic Chemistry. In the second chapter, we studied a new type of chemical bond between a pair of identical or similar hydrogen atoms that are close to electrical neutrality, known as hydrogen-hydrogen (H-H) bond. In this study performed with complexed alkanes, provides new and important information about their stability involving this type of interaction. We show that the H-H bond playing a secondary role in the stability of branched alkanes in comparison with linear or less branched isomers. In the third chapter, we study the electronic structure and the stability of tetrahedrane, substituted tetrahedranes and silicon and germanium parents, it was evaluated the substituent effect on the carbon cage in the tetrahedrane derivatives and the results indicate that stronger electron withdrawing groups (EWG) makes the tetrahedrane cage slightly unstable while slight EWG causes a greater instability in the tetrahedrane cage. We showed that the sigma aromaticity EWG and electron donating groups (EDG) results in decrease and increase, respectively, of NICS and D3BIA aromaticity indices. In addition, another factor can be utilized to explain the stability of tetra-tert-butyltetrahedrane as well as HH bond. GVB and ADMP were also used to explain the stability effect of the substituents bonded to the carbon of the tetrahedrane cage. In the fourth chapter, we performed a theoretical investigation of the inhibitory effect of the drug abiraterone (ABE), used in the prostate cancer treatment as CYP17 inhibitor, comparing the interaction energies and electron density of the ABE with the natural substrate, pregnenolone (PREG). Molecular dynamics and docking were used to obtain the CYP1ABE and CYP17-PREG complexes. From molecular dynamics was obtained that the ABE has higher diffusion trend water CYP17 binding site compared to the PREG. With the ONIOM (B3LYP:AMBER) method, we find that the interaction electronic energy of ABE is 21.38 kcal mol-1 more stable than PREG. The results obtained by QTAIM indicate that such stability is due a higher electronic density of interactions between ABE and CYP17

Positive Encounters: An Investigative Study on the Effects of Interactions between Cations and Quantum Dots

Thesis (Ph.D.)--University of Washington, 2016-06

Efeito do substituinte “-F” nas propriedades electrónicas e estruturais de fluoroquinolonas

Este estudo incide sobre as características que a presença do ião flúor em moléculas concede. Mais concretamente em fluoroquinolonas, antibióticos que cada vez são mais utilizados. Fez-se uma analise de vários parâmetros para obtermos informação sobre a interação fármaco-receptor nas fluoroquinolonas. Sendo para isso utilizadas técnicas de caracterização química computacional para conseguirmos caracterizar eletronicamente e estruturalmente (3D) as fluoroquinolonas em complemento aos métodos semi-empíricos utilizados inicialmente. Como é sabido, a especificidade e a afinidade para o sitio alvo, é essencial para eficácia de um fármaco. As fluoroquinolonas sofreram um grande desenvolvimento desde a primeira quinolona sintetizada em 1958, sendo que desde ai foram sintetizadas inúmeros derivados da mesma. Este facto deve-se a serem facilmente manipuladas, derivando fármacos altamente potentes, espectro alargado, factores farmacocinéticos optimizados e efeitos adversos reduzidos. A grande alteração farmacológica para o aumento do interesse neste grupo, foi a substituição em C6 de um átomo de flúor em vez de um de hidrogénio. Para obtermos as informações sobre a influência do ião flúor sobre as propriedades estruturais e electrónicas das fluoroquinolonas, foi feita uma comparação entre a fluoroquinolona com flúor em C6 e com hidrogénio em C6. As quatro fluoroquinolonas presentes neste estudo foram: ciprofloxacina, moxiflocacina, sparfloxacina e pefloxacina. As informações foram obtidas por programas informáticos de mecânica quântica e molecular. Concluiu-se que a presença de substituinte flúor não modificava de forma significativa a geometria das moléculas mas sim a distribuição da carga no carbono vicinal e nos átomos em posição alfa, beta e gama relativamente a este. Esta modificação da distribuição electrónica pode condicionar a ligação do fármaco ao receptor, modificando a sua actividade farmacológica.

Multiple sulfur isotope fractionations in inorganic aqueous systems

New constraints on isotope fractionation factors in inorganic aqueous sulfur systems based on theoretical and experimental techniques relevant to studies of the sulfur cycle in modern environments and the geologic rock record are presented in this dissertation. These include theoretical estimations of equilibrium isotope fractionation factors utilizing quantum mechanical software and a water cluster model approach for aqueous sulfur compounds that span the entire range of oxidation state for sulfur. These theoretical calculations generally reproduce the available experimental determinations from the literature and provide new constraints where no others are available. These theoretical calculations illustrate in detail the relationship between sulfur bonding environment and the mass dependence associated with equilibrium isotope exchange reactions involving all four isotopes of sulfur. I additionally highlight the effect of isomers of protonated compounds (compounds with the same chemical formula but different structure, where protons are bound to either sulfur or oxygen atoms) on isotope partitioning in the sulfite (S4+) and sulfoxylate (S2+) systems, both of which are key intermediates in oxidation-reduction processes in the sulfur cycle. I demonstrate that isomers containing the highest degree of coordination around sulfur (where protonation occurs on the sulfur atom) have a strong influence on isotopic fractionation factors, and argue that isomerization phenomenon should be considered in models of the sulfur cycle. Additionally, experimental results of the reaction rates and isotope fractionations associated with the chemical oxidation of aqueous sulfide are presented. Sulfide oxidation is a major process in the global sulfur cycle due largely to the sulfide-producing activity of anaerobic microorganisms in organic-rich marine sediments. These experiments reveal relationships between isotope fractionations and reaction rate as a function of both temperature and trace metal (ferrous iron) catalysis that I interpret in the context of the complex mechanism of sulfide oxidation. I also demonstrate that sulfide oxidation is a process associated with a mass dependence that can be described as not conforming to the mass dependence typically associated with equilibrium isotope exchange. This observation has implications for the inclusion of oxidative processes in environmental- and global-scale models of the sulfur cycle based on the mass balance of all four isotopes of sulfur. The contents of this dissertation provide key reference information on isotopic fractionation factors in aqueous sulfur systems that will have far-reaching applicability to studies of the sulfur cycle in a wide variety of natural settings.

Air oxygenation chemistry of 4-TBC catalyzed by chloro bridged dinuclear copper(II) complexes of pyrazole based tridentate ligands: synthesis, structure, magnetic and computational studies

Four dinuclear bis(mu-Cl) bridged copper(II) complexes, Cu-2(mu-Cl)(2)(L-X)(2)](ClO4)(2) (L-X = N,N-bis(3,5-dimethylpyrazole-1-yl)-methyl]benzylamine with X = H(1), OMe(2), Me(3) and Cl(4)), have been synthesized and characterized by the single crystal X-ray diffraction method. In these complexes, each copper(II) center is penta-coordinated with square-pyramidal geometry. In addition to the tridentate L-X ligand, a chloride ion occupies the last position of the square plane. This chloride ion is also bonded to the neighboring Cu(II) site in its axial position forming an SP-I dinuclear Cu(II) unit that exhibits small intramolecular ferromagnetic interactions and supported by DFT calculations. The complexes 1-3 exhibit methylmonooxygenase (pMMO) behaviour and oxidise 4-tert-butylcatechol (4-TBCH2) with molecular oxygen in MeOH or MeCN to 4-tert-butyl-benzoquinone (4-TBQ), 5-methoxy-4-tert-butyl-benzoquinone (5-MeO-4-TBQ) as the major products along with 6,6'-Bu-t-biphenyl-3,4,3',4'-tetraol and others as minor products. These are further confirmed by ESI- and FAB-mass analyses. A tentative catalytic cycle has been framed based on the mass spectral analysis of the products and DFT calculations on individual intermediates that are energetically feasible.

Computational investigations on periodic trends in Group 14 chemistry

Describes the use of computer-aided molecular modelling to investigate trends in the chemistry of the Group 14 elements, namely carbon, silicon, germanium, tin and lead. The chemical behaviour of two classes of molecules containing Group 14 elements was related to trends in the fundamental properties of these elements.

Computational Methods in Biophysics and Medicinal Chemistry: Applications and Challenges

In this thesis I described the theory and application of several computational methods in solving medicinal chemistry and biophysical tasks. I pointed out to the valuable information which could be achieved by means of computer simulations and to the possibility to predict the outcome of traditional experiments. Nowadays, computer represents an invaluable tool for chemists. In particular, the main topics of my research consisted in the development of an automated docking protocol for the voltage-gated hERG potassium channel blockers, and the investigation of the catalytic mechanism of the human peptidyl-prolyl cis-trans isomerase Pin1.

Influence of reaction atmosphere (H2O, N2, H2, CO2, CO) on fluidized-bed fast pyrolysis of biomass using detailed tar vapor chemistry in computational fluid dynamics

Secondary pyrolysis in fluidized bed fast pyrolysis of biomass is the focus of this work. A novel computational fluid dynamics (CFD) model coupled with a comprehensive chemistry scheme (134 species and 4169 reactions, in CHEMKIN format) has been developed to investigate this complex phenomenon. Previous results from a transient three-dimensional model of primary pyrolysis were used for the source terms of primary products in this model. A parametric study of reaction atmospheres (H2O, N2, H2, CO2, CO) has been performed. For the N2 and H2O atmosphere, results of the model compared favorably to experimentally obtained yields after the temperature was adjusted to a value higher than that used in experiments. One notable deviation versus experiments is pyrolytic water yield and yield of higher hydrocarbons. The model suggests a not overly strong impact of the reaction atmosphere. However, both chemical and physical effects were observed. Most notably, effects could be seen on the yield of various compounds, temperature profile throughout the reactor system, residence time, radical concentration, and turbulent intensity. At the investigated temperature (873 K), turbulent intensity appeared to have the strongest influence on liquid yield. With the aid of acceleration techniques, most importantly dimension reduction, chemistry agglomeration, and in-situ tabulation, a converged solution could be obtained within a reasonable time (∼30 h). As such, a new potentially useful method has been suggested for numerical analysis of fast pyrolysis.

Computational fluid dynamics analysis of a flat plate photocatalytic reactor for storm and wastewater reuse

A computational fluid dynamics (CFD) analysis has been performed for a flat plate photocatalytic reactor using CFD code FLUENT. Under the simulated conditions (Reynolds number, Re around 2650), a detailed time accurate computation shows the different stages of flow evolution and the effects of finite length of the reactor in creating flow instability, which is important to improve the performance of the reactor for storm and wastewater reuse. The efficiency of a photocatalytic reactor for pollutant decontamination depends on reactor hydrodynamics and configurations. This study aims to investigate the role of different parameters on the optimization of the reactor design for its improved performance. In this regard, more modelling and experimental efforts are ongoing to better understand the interplay of the parameters that influence the performance of the flat plate photocatalytic reactor.

Computational design of cyclic nitroxides as efficient redox mediators for dye-sensitized solar cells

Cyclic nitroxide radicals represent promising alternatives to the iodine-based redox mediator commonly used in dye-sensitized solar cells (DSSCs). To date DSSCs with nitroxide-based redox mediators have achieved energy conversion efficiencies of just over 5 % but efficiencies of over 15 % might be achievable, given an appropriate mediator. The efficacy of the mediator depends upon two main factors: it must reversibly undergo one-electron oxidation and it must possess an oxidation potential in a range of 0.600-0.850 V (vs. a standard hydrogen electrode (SHE) in acetonitrile at 25 °C). Herein, we have examined the effect that structural modifications have on the value of the oxidation potential of cyclic nitroxides as well as the reversibility of the oxidation process. These included alterations to the N-containing skeleton (pyrrolidine, piperidine, isoindoline, azaphenalene, etc.), as well as the introduction of different substituents (alkyl-, methoxy-, amino-, carboxy-, etc.) to the ring. Standard oxidation potentials were calculated using high-level ab initio methodology that was demonstrated to be very accurate (with a mean absolute deviation from experimental values of only 16 mV). An optimal value of 1.45 for the electrostatic scaling factor for UAKS radii in acetonitrile solution was obtained. Established trends in the values of oxidation potentials were used to guide molecular design of stable nitroxides with desired E° ox and a number of compounds were suggested for potential use as enhanced redox mediators in DSSCs. Copyright © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Computational fluid dynamic analysis of intracranial aneurysm bleb Formation (NS26P)

Purpose: The management of unruptured aneurysms remains controversial as treatment infers potential significant risk to the currently well patient. The decision to treat is based upon aneurysm location, size and abnormal morphology (e.g. bleb formation). A method to predict bleb formation would thus help stratify patient treatment. Our study aims to investigate possible associations between intra-aneurysmal flow dynamics and bleb formation within intracranial aneurysms. Competing theories on aetiology appear in the literature. Our purpose is to further clarify this issue. Methodology: We recruited data from 3D rotational angiograms (3DRA) of 30 patients with cerebral aneurysms and bleb formation. Models representing aneurysms pre-bleb formation were reconstructed by digitally removing the bleb, then computational fluid dynamics simulations were run on both pre and post bleb models. Pulsatile flow conditions and standard boundary conditions were imposed. Results: Aneurysmal flow structure, impingement regions, wall shear stress magnitude and gradients were produced for all models. Correlation of these parameters with bleb formation was sought. Certain CFD parameters show significant inter patient variability, making statistically significant correlation difficult on the partial data subset obtained currently. Conclusion: CFD models are readily producible from 3DRA data. Preliminary results indicate bleb formation appears to be related to regions of high wall shear stress and direct impingement regions of the aneurysm wall.