| 1. |
- Aidas, Kestutis, et al.
(författare)
-
A quantum mechanics/molecular dynamics study of electric field gradient fluctuations in the liquid phase. The case of Na+ in aqueous solution
- 2013
-
Ingår i: Physical Chemistry, Chemical Physics - PCCP. - : Royal Society of Chemistry (RSC). - 1463-9076 .- 1463-9084. ; 15:5, s. 1621-1631
-
Tidskriftsartikel (refereegranskat)abstract
- The Na-23 quadrupolar coupling constant of the Na+ ion in aqueous solution has been predicted using molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics methods for the calculation of electric field gradients. The developed computational approach is generally expected to provide reliable estimates of the quadrupolar coupling constants of monoatomic species in condensed phases, and we show here that intermolecular polarization and non-electrostatic interactions are of crucial importance as they result in a 100% increased quadrupolar coupling constant of the ion as compared to a simpler pure electrostatic picture. These findings question the reliability of the commonly applied classical Sternheimer approximation for the calculations of the electric field gradient. As it can be expected from symmetry considerations, the quadrupolar coupling constants of the 5- and 6-coordinated Na+ ions in solution are found to differ significantly.
|
|
| 2. |
- Aidas, Kestutis, et al.
(författare)
-
Gauge-origin independent magnetizabilities from hybrid quantum mechanics/molecular mechanics models: Theory and applications to liquid water
- 2007
-
Ingår i: Chemical Physics Letters. - : Elsevier BV. - 0009-2614. ; 442:4-6, s. 322-328
-
Tidskriftsartikel (refereegranskat)abstract
- The theory of a hybrid quantum mechanics/molecular mechanics (QM/MM) approach for gauge-origin independent calculations of the molecular magnetizability using Hartree-Fock or Density Functional Theory is presented. The method is applied to liquid water using configurations generated from classical Molecular Dynamics simulation to calculate the statistical averaged magnetizability. Based on a comparison with experimental data, treating only one water molecule quantum mechanically appears to be insufficient, while a quantum mechanical treatment of also the first solvation shell leads to good agreement between theory and experiment. This indicates that the gas-to-liquid phase shift for the molecular magnetizability is to a large extent of non-electrostatic nature. (c) 2007 Elsevier B.V. All rights reserved.
|
|
| 3. |
- Aidas, Kestutis, et al.
(författare)
-
On the performance of quantum chemical methods to predict solvatochromic effects: The case of acrolein in aqueous solution.
- 2008
-
Ingår i: Journal of Chemical Physics. - : AIP Publishing. - 0021-9606 .- 1089-7690. ; 128:19, s. 1-194503
-
Tidskriftsartikel (refereegranskat)abstract
- The performance of the Hartree-Fock method and the three density functionals B3LYP, PBE0, and CAM-B3LYP is compared to results based on the coupled cluster singles and doubles model in predictions of the solvatochromic effects on the vertical n-->pi(*) and pi-->pi(*) electronic excitation energies of acrolein. All electronic structure methods employed the same solvent model, which is based on the combined quantum mechanics/molecular mechanics approach together with a dynamical averaging scheme. In addition to the predicted solvatochromic effects, we have also performed spectroscopic UV measurements of acrolein in vapor phase and aqueous solution. The gas-to-aqueous solution shift of the n-->pi(*) excitation energy is well reproduced by using all density functional methods considered. However, the B3LYP and PBE0 functionals completely fail to describe the pi-->pi(*) electronic transition in solution, whereas the recent CAM-B3LYP functional performs well also in this case. The pi-->pi(*) excitation energy of acrolein in water solution is found to be very dependent on intermolecular induction and nonelectrostatic interactions. The computed excitation energies of acrolein in vacuum and solution compare well to experimental data.
|
|
| 4. |
- Aidas, Kestutis, et al.
(författare)
-
Photoabsorption of Acridine Yellow and Proflavin Bound to Human Serum Albumin Studied by Means of Quantum Mechanics/Molecular Dynamics
- 2013
-
Ingår i: Journal of Physical Chemistry B. - : American Chemical Society (ACS). - 1520-6106 .- 1520-5207. ; 117:7, s. 2069-2080
-
Tidskriftsartikel (refereegranskat)abstract
- Attempting to unravel mechanisms in optical probing of proteins, we have performed pilot calculations of two cationic chromophores-acridine yellow and proflavin-located at different binding sites within human serum albumin, including the two primary drug binding sites as well as a heme binding site. The computational scheme adopted involves classical molecular dynamics simulations of the ligands bound to the protein and subsequent linear response polarizable embedding density functional theory calculations of the excitation energies. A polarizable embedding potential consisting of point charges fitted to reproduce the electrostatic potential and isotropic atomic polarizabilities computed individually for every residue of the protein was used in the linear response calculations. Comparing the calculated aqueous solution-to-protein shifts of maximum absorption energies to available experimental data, we concluded that the cationic proflavin chromophore is likely not to bind albumin at its drug binding site I nor at its heme binding site. Although agreement with experimental data could only be obtained in qualitative terms, our results clearly indicate that the difference in optical response of the two probes is due to deprotonation, and not, as earlier suggested, to different binding sites. The ramifications of this finding for design of molecular probes targeting albumin or other proteins is briefly discussed.
|
|
| 5. |
- Aidas, Kestutis, et al.
(författare)
-
Solvent effects on NMR isotropic shielding constants. A comparison between explicit polarizable discrete and continuum approaches
- 2007
-
Ingår i: The Journal of Physical Chemistry Part A: Molecules, Spectroscopy, Kinetics, Environment and General Theory. - : American Chemical Society (ACS). - 1520-5215. ; 111:20, s. 4199-4210
-
Tidskriftsartikel (refereegranskat)abstract
- The gas-to-aqueous solution shifts of the O-17 and C-13 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.
|
|
| 6. |
- Aidas, Kestutis, et al.
(författare)
-
The Dalton quantum chemistry program system
- 2014
-
Ingår i: Wiley Interdisciplinary Reviews. Computational Molecular Science. - : Wiley. - 1759-0876. ; 4:3, s. 269-284
-
Tidskriftsartikel (refereegranskat)abstract
- Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, MOller-Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from for a number of UNIX platforms.
|
|
| 7. |
- Eriksen, Janus J., et al.
(författare)
-
Computational Protocols for Prediction of Solute NMR Relative Chemical Shifts. A Case Study of L-Tryptophan in Aqueous Solution
- 2011
-
Ingår i: Journal of Computational Chemistry. - : Wiley. - 0192-8651 .- 1096-987X. ; 32:13, s. 2853-2864
-
Tidskriftsartikel (refereegranskat)abstract
- In this study, we have applied two different spanning protocols for obtaining the molecular conformations of L-tryptophan in aqueous solution, namely a molecular dynamics simulation and a molecular mechanics conformational search with subsequent geometry re-optimization of the stable conformers using a quantum mechanically based method. These spanning protocols represent standard ways of obtaining a set of conformations on which NMR calculations may be performed. The results stemming from the solute-solvent configurations extracted from the MD simulation at 300 K are found to be inferior to the results stemming from the conformations extracted from the MM conformational search in terms of replicating an experimental reference as well as in achieving the correct sequence of the NMR relative chemical shifts of L-tryptophan in aqueous solution. We find this to be due to missing conformations visited during the molecular dynamics run as well as inaccuracies in geometrical parameters generated from the classical molecular dynamics simulations.
|
|
| 8. |
- Frecus, Bogdan, et al.
(författare)
-
EPR spin Hamiltonian parameters of encapsulated spin-labels : impact of the hydrogen bonding topology
- 2013
-
Ingår i: Physical Chemistry, Chemical Physics - PCCP. - : Royal Society of Chemistry (RSC). - 1463-9076 .- 1463-9084. ; 15:7, s. 2427-2434
-
Tidskriftsartikel (refereegranskat)abstract
- Encapsulation of spin-labels into "host'' compounds, like cucurbit[n]urils or cyclodextrins, in solutions has profound effects on the EPR spin Hamiltonian parameters of the spin-labels. In this work we study the microscopic origin of the EPR spin Hamiltonian parameters of spin-labels enclosed in hydrophobic cavities. We focus on the dependence of the EPR properties of encapsulated spin-labels on the hydrogen bonding topologies that occur upon encapsulation, and quantize various contributions to these parameters according to specific hydrogen bonding patterns. The obtained results provide refined insight into the role of the hydrogen bonding induced encapsulation shifts of EPR spin Hamiltonian parameters in solvated "spin-label@host compound'' complexes.
|
|
| 9. |
- Genheden, Samuel, et al.
(författare)
-
Accurate Predictions of Nonpolar Solvation Free Energies Require Explicit Consideration of Binding-Site Hydration
- 2011
-
Ingår i: Journal of the American Chemical Society. - : American Chemical Society (ACS). - 1520-5126 .- 0002-7863. ; 133:33, s. 13081-13092
-
Tidskriftsartikel (refereegranskat)abstract
- Continuum solvation methods are frequently used to increase the efficiency of computational methods to estimate free energies. In this paper, we have evaluated how well such methods estimate the nonpolar solvation free-energy change when a ligand binds to a protein. Three different continuum methods at various levels of approximation were considered, viz., the polarized continuum model (PCM), a method based on cavity and dispersion terms (CD), and a method based on a linear relation to the solvent-accessible surface area (SASA). Formally rigorous double-decoupling thermodynamic integration was used as a benchmark for the continuum methods. We have studied four protein-ligand complexes with binding sites of varying solvent exposure, namely the binding of phenol to ferritin, a biotin analogue to avidin, 2-aminobenzimidazole to trypsin, and a substituted galactoside to galectin-3. For ferritin and avidin, which have relatively hidden binding sites, rather accurate nonpolar solvation free energies could be obtained with the continuum methods if the binding site is prohibited to be filled by continuum water in the unbound state, even though the simulations and experiments show that the ligand replaces several water molecules upon binding. For the more solvent exposed binding sites of trypsin and galectin-3, no accurate continuum estimates could be obtained, even if the binding site was allowed or prohibited to be filled by continuum water. This shows that continuum methods fail to give accurate free energies on a wide range of systems with varying solvent exposure because they lack a microscopic picture of binding-site hydration as well as information about the entropy of water molecules that are in the binding site before the ligand binds. Consequently, binding affinity estimates based upon continuum solvation methods will give absolute binding energies that may differ by up to 200 kJ/mol depending on the method used. Moreover, even relative energies between ligands with the same scaffold may differ by up to 75 kJ/mol. We have tried to improve the continuum solvation methods by adding information about the solvent exposure of the binding site or the hydration of the binding site, and the results are promising at least for this small set of complexes.
|
|
| 10. |
- Genheden, Samuel, et al.
(författare)
-
Nonpolar Solvation Free Energies of Protein-Ligand Complexes
- 2010
-
Ingår i: Journal of Chemical Theory and Computation. - : American Chemical Society (ACS). - 1549-9618 .- 1549-9626. ; 6:11, s. 3558-3568
-
Tidskriftsartikel (refereegranskat)abstract
- Recent investigations have indicated that different solvation methods give qualitatively different results for the nonpolar solvation contribution to ligand-binding affinities. Therefore, we have calculated the nonpolar solvation contribution to the free energy of benzene binding to the T4 lysozyme Leu99Ala mutant using thermodynamic integration (TI) and three approximate methods. The total binding free energy was calculated with TI and then decomposed into contributions from the solvent and the solute. The nonpolar contribution from the solute was compared to approximate methods within the framework of the molecular-mechanics and generalized Born with surface area method (MM/GBSA). First, the nonpolar solvation energy was calculated with a linear relation to the solvent-accessible surface area (SASA). Second, a recent approach that divides the nonpolar solvation energy into cavity and dispersion parts was used, and third, the nonpolar solvation energy was calculated with the polarized continuum model (PCM). Surprisingly, the simple SASA estimate reproduces the TI results best. However, the reason for this is that all continuum methods assume that the benzene cavity is filled with water for the free protein, contrary to both experimental and simulation results. We present a method to avoid this assumption and then, PCM provides results that are closest to the results obtained with TI.
|
|
