Molecular dynamics simulation of liquid trimethylphosphine

Structural and dynamical properties of liquid trimethylphosphine (TMP), (CH(3))(3)P, as a function of temperature is investigated by molecular dynamics (MD) simulations. The force field used in the MD simulations, which has been proposed from molecular mechanics and quantum chemistry calculations, is able to reproduce the experimental density of liquid TMP at room temperature. Equilibrium structure is investigated by the usual radial distribution function, g(r), and also in the reciprocal space by the static structure factor, S(k). On the basis of center of mass distances, liquid TMP behaves like a simple liquid of almost spherical particles, but orientational correlation due to dipole-dipole interactions is revealed at short-range distances. Single particle and collective dynamics are investigated by several time correlation functions. At high temperatures, diffusion and reorientation occur at the same time range as relaxation of the liquid structure. Decoupling of these dynamic properties starts below ca. 220 K, when rattling dynamics of a given TMP molecules due to the cage effect of neighbouring molecules becomes important. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624408]

SwissSidechain: a molecular and structural database of non-natural sidechains.

Amino acids form the building blocks of all proteins. Naturally occurring amino acids are restricted to a few tens of sidechains, even when considering post-translational modifications and rare amino acids such as selenocysteine and pyrrolysine. However, the potential chemical diversity of amino acid sidechains is nearly infinite. Exploiting this diversity by using non-natural sidechains to expand the building blocks of proteins and peptides has recently found widespread applications in biochemistry, protein engineering and drug design. Despite these applications, there is currently no unified online bioinformatics resource for non-natural sidechains. With the SwissSidechain database (http://www.swisssidechain.ch), we offer a central and curated platform about non-natural sidechains for researchers in biochemistry, medicinal chemistry, protein engineering and molecular modeling. SwissSidechain provides biophysical, structural and molecular data for hundreds of commercially available non-natural amino acid sidechains, both in l- and d-configurations. The database can be easily browsed by sidechain names, families or physico-chemical properties. We also provide plugins to seamlessly insert non-natural sidechains into peptides and proteins using molecular visualization software, as well as topologies and parameters compatible with molecular mechanics software.

Expanding molecular modeling and design tools to non-natural sidechains.

Protein-protein interactions encode the wiring diagram of cellular signaling pathways and their deregulations underlie a variety of diseases, such as cancer. Inhibiting protein-protein interactions with peptide derivatives is a promising way to develop new biological and therapeutic tools. Here, we develop a general framework to computationally handle hundreds of non-natural amino acid sidechains and predict the effect of inserting them into peptides or proteins. We first generate all structural files (pdb and mol2), as well as parameters and topologies for standard molecular mechanics software (CHARMM and Gromacs). Accurate predictions of rotamer probabilities are provided using a novel combined knowledge and physics based strategy. Non-natural sidechains are useful to increase peptide ligand binding affinity. Our results obtained on non-natural mutants of a BCL9 peptide targeting beta-catenin show very good correlation between predicted and experimental binding free-energies, indicating that such predictions can be used to design new inhibitors. Data generated in this work, as well as PyMOL and UCSF Chimera plug-ins for user-friendly visualization of non-natural sidechains, are all available at http://www.swisssidechain.ch. Our results enable researchers to rapidly and efficiently work with hundreds of non-natural sidechains.

Metodologia para obtenção de parâmetros de mecânica molecular aplicados a compostos de coordenação

A methodology is presented to obtain force field parameters to be used in molecular mechanics. The case of Ru(II) is investigated and the parameters obtained, specially its covalent radii, are employed to model Ru(II) coordination compound. The combined use of molecular mechanics with ab initio methods allowed us to predict the metal-ligand stretching force constant for Ru(II) coordination compounds.

Estudo de modelagem molecular de complexos ferriprotoporfirina-IX e quinolinocarbinolaminas antimaláricas: proposta de um farmacóforo

Quinine and quinidine are well-known 4-quinolinecarbinolamines that exhibit antimalarial activity, but, in contrast, their epimers 9-epiquinine and 9-epiquinidine are almost inactive. Literature data are conflicting in describing the 4-quinolinecarbinolamine interaction mode with the molecular target, the ferriprotoporphyrin IX [Fe(III)PPIX]. In the present paper, a pharmacophore is proposed based on the binding of the non-aromatic nitrogen to the iron atom. The 4-quinolinecarbinolamine antimalarials were superimposed on the pharmacophore under consideration and complexes with Fe(III)PPIX were constructed. Conformational analyses of the complexes were performed applying the MM+ molecular mechanics method. The analysis of the complexes showed that the proposed ligand mode is possible although it does not explain the activity differences between epimers. A discussion of the structural aspects is also provided.

Studies in molecular spectroscopy

This thesis is divided into two parts. The first part deals with some studies in molecular mechanics Using spectroscopic data and has four chapters in it. Certain approximation methods for the evaluation of molecular force fields are herein developed The second part, which consists of the last two chaptcrs, deals with infrared spectral studies of ternary liquid systems and a polymer film prepared by glow discharge method.

Metal complexes of a tetraazacyclophane: solution and molecular modelling studies

An alternative synthetic approach to yield the compound 2,3,5,6,8,9,11,14-octahydrobenzo[1][ 1,4,7,10]tetraazacyclotetradecine (bz[14]N-4) is presented. The protonation constants of bz[14]N-4 and the stability constants of its complexes with Ni2+, Cu2+, Zn2+, Cd2+, Pb2+ were determined in H2O at 25degreesC with ionic strength 0.10 mol dm(-3) in KNO3 and they were compared with structurally related macrocycles cyclam (1,4,8,11-tetraazacyclotetradecane) and cyclen (1,4,7,10-tetraazacyclododecane). These studies indicate that only 1 : 1 ( M : L) species are formed in solution, and the ligand exhibits a high affinity for larger ions such as Cd2+ and Pb2+. The X-ray study of [bz[14]N4H3](3+) shows that an inclusion compound with a chloride counter-anion is formed through NH...Cl hydrogen bonds. Spectroscopic data in solution ( electronic and NMR spectra) showed that the macrocycle adopts a planar arrangement upon metal complexation. Molecular mechanics calculations reveal that in spite of the presence of the benzene ring in the macrocyclic framework this ligand can encapsulate metal ions with different stereo-electronic sizes in square planar arrangements. Our results indicate that the presence of the benzene ring in the backbone of the bz[14]N-4 confers a coordination behaviour intermediate between that of cyclam and cyclen.

Molecular dynamics simulations of proteins: can the explicit water model be varied?

In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties.

Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

Use of grape polyphenols against carcinogenesis: Putative molecular mechanisms of action using in vitro and in vivo test systems

Polyphenols are present in foods and beverages and are related to sensorial qualities such as color, bitterness, and astringency, which are relevant in wine, tea, grape juice, and other products. These compounds occur naturally in forms varying from simple phenolic acids to complex polymerized tannins. Thus, it is reasonable to expect that grape-derived products elaborated in the presence of skins and seeds, such as wine and grape juice, are natural sources of flavonoids in the diet. Carcinogenesis is a multistep process that is characterized by genetic, epigenetic, and phenotypic changes. With increasing knowledge of these mechanisms, and the conclusion that most cases of cancer are preventable, efforts have focused on identifying the agents with potential anticancer properties. The use of grape polyphenols against the carcinogenesis process seems to be a suitable alternative for either prevention and/or therapeutic purposes. The aim of this article is to show the molecular data generated from the use of grape polyphenols against carcinogenesis using in vivo and in vitro test systems. © Mary Ann Liebert, Inc. and Korean Society of Food Science and Nutrition.

Spectroscopic properties from hybrid QM, MM molecular dynamics simulation

The central aim of this thesis work is the application and further development of a hybrid quantum mechanical/molecular mechanics (QM/MM) based approach to compute spectroscopic properties of molecules in complex chemical environments from electronic structure theory. In the framework of this thesis, an existing density functional theory implementation of the QM/MM approach is first used to calculate the nuclear magnetic resonance (NMR) solvent shifts of an adenine molecule in aqueous solution. The findings show that the aqueous solvation with its strongly fluctuating hydrogen bond network leads to specific changes in the NMR resonance lines. Besides the absolute values, also the ordering of the NMR lines changes under the influence of the solvating water molecules. Without the QM/MM scheme, a quantum chemical calculation could have led to an incorrect assignment of these lines. The second part of this thesis describes a methodological improvement of the QM/MM method that is designed for cases in which a covalent chemical bond crosses the QM/MM boundary. The development consists in an automatized protocol to optimize a so-called capping potential that saturates the electronic subsystem in the QM region. The optimization scheme is capable of tuning the parameters in such a way that the deviations of the electronic orbitals between the regular and the truncated (and "capped") molecule are minimized. This in turn results in a considerable improvement of the structural and spectroscopic parameters when computed with the new optimized capping potential within the QM/MM technique. This optimization scheme is applied and benchmarked on the example of truncated carbon-carbon bonds in a set of small test molecules. It turns out that the optimized capping potentials yield an excellent agreement of NMR chemical shifts and protonation energies with respect to the corresponding full molecules. These results are very promising, so that the application to larger biological complexes will significantly improve the reliability of the prediction of the related spectroscopic properties.

Molecular Machanics of Water Monomer-Maple Code

This is the Maple code to support the molecular dynamics of a water monomer molecule, allowing investigation of the classical vibrations of this molecule.

Cyclooligomerisation of azido-alkyne-functionalised sugars: synthesis of 1,6-linked cyclic pseudo-galactooligosaccharides and assessment of their sialylation by Trypanosoma cruzi trans-sialidase

Cyclic pseudo-galactooligosaccharides were synthesized by cyclooligomerisation of isomeric azido-alkyne derivatives of beta-D-galactopyranose under Cu(I)-catalysed azide-alkyne 1,3-dipolar cycloaddition reaction conditions. The principal products isolated were cyclic dimers and trimers, with lower amounts of cyclic tetramer and pentamer also evident in some cases. Molecular mechanics calculations suggest very compact but flexible structures for the cyclic trimers, with secondary OH groups exposed outside the macrocycle and available for enzymatic glycosylation. The cyclic dimers and trimers represent a new type of acceptor substrate for Trypanosoma cruzi trans-sialidase, giving rise to doubly and triply sialylated glycomacrocycles, respectively.

Detailed H-1 and C-13 NMR structural assignment and relative stereochemistry determination for three closely related sesquiterpene lactones

A complete analysis of H-1 and C-13 NMR spectra of the trypanocidal sesquiterpene lactone eremantholide C and two of its analogues is described. These structurally similar sesquiterpene lactones were submitted to H-1 NMR, C-13 (H-1) NMR, gCOSY, gHSQC, gHMBC, J-resolved and DPFGSE-NOE NMR techniques. The detailed analysis of those results, correlated to some computational calculations (molecular mechanics), led to the total and unequivocal assignment of all H-1 and C-13 NMR data. The determination of all H-1/H-1 coupling constants and all signal multiplicities, together with the elimination of previous ambiguities were also achieved. Copyright (C) 2008 John Wiley & Sons, Ltd.

Complexes of 1-thia-4,7,10-triazacyclododecane ([12]aneN(3)S) - Synthesis, structure and stability constants

The 12-membered macrocyclic ligand 1-thia-4,7, 10-triazacyclododecane ([12]aneN(3)S) has been synthesised, although upon crystallization from acetonitrile a product in which carbon dioxide had added to one secondary amine in the macrocyclic ring (H[12]aneN(3)SCO(2). H2O) was isolated and subsequently characterised by X-ray crystallography. The protonation constants for [12]aneN(3)S and stability constants with Zn(II), Pb(II), Cd(II) and Cu(II) have been determined either potentiometrically or spectrophotometrically in aqueous solution, and compared with those measured or reported for the ligands 1-oxa-4,7,10-triazacyclododecane ([12]aneN(3)O) and 1,4,7,10-tetraazacyclododecane ([12]aneN(4)). The magnitudes of the stability constants are consistent with trends observed previously for macrocyclic ligands as secondary amine donors are replaced with oxygen and thioether donors although the stability constant for the [Hg([12]aneN(4))](2+) complex has been estimated from an NMR experiment to be at least three orders of magnitude larger than reported previously. Zinc(II), mercury(II), lead(II), copper(II) and nickel(II) complexes of [12]aneN(3)S have been isolated and characterised by X-ray crystallography. In the case of copper(II), two complexes [Cu([12]aneN(3)S)(H2O)](ClO4)(2) and [Cu-2([12]aneN(3)S)(2)(OH)(2)](ClO4)(2) were isolated, depending on the conditions employed. Molecular mechanics calculations have been employed to investigate the relative metal ion size preferences of the [3333], asym-[2424] and sym-[2424] conformation isomers. The calculations predict that the asym-[2424] conformer is most stable for M-N bond lengths in the range 2.00-2.25 Angstrom whilst for the larger metal ions the [3333] conformer is dominant. The disorder seen in the structure of the [Zn([12]aneN(3)S)(NO3)](+) complex is also explained by the calculations. (C) 1999 Elsevier Science Ltd. All rights reserved.