Comparison of Implicit and Explicit Solvation Models for Iota-Cyclodextrin Conformation Analysis from Replica Exchange Molecular Dynamics

Large ring cyclodextrins have become increasingly important for drug delivery applications. In this work, we have performed replica-exchange molecular dynamics simulations using both implicit and explicit water solvation models to study the conformational diversity of iota-cyclodextrin containing 14 α-1,4 glycosidic linked d-glucopyranose units (CD14). The new quantifiable calculation methods are proposed to analyze the openness, bending, and twisted conformation of CD14 in terms of circularity, biplanar angle, and one-directional conformation (ODC). CD14 in GB implicit water model (Igb5) was found mostly in an opened conformation with average circularity of 0.39 ± 0.16 and a slight bend with average biplanar angle of 145.5 ± 16.0°. In contrast, CD14 in TIP3P explicit water solvation is significantly twisted with average circularity of 0.16 ± 0.10, while 29.1% are ODCs. In addition, classification of CD14 conformations using a Gaussian mixture model (GMM) shows that 85.0% of all CD14 in implicit water at 300 K correspond to the elliptical conformation, in contrast to 82.3% in twisted form in explicit water. GMM clustering also reveals minority conformations of CD14 such as the 8-shape, boat-form, and twisted conformations. This work provides fundamental insights into CD14 conformation, influence of solvation models, and also proposes new quantifiable analysis techniques for molecular conformation studies in the future.

[1]  Walter Rocchia,et al.  Application of conformational clustering in protein-ligand docking. , 2012, Methods in molecular biology.

[2]  K. Harata,et al.  The Structure of the Cyclodextrin Complex. XX. Crystal Structure of Uncomplexed Hydrated γ-Cyclodextrin , 1987 .

[3]  W. Saenger,et al.  An orthorhombic crystal form of cyclohexaicosaose, CA26.32.59 H(2)O: comparison with the triclinic form. , 2001, Carbohydrate research.

[4]  P. Ivanov Conformations of some lower-size large-ring cyclodextrins derived from conformational search with molecular dynamics and principal component analysis , 2012 .

[5]  Florian Sittel,et al.  Robust Density-Based Clustering To Identify Metastable Conformational States of Proteins. , 2016, Journal of chemical theory and computation.

[6]  W. Saenger,et al.  Topography of cyclodextrin inclusion complexes. 15. Crystal and molecular structure of the cyclohexaamylose-7.57 water complex, form III. Four- and six-membered circular hydrogen bonds , 1981 .

[7]  D. van der Spoel,et al.  A temperature predictor for parallel tempering simulations. , 2008, Physical chemistry chemical physics : PCCP.

[8]  K. Harata,et al.  X-Ray Structure of i-Cyclodextrin , 1998 .

[9]  Amedeo Caflisch,et al.  Weighted Distance Functions Improve Analysis of High-Dimensional Data: Application to Molecular Dynamics Simulations. , 2015, Journal of chemical theory and computation.

[10]  H. Ueda Physicochemical Properties and Complex Formation Abilities of Large-Ring Cyclodextrins , 2002 .

[11]  F. Momany,et al.  DFT energy optimization of a large carbohydrate: cyclomaltohexaicosaose (CA-26). , 2012, Journal of Physical Chemistry B.

[12]  V. Hornak,et al.  Modified replica exchange simulation methods for local structure refinement. , 2005, The journal of physical chemistry. B.

[13]  T. Rungrotmongkol,et al.  Conformation study of ɛ-cyclodextrin: Replica exchange molecular dynamics simulations. , 2016, Carbohydrate polymers.

[14]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[15]  U. Hansmann Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[16]  T. Fujiwara,et al.  Structure of δ-Cyclodextrin 13.75H2O , 1990 .

[17]  W. Zimmermann,et al.  Effect of ethanol on the synthesis of large-ring cyclodextrins by cyclodextrin glucanotransferases , 2007 .

[18]  G. Sheldrick,et al.  V-Amylose at atomic resolution: X-ray structure of a cycloamylose with 26 glucose residues (cyclomaltohexaicosaose). , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[19]  Christian Betzel,et al.  Topography of cyclodextrin inclusion complexes, part 20. Circular and flip-flop hydrogen bonding in .beta.-cyclodextrin undecahydrate: a neutron diffraction study , 1984 .

[20]  Wolfram Saenger,et al.  Strain-Induced "Band Flips" in Cyclodecaamylose and Higher Homologues. , 1998, Angewandte Chemie.

[21]  W. Zimmermann,et al.  Effect of the reaction temperature on the transglycosylation reactions catalyzed by the cyclodextrin glucanotransferase from Bacillus macerans for the synthesis of large-ring cyclodextrins , 2004 .

[22]  W. Saenger,et al.  Band-flip and kink as novel structural motifs in α-(1→4)-d-glucose oligosaccharides. Crystal structures of cyclodeca- and cyclotetradecaamylose , 1999 .

[23]  C. Jaime,et al.  Computational study on the intramolecular self-organization of the macrorings of some 'giant' cyclodextrins (CD(n), n = 40, 70, 85, 100). , 2015, Organic & biomolecular chemistry.

[24]  C. Jaime,et al.  Computational study on the conformations of CD38 and inclusion complexes of some lower-size large-ring cyclodextrins , 2014 .

[25]  Hugh Nymeyer,et al.  Atomic Simulations of Protein Folding, Using the Replica Exchange Algorithm , 2004, Numerical Computer Methods, Part D.