Quantum Mechanical Region Research Articles

Quantum Mechanics/Extremely Localized Molecular Orbital Embedding Strategy for Excited States: Coupling to Time-Dependent Density Functional Theory and Equation-of-Motion Coupled Cluster.

The QM/ELMO (quantum mechanics/extremely localized molecular orbital) method is a recently developed embedding technique in which the most important region of the system under examination is treated at fully quantum mechanical level, while the rest is described by means of transferred and frozen extremely localized molecular orbitals. In this paper, we propose the first application of the QM/ELMO approach to the investigation of excited states, and, in particular, we present the coupling of the QM/ELMO philosophy with Time-Dependent Density Functional Theory (TDDFT) and Equation-of-Motion Coupled Cluster with single and double substitutions (EOM-CCSD). The proposed TDDFT/ELMO and EOM-CCSD/ELMO strategies underwent a series of preliminary tests that were already considered for the validation of other embedding methods for excited states. The obtained results showed that the novel techniques allow the accurate description of localized excitations in large systems by only including a relatively small number of atoms in the region treated at fully quantum chemical level. Furthermore, for TDDFT/ELMO, it was also observed that (i) the method enables to avoid the presence of artificial low-lying charge-transfer states that may affect traditional TDDFT calculations, even using functionals that do not take into account long-range corrections, and (ii) the novel approach can be also successfully exploited to investigate local electronic transitions in quite large systems (e.g., reduced model of the Green Fluorescent Protein), and the accuracy of the results can be improved by including a sufficient number of chemically crucial fragments/residues in the quantum mechanical region. Finally, concerning EOM-CCSD/ELMO, it was also seen that, despite the quite crude approximation of an embedding potential given by frozen extremely localized molecular orbitals, the new strategy is able to satisfactorily account for the effects of the environment. This work paves the way to further extensions of the QM/ELMO philosophy for the study of local excitations in extended systems, suggesting the coupling of the QM/ELMO approach with other quantum chemical strategies for excited states, from the simplest ΔSCF techniques to the most advanced and computationally expensive multireferences methods.

Quantum refinement with multiple conformations: application to the P-cluster in nitrogenase.

X-ray crystallography is the main source of atomistic information on the structure of proteins. Normal crystal structures are obtained as a compromise between the X-ray scattering data and a set of empirical restraints that ensure chemically reasonable bond lengths and angles. However, such restraints are not always available or accurate for nonstandard parts of the structure, for example substrates, inhibitors and metal sites. The method of quantum refinement, in which these empirical restraints are replaced by quantum-mechanical (QM) calculations, has previously been suggested for small but interesting parts of the protein. Here, this approach is extended to allow for multiple conformations in the QM region by performing separate QM calculations for each conformation. This approach is shown to work properly and leads to improved structures in terms ofelectron-density maps and real-space difference density Z-scores. It is also shown that the quality of the structures can be gauged using QM strain energies. The approach, called ComQumX-2QM, is applied to the P-cluster in two different crystal structures of the enzyme nitrogenase, i.e. an Fe8S7Cys6 cluster, used for electron transfer. One structure is at a very high resolution (1.0 Å) and shows a mixture of two different oxidation states, the fully reduced PN state (Fe82+, 20%) and the doubly oxidized P2+ state (80%). In the original crystal structure the coordinates differed for only two iron ions, but here it is shown that the two states also show differences in other atoms of up to 0.7 Å. The second structure is at a more modest resolution, 2.1 Å, and was originally suggested to show only the one-electron oxidized state, P1+. Here, it is shown that it is rather a 50/50% mixture of the P1+ and P2+ states and that many of the Fe-Fe and Fe-S distances in the original structure were quite inaccurate (by up to 0.8 Å). This shows that the new ComQumX-2QM approach can be used to sort out what is actually seen in crystal structures with dual conformations and to give locally improved coordinates.

On the Catalytic Activity of the Engineered Coiled-Coil Heptamer Mimicking the Hydrolase Enzymes: Insights from a Computational Study.

Recently major advances were gained on the designed proteins aimed to generate biomolecular mimics of proteases. Although such enzyme-like catalysts must still suffer refinements for improving the catalytic activity, at the moment, they represent a good example of artificial enzymes to be tested in different fields. Herein, a de novo designed homo-heptameric peptide assembly (CC-Hept) where the esterase activity towards p-nitro-phenylacetate was obtained for introduction of the catalytic triad (Cys-His-Glu) into the hydrophobic matrix, is the object of the present combined molecular dynamics and quantum mechanics/molecular mechanics investigation. Constant pH Molecular Dynamics simulations on the apoform of CC-Hept suggested that the Cys residues are present in the protonated form. Molecular dynamics (MD) simulations of the enzyme–substrate complex evidenced the attitude of the enzyme-like system to retain water molecules, necessary in the hydrolytic reaction, in correspondence of the active site, represented by the Cys-His-Glu triad on each of the seven chains, without significant structural perturbations. A detailed reaction mechanism of esterase activity of CC-Hept-Cys-His-Glu was investigated on the basis of the quantum mechanics/molecular mechanics calculations employing a large quantum mechanical (QM) region of the active site. The proposed mechanism is consistent with available esterases kinetics and structural data. The roles of the active site residues were also evaluated. The deacylation phase emerged as the rate-determining step, in agreement with esterase activity of other natural proteases.

Influence of Electronic Polarization on the Spectral Density.

Accurate spectral densities are necessary for computing realistic exciton dynamics and nonlinear optical spectra of chromophores in condensed-phase environments, including multichromophore pigment-protein systems. However, due to the significant computational cost of computing spectral densities from first principles, requiring many thousands of excited-state calculations, most simulations of realistic systems rely on treating the environment as fixed-point charges. Here, using a number of representative systems ranging from solvated chromophores to the photoactive yellow protein (PYP), we demonstrate that the quantum mechanical (QM) electronic polarization of the environment is key to obtaining accurate spectral densities and line shapes within the cumulant framework. We show that the QM environment can enhance or depress the coupling of fast chromophore degrees of freedom to the energy gap, altering the electronic-vibrational coupling and the resulting vibronic progressions in the absorption spectrum. In analyzing the physical origin of peaks in the spectral density, we identify vibrational modes that couple the electron and the hole as being particularly sensitive to the QM screening of the environment. For PYP, we reveal the need for careful determination of the appropriate QM region to obtain reliable spectral densities. Our results indicate that the QM polarization of the environment can be crucial not just for excitation energies but also for electronic-vibrational coupling in complex systems with implications for the correct modeling of linear and nonlinear optical spectroscopy in the condensed phase as well as energy transfer in pigment-protein complexes.

Coupled Cluster Theory with Induced Dipole Polarizable Embedding for Ground and Excited States.

In this work, we present the theory and implementation of the coupled cluster single and double excitations (CCSD) method combined with a classical polarizable molecular mechanics force field (MMPol) based on the induced dipole model. The method is developed to compute electronic excitation energies within the state specific (SS) and linear response (LR) formalisms for the interaction of the quantum mechanical and classical regions. Furthermore, we consider an approximate expression of the correlation energy, originally developed for CCSD with implicit solvation models, where the interaction term is linear in the coupled cluster density. This approximation allows us to include the explicit contribution of the environment to the CC equations without increasing the computational effort. The test calculations on microsolvated systems, where the CCSD/MMPol method is compared to full CCSD calculations, demonstrates the reliability of this computational protocol for all interaction schemes (errors < 2%). We also show that it is important to include induced dipoles on all atom centers of the classical region and that too diffuse functions in the basis set may be problematic due to too strong interaction with the environment.

Quantum region size effect on key properties of the metal site in mutated zinc finger

Abstract We inspect the dependence of quantum mechanics/molecular mechanics hybrid calculations (QM/MM) on the quantum mechanics region size using mutated structures of the zinc finger fragment as test cases. We focus on the effect of removing the Alanine residue from the QM region on the charge transfer as well as on the stabilization energy of the metal site. All Calculations were carried out using the hybrid Local Self-Consistent Field approach (LSCF/MM) at the B3LYP-LANL2DZ/CHARMM27 level of theory. Results show that metal-ligand charge transfer is not very sensitive to the theoretical treatment of the Alanine residue when this latter is included or excluded from the QM region, and the change in the stabilization energy can give us a picture on the proximity of the Alanine to the metal ion, as well as on the most favourable region in which this residue must be included.

Effect of arsenate substitution on phosphate repository of cell: a computational study.

The structural analogy with phosphate derives arsenate into various metabolic processes associated with phosphate inside the organisms. But it is difficult to evaluate the effect of arsenate substitution on the stability of individual biological phosphate species, which span from a simpler monoester form like pyrophosphate to a more complex phosphodiester variant like DNA. In this study, we have classified the physiological phosphate esters into three different classes on the basis of their structural differences. This classification has helped us to present a concise theoretical study on the kinetic stability of phosphate analogue species of arsenate against hydrolysis. All the calculations have been carried out using QM/MM methods of our Own N-layer Integrated molecular Orbital molecular Mechanics (ONIOM). For quantum mechanical region, we have used M06-2X density functional with 6-31+G(2d,2p) basis set and for molecular mechanics we have used the AMBER force field. The calculated rate constants for hydrolysis show that none of the phosphate analogue species of arsenate has a reasonable stability against hydrolysis.

Exchange-repulsion energy in QM/EFP.

The effective fragment potential (EFP) is a quantum mechanics (QM)-based model designed to accurately describe intermolecular interactions. Hybrid QM/EFP calculations combine quantum mechanical methods with an EFP embedding to study complex systems in which many-body effects are relevant. As in EFP-only calculations, non-bonded interactions between the QM region and EFP fragments are computed as a sum of electrostatic, polarization, dispersion, and exchange-repulsion energies. The exchange-repulsion term is a computational bottleneck of the EFP calculations. Here, we present a general procedure for computing the QM/EFP exchange-repulsion interactions based on one-electron contributions to the QM Hamiltonian, by using Gaussian functions to represent localized molecular orbitals of the effective fragments. The accuracy of the exchange-repulsion and total QM/EFP interaction energies is evaluated on a diverse set of dimers, including complexes from the S22 dataset of non-covalent interactions. In most cases, the QM/EFP energies are at least as accurate as corresponding EFP energies. A simple and computationally efficient form of the introduced QM/EFP exchange-repulsion term will facilitate further developments and applications of QM/EFP methods.

QM/MM study of the reaction mechanism of sulfite oxidase

Sulfite oxidase is a mononuclear molybdenum enzyme that oxidises sulfite to sulfate in many organisms, including man. Three different reaction mechanisms have been suggested, based on experimental and computational studies. Here, we study all three with combined quantum mechanical (QM) and molecular mechanical (QM/MM) methods, including calculations with large basis sets, very large QM regions (803 atoms) and QM/MM free-energy perturbations. Our results show that the enzyme is set up to follow a mechanism in which the sulfur atom of the sulfite substrate reacts directly with the equatorial oxo ligand of the Mo ion, forming a Mo-bound sulfate product, which dissociates in the second step. The first step is rate limiting, with a barrier of 39–49 kJ/mol. The low barrier is obtained by an intricate hydrogen-bond network around the substrate, which is preserved during the reaction. This network favours the deprotonated substrate and disfavours the other two reaction mechanisms. We have studied the reaction with both an oxidised and a reduced form of the molybdopterin ligand and quantum-refinement calculations indicate that it is in the normal reduced tetrahydro form in this protein.

The Quality of the Embedding Potential Is Decisive for Minimal Quantum Region Size in Embedding Calculations: The Case of the Green Fluorescent Protein.

The calculation of spectral properties for photoactive proteins is challenging because of the large cost of electronic structure calculations on large systems. Mixed quantum mechanical (QM) and molecular mechanical (MM) methods are typically employed to make such calculations computationally tractable. This study addresses the connection between the minimal QM region size and the method used to model the MM region in the calculation of absorption properties-here exemplified for calculations on the green fluorescent protein. We find that polarizable embedding is necessary for a qualitatively correct description of the MM region, and that this enables the use of much smaller QM regions compared to fixed charge electrostatic embedding. Furthermore, absorption intensities converge very slowly with system size and inclusion of effective external field effects in the MM region through polarizabilities is therefore very important. Thus, this embedding scheme enables accurate prediction of intensities for systems that are too large to be treated fully quantum mechanically.

Improved Electrostatic Embedding for Fragment-Based Chemical Shift Calculations in Molecular Crystals.

Fragment-based methods predict nuclear magnetic resonance (NMR) chemical shielding tensors in molecular crystals with high accuracy and computational efficiency. Such methods typically employ electrostatic embedding to mimic the crystalline environment, and the quality of the results can be sensitive to the embedding treatment. To improve the quality of this embedding environment for fragment-based molecular crystal property calculations, we borrow ideas from the embedded ion method to incorporate self-consistently polarized Madelung field effects. The self-consistent reproduction of the Madelung potential (SCRMP) model developed here constructs an array of point charges that incorporates self-consistent lattice polarization and which reproduces the Madelung potential at all atomic sites involved in the quantum mechanical region of the system. The performance of fragment- and cluster-based 1H, 13C, 14N, and 17O chemical shift predictions using SCRMP and density functionals like PBE and PBE0 are assessed. The improved embedding model results in substantial improvements in the predicted 17O chemical shifts and modest improvements in the 15N ones. Finally, the performance of the model is demonstrated by examining the assignment of the two oxygen chemical shifts in the challenging γ-polymorph of glycine. Overall, the SCRMP-embedded NMR chemical shift predictions are on par with or more accurate than those obtained with the widely used gauge-including projector augmented wave (GIPAW) model.

Coupling Real-Time Time-Dependent Density Functional Theory with Polarizable Force Field.

Real-time time-dependent density functional theory (RT-TDDFT) is a powerful tool for obtaining spectroscopic observables and understanding complex, time-dependent properties. Currently, performing RT-TDDFT calculations on large, fully quantum mechanical systems is not computationally feasible. Previously, polarizable mixed quantum mechanical and molecular mechanical (QM/MMPol) models have been successful in providing accurate, yet efficient, approximations to a fully quantum mechanical system. Here we develop a coupling scheme between induced dipole based QM/MMPol and RT-TDDFT. Our approach is validated by comparing calculated spectra with both real-time and linear-response TDDFT calculations. The model developed within provides an accurate method for performing RT-TDDFT calculations on extended systems while accounting for mutual polarization between the quantum mechanical and molecular mechanical regions.

Efficient Strategy for the Calculation of Solvation Free Energies in Water and Chloroform at the Quantum Mechanical/Molecular Mechanical Level

The partitioning of solute molecules between immiscible solvents with significantly different polarities is of great importance. The polarization between the solute and solvent molecules plays an essential role in determining the solubility of the solute, which makes computational studies utilizing molecular mechanics (MM) rather difficult. In contrast, quantum mechanics (QM) can provide more reliable predictions. In this work, the partition coefficients of the side chain analogs of some amino acids between water and chloroform were computed. The QM solvation free energies were calculated indirectly via a series of MM states using the multistate Bennett acceptance ratio (MBAR) and the MM-to-QM corrections were applied at the two endpoints using thermodynamic perturbation (TP). Previously, it has been shown (Jia et al. J. Chem. Theory Comput. 2016, 12, 499-511) that this method provides the minimal variance in the results without running QM simulations. However, if there is insufficient overlap in phase space between the MM and QM Hamiltonians, this method fails. In this work, we propose, for the first time, a quantity termed the reweighting entropy that serves as a metric for the reliability of the TP calculations. If the reweighting entropy is below a certain threshold (0.65 for the solvation free energy calculations in this work), this MM-to-QM correction should be avoided and two alternative methods can be employed by either introducing a semiempirical state or conducting nonequilibrium simulations. However, the results show that the QM methods are not guaranteed to yield better results than the MM methods. Further improvement of the QM methods are imperative, especially the treatment of the van der Waals and the electrostatic interactions between the QM region and the MM region in the first shell. We also propose a scheme for the calculation of the van der Waals parameters for the solute molecules in nonaqueous solvent, which improves the quality of the computed thermodynamic properties. Furthermore, the force field parameters for the sulfur-containing molecules are also optimized.

Quantum speed limit for a relativistic electron in the noncommutative phase space

The influence of the noncommutativity on the average speed of a relativistic electron interacting with a uniform magnetic field within the minimum evolution time is investigated. We find that it is possible for the wave packet of the electron to travel faster than the speed of light in vacuum because of the noncommutativity. It is a clear signature of violating Lorentz invariance in the noncommutative relativistic quantum mechanical region.

A Multidimensional B-Spline Correction for Accurate Modeling Sugar Puckering in QM/MM Simulations.

The computational efficiency of approximate quantum mechanical methods allows their use for the construction of multidimensional reaction free energy profiles. It has recently been demonstrated that quantum models based on the neglect of diatomic differential overlap (NNDO) approximation have difficulty modeling deoxyribose and ribose sugar ring puckers and thus limit their predictive value in the study of RNA and DNA systems. A method has been introduced in our previous work to improve the description of the sugar puckering conformational landscape that uses a multidimensional B-spline correction map (BMAP correction) for systems involving intrinsically coupled torsion angles. This method greatly improved the adiabatic potential energy surface profiles of DNA and RNA sugar rings relative to high-level ab initio methods even for highly problematic NDDO-based models. In the present work, a BMAP correction is developed, implemented, and tested in molecular dynamics simulations using the AM1/d-PhoT semiempirical Hamiltonian for biological phosphoryl transfer reactions. Results are presented for gas-phase adiabatic potential energy surfaces of RNA transesterification model reactions and condensed-phase QM/MM free energy surfaces for nonenzymatic and RNase A-catalyzed transesterification reactions. The results show that the BMAP correction is stable, efficient, and leads to improvement in both the potential energy and free energy profiles for the reactions studied, as compared with ab initio and experimental reference data. Exploration of the effect of the size of the quantum mechanical region indicates the best agreement with experimental reaction barriers occurs when the full CpA dinucleotide substrate is treated quantum mechanically with the sugar pucker correction.

Convergence of Computed Aqueous Absorption Spectra with Explicit Quantum Mechanical Solvent.

For reliable condensed phase simulations, an accurate model that includes both short- and long-range interactions is required. Short- and long-range interactions can be particularly strong in aqueous solution, where hydrogen-bonding may play a large role at short range and polarization may play a large role at long range. Although short-range solute-solvent interactions such as charge transfer, hydrogen bonding, and solute-solvent polarization can be taken into account with a quantum mechanical (QM) treatment of the solvent, it is unclear how much QM solvent is necessary to accurately model interactions with different solutes. In this work, we investigate the effect of explicit QM solvent on absorption spectra computed for a series of solutes with decreasing polarity. By adjusting the boundary between QM and classical molecular mechanical solvent to include up to 400 QM water molecules, convergence of the calculated absorption spectra with respect to the size of the QM region is achieved. We find that the rate of convergence does not correlate with solute polarity when excitation energies are calculated using time-dependent density functional theory with a range-separated hybrid functional, but it does correlate with solute polarity when using configuration interaction singles. We also find that larger basis sets converge the computed spectrum with fewer QM solvent molecules. To optimize the computational cost with respect to convergence, we test a mixed basis set with more basis functions for atoms of the chromophore and the solvent molecules that are nearest to it and fewer basis functions for the atoms of the remaining solvent molecules in the QM region. Our results show that using a mixed basis set is a potentially effective way to significantly lower the computational cost while reproducing the results computed with larger basis sets.

Automated Fragmentation Polarizable Embedding Density Functional Theory (PE-DFT) Calculations of Nuclear Magnetic Resonance (NMR) Shielding Constants of Proteins with Application to Chemical Shift Predictions.

Full-protein nuclear magnetic resonance (NMR) shielding constants based on ab initio calculations are desirable, because they can assist in elucidating protein structures from NMR experiments. In this work, we present NMR shielding constants computed using a new automated fragmentation (J. Phys. Chem. B 2009, 113, 10380-10388) approach in the framework of polarizable embedding density functional theory. We extend our previous work to give both basis set recommendations and comment on how large the quantum mechanical region should be to successfully compute 13C NMR shielding constants that are comparable with experiment. The introduction of a probabilistic linear regression model allows us to substantially reduce the number of snapshots that are needed to make comparisons with experiment. This approach is further improved by augmenting snapshot selection with chemical shift predictions by which we can obtain a representative subset of snapshots that gives the smallest predicted error, compared to experiment. Finally, we use this subset of snapshots to calculate the NMR shielding constants at the PE-KT3/pcSseg-2 level of theory for all atoms in the protein GB3.

QM/MM Methods for Crystalline Defects. Part 2: Consistent Energy and Force-Mixing

We develop and analyze QM/MM (quantum/classic) hybrid methods for crystalline defects within the context of the tight-binding model. QM/MM methods employ accurate quantum mechanics (QM) models only in regions of interest (defects) and switch to computationally cheaper interatomic potential molecular mechanics (MM) models to describe the crystalline bulk. We propose new energy-based and force-based QM/MM methods, building on two principles: (i) locality of the QM model; and (ii) constructing the MM model as an explicit and controllable approximation of the QM model. This approach enables us to rigorously establish convergence rates in terms of the size of the QM region.

Quantum Speed Limit for Relativistic Spin-0 and Spin-1 Bosons on Commutative and Noncommutative Planes

Quantum speed limits of relativistic charged spin-0 and spin-1 bosons in the background of a homogeneous magnetic field are studied on both commutative and noncommutative planes. We show that, on the commutative plane, the average speeds of wave packets along the radial direction during the interval in which a quantum state is evolving from an initial state to the orthogonal final one can not exceed the speed of light, regardless of the intensities of the magnetic field. However, due to the noncommutativity, the average speeds of the wave packets on noncommutative plane will exceed the speed of light in vacuum provided the intensity of the magnetic field is strong enough. It is a clear signature of violating Lorentz invariance in the relativistic quantum mechanics region.

Convergence of Excitation Energies in Mixed Quantum and Classical Solvent: Comparison of Continuum and Point Charge Models

Mixed quantum mechanical (QM)/classical methods provide a computationally efficient approach to modeling both ground and excited states in the condensed phase. To accurately model short-range interactions, some amount of the environment can be included in the QM region, whereas a classical model can treat long-range interactions to maintain computational affordability. The best computational protocol for these mixed QM/classical methods can be determined by examining convergence of molecular properties. Here, we compare molecular mechanical (MM) fixed point charges to a polarizable continuum model (PCM) for computing electronic excitations in solution. We computed the excitation energy of three pairs of neutral/anionic molecules in aqueous solvent, including up to 250 water molecules in the QM region. Interestingly, the convergence is similar for MM point charges and a PCM, with convergence achieved when at least one full solvation shell is treated with QM. Although the van der Waals (VDW) definition of the PCM cavity is adequate with small amounts of QM solvent, larger QM solvent layers had gaps in the VDW PCM cavity, leading to asymptotically incorrect excitation energies. Given that the VDW cavity leads to unphysical solute-solvent interactions, we advise using a solvent-excluded surface cavity for QM/PCM calculations that include QM solvent.