Structure and thermodynamics of alpha-, beta-, and gamma-cyclodextrin dimers. Molecular dynamics studies of the solvent effect and free binding energies

The alpha-, beta-, and gamma-cyclodextrin (CyDs) dimers were studied by molecular dynamics (MD) simulations in water as an explicit solvent. The relative stability of dimers and the involved molecular interactions were determined. Three possible starting orientations were considered for the dimers: head-to-head, head-to-tail, and tail-to-tail. MD simulations were performed over a period of 5 ns to ensure the stability of the system for both the CyD dimers and monomers. The MM-PBSA methodology was used to obtain the free binding energy of the dimers and to determine the most stable arrangement for each solvated CyD. In a vacuum, MD simulations provided the head-to-head orientation as the most stable orientation for the three CyDs, while in aqueous solution the, the head-to-tail orientation was found to be the most stable for the alpha-CyD dimer and the tail-to-tail orientation the most stable for the beta- and gamma-CyD dimers.     

Structure and ultrafast dynamics of liquid water: a quantum mechanics/molecular mechanics molecular dynamics simulations study

A quantum mechanics/molecular mechanics molecular dynamics simulation was performed for liquid water to investigate structural and dynamical properties of this peculiar liquid. The most important region containing a central reference molecule and all nearest surrounding molecules (first coordination shell) was treated by Hartree-Fock (HF), post-Hartree-Fock [second-order Moller-Plesset perturbation theory (MP2)], and hybrid density functional B3LYP [Becke's three parameter functional (B3) with the correlation functional of Lee, Yang, and Parr (LYP)] methods. In addition, another HF-level simulation (2HF) included the full second coordination shell. Site to site interactions between oxygen-oxygen, oxygen-hydrogen, and hydrogen-hydrogen atoms of all ab initio methods were compared to experimental data. The absence of a second peak and the appearance of a shoulder instead in the gO-O graph obtained from the 2HF simulation is notable, as this feature has been observed so far only for pressurized or heated water. Dynamical data show that the 2HF procedure compensates some of the deficiency of the HF one-shell simulation, reducing the difference between correlated (MP2) and HF results. B3LYP apparently leads to too rigid structures and thus to an artificial slow down of the dynamics.     

Glucose orientation and dynamics in alpha-, beta-, and gamma-cyclodextrins

We investigate, using molecular dynamics (MD) computer simulations, the conformational behavior of alpha-, beta-, and gamma-cyclodextrins (CDs). Our analysis of a 30 ns trajectory of CD solution dynamics reveals the underlying conformational behaviours of the CDs that explain their relative flexibility. The distributions of the torsion angles related to the glycosidic linkages, P(phi,psi) were calculated for the three CDs. Most noticeable is the limited range in phi torsion rotations compared with psi rotations for all the CDs. This difference between the three CDs is amplified in the motion and dynamics of their glucose monomers when we monitor their orientational and librational positions relative to the macrocyclic mean plane. The relaxation times of the monomers to their equilibrium orientations follow the pattern gamma-CD > alpha-CD > beta-CD. The root-mean-square deviations of the motion of the monomer centers of mass from the mean macrocyclic planes exhibit the same trend.     

Molecular dynamics simulations of structure and dynamics of organic molecular crystals

A set of model compounds covering a range of polarity and flexibility have been simulated using GAFF, CHARMM22, OPLS and MM3 force fields to examine how well classical molecular dynamics simulations can reproduce structural and dynamic aspects of organic molecular crystals. Molecular structure, crystal structure and thermal motion, including molecular reorientations and internal rotations, found from the simulations have been compared between force fields and with experimental data. The MM3 force field does not perform well in condensed phase simulations, while GAFF, CHARMM and OPLS perform very similarly. Generally molecular and crystal structure are reproduced well, with a few exceptions. The atomic displacement parameters (ADPs) are mostly underestimated in the simulations with a relative error of up to 70%. Examples of molecular reorientation and internal rotation, observed in the simulations, include in-plane reorientations of benzene, methyl rotations in alanine, decane, isopropylcyclohexane, pyramidal inversion of nitrogen in amino group and rotation of the whole group around the C-N bond. Frequencies of such dynamic processes were calculated, as well as thermodynamic properties for reorientations in benzene and alanine. We conclude that MD simulations can be used for qualitative analysis, while quantitative results should be taken with caution. It is important to compare the outcomes from simulations with as many experimental quantities as available before using them to study or quantify crystal properties not available from experiment.     

Computational studies on pseudorotaxanes by molecular dynamics and free energy perturbation simulations

A computational scheme that comprises the utilization of the AMBER force field with RESP charges and an explicit solvent model for acetonitrile proved to be useful for studying the structures and energetics of pseudorotaxanes of benzidine and 4,4'-biphenol with cyclobis(paraquat-p-phenylene). The scheme can be further utilized for modeling [2]rotaxanes.     

Molecular electrostatic potentials and hydrogen bonding in alpha-, beta-, and gamma-cyclodextrins

Cyclodextrins (CDs) are cyclic oligomers of glucose having the toroid of sugars elaborating a central cavity of varying size depending on the number of glucoses. The central hydrophobic cavity of CD shows a binding affinity toward different guest molecules, which include small substituted benzenes to long chain surfactant molecules leading to a variety of inclusion complexes when the size and shape complementarity of host and guest are compatible. Further, interaction of guest molecules with the outer surface of alpha-CD has also been observed. Primarily it is the electrostatic interactions that essentially constitute a driving force for the formation of inclusion complexes. To gain insights for these interactions, the electronic structure and the molecular electrostatic potentials in alpha-, beta-, and gamma-CDs are derived using the hybrid density functional theory employing the three-parameter exchange correlation functional due to Becke, Lee, Yang, and Parr (B3LYP). The present work demonstrates how the topography of the molecular electrostatic potential (MESP) provides a measure of the cavity dimensions and understanding of the hydrogen-bonded interactions involving primary and secondary hydroxyl groups. In alpha-CD, hydrogen-bonded interactions between primary -OH groups engender a "cone-like" structure, while in beta- or gamma-CD the interactions from the primary -OH with ether oxygen in glucose ring facilitates a "barrel-like" structure. Further, the strength of hydrogen-bonded interactions of primary -OH groups follows the rank order alpha-CD > beta-CD > gamma-CD, while the secondary hydrogen-bonded interactions exhibit a reverse trend. Thus weak hydrogen-bonded interactions prevalent in gamma-CD manifest in shallow MESP minima near hydroxyl oxygens compared to those in alpha- or beta-CD. Furthermore, electrostatic potential topography reveals that the guest molecule tends to penetrate inside the cavity forming the inclusion complex in beta- or gamma-CD.     

Experimental (NMR) and Computational (MD) Studies on the Inclusion Complexes of 1-Bromoadamantane with alpha-, beta-, and gamma-Cyclodextrin

The inclusion complexes between the most commonly used cyclodextrins (alpha-, beta-, and gamma-CD) and 1-bromoadamantane were prepared and studied experimentally by NMR methods and by molecular dynamics simulations (AMBER force field) with solvation. The NMR results suggest host/guest ratios of 2:1, 1:1, and 1:1 for the complexes with alpha-, beta-, and gamma-cyclodextrin, respectively, as well as defined geometries for the complexes. Averaged geometrical data from the molecular dynamics simulations agree with the complexation geometries deduced experimentally.     

Molecular motions of alpha-cyclodextrin on a dodecyl chain studied by molecular dynamics simulations

Motions of an alpha-cyclodextrin (alpha-CD) molecule on a dodecyl chain adopting the all-trans conformation were investigated in the presence of water by molecular dynamics simulations with CVFF force fields, where the trimethylammonium group of dodecyltrimethylammonium bromide (DTAB) is protruded outside the secondary hydroxyl rim of alpha-CD (the secondary-in structure). The alpha-CD molecule shuttled rapidly on the chain without decomplexation. This rapid motion is consistent with the NMR data. The plane formed by 6 O4 atoms of alpha-CD is most populated between the C6 and C7 atoms of DTAB. This structure is very close to that estimated by NMR. The alpha-CD molecule underwent a restricted rotation in a range of 60 degrees with regard to the plane of the dodecyl chain: this plane at the most population is middle between the two diagonal lines of the normal hexagon formed by 6 O4 atoms of alpha-CD. The published NMR data were reanalyzed in terms of the rotation angle, and a slightly better structure was obtained. The distortion of the alpha-CD cavity from the normal hexagon was decreased upon complex formation with DTAB. The deviation of the center of alpha-CD from the center of the dodecyl chain predicted by molecular dynamics simulations is consistent with the NMR data. The secondary-in structure is energetically more stable than the primary-in structure, as calculated by molecular mechanics with CVFF and Amber force fields. This result is consistent with the NMR data. Molecular dynamics simulations were also carried out for the primary-in structure. Some of the results are close to those of the secondary-in structure.     

Atomistic modeling of crystal-to-amorphous transition and associated kinetics in the Ni-Nb system by molecular dynamics simulations

With the aid of ab initio calculations, an n-body potential of the Ni-Nb system is constructed under the Finnis-Sinclair formalism and the constructed potential is capable of not only reproducing some static physical properties but also revealing the atomistic mechanism of crystal-to-amorphous transition and associated kinetics. With application of the constructed potential, molecular dynamics simulations using the solid solution models reveal that the physical origin of crystal-to-amorphous transition is the crystalline lattice collapsing while the solute atoms are exceeding the critical solid solubilities, which are determined to be 19 atom % Ni and 13 atom % Nb for the Nb- and Ni-based solid solutions, respectively. It follows that an intrinsic glass-forming ability of the Ni-Nb system is within 19-87 atom % Ni, which matches well with that observed in ion beam mixing/solid-state reaction experiments. Simulations using the Nb/Ni/Nb (Ni/Nb/Ni) sandwich models indicate that the amorphous layer at the interfaces grows in a layer-by-layer mode and that, upon dissolving solute atoms, the Ni lattice approaches and exceeds its critical solid solubility faster than the Nb lattice, revealing an asymmetric behavior in growth kinetics. Moreover, an energy diagram is obtained by computing the energetic sequence of the Ni(x)Nb(100)(-)(x) alloy in fcc, bcc, and amorphous structures, respectively, over the entire composition range, and the diagram could serve as a guide for predicting the metastable alloy formation in the Ni-Nb system.     

Self-assembling cyclic peptides: molecular dynamics studies of dimers in polar and nonpolar solvents

The self-assembly of cyclic D,L-alpha-peptides into hollow nanotubes is a crucial mechanistic step in their application as antibacterial and drug-delivery agents. To understand this process, molecular dynamics (MD) simulations were performed on dimers of cyclic peptides formed from cyclo [(-L-Trp-D-N-MeLeu-)4-]2 and cyclo [(-L-Trp-D-Leu-)4-]2 subunits in nonpolar (nonane) and polar (water) solvent. The dimers were observed to be stable only in nonpolar solvent over the full 10 ns length of the MD trajectory. The behavior of the dimers in different solvents is rationalized in terms of the intersubunit hydrogen bonding, hydrogen bonding with the solvent, and planarity of the rings. It is shown that the phi and psi dihedral angles of a single uncapped ring in nonane lie in the beta-sheet region of the Ramachandran plot, and the ring stays in a flat conformation. Steered MD (SMD) simulations based on Jarzynski's equality were performed to obtain the potential of mean force as a function of the distance between the two rings of the capped dimer in nonane. It is also shown that a single peptide subunit prefers to reside close to the nonane/water interface rather than in bulk solvent because of the amphiphilic character of the peptide ring. The present MD results build the foundation for using MD simulations to study the mechanism of the formation of cyclic peptide nanotubes in lipid bilayers.     

Thermodiffusion in model nanofluids by molecular dynamics simulations

In this work, a new algorithm is proposed to compute single particle (infinite dilution) thermodiffusion using nonequilibrium molecular dynamics simulations through the estimation of the thermophoretic force that applies on a solute particle. This scheme is shown to provide consistent results for model nanofluids in the liquid state (spherical nonmetallic nanoparticles+Lennard-Jones fluid) where it appears that thermodiffusion amplitude, as well as thermal conductivity, decreases with nanoparticle concentration. Then, by changing the nature of the nanoparticle (size, mass, and internal stiffness) and that of the solvent (quality and viscosity), various trends are exhibited. In all cases, the single particle thermodiffusion is positive, i.e., the nanoparticle tends to migrate toward the cold area. The single particle thermal diffusion coefficient is shown to be independent of the size of the nanoparticle (diameter of 0.8-4 nm), whereas it increases with the quality of the solvent and is inversely proportional to the viscosity of the fluid. In addition, this coefficient is shown to be independent of the mass of the nanoparticle and to increase with the stiffness of the nanoparticle internal bonds. Besides, for these configurations, the mass diffusion coefficient behavior appears to be consistent with a Stokes-Einstein-like law.     

Atomic and molecular quantum mechanics by the path integral molecular dynamics method

The quantum path integral molecular dynamics method was applied to studies of excess electron localization by a Na⁺ ion and by a NaCl molecule. Spatial and energetic characterization of the ground state of the excess electron compare favorably with results of model potential calculations for Na and with SCF Cl calculations for NaCl⁻.     

Molecular modeling studies of the reactions of phenoxy radical dimers: pathways to dibenzofurans

The B3LYP hybrid density functional computational technique was applied to describe the sequence of phenoxy radicals coupling reactions leading to the formation of dibenzofurans. Reaction kinetic parameters were estimated for key reactions. Aromatization of bis-keto dimers of phenoxy radicals followed by intermediate dehydration or dehydroxylation was demonstrated to be a strongly stereoselective process. While the S,S-diastereomer of the ortho-C//ortho-C keto dimer forms (o,o')-dihydroxybiphenyl, a known dibenzofuran intermediate, via inter-ring hydrogen transfer reaction, the less stable R,S-stereoisomer can easily be transformed into another 5-hydroxyl-4,5-cyclohexadiene-2,3-benzofuran intermediate that provides an energetically more favorable pathway for formation of dibenzofuran. The possible channels of radical-chain processes that convert these intermediates to dibenzofuran and polychlorinated dibenzofurans are discussed.     

A new set of molecular mechanics parameters for hydroxyproline and its use in molecular dynamics simulations of collagen-like peptides

Recently, the importance of proline ring pucker conformations in collagen has been suggested in the context of hydroxylation of prolines. The previous molecular mechanics parameters for hydroxyproline, however, do not reproduce the correct pucker preference. We have developed a new set of parameters that reproduces the correct pucker preference. Our molecular dynamics simulations of proline and hydroxyproline monomers as well as collagen-like peptides, using the new parameters, support the theory that the role of hydroxylation in collagen is to stabilize the triple helix by adjusting to the right pucker conformation (and thus the right phi angle) in the Y position.     

Investigation of chiral recognition by molecular micelles with molecular dynamics simulations

Molecular dynamics simulations were used to characterize the binding of the chiral drugs chlorthalidone and lorazepam to the molecular micelle poly-(sodium undecyl-(L)-leucine-valine). The project’s goal was to characterize the nature of chiral recognition in capillary electrophoresis separations that use molecular micelles as the chiral selector. The shapes and charge distributions of the chiral molecules investigated, their orientations within the molecular micelle chiral binding pockets, and the formation of stereoselective intermolecular hydrogen bonds with the molecular micelle were all found to play key roles in determining where and how lorazepam and chlorthalidone enantiomers interacted with the molecular micelle.     

Molecular dynamics simulations of porphyrin-dendrimer systems: toward modeling electron transfer in solution

We have performed computational simulations of porphyrin-dendrimer systems--a cationic porphyrin electrostatically associated to a negatively charged dendrimer--using the method of classical molecular dynamics (MD) with an atomistic force field. Previous experimental studies have shown a strong quenching effect of the porphyrin fluorescence that was assigned to electron transfer (ET) from the dendrimer's tertiary amines (Paulo, P. M. R.; Costa, S. M. B. J. Phys. Chem. B 2005, 109, 13928). In the present contribution, we evaluate computationally the role of the porphyrin-dendrimer conformation in the development of a statistical distribution of ET rates through its dependence on the donor-acceptor distance. We started from simulations without explicit solvent to obtain trajectories of the donor-acceptor distance and the respective time-averaged distributions for two dendrimer sizes and different initial configurations of the porphyrin-dendrimer pair. By introducing explicit solvent (water) in our simulations, we were able to estimate the reorganization energy of the medium for the systems with the dendrimer of smaller size. The values obtained are in the range 0.6-1.5 eV and show a linear dependence with the inverse of the donor-acceptor distance, which can be explained by a two-phase dielectric continuum model taking into account the medium heterogeneity provided by the dendrimer organic core. Dielectric relaxation accompanying ET was evaluated from the simulations with explicit solvent showing fast decay times of some tens of femtoseconds and slow decay times in the range of hundreds of femtoseconds to a few picoseconds. The variations of the slow relaxation times reflect the heterogeneity of the dendrimer donor sites which add to the complexity of ET kinetics as inferred from the experimental fluorescence decays.     

The cationic halochromism of phenolate betaines: Molecular dynamics and quantum mechanics studies

The cationic halochromism of phenolate betaines was reproduced with the aid of a simple theoretical model, by calculation of the longest wavelength transition energies of supermolecules obtained by positioning a cation Mⁿ⁺ at a variable distance from the oxygen atom of the dye. The theoretical results were compared with experimental data for three systems, Reichardt’s betaine 1, Brooker’s merocyanine 2 and the N-methyl-8-oxyquinolinium dye 3. The model was validated by molecular dynamics simulations of solutions of dye 3, in methanol and DMSO, in the presence of variable concentrations of Na⁺.