Modelling of the sodium complex of a calixarene tetraester in the 1,3-alternate conformation

This paper describes the use of molecular mechanics to model the geometry of the sodium complex of a calix[4] arene tetraester, in the 1,3-alternate conformation 1. Partial charges were assigned to the calixarene on the basis of semi-empirical (AM1, PM3, MNDO, INDO, CNDO and ZINDO) calculations and the binding of the sodium ion to the calixarene was modelled using molecular mechanics. Agreement between the optimised and X-ray structures of the complex was very good. The effect of placing the cation in different starting positions on the energy-minimised geometry of the complex is described.

Structural and kinetic studies of spin crossover in an Iron(II) complex with a novel tripodal ligand

Configurational and ligand conformational influences on the kinetics of (1)A(1) reversible arrow T-5(2) spin crossover in the Fe(II) complex with the novel tripodal ligand, 1,1,1-tris((N-(2-pyridylmethyl)-N-methylamino)methyl)ethane (tptMetame), have been explored. Despite having six chelate rings and three chiral nitrogen atoms, only one enantiomeric pair of isomers, Delta, SSS, and Lambda, RRR, of the complex ion is observed. The conformation of the three rings forming the upper ''cap'' of the complex structure can be assigned delta or lambda with respect to the 3-fold molecular axis. X-ray data at 300 and 153 K, above and below the critical temperature for the spin transition, show that the conformation of the ligand ''cap'' is the same as the absolute configuration of the complex, with the same Lambda lambda(CAP)(or Delta delta(CAP)) combination prevailing for both the LS ((1)A(1)) and HS (T-5(2)) isomers. Molecular mechanics calculations further show that the ligand energy remains lowest for this Lambda lambda(CAP) (or Delta delta(CAP)) combination at all Fe-N distances over the range spanning the LS and HS isomers. Measurements of the spin crossover relaxation time have been carried out in solution over the temperature range 293-170 K. The observed monophasic relaxation traces are also consistent with the absolute configuration of the complex remaining unaltered during the spin crossover.

The role of the active site residues in human galactokinase: implications for the mechanisms of GHMP kinases.

Galactokinase catalyses the phosphorylation of galactose at the expense of ATP. Like other members of the GHMP family of kinases it is postulated to function through an active site base mechanism in which Asp-186 abstracts a proton from galactose. This asparate residue was altered to alanine and to asparagine by site-directed mutagenesis of the corresponding gene. This resulted in variant enzyme with no detectable galactokinase activity. Alteration of Arg-37, which lies adjacent to Asp-186 and is postulated to assist the catalytic base, to lysine resulted in an active enzyme. However, alteration of this residue to glutamate abolished activity. All the variant enzymes, except the arginine to lysine substitution, were structurally unstable (as judged by native gel electrophoresis in the presence of urea) compared to the wild type. This suggests that the lack of activity results from this structural instability, in addition to any direct effects on the catalytic mechanism. Computational estimations of the pK(a) values of the arginine and aspartate residues, suggest that Arg-37 remains protonated throughout the catalytic cycle whereas Asp-186 has an abnormally high pK(a) value (7.18). Quantum mechanics/molecular mechanics (QM/MM) calculations suggest that Asp-186 moves closer to the galactose molecule during catalysis. The experimental and theoretical studies presented here argue for a mechanism in which the C-1-OH bond in the sugar is weakened by the presence of Asp-186 thus facilitating nucleophilic attack by the oxygen atom on the gamma-phosphorus of ATP.

Modéliser la polarisation électronique par un continuum diélectrique intramoléculaire vers un champ de force polarisable pour la chimie bioorganique

Thèse numérisée par la Division de la gestion de documents et des archives de l'Université de Montréal

Étude des interactions supramoléculaires par modélisation moléculaire

L’avancée des infrastructures informatiques a permis l’émergence de la modélisation moléculaire. À cet effet, une multitude de modèles mathématiques sont aujourd’hui disponibles pour simuler différents systèmes chimiques. À l’aide de la modélisation moléculaire, différents types d’interactions chimiques ont été observés. À partir des systèmes les plus simples permettant l’utilisation de modèles quantiques rigoureux, une série d’approximations a été considérée pour rendre envisageable la simulation de systèmes moléculaires de plus en plus complexes. En premier lieu, la théorie de la fonctionnelle de densité dépendante du temps a été utilisée pour simuler les énergies d’excitation de molécules photoactives. De manière similaire, la DFT indépendante du temps a permis la simulation du pont hydrogène intramoléculaire de structures analogues au 1,3,5-triazapentadiène et la rationalisation de la stabilité des états de transition. Par la suite, la dynamique moléculaire et la mécanique moléculaire ont permis de simuler les interactions d’un trimère d’acide cholique et d’un pyrène dans différents solvants. Cette même méthodologie a été utilisée pour simuler les interactions d’un rotaxane-parapluie à l’interface d’un système biphasique. Finalement, l’arrimage moléculaire et les fonctions de score ont été utilisés pour simuler les interactions intermoléculaires entre une protéine et des milliers de candidats moléculaires. Les résultats ont permis de mettre en place une stratégie de développement d’un nouvel inhibiteur enzymatique.

Conception, synthèse et caractérisation de nouvelles macromolécules branchées biocompatibles pour encapsuler des principes actifs hydrophobes

La vectorisation des médicaments est une approche très prometteuse tant sur le plan médical qu’économique pour la livraison des substances actives ayant une faible biodisponibilité. Dans ce contexte, les polymères en étoile et les dendrimères, macromolécules symétriques et branchées, semblent être les solutions de vectorisation les plus attrayantes. En effet, ces structures peuvent combiner efficacement une stabilité élevée dans les milieux biologiques à une capacité d’encapsulation des principes actifs. Grâce à leur architecture bien définie, ils permettent d’atteindre un très haut niveau de reproductibilité de résultats, tout en évitant le problème de polydispersité. Bien que des nombreuses structures dendritiques aient été proposées ces dernières années, il est cependant à noter que la conception de nouveaux nanovecteurs dendritiques efficaces est toujours d’actualité. Ceci s’explique par des nombreuses raisons telles que celles liées à la biocompatibilité, l’efficacité d’encapsulation des agents thérapeutiques, ainsi que par des raisons économiques. Dans ce projet, de nouvelles macromolécules branchées biocompatibles ont été conçues, synthétisées et évaluées. Pour augmenter leur efficacité en tant qu’agents d’encapsulations des principes actifs hydrophobes, les structures de ces macromolécules incluent un coeur central hydrophobe à base de porphyrine, décanediol ou trioléine modifié et, également, une couche externe hydrophile à base d’acide succinique et de polyéthylène glycol. Le choix des éléments structuraux de futures dendrimères a été basé sur les données de biocompatibilité, les résultats de nos travaux de synthèse préliminaires, ainsi que les résultats de simulation in silico réalisée par une méthode de mécanique moléculaire. Ces travaux ont permis de choisir des composés les plus prometteurs pour former efficacement et d’une manière bien contrôlable des macromolécules polyesters. Ils ont aussi permis d’évaluer au préalable la capacité de futurs dendrimères de capter une molécule médicamenteuse (itraconazole). Durant cette étape, plusieurs nouveaux composés intermédiaires ont été obtenus. L’optimisation des conditions menant à des rendements réactionnels élevés a été réalisée. En se basant sur les travaux préliminaires, l’assemblage de nouveaux dendrimères de première et de deuxième génération a été effectué, en utilisant les approches de synthèse divergente et convergente. La structure de nouveaux composés a été prouvée par les techniques RMN du proton et du carbone 13C, spectroscopie FTIR, UV-Vis, analyse élémentaire, spectrométrie de masse et GPC. La biocompatibilité de produits a été évaluée par les tests de cytotoxicité avec le MTT sur les macrophages murins RAW-262.7. La capacité d’encapsuler les principes actifs hydrophobes a été étudiée par les tests avec l’itraconazole, un antifongique puissant mais peu biodisponible. La taille de nanoparticules formées dans les solutions aqueuses a été mesurée par la technique DLS. Ces mesures ont montré que toutes les structures dendritiques ont tendance à former des micelles, ce qui exclue leurs applications en tant que nanocapsules unimoléculaires. L’activité antifongique des formulations d’itraconazole encapsulé avec les dendrimères a été étudiée sur une espèce d’un champignon pathogène Candida albicans. Ces tests ont permis de conclure que pour assurer l’efficacité du traitement, un meilleur contrôle sur le relargage du principe actif était nécessaire.

Algorithm for generating defective graphene sheets

An algorithm is presented for the generation of molecular models of defective graphene fragments, containing a majority of 6-membered rings with a small number of 5- and 7-membered rings as defects. The structures are generated from an initial random array of points in 2D space, which are then subject to Delaunay triangulation. The dual of the triangulation forms a Voronoi tessellation of polygons with a range of ring sizes. An iterative cycle of refinement, involving deletion and addition of points followed by further triangulation, is performed until the user-defined criteria for the number of defects are met. The array of points and connectivities are then converted to a molecular structure and subject to geometry optimization using a standard molecular modeling package to generate final atomic coordinates. On the basis of molecular mechanics with minimization, this automated method can generate structures, which conform to user-supplied criteria and avoid the potential bias associated with the manual building of structures. One application of the algorithm is the generation of structures for the evaluation of the reactivity of different defect sites. Ab initio electronic structure calculations on a representative structure indicate preferential fluorination close to 5-ring defects.

Dinuclear copper and zinc complexes of a hexaazamacrocycle containing p-xylyl spacers and bridging anions: theoretical and spectroscopic studies

The hexaazamacrocycle 7,22-dimethyl-3,7,11,18,22,26-hexaazatricyclo[26.2.2.2(13,16)] tetratriaconta-1(30), 13,15,28,31,33- hexaene (Me-2[30] pbz(2)N(6)) was synthesized and characterised by single crystal X-ray diffraction. The macrocycle adopts a conformation with the two aromatic rings almost parallel at a distance of ca. 4.24 Angstrom, but displaced relative to each other by ca. 1.51 Angstrom. The protonation constants of this compound and the stability constants of its complexes with Cu2+ and Zn2+, were determined in water - methanol (9 : 1 v/v) at 25 degreesC with ionic strength 0.10 mol dm(-3) in KCl. The potentiometric and spectroscopic studies (NMR of zinc, cadmium and lead complexes, and EPR of the copper complexes) indicate the formation of only dinuclear complexes. The association constants of the dinuclear copper complex with anions ( thiocyanate, terephthalate and glyphosate) and neutral molecules (1,4-benzenedimethanol, p-xylylenediamine and terephthalic acid) were determined at 20 degreesC in methanol. The structural preferences of this ligand and of its dinuclear copper(II) complex with a variety of bridging ligands were evaluated theoretically by molecular mechanics calculations (MM) and molecular dynamics (MD) using quenching techniques.

Metal complexes of a dipyridine octaazamacrocycle: stability constants, structural and modelling studies

Two 28-membered octaazamacrocycles, [28]py(2)N(6) and Me-2[28]py(2)N(6), have been synthesized. The protonation constants of the N-methyl. derivative and the stability constants of its complexes with Ni2+, Cu2+, Zn2+, Cd2+, and Pb2+ were determined at 25degreesC in 0.10 mol dm(-3) KNO3. The high overall basicity of Me-2[28]py(2)N(6) is ascribed to the weaker repulsion between protonated contiguous charged ammonium sites separated by propyl chains. These studies together with NMR, UV-vis and EPR spectroscopies indicated the presence of mono- and di-nuclear species, The single crystal structure of the complex [Ni-2([28]py(2)N(6))(H2O)(4)]Cl-4.3H(2)O was determined, and showed each nickel centre in a distorted octahedral co-ordination environment. The nickel centres are held within the macrocycle at a large distance of 6.991(g) Angstrom from each other. The formation of mononuclear complexes was evaluated theoretically via molecular mechanics (MM) and molecular dynamics (MD) calculations and showed that these large macrocycles have sufficient flexibility to encapsulate metal ions with different stereo-electronic sizes. Structures for small and large metal ions are proposed.

Mapping antimalarial pharmacophores as a useful tool for the rapid discovery of drugs effective in vivo: design, construction, characterization, and pharmacology of metaquine

Resistant strains of Plasmodium falciparum and the unavailability of useful antimalarial vaccines reinforce the need to develop new efficacious antimalarials. This study details a pharmacophore model that has been used to identify a potent, soluble, orally bioavailable antimalarial bisquinoline, metaquine (N,N'-bis(7-chloroquinolin-4-yl)benzene-1,3-diamine) (dihydrochloride), which is active against Plasmodium berghei in vivo (oral ID50 of 25 mu mol/kg) and multidrug-resistant Plasmodium falciparum K1 in vitro (0.17 mu M). Metaquine shows strong affinity for the putative antimalarial receptor, heme at pH 7.4 in aqueous DMSO. Both crystallographic analyses and quantum mechanical calculations (HF/6-31+G*) reveal important regions of protonation and bonding thought to persist at parasitic vacuolar pH concordant with our receptor model. Formation of drug-heme adduct in solution was confirmed using high-resolution positive ion electrospray mass spectrometry. Metaquine showed strong binding with the receptor in a 1: 1 ratio (log K = 5.7 +/- 0.1) that was predicted by molecular mechanics calculations. This study illustrates a rational multidisciplinary approach for the development of new 4-aminoquinoline antimalarials, with efficacy superior to chloroquine, based on the use of a pharmacophore model.

Development of novel brush-type chiral stationary phases based on terpenoid selectors: HPLC evaluation and theoretical investigation of enantioselective binding interactions

The terpenoid chiral selectors dehydroabietic acid, 12,14-dinitrodehydroabietic acid and friedelin have been covalently linked to silica gel yielding three chiral stationary phases CSP 1, CSP 2 and CSP 3, respectively. The enantiodiscriminating capability of each one of these phases was evaluated by HPLC with four families of chiral aromatic compounds composed of alcohols, amines, phenylalanine and tryptophan amino acid derivatives and beta-lactams. The CSP 3 phase, containing a selector with a large friedelane backbone is particularly suitable for resolving free alcohols and their derivatives bearing fluorine substituents, while CSP 2 with a dehydroabietic architecture is the only phase that efficiently discriminates 1, 1'-binaphthol atropisomers. CSP 3 also gives efficient resolution of the free amines. All three phases resolve well the racemates of N-trifluoracetyl and N-3,5-dinitrobenzoyl phenylalanine amino acid ester derivatives. Good enantioseparation of beta-lactams and N-benzoyl tryptophan amino acid derivatives was achieved on CSP 1. In order to understand the structural factors that govern the chiral molecular recognition ability of these phases, molecular dynamics simulations were carried out in the gas phase with binary diastereomeric complexes formed by the selectors of CSP 1 and CSP 2 and several amino acid derivatives. Decomposition of molecular mechanics energies shows that van der Waals interactions dominate the formation of the diastereomeric transient complexes while the electrostatic binding interactions are primarily responsible for the enantioselective binding of the (R)- and (S)-analytes. Analysis of the hydrogen bonds shows that electrostatic interactions are mainly associated with the formation of N-(HO)-O-...=C enantio selective hydrogen bonds between the amide binding sites from the selectors and the carbonyl groups of the analytes. The role of mobile phase polarity, a mixture of n-hexane and propan-2-ol in different ratios, was also evaluated through molecular dynamics simulations in explicit solvent. (c) 2006 Elsevier Ltd. All rights reserved.

Synthesis and properties of polyaromatic dendrimers possessing a repetitive amide-ester coupling sequence

A series of novel polyaromatic dendrimers that feature tris-(2-ethylamino)amine as the central core unit has been synthesized up to the third generation by employing a convergent growth strategy. The building blocks 1,3-diamino-2-hydroxypropane and 4-carboxybenzaldehyde were used for dendron construction, a process that involved the cyclic repetition of esterification, oxidation and selective amidation steps. Molecular modelling of this class of dendrimers has been used to predict potential solution state conformations employing molecular mechanics and molecular dynamic simulations. In addition, the results of preliminary metal binding studies using the first generation dendritic system are also outlined. (C) 2003 Elsevier Science Ltd. All rights reserved.

Computer simulations of structures, energetics and dynamics of myoglobin center dot center dot center dot ligand complexes

Myoglobin has been studied in considerable detail using different experimental and computational techniques over the past decades. Recent developments in time-resolved spectroscopy have provided experimental data amenable to detailed atomistic simulations. The main theme of the present review are results on the structures, energetics and dynamics of ligands ( CO, NO) interacting with myoglobin from computer simulations. Modern computational methods including free energy simulations, mixed quantum mechanics/molecular mechanics simulations, and reactive molecular dynamics simulations provide insight into the dynamics of ligand dynamics in confined spaces complementary to experiment. Application of these methods to calculate and understand experimental observations for myoglobin interacting with CO and NO are presented and discussed.

A complete atomistic model of molten polyethylene from neutron scattering data: a new methodology for polymer structure

A new approach to the study of the local organization in amorphous polymer materials is presented. The method couples neutron diffraction experiments that explore the structure on the spatial scale 1–20 Å with the reverse Monte Carlo fitting procedure to predict structures that accurately represent the experimental scattering results over the whole momentum transfer range explored. Molecular mechanics and molecular dynamics techniques are also used to produce atomistic models independently from any experimental input, thereby providing a test of the viability of the reverse Monte Carlo method in generating realistic models for amorphous polymeric systems. An analysis of the obtained models in terms of single chain properties and of orientational correlations between chain segments is presented. We show the viability of the method with data from molten polyethylene. The analysis derives a model with average C-C and C-H bond lengths of 1.55 Å and 1.1 Å respectively, average backbone valence angle of 112, a torsional angle distribution characterized by a fraction of trans conformers of 0.67 and, finally, a weak interchain orientational correlation at around 4 Å.

Theoretical Analysis of Aggregation in Block-Copolymer Films: The Optical Signature

We focus this work on the theoretical investigation of the block-copolymer poly [oxyoctyleneoxy-(2,6-dimethoxy-1,4phenylene-1,2-ethinylene-phenanthrene-2,4diyl) named as LaPPS19, recently proposed for optoelectronic applications. We used for that a variety of methods, from molecular mechanics to quantum semiempirical techniques (AMI, ZINDO/S-CIS). Our results show that as expected isolated LaPPS19 chains present relevant electron localization over the phenanthrene group. We found, however, that LaPPS19 could assemble in a pi-stacked form, leading to impressive interchain interaction; the stacking induces electronic delocalization between neighbor chains and introduces new states below the phenanthrene-related absorption; these results allowed us to associate the red-shift of the absorption edge, seen in the experimental results, to spontaneous pi-stack aggregation of the chains. (C) 2009 Wiley Periodicals, Inc. Int J Quantum Chem 110: 885-892, 2010