Benchmarks of the density functional tight-binding method for redox, protonation and electronic properties of quinones.

Organic materials with controllable molecular design and sustainable resources are promising electrode materials. Crystalline quinones have been investigated in a variety of rechargeable battery chemistries due to their ubiquitous nature, voltage tunability and environmental friendliness. In acidic electrolytes, quinone crystals can undergo proton-coupled electron transfer (PCET), resulting in charge storage. However, the detailed mechanism of this phenomenon remains elusive. To model PCET in crystalline quinones, force field-based methods are not viable due to variable redox states of the quinone molecules during battery operation and computationally efficient quantum mechanical methods are strongly desired. The semi-empirical density functional tight-binding (DFTB) method has been widely used to study inorganic crystalline systems and biological systems but has not been comprehensively benchmarked for studying charge transport in quinones. In this work, we benchmark the third order variant of DFTB (DFTB3) for the reduction potential of quinones in aqueous solution, energetics of proton transfer between quinones and between quinones and water, and structural and electronic properties of crystalline quinones. Our results reveal the deficiencies of the DFTB3 method in describing the proton affinity of quinones and the structural and electronic properties of crystalline quinones, and highlight the need for further development of the DFTB method for describing charge transport in crystalline quinones.

Role of the CarH photoreceptor protein environment in the modulation of cobalamin photochemistry.

The photochemistry of cobalamins has recently been found to have biological importance, with the discovery of bacterial photoreceptor proteins, such as CarH and AerR. CarH and AerR, are involved in the light regulation of carotenoid biosynthesis and bacteriochlorophyll biosynthesis, respectively, in bacteria. Experimental transient absorption spectroscopic studies have indicated unusual photochemical behavior of 5'-deoxy-5'-adenosylcobalamin (AdoCbl) in CarH, with excited-state charge separation between cobalt and adenosyl and possible heterolytic cleavage of the Co-adenosyl bond, as opposed to the homolytic cleavage observed in aqueous solution and in many AdoCbl-based enzymes. We employ molecular dynamics and hybrid quantum mechanical/molecular mechanical calculations to obtain a microscopic understanding of the modulation of the excited electronic states of AdoCbl by the CarH protein environment, in contrast to aqueous solution and AdoCbl-based enzymes. Our results indicate a progressive stabilization of the electronic states involving charge transfer (CT) from cobalt/corrin to adenine on changing the environment from gas phase to water to solvated CarH. The solvent exposure of the adenosyl ligand in CarH, the π-stacking interaction between a tryptophan and the adenine moiety, and the hydrogen-bonding interaction between a glutamate and the lower axial ligand of cobalt are found to contribute to the stabilization of the states involving CT to adenine. The combination of these three factors, the latter two of which can be experimentally tested via mutagenesis studies, is absent in an aqueous solvent environment and in AdoCbl-based enzymes. The favored CT from metal and/or corrin to adenine in CarH may promote heterolytic cleavage of the cobalt-adenosyl bond proposed by experimental studies. Overall, this work provides novel, to our knowledge, physical insights into the mechanism of CarH function and directions for future experimental investigations. The fundamental understanding of the mechanism of CarH functioning will serve the development of optogenetic tools based on the new class of B12-dependent photoreceptors.

Mechanism and potential sites of potassium interaction with glutamate transporters.

In the mammalian glutamate transporters, countertransported intracellular K+ is essential for relocating the glutamate binding site to the extracellular side of the membrane. This K+-dependent process is believed to be rate limiting for the transport cycle. In contrast, extracellular K+ induces glutamate release upon transporter reversal. Here, we analyzed potential K+ binding sites using molecular dynamics (MD) simulations and site-directed mutagenesis. Two candidate sites were identified by spontaneous K+ binding in MD simulations, one site (K1 site) overlapping with the Na1 Na+ binding site and the K2 site being localized under hairpin loop 2 (HP2). Mutations to conserved amino acid residues in these sites resulted in several transporters that were defective in K+-induced reverse transport and which bound K+ with reduced apparent affinity compared with the wild-type transporter. However, external K+ interaction was abolished in only one mutant transporter EAAC1D454A in the K1 site. Our results, for the first time, directly demonstrate effects of K1-site mutations on K+ binding, in contrast to previous reports on K+ binding sites based on indirect evidence. We propose that K+ binding to the K1 site is responsible for catalyzing the relocation step, whereas binding to the K2 site may have an as-of-yet unidentified regulatory function.

Coiled-coil registry shifts in the F684I mutant of Bicaudal D result in cargo-independent activation of dynein motility.

The dynein adaptor Drosophila Bicaudal D (BicD) is auto-inhibited and activates dynein motility only after cargo is bound, but the underlying mechanism is elusive. In contrast, we show that the full-length BicD/F684I mutant activates dynein processivity even in the absence of cargo. Our X-ray structure of the C-terminal domain of the BicD/F684I mutant reveals a coiled-coil registry shift; in the N-terminal region, the two helices of the homodimer are aligned, whereas they are vertically shifted in the wild-type. One chain is partially disordered and this structural flexibility is confirmed by computations, which reveal that the mutant transitions back and forth between the two registries. We propose that a coiled-coil registry shift upon cargo-binding activates BicD for dynein recruitment. Moreover, the human homolog BicD2/F743I exhibits diminished binding of cargo adaptor Nup358, implying that a coiled-coil registry shift may be a mechanism to modulate cargo selection for BicD2-dependent transport pathways.

Role of Coiled-Coil Registry Shifts in the Activation of Human Bicaudal D2 for Dynein Recruitment upon Cargo Binding.

Dynein adaptors such as Bicaudal D2 (BicD2) recognize cargo for dynein-dependent transport, and cargo-bound adaptors are required to activate dynein for processive transport, but the mechanism of action is unknown. Here we report the X-ray structure of the cargo-binding domain of human BicD2 and investigate the structural dynamics of the coiled-coil. Our molecular dynamics simulations support the fact that BicD2 can switch from a homotypic coiled-coil registry, in which both helices of the homodimer are aligned, to an asymmetric registry, where a portion of one helix is vertically shifted, as both states are similarly stable and defined by distinct conformations of F743. The F743I variant increases dynein recruitment in the Drosophila homologue, whereas the human R747C variant causes spinal muscular atrophy. We report spontaneous registry shifts for both variants, which may be the cause for BicD2 hyperactivation and disease. We propose that a registry shift upon cargo binding may activate autoinhibited BicD2 for dynein recruitment.

Sterol A-ring plasticity in hedgehog protein cholesterolysis supports a primitive substrate selectivity mechanism.

Cholesterolysis of Hedgehog family proteins couples endoproteolysis to protein C-terminal sterylation. The transformation is self-catalyzed by HhC, a partially characterized enzymatic domain found in precursor forms of Hedgehog. Here we explore spatial ambiguity in sterol recognition by HhC, using a trio of derivatives where the sterol A-ring is contracted, fused, or distorted. Sterylation assays indicate that these geometric variants react as substrates with relative activity: cholesterol, 1.000 > A-ring contracted, 0.100 > A-ring fused, 0.020 > A-ring distorted, 0.005. Experimental results and computational sterol docking into the first HhC homology model suggest a partially unstructured binding site with substrate recognition governed in large part by hydrophobic interactions.

Role of active site conformational changes in photocycle activation of the AppA BLUF photoreceptor.

Blue light using flavin adenine dinucleotide (BLUF) proteins are essential for the light regulation of a variety of physiologically important processes and serve as a prototype for photoinduced proton-coupled electron transfer (PCET). Free-energy simulations elucidate the active site conformations in the AppA (activation of photopigment and puc expression) BLUF domain before and following photoexcitation. The free-energy profile for interconversion between conformations with either Trp104 or Met106 closer to the flavin, denoted Trpin/Metout and Trpout/Metin, reveals that both conformations are sampled on the ground state, with the former thermodynamically favorable by ∼3 kcal/mol. These results are consistent with the experimental observation of both conformations. To analyze the proton relay from Tyr21 to the flavin via Gln63, the free-energy profiles for Gln63 rotation were calculated on the ground state, the locally excited state of the flavin, and the charge-transfer state associated with electron transfer from Tyr21 to the flavin. For the Trpin/Metout conformation, the hydrogen-bonding pattern conducive to the proton relay is not thermodynamically favorable on the ground state but becomes more favorable, corresponding to approximately half of the configurations sampled, on the locally excited state. The calculated energy gaps between the locally excited and charge-transfer states suggest that electron transfer from Tyr21 to the flavin is more facile for configurations conducive to proton transfer. When the active site conformation is not conducive to PCET from Tyr21, Trp104 can directly compete with Tyr21 for electron transfer to the flavin through a nonproductive pathway, impeding the signaling efficiency.

Assessing the Potential Effects of Active Site Mg2+ Ions in the glmS Ribozyme-Cofactor Complex.

Ribozymes employ diverse catalytic strategies in their self-cleavage mechanisms, including the use of divalent metal ions. This work explores the effects of Mg2+ ions in the active site of the glmS ribozyme-GlcN6P cofactor complex using computational methods. Deleterious and potentially beneficial effects of an active site Mg2+ ion on the self-cleavage reaction were identified. The presence of a Mg2+ ion near the scissile phosphate oxygen atoms at the cleavage site was determined to be deleterious, and thereby anticatalytic, due to electrostatic repulsion of the cofactor, disruption of key hydrogen-bonding interactions, and obstruction of nucleophilic attack. On the other hand, the presence of a Mg2+ ion at another position in the active site, the Hoogsteen face of the putative base, was found to avoid these deleterious effects and to be potentially catalytically favorable owing to the stabilization of negative charge and pKa shifting of the guanine base.

Effect of Protonation upon Electronic Coupling in the Mixed Valence and Mixed Protonated Complex, [Ni(2,3-pyrazinedithiol)2].

We demonstrate that protonation of a mixed valence molecule, generating a mixed valence mixed protonated (MVMP) state, results in a severe reduction in the electronic coupling intimately connected with electron transfer kinetics. This phenomenon is illustrated by synthesizing a mixed valence molecule, [Ni(2,3-pyrazinedithiol)2], that can be asymmetrically protonated, rendering the MVMP state. We characterize the structural, electronic, vibrational, and magnetic properties of this complex in five different states, including the mixed valence and MVMP states, and then analyze the intervalence charge transfer (IVCT) band to demonstrate a five-fold reduction in electronic coupling upon protonation. We conclude that the reduction in electronic coupling is a result of the asymmetry of the electronic orbitals of the redox sites that results from the asymmetric protonation. This conclusion suggests that many systems designed to link electron and proton transfer will also exhibit a decrease in electronic coupling upon protonation as the strength of the interaction between redox and protonation sites is increased.

Proton Quantization and Vibrational Relaxation in Nonadiabatic Dynamics of Photoinduced Proton-Coupled Electron Transfer in a Solvated Phenol-Amine Complex.

Nonadiabatic dynamics simulations of photoinduced proton-coupled electron transfer (PCET) in a phenol-amine complex in solution were performed. The electronic potential energy surfaces were generated on-the-fly with a hybrid quantum mechanical/molecular mechanical approach that described the solute with a multiconfigurational method in a bath of explicit solvent molecules. The transferring hydrogen nucleus was represented as a quantum mechanical wave function calculated with grid-based methods, and surface hopping trajectories were propagated on the adiabatic electron-proton vibronic surfaces. Following photoexcitation to the excited S1 electronic state, the overall decay to the ground vibronic state was found to be comprised of relatively fast decay from a lower proton vibrational state of S1 to a highly excited proton vibrational state of the ground S0 electronic state, followed by vibrational relaxation within the S0 state. Proton transfer can occur either on the highly excited proton vibrational states of S0 due to small environmental fluctuations that shift the delocalized vibrational wave functions or on the low-energy proton vibrational states of S1 due to solvent reorganization that alters the asymmetry of the proton potential and reduces the proton transfer barrier. The isotope effect arising from replacing the transferring hydrogen with deuterium is predicted to be negligible because hydrogen and deuterium behave similarly in both types of proton transfer processes. Although an isotope effect could be observed for other systems, in general the absence of an isotope effect does not imply the absence of proton transfer in photoinduced PCET systems. This computational approach is applicable to a wide range of other photoinduced PCET processes.

Copper Oxidation/Reduction in Water and Protein: Studies with DFTB3/MM and VALBOND Molecular Dynamics Simulations.

We apply two recently developed computational methods, DFTB3 and VALBOND, to study copper oxidation/reduction processes in solution and protein. The properties of interest include the coordination structure of copper in different oxidation states in water or in a protein (plastocyanin) active site, the reduction potential of the copper ion in different environments, and the environmental response to copper oxidation. The DFTB3/MM and VALBOND simulation results are compared to DFT/MM simulations and experimental results whenever possible. For a copper ion in aqueous solution, DFTB3/MM results are generally close to B3LYP/MM with a medium basis, including both solvation structure and reduction potential for Cu(II); for Cu(I), however, DFTB3/MM finds a two-water coordination, similar to previous Born-Oppenheimer molecular dynamics simulations using BLYP and HSE, whereas B3LYP/MM leads to a tetrahedron coordination. For a tetraammonia copper complex in aqueous solution, VALBOND and DFTB3/MM are consistent in terms of both structural and dynamical properties of solvent near copper for both oxidation states. For copper reduction in plastocyanin, DFTB3/MM simulations capture the key properties of the active site, and the computed reduction potential and reorganization energy are in fair agreement with experiment, especially when the periodic boundary condition is used. Overall, the study supports the value of VALBOND and DFTB3(/MM) for the analysis of fundamental copper redox chemistry in water and protein, and the results also help highlight areas where further improvements in these methods are desirable.

DFTB3 Parametrization for Copper: The Importance of Orbital Angular Momentum Dependence of Hubbard Parameters.

We report the parametrization of a density functional tight binding method (DFTB3) for copper in a spin-polarized formulation. The parametrization is consistent with the framework of 3OB for main group elements (ONCHPS) and can be readily used for biological applications that involve copper proteins/peptides. The key to our parametrization is to introduce orbital angular momentum dependence of the Hubbard parameter and its charge derivative, thus allowing the 3d and 4s orbitals to adopt different sizes and responses to the change of charge state. The parametrization has been tested by applying to a fairly broad set of molecules of biological relevance, and the properties of interest include optimized geometries, ligand binding energies, and ligand proton affinities. Compared to the reference QM level (B3LYP/aug-cc-pVTZ, which is shown here to be similar to the B97-1 and CCSD(T) results, in terms of many properties of interest for a set of small copper containing molecules), our parametrization generally gives reliable structural properties for both Cu(I) and Cu(II) compounds, although several exceptions are also noted. For energetics, the results are more accurate for neutral ligands than for charged ligands, likely reflecting the minimal basis limitation of DFTB3; the results generally outperform NDDO based methods such as PM6 and even PBE with the 6-31+G(d,p) basis. For all ligand types, single-point B3LYP calculations at DFTB3 geometries give results very close (∼1-2 kcal/mol) to the reference B3LYP values, highlighting the consistency between DFTB3 and B3LYP structures. Possible further developments of the DFTB3 model for a better treatment of transition-metal ions are also discussed. In the current form, our first generation of DFTB3 copper model is expected to be particularly valuable as a method that drives sampling in systems that feature a dynamical copper binding site.

Role of Solvent Dynamics in Photoinduced Proton-Coupled Electron Transfer in a Phenol-Amine Complex in Solution.

Photoinduced proton-coupled electron transfer (PCET) plays an essential role in a wide range of energy conversion processes. Previous experiments on a phenol-amine complex in solution provided evidence of an electron-proton transfer (EPT) excited state characterized by both intramolecular charge transfer and proton transfer from the phenol to the amine. Herein we analyze hundreds of surface hopping trajectories to investigate the role of solvent dynamics following photoexcitation to the EPT state. This solvent dynamics leads to a significant decrease in the energy gap between the ground and EPT states, thereby facilitating decay to the ground state, and generates an electrostatic environment conducive to proton transfer on the EPT state. In addition to solvent reorganization, the geometrical properties at the hydrogen-bonding interface must be suitable to allow proton transfer. These mechanistic insights elucidate the underlying fundamental physical principles of photoinduced PCET processes.

Microscopic basis for kinetic gating in Cytochrome c oxidase: insights from QM/MM analysis.

Understanding the mechanism of vectorial proton pumping in biomolecules requires establishing the microscopic basis for the regulation of both thermodynamic and kinetic features of the relevant proton transfer steps. For the proton pump cytochrome c oxidase, while the regulation of thermodynamic driving force for key proton transfers has been discussed in great detail, the microscopic basis for the control of proton transfer kinetics has been poorly understood. Here we carry out extensive QM/MM free energy simulations to probe the kinetics of relevant proton transfer steps and analyze the effects of local structure and hydration level. We show that protonation of the proton loading site (PLS, taken to be a propionate of heme a3) requires a concerted process in which a key glutamic acid (Glu286H) delivers the proton to the PLS while being reprotonated by an excess proton coming from the D-channel. The concerted nature of the mechanism is a crucial feature that enables the loading of the PLS before the cavity containing Glu286 is better hydrated to lower its pK a to experimentally measured range; the charged rather than dipolar nature of the process also ensures a tight coupling with heme a reduction, as emphasized by Siegbahn and Blomberg. In addition, we find that rotational flexibility of the PLS allows its protonation before that of the binuclear center (the site where oxygen gets reduced to water). Together with our recent study (P. Goyal, et al., Proc. Natl. Acad. Sci. USA, 110:18886-18891, 2013) that focused on the modulation of Glu286 pK a , the current work suggests a mechanism that builds in a natural sequence for the protonation of the PLS prior to that of the binuclear center. This provides microscopic support to the kinetic constraints revealed by kinetic network analysis as essential elements that ensure an efficient vectorial proton transport in cytochrome c oxidase.

Role of the active site guanine in the glmS ribozyme self-cleavage mechanism: quantum mechanical/molecular mechanical free energy simulations.

The glmS ribozyme catalyzes a self-cleavage reaction at the phosphodiester bond between residues A-1 and G1. This reaction is thought to occur by an acid-base mechanism involving the glucosamine-6-phosphate cofactor and G40 residue. Herein quantum mechanical/molecular mechanical free energy simulations and pKa calculations, as well as experimental measurements of the rate constant for self-cleavage, are utilized to elucidate the mechanism, particularly the role of G40. Our calculations suggest that an external base deprotonates either G40(N1) or possibly A-1(O2'), which would be followed by proton transfer from G40(N1) to A-1(O2'). After this initial deprotonation, A-1(O2') starts attacking the phosphate as a hydroxyl group, which is hydrogen-bonded to deprotonated G40, concurrent with G40(N1) moving closer to the hydroxyl group and directing the in-line attack. Proton transfer from A-1(O2') to G40 is concomitant with attack of the scissile phosphate, followed by the remainder of the cleavage reaction. A mechanism in which an external base does not participate, but rather the proton transfers from A-1(O2') to a nonbridging oxygen during nucleophilic attack, was also considered but deemed to be less likely due to its higher effective free energy barrier. The calculated rate constant for the favored mechanism is in agreement with the experimental rate constant measured at biological Mg(2+) ion concentration. According to these calculations, catalysis is optimal when G40 has an elevated pKa rather than a pKa shifted toward neutrality, although a balance among the pKa's of A-1, G40, and the nonbridging oxygen is essential. These results have general implications, as the hammerhead, hairpin, and twister ribozymes have guanines at a similar position as G40.

Nonadiabatic dynamics of photoinduced proton-coupled electron transfer in a solvated phenol-amine complex.

Photoinduced concerted electron-proton transfer (EPT), denoted photo-EPT, is important for a wide range of energy conversion processes. Transient absorption and Raman spectroscopy experiments on the hydrogen-bonded p-nitrophenylphenol-t-butylamine complex, solvated in 1,2-dichloroethane, suggested that this complex may undergo photo-EPT. The experiments probed two excited electronic states that were interpreted as an intramolecular charge transfer (ICT) state and an EPT state. Herein mixed quantum mechanical/molecular mechanical nonadiabatic surface hopping dynamics is used to investigate the relaxation pathways following photoexcitation. The potential energy surface is generated on the fly with a semiempirical floating occupation molecular orbital complete active space configuration interaction method for the solute molecule and a molecular mechanical force field for the explicit solvent molecules. The free energy curves along the proton transfer coordinate illustrate that proton transfer is thermodynamically and kinetically favorable on the lower-energy excited state but not on the higher-energy excited state, supporting the characterization of these states as EPT and ICT, respectively. The nonadiabatic dynamics simulations indicate that the population decays from the ICT state to the EPT state in ∼100 fs and from the EPT state to the ground state on the slower time scale of ∼1 ps, qualitatively consistent with the experimental measurements. For ∼54% of the trajectories, the proton transfers from the phenol to the amine in ∼400 fs on the EPT state and then transfers back to the phenol rapidly upon decay to the ground state. Thus, these calculations augment the original interpretation of the experimental data by providing evidence of proton transfer on the EPT state prior to decay to the ground state. The fundamental insights obtained from these simulations are also relevant to other photo-EPT processes.

Molecular simulation of water and hydration effects in different environments: challenges and developments for DFTB based models.

We discuss the description of water and hydration effects that employs an approximate density functional theory, DFTB3, in either a full QM or QM/MM framework. The goal is to explore, with the current formulation of DFTB3, the performance of this method for treating water in different chemical environments, the magnitude and nature of changes required to improve its performance, and factors that dictate its applicability to reactions in the condensed phase in a QM/MM framework. A relatively minor change (on the scale of kBT) in the O-H repulsive potential is observed to substantially improve the structural properties of bulk water under ambient conditions; modest improvements are also seen in dynamic properties of bulk water. This simple change also improves the description of protonated water clusters, a solvated proton, and to a more limited degree, a solvated hydroxide. By comparing results from DFTB3 models that differ in the description of water, we confirm that proton transfer energetics are adequately described by the standard DFTB3/3OB model for meaningful mechanistic analyses. For QM/MM applications, a robust parametrization of QM-MM interactions requires an explicit consideration of condensed phase properties, for which an efficient sampling technique was developed recently and is reviewed here. The discussions help make clear the value and limitations of DFTB3 based simulations, as well as the developments needed to further improve the accuracy and transferability of the methodology.

Changing hydration level in an internal cavity modulates the proton affinity of a key glutamate in cytochrome c oxidase.

Cytochrome c oxidase contributes to the transmembrane proton gradient by removing two protons from the high-pH side of the membrane each time the binuclear center active site is reduced. One proton goes to the binuclear center, whereas the other is pumped to the low-pH periplasmic space. Glutamate 286 (Glu286) has been proposed to serve as a transiently deprotonated proton donor. Using unrestrained atomistic molecular dynamics simulations, we show that the size of and water distribution in the hydrophobic cavity that holds Glu286 is controlled by the protonation state of the propionic acid of heme a3, a group on the proton outlet pathway. Protonation of the propionate disrupts hydrogen bonding to two side chains, allowing a loop to swing open. Continuum electrostatics and atomistic free-energy perturbation calculations show that the resultant changes in hydration and electrostatic interactions lower the Glu proton affinity by at least 5 kcal/mol. These changes in the internal hydration level occur in the absence of major conformational transitions and serve to stabilize needed transient intermediates in proton transport. The trigger is not the protonation of the Glu of interest, but rather the protonation of a residue ∼10 Šaway. Thus, unlike local water penetration to stabilize a new charge, this finding represents a specific role for water molecules in the protein interior, mediating proton transfers and facilitating ion transport.

Proton storage site in bacteriorhodopsin: new insights from quantum mechanics/molecular mechanics simulations of microscopic pK(a) and infrared spectra.

Identifying the group that acts as the proton storage/loading site is a challenging but important problem for understanding the mechanism of proton pumping in biomolecular proton pumps, such as bacteriorhodopsin (bR) and cytochrome c oxidase. Recent experimental studies of bR propelled the idea that the proton storage/release group (PRG) in bR is not an amino acid but a water cluster embedded in the protein. We argue that this idea is at odds with our knowledge of protein electrostatics, since invoking the water cluster as the PRG would require the protein to raise the pK(a) of a hydronium by almost 11 pK(a) units, which is difficult considering known cases of pK(a) shifts in proteins. Our recent quantum mechanics/molecular mechanics (QM/MM) simulations suggested an alternative "intermolecular proton bond" model in which the stored proton is shared between two conserved Glu residues (194 and 204). Here we show that this model leads to microscopic pK(a) values consistent with available experimental data and the functional requirement of a PRG. Extensive QM/MM simulations also show that, independent of a number of technical issues, such as the influence of QM region size, starting X-ray structure, and nuclear quantum effects, the "intermolecular proton bond" model is qualitatively consistent with available spectroscopic data. Potential of mean force calculations show explicitly that the stored proton strongly prefers the pair of Glu residues over the water cluster. The results and analyses help highlight the importance of considering protein electrostatics and provide arguments for why the "intermolecular proton bond" model is likely applicable to the PRG in biomolecular proton pumps in general.

Application of the SCC-DFTB method to neutral and protonated water clusters and bulk water.

The self-consistent charge density functional tight-binding (SCC-DFTB) method has been actively employed to study proton-transfer processes in biological systems. Recent studies in the literature employing SCC-DFTB reported that the method favors the Zundel form of the hydrated proton over the Eigen form, both in gas-phase water clusters and in bulk water, in disagreement with both higher-level calculations and experimental data. In this work, we explore the performance of SCC-DFTB for protonated gas-phase water clusters and bulk water (the latter both with and without an excess proton) with a modified O-H repulsive potential reported in our earlier work and with on-site third-order expansion of the DFT energy. Our results show that, with the proper set of published parameters, SCC-DFTB does correctly favor the Eigen form of the hydrated proton as compared to the Zundel form, both in gas-phase clusters and in the bulk; the amphiphilic character of the hydrated proton discussed in the literature has also been observed. The analyses do, however, bring forth remaining limitations in terms of the solvation structure around the hydrated proton as well as the structure of bulk water, which can guide future improvements of the method.