400 likes | 521 Views
Theoretical Inorganic Chemistry Group. Structure and dynamics of b -cyclodextrin and glycine at quantum mechanical level. Hélio A. Duarte , Hélio F. Dos Santos, Thomas Heine, Serguei Patchkovskii duarteh@ufmg.br Department of Chemistry - ICEx,
E N D
Theoretical Inorganic Chemistry Group Structure and dynamics of b-cyclodextrin and glycine at quantum mechanical level Hélio A. Duarte, Hélio F. Dos Santos, Thomas Heine, Serguei Patchkovskii duarteh@ufmg.br Department of Chemistry - ICEx, Federal University of Minas Gerais - UFMG ACS 232nd National Meeting San Francisco, CA - USA
Outline • Motivation • Spironolactone and its Complexes with b-cyclodextrin • b-cyclodextrine in aqueous solution – molecular dynamics using DFTB/MM approach. • Glycine in aqueous solution – molecular dynamics using full DFTB.
b-cyclodextrine Consists of 7 D-glucose linked by a (1-4) interglucose bonds.
b-Cyclodextrine • Inclusion compounds • Drug Delivery Systems • Improved molecular switches • Artificial enzymes • Rotaxamers • Nanoreactors • Self-assembling systems
Spironolactone and its Complexes with b-cyclodextrin Lula, Gomes, Piló-Veloso, De Noronha, Duarte, Santos, Sinisterra, J. Inclusion Phenon. Macroc. Chem., (2006).
Spironolactone : b-cyclodextrin • Complexes 1:1 and 1:2 are formed and well characterized by ROESY-NMR. • The rings A and DE are involved in the inclusion process. Simulation at gas phase: DC-SCC-DFTB *Zhechkov, L.; Heine, T.; Patchkovskii, S.; Seifert, G.; Duarte, H. A. JCTC2005, 1, 841. *Elstner, et al., Phys. Rev. B, 1998, 58, 7260. *Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, R. Physical Review B1995, 51, 12947
DC-SCC-DFTB calculations of 1:1 complexes at gas phase A-Head
DC-SCC-DFTB calculations of 1:2 complexes at gas phase Head-Head arrangement.
Inclusion process: guest:host according to Rekharsky and Inoue, Chem. Rev., 98, 1875. • penetration of the hydrophobic part of the guest molecule into the cylodextrin cavity • dehydration of the organic guest. • hydrogen bonding interactions • release of the water molecules to bulk water • conformational changes or strain release of the CyD upon complexation • how many water molecules are inside of the cavity before and after complexation.
Methodology • Born-Oppenheimer Molecular Dynamics • QM/MM calculations • QM : DC-SCC-DFTB* method • MM: employs Rappé’s universal force field (UFF). • Cubic box with a lattice vector length of 34.92 Å. • 1385 water molecules and b-CyD. • Microcanonical NVE ensemble. • MD run: 160 ps with a time step of 0.5 fs. • Program: deMon program (NRC-2004, Canada) *Zhechkov, L.; et al. JCTC2005, 1, 841. *Elstner, et al., Phys. Rev. B, 1998, 58, 7260. *Porezag, D. et al. Physical Review B1995, 51, 12947
Setup of the simulation: The periodic simulation box is given. b-CyD, given in bold, is treated quantum mechanically. The surrounding waters (wireframe model) and all solute-solvent interactions are approximated with the universal force field (UFF) employing TIP3P partial charges on water.
Dihedral angles: C2C3C4C5O4O4’O4’’O4’’’ Angles: C1O4’C4’O4O4’O4’’
The root mean square deviations (RMSD) of the coordinates between two snapshots of a MD trajectory provides information about the flexibility of the b-CyD. The water surrounding the b-CyD acts as a cushion, decreasing its free motion. Figure 2. RMSQ for b-CyD in gas phase (dashed) and in solution (full).
Figure 3. Configuration space taken by b-CyD in aqueous solution.
In the radial distribution function (RDF), the range of r below 4.2 A corresponds to the encapsulated water molecules and integrates to 7.9. This is in agreement with X-ray and neutron diffraction studies, which arrived at 7 water molecules. Figure 4. RDF with respect to the distance between the centres of mass of b-CyD and water molecules.
Motion of the water molecules in the cavity of b-cyclodextrine. Solvent water molecules were removed for better viewing.
Figure 5. Configurational space taken by the water molecules encapsulated in b-CyD. For sake of clarity, only the initial structure of b-CyD is shown.
Table I. Average number and oxygen-oxygen distance of hydrogen bonds between water and -CyD. 91% of the HBs formed with the glycosidic (O4) and 64.8 % of the pyranoid (O5) oxygens are due to the encapsulated water molecules. For the primary (O6) and secondary (O2,O3) hydroxyls, 96% of HBs are due to outer solvent.
65% 90%
Dwell time of water molecules in the cavity A=33.2A2 • No preferential side for the water molecules to enter the cavity. • Roughly 50% of the water molecules come inside and get out through the top side. A=28.3A2
Fig. 6. Dwell time distribution of the water molecules. There is strong peak at 70 fs dwell time of the encapsulated water molecules. Much longer dwell times are possible, up to several ps.
Angiotensine(1-7):Cyclodextrine Preliminary Results NOESY-NMR TYR (H3/H5 and H2/H6) and b-CyD (H3 and H5) The chemical structure of angiotensin (1-7), [AspArgValTyrIleHisPro]
Methodology • Born-Oppenheimer Molecular Dynamics • QM/MM calculations • QM : DC-DFTB* method • MM: employs Rappé’s universal force field (UFF). • Cubic box with a lattice vector length of 61.0 Å. • 7381 water molecules and Ang(1-7):b-CyD. • Microcanonical NVE ensemble. • MD run: with a time step of 0.5 fs. • Program: deMon program (NRC-2004, Canada) *Zhechkov, L.; et al. JCTC2005, 1, 841. * Elstner, et al., Phys. Rev. B, 1998, 58, 7260. *Porezag, D. et al. Physical Review B1995, 51, 12947
ANG:CYD → Preliminary Results Structural parameters of the ang(1-7):b-CyD
Glycine in Aqueous Solution Progress report Neutral form Zwitterionic form
Methodology • Born-Oppenheimer Molecular Dynamics • QM : DC-DFTB* method • Cubic box with a lattice vector length of 16.0 Å. • 129 water molecules and glycine. • Microcanonical NVE ensemble. • MD run: 100 ps with a time step of 0.5 fs. • Program: deMon program (NRC-2004, Canada) *Zhechkov, L.; et al. JCTC2005, 1, 841. * Elstner, et al., Phys. Rev. B, 1998, 58, 7260. *Porezag, D. et al. Physical Review B1995, 51, 12947
RDF with respect to the distance between the centres of mass of glycine and water. 22 water molecules in the first solvation shell.
Geometrical Properties of glycine 1 5 2 4 3
Thermodynamical properties Neutral Zwitterion DFTB-MD: DENVE = -25.5 kcal/mol PBE/TZVP/UAHF-PCM: DG = -23.4 kcal/mol Exp*. : DH = -10.3 kcal/mol DG = -7.2 kcal/mol * Quoted from Wada et al., Bull. Chem.Soc. Jpn, 55, 3064 (1992).
Grupo de Pesquisa em Química Inorgânica Teórica - GPQIT • Collaborators: • Prof. Ruben Sinisterra (DQ-UFMG) • Prof. Hélio F. Dos Santos (DQ-UFJF) • Prof. Gotthard Seifert (TU-Dresden) • Prof. Thomas Heine (TU-Dresden) • Dr. Serguei Patchkovskii (NRC-Canada)
Grupo de Pesquisa em Química Inorgânica Teórica - GPQIT • Team: • Dr. Heitor Avelino de Abreu (CNPq) • Antonio Noronha (PhD Student) • Augusto Faria Oliveira (PhD Student) • Luciana Guimarães (PhD Student) • Guilherme Ferreira (IC) • Conny Cerai (IC) • Danniel Brandão (IC) • Leonardo R. R. de Oliveira (IC)
Support • UFMG • Instituto do Milênio: Água - Uma Visão Mineral(PADCT/CNPq) • CNPq • CAPES • FAPEMIG • PRONEX