共查询到20条相似文献,搜索用时 15 毫秒
1.
In this article, we propose a new multigrid-based algorithm for solving integral equations of the reference interactions site model (RISM). We also investigate the relationship between the parameters of the algorithm and the numerical accuracy of the hydration free energy calculations by RISM. For this purpose, we analyzed the performance of the method for several numerical tests with polar and nonpolar compounds. The results of this analysis provide some guidelines for choosing an optimal set of parameters to minimize computational expenses. We compared the performance of the proposed multigrid-based method with the one-grid Picard iteration and nested Picard iteration methods. We show that the proposed method is over 30 times faster than the one-grid iteration method, and in the high accuracy regime, it is almost seven times faster than the nested Picard iteration method. 相似文献
2.
Masahiro Kinoshita Yuko Okamoto Fumio Hirata 《Journal of computational chemistry》1997,18(10):1320-1326
We have developed an algorithm for solving the reference interaction site model (RISM) equations for water near a solute molecule with many atomic sites (interaction sites). It is a hybrid of the Newton–Raphson and Picard methods and is judiciously constructed. Various considerations are given so that the computer time can be saved as much as possible. The robustness and high efficiency of the algorithm has been demonstrated for calculating hydration free energies of Met-enkephalin (a peptide with 75 sites) with different conformations. The Jacobian matrix is treated as part of the input data, and it has been found that the same matrix can be used for a considerably large set of different conformations of the solute molecule. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1320–1326, 1997 相似文献
3.
Prof. Ran Friedman 《Chemphyschem》2023,24(2):e202200516
The use of actinides for medical, scientific and technological purposes has gained momentum in the recent years. This creates a need to understand their interactions with biomolecules, both at the interface and as they become complexed. Calculation of the Gibbs binding energies of the ions to biomolecules, i. e., the Gibbs energy change associated with a transfer of an ion from the water phase to its binding site, could help to understand the actinides’ toxicities and to design agents that bind them with high affinities. To this end, there is a need to obtain accurate reference values for actinide hydration, that for most actinides are not available from experiment. In this study, a set of ionic radii is developed that enables future calculations of binding energies for Pu3+ and five actinides with renewed scientific and technological interest: Ac3+, Am3+, Cm3+, Bk3+ and Cf3+. Reference hydration energies were calculated using quantum chemistry and ion solvation theory and agree well for all ions except Ac3+, where ion solvation theory seems to underestimate the magnitude of the Gibbs hydration energy. The set of radii and reference energies that are presented here provide means to calculate binding energies for actinides and biomolecules. 相似文献
4.
Electrostatic free energies of solvation for 15 neutral amino acid side chain analogs are computed. We compare three methods of varying computational complexity and accuracy for three force fields: free energy simulations, Poisson-Boltzmann (PB), and linear response approximation (LRA) using AMBER, CHARMM, and OPLS-AA force fields. We find that deviations from simulation start at low charges for solutes. The approximate PB and LRA produce an overestimation of electrostatic solvation free energies for most of molecules studied here. These deviations are remarkably systematic. The variations among force fields are almost as large as the variations found among methods. Our study confirms that success of the approximate methods for electrostatic solvation free energies comes from their ability to evaluate free energy differences accurately. 相似文献
5.
Hydration free energy (HFE) is generally used for evaluating molecular solubility, which is an important property for pharmaceutical and chemical engineering processes. Accurately predicting HFE is also recognized as one fundamental capability of molecular mechanics force field. Here, we present a systematic investigation on HFE calculations with AMOEBA polarizable force field at various parameterization and simulation conditions. The HFEs of seven small organic molecules have been obtained alchemically using the Bennett Acceptance Ratio method. We have compared two approaches to derive the atomic multipoles from quantum mechanical calculations: one directly from the new distributed multipole analysis and the other involving fitting to the electrostatic potential around the molecules. Wave functions solved at the MP2 level with four basis sets (6-311G*, 6-311++G(2d,2p), cc-pVTZ, and aug-cc-pVTZ) are used to derive the atomic multipoles. HFEs from all four basis sets show a reasonable agreement with experimental data (root mean square error 0.63 kcal/mol for aug-cc-pVTZ). We conclude that aug-cc-pVTZ gives the best performance when used with AMOEBA, and 6-311++G(2d,2p) is comparable but more efficient for larger systems. The results suggest that the inclusion of diffuse basis functions is important for capturing intermolecular interactions. The effect of long-range correction to van der Waals interaction on the hydration free energies is about 0.1 kcal/mol when the cutoff is 12?, and increases linearly with the number of atoms in the solute/ligand. In addition, we also discussed the results from a hybrid approach that combines polarizable solute with fixed-charge water in the HFE calculation. 相似文献
6.
Cristiano Ruch Werneck Guimares Ricardo Bicca de Alencastro 《International journal of quantum chemistry》2001,85(6):713-726
Thrombin is a serine protease involved in blood coagulation. Since thrombin inhibitors appear to be effective in the treatment and prevention of thrombotic and embolic disorders, considerable attention has been focused on the structure and interactions of the enzyme. In this work, we calculated the relative free energies of hydration of the new thrombin inhibitor candidates, p‐substituted derivatives of benzamidine, a well‐known noncovalent thrombin inhibitor. We used molecular dynamics and the finite difference thermodynamic integration (FDTI) algorithm within the Discover program of MSI. We have shown that the orthogonality problem that occurs in the calculation of intraperturbed‐group contributions to the free energy is treated adequately by the FDTI method. We have also shown that problems of singularity and convergence in free energy calculations can be properly solved using this method. To conclude, the calculated free energies of hydration gave the following order of solvation for the candidates: p‐(2‐oxo‐1‐propyl)benzamidine > p‐methylbenzamidine > p‐ethylbenzamidine > p‐(1‐propyl)benzamidine > benzamidine. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001 相似文献
7.
We calculated the free energy of solvation of the neutral analogs of 18 amino acid side-chains (not including glycine and proline) using the OPLS all-atom force field in TIP4P water, SPC water, and cyclohexane by molecular dynamics simulation and thermodynamic integration. The average unsigned errors in the free energies of solvation in TIP4P, SPC, and cyclohexane are 4.4, 4.9, and 2.1 kJ/mol respectively. Most of the calculated hydration free energies are not favorable enough compared to experiment. The largest errors are found for tryptophan, histidine, glutamic acid, and glutamine. The average unsigned errors in the free energy of transfer from TIP4P to cyclohexane and from SPC to cyclohexane are 4.0 and 4.1 kJ/mol, respectively. The largest errors, of more than 7.5 kJ/mol, are found for histidine, glutamine, and glutamatic acid. 相似文献
8.
Present paper discusses the behaviour of transfer free energy of some amino acids from water to aqueous solution of Ethanol, 2‐PrOH ad t‐BuOH at different compositions. Dissection of transfer free energy into cavity term, interaction term and electrical term reveals that cavity forming free energy of transfer ΔGt0 (Cav.) plays an important role in dictating actual interaction of amino acids in these mixed solvents. Cavity forming free energy of transfer has been estimated by using Scaled Particle Theory (SPT). 相似文献
9.
Summary A theoretical conformational study was performed on leu-enkephalin in its zwitterionic form, both in vacuo and in the presence of a number, n, of up to 13 water molecules saturating its first hydration shell. The intramolecular energy of enkephalin as well as the intermolecular enkephalin-water and water-water interaction energies were computed with the SIBFA procedure (Sum of Interactions Between Fragments Ab initio computed), which uses additive ab initio multipole systematics and analytical formulas grounded on ab initio SCF computations. Energy minimizations were performed with a polyvalent minimizer, Merlin, with which three distinct derivative and three distinct nonderivative minimizers can be activated in a sequential fashion.Eight different candidate conformations of enkephalin were used as starting points. These conformations are either those found in distinct X-ray structures, or those proposed on the basis of theoretical computations by other authors. In the absence of hydration, they converged towards distinct folded energy-minima, the best four ones being separated by an energy gap of 8.7 kcal/mol. In marked contrast, with up to n=13, the energetical separation between the six best conformers narrowed down to 4 kcal/mol. They can be characterized by: (a) either a direct or a water-mediated ammonium-carboxylate interaction; b) either a close proximity (as in morphine) or a large separation between the aromatic rings of Tyr and Phe (intercenter separations of 4.5 Å and 10.5 Å, respectively), with each of the four mutual combinations of (a) and (b) being represented. 相似文献
10.
Masahiro Kinoshita Yuko Okamoto Fumio Hirata 《Journal of computational chemistry》1998,19(15):1724-1735
We developed a robust, highly efficient algorithm for solving the full reference interaction site model (RISM) equations for salt solutions near a solute molecule with many atomic sites. It was obtained as an extension of our previously reported algorithm for pure water near the solute molecule. The algorithm is a judicious hybrid of the Newton–Raphson and Picard methods. The most striking advantage is that the Jacobian matrix is just part of the input data and need not be recalculated at all. To illustrate the algorithm, we solved the full RISM equations for a dipeptide (NH2(SINGLE BOND)CHCH3(SINGLE BOND)CONH(SINGLE BOND)CHCH3(SINGLE BOND)COOH) in a 1 M NaCl solution. The extended simple point charge (SPC/E) model was employed for water molecules. Two different conformations of the dipeptide were considered. It was assumed for each conformation that the dipeptide was present either as an un-ionized form or as a zwitterion. The structure of the salt solution near the dipeptide and salt effects on the solvation free energy were also discussed. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1724–1735, 1998 相似文献
11.
Kyu‐Kwang Han Ki Hae Kim Byung Jin Mhin Hyeon S. Son 《Journal of computational chemistry》2001,22(10):1004-1009
We demonstrate the feasibility of using multiensemble sampling method (MESM) to determine the free energy difference between two far states for which the configurational distributions do not overlap at all. The MESM is a recently developed non‐Boltzmann sampling technique. The free energy of charging a sodium ion in water is accurately calculated in a single simulation, introducing nine intermediate ionic states. This is due to the ability of the method to explore the relevant parts of configuration space equally for every state, and this ability comes from the universality of weighting function W and the simplicity in adjusting its parameters. Detailed procedures of adjusting the parameters are presented. The comparison with a free energy perturbation method (FEPM) shows that the MESM is more reliable and efficient. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1004–1009, 2001 相似文献
12.
本文根据静电理论和晶体场理论将金属离子的水化分为水合和水溶剂化两个过程,水化能主要决定于这两个过程的能量变化,并由此推得水化能的计算公式。计算结果与实验值基本相符,从而证明这一模型是合理的。 相似文献
13.
Modern classical force fields have been traditionally parameterized by attempting to maximize agreement to any number of experimental and/or quantum mechanical target properties. As these force fields are pushed towards obtaining quantitative estimates of often subtle energetic differences, stringent and consistent parameterization criteria, particularly in regard to charge distributions, are required to ensure that systematic errors cancel, that parameters are transferable between molecules, and that performance does not significantly deteriorate when using more approximate methods, such as with continuum solvent models. Relative free energies of hydration are presented here for 40 mono- and disubstituted benzenes modeled with the OPLS-AA force field; heats of vaporization and pure liquid densities at standard conditions are presented when experimental data is available. Overall agreement between OPLS-AA and experiment is remarkable (average error = 0.5 kcal/mol for DeltaDeltaG(hydration), 1.0 kcal/mol for DeltaH(vap) (0), 0.02 g/mL for densities), yet several functional groups are identified as having consistent and correctable errors (alkyl-, nitro-, and thiobenzenes). Relative free energies of hydration obtained with rigorous free energy perturbations using explicit solvent are also compared with energies from minimizations using a generalized Born model (GB). There is high correlation between these estimates (R = 0.99), and as demonstrated here, reparameterization of the aforementioned groups can be guided with rapid GB calculations. 相似文献
14.
We develop a new method for calculating the hydration free energy (HFE) of a protein with any net charge. The polar part of the energetic component in the HFE is expressed as a linear combination of four geometric measures (GMs) of the protein structure and the generalized Born (GB) energy plus a constant. The other constituents in the HFE are expressed as linear combinations of the four GMs. The coefficients (including the constant) in the linear combinations are determined using the three‐dimensional reference interaction site model (3D‐RISM) theory applied to sufficiently many protein structures. Once the coefficients are determined, the HFE and its constituents of any other protein structure are obtained simply by calculating the four GMs and GB energy. Our method and the 3D‐RISM theory give perfectly correlated results. Nevertheless, the computation time required in our method is over four orders of magnitude shorter. 相似文献
15.
The development and parameterization of a solvent potential of mean force designed to reproduce the hydration thermodynamics of small molecules and macromolecules aimed toward applications in conformation prediction and ligand binding free energy prediction is presented. The model, named SGB/NP, is based on a parameterization of the Surface Generalized Born continuum dielectric electrostatic model using explicit solvent free energy perturbation calculations and a newly developed nonpolar hydration free energy estimator motivated by the results of explicit solvent simulations of the thermodynamics of hydration of hydrocarbons. The nonpolar model contains, in addition to the more commonly used solvent accessible surface area term, a component corresponding to the attractive solute-solvent interactions. This term is found to be important to improve the accuracy of the model, particularly for cyclic and hydrogen bonding compounds. The model is parameterized against the experimental hydration free energies of a set of small organic molecules. The model reproduces the experimental hydration free energies of small organic molecules with an accuracy comparable or superior to similar models employing more computationally demanding estimators and/or a more extensive set of parameters. 相似文献
16.
The evaluation of water binding free energies around solute molecules is important for the thermodynamic characterization of hydration or association processes. Here, a rapid approximate method to estimate water binding free energies around (bio)macromolecules from a single molecular dynamics simulation is presented. The basic idea is that endpoint free‐energy calculation methods are applied and the endpoint quantities are monitored on a three‐dimensional grid around the solute. Thus, a gridded map of water binding free energies around the solute is obtained, that is, from a single short simulation, a map of favorable and unfavorable water binding sites can be constructed. Among the employed free‐energy calculation methods, approaches involving endpoint information pertaining to actual thermodynamic integration calculations or endpoint information as exploited in the linear interaction energy method were examined. The accuracy of the approximate approaches was evaluated on the hydration of a cage‐like molecule representing either a nonpolar, polar, or charged water binding site and on α‐ and β‐cyclodextrin molecules. Among the tested approaches, the linear interaction energy method is considered the most viable approach. Applying the linear interaction energy method on the grid around the solute, a semi‐quantitative thermodynamic characterization of hydration around the whole solute is obtained. Disadvantages are the approximate nature of the method and a limited flexibility of the solute. © 2016 Wiley Periodicals, Inc. 相似文献
17.
18.
A protein energy surface is constructed. Validation is through applications of global energy minimization to surface loops of protein crystal structures. For 9 of 10 predictions, the native backbone conformation is identified correctly. Electrostatic energy is modeled as a pairwise sum of interactions between anisotropic atomic charge densities. Model repulsion energy has a softness similar to that seen in ab initio data. Intrinsic torsional energy is modeled as a sum over pairs of adjacent torsion angles of 2-dimensional Fourier series. Hydrophobic energy is that of a hydration shell model. The remainder of hydration free energy is obtained as the energetic effect of a continuous dielectric medium. Parameters are adjusted to reproduce the following data: a complete set of ab initio energy surfaces, meaning one for each pair of adjacent torsion angles of each blocked amino acid; experimental crystal structures and sublimation energies for nine model compounds; ab initio energies over 1014 conformations of 15 small-molecule dimers; and experimental hydration free energies for 48 model compounds. All ab initio data is at the Hartree–Fock/6–31G* level. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 548–573, 1998 相似文献
19.
Yvonne Paterson George Némethy Harold A. Scheraga 《Journal of solution chemistry》1982,11(12):831-856
A hydration-shell model has been developed for calculating the interaction energy between ions in water. The hydration shell around each ion contains a few tightly bound water molecules and a larger number of less tightly bound molecules. The energies of their interaction with the ion and the size of the hydration shell have been derived from published experimental data for ion-water clusters in the gas phase. An expression derived for the interaction energy of two univalent ions in water incorporates the following effects: a Lennard-Jones 6–12 interaction, a Coulomb interaction between the charges, the polarization of the hydration shells by a neighboring ion, and an energy term for the removal of water from the hydration shells when the hydration shells of two ions overlap. The effective dielectric constant at small ion-ion distances is the only adjustable parameter. Computed interaction energies for aqueous solutions of twelve alkali halides match experimental values, derived from conductimetric measurements, with an average error of ±14%. 相似文献
20.
We present results showing the importance of appropriate treatment of atomic masses in molecular dynamics (MD)-based single topology free-energy perturbations (FEPs) on small molecule systems. The reversibility of gas phase simulations is significantly improved by scaling the atomic mass of mutated atoms with the lambda variable normally used for the scaling of energy terms. Because this effect is less pronounced for solvated systems, it will not cancel in estimates of the relative hydration free energy difference. The advantage of mass scaling is demonstrated by a null mutation of ethane to ethane and the calculation of the relative hydration free energy difference between ethane and n-propane. Furthermore, it is found that the simulation time necessary for converged MD/FEPs is prohibitively large for relative hydration free energy calculations on cyclic alkanes. Therefore, we explore an alternative free energy pathway including strongly constrained conformations to improve convergence in FEP simulations of flexible molecules. 相似文献