No description available.

Doctor of Philosophy / Department of Chemistry / Paul E. Smith / Molecular Dynamics (MD) simulations have played an important role in providing detailed atomic information for the study of biological systems. The quality of an MD simulation depends on both the degree of sampling and the accuracy of force field. Kirkwood-Buff (KB) theory provides a relationship between species distributions from simulation results and thermodynamic properties from experiments. Recently, it has been used to develop new, hopefully improved, force fields and to study preferential interactions. Here we combine KB theory and MD simulations to study a variety of intermolecular interactions in solution. Firstly, we present a force field for neutral amines and carboxylic acids. The parameters were developed to reproduce the composition dependent KB integrals obtained from an analysis of the experimental data, allowing for accurate descriptions of activities involved with uncharged N-terminus and lysine residues, as well as the protonated states for the C-terminus and both aspartic and glutamic acids. Secondly, the KB force fields and KB theory are used to investigate the urea cosolvent effect on peptide aggregation behavior by molecular dynamics simulation. Neo-pentane, benzene, glycine and methanol are selected to represent different characteristics of proteins. The chemical potential derivatives with respect to the cosolvent concentrations are calculated and analyzed, and the four solutes exhibit large differences. Finally, the contributions from the vibrational partition function to the total free energy and enthalpy changes are investigated for several systems and processes including: the enthalpy of evaporation, the free energy of solvation, the activity of a solute in solution, protein folding, and the enthalpy of mixing. The vibrational frequencies of N-methylacetamide, acetone and water are calculated using density functional theory and MD simulations. We argue that the contributions from the vibrational partition function are large and in classical force fields these contributions should be implicitly included by the use of effective intermolecular interactions.

Master of Science / Department of Biochemistry / Paul E. Smith / We are developing a force field (FF) for molecular dynamics (MD) simulations of peptides and small proteins that is grounded in the Kirkwood-Buff theory of solutions. Here we present the Kirkwood-Buff Force Field (KBFF) parameters for the aromatic amino acids, based upon simulations of binary mixtures of small molecules representative of these amino acids over their entire composition ranges (excluding Histidine). Many aromatics are not fully soluble in water, so they have instead been studied in solvents of methanol or toluene. The parameters were developed by studying the following binary solutions: Phenylalanine − benzene + methanol, toluene + methanol, and toluene + benzene; Tyrosine − toluene + phenol and toluene + p-Cresol; Tryptophan − pyrrole + methanol and indole + methanol; Histidine − pyrrole + methanol, pyridine + methanol, pyridine + water, histidine + water (at 0.25 molal), and histidine monohydrochloride + water (at 0.3 molal and 0.6 molal). Our simulations reproduce the Kirkwood-Buff integrals, which guarantees that the KBFF provides an adequate balance of solute-solvent, solute-solute, and solvent-solvent interactions. Additionally, we show that the KBFF does not sacrifice reproduction of other solution properties in order to achieve this improved description of intermolecular interactions. We present these results as validating evidence for the future use of the KBFF in simulations of peptides and small proteins.

Doctor of Philosophy / Department of Chemistry / Paul E. Smith / The force fields used in computer simulations play an important role in describing a particular system. In order to estimate the accuracy of a force field, physical or thermodynamic properties are usually compared with simulation results. Recently, we have been developing a force field which is called the Kirkwood-Buff Force Field (KBFF). This force field is established by transforming experimental data into Kirkwood-Buff (KB) integrals and then attempting to reproduce those KBIs with molecular dynamic (MD) simulations. Here we investigate a variety of intermolecular interactions in aqueous solutions through KB theory and molecular simulations. First, we describe a force field for the simulation of alkali halide aqueous solutions. These models are developed specifically to reproduce the experimentally determined Kirkwood-Buff integrals and solution activities as a function of molality. Additionally, other experimentally known properties including ion diffusion constants, relative permittivities, the densities and heats of mixing are reproduced by these models. Second, In an effort to understand the interactions which occur between amino acids in solution we have developed new force fields for simple amino acids and their analogs including glycine, betaine, β-alanine, dl-alanine, NH4Cl, NH4Br, N(CH3)4Cl, N(CH3)4Br, CH3NH3Cl, and CH3COONa. The new force fields reproduce the experimental Kirkwood-Buff integrals which describe the relative distribution of all the species in a solution mixture. Furthermore, it is shown that these simple amino acids can be understood in terms of the interactions of their functional groups and that, to a very good approximation, the transferability and additivity usually assumed in the development of biomolecular force fields appear to hold true. Third, an analysis of the effect of a cosolvent on the association of a solute in solution is presented by using the Kirkwood-Buff theory of solutions. The derived expressions provide a foundation for the investigation of cosolvent effects on molecular and biomolecular equilibria, including protein association, aggregation, and cellular crowding. Finally, in an effort to understand peptide aggregation at the atomic level we have performed simulations of polyglycine ((gly)n) using our recently developed force fields. Experimentally, the association of glycine polypeptides increases with n. Our force fields reproduce this behavior, and we investigated the reasons behind this trend. In addition to studying closed ensembles, we also simulate these systems in a semi-open ensemble that was designed to mimic cellular environments typically open to water, using a simple direct approach. The differences between the two ensembles are investigated and compared with our recent theoretical descriptions of aggregating systems using Kirkwood-Buff theory.

By developing strong, self-sustaining individuals in organisations, employees will need to withstand the forces of stressors and be able to cope with continuous changes within their working environment. To move the equilibrium level from the current to the desired condition, the field of forces must be altered - by adding driving forces and by removing restraining forces. An individual's sense of coherence may either alleviate or aggravate reactions to a stressor and moderate the impact of occupational stressors on the individual's affective outcomes. The objectives of this study were to validate the 13-item version of the Orientation to Life Questionnaire (OLQ) (Antonovsky, 1987) and to determine the experience of employees with high levels of sense of coherence regarding helping and restraining factors within the workplace (compared to that of those with a low sense of coherence). A cross-sectional s w e y design was used. The total population (N = 2 678) of employees in a financial institution in Gauteng was used in this study. Random samples of groups with a strong (n = 300) and low (n = 300) sense of coherence were taken for purposes of the qualitative study. The OLQ (Antonovsky, 1987) and a biographical questionnaire were administered. The scale showed acceptable reliability and construct validity. The study set out to determine the applicability of the theoretical model of sense of coherence to employees in a financial institution. Reliability analysis revealed that the three subscales of sense of coherence were sufficiently internally consistent. The reliability of the measuring instrument were assessed with the use of Cronbach alpha coefficients. Descriptive statistics (e.g. mean and standard deviations) were used to analyse the data. By using the structural equation modelling approach, a one-dimensional factor structure for sense of coherence amongst employees in a financial institution emerged. The reliability analysis revealed that the alpha scores were acceptable. It can therefore be concluded that the 13-item version of the OLQ is a reliable and valid measuring instrument. The results showed that although employee groups with high and low levels of sense of coherence are experiencing similar helping and restraining factors within their working environment, helping factors are being experienced with a higher frequency by groups with high levels of sense of coherence and restraining factors with a higher frequency by groups with low levels of sense of coherence. Recommendations were made for future research. / Thesis (M.A. (Industrial Psychology))--North-West University, Potchefstroom Campus, 2007

Doctor of Philosophy / Department of Chemistry / Paul E. Smith / The impact of computer simulations has become quite significant especially with the development of supercomputers during the last couple of decades. They are used in a wide range of purposes such as exploring experimentally inaccessible phenomena and providing an alternative when experiments are expensive, dangerous, time consuming, difficult and controversial. In terms of applications in biological systems molecular modeling techniques can be used in rational drug design, predicting structures of proteins and circumstances where the atomic level descriptions provided by them are valuable for the understanding of the systems of interest. Hence, the potential of computer simulations of biomolecular systems is undeniable. Irrespective of the promising uses of computer simulations, it cannot be guaranteed that the results will be realistic. The precision of a molecular simulation depends on the degree of sampling achieved during the simulation while the accuracy of the results depends on the satisfactory description of intramolecular and intermolecular interactions in the system, i.e. the force field. Recently, we have been developing a force field for molecular dynamics simulations of biological systems based on the Kirkwood Buff (KB) theory of solutions, not only with an emphasis on the accurate description of intermolecular interactions, but also by reproducing several physical properties such as partial molar volume, compressibility and composition dependent chemical potential derivatives to match with respective experimental values. In this approach simulation results in terms of KB integrals can be directly compared with experimental data through a KB analysis of the solution properties and therefore it provides a simple and clear method to test the capability of the KB derived force field. Initially, we have provided a rigorous framework for the analysis of experimental and simulation data concerning open and closed multicomponent systems using the KB theory of solutions. The results are illustrated using computer simulations for various concentrations of the solutes Gly, Gly₂ and Gly₃ in both open and closed systems, and in the absence or presence of NaCl as a cosolvent. Then, we have attempted to quantify the interactions between amino acids in aqueous solutions using the KB theory of solutions. The results are illustrated using computer simulations for various concentrations of the twenty zwitterionic amino acids at ambient temperature and pressure. Next, several amino acids were also studied at higher temperatures and pressures and the results are discussed in terms of the preferential (solute over solvent) interactions between the amino acids. Finally, we have described our most recent efforts towards a complete force field for peptides and proteins. The results are illustrated using molecular dynamics simulations of several tripeptides, selected peptides and selected globular proteins at ambient temperature and pressure followed by replica exchange molecular dynamics simulations of a few selected peptides.

One of the central tasks for biomolecular modeling is to develop accurate and computationally cheap methods. In this dissertation, we present the development of a brand new polarizable force field—Polarizable Simulations with Second order Interaction Model (POSSIM) involving electrostatic polarization. The POSSIM framework combines accuracy of a polarizable force field and computational efficiency of the second-order approximation of the full-scale induced point dipole polarization formalism. POSSIM force field has been extended to include parameters for small molecules serving as models for peptide and protein side-chains. Parameters have been fitted to permit reproducing many-body energies, gas-phase dimerization energies and geometries and liquid-phase heats of vaporization and densities. Quantum mechanical and experimental data have been used as the target for the fitting. The resulting parameters can be used for simulations of the parameterized molecules themselves or their analogues. In addition to this, these force field parameters have been employed in further development of the POSSIM fast polarizable force field for proteins. The POSSIM framework has been expanded to include a complete polarizable force field for proteins. Most of the parameter fitting was done to high-level quantum mechanical data. Conformational geometries and energies for dipeptides have been reproduced within average errors of ca. 0.5 kcal/mol for energies of the conformers (for the electrostatically neutral residues) and 9.7º for key dihedral angles. We have also validated this force field by simulating an elastin-like polypeptide GVG(VPGVG)3 in aqueous solution. Elastin-like peptides with the (VPGVG)n motif are known to exhibit anomalous behavior of their radius of gyration that increases when temperature is lowered (the so called inverse temperature transition). We have simulated the system with the OPLS-AA and POSSIM force fields and demonstrated that our newly developed polarizable POSSIM parameters permit to capture the experimentally observed decrease of the radius of gyration with increasing temperature, while the fixed-charges OPLS-AA ones do not. Furthermore, our fitting of the force field parameters for the peptides and proteins has been streamlined compared with the previous generation of the complete polarizable force field and relied more on transferability of parameters for non-bonded interactions (including the electrostatic component). The resulting deviations from the quantum mechanical data are similar to those achieved with the previous generation, thus the technique is robust and the parameters are transferable. At the same time, the number of parameters used in this work was noticeably smaller than that of the previous generation of our complete polarizable force field for proteins, thus the transferability of this set can be expected to be greater and the danger of force field fitting artifacts is lower. Therefore, we believe that this force field can be successfully applied in a wide variety of applications to proteins and protein-ligand complexes.

Simulações de Dinâmica Molecular tornaram-se uma ferramenta indispensável no estudo de sistemas em fase condensada, incluindo sistemas líquido-Cristalinos, e na predição de propriedades dinâmicas. Cristais líquidos possuem um leque enorme de aplicações, mas o estudo teórico destes sistemas torna-se complicado devido ao seu tamanho, ao método utilizado para obtenção dos parâmetros do campo de força e, principalmente, à transferibilidade dos parâmetros para outro estado termodinâmico. Tendo isso em vista, este trabalho propõe uma metodologia para parametrizar campos de força derivados de cálculos mecânico-quânticos que possuam um grau de transferibilidade confiável. O sistema escolhido neste trabalho foi o 4’-Pentil- 4-Carbonitrila, também conhecido como 5CB, pois o mesmo já é utilizado em diversos aparelhos óptico-eletrônicos; a parametrização intramolecular foi feita com o programa JOYCE, e a intermolecular, com o programa PICKY. Os parâmetros intramoleculares obtidos mostraram uma boa descrição da geometria interna do sistema, contribuindo para a parametrização intermolecular, a qual obteve uma excelente descrição de propriedades termodinâmicas. Este trabalho corrobora para a hipótese de que campos de força derivados de cálculos mecânico- quânticos podem descrever diferentes fases termodinâmicas com um alto grau de confiabilidade. / Simulations using Molecular dynamics have become a powerful tool in the study of systems in condensed phase, including liquid-crystalline systems, and in the prediction of dynamical properties. Liquid Crystals have many applications, but the theoretical study of these systems is more complex because of their size, the method that was used in the force field parametrization and, mainly, because in most of the cases parameters cannot be transferred to another thermodynamical state. With that in mind, this work propose a methodology to parametrize force fields derived from quantummechanical calculations and which can be transferred to other thermodynamical state without losing important information. The chosen system in this work was the 4-Cyano-4’-Pentylbiphenyl, also known as 5CB, which have been used in many optical-electronic device and the intramolecular parametrization was done with the JOYCE program and the program PICKY was used in the intermolecular parametrization. The intramolecular parameters obtained show a good description of the internal geometry, contributes to the intermolecular parametrization, with we obtained an accurate description of the thermodynamical and physical chemical properties. This work corroborate to the hypothesis that force field derived from quantum-mechanical calculations can describe different thermodynamical states without losing important information.

In the automotive industry, a kind of fuel additives, known as surfactant, is used to protect metallic surfaces. Its efficiency strongly depends on factors such as temperature, solvent properties and the presence of other surfactants in the system. In this thesis, both classical and ReaxFF molecular dynamics (MD) simulations are used in studying the impacts of these factors on the adsorption of organic surfactants at the fluid-solid interface. Firstly, a classical MD simulation study of competitive adsorption is carried out on a multi-functional phenol and amine surfactant model with ethanol at the oil/iron oxide interface. As the concentration of ethanol increases, the ethanol molecules effectively compete for the adsorption sites on the iron oxide surface. This observation concurs with the experimental findings of similar oil/iron oxide systems. Unlike most MD interfacial studies, ReaxFF MD uses a fully flexible and polarizable solid surface. The second part of the thesis includes a study on the effect of polarity of organic molecules on the structure of iron oxide using ReaxFF-based MD simulations. The simulation results suggest that care must be taken when parameterising empirical and transferable force fields because the fixed charges on a solid slab may not be a perfect representation of the real system, especially when the solid is in contact with polar compounds. Lastly, but not the least, missing ReaxFF interaction parameters for Fe/N have been developed to simulate the adsorption of amine based surfactant on iron oxide. The parameterisation of the force field is done by fitting these interaction parameters to a set of quantum mechanical data involving iron-based clusters. These newly developed parameters are able to capture chemisorption and proton transfer between hexadecylamine and iron oxide.

Métodos computacionais assumiram a partir de 1980, um papel de destaque no planejamento de novos fármacos, oferecendo abordagens racionais para reduzir o grau de incerteza na geração de novos compostos bioativos. Dentre estes métodos, destacam-se aqueles dependentes de campos de força, como o atracamento e a dinâmica molecular. Infelizmente, estes métodos exigem estratégias de parametrização capazes de lidar com a diversidade química associada ao planejamento de fármacos. Os esforços atuais neste sentido são focados em fase gasosa, como no caso do GAFF e MMFF94, Assim, o presente trabalho busca explorar a capacidade de um campo de força baseado em fase condensada, o GROMOS, na reprodução de propriedades fisico-químicas de anéis aromáticos comumente encontrados em fármacos. Assim, parâmetros ligados e de Lennard-Jones do GROMOS53a6 foram utilizados para a construção das topologias destas moléculas orgânicas, enquanto novos parâmetros coulômbicos e torsionais foram gerados. Em seguida, suas propriedades físico-químicas foram simuladas e comparadas aos respectivos valores experimentais, permitindo a determinação da qualidade de cada topologia. Até o momento, 41 moléculas foram parametrizadas com sucesso, levando a erros absolutos abaixo de 15% para densidade, entalpia de vaporização e capacidade térmica isobárica. A partir desta etapa de validação, os parâmetros obtidos foram aplicados no estudo de hexafirinas sintéticas em diferentes solventes e íons, acessando com sucesso a conformação e coordenação das moléculas envolvidas. Desta forma, os dados obtidos constituem-se em um benchmark para futuros estudos baseados no campo de força GROMOS, sugerindo seu potencial para estudos de moléculas orgânicas sintéticas em fase condensada. Abre-se, ainda, a perspectiva de emprego de técnicas de dinâmica molecular, com estes parâmetros, no estudo do perfil conformacional, dinâmica e flexibilidade de fármacos em solução. / Since 1980, computational methods took a major role in drug design, offering rational approaches to mitigate the uncertainty in the generation of new bioactive compounds. Among these methods, it is possible to highlight the force field dependents, like molecular docking and molecular dynamics. Unfortunately, these methods demand parametrization strategies capable to to deal with such chemical diversity associated with drug planning. The efforts currently available in this sense are focused on gas phase, like GAFF and MMF94, Therefore, the present work aims to explore the capacity of a force field based on condensed phase, GROMOS, on reproducing physical-chemical properties of aromatic rings commonly found on drugs. Thus, bonded and Lennard-Jones parameters of GROMOS53a6 was used to build topologies for these organic molecules, while new coulombic and torsional parameters were generated. Next, their physical-chemical properties were simulated and compared to respective experimental data, allowing a quality determination of each topology. So far, 41 molecules were successfully parameterized, leading to absolute errors below 15% for density, enthalpy of vaporization and isobaric heat capacity. From this validation step, the obtained parameters were applied on the synthetic hexaphyrin studies in different solvents and ions, successfully accessing the conformation and coordination of the participating molecules. Therefore, the obtained data constitute a benchmark for future studies based on GROMOS force field, suggesting its potential to carry out studies on synthetic organic molecules in condensed phase. Yet, it opens the perspective of applying molecular dynamics techniques with these parameters on the study of conformational profile, dynamics and flexibility of drugs in solution.