Glycosidic linkage conformation of methyl-alpha-mannopyranoside

We study the preferred conformation of the glycosidic linkage of methyl-alpha-mannopyranoside in the gas phase and in aqueous solution. Results obtained utilizing Car-Parrinello molecular dynamics (CPMD) simulations are compared to those obtained from classical molecular dynamics (MD) simulations. We describe classical simulations performed with various water potential functions to study the impact of the chosen water potential on the predicted conformational preference of the glycosidic linkage of the carbohydrate in aqueous solution. In agreement with our recent studies, we find that results obtained with CPMD simulations differ from those obtained from classical simulations. In particular, this study shows that the trans (t) orientation of the glycosidic linkage of methyl-alpha-mannopyranoside is preferred over its gauche anticlockwise (g-) orientation in aqueous solution. CPMD simulations indicate that this preference is due to intermolecular hydrogen bonding with surrounding water molecules, whereas no such information could be demonstrated by classical MD simulations. This study emphasizes the importance of ab initio MD simulations for studying the structural properties of carbohydrates in aqueous solution.     

Competition of Intra- and Intermolecular Forces in Anthraquinone and Its Selected Derivatives

Intra- and intermolecular forces competition was investigated in the 9,10-anthraquinone (1) and its derivatives both in vacuo and in the crystalline phase. The 1,8-dihydroxy-9,10-anthraquinone (2) and 1,8-dinitro-4,5-dihydroxy-anthraquinone (3) contain Resonance-Assisted Hydrogen Bonds (RAHBs). The intramolecular hydrogen bonds properties were studied in the electronic ground and excited states employing Møller-Plesset second-order perturbation theory (MP2), Density Functional Theory (DFT) method in its classical formulation as well as its time-dependent extension (TD-DFT). The proton potential functions were obtained via scanning the OH distance and the dihedral angle related to the OH group rotation. The topological analysis was carried out on the basis of theories of Atoms in Molecules (AIM—molecular topology, properties of critical points, AIM charges) and Electron Localization Function (ELF—2D maps showing bonding patterns, calculation of electron populations in the hydrogen bonds). The Symmetry-Adapted Perturbation Theory (SAPT) was applied for the energy decomposition in the dimers. Finally, Car–Parrinello molecular dynamics (CPMD) simulations were performed to shed light onto bridge protons dynamics upon environmental influence. The vibrational features of the OH stretching were revealed using Fourier transformation of the autocorrelation function of atomic velocity. It was found that the presence of OH and NO2 substituents influenced the geometric and electronic structure of the anthraquinone moiety. The AIM and ELF analyses showed that the quantitative differences between hydrogen bonds properties could be neglected. The bridged protons are localized on the donor side in the electronic ground state, but the Excited-State Intramolecular Proton Transfer (ESIPT) was noticed as a result of the TD-DFT calculations. The hierarchy of interactions determined by SAPT method indicated that weak hydrogen bonds play modifying role in the organization of these crystal structures, but primary ordering factor is dispersion. The CPMD crystalline phase results indicated bridged proton-sharing in the compound 2.     

Fermi resonance and strong anharmonic effects in the absorption spectra of the ν-OH (ν-OD) vibration of solid H- and D-benzoic acid

The vibrational spectra of polycrystalline benzoic acid (BA) and its deuterated derivative were studied over the wide frequency region 4000–10 cm⁻¹ by IR and Raman methods. A theoretical analysis of the hydrogen bond frequency region and calculations at the B3LYP/6-311++G(2d, 2p) level for the benzoic acid cyclic dimer in the gas phase were made. In order to study the dynamics of proton transfer two formalisms were applied: Car–Parrinello Molecular Dynamics (CPMD) and Path Integrals Molecular Dynamics (PIMD). It was shown that the experimentally observed very broad ν-OH band absorption is the result of complex anharmonic interaction: Fermi resonance between the OH-stretching and bending vibrations and strong interaction of the ν-OH stretching with the low frequency phonons. The theoretical analysis in the framework of such an approach gave a good correlation with experiment. From the CPMD calculations it was confirmed that the O–HO bridge is not rigid, with the OO distance being described by a large amplitude motion. For the benzoic acid dimer we observed stepwise (asynchronous) proton transfer.     

Revealing Intra- and Intermolecular Interactions Determining Physico-Chemical Features of Selected Quinolone Carboxylic Acid Derivatives

The intra- and intermolecular interactions of selected quinolone carboxylic acid derivatives were studied in monomers, dimers and crystals. The investigated compounds are well-recognized as medicines or as bases for further studies in drug design. We employed density functional theory (DFT) in its classical formulation to develop gas-phase and solvent reaction field (PCM) models describing geometric, energetic and electronic structure parameters for monomers and dimers. The electronic structure was investigated based on the atoms in molecules (AIM) and natural bond orbital (NBO) theories. Special attention was devoted to the intramolecular hydrogen bonds (HB) present in the investigated compounds. The characterization of energy components was performed using symmetry-adapted perturbation theory (SAPT). Finally, the time-evolution methods of Car–Parrinello molecular dynamics (CPMD) and path integral molecular dynamics (PIMD) were employed to describe the hydrogen bond dynamics as well as the spectroscopic signatures. The vibrational features of the O-H stretching were studied using Fourier transformation of the autocorrelation function of atomic velocity. The inclusion of quantum nuclear effects provided an accurate depiction of the bridged proton delocalization. The CPMD and PIMD simulations were carried out in the gas and crystalline phases. It was found that the polar environment enhances the strength of the intramolecular hydrogen bonds. The SAPT analysis revealed that the dispersive forces are decisive factors in the intermolecular interactions. In the electronic ground state, the proton-transfer phenomena are not favourable. The CPMD results showed generally that the bridged proton is localized at the donor side, with possible proton-sharing events in the solid-phase simulation of stronger hydrogen bridges. However, the PIMD enabled the quantitative estimation of the quantum effects inclusion—the proton position was moved towards the bridge midpoint, but no qualitative changes were detected. It was found that the interatomic distance between the donor and acceptor atoms was shortened and that the bridged proton was strongly delocalized.     

Elucidating the vibrational spectra of hydrogen-bonded aggregates in solution: electronic structure calculations with implicit solvent and first-principles molecular dynamics simulations with explicit solvent for 1-hexanol in n-hexane

Fourier transform infrared spectroscopy is a popular method for the experimental investigation of hydrogen-bonded aggregates, but linking spectral information to microscopic information on aggregate size distribution and aggregate architecture is an arduous task. Static electronic structure calculations with an implicit solvent model, Car-Parrinello molecular dynamics (CPMD) using the Becke-Lee-Yang-Parr (BLYP) exchange and correlation energy functionals and classical molecular dynamics simulations for the all-atom version of the optimized parameters for liquid simulations (OPLS-AA) force field were carried out for an ensemble of 1-hexanol aggregates solvated in n-hexane. The initial configurations for these calculations were size-selected from a distribution of aggregates obtained from a large-scale Monte Carlo simulation. The vibrational spectra computed from the static electronic structure calculations for monomers and dimers and from the CPMD simulations for aggregates up to pentamers demonstrate the extent of the contribution of dangling or nondonating hydroxyl groups found in linear and branched aggregates to the "monomeric" peak. Furthermore, the computed spectra show that there is no simple relationship between peak shift and aggregate size nor architecture, but the effect of hydrogen-bond cooperativity is shown to differentiate polymer-like (cooperative) and dimer-like (noncooperative) hydrogen bonds in the vibrational spectrum. In contrast to the static electronic structure calculations and the CPMD simulations, the classical molecular dynamics simulations greatly underestimate the vibrational peak shift due to hydrogen bonding.     

Car-Parrinello simulation of an O-H stretching envelope and potential of mean force of an intramolecular hydrogen bonded system: application to a Mannich base in solid state and in vacuum

Car-Parrinello molecular dynamics (CPMD) study was performed for an anharmonic system-an intramolecularly hydrogen bonded Mannich-base-type compound, 4,5-dimethyl-2(N,N-dimethylaminemethyl)phenol, to investigate the vibrational spectrum associated with the O-H stretching. Calculations were carried out for the solid state and for an isolated molecule. The classical CPMD simulation was performed and then the proton potential snapshots were extracted from the trajectory. The vibrational Schrodinger equation for the snapshots was solved numerically, and the (O-H) envelope was calculated as a superposition of the 0-->1 transitions. The potential of mean force for the proton stretching mode was calculated from the proton vibrational eigenfunctions and eigenvalues incorporating statistical sampling, nuclear quantum effects, and effects of the environment. Perspectives for application of the presented methodology in the computational support of biocatalysis are given in the study.     

Preferred conformation of the glycosidic linkage of methyl-beta-mannose

The conformational preference of the glycosidic linkage of methyl-beta-mannose was studied in the gas phase and in aqueous solution by ab initio calculations, and by molecular dynamics (MD) and Car-Parrinello molecular dynamics (CPMD) simulations. MD simulations were performed with various water potential functions to study the impact of the chosen water potential on the predicted conformational preference of the glycosidic linkage of the carbohydrate in solution. This study shows that the trans (t) orientation of the glycosidic linkage of methyl-beta-mannose is preferred over its gauche clockwise (g+) orientation in solution. CPMD simulations clearly indicate that this preference is due to intermolecular hydrogen bonding with surrounding water molecules, whereas no such information could be demonstrated by MD simulations. This study demonstrates the importance of ab initio molecular dynamics simulations in studying the structural properties of carbohydrate-water interactions.     

Combined ab initio and classical molecular dynamics simulations of alkyl-lithium aggregates in ethereal solutions

Molecular dynamics simulations of organolithium aggregates in solution are reported for the first time. We use a combined quantum/classical force field (the so-called QM/MM approach) and study ethyl-lithium aggregates in dimethyl ether (DME) solvent. The solutes are described at the Density Functional Theory level while solvent molecules are described using molecular mechanics. NVT Molecular Dynamics simulations at 200 K are carried out in the Born–Oppenheimer approximation. After equilibration, the production phase was run for 80 ps (monomer), 40 ps (dimer) and 26 ps (tetramer). The analysis of the results focuses on Li coordination as a function of aggregate size and we show that the total Li coordination number is always 4. No decoordination has been observed along the simulations. Fluctuations of the structures are predicted to be large in some cases and possible implications on reactivity are discussed.     

Car-Parrinello molecular dynamics simulation of Fe 3+ (aq)

The optimized geometry and energetic properties of Fe(D2O)n 3+ clusters, with n = 4 and 6, have been studied with density-functional theory calculations and the BLYP functional, and the hydration of a single Fe 3+ ion in a periodic box with 32 water molecules at room temperature has been studied with Car-Parrinello molecular dynamics and the same functional. We have compared the results from the CPMD simulation with classical MD simulations, using a flexible SPC-based water model and the same number of water molecules, to evaluate the relative strengths and weaknesses of the two MD methods. The classical MD simulations and the CPMD simulations both give Fe-water distances in good agreement with experiment, but for the intramolecular vibrations, the classical MD yields considerably better absolute frequencies and ion-induced frequency shifts. On the other hand, the CPMD method performs considerably better than the classical MD in describing the intramolecular geometry of the water molecule in the first hydration shell and the average first shell...second shell hydrogen-bond distance. Differences between the two methods are also found with respect to the second-shell water orientations. The effect of the small box size (32 vs 512 water molecules) was evaluated by comparing results from classical simulations using different box sizes; non-negligible effects are found for the ion-water distance and the tilt angles of the water molecules in the second hydration shell and for the O-D stretching vibrational frequencies of the water molecules in the first hydration shell.     

Ab initio molecular dynamics study of 45S5 bioactive silicate glass

Bioglass 45S5, the prototype of bioactive melt-quenched silicate glasses, was modeled by means of Car-Parrinello molecular dynamics (CPMD) simulations. Although long-range structural properties cannot be modeled by using this ab initio approach, the accuracy of CPMD simulations is exploited here to provide insight into the short-range structure and to analyze vibrational and electronic properties of this biomaterial. Detailed structural analysis in the short-range scale provided insight into the local environment of modifier Na and Ca ions: a possible key role of these cations in organizing the glass network by connecting different chains and fragments into specific, rather flexible geometries was proposed. The individual contributions of different species to the vibrational density of states were separated and discussed, allowing the identification of specific features in the vibrational spectrum, such as those related to phosphate groups. The components of the electronic density of states were also analyzed, enabling us to identify correlations between the electronic structure and the structural properties, such as the different bonding character of Si-O bonds involving bridging or nonbridging oxygen atoms.     

Investigation of structure of liquid 2,2,2 trifluoroethanol: neutron diffraction, molecular dynamics, and ab initio quantum chemical study

The molecular conformation and intermolecular H bonding in liquid 2,2,2 trifluoroethanol (TFE) have been studied by neutron diffraction with hydrogen/deuterium isotopic substitution at room temperature. For comparison, conformations of molecules and their dimers in the gas phase have also been calculated, based on the density functional theory. Energies, geometry, and vibrational frequencies of dimers were analyzed. Diffraction data analyzed by the "Monte Carlo determination of g(r)" (MCGR) method resulted in a molecular structure in agreement with the findings from gas phase electron diffraction experiments and density functional calculations. The intermolecular structure functions were compared to the same functions obtained from a molecular dynamics simulation. All of the composite radial distribution functions are in good agreement with the simulation results. According to our calculation the hydrogen-bonded aggregation size is smaller in pure liquid TFE than in pure liquid ethanol.     

Identification of active sites of biomolecules. 1. Methyl-alpha-mannopyranoside and Fe(III)

Car-Parrinello molecular dynamics (CPMD) simulations, DFT chemical reactivity index calculations, and mass spectrometric measurements are combined in an integrated effort to elucidate the details of the coordination of a transition-metal ion to a carbohydrate. The impact of the interaction with the FeIII ion on the glycosidic linkage conformation of methyl-alpha-d-mannopyranoside is studied by classical molecular dynamics (MD) and CPMD simulations. This study shows that FeIII interacts with specific hydroxyl oxygen atoms of the carbohydrate, affecting the ground state carbohydrate conformation. These conformational details are discussed in terms of a set of supporting experiments involving electrospray ionization mass spectrometry, and CPMD simulations clearly indicate that the specific conformational preference is due to intramolecular hydrogen bonding. Classical MD simulations proved insensitive to these important chemical properties. Thus, we demonstrate the importance of chemical reactivity calculations and CPMD simulations in predicting the active sites of biological molecules toward metal cations.     

Computational study on the properties and structure of methyl lactate

A theoretical study on the properties and molecular level structure of the very important green solvent methyl lactate is carried out in the gas phase and methanol and water solutions, with the solvent treated both explicitly and as a continuum. Torsional barriers giving rise to different conformers by rotation of the hydroxyl and methyl groups were analyzed using density functional theory (DFT) to establish the most stable conformer both in gas phase and solution. DFT computations on lactate dimers were also done to study short-range features, and the effect of the surrounding solvent on intra- and intermolecular hydrogen bonding was analyzed according to the polarizable continuum model approach. We have also studied lactate/water and lactate/methanol small clusters together with the corresponding binding energies. Moreover, classical molecular dynamics simulations (MD) were carried out to study medium- and large-range effects at lower computational cost. MD simulations at different pressure and temperature conditions on pure lactate were carried out, and mixtures with water and methanol of different compositions were also studied. Structural information, analyzed through the radial distribution functions, together with dynamic aspects of pure and mixed fluids were considered. The intramolecular hydrogen bonding ability of methyl lactate together with the possibility of homo- and hetero-intermolecular association determines the behavior of this molecule in pure fluids or in mixed.     

Brownian Dynamics, Molecular Dynamics, and Monte Carlo modeling of colloidal systems

This paper serves as an introductory review of Brownian Dynamics (BD), Molecular Dynamics (MD), and Monte Carlo (MC) modeling techniques. These three simulation methods have proven to be exceptional investigative solutions for probing discrete molecular, ionic, and colloidal motions at their basic microscopic levels. The review offers a general study of the classical theories and algorithms that are foundational to Brownian Dynamics, Molecular Dynamics, and Monte Carlo simulations. Important topics of interest include fundamental theories that govern Brownian motion, the Langevin equation, the Verlet algorithm, and the Metropolis method. Brownian Dynamics demonstrates advantages over Molecular Dynamics as pertaining to the issue of time-scale separation. Monte Carlo methods exhibit strengths in terms of ease of implementation. Hybrid techniques that combine these methods and draw from these efficacies are also presented. With their rigorous microscopic approach, Brownian Dynamics, Molecular Dynamics, and Monte Carlo methods prove to be especially viable modeling methods for problems with challenging complexities such as high-level particle concentration and multiple particle interactions. These methods hold promising potential for effective modeling of transport in colloidal systems.     

Modeling fluid diffusion using the lattice density functional theory approach: counterdiffusion in an external field

The dependence of the diffusivity on temperature, pressure, and composition is not understood well; consequently, data is preferred significantly over correlations in most practical situations. Even in dilute gases, the contributions of attractions and repulsions to the diffusivity are difficult to understand on a molecular level without performing simulations. We have developed a Lattice Density Functional Theory (LDFT) approach for modeling diffusion to supplement existing methods that are very rigorous but computationally demanding. The LDFT approach is analogous to the van der Waals equation in how it accounts for molecular interactions in that it has first-order corrections to ideal behavior; it is an extension of the Equilibrium LDFT for adsorption and phase behavior. In this work, the LDFT approach is presented and demonstrated by modeling the problem of color counterdiffusion in an externally-applied potential field. This potential field, in combination with the intermolecular potential function, creates a diffusion regime in which repulsions cause oscillations in the density profile. Using the LDFT approach, the oscillations were described and attributed to nearest-neighbor and next nearest-neighbor interactions. The LDFT approach gives qualitative and quantitative agreement with dual control-volume Grand Canonical Molecular Dynamics simulations.     

Intermolecular vibrations and diffusive orientational dynamics of Cs condensed ring aromatic molecular liquids

The ultrafast dynamics, including the intermolecular vibrations and the diffusive orientational dynamics, of the neat C(s) symmetry condensed ring aromatic molecular liquids benzofuran, 1-fluoronaphtalene, and quinoline were investigated for the first time by means of femtosecond Raman-induced Kerr effect spectroscopy. To understand the features of these C(s) condensed ring aromatic molecular liquids, reference singular aromatic molecular liquids, furan, fluorobenzene, pyridine, and benzene, were also studied. High quality low-frequency Kerr spectra of the aromatic molecular liquids were obtained by Fourier-transform deconvolution analysis of the measured Kerr transients. The Kerr spectra of the C(s) condensed ring aromatic molecular liquids are bimodal, as are those of the reference singular aromatic molecular liquids. The first moment of the intermolecular vibrational spectrum and the peak frequencies of the high- and low-frequency components in the broad spectrum band were compared with their molecular properties such as the rotational constants, molecular weight, and intermolecular (bimolecular) force. The comparisons show that the molecular volume (related to molecular weight and rotational constants) is a dominant property for the characteristic frequency of the entire intermolecular vibrational spectrum. The observed intramolecular vibrational modes in the Kerr spectra of the aromatic molecular liquids were also assigned on the basis of the ab initio quantum chemical calculation results. In their picosecond diffusive orientational dynamics, the slowest relaxation time constant for both the condensed ring and singular aromatic molecular liquids can be accounted for by the simple Stokes-Einstein-Debye hydrodynamic model.     

Spin-orbit density functional and ab initio study of HgX(n) (X=F, Cl, Br, and I; n=1, 2, and 4)

Quantum chemical calculations of HgX(n) (X=F, Cl, Br, and I; n=1, 2, and 4) in the gas phase are performed using the density functional theory (DFT), two-component spin-orbit (SO) DFT, and high-level ab initio method with relativistic effective core potentials (RECPs). Molecular geometries, vibrational frequencies, and various thermochemical energies are calculated and compared with available experimental results. We assess the performances of DFT functionals for calculating various molecular properties. The PBE0 functional is generally reasonable for the molecular geometries and the vibrational frequencies, but the M06 functional is more appropriate for estimating thermochemical energies. Both shape-consistent and energy-consistent RECPs correctly describe the SO effect.