Yes / A reliable description of ion pair interactions for biological systems, particularly those involving polyatomic ions such as carboxylate and divalent ions such as Ca2+, using biomolecular force-fields is essential for making useful predictions for a range of protein functions. In particular, the interaction of divalent ions with the double carboxylate group present in γ-carboxyglutamic acid (Gla), relevant to the function of many proteins, is relatively understudied using biomolecular force-fields. Using force-field based metadynamics simulations to predict the free energy of binding between Ca2+ and the carboxylate group in liquid water, we show that a widely-used biomolecular force-field, CHARMM22*, substantially over-estimates the binding strength between Ca2+ and the side-chains of both glutamic acid (Glu) and Gla, compared with experimental data obtained for the analogous systems of aqueous calcium–acetate and calcium–malonate. To correct for this, we propose and test a range of modifications to the σ value of the heteroatomic Lennard–Jones interaction between Ca2+ and the oxygen of the carboxylate group. Our revised parameter set can recover the same three association modes of this aqueous ion pair as the standard parameter set, and yields free energies of binding for the carboxylate–Ca2+ interaction in good agreement with experimental data. The revised parameter set recovers other structural properties of the ion pair in agreement with the standard CHARMM22* parameter set.

No description available.

A great interest has been shown for self-assembled organic nano-structures that can be used in a variety of optoelectronic applications, from element detection to home electronics. It is known from experimental research that sexiphenyl (6P) grown on muscovite mica substrate form uniaxially self-assembled nanofibers which together with sexithiophene (6T) deposited on top gives the possibility to tune their polarized emission. A key to continue develop and explore the full potential of this technique is to understand the mechanisms behind the growth. This thesis investigate the initial growth of 6P and 6T on a 6Pˆ<img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cleft(11%20%5Cbar%7B1%7D%20%5Cright)%20%20%20" /> nanofiber substrate through Molecular Dynamics (MD) simulations. The adsorption of the molecules has been simulated with Simulated Annealing (SA) where 6P align perfectly with the substrate for all coverage while 6T starts to align after a certain amount of coverage. Both molecules show a monotonic increase in the adsorption energy per molecule with an increasing coverage. The surface diffusion of the molecules has been studied and shows a higher movement for both in the direction of the longmolecular axis. / Project P25154-N20 "Hetero-epitaxy of organic-organic nanofibers"

Fur (Ferric Uptake Regulator) est un régulateur transcriptionnel spécifique des bactéries qui intervient dans le contrôle de l'homéostasie du fer, ce qui en fait une cible antibactérienne intéressante. Avant mon arrivée au laboratoire, quatre inhibiteurs interagissant spécifiquement avec Fur avaient été isolés. La partie active de ces inhibiteurs consiste en des peptides de 13 acides aminés. Au cours de cette thèse, j'ai utilisé une double-approche : théorique et expérimentale pour étudier l'interaction de ces peptides avec Fur afin de comprendre le mécanisme d'inhibition. J'ai synthétisé plusieurs séquences peptidiques, montré par des tests biochimiques que certaines inhibaient Fur et déterminé les interactions importantes à l'activité inhibitrice. J'ai obtenu des modèles théoriques des complexes Fur/peptides par amarrage moléculaire, cohérents avec les résultats expérimentaux, qui ont mis en évidence une zone d'inhibition de Fur. Des criblages in silico dans cette zone ont permis de sélectionner de petites molécules, inhibitrices potentielles de Fur et donc intéressantes pour des applications thérapeutiques. / Fur (Ferric Uptake Regulator) is a transcriptional regulator involved in the control of iron homeostasis. Specific to bacteria, Fur is an attractive antibacterial target. Before my arrival in the laboratory, four inhibitors interacting specifically with Fur had been isolated. The active part of these inhibitors consists of peptides of 13 amino acids. In this thesis I have used both theoretical and experimental approaches to study interactions of these peptides with Fur in order to understand the inhibition mechanism. I have synthesized several peptide sequences, shown through biochemical assays that some of them could inhibit Fur and I have identified residues important to the inhibitory activity. I‘ve obtained theoretical models of Fur/peptide complexes consistent with experimental results, which reveal an inhibition pocket in Fur. Small molecules have then been selected though In silico screening of this pocket, that could potentially inhibit Fur, and thus be interesting for therapeutic applications.

The aim of the study was to analyse where different organic molecules situated themselves in relation to the water surface of a water microdroplet and use the resulting data to compare three different forcefields in the simulation package GROMACS. The forcefields used were: General AMBER forcefield (GAFF), Optimized potentials for liquid simulations - all atoms (OPLS-AA), and CHARMM general force field (CGenFF). A library of 146 molecules were simulated using molecular dynamics. Out of the 146 molecules only 65 resulted in useful data for the comparison of the forcefields. The molecules were placed in the centre of a water microdroplet and their movements were simulated for a duration of 1 ns. The trajectories and positions of the molecules were stored and from each simulation a density profile was generated, showing where the molecules situated themselves. The distance from the peak of the density profile to the water surface was calculated and compared between the different forcefields. To analyse the data further some of the molecules were divided into subsets based on their functional groups to see if any trends were visible. Although inconclusive, the data suggested that different forcefields were more or less agreeable depending on the functional group of the molecules, for example OPLS-AA differed from CGenFF and GAFF in the case of alcohols.

Proteins that have been modified by attaching them to a surface or to a polyethylene glycol (PEG) molecule can see many uses in therapeutics and diagnostics -- these unique proteins are called protein devices. Current techniques can perform these functionalizations at a specific residue on the protein, but what remains is identifying what happens to protein structure when mutated, and where to perform the attachment. Both of these issues can be examined using molecular dynamic (MD) simulations. Currently, simulations of the unnatural amino acid (uAA) mutations necessary for protein device functionalization cannot be executed, and full-protein screens of all possible protein device models have never been attempted. Results from this dissertation first employs a new model for simulating PEGylated protein devices building off of previous studies that explore where to attach functional groups. Next, many current assumptions in the community regarding ideal attachment sites are examined. Some of these factors include primary chain location, amino acid type, solvent accessibility, and secondary structure. The focus then turns to novel tertiary structure factors that could influence how well attachment locations affect overall protein device stability. The usefulness of each factor is analyzed to show what factors provided the best predictive power for a site's performance in the screen. A general heuristic is given that could aid in future screens of other protein devices to reduce compute time and quickly identify sites for experimental examination. To explore uAA mutation effects on protein structure, parameters are developed for linear moiety R-groups present in these novel amino acids. The CHARMM and CGenFF force fields currently lack parameters for most linear-angle molecular moieties. This work proposes a method that (1) develops CHARMM parameters for four small molecules that contain terminal azido and alkynyl groups using ffTK, (2) addresses linearity issues, and (3) validates ffTK results via in silico MD simulation. Dihedral analysis examines the linear-angle-containing dihedrals and compares methods for the moiety parameterization. Next, the small molecule parameters are combined with CGenFF to generate parameters for unnatural amino acid MD simulation in a protein. Finally, validation confirms the parameters derived in this work to appropriately simulate unnatural amino acids and small molecules with azido and alkynyl groups.

Proteins that have been modified by attaching them to a surface or to a polyethylene glycol (PEG) molecule can see many uses in therapeutics and diagnostics -- these unique proteins are called protein devices. Current techniques can perform these functionalizations at a specific residue on the protein, but what remains is identifying what happens to protein structure when mutated, and where to perform the attachment. Both of these issues can be examined using molecular dynamic (MD) simulations. Currently, simulations of the unnatural amino acid (uAA) mutations necessary for protein device functionalization cannot be executed, and full-protein screens of all possible protein device models have never been attempted. Results from this dissertation first employs a new model for simulating PEGylated protein devices building off of previous studies that explore where to attach functional groups. Next, many current assumptions in the community regarding ideal attachment sites are examined. Some of these factors include primary chain location, amino acid type, solvent accessibility, and secondary structure. The focus then turns to novel tertiary structure factors that could influence how well attachment locations affect overall protein device stability. The usefulness of each factor is analyzed to show what factors provided the best predictive power for a site's performance in the screen. A general heuristic is given that could aid in future screens of other protein devices to reduce compute time and quickly identify sites for experimental examination. To explore uAA mutation effects on protein structure, parameters are developed for linear moiety R-groups present in these novel amino acids. The CHARMM and CGenFF force fields currently lack parameters for most linear-angle molecular moieties. This work proposes a method that (1) develops CHARMM parameters for four small molecules that contain terminal azido and alkynyl groups using ffTK, (2) addresses linearity issues, and (3) validates ffTK results via in silico MD simulation. Dihedral analysis examines the linear-angle-containing dihedrals and compares methods for the moiety parameterization. Next, the small molecule parameters are combined with CGenFF to generate parameters for unnatural amino acid MD simulation in a protein. Finally, validation confirms the parameters derived in this work to appropriately simulate unnatural amino acids and small molecules with azido and alkynyl groups.

Yes / Materials-binding peptides offer promising routes to the production of tailored Pd nanomaterials in aqueous media, enabling the optimization of catalytic properties. However, the atomic-scale details needed to make these advances are relatively scarce and challenging to obtain. Molecular simulations can provide key insights into the structure of peptides adsorbed at the aqueous Pd interface, provided that the force-field can appropriately capture the relevant bio-interface interactions. Here, we introduce and apply a new polarizable force field, PdP-CHARMM, for the simulation of biomolecule–Pd binding under aqueous conditions. PdP-CHARMM was parametrized with density functional theory (DFT) calculations, using a process compatible with similar polarizable force-fields created for Ag and Au surfaces, ultimately enabling a direct comparison of peptide binding modes across these metal substrates. As part of our process for developing PdP-CHARMM, we provide an extensive study of the performance of ten different dispersion-inclusive DFT functionals in recovering biomolecule–Pd(111) binding. We use the functional with best all-round performance to create PdP-CHARMM.We then employ PdP-CHARMM and metadynamics simulations to estimate the adsorption free energy for a range of amino acids at the aqueous Pd(111) interface. Our findings suggest that only His and Met favor direct contact with the Pd substrate, which we attribute to a remarkably robust interfacial solvation layering. Replica-exchange with solute tempering molecular dynamics simulations of two experimentally-identified Pd-binding peptides also indicate surface contact to be chiefly mediated by His and Met residues at aqueous Pd(111). Adsorption of these two peptides was also predicted for the Au(111) interface, revealing distinct differences in both the solvation structure and modes of peptide adsorption at the Au and Pd interfaces. We propose that this sharp contrast in peptide binding is largely due to the differences in interfacial solvent structuring. / Air Force Office for Scientfi c Research (Grant #FA9550-12-1-0226)

Fur (Ferric Uptake Regulator) est un régulateur transcriptionnel spécifique des bactéries qui intervient dans le contrôle de l'homéostasie du fer, ce qui en fait une cible antibactérienne intéressante. Avant mon arrivée au laboratoire, quatre inhibiteurs interagissant spécifiquement avec Fur avaient été isolés. La partie active de ces inhibiteurs consiste en des peptides de 13 acides aminés. Au cours de cette thèse, j'ai utilisé une double-approche : théorique et expérimentale pour étudier l'interaction de ces peptides avec Fur afin de comprendre le mécanisme d'inhibition. J'ai synthétisé plusieurs séquences peptidiques, montré par des tests biochimiques que certaines inhibaient Fur et déterminé les interactions importantes à l'activité inhibitrice. J'ai obtenu des modèles théoriques des complexes Fur/peptides par amarrage moléculaire, cohérents avec les résultats expérimentaux, qui ont mis en évidence une zone d'inhibition de Fur. Des criblages in silico dans cette zone ont permis de sélectionner de petites molécules, inhibitrices potentielles de Fur et donc intéressantes pour des applications thérapeutiques.

Indiana University-Purdue University Indianapolis (IUPUI) / The multi-protein complex of F0F1 ATP synthase has been of great interest in the fields of microbiology and biochemistry, due to the ubiquitous use of ATP as a biological energy source. Efforts to better understand this complex have been made through structural determination of segments based on NMR and crystallographic data. Some experiments have provided useful data, while others have brought up more questions, especially when structures and functions are compared between bacteria and species with chloroplasts or mitochondria. The epsilon subunit is thought to play a signi cant role in the regulation of ATP synthesis and hydrolysis, yet the exact pathway is unknown due to the experimental difficulty in obtaining data along the transition pathway. Given starting and end point protein crystal structures, the transition pathway of the epsilon subunit was examined through computer simulation.The purpose of this investigation is to determine the likelihood of one such proposed mechanism for the involvement of the epsilon subunit in ATP regulation in bacterial species such as E. coli.