Determination of Amino Acids' pKa: Importance of Cavity Scaling within Implicit Solvation Models and Choice of DFT Functionals.

Protonation states of molecules significantly influence the thermodynamics and kinetics of chemical reactions. This is especially important in biochemical processes, where appropriate protonation states of amino acids control the exo/endoergicity of practically all biochemical cycles. This paper is focused on appraisal of the impact of DFT functionals and PCM solvation models on the accuracy of pKa evaluations for all proteinogenic amino acids. Eight functionals (B3LYP, PBE0, revPBE0, M06-2X, M11, M11-L, TPSSh, and ωB97X-D) and four basis sets are considered, together with four kinds of implicit solvation models when additional attention is paid to a cavity construction. An influence of nonelectrostatic contributions and Wertz's corrections on Gibbs free energy is investigated together with accuracy of provided proton solvation energy. The best model is based on the M06-2X/6-311++G**/D-PCM/UAKS computational level. The fitting procedure is utilized to improve the accuracy of the evaluated models. All of these results are also compared with values obtained from the COSMOtherm program and CCSD(T) calculations. Results for cysteine and histidine are discussed individually, as they can be found in different protonation states at neutral pH.

Tryptophan to Tryptophan Hole Hopping in an Azurin Construct.

Electron transfer (ET) between neutral and cationic tryptophan residues in the azurin construct [ReI(H126)(CO)3(dmp)](W124)(W122)CuI (dmp = 4,7-Me2-1,10-phenanthroline) was investigated by Born-Oppenheimer quantum-mechanics/molecular mechanics/molecular dynamics (QM/MM/MD) simulations. We focused on W124•+ ← W122 ET, which is the middle step of the photochemical hole-hopping process *ReII(CO)3(dmp•-) ← W124 ← W122 ← CuI, where sequential hopping amounts to nearly 10,000-fold acceleration over single-step tunneling (ACS Cent. Sci. 2019, 5, 192-200). In accordance with experiments, UKS-DFT QM/MM/MD simulations identified forward and reverse steps of W124•+ ↔ W122 ET equilibrium, as well as back ET ReI(CO)3(dmp•-) → W124•+ that restores *ReII(CO)3(dmp•-). Strong electronic coupling between the two indoles (≥40 meV in the crossing region) makes the productive W124•+ ← W122 ET adiabatic. Energies of the two redox states are driven to degeneracy by fluctuations of the electrostatic potential at the two indoles, mainly caused by water solvation, with contributions from the protein dynamics in the W122 vicinity. ET probability depends on the orientation of Re(CO)3(dmp) relative to W124 and its rotation diminishes the hopping yield. Comparison with hole hopping in natural systems reveals structural and dynamics factors that are important for designing efficient hole-hopping processes.

Photoinduced hole hopping through tryptophans in proteins.

Hole hopping through tryptophan/tyrosine chains enables rapid unidirectional charge transport over long distances. We have elucidated structural and dynamical factors controlling hopping speed and efficiency in two modified azurin constructs that include a rhenium(I) sensitizer, Re(His)(CO)3(dmp)+, and one or two tryptophans (W1, W2). Experimental kinetics investigations showed that the two closely spaced (3 to 4 Å) intervening tryptophans dramatically accelerated long-range electron transfer (ET) from CuI to the photoexcited sensitizer. In our theoretical work, we found that time-dependent density-functional theory (TDDFT) quantum mechanics/molecular mechanics/molecular dynamics (QM/MM/MD) trajectories of low-lying triplet excited states of ReI(His)(CO)3(dmp)+-W1(-W2) exhibited crossings between sensitizer-localized (*Re) and charge-separated [ReI(His)(CO)3(dmp•-)/(W1•+ or W2•+)] (CS1 or CS2) states. Our analysis revealed that the distances, angles, and mutual orientations of ET-active cofactors fluctuate in a relatively narrow range in which the cofactors are strongly coupled, enabling adiabatic ET. Water-dominated electrostatic field fluctuations bring *Re and CS1 states to a crossing where *Re(CO)3(dmp)+←W1 ET occurs, and CS1 becomes the lowest triplet state. ET is promoted by solvation dynamics around *Re(CO)3(dmp)+(W1); and CS1 is stabilized by Re(dmp•-)/W1•+ electron/hole interaction and enhanced W1•+ solvation. The second hop, W1•+←W2, is facilitated by water fluctuations near the W1/W2 unit, taking place when the electrostatic potential at W2 drops well below that at W1•+ Insufficient solvation and reorganization around W2 make W1•+←W2 ET endergonic, shifting the equilibrium toward W1•+ and decreasing the charge-separation yield. We suggest that multiscale TDDFT/MM/MD is a suitable technique to model the simultaneous evolution of photogenerated excited-state manifolds.

QM and QM/MM umbrella sampling MD study of the formation of Hg(II)-thymine bond: Model for evaluation of the reaction energy profiles in solutions with constant pH.

The formation of the Hg-N3(T) bond between the 1-methylthymine (T) molecule and the hydrated Hg2+ cation was explored with the combined quantum mechanics/molecular mechanics (QM/MM) method including Free Energy Perturbation corrections. The thermodynamic properties were determined in the whole pH range, when these systems were explicitly investigated and considered as the QM part: (1) T + [Hg(H2 O)6 ]2+ , (2) T + [Hg(H2 O)5 (OH)]+ , (3) T + Hg(H2 O)4 (OH)2 , and (4) N3-deprotonated T + Hg(H2 O)4 (OH)2 . The MM part contained only solvent molecules and counterions. As a result, the dependence of Gibbs-Alberty reaction free energy on pH was obtained along the reaction coordinate. We found that an endoergic reaction in acidic condition up to pH < 4-5 becomes exoergic for a higher pH corresponding to neutral and basic solutions. The migration of the Hg2+ cation between N3 and O4/2 positions in dependence on pH is discussed as well. For the verification, DFT calculations of stationary points were performed confirming the qualitative trends of QM/MM MD simulations and NMR parameters were determined for them.

Protein environment affects the water-tryptophan binding mode. MD, QM/MM, and NMR studies of engrailed homeodomain mutants.

Water molecules can interact with aromatic moieties using either their O-H bonds or their lone-pairs of electrons. In proteins, water-π interactions have been reported to occur with tryptophan and histidine residues, and dynamic exchange between O-Hπ hydrogen bonding and lone-pairπ interactions was suggested to take place, based on ab initio calculations. Here we used classical and QM/MM molecular dynamics simulations, complemented with an NMR study, to examine a specific water-indole interaction observed in the engrailed homeodomain and in its mutants. Our simulations indicate that the binding mode between water and indole can adapt to the potential created by the surrounding amino acids (and by the residues at the DNA surface in protein-DNA complexes), and support the model of dynamic switching between the O-Hπ hydrogen bonding and lone-pairπ binding modes.

Redox Potentials for Tetraplatin, Satraplatin, Its Derivatives, and Ascorbic Acid: A Computational Study.

Redox potentials of the Pt(IV) complexes, such as satraplatin, tetraplatin, and several others, are determined at the density functional theory (DFT) level (with B3LYP, ω-B97XD, PBE1PBE, TPSSTPSS, M06-L, M11-L, and MN12-L functionals) and compared with post-Hartree-Fock methods MP2 and CCSD(T). Calculations are performed in water solution employing an implicit solvation model. The impact of replacement of a chloro ligand by a water molecule (hydration in the equatorial plane of the complexes) is also explored. Furthermore, an influence of solvent pH on the magnitude of the redox potentials is discussed for such hydrated complexes. The obtained results are compared with available experimental data leading to a root-mean-square deviation (RMSD) of ca. 0.23 V, using the CCSD(T)/6-31+G(d)/IEF-PCM/scaled-UAKS level. Distribution of the electron density is analyzed at the B3LYP/6-311++G(2df,2pd) level. Also, a correlation between binding energies of axial ligands and the redox potential is demonstrated. Since the Pt(IV) complexes are considered in the framework of anticancer treatment, possible reducing agents in bioenvironment are searched. From this reason, the reduction potential of different protonation states of ascorbic acid is also presented.

The influence of the metal cations and microhydration on the reaction trajectory of the N3 ↔ O2 thymine proton transfer: Quantum mechanical study.

This study involves the intramolecular proton transfer (PT) process on a thymine nucleobase between N3 and O2 atoms. We explore a mechanism for the PT assisted by hexacoordinated divalent metals cations, namely Mg2+ , Zn2+ , and Hg2+ . Our results point out that this reaction corresponds to a two-stage process. The first involves the PT from one of the aqua ligands toward O2. The implications of this stage are the formation of a hydroxo anion bound to the metal center and a positively charged thymine. To proceed to the second stage, a structural change is needed to allow the negatively charged hydroxo ligand to abstract the N3 proton, which represents the final product of the PT reaction. In the presence of the selected hexaaqua cations, the activation barrier is at most 8 kcal/mol. © 2017 Wiley Periodicals, Inc.

Study on electronic properties, thermodynamic and kinetic parameters of the selected platinum(II) derivatives interacting with guanine.

Interaction of hydrated forms of several potential anticancer agents (PtCl2(diaminocyclohexane), trans-[PtCl2(NH3)(thiazole)], cis-[PtCl2(NH3)(piperidine)], and cis-PtCl2(NH3)(cyclohexylamine) complexes) with guanine are explored and compared with an analogous interaction of cisplatin. Basic electronic properties, binding and stabilization energies are determined and energy profiles for the aquation reaction are estimated at the B3LYP/6-311++G(2df,2pd) level of theory. It is found that the substitution reaction is an exothermic and exergonic process with ΔG slightly less negative than -20kcal/mol. The largest energy release occurs for PtCl(H2O)(diaminocyclohexane) complex. The rate constants for the Pt(II) complexes in the chloro- and hydroxo-form are compared and an impact of the ligand in the trans position to water is discussed.

Side Reactions with an Equilibrium Constraint: Detailed Mechanism of the Substitution Reaction of Tetraplatin with dGMP as a Starting Step of the Platinum(IV) Reduction Process.

Two possible pathways of the substitution reaction within the reduction process of the PtIV(DACH)Cl4 by dGMP are compared: associative reaction course and autocatalytic Basolo-Pearson mechanisms. Since two forms: single-protonated and fully deprotonated phosphate group of dGMP are present in equilibrium at neutral and mildly acidic solutions, consideration of a side reactions scheme with acido-basic equilibrium-constraint is a very important model for obtaining reliable results. The examined complexes are optimized at the B3LYP-GD3BJ/6-31G(d) level with the COSMO implicit solvation model and Klamt's radii used for cavity construction. Energy characteristics and thermodynamics for all reaction branches are determined using the B3LYP-GD3BJ/6-311++G(2df,2pd)/IEF-PCM/scaled-UAKS level with Wertz's entropy corrections. Rate constants are estimated for each individual branch according to Eyring's transition state theory (TST), averaged according to equilibrium constraint and compared with available experimental data. The determined reaction barriers of the autocatalytic pathway fairly correspond with experimental values. Furthermore, autocatalytic reaction of tetraplatin and its two analogues complexes [PtIV(en)Cl4 and PtIV(NH3)2Cl4] are explored and compared with measured data in order to examined general reaction descriptors.

Estimation of Transition-Metal Empirical Parameters for Molecular Mechanical Force Fields.

Force-field parameters of the first row transition metals together with a few additional common elements such as those from the second (Rh, Ru) and third (Hg, Pt) rows of elements in ligated forms were determined based on the density functional theory calculations. Bonding characteristics were determined by averaging metal-ligand force constants in optimal geometries from several chosen complexes of each metal in the most common oxidation numbers and structural arrangements. Parameters of Lennard-Jones potential were determined based on a supermolecular model. Our determined molecular mechanical parameters are compared with presently available parameters published by other groups. We performed two different kinds of testing in order to demonstrate the reliability of these parameters in the case of ligated metallo complexes. First, the nonbonding potential was constructed for an additional set of 19 larger systems containing common complexes with organic molecules. The second test compares the Pt-O and Pt-H radial distribution functions for cisplatin in a box of TIP3P water with lately published studies.

Reduction Process of Tetraplatin in the Presence of Deoxyguanosine Monophosphate (dGMP): A Computational DFT Study.

The reduction mechanism of [Pt(IV) (dach)Cl4 ] (dach=diaminocyclohexyl) in the presence of dGMP was studied. The first step is substitution of a chloro ligand by dGMP, followed by nucleophilic attack of a phosphate or sugar oxygen atom to the C8-position of guanine. Subsequent reduction forms the [Pt(II) (dach)Cl2 ] complex. The whole process is completed by a hydrolysis. Two different pathways for the substitution reaction were examined: a direct associative and a Basolo-Pearson autocatalytic mechanism. All the explored structures were optimized at the B3LYP-D3/6-31G(d) level and by using the COSMO solvation model with Klamt's radii. Single-point energetics was determined at the B3LYP-GD3BJ/6-311++G(2df,2pd)/PCM/scaled-UAKS level. Activation barriers were used for an estimation of the rate constants and these were compared with experimental values. It was found that the rate-determining step is the nucleophilic attack with a slightly faster performance in the 3'-dGMP branch than in the case of 5'-dGMP with activation barriers of 21.1 and 20.4 kcal mol(-1) (experimental: 23.8 and 23.2 kcal mol(-1) ). The reduction reaction is connected with an electron flow from guanine. The product of the reduction reaction is a chelate structure, which dissociates within the last reaction step, that is, a hydrolysis reaction. The whole redox process (substitution, reduction, and hydrolysis) is exergonic by 34 and 28 kcal mol(-1) for 5'-dGMP and 3'-dGMP, respectively.