Simulations of vibrational spectra from classical trajectories: calibration with ab initio force fields

An algorithm allowing simulating vibrational spectra from classical time-dependent trajectories was applied for infrared absorption, vibrational circular dichroism, Raman, and Raman optical activity of model harmonic systems. The implementation of the theory within the TINKER molecular dynamics (MD) program package was tested with ab initio harmonic force fields in order to determine the feasibility for more extended MD simulations. The results suggest that sufficiently accurate frequencies can be simulated with integration time steps shorter than about 0.5 fs. For a given integration time step, lower vibrational frequencies ( approximately 0-2000 cm(-1)) could be reproduced with a higher accuracy than higher-frequency vibrational modes (e.g., O-H and C-H stretching). In principle, the algorithm also provides correct intensities for ideal systems. In applied simulations, however, the intensity profiles are affected by an unrealistic energy distribution between normal modes and a slow energy relaxation. Additionally, the energy fluctuations may cause weakening of the intensities on average. For ab initio force fields, these obstacles could be overcome by an arbitrary normal mode energy correction. For general MD simulations, averaging of many shorter MD trajectories started with randomly distributed atomic velocities provided the best spectral shapes. alpha-pinene, D-gluconic acid, formaldehyde dimer, and the acetylprolineamide molecule were used in the tests.     

Nitrile and thiocyanate IR probes: molecular dynamics simulation studies

Nitrile- and thiocyanate-derivatized amino acids have been found to be useful IR probes for investigating their local electrostatic environments in proteins. To shed light on the CN stretch frequency shift and spectral lineshape change induced by interactions with hydrogen-bonding solvent molecules, we carried out both classical and quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations for MeCN and MeSCN in water. These QM/MM and conventional force field MD simulation results were found to be inconsistent with the experimental results as well as with the high-level ab initio calculation results of MeCN-water and MeSCN-water potential energies. Thus, a new set of atomic partial charges of MeCN and MeSCN is obtained. By using the MD simulation trajectories and the electrostatic potential model recently developed, the CN and SCN stretching mode frequency trajectories were obtained and used to simulate the IR spectra. The C[Triple Bond]N frequency blueshifts of MeCN and MeSCN in water are estimated to be 9.0 and 1.9 cm(-1), respectively, in comparison with those of gas phase values. These values are found to be in reasonable agreement with the experimentally measured IR spectra of MeCN, MeSCN, beta-cyano-L-alanine, and cyanylated cysteine in water and other polar solvents.     

Experimental probing and modeling of key sorbent-solute interactions of norephedrine enantiomers with polysaccharide-based chiral stationary phases

The key interactions of a chiral solute, norephedrine or 2-amino-1-phenyl-1-propanol (PPA), with three commercially important polysaccharide-based chiral stationary phases, amylose Tris(3,5-dimethylphenylcarbamate) (ADMPC), amylose Tris((S)-alpha-methylbenzylcarbamate) (ASMBC) and cellulose Tris(3,5-dimethylphenylcarbamate) (CDMPC) are studied in detail using different experimental techniques and molecular simulations. The HPLC retention factors of the enantiomers of PPA in n-hexane/2-propanol (90/10, v/v) at 298 K vary significantly with these sorbents. The enantioselectivities of -PPA versus +PPA are 2.4, 1.0, and 0.8 (reversal in the elution order), respectively. The observed changes in the wavenumbers and the intensities of the amide bands of these polymers in the attenuated total reflection-infrared spectroscopy (ATR-IR) spectra upon absorption of each enantiomer are different. The IR wavenumbers, and the H-bonding interaction energies of the polymer side chains with each enantiomer (polymer-solute) in four different binding configurations are estimated and ranked using the density functional theory (the DFT/B3LYP/6-311+g(d,p) level of theory). X-ray diffraction (XRD) results show that the polymer crystallinity increases significantly upon absorption of each enantiomer. The helical pitches and the inter-rod packing for these polymers are inferred from the XRD results and incorporated into the molecular dynamics (MD) simulations. The elution orders predicted for the enantiomers of PPA using MD simulations of the polymer-PPA binary systems are consistent with the chromatography results. The enantioselectivity observed in ADMPC is hypothesized to be due to having three simultaneous interactions (two H-bond and one pi-pi) of the polymer with -PPA versus two interactions (one H-bond and one pi-pi) with +PPA.     

Shielding of ionic interactions by sulfur dioxide in an ionic liquid

The effect of adding SO2 on the structure and dynamics of 1-butyl-3-methylimidazolium bromide (BMIBr) was investigated by low-frequency Raman spectroscopy and molecular dynamics (MD) simulations. The MD simulations indicate that the long-range structure of neat BMIBr is disrupted resulting in a liquid with relatively low viscosity and high conductivity, but strong correlation of ionic motion persists in the BMIBr-SO2 mixture due to ionic pairing. Raman spectra within the 5相似文献   

In silico models for the human alpha4beta2 nicotinic acetylcholine receptor

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is one of the most widely expressed nAChR subtypes in the brain. Its subunits have high sequence identity (54 and 46% for alpha4 and beta2, respectively) with alpha and beta subunits in Torpedo nAChR. Using the known structure of the Torpedo nAChR as a template, the closed-channel structure of the alpha4beta2 nAChR was constructed through homology modeling. Normal-mode analysis was performed on this closed structure and the resulting lowest frequency mode was applied to it for a "twist-to-open" motion, which increased the minimum pore radius from 2.7 to 3.4 A and generated an open-channel model. Nicotine could bind to the predicted agonist binding sites in the open-channel model but not in the closed one. Both models were subsequently equilibrated in a ternary lipid mixture via extensive molecular dynamics (MD) simulations. Over the course of 11 ns MD simulations, the open channel remained open with filled water, but the closed channel showed a much lower water density at its hydrophobic gate comprised of residues alpha4-V259 and alpha4-L263 and their homologous residues in the beta2 subunits. Brownian dynamics simulations of Na+ permeation through the open channel demonstrated a current-voltage relationship that was consistent with experimental data on the conducting state of alpha4beta2 nAChR. Besides establishment of the well-equilibrated closed- and open-channel alpha4beta2 structural models, the MD simulations on these models provided valuable insights into critical factors that potentially modulate channel gating. Rotation and tilting of TM2 helices led to changes in orientations of pore-lining residue side chains. Without concerted movement, the reorientation of one or two hydrophobic side chains could be enough for channel opening. The closed- and open-channel structures exhibited distinct patterns of electrostatic interactions at the interface of extracellular and transmembrane domains that might regulate the signal propagation of agonist binding to channel opening. A potential prominent role of the beta2 subunit in channel gating was also elucidated in the study.     

Effect of nanotube-length on the transport properties of single-file water molecules: transition from bidirectional to unidirectional

We use molecular dynamics (MD) simulations to study the transport of single-file water molecules through carbon nanotubes (CNTs) with various lengths in an electric field. Most importantly, we find that even the water dipoles inside the CNT are maintained along the field direction, a large amount of water molecules can still transport against the field direction for short CNTs, leading to a low unidirectional transport efficiency (η). As the CNT length increases, the efficiency η will increase remarkably, and achieves the maximum value of 1.0 at or exceeding a critical CNT length. Consequently, the transition from bidirectional to unidirectional transport is observed and is found to be relevant to thermal fluctuations of the two reservoirs, which is explored by the interaction between water molecules inside and outside the CNT. We also find that the water flow vs CNT length follows an exponential decay of f ～ exp?(- L/L(0)), and the average translocation time of individual water molecules yields to a power law of τ(trans) ～ L(υ), where L(0) and ν are constant and slightly depend on the field strength. We further compare our results with the continuous-time random-walk (CTRW) model and find that the water flow can also be described by a power law of f ～ L(-μ) modified from CTRW. Our results provide some new physical insights into the biased transport of single-file water molecules, which show the feasibility of using CNTs with any length to pump water in an electric field. The mechanism is important for designing efficient nanofluidic apparatuses.     

Molecular dynamics integration and molecular vibrational theory. III. The infrared spectrum of water

The new symplectic molecular dynamics (MD) integrators presented in the first paper of this series were applied to perform MD simulations of water. The physical properties of a system of flexible TIP3P water molecules computed by the new integrators, such as diffusion coefficients, orientation correlation times, and infrared (IR) spectra, are in good agreement with results obtained by the standard method. The comparison between the new integrators' and the standard method's integration time step sizes indicates that the resulting algorithm allows a 3.0 fs long integration time step as opposed to the standard leap-frog Verlet method, a sixfold simulation speed-up. The accuracy of the method was confirmed, in particular, by computing the IR spectrum of water in which no blueshifting of the stretching normal mode frequencies is observed as occurs with the standard method.     

Using MD snapshots in ab initio and DFT calculations: OH vibrations in the first hydration shell around Li+(aq)

The average OH stretching vibrational frequency for the water molecules in the first hydration shell around a Li(+) ion in a dilute aqueous solution was calculated by a hybrid molecular dynamics + quantum-mechanical ("MD + QM") approach. Using geometry configurations from a series of snapshots from an MD simulation, the anharmonic, uncoupled OH stretching frequencies were calculated for 100 first-shell OH oscillators at the B3LYP and HF/6-31G(d,p) levels of theory, explicitly including the first shell and the relevant second shell water molecules into charge-embedded supermolecular QM calculations. Infrared intensity-weighting of the density-of-states (DOS) distributions by means of the squared dipole moment derivatives (which vary by a factor of 20 over the OH stretching frequency band at the B3LYP level), changes the downshift from approximately -205 to -275 cm(-1) at the B3LYP level. Explicit inclusion of relevant third-shell water molecules in the supermolecular cluster leads to a further downshift by approximately -30 cm(-1). Our final estimated average downshift is approximately -305 cm(-1). The experimental value lies somewhere in the range between -290 and -420 cm(-1). Also, the absolute nu(OH) frequency is well reproduced in our calculations. "In-liquid" instantaneous correlation curves between nu(OH) and various typical H-bond strength parameters such as R(O...O), R(H...O), the intramolecular OH bond length, and the IR intensity are presented. Some of these correlations are robust and persist also for the rather distorted instantaneous geometries in the liquid; others are less so.     

Effects of ionic strength on SAXS data for proteins revealed by molecular dynamics simulations

The combination of small-angle X-ray solution scattering (SAXS) experiments and molecular dynamics (MD) simulations is now becoming a powerful tool to study protein conformations in solution at an atomic resolution. In this study, we investigated effects of ionic strength on SAXS data theoretically by using MD simulations of hen egg white lysozyme at various NaCl concentrations from 0 to 1 M. The calculated SAXS excess intensities showed a significant dependence on ion concentration, which originates from the different solvent density distributions in the presence and absence of ions. The addition of ions induced a slow convergence of the SAXS data, and a ～20 ns simulation is required to obtain convergence of the SAXS data with the presence of ions whereas only a 0.2 ns simulation is sufficient in the absence of ions. To circumvent the problem of the slow convergence in the presence of ions, we developed a novel method that reproduces the SAXS excess intensities with the presence of ions from short MD trajectories in pure water. By applying this method to SAXS data for the open and closed forms of transferrin at 1 M ion concentration, the correct form could be identified by simply using short MD simulations of the protein in pure water for 0.2 ns.     

Comparison of the hydration and diffusion of protons in perfluorosulfonic acid membranes with molecular dynamics simulations. scui@utk.edu

Classical molecular dynamics (MD) simulations were performed to determine the hydrated morphology and hydronium ion diffusion coefficients in two different perfluorosulfonic acid (PFSA) membranes as functions of water content. The structural and transport properties of 1143 equivalent weight (EW) Nafion, with its relatively long perfluoroether side chains, are compared to the short-side-chain (SSC) PFSA ionomer at an EW of 977. The separation of the side chains was kept uniform in both ionomers consisting of -(CF 2) 15- units in the backbone, and the degree of hydration was varied from 5 to 20 weight % water. The MD simulations indicated that the distribution of water clusters is more dispersed in the SSC ionomer, which leads to a more connected water-channel network at the low water contents. This suggests that the SSC ionomer may be more inclined to form sample-spanning aqueous domains through which transport of water and protons may occur. The diffusion coefficients for both hydronium ions and water molecules were calculated at hydration levels of 4.4, 6.4, 9.6, and 12.8 H 2O/SO 3H for each ionomer. When compared to experimental proton diffusion coefficients, this suggests that as the water content is increased the contribution of proton hopping to the overall proton diffusion increases.     

The structure,thermodynamic properties and infrared spectra of liquid water and ice

Monte Carlo simulations are used, together with models of the intramolecular and intermolecular potential surfaces, to model liquid water and several phases of ice. Intramolecular relaxation makes important contributions to both thermodynamic and structural properties. A quantum local mode analysis of the Monte Carlo configurations is used to predict the density of states and infrared absorption intensities for the intramolecular bending and stretching vibrations. The large shifts from the gas phase OH stretch frequencies observed experimentally in the liquid and solid phases are due to anharmonic terms in the intramolecular surface rather than to harmonic intermolecular coupling. A significant contribution to observed changes in IR intensity on condensation arises from the large molecular polarisability.     

Dependences of Water Permeation through Cyclic Octa-peptide Nanotubes on Channel Length and Membrane Thickness

Effects of the channel length and membrane thickness on the water permeation through the transmembrane cyclic octa-peptide nanotubes (octa-PNTs) have been studied by molecular dynamics (MD) simulations. The water osmotic permeability (p(f)) through the PNTs of k × (WL)(4)/POPE (1-palmitoyl-2-oleoyl-glycerophosphoethanolamine; k = 6, 7, 8, 9, and 10) was found to decay with the channel length (L) along the axis (～L(-2.0)). Energetic analysis showed that a series of water binding sites exist in these transmembrane PNTs, with the barriers of ～3k(B)T, which elucidates the tendency of p(f) well. Water diffusion permeability (p(d)) exhibits a relationship of ～L(-1.8), which results from the novel 1-2-1-2 structure of water chain in such confined nanolumens. In the range of simulation accuracy, the ratio (p(f)/p(d)) of the water osmotic and diffusion permeability is approximately a constant. MD simulations of water permeation through the transmembrane PNTs of 8 × (WL)(4)/octane with the different octane membrane thickness revealed that the water osmotic and diffusion permeability (p(f) and p(d)) are both independent of the octane membrane thickness, confirmed by the weak and nearly same interactions between the channel water and octane membranes with the different thickness. The results may be helpful for revealing the permeation mechanisms of biological water channels and designing artificial nanochannels.     

Comparison of the glass transition temperature and fragility parameter of isomalto-olygomer predicted by molecular dynamics simulations with those measured by differential scanning calorimetry

The purpose of this study is to examine whether molecular dynamics (MD) simulations using a commercially available software for personal computers can estimate the glass transition temperature (Tg) of amorphous systems containing pharmaceutically-relevant excipients. MD simulations were carried out with an amorphous matrix model constructed from isomaltoheptaose, and the Tg estimated from the calculated density versus temperature profile was compared with the Tg measured by differential scanning calorimetry (DSC) for freeze-dried isomalto-oligomer having an average molecular weight close to that of isomaltoheptaose. The Tg values determined by DSC were lower by 10 to 20 K than those extrapolated from the Tg values estimated by MD simulation. Fragility parameter was estimated to be 56 and 51 from MD simulation and from DSC measurement, respectively. Thus, the results suggest that MD simulation can provide approximate estimates for the Tg and fragility parameter of amorphous formulations. However, a reduction of the cooling rate, achievable by sufficiently elongating the simulation duration, is necessary for more accurate estimation.     

Human aquaporin 4 gating dynamics in dc and ac electric fields: a molecular dynamics study

Water self-diffusion within human aquaporin 4 has been studied using molecular dynamics (MD) simulations in the absence and presence of external ac and dc electric fields. The computed diffusive (p(d)) and osmotic (p(f)) permeabilities under zero-field conditions are (0.718 ± 0.24) × 10(-14) cm(3) s(-1) and (2.94 ± 0.47) × 10(-14) cm(3) s(-1), respectively; our p(f) agrees with the experimental value of (1.50 ± 0.6) × 10(-14) cm(3) s(-1). A gating mechanism has been proposed in which side-chain dynamics of residue H201, located in the selectivity filter, play an essential role. In addition, for nonequilibrium MD in external fields, it was found that water dipole orientation within the constriction region of the channel is affected by electric fields (e-fields) and that this governs the permeability. It was also found that the rate of side-chain flipping motion of residue H201 is increased in the presence of e-fields, which influences water conductivity further.     

Isotope effects in liquid water by infrared spectroscopy. IV. No free OH groups in liquid water

The presence of free OH (OH not H-bonded) in bulk water is a key element for the determination of its molecular structure. The OH covalent bond infrared (IR) absorption is highly sensitive to the molecular environment. For this reason, IR spectroscopy is used for the determination of free OH. A workable definition of this is obtained with methanol (MeOH) in hexane where minute quantities of free OH are present. These absorb at 3654?cm(-1) (a 27?cm(-1) redshift from the gas position) with a full width at half height of 35?cm(-1). The IR spectrum of water between room temperature and 95?°C does not display such a band near 3650?cm(-1). This indicates that we do not see, in the IR spectra, the "free" OH group. From this we conclude that it is not present in liquid water at least down to the 1000 ppm level which is the limit of detectivity of our spectrometer. Other spectroscopic considerations of methanol and water in acetonitrile solutions indicate that weak H-bonds are also not present in liquid water.     

Al3+, Ca2+, Mg2+, and Li+ in aqueous solution: calculated first-shell anharmonic OH vibrations at 300 K

The anharmonic OH stretching vibrational frequencies, ν(OH), for the first-shell water molecules around the Li(+), Ca(2+), Mg(2+), and Al(3+) ions in dilute aqueous solutions have been calculated based on classical molecular dynamics (MD) simulations and quantum-mechanical (QM) calculations. For Li(+)(aq), Ca(2+)(aq), Mg(2+)(aq), and Al(3+)(aq), our calculated IR frequency shifts, Δν(OH), with respect to the gas-phase water frequency, are about -300, -350, -450, and -750?cm(-1), compared to -290, -290, -420, and -830?cm(-1) from experimental infrared (IR) studies. The agreement is thus quite good, except for the order between Li(+) and Ca(2+). Given that the polarizing field from the Ca(2+) ion ought to be larger than that from Li(+)(aq), our calculated result seems reasonable. Also the absolute OH frequencies agree well with experiment. The method we used is a sequential four-step procedure: QM(electronic) to make a force field+MD simulation+QM(electronic) for point-charge-embedded M(n+) (H(2)O)(y) (second?shell) (H(2)O)(z) (third?shell) clusters+QM(vibrational) to yield the OH spectrum. The many-body Ca(2+)-water force-field presented in this paper is new. IR intensity-weighting of the density-of-states frequency distributions was carried out by means of the squared dipole moment derivatives.     

Structured water in polyelectrolyte dendrimers: understanding small angle neutron scattering results through atomistic simulation

Based on atomistic molecular dynamics (MD) simulations, the small angle neutron scattering (SANS) intensity behavior of a single generation-4 polyelectrolyte polyamidoamine starburst dendrimer is investigated at different levels of molecular protonation. The SANS form factor, P(Q), and Debye autocorrelation function, γ(r), are calculated from the equilibrium MD trajectory based on a mathematical approach proposed in this work. The consistency found in comparison against previously published experimental findings (W.-R. Chen, L. Porcar, Y. Liu, P. D. Butler, and L. J. Magid, Macromolecules 40, 5887 (2007)) leads to a link between the neutron scattering experiment and MD computation, and fresh perspectives. The simulations enable scattering calculations of not only the hydrocarbons but also the contribution from the scattering length density fluctuations caused by structured, confined water within the dendrimer. Based on our computational results, we explore the validity of using radius of gyration R(G) for microstructure characterization of a polyelectrolyte dendrimer from the scattering perspective.     

Odd-even variations in the wettability of n-alkanethiolate monolayers on gold by water and hexadecane: a molecular dynamics simulation study

Molecular dynamics (MD) simulations were performed to investigate odd-even chain length dependencies in the wetting properties of self-assembled monolayers (SAMs) of n-alkanethiols [CH3(CH2)n-1SH] on gold by water and hexadecane. Experimentally, the contact angle of hexadecane on the SAMs depends on whether n is odd or even, while contact angles for water show no odd-even dependence. Our MD simulations of this system included a microscopic droplet of either 256 water molecules or 60 hexadecane molecules localized on an n-alkanethiolate SAM on gold with either an even or odd chain length. Contact angles calculated for these nanoscopic droplets were consistent with experimentally observed macroscopic trends in wettability, namely, that hexadecane is sensitive to structural differences between odd- and even-chained SAMs while water is not. Structural properties for the SAMs (including features such as chain tilt, chain twist, and terminal methyl group tilt) were calculated during the MD simulations and used to generate IR spectra of these films that compared favorably with experimental spectra. MD simulations of SAMs in contact with slabs of water and hexadecane revealed that the effects of these solvents on the structure of the SAM was restricted to the chain terminus and had no effect on the inner structure of the SAM. The density profiles for water and hexadecane on the SAMs were different in that water displayed a significant depletion in its density at the liquid/SAM interface from its bulk value, while no such depletion occurred for hexadecane. This difference in contact may explain the lack of an odd-even variation in the wetting characteristics of water on these surfaces, because the water molecules are positioned further away from the surface and, therefore, are not sensitive to the structural differences in the average orientations for the terminal methyl groups in odd- and even-chained SAMs. In contrast, the differences in the wetting properties of hexadecane on the odd- and even-chained SAMs may reflect the closer proximity of these molecules to the SAM surface and a resulting greater sensitivity to the differences in the terminal methyl group orientations in the SAMs. SAM-solvent interaction energies were calculated during the MD simulations, yielding interaction energies that differed on the even- and odd-chained surfaces by approximately 10% for hexadecane and negligibly for water, in accord with estimates using experimental wetting results.     

Demonstrating accuracy of the already proposed protocol for structure elucidation of cyclodextrin inclusion complexes by validation using quantitative ROESY analysis

In this study, we attempt to ascertain the accuracy of the structures determined using our previously developed method and hence the accuracy of our method. In the present report, we have taken the guest molecule cetirizine (CTZ) and the host molecules are α-cyclodextrin (α-CD) and β-cyclodextrin (β-CD). Structures with good accuracy were elucidated using a productive fusion of experimental and computational methods. We performed molecular mechanics studies (MM) in light of experimental ROESY studies, followed by molecular dynamics studies (MD). The results from these studies were analyzed using quantitative ROESY analysis to determine the final accurate structures of the complexes. The accuracy of these structures was tested using density functional theory (DFT) that is an accurate method for structure determination. DFT studies were carried out using the functionals B3LYP and M06L with def-TZVP basis set and similarly quantitative ROESY analysis was performed for the obtained structures. The ROESY intensities of the structures obtained from MM and MD studies, were compared with ROESY intensities of the structures obtained from DFT studies. Calculated ROESY intensities of the structures obtained from B3LYP functional are comparable, with intensities of structures obtained from MM and MD studies, but M06L functional showed poor results. In addition to the accuracy of MM and MD studies, low computational cost and less time input make it good method for structural studies for CD inclusion complexes.

     

IR spectra of N-methylacetamide in water predicted by combined quantum mechanical/molecular mechanical molecular dynamics simulations

We applied the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulation method in assessing IR spectra of N-methylacetamide and its deuterated form in aqueous solutions. The model peptide is treated at the Austin Model 1 (AM1) level and the induced dipole effects by the solvent are incorporated in fluctuating solute dipole moments, which are calculated using partial charges from Mulliken population analyses without resorting to any available high-level ab initio dipole moment data. Fourier transform of the solute dipole autocorrelation function produces in silico IR spectra, in which the relative peak intensities and bandwidths of major amide bands are quantitatively compatible with experimental results only when both geometric and electronic polarizations of the peptide by the solvent are dealt with at the same quantum-mechanical level. We cast light on the importance of addressing dynamic charge fluctuations of the solute in calculating IR spectra by comparing classical and QM/MM MD simulation results. We propose the adjustable scaling factors for each amide mode to be directly compared with experimental data.