Cation complexation by chemically modified calixarenes. 11. Complexation and extraction of alkali cations by calix[5]- and -[6]arene ketones. Crystal and molecular structures of calix[5]arene ketones and Na+ and Rb+ complexes

A series of four calix[5]arenes and three calix[6]arenes (R-calixarene-OCH2COR1) (R = H or Bu-t) with alkyl ketone residues (R-1 = Me or Bu-t) on the lower rim have been synthesized, and their affinity for complexation of alkali cations has been assessed through phase-transfer experiments and stability constant measurements. The conformations of these ketones have been probed by H-1 NMR and X-ray diffraction analysis, and by molecular mechanics calculations. Pentamer 3 (R R-1 = Bu-t) possesses a symmetrical cone conformation in solution and a very distorted cone conformation in the solid state. Pentamer 5 (R = H, R-1 = Bu-t) exists in a distorted 1,2-alternate conformation in the solid state, but in solution two slowly interconverting conformations, one a cone and the other presumed to be 1,2-alternate, can be detected. X-ray structure analysis of the sodium and rubidium perchlorate complexes of 3 reveal the cations deeply encapsulated by the ethereal and carbonyl oxygen atoms in distorted cone conformations which can be accurately reproduced by molecular mechanics calculations. The phase-transfer and stability constant data reveal that the extent of complexation depends on calixarene size and the nature of the alkyl residues adjacent to the ketonic carbonyls with tert-butyl much more efficacious than methyl.

Molecular modeling of calixarenes with group I metal ions

Molecular mechanics calculations have been used to model the geometries of the complexes of Group I metal ions with calix[n]arenes (n = 4,5). A simple procedure in which the calixarene atoms are assigned partial charges on the basis of AM1 calculations and the metal ions are allowed to bind electrostatically to the calixarenes produces surprising good results when the resulting structures are compared to known crystallographic data on the complexes. Encapsulated solvent molecules and/or counterions can be included in the calculations and, indeed, are necessary to reproduce the X-ray data.

Binding modes of diketo-acid inhibitors of HIV-1 integrase: A comparative molecular dynamics simulation study

HIV-1 integrase (IN) has become an attractive target since drug resistance against HIV-1 reverse transcriptase (RT) and protease (PR) has appeared. Diketo acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN: however the action mechanism is not well understood. Here, to study the inhibition mechanism of DKAs we performed 10 ns comparative molecular dynamics simulations on HIV-1 IN bound with three most representative DMA inhibitors: Shionogi inhibitor, S-1360 and two Merck inhibitors L-731,988 and L-708,906. Our simulations show that the acidic part of S-1360 formed salt bridge and cation-pi interactions with Lys159. In addition, the catalytic Glu152 in S-1360 was pushed away from the active site to form an ion-pair interaction with Arg199. The Merck inhibitors can maintain either one or both of these ion-pair interaction features. The difference in potencies of the DMA inhibitors is thus attributed to the different binding modes at the catalytic site. Such structural information at atomic level, not only demonstrates the action modes of DMA inhibitors but also provides a novel starting point for structural-based design of HIV-1 IN inhibitors.

The Plasmodium falciparum Malaria M1 Alanyl Aminopeptidase (PfA-M1): Insights of Catalytic Mechanism and Function from MD Simulations

Malaria caused by several species of Plasmodium is major parasitic disease of humans, causing 1-3 million deaths worldwide annually. The widespread resistance of the human parasite to current drug therapies is of major concern making the identification of new drug targets urgent. While the parasite grows and multiplies inside the host erythrocyte it degrades the host cell hemoglobin and utilizes the released amino acids to synthesize its own proteins. The P. falciparum malarial M1 alanyl-aminopeptidase (PfA-M1) is an enzyme involved in the terminal stages of hemoglobin digestion and the generation of an amino acid pool within the parasite. The enzyme has been validated as a potential drug target since inhibitors of the enzyme block parasite growth in vitro and in vivo. In order to gain further understanding of this enzyme, molecular dynamics simulations using data from a recent crystal structure of PfA-M1 were performed. The results elucidate the pentahedral coordination of the catalytic Zn in these metallo-proteases and provide new insights into the roles of this cation and important active site residues in ligand binding and in the hydrolysis of the peptide bond. Based on the data, we propose a two-step catalytic mechanism, in which the conformation of the active site is altered between the Michaelis complex and the transition state. In addition, the simulations identify global changes in the protein in which conformational transitions in the catalytic domain are transmitted at the opening of the N-terminal 8 angstrom-long channel and at the opening of the 30 angstrom-long C-terminal internal chamber that facilitates entry of peptides to the active site and exit of released amino acids. The possible implications of these global changes with regard to enzyme function are discussed.

The molecular dynamics of Trypanosoma brucei UDP-galactose 4'-epimerase:a drug target for African sleeping sickness

During the past century, several epidemics of human African trypanosomiasis, a deadly disease caused by the protist Trypanosoma brucei, have afflicted sub-Saharan Africa. Over 10 000 new victims are reported each year, with hundreds of thousands more at risk. As current drug treatments are either highly toxic or ineffective, novel trypanocides are urgently needed. The T. brucei galactose synthesis pathway is one potential therapeutic target. Although galactose is essential for T. brucei survival, the parasite lacks the transporters required to intake galactose from the environment. UDP-galactose 4'-epimerase (TbGalE) is responsible for the epimerization of UDP-glucose to UDP-galactose and is therefore of great interest to medicinal chemists. Using molecular dynamics simulations, we investigate the atomistic motions of TbGalE in both the apo and holo states. The sampled conformations and protein dynamics depend not only on the presence of a UDP-sugar ligand, but also on the chirality of the UDP-sugar C4 atom. This dependence provides important insights into TbGalE function and may help guide future computer-aided drug discovery efforts targeting this protein.

A Neutron Diffraction and Molecular Dynamics Investigation of Acetate-based Ionic Liquids as Solvents for Glucose

The liquid state structure of the ionic liquid, 1-ethyl-3-methylimidazolium acetate, and the solute/solvent structure of glucose dissolved in the ionic liquid at a 1: 6 molar ratio have been investigated at 323 K by molecular dynamics simulations and neutron diffraction experiments using H/D isotopically substituted materials. Interactions between hydrogen-bond donating cation sites and polar, directional hydrogen-bond accepting acetate anions are examined. Ion-ion radial distribution functions for the neat ionic liquid, calculated from both MD and derived from the empirical potential structure refinement model to the experimental data, show the alternating shell-structure of anions around the cation, as anticipated. Spatial probability distributions reveal the main anion-to-cation features as in-plane interactions of anions with imidazolium ring hydrogens and cation-cation planar stacking. Interestingly, the presence of the polarised hydrogen-bond acceptor anion leads to increased anion-anion tail-tail structuring within each anion shell, indicating the onset of hydrophobic regions within the anion regions of the liquid.

PLUMED 2: New feathers for an old bird

Enhancing sampling and analyzing simulations are central issues in molecular simulation. Recently, we introduced PLUMED, an open-source plug-in that provides some of the most popular molecular dynamics (MD) codes with implementations of a variety of different enhanced sampling algorithms and collective variables (CVs). The rapid changes in this field, in particular new directions in enhanced sampling and dimensionality reduction together with new hardware, require a code that is more flexible and more efficient. We therefore present PLUMED 2 here a,complete rewrite of the code in an object-oriented programming language (C++). This new version introduces greater flexibility and greater modularity, which both extends its core capabilities and makes it far easier to add new methods and CVs. It also has a simpler interface with the MD engines and provides a single software library containing both tools and core facilities. Ultimately, the new code better serves the ever-growing community of users and contributors in coping with the new challenges arising in the field.

Program summary

Program title: PLUMED 2

Catalogue identifier: AEEE_v2_0

Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEEE_v2_0.html

Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland

Licensing provisions: Yes

No. of lines in distributed program, including test data, etc.: 700646

No. of bytes in distributed program, including test data, etc.: 6618136

Distribution format: tar.gz

Programming language: ANSI-C++.

Computer: Any computer capable of running an executable produced by a C++ compiler.

Operating system: Linux operating system, Unix OSs.

Has the code been vectorized or parallelized?: Yes, parallelized using MPI.

RAM: Depends on the number of atoms, the method chosen and the collective variables used.

Classification: 3, 7.7, 23. Catalogue identifier of previous version: AEEE_v1_0.

Journal reference of previous version: Comput. Phys. Comm. 180 (2009) 1961.

External routines: GNU libmatheval, Lapack, Bias, MPI. (C) 2013 Elsevier B.V. All rights reserved.

Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems

The generalized Langevin equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general nonequilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e. g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly nontrivial, especially if the memory kernels decay nonexponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this paper, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.

Protection of DNA against low-energy electrons by amino acids: A first-principles molecular dynamics study

Using first-principles molecular dynamics simulations, we have investigated the notion that amino acids can play a protective role when DNA is exposed to excess electrons produced by ionizing radiation. In this study we focus on the interaction of glycine with the DNA nucleobase thymine. We studied thymine-glycine dimers and a condensed phase model consisting of one thymine molecule solvated in amorphous glycine. Our results show that the amino acid acts as a protective agent for the nucleobase in two ways. If the excess electron is initially captured by the thymine, then a proton is transferred in a barrier-less way from a neighboring hydrogen-bonded glycine. This stabilizes the excess electron by reducing the net partial charge on the thymine. In the second mechanism the excess electron is captured by a glycine, which acts as a electron scavenger that prevents electron localization in DNA. Both these mechanisms introduce obstacles to further reactions of the excess electron within a DNA strand, e.g. by raising the free energy barrier associated with strand breaks.

The diamond pyramid structure in electroless copper deposit, its atomic model and molecular dynamics simulation

In this seminar, I will talk about the discovery of the diamond pyramid structures in the electroless copper deposits on both epoxy and stainless steel substrates. The surface morphology of the structure was characterized with scanning electron microscopy (SEM). According to the morphological feature of the structure, an atom model was brought forward in order to describe the possible mechanism of forming such structure. Molecular dynamics simulations were then carried out to investigate the growing process of the diamond pyramid structure. The final structures of the simulation were compared with the SEM images and the atomic model. The radial distribution function of the final structures of the simulation was compared with that calculated from the X-ray diffraction pattern of the electroless copper deposit sample.

Molecular dynamics simulation (MDS) to study nanoscale machining processes

Molecular Dynamics Simulations (MDS) are constantly being used to make important contributions to our fundamental understanding of material behaviour, at the atomic scale, for a variety of thermodynamic processes. This chapter shows that molecular dynamics simulation is a robust numerical analysis tool in addressing a range of complex nanofinishing (machining) problems that are otherwise difficult or impossible to understand using other methods. For example the mechanism of nanometric cutting of silicon carbide is influenced by a number of variables such as machine tool performance, machining conditions, material properties, and cutting tool performance (material microstructure and physical geometry of the contact) and all these variables cannot be monitored online through experimental examination. However, these could suitably be studied using an advanced simulation based approach such as MDS. This chapter details how MD simulation can be used as a research and commercial tool to understand key issues of ultra precision manufacturing research problems and a specific case was addressed by studying diamond machining of silicon carbide. While this is appreciable, there are a lot of challenges and opportunities in this fertile area. For example, the world of MD simulations is dependent on present day computers and the accuracy and reliability of potential energy functions [109]. This presents a limitation: Real-world scale simulation models are yet to be developed. The simulated length and timescales are far shorter than the experimental ones which couples further with the fact that contact loading simulations are typically done in the speed range of a few hundreds of m/sec against the experimental speed of typically about 1 m/sec [17]. Consequently, MD simulations suffer from the spurious effects of high cutting speeds and the accuracy of the simulation results has yet to be fully explored. The development of user-friendly software could help facilitate molecular dynamics as an integral part of computer-aided design and manufacturing to tackle a range of machining problems from all perspectives, including materials science (phase of the material formed due to the sub-surface deformation layer), electronics and optics (properties of the finished machined surface due to the metallurgical transformation in comparison to the bulk material), and mechanical engineering (extent of residual stresses in the machined component) [110]. Overall, this chapter provided key information concerning diamond machining of SiC which is classed as hard, brittle material. From the analysis presented in the earlier sections, MD simulation has helped in understanding the effects of crystal anisotropy in nanometric cutting of 3C-SiC by revealing the atomic-level deformation mechanisms for different crystal orientations and cutting directions. In addition to this, the MD simulation revealed that the material removal mechanism on the (111) surface of 3C-SiC (akin to diamond) is dominated by cleavage. These understandings led to the development of a new approach named the “surface defect machining” method which has the potential to be more effective to implement than ductile mode micro laser assisted machining or conventional nanometric cutting.

Estudo de receptores sintéticos no reconhecimento molecular de substratos aniónicos e iões metálicos

O trabalho descrito insere-se no âmbito da Química Supramolecular e consistiu no desenvolvimento de receptores artificiais, na forma protonada ou complexada, para o reconhecimento molecular de substratos aniónicos derivados de ácidos carboxílicos incluindo os herbicidas PMG2-, ATCP- e 2,4-D-. Foram investigadas duas séries de aniões, uma alifática (ox2-, mal2-, suc2-, glu2-, adip2-, cit3- e cta3-) e outra aromática (bzc-, naphc-, anthc-, pyrc-, ph2-, iph2-, tph2-, btc3-, dihyac2-, 4,4-dibzc2-, 3-nitrobzc- e 4-nitrobzc-). Foram sintetizados sete novos ligandos macrocíclicos simétricos constituídos por anéis aromáticos piridina ou fenantrolina ligados por cadeias de poliaminas saturadas. O comportamento ácido-base destes macrociclos foi investigado em solução aquosa e as constantes de protonação correspondentes determinadas por métodos potenciométricos e de RMN de 1H. As propriedades de complexação destes ligandos com os iões metálicos Ni2+, Cu2+, Zn2+, Cd2+ e Pb2+ foram também estudadas por métodos potenciométricos nas mesmas condições experimentais, tendo revelado que os macrociclos de dimensão média são capazes de acomodar um ou dois iões metálicos. O complexo dinuclear de Cu(II) derivado do macrociclo com dois grupos piridina foi utilizado como receptor de aniões carboxilato originando complexos ternários. Todos os complexos foram caracterizados em solução por espectroscopias de UV/vis/IVpróx e de RMN. As espécies paramagnéticas foram também caracterizadas por espectroscopia de RPE. A formação de espécies ternárias foi ainda investigada por espectrometria de massa ESI-MS e ESI-MS/MS. As estruturas cristalinas de alguns dos complexos foram determinadas por difracção de raios-X. As formas protonadas dos macrociclos foram utilizadas como receptores de uma grande variedade de aniões carboxilato. O reconhecimento molecular entre os receptores e os substratos aniónicos foi investigado em solução por métodos potenciométricos e de espectroscopia de RMN com determinação das constantes de associação. Os agregados supramoleculares foram caracterizados no estado sólido por difracção de raios-X. Finalmente as associações supramoleculares foram estudadas em solução por métodos de dinâmica molecular com determinação dos termos entrópicos e entálpicos das energias livres de ligação. Em resumo, nesta tese apresentam-se os resultados de estudos realizados com duas famílias de macrociclos: desde a síntese dos compostos, passando por estudos em solução e finalizando com simulação molecular. Este estudo sistemático através da conjugação de metodologias complementares permitiu caracterizar ao nível macroscópico e microscópico as associações moleculares.

Optimization of molecular dynamics simulation code and applications to biomolecular systems

Tese de doutoramento, Bioquimica, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2015

Assessing the Reproducibility of Clustering of Molecular Dynamics Conformations on Self-Organizing Maps

The goal of most clustering algorithms is to find the optimal number of clusters (i.e. fewest number of clusters). However, analysis of molecular conformations of biological macromolecules obtained from computer simulations may benefit from a larger array of clusters. The Self-Organizing Map (SOM) clustering method has the advantage of generating large numbers of clusters, but often gives ambiguous results. In this work, SOMs have been shown to be reproducible when the same conformational dataset is independently clustered multiple times (~100), with the help of the Cramérs V-index (C_v). The ability of C_v to determine which SOMs are reproduced is generalizable across different SOM source codes. The conformational ensembles produced from MD (molecular dynamics) and REMD (replica exchange molecular dynamics) simulations of the penta peptide Met-enkephalin (MET) and the 34 amino acid protein human Parathyroid Hormone (hPTH) were used to evaluate SOM reproducibility. The training length for the SOM has a huge impact on the reproducibility. Analysis of MET conformational data definitively determined that toroidal SOMs cluster data better than bordered maps due to the fact that toroidal maps do not have an edge effect. For the source code from MATLAB, it was determined that the learning rate function should be LINEAR with an initial learning rate factor of 0.05 and the SOM should be trained by a sequential algorithm. The trained SOMs can be used as a supervised classification for another dataset. The toroidal 10×10 hexagonal SOMs produced from the MATLAB program for hPTH conformational data produced three sets of reproducible clusters (27%, 15%, and 13% of 100 independent runs) which find similar partitionings to those of smaller 6×6 SOMs. The χ^2 values produced as part of the C_v calculation were used to locate clusters with identical conformational memberships on independently trained SOMs, even those with different dimensions. The χ^2 values could relate the different SOM partitionings to each other.

Elementary steps in aqueous proton transfer reactions : a first principles molecular dynamics study

La nature des acides dans un environnement aqueux est primordiale dans de nombreux aspects de la chimie et de la biologie. La caractéristique principale d'un acide est sa capacité à transférer un proton vers une molécule d'eau ou vers n'importe quelle base, mais ce procédé n'est pas aussi simple qu'il y paraît. Il peut au contraire être extrêmement complexe et dépendre de manière cruciale de la solvatation des différents intermédiaires de réaction impliqués. Cette thèse décrit les études computationnelles basées sur des simulations de dynamique moléculaire ab initio qui ont pour but d'obtenir une description à l'échelle moléculaire des divers procédés de transferts de proton entre acide et bases dans un milieu aqueux. Pour cela, nous avons étudié une serie de système, dont l'acide hydrofluorique aqueux, l'acide trifluoroacétique aqueux, et un système modèle constitué d'un phénol et d'une entité carboxylate reliés entre eux par une molécule d'eau en solution aqueuse. Deux états intermédiaires ont été identifiés pour le transfert d'un proton depuis un acide. Ces intermédiaires apparaissent stabilisés par un motif local de solvatation via des ponts H. Leurs signatures spectroscopiques ont été caractérisées au moyen de la spectroscopie infrarouge, en utilisant le formalisme de la dynamique moléculaire ab initio, qui inclut l'effet quantique nucléaire de manière explicite. Cette étude a aussi identifié trois chemins de réaction élémentaire, qui sont responsable pour le transfert d'un proton d'un acide à une base, ainsi que leurs échelles de temps caractéristiques. Les conclusions tirées de ces études sont discutées dans les détails, au niveau moléculaire, avec une emphase sur les comparaisons entre les résultats théoriques et les mesures expérimentales obtenues dans a littérature ou via des collaborateurs.