Probing the transition from hydrophilic to hydrophobic solvation with atomic scale resolution

Picosecond and femtosecond X-ray absorption spectroscopy is used to probe the changes of the solvent shell structure upon electron abstraction of aqueous iodide using an ultrashort laser pulse. The transient L(1,3) edge EXAFS at 50 ps time delay points to the formation of an expanded water cavity around the iodine atom, in good agreement with classical and quantum mechanical/molecular mechanics (QM/MM) molecular dynamics (MD) simulations. These also show that while the hydrogen atoms pointed toward iodide, they predominantly point toward the bulk solvent in the case of iodine, suggesting a hydrophobic behavior. This is further confirmed by quantum chemical (QC) calculations of I(-)/I(0)(H(2)O)(n=1-4) clusters. The L(1) edge sub-picosecond spectra point to the existence of a transient species that is not present at 50 ps. The QC calculations and the QM/MM MD simulations identify this transient species as an I(0)(OH(2)) complex inside the cavity. The simulations show that upon electron abstraction most of the water molecules move away from iodine, while one comes closer to form the complex that lives for 3-4 ps. This time is governed by the reorganization of the main solvation shell, basically the time it takes for the water molecules to reform an H-bond network. Only then is the interaction with the solvation shell strong enough to pull the water molecule of the complex toward the bulk solvent. Overall, much of the behavior at early times is determined by the reorientational dynamics of water molecules and the formation of a complete network of hydrogen bonded molecules in the first solvation shell.     

Simulations of the temperature dependence of amide I vibration

For spectroscopic studies of peptide and protein thermal denaturation it is important to single out the contribution of the solvent to the spectral changes from those originated in the molecular structure. To obtain insights into the origin and size of the temperature solvent effects on the amide I spectra, combined molecular dynamics and density functional simulations were performed with the model N-methylacetamide molecule (NMA). The computations well reproduced frequency and intensity changes previously observed in aqueous NMA solutions. An empirical correction of vacuum frequencies in single NMA molecule based on the electrostatic potential of the water molecules provided superior results to a direct density functional average obtained for a limited number of solute-solvent clusters. The results thus confirm that the all-atom quantum and molecular mechanics approach captures the overall influence of the temperature dependent solvent properties on the amide I spectra and can improve the accuracy and reliability of molecular structural studies.     

On the solvation of the Zn2+ ion in methanol: a combined quantum mechanics, molecular dynamics, and EXAFS approach

The solvation properties of the Zn(2+) ion in methanol solution have been investigated using a combined approach based on molecular dynamics (MD) simulations and extended X-ray absorption fine structure (EXAFS) experimental results. The quantum mechanical potential energy surface for the interaction of the Zn(2+) ion with a methanol molecule has been calculated taking into account the effect of bulk solvent by the polarizable continuum model (PCM). The effective Zn-methanol interactions have been fitted by suitable analytical potentials, and have been utilized in the MD simulation to obtain the structural properties of the solution. The reliability of the whole procedure has been assessed by comparing the theoretical structural results with the EXAFS experimental data. The structural parameters of the first solvation shells issuing from the MD simulations provide an effective complement to the EXAFS experiments.     

Development of the force field parameters for phosphoimidazole and phosphohistidine

Phosphorylation of histidine-containing proteins is a key step in the mechanism of many phosphate transfer enzymes (kinases, phosphatases) and is the first stage in a wide variety of signal transduction cascades in bacteria, yeast, higher plants, and mammals. Studies of structural and dynamical aspects of such enzymes in the phosphorylated intermediate states are important for understanding the intimate molecular mechanisms of their functioning. Such information may be obtained via molecular dynamics and/or docking simulations, but in this case appropriate force field parameters for phosphohistidine should be explicitly defined. In the present article we describe development of the GROMOS96 force field parameters for phosphoimidazole molecule--a realistic model of the phosphohistidine side chain. The parameterization is based on the results of ab initio quantum chemical calculations with subsequent refinement and testing using molecular mechanics and molecular dynamics simulations. The set of force constants and equilibrium geometry is employed to derive force field for the phosphohistidine moiety. Resulting parameters and topology are incorporated into the molecular modeling package GROMACS and used in molecular dynamics simulations of a phosphohistidine-containing protein in explicit solvent.     

力场/溶剂水模型对EGFRvⅢ抗原-MR1(scFv)抗体复合物的MM-PBSA计算精度的影响与实验验证

以表皮生长因子Ⅲ型突变体(EGFRvⅢ)抗原多肽与其抗体(MR1)及其人源化突变体的复合物结构为出发点,采用分子动力学中的6种常用力场及3种常用溶剂水模型,分别对上述抗原-抗体复合物进行100ns的分子动力学模拟与分子力学和连续介质模型计算自由能(MM-PBSA),并在实验上利用等温滴定量热(ITC)仪测定了抗原和抗体相互作用的热力学参数.通过在结构变化、能量变化及野生型与突变体比较等几个方面进行综合分析,给出了最佳的计算模型.对不同力场及水模型计算精度等相关问题进行了探讨.     

Solvent effects on peroxynitrite structure and properties from QM/MM simulations

We have investigated the structure and the vibrational spectrum of peroxynitrite anion in aqueous solution by means of combined quantum-classical (QM/MM) molecular dynamics simulations. In our QM/MM scheme, the reactant was modeled using density functional theory with a Gaussian basis set and the solvent was described using the mean-field TIP4P and the polarizable TIP4P-FQ force fields. The choice of basis sets, functionals and force field parameters has been validated by performing calculations on isolated peroxynitrite and on small peroxynitrite-water complexes. Poor values for isolated peroxynitrite structural properties and vibrational frequencies are found for most ab initio methods, particularly regarding the ON-OO(-) bond distance and the harmonic stretching frequency, probably due to the singlet-triplet instability found in the HF wave function. On the other hand, DFT methods yield results in better agreement with high level CCSD(T) ab initio calculations. We have computed the vibrational spectrum for aqueous peroxynitrite by calculating the Fourier transform of the velocity autocorrelation function extracted from the QM-MM molecular dynamics simulations. Our computational scheme, which allows for the inclusion of both anharmonicity and solvent effects, is able to clarify previous discrepancies regarding the experimental spectra assignments and to shed light on the subtle interplay between solvation and peroxynitrite structure and properties.     

Implementation of the analytic energy gradient for the combined time-dependent density functional theory/effective fragment potential method: application to excited-state molecular dynamics simulations

Excited-state quantum mechanics/molecular mechanics molecular dynamics simulations are performed, to examine the solvent effects on the fluorescence spectra of aqueous formaldehyde. For that purpose, the analytical energy gradient has been derived and implemented for the linear-response time-dependent density functional theory (TDDFT) combined with the effective fragment potential (EFP) method. The EFP method is an efficient ab initio based polarizable model that describes the explicit solvent effects on electronic excitations, in the present work within a hybrid TDDFT/EFP scheme. The new method is applied to the excited-state MD of aqueous formaldehyde in the n-π* state. The calculated π*→n transition energy and solvatochromic shift are in good agreement with other theoretical results.     

On the interpretation of continuous wave electron spin resonance spectra of tempo-palmitate in 5-cyanobiphenyl

Electron spin resonance (ESR) measurements are highly informative on the dynamic behavior of molecules, which is of fundamental importance to understand their stability, biological functions and activities, and catalytic action. The wealth of dynamic information which can be extracted from a continuous wave electron spin resonance (cw-ESR) spectrum can be inferred by a basic theoretical approach defined within the stochastic Liouville equation formalism, i.e., the direct inclusion of motional dynamics in the form of stochastic (Fokker-Planck/diffusive) operators in the super Hamiltonian H governing the time evolution of the system. Modeling requires the characterization of magnetic parameters (e.g., hyperfine and Zeeman tensors) and the calculation of ESR observables in terms of spectral densities. The magnetic observables can be pursued by the employment of density functional theory which is apt, provided that hybrid functionals are employed, for the accurate computation of structural properties of molecular systems. Recently, an ab initio integrated computational approach to the in silico interpretation of cw-ESR spectra of multilabeled systems in isotropic fluids has been discussed. In this work we present the extension to the case of nematic liquid crystalline environments by performing simulations of the ESR spectra of the prototypical nitroxide probe 4-(hexadecanoyloxy)-2,2,6,6-tetramethylpiperidine-1-oxy in isotropic and nematic phases of 5-cyanobiphenyl. We first discuss the basic ingredients of the integrated approach, i.e., (1) determination of geometric and local magnetic parameters by quantum-mechanical calculations, taking into account the solvent and, when needed, the vibrational averaging contributions; (2) numerical solution of a stochastic Liouville equation in the presence of diffusive rotational dynamics, based on (3) parameterization of diffusion rotational tensor provided by a hydrodynamic model. Next we present simulated spectra with minimal resorting to fitting procedures, proving that the combination of sensitive ESR spectroscopy and sophisticated modeling can be highly helpful in providing three-dimensional structural and dynamic information on molecular systems in anisotropic environments.     

Protonation effects on the resonance Raman properties of a novel (terpyridine)Ru(4H-imidazole) complex: an experimental and theoretical case study

The optically active states in a novel (terpyridine)Ru(4H-imidazole) complex displaying an unusually broad and red-shifted absorption in the visible range are investigated experimentally and theoretically. Since this property renders the complex promising for an application as sensitizer in dye-sensitized solar cells, a detailed knowledge on the correlation between features in the absorption spectrum and structural elements is indispensable in order to develop strategies for spectroscopy/theory-guided design of such molecular components. To this aim, time-dependent density functional theory calculations, including solvent effects, are employed to analyze the experimental UV-vis absorption and resonance Raman (RR) spectra of the unprotonated and protonated forms of the complex. This provides a detailed photophysical picture for a complex belonging to a novel class of Ru-polypyridine black absorbers, which can be tuned by external pH stimuli. The complex presents two absorption maxima in the visible region, which are assigned by the calculations to metal-to-ligand charge transfer (MLCT) and intra-ligand states, respectively. RR simulations are performed in resonance with both bands and are found to correctly reproduce the observed effects of protonation. Finally, the examination of the molecular orbitals and of the RR spectra for the MLCT state shows that protonation favors a charge transfer excitation to the 4H-imidazole ligand.     

Development and validation of an integrated computational approach for the study of ionic species in solution by means of effective two-body potentials. The case of Zn2+, Ni2+, and Co2+ in aqueous solutions

In this paper we have developed an effective computational procedure for the structural and dynamical investigation of ions in aqueous solutions. Quantum mechanical potential energy surfaces for the interaction of a transition metal ion with a water molecule have been calculated taking into account the effect of bulk solvent by the polarizable continuum model (PCM). The effective ion-water interactions have been fitted by suitable analytical potentials, and have been utilized in molecular dynamics (MD) simulations to obtain structural and dynamical properties of the ionic aqueous solutions. This procedure has been successfully applied to the Co2+-H2O open-shell system and, for the first time, Co-oxygen and Co-hydrogen pair potential functions have been determined and employed in MD simulations. The reliability of the whole procedure has been assessed by applying it also to the Zn2+ and Ni2+ aqueous solutions, and the structural and dynamical properties of the three systems have been calculated by means of MD simulations and have been found to be in very good agreement with experimental results. The structural parameters of the first solvation shells issuing from the MD simulations provide an effective complement to extended X-ray absorption fine structure (EXAFS) experiments.     

Solvent effects on NMR isotropic shielding constants. a comparison between explicit polarizable discrete and continuum approaches

The gas-to-aqueous solution shifts of the 17O and 13C NMR isotropic shielding constants for the carbonyl chromophore in formaldehyde and acetone are investigated. For the condensed-phase problem, we use the hybrid density functional theory/molecular mechanics approach in combination with a statistical averaging over an appropriate number of solute-solvent configurations extracted from classical molecular dynamics simulations. The PBE0 exchange-correlation functional and the 6-311++G(2d,2p) basis set are used for the calculation of the shielding constants. London atomic orbitals are employed to ensure gauge-origin independent results. The effects of the bulk solvent molecules are found to be crucial in order to calculate accurate solvation shifts of the shielding constants. Very good agreement between the computed and experimental solvation shifts is obtained for the shielding constants of acetone when a polarizable water potential is used. Supermolecular results based on geometry-optimized molecular structures are presented. We also compare the results obtained with the polarizable continuum model to the results obtained using explicit MM molecules to model the bulk solvent effect.     

Ferric complexes of 3-hydroxy-4-pyridinones characterized by density functional theory and Raman and UV-vis spectroscopies

Deferiprone and other 3-hydroxy-4-pyridinones are used in metal chelation therapy of iron overload. To investigate the structure and stability of these compounds in the natural aqueous environment, ferric complexes of deferiprone and amino acid maltol conjugates were synthesized and studied by computational and optical spectroscopic methods. The complexation caused characteristic intensity changes, a 300× overall enhancement of the Raman spectrum, and minor changes in UV-vis absorption. The spectra were interpreted on the basis of density functional theory (DFT) calculations. The CAM-B3LYP and ωB97XD functionals with CPCM solvent model were found to be the most suitable for simulations of the UV-vis spectra, whereas B3LYP, B3LYPD, B3PW91, M05-2X, M06, LC-BLYP, ωB97XD, and CAM-B3LYP functionals were all useful for simulation of the Raman scattering. Characteristic Raman band frequencies for 3-hydroxy-4-pyridinones were assigned to molecular vibrations. The computed conformer energies consistently suggest the presence of another isomer of the deferiprone-ferric complex in solution, in addition to that found previously by X-ray crystallography. However, the UV-vis and Raman spectra of the two species are similar and could not be resolved. In comparison to UV-vis, the Raman spectra and their combination with calculations appear more promising for future studies of iron sequestrating drugs and artificial metalloproteins as they are more sensitive to structural details.     

Solution structure of papain as studied by molecular mechanics and molecular dynamics techniques

In this paper, we report the molecular mechanics and molecular dynamics studies of the hydration of papain using the AMBER and CHARMm programs. We studied papain in an environment with minimal hydration involving only the X-ray waters and also in the hydrated environment by adding an extra 9 Å layer of water around the residues. The effect of nonbond cutoff was studied by performing minimizations with 8 Å and 15 Å nonbond cutoffs using the program AMBER. Two different solvent models—a constant dielectric and a distance-dependent dielectric—were considered. The AMBER-minimized structure and the average structure obtained from the CHARMm simulations for papain solvated with X-ray waters are presented and compared with the X-ray crystal structure results. Results of a similar comparison of the hydrated structures were also presented. The calculated root mean square deviation between the minimized structure and the X-ray structure is smaller for the hydrated system than for the system hydrated with only the X-ray waters. Results of the molecular mechanics and molecular dynamics simulations were presented for the various regions of papain. The hydration of the active site of papain and the effect of hydration on the torsional motion of the active site residues are presented. © 1996 by John Wiley & Sons, Inc.     

New AMBER force field parameters of heme iron for cytochrome P450s determined by quantum chemical calculations of simplified models

The heme protein, cytochrome P450, is an oxidoreductase that plays an important role in drug metabolism. To model P450s using molecular mechanics methods and classical molecular dynamics simulations, force field parameters and atomic charges are required. Because these parameters are generally obtained by quantum chemical methods, an appropriate simplified model for the iron-porphyrin system was needed. In this study, two models with a five-coordinated Fe(III) mimicking the sextet spin state of P450s are proposed, which are optimized by semiempirical and ab initio unrestricted Hartree-Fock methods. The results produced using the simpler of the two models were similar to those of the more complex model; therefore, the more simplified model of P450 can be used without a loss of accuracy. Furthermore, several quantum chemical calculations were carried out on the simpler model to investigate which method was most suitable for iron-porphyrin systems. The results calculated by hybrid density functional theory (DFT), with the MIDI basis set for iron, reproduced the three-dimensional structures determined by X-ray diffraction and extended X-ray absorption fine-structure experiments. From these results, atomic charges and force-field parameters for molecular mechanics and molecular dynamics calculations were obtained.     

O?+?C2H4 potential energy surface: lowest-lying singlet at the multireference level

Extensive density functional theory (DFT) calculations have been performed to develop a force field for the classical molecular dynamics (MD) simulations of various azobenzene derivatives. Besides azobenzene, we focused on a thiolated azobenzene’s molecular rod (4′-{[(1,1′-biphenyl)-4-yl]diazenyl}-(1,1′-biphenyl)-4-thiol) that has been previously demonstrated to photoisomerize from trans to cis with high yields on surfaces. The developed force field is an extension of OPLS All Atoms, and key bonding parameters are parameterized to reproduce the potential energy profiles calculated by DFT. For each of the parameterized molecule, we propose three sets of parameters: one best suited for the trans configuration, one for the cis configuration, and finally, a set able to describe both at a satisfactory degree. The quality of the derived parameters is evaluated by comparing with structural and vibrational experimental data. The developed force field opens the way to the classical MD simulations of self-assembled monolayers (SAMs) of azobenzene’s molecular rods, as well as to the quantum mechanics/molecular mechanics study of photoisomerization in SAMs.     

XAFS spectral analysis of the cadmium coordination geometry in cadmium thiolate clusters in metallothionein

We report the combination of measurement and prediction of X-ray absorption fine structure (XAFS) data, where the term XAFS refers to the overall spectrum that encompasses both the X-ray Absorption Near Edge Structure (XANES) region as well as the Extended X-ray Absorption Fine Structure (EXAFS) region, to evaluate the cadmium thiolate cluster structures in the metalloprotein metallothionein. XAFS spectra were simulated using coordinates from molecular models of the protein calculated by molecular mechanics/molecular dynamics (MM3/MD), from NMR analyses, and from analysis of X-ray diffraction data. XAFS spectra were also simulated using the coordinates from X-ray crystallographic data for [Cd(SPh)4]2-, CdS, [Cd2(mu-SPh)2(SPh)4]2-, and [Cd4(mu-SPh)6(SPh)4]2-. The simulated XAFS data that were calculated using the FEFF8 program closely resemble the experimental data reported for [Cd(SPh)4]2-, CdS, [Cd2(mu-SPh)2(SPh)4]2-, [Cd4(mu-SPh)6(SPh)4]2-, rabbit liver metallothionein cadmium alpha-domain (Cd4-alpha MT), and cadmium rabbit liver betaalpha metallothionein (Cd7-betaalpha MT). MM3 force field parameters were modified to include cadmium-sulfur bonding and were initially set to values derived from published X-ray diffraction and EXAFS experimental data. The force field was further calibrated and adjusted through comparison between experimental spectra taken from the literature and simulated XAFS spectra calculated using the FEFF8 program in combination with atomic coordinates from MM3/MD energy minimization. MM3/MD techniques were used with the calibrated force field to predict the high-resolution structure of the metal clusters in rabbit liver Cd7-MT. Structures for Cd3S9 (beta) MT and Cd4S11 (alpha) MT domains from MM3/MD calculations and those previously reported for Cd7-MT on the basis of 1H and 113Cd NMR data were compared. Structural differences between the different models for these cadmium thiolate clusters were evident. Combining the measurement and simulation of XAFS data provides an excellent method of assessing, modeling, and predicting metal-binding sites in metalloproteins when X-ray absorption spectroscopy (XAS) data are available.     

A theoretical study of the mechanism of the nucleotidyl transfer reaction catalyzed by yeast RNA polymerase II

The mechanism of the nucleotidyl transfer reaction catalyzed by yeast RNA polymerase II has been investigated using molecular mechanics and quantum mechanics methods.Molecular dynamics(MD) simulations were carried out using the TIP3 water model and generalized solvent boundary potential(GSBP) by CHARMM based on the X-ray crystal structure.Two models of the ternary elongation complex were constructed based on CHARMM MD calculations.All the species including reactants,transition states,intermediates,and products were optimized using the DFT-PBE method coupled with the basis set DZVP and the auxiliary basis set GEN-A2.Three pathways were explored using the DFT method.The most favorable reaction pathway involves indirect proton migration from the RNA primer 3’-OH to the oxygen atom of-phosphate via a solvent water molecule,proton rotation from the oxygen atom of-phosphate to the-phosphate side,the RNA primer 3’-O nucleophilic attack on the-phosphorus atom,and P-O bond breakage.The corresponding reaction potential profile was obtained.The rate limiting step,with a barrier height of 21.5 kcal/mol,is the RNA primer 3’-O nucleophilic attack,rather than the commonly considered proton transfer process.A high-resolution crystal structure including crystallographic water molecules is required for further studies.     

99Tc NMR as a promising technique for structural investigation of biomolecules: theoretical studies on the solvent and thermal effects of phenylbenzothiazole complex

The phenylbenzothiazole compounds show antitumor properties and are highly selective. In this paper, the ⁹⁹Tc chemical shifts based on the (^99mTc)(CO)₃(NNO) complex conjugated to the antitumor agent 2‐(4′‐aminophenyl)benzothiazole are reported. Thermal and solvent effects were studied computationally by quantum‐chemical methods, using the density functional theory (DFT) (DFT level BPW91/aug‐cc‐pVTZ for the Tc and BPW91/IGLO‐II for the other atoms) to compute the NMR parameters for the complex. We have calculated the ⁹⁹Tc NMR chemical shifts of the complex in gas phase and solution using different solvation models (polarizable continuum model and explicit solvation). To evaluate the thermal effect, molecular dynamics simulations were carried, using the atom‐centered density matrix propagation method at the DFT level (BP86/LanL2dz). The results highlight that the ⁹⁹Tc NMR spectroscopy can be a promising technique for structural investigation of biomolecules, at the molecular level, in different environments. Copyright © 2014 John Wiley & Sons, Ltd.